Synchronization Properties of Network Motifs

We address the problem of understanding the variable abundance of 3-node and 4-node subgraphs (motifs) in complex networks from a dynamical point of view. As a criterion in the determination of the functional significance of a n-node subgraph, we propose an analytic method to measure the stability of the synchronous state (SSS) the subgraph displays. We show that, for undirected graphs, the SSS is correlated with the relative abundance, while in directed graphs the correlation exists only for some specific motifs.Comment: 7 pages, 3 figure

Failed tidal disruption events and X-ray flares from the Galactic Centre

The process of tidal disruption of stars by a supermassive black hole provides luminous UV and soft X-ray flares with peak luminosities of 481046 erg\u2009s 121 and duration of a few months. As part of a wider exploration of the effects of stellar rotation on the outcome of a TDE, we have performed hydrodynamical simulations of the disruption of a rotating star whose spin axis is opposite to the orbital axis. Such a retrograde rotation makes the star more resilient to tidal disruption, so that, even if its orbit reaches the formal tidal radius, it actually stays intact after the tidal encounter. However, the outer layers of the star are initially stripped away from the core, but then fall back on to the star itself, producing a newly formed accretion disc around the star. We estimate that the accretion rate on to the star would be strongly super-Eddington (for the star) and would result in an X-ray flare with luminosity of the order of 481040 erg\u2009s 121 and duration of a few months. We speculate that such events might be responsible for the known X-ray flares from Sgr A* in the recent past

Grain growth in the envelopes and disks of Class I protostars

We present new 3 mm ATCA data of two Class I Young Stellar Objects in the Ophiucus star forming region: Elias29 and WL12. For our analysis we compare them with archival 1.1 mm SMA data. In the (u,v) plane the two sources present a similar behavior: a nearly constant non-zero emission at long baselines, which suggests the presence of an unresolved component and an increase of the fluxes at short baselines, related to the presence of an extended envelope. Our data analysis leads to unusually low values of the spectral index $\alpha_{\rm 1.1-3mm}$, which may indicate that mm-sized dust grains have already formed both in the envelopes and in the disk-like structures at such early stages. To explore the possible scenarios for the interpretation of the sources we perform a radiative transfer modeling using a Monte Carlo code, in order to take into account possible deviations from the Rayleigh-Jeans and optically thin regimes. Comparison between the model outputs and the observations indicates that dust grains may form aggregates up to millimeter size already in the inner regions of the envelopes of Class I YSOs. Moreover, we conclude that the embedded disk-like structures in our two Class Is are probably very compact, in particular in the case of WL12, with outer radii down to tens of AU.Comment: 12 pages, 8 figures, Accepted for publication in A&

Multiple tidal disruption flares in the active galaxy IC 3599

Tidal disruption events occur when a star passes too close to a massive black hole and it is totally ripped apart by tidal forces. It may also happen that the star is not close enough to the black hole to be totally disrupted and a less dramatic event might happen. If the stellar orbit is bound and highly eccentric, just like some stars in the centre of our own Galaxy, repeated flares should occur. When the star approaches the black hole tidal radius at periastron, matter might be stripped resulting in lower intensity outbursts recurring once every orbital period. We report on Swift observations of a recent bright flare from the galaxy IC 3599 hosting a middle-weight black hole, where a possible tidal disruption event was observed in the early 1990s. By light curve modelling and spectral fitting we can consistently account for the events as the non-disruptive tidal stripping of a star into a highly eccentric orbit. The recurrence time is 9.5 yr. IC 3599 is also known to host a low-luminosity active galactic nucleus. Tidal stripping from this star over several orbital passages might be able to spoon-feed also this activity.Comment: Accepted for publication to Astronomy & Astrophysic

Serum creatine kinase isoenzymes in children with osteogenesis imperfecta

This study evaluates serum creatine kinase isoenzyme activity in children with osteogenesis imperfecta to determine its usefulness as a biochemical marker during treatment with bisphosphonate. The changes of creatine kinase (CK) isoenzyme activity during and after discontinuation therapy were observed. These results could be useful in addressing over-treatment risk prevention. Introduction The brain isoenzyme of creatine kinase (CKbb) is highly expressed in mature osteoclasts during osteoclastogenesis, thus plays an important role in bone resorption. We previously identified high serum CKbb levels in 18 children with osteogenesis imperfect (OI) type 1 treated for 1 year with bisphosphonate (neridronate). In the present study, serum CK isoenzymes were evaluated in the same children with continuous versus discontinued neridronate treatment over a further 2-year follow-up period. Methods This study included 18 children with OI type 1, 12 with continued (group A) and 6 with ceased (group B) neridronate treatment. Auxological data, serum biochemical markers of bone metabolism, bone mineral density z-score, and serum total CK and isoenzyme activities were determined in both groups. Results Serum CKbb was progressively and significantly increased in group A (p < 0.004) but rapidly decreased to undetectable levels in group B. In both groups, the cardiac muscle creatine kinase isoenzyme (CKmb) showed a marked decrease, while serum C-terminal telopeptide (CTx) levels were almost unchanged. Conclusions This study provides evidence of the cumulative effect of neridronate administration in increasing serum CKbb levels and the reversible effect after its discontinuation. This approach could be employed for verifying the usefulness of serum CKbb as a biochemical marker in patients receiving prolonged bisphosphonate treatment. Moreover, the decreased serum CKmb levels suggest a systemic effect of these drugs

Black hole mergers: the first light

The coalescence of supermassive black hole binaries occurs via the emission of gravitational waves, that can impart a substantial recoil to the merged black hole. We consider the energy dissipation, that results if the recoiling black hole is surrounded by a thin circumbinary disc. Our results differ significantly from those of previous investigations. We show analytically that the dominant source of energy is often potential energy, released as gas in the outer disc attempts to circularize at smaller radii. Thus, dimensional estimates, that include only the kinetic energy gained by the disc gas, underestimate the real energy loss. This underestimate can exceed an order of magnitude, if the recoil is directed close to the disc plane. We use three dimensional Smooth Particle Hydrodynamics (SPH) simulations and two dimensional finite difference simulations to verify our analytic estimates. We also compute the bolometric light curve, which is found to vary strongly depending upon the kick angle. A prompt emission signature due to this mechanism may be observable for low mass (10^6 Solar mass) black holes whose recoil velocities exceed about 1000 km/s. Emission at earlier times can mainly result from the response of the disc to the loss of mass, as the black holes merge. We derive analytically the condition for this to happen.Comment: 16 pages, accepted by MNRAS. Animations of the simulations are available at http://jilawww.colorado.edu/~pja/recoil.htm

Stellar disruption by a supermassive black hole: is the light curve really proportional to t−5/3t^{-5/3}?

In this paper we revisit the arguments for the basis of the time evolution of the flares expected to arise when a star is disrupted by a supermassive black hole. We present a simple analytic model relating the lightcurve to the internal density structure of the star. We thus show that the standard lightcurve proportional to $t^{-5/3}$ only holds at late times. Close to the peak luminosity the lightcurve is shallower, deviating more strongly from $t^{-5/3}$ for more centrally concentrated (e.g. solar--type) stars. We test our model numerically by simulating the tidal disruption of several stellar models, described by simple polytropic spheres with index $\gamma$. The simulations agree with the analytical model given two considerations. First, the stars are somewhat inflated on reaching pericentre because of the effective reduction of gravity in the tidal field of the black hole. This is well described by a homologous expansion by a factor which becomes smaller as the polytropic index becomes larger. Second, for large polytropic indices wings appear in the tails of the energy distribution, indicating that some material is pushed further away from parabolic orbits by shocks in the tidal tails. In all our simulations, the $t^{-5/3}$ lightcurve is achieved only at late stages. In particular we predict that for solar type stars, this happens only after the luminosity has dropped by at least two magnitudes from the peak. We discuss our results in the light of recent observations of flares in otherwise quiescent galaxies and note the dependence of these results on further parameters, such as the star/hole mass ratio and the stellar orbit.Comment: 10 pages, 10 figures, MNRAS accepte

Obscuring and feeding supermassive black holes with evolving nuclear star clusters

Recently, high resolution observations with the help of the near-infrared adaptive optics integral field spectrograph SINFONI at the VLT proved the existence of massive and young nuclear star clusters in the centres of a sample of Seyfert galaxies. With the help of high resolution hydrodynamical simulations with the PLUTO-code, we follow the evolution of such clusters, especially focusing on mass and energy feedback from young stars. This leads to a filamentary inflow of gas on large scales (tens of parsec), whereas a turbulent and very dense disc builds up on the parsec scale. Here, we concentrate on the long-term evolution of the nuclear disc in NGC 1068 with the help of an effective viscous disc model, using the mass input from the large scale simulations and accounting for star formation in the disc. This two-stage modelling enables us to connect the tens of parsec scale region (observable with SINFONI) with the parsec scale environment (MIDI observations). At the current age of the nuclear star cluster, our simulations predict disc sizes of the order of 0.8 to 0.9 pc, gas masses of 1.0e6 Msun and mass transfer rates through the inner boundary of 0.025 Msun/yr in good agreement with values derived from observations.Comment: 6 pages, 3 figures, to appear in the proceedings of the IAU General Assembly 2009, Rio de Janeiro, S267 Co-evolution of Central Black Holes and Galaxie

The puzzling source IGR J17361-4441 in NGC 6388: a possible planetary tidal disruption event

On 2011 August 11, INTEGRAL discovered the hard X-ray source IGR J17361-4441 near the centre of the globular cluster NGC 6388. Follow up observations with Chandra showed the position of the transient was inconsistent with the cluster dynamical centre, and thus not related to its possible intermediate mass black hole. The source showed a peculiar hard spectrum (Gamma \approx 0.8) and no evidence of QPOs, pulsations, type-I bursts, or radio emission. Based on its peak luminosity, IGR J17361-4441 was classified as a very faint X-ray transient, and most likely a low-mass X-ray binary. We re-analysed 200 days of Swift/XRT observations, covering the whole outburst of IGR J17361-4441 and find a t^{-5/3} trend evident in the light curve, and a thermal emission component that does not evolve significantly with time. We investigate whether this source could be a tidal disruption event, and for certain assumptions find an accretion efficiency epsilon \approx 3.5E-04 (M_{Ch}/M) consistent with a massive white dwarf, and a disrupted minor body mass M_{mb}=1.9E+27(M/M_{Ch}) g in the terrestrial-icy planet regime. These numbers yield an inner disc temperature of the order kT_{in} \approx 0.04 keV, consistent with the blackbody temperature of kT_{in} \approx 0.08 keV estimated by spectral fitting. Although the density of white dwarfs and the number of free-floating planets are uncertain, we estimate the rate of planetary tidal disruptions in NGC 6388 to be in the range 3E-06 to 3E-04 yr^{-1}. Averaged over the Milky Way globular clusters, the upper limit value corresponds to 0.05 yr^{-1}, consistent with the observation of a single event by INTEGRAL and Swift.Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal on 2014 July 16; 9 pages, 5 figures. Added references; corrected typo

Effects of photoevaporation on protoplanetary disc &#8216;isochrones&#8217;

Protoplanetary discs are the site of star and planet formation, and their evolution and consequent dispersal deeply affect the formation of planetary systems. In the standard scenario they evolve on time-scales similar to Myr due to the viscous transport of angular momentum. The analytical self-similar solution for their evolution predicts specific disc isochrones in the accretion rate-disc mass plane. However, photoevaporation by radiation emitted by the central star is likely to dominate the gas disc dispersal of the innermost region, introducing another (shorter) time-scale for this process. In this paper, we include the effect of internal (X and EUV) photoevaporation on the disc evolution, finding numerical solutions for a population of protoplanetary discs. Our models naturally reproduce the expected quick dispersal of the inner region of discs when their accretion rates match the rate of photoevaporative mass loss, in line with previous studies. We find that photoevaporation preferentially removes the lightest discs in the sample. The net result is that, counter-intuitively, photoevaporation increases the average disc mass in the sample, by dispersing the lightest discs. At the same time, photoevaporation also reduces the mass accretion rate by cutting the supply of material from the outer to the inner disc. In a purely viscous framework, this would be interpreted as the result of a longer viscous evolution, leading to an overestimate of the disc age. Our results thus show that photoevaporation is a necessary ingredient to include when interpreting observations of large disc samples with measured mass accretion rates and disc masses. Photoevaporation leaves a characteristic imprint on the shape of the isochrone. Accurate data of the accretion rate-disc mass plane in the low disc mass region therefore give clues on the typical photoevaporation rate