BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jorge Cuadra

Jorge Cuadra appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.EPastro-ph.HEastro-ph.SRastro-ph.COgr-qcastro-ph.IMcond-mat.mes-hallphysics.optics

Catalog footprint

What is connected

26works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2022arXiv

Galactic center gamma-ray production by cosmic rays from stellar winds and Sgr A East

The High Energy Stereoscopic System (HESS), the Major Atmospheric Gamma-ray Imaging Cherenkov Telescope (MAGIC), and the Very Energetic Radiation Imaging Telescope Array System (VERITAS) have observed diffuse gamma-ray emission strongly correlated with the central molecular zone in the Galactic center. The most accepted scenario to generate this emission is via a hadronic interaction between cosmic rays (CRs) and ambient gas, where CRs are accelerated from a central and continuous source of 1 PeV protons (PeVatron). We explore the influence of the three-dimensional (3D) shape of the central molecular zone on the indirect observation of the CR energy density via gamma-ray detection. We simulated synthetic gamma-ray maps using a CR diffusion model with spherical injection, one isotropic diffusion coefficient, no advection, and mono-energetic particles of 1 PeV. Also, we used two different 3D gas distributions considering the observed gas column density, both with and without an inner cavity. We find that when using a persistent CR source, a disk-like gas distribution is needed to reproduce the existing CR indirect observations. This is in agreement with the continuous gas distribution implied by some dynamical models and studies based on the comparison of emission and absorption molecular lines. However, it contradicts several models of the central molecular zone, which imply that this structure has a significant inner cavity. This tension can be reconciled by an additional, impulsive CR injection. If the central molecular zone has a cavity, a composite CR population, coming from the stellar winds of the Wolf-Rayet stars in the central 0.5 pc and the supernova Sgr A East, produces a good match to the observed gamma-ray morphology in the Galactic center.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

New self-consistent wind parameters fitting optical spectra of O-type stars observed with HERMES spectrograph

We perform spectral fitting for a set of O-type stars based on self-consistent wind solutions, which provide mass-loss rate and velocity profiles directly derived from the initial stellar parameters. The great advantage of this self-consistent spectral fitting is therefore the reduction of the number of free parameters to be tuned. Self-consistent values for the line-force parameters (k,alpha,delta) and subsequently for the mass-loss rate and terminal velocity are provided by the m-CAK prescription introduced in Paper I, updated in this work with improvements such as a temperature structure for the wind, self-consistently evaluated from the line-acceleration. Synthetic spectra are calculated using the radiative transfer code FASTWIND, replacing the classical beta-law for our new calculated velocity profiles. We found that self-consistent m-CAK solutions provide values for theoretical mass-loss rates on the order of the most recent predictions of other studies. From here, we generate synthetic spectra with self-consistent hydrodynamics to fit and obtain a new set of stellar and wind parameters for our sample of O-type stars whose spectra was taken with the high resolution echelle spectrograph HERMES (R=85000). We find a satisfactory global fit for our observations, with good accuracy for photospheric He I and He II lines and a quite acceptable fit for H lines. Although this self-consistent spectral analysis is currently constrained in the optical wavelength range only, this is an important step towards the determination of stellar and wind parameters without using a beta-law. Given these results, we expect that the values introduced here should be helpful for future studies about the stars constituting this sample, together with the prospective that the m-CAK self-consistent prescription be extended to numerous studies about massive stars in future.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

3D simulations of clump formation in stellar wind collisions

The inner parsec of our Galaxy contains tens of Wolf-Rayet stars whose powerful outflows are constantly interacting while filling the region with hot, diffuse plasma. Theoretical models have shown that, in some cases, the collision of stellar winds can generate cold, dense material in the form of clumps. However, their formation process and properties are not well understood yet. In this work we present, for the first time, a statistical study of the clump formation process in unstable wind collisions. We study systems with dense outflows (${\sim}10^{-5}\rm\ M_{\odot}\ yr^{-1}$), wind speeds of $500$-$1500\rm\ km\ s^{-1}$, and stellar separations of ${\sim}20$-$200\rm\ au$. We develop 3D high resolution hydrodynamical simulations of stellar wind collisions with the adaptive-mesh refinement grid-based code Ramses. We aim to characterise the initial properties of clumps that form through hydrodynamic instabilities, mostly via the non-linear thin shell instability (NTSI). Our results confirm that more massive clumps are formed in systems whose winds are close to the transition between the radiative and adiabatic regimes. Increasing either the wind speed or the degree of asymmetry increases the dispersion of the clump mass and ejection speed distributions. Nevertheless, the most massive clumps are very light (${\sim}10^{-3}$-$10^{-2}\rm\ M_{\oplus}$), about three orders of magnitude less massive than theoretical upper limits. Applying these results to the Galactic Centre we find that clumps formed through the NTSI should not be heavy enough either to affect the thermodynamic state of the region or to survive for long enough to fall onto the central super-massive black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Binary-induced spiral arms inside the disc cavity of AB Aurigae

In this work we demonstrate that the inner spiral structure observed in AB Aurigae can be created by a binary star orbiting inside the dust cavity. We find that a companion with a mass-ratio of 0.25, semi-major axis of 40 au, eccentricity of 0.5, and inclination of 90° produces gaseous spirals closely matching the ones observed in $^{12}$CO (2-1) line emission. Based on dust dynamics in circumbinary discs (Poblete, Cuello, and Cuadra 2019), we constrain the inclination of the binary with respect to the circumbinary disc to range between 60° and 90°. We predict that the stellar companion is located roughly 0.18 arcsec from the central star towards the east-southeast, above the plane of the disc. Should this companion be detected in the near future, our model indicates that it should be moving away from the primary star at a rate of 6 mas/yr on the plane of the sky. Since our companion is inclined, we also predict that the spiral structure will appear to change with time, and not simply co-rotate with the companion.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Feedback-limited Accretion: Variable Luminosity from Growing Planets

Planets form in discs of gas and dust around stars, and continue to grow by accretion of disc material while available. Massive planets clear a gap in their protoplanetary disc, but can still accrete gas through a circumplanetary disk. For high enough accretion rates the planet should be detectable at infrared wavelengths. As the energy of the gas accreted on to the planet is released, the planet surroundings heat up in a feedback process. We aim to test how this planet feedback affects the gas in the coorbital region and the accretion rate itself. We modified the 2D code FARGO-AD to include a prescription for the accretion and feedback luminosity of the planet and use it to model giant planets on 10 au circular and eccentric orbits around a solar mass star. We find that this feedback reduces but does not halt the accretion on to the planet, although this result might depend on the near-coincident radial ranges where both recipes are implemented. Our simulations also show that the planet heating gives the accretion rate a stochastic variability with an amplitude $Δ\dot{M}_p \sim 0.1 \dot{M}_p$. A planet on an eccentric orbit ($e=0.1)$ presents a similar variability amplitude, but concentrated on a well-defined periodicity of half the orbital period and weaker broadband noise, potentially allowing observations to discriminate between both cases. Finally, we find that the heating of the coorbital region by the planet feedback alters the gas dynamics, reducing the difference between its orbital velocity and the Keplerian motion at the edge of the gap, which can have important consequences for the formation of dust rings.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How Jupiters save or destroy inner Neptunes around evolved stars

In about 6 Giga years our Sun will evolve into a red giant and finally end its life as a white dwarf. This stellar metamorphosis will occur to virtually all known host stars of exo-planetary systems and is therefore crucial for their final fate. It is clear that the innermost planets will be engulfed and evaporated during the giant phase and that planets located farther out will survive. However, the destiny of planets in-between, at ~ 1 - 10 au, has not yet been investigated with a multi-planet tidal treatment. We here combine for the first time multi-planet interactions, stellar evolution, and tidal effects in an N-body code to study the evolution of a Neptune-Jupiter planetary system. We report that the fate of the Neptune-mass planet, located closer to the star than the Jupiter-mass planet, can be very different from the fate of a single Neptune. The simultaneous effects of gravitational interactions, mass loss and tides can drive the planetary system towards mean motion resonances. Crossing these resonances affects particularly the eccentricity of the Neptune and thereby also its fate, which can be engulfment, collision with the Jupiter-mass planet, ejection from the system, or survival at a larger separation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Colliding Winds in and around the Stellar Group IRS 13E at the Galactic Center

IRS~13E is an enigmatic compact group of massive stars located in projection only 3.6 arcseconds away from Sgr A*. This group has been suggested to be bounded by an intermediate-mass black hole (IMBH). We present a multi-wavelength study of the group and its interplay with the environment. Based on Chandra observations, we find the X-ray spectrum of IRS~13E can be well characterized by an optically thin thermal plasma. The emission peaks between two strongly mass-losing Wolf-Rayet stars of the group. These properties can be reasonably well reproduced by simulated colliding winds of these two stars. However, this scenario under-predicts the X-ray intensity in outer regions. The residual emission likely results from the ram-pressure confinement of the IRS~13E group wind by the ambient medium and is apparently associated with a shell-like warm gas structure seen in Pa-alpha and in ALMA observations. These latter observations also show strongly peaked thermal emission with unusually large velocity spread between the two stars. These results indicate that the group is colliding with the bar of the dense cool gas mini-spiral around Sgr A*. The extended X-ray morphology of IRS~13E and its association with the bar further suggest that the group is physically much farther away than the projected distance from Sgr A*. The presence of an IMBH, while favorable to keep the stars bound together, is not necessary to explain the observed stellar and gas properties of IRS~13E.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Stellar winds pump the heart of the Milky Way

The central super-massive black hole of the Milky Way, Sgr A*, accretes at a very low rate making it a very underluminous galactic nucleus. Despite the tens of Wolf-Rayet stars present within the inner parsec supplying ${\sim}10^{-3}\rm\ M_{\odot}\ yr^{-1}$ in stellar winds, only a negligible fraction of this material ($<10^{-4}$) ends up being accreted onto Sgr A*. The recent discovery of cold gas (${\sim}10^4\rm\ K$) in its vicinity raised questions about how such material could settle in the hostile (${\sim}10^7\rm\ K$) environment near Sgr A*. In this work we show that the system of mass-losing stars blowing winds can naturally account for both the hot, inefficient accretion flow, as well as the formation of a cold disk-like structure. We run hydrodynamical simulations using the grid-based code Ramses starting as early in the past as possible to observe the state of the system at the present time. Our results show that the system reaches a quasi-steady state in about ${\sim}500\rm\ yr$ with material being captured at a rate of ${\sim}10^{-6}\rm\ M_{\odot}\ yr^{-1}$ at scales of ${\sim}10^{-4}\rm\ pc$, consistent with the observations and previous models. However, on longer timescales ($\gtrsim3000\rm\ yr$) the material accumulates close to the black hole in the form of a disk. Considering the duration of the Wolf-Rayet phase (${\sim}10^5\rm\ yr$), we conclude that this scenario likely has already happened, and could be responsible for the more active past of Sgr A*, and/or its current outflow. We argue that the hypothesis of the mass-losing stars being the main regulator of the activity of the black hole deserves further consideration.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Visualizing plasmon-exciton polaritons at the nanoscale using electron microscopy

Polaritons are compositional light-matter quasiparticles that have recently enabled remarkable breakthroughs in quantum and nonlinear optics, as well as in material science. Despite the enormous progress, however, a direct nanometer-scale visualization of polaritons has remained an open challenge. Here, we demonstrate that plasmon-exciton polaritons, or plexcitons, generated by a hybrid system composed of an individual silver nanoparticle and a few-layer transition metal dichalcogenide can be spectroscopically mapped with nanometer spatial resolution using electron energy loss spectroscopy in a scanning transmission electron microscope. Our experiments reveal important insights about the coupling process, which have not been reported so far. These include nanoscale variation of Rabi splitting and plasmon-exciton detuning, as well as absorption-dominated extinction signals, which in turn provide the ultimate evidence for the plasmon-exciton hybridization in the strong coupling regime. These findings pioneer new possibilities for in-depth studies of polariton-related phenomena with nanometer spatial resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Inner mean-motion resonances with eccentric planets: A possible origin for exozodiacal dust clouds

High levels of dust have been detected in the immediate vicinity of many stars, both young and old. A promising scenario to explain the presence of this short-lived dust is that these analogues to the Zodiacal cloud (or exozodis) are refilled in situ through cometary activity and sublimation. As the reservoir of comets is not expected to be replenished, the presence of these exozodis in old systems has yet to be adequately explained. It was recently suggested that mean-motion resonances (MMR) with exterior planets on moderately eccentric ($\mathrm{e_p}\gtrsim 0.1$) orbits could scatter planetesimals on to cometary orbits with delays of the order of several 100 Myr. Theoretically, this mechanism is also expected to sustain continuous production of active comets once it has started, potentially over Gyr-timescales. We aim here to investigate the ability of this mechanism to generate scattering on to cometary orbits compatible with the production of an exozodi on long timescales. We combine analytical predictions and complementary numerical N-body simulations to study its characteristics. We show, using order of magnitude estimates, that via this mechanism, low mass discs comparable to the Kuiper Belt could sustain comet scattering at rates compatible with the presence of the exozodis which are detected around Solar-type stars, and on Gyr timescales. We also find that the levels of dust detected around Vega could be sustained via our proposed mechanism if an eccentric Jupiter-like planet were present exterior to the system's cold debris disc.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Resolving the Planetesimal Belt of HR 8799 with ALMA

The star HR 8799 hosts one of the largest known debris discs and at least four giant planets. Previous observations have found evidence for a warm belt within the orbits of the planets, a cold planetesimal belt beyond their orbits and a halo of small grains. With the infrared data, it is hard to distinguish the planetesimal belt emission from that of the grains in the halo. With this in mind, the system has been observed with ALMA in band 6 (1.34 mm) using a compact array format. These observations allow the inner edge of the planetesimal belt to be resolved for the first time. A radial distribution of dust grains is fitted to the data using an MCMC method. The disc is best fit by a broad ring between $145^{+12}_{-12}$ AU and $429^{+37}_{-32}$ AU at an inclination of $40^{+5}_{-6}$° and a position angle of $51^{+8}_{-8}$°. A disc edge at ~145 AU is too far out to be explained simply by interactions with planet b, requiring either a more complicated dynamical history or an extra planet beyond the orbit of planet b.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

SDSS J0159+0105: A Radio-Quiet Quasar with a Centi-Parsec Supermassive Black Hole Binary Candidate

We report a candidate centi-parsec supermassive black hole binary (SMBHB) in the radio-quiet quasar SDSS J0159+0105 at z=0.217. With a modified lomb-scargle code GLSdeDRW and the auto-correlation analysis ACF, we detect two significant (at P>99%) periodic signals at ~741 day and ~1500 day from the 8.1-year Catalina V-band light curve of this quasar. The period ratio, which is close to 1:2, is typical of a black-hole binary system with a mass ratio of 0.05<q<0.8 according to recent numerical simulations. SDSS J0159+0105 has two SDSS spectroscopic observations separated by ~10 years. There is a significant change in the broad H-beta profile between the two epochs, which can be explained by a single broad-line region (BLR) around the binary system illuminated by the aforementioned mini-disks, or a stream of gas flowing from the circumbinary disk to one of the SMBHs. From the single BLR assumption and the orbital period t_orb ~1500 day, we estimate the total virial masses of M_SMBHB ~ 1.3x10^8 M_sun, the average distances of BLR of ~0.04pc (~50 light-day, with +/-0.3 dex uncertainty), and a SMBHB separation of d= (0.01pc)M_{8,tot}^{1/3} (T_rest/3.3yr)^{2/3} ~ 0.013 pc (15 light-day). Based on analytical work, the postulated circumbinary disk has an inner radius of 2d = 0.026 pc (30 light-day). SDSS J0159+0105 also displays unusual spectral energy distribution. The unique properties of SDSS J0159+0105 are consistent with it being a centi-parsec SMBHB. The GLSdeDRW code link: http://butler.lab.asu.edu/qso_period/

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spiral waves triggered by shadows in transition disks

Circumstellar asymmetries such as central warps have recently been shown to cast shadows on outer disks. We investigate the hydrodynamical consequences of such variable illumination on the outer regions of a transition disk, and the development of spiral arms. Using 2D simulations, we follow the evolution of a gaseous disk passively heated by the central star, under the periodic forcing of shadows with an opening angle of $\sim$28$^\circ$. With a lower pressure under the shadows, each crossing results in a variable azimuthal acceleration, which in time develops into spiral density waves. Their pitch angles evolve from $Π\sim 15^\circ-22^\circ$ at the onset, to $\sim$11$^\circ$-14$^\circ$, over $\sim$65~AU to 150~AU. Self-gravity enhances the density contrast of the spiral waves, as also reported previously for spirals launched by planets. Our control simulations with unshadowed irradiation do not develop structures, except for a different form of spiral waves seen at later times only in the gravitationally unstable control case. Scattered light predictions in the $H$-band show that such illumination spirals should be observable. We suggest that spiral arms in the case-study transition disk HD~142527 could be explained as a result of shadowing from the tilted inner disk.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accretion kinematics through the warped transition disk in HD142527 from resolved CO(6-5) observations

The finding of residual gas in the large central cavity of the HD142527 disk motivates questions on the origin of its non-Keplerian kinematics, and possible connections with planet formation. We aim to understand the physical structure that underlies the intra-cavity gaseous flows, guided by new molecular-line data in CO(6-5) with unprecedented angular resolutions. Given the warped structure inferred from the identification of scattered-light shadows cast on the outer disk, the kinematics are consistent, to first order, with axisymmetric accretion onto the inner disk occurring at all azimuth. A steady-state accretion profile, fixed at the stellar accretion rate, explains the depth of the cavity as traced in CO isotopologues. The abrupt warp and evidence for near free-fall radial flows in HD 142527 resemble theoretical models for disk tearing, which could be driven by the reported low mass companion, whose orbit may be contained in the plane of the inner disk. The companion's high inclination with respect to the massive outer disk could drive Kozai oscillations over long time-scales; high-eccentricity periods may perhaps account for the large cavity. While shadowing by the tilted disk could imprint an azimuthal modulation in the molecular-line maps, further observations are required to ascertain the significance of azimuthal structure in the density field inside the cavity of HD142527.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hydrodynamic and radiative transfer modeling of X-ray emission from colliding WR winds: WR 140 & the Galactic center

Colliding Wolf-Rayet (WR) winds produce thermal X-ray emission widely observed by X-ray telescopes. In wide WR+O binaries, such as WR 140, the X-ray flux is tied to the orbital phase, and is a direct probe of the winds' properties. In the Galactic center, $\sim$30 WRs orbit the super massive black hole (SMBH) within $\sim$10", leading to a smorgasbord of wind-wind collisions. To model the X-ray emission of WR 140 and the Galactic center, we perform 3D hydrodynamic simulations to trace the complex gaseous flows, and then carry out 3D radiative transfer calculations to compute the variable X-ray spectra. The model WR 140 RXTE light curve matches the data well for all phases except the X-ray minimum associated with periastron, while the model spectra agree with the RXTE hardness ratio and the shape of the Suzaku observations throughout the orbit. The Galactic center model of the Chandra flux and spectral shape match well in the region r$<$3", but the model flux falls off too rapidly beyond this radius.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Precession and accretion in circumbinary discs: The case of HD 104237

We present the results of smoothed particle hydrodynamics (SPH) simulations of the disc around the young, eccentric stellar binary HD 104237. We find that the binary clears out a large cavity in the disc, driving a significant eccentricity at the cavity edge. This then precesses around the binary at a rate of $\dot{\varpi} = 0.48^{\circ}T_{\mathrm{b}}^{-1}$, which for HD 104237 corresponds to 40 years. We find that the accretion pattern into the cavity and onto the binary changes with this precession, resulting in a periodic accretion variability driven purely by the physical parameters of the binary and its orbit. For each star we find that this results in order of magnitude changes in the accretion rate. We also find that the accretion variability allows the primary to accrete gas at a higher rate than the secondary for approximately half of each precession period. Using a large number of 3-body integrations of test particles orbiting different binaries, we find good agreement between the precession rate of a test particle and our SPH disc precession. These rates also agree very well with the precession rates predicted by the analytic theory of Leung & Lee (2013), showing that their prescription can be accurately used to predict long-term accretion variability timescales for eccentric binaries accreting from a disc. We discuss the implications of our result, and suggest that this process provides a viable way of preserving unequal mass ratios in accreting eccentric binaries in both the stellar and supermassive black hole regimes

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Probing dark matter crests with white dwarfs and IMBHs

White dwarfs (WDs) are the most promising captors of dark matter (DM) particles in the crests that are expected to build up in the cores of dense stellar clusters. The DM particles could reach sufficient densities in WD cores to liberate energy through self-annihilation. The extinction associated with our Galactic Centre, the most promising region where to look for such effects, makes it impossible to detect the potential associated luminosity of the DM-burning WDs. However, in smaller stellar systems which are close enough to us and not heavily extincted, such as $ω-$Cen, we may be able to detect DM-burning WDs. We investigate the prospects of detection of DM-burning WDs in a stellar cluster harbouring an IMBH, which leads to higher densities of DM at the centre as compared with clusters without one. We calculate the capture rate of WIMPs by a WD around an IMBH and estimate the luminosity that a WD would emit depending on its distance to the center of the cluster. Direct-summation $N-$body simulations of $ω-$Cen yield a non-negligible number of WDs in the range of radii of interest. We apply our assumption to published HST/ACS observations of stars in the center of $ω-$Cen to search for DM burning WDs and, although we are not able to identify any evident candidate because of crowding and incompleteness, we proof that their bunching up at high luminosities would be unique. We predict that DM burning will lead to a truncation of the cooling sequence at the faint end. The detection of DM burning in future observations of dense stellar clusters, such as globular clusters or ultra-compact dwarf galaxies could allow us to probe different models of DM distributions and characteristics such as the DM particle scattering cross section on nucleons. On the other hand, if DM-burning WDs really exist, their number and properties could give hints to the existence of IMBHs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Protoplanetary disks including radiative feedback from accreting planets

While recent observational progress is converging on the detection of compact regions of thermal emission due to embedded protoplanets, further theoretical predictions are needed to understand the response of a protoplanetary disk to the planet formation radiative feedback. This is particularly important to make predictions for the observability of circumplanetary regions. In this work we use 2D hydrodynamical simulations to examine the evolution of a viscous protoplanetary disk in which a luminous Jupiter-mass planet is embedded. We use an energy equation which includes the radiative heating of the planet as an additional mechanism for planet formation feedback. Several models are computed for planet luminosities ranging from $10^{-5}$ to $10^{-3}$ Solar luminosities. We find that the planet radiative feedback enhances the disk's accretion rate at the planet's orbital radius, producing a hotter and more luminous environement around the planet, independently of the prescription used to model the disk's turbulent viscosity. We also estimate the thermal signature of the planet feedback for our range of planet luminosities, finding that the emitted spectrum of a purely active disk, without passive heating, is appreciably modified in the infrared. We simulate the protoplanetary disk around HD 100546 where a planet companion is located at about 68 au from the star. Assuming the planet mass is 5 Jupiter masses and its luminosity is $\sim 2.5 \times 10^{-4} \, L_\odot$, we find that the radiative feedback of the planet increases the luminosity of its $\sim 5$ au circumplanetary disk from $10^{-5} \, \rm L_\odot$ (without feedback) to $10^{-3} \, \rm L_\odot$, corresponding to an emission of $\sim 1 \, \rm mJy$ in $L^\prime$ band after radiative transfer calculations, a value that is in good agreement with HD 100546b observations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stability of Gas Clouds in Galactic Nuclei: An Extended Virial Theorem

Cold gas entering the central $1$ to $10^2$ pc of a galaxy fragments and condenses into clouds. The stability of the clouds determines whether they will be turned into stars or can be delivered to the central supermassive black hole (SMBH) to turn on an active galactic nucleus (AGN). The conventional criteria to assess the stability of these clouds, such as the Jeans criterion and Roche (or tidal) limit, are insufficient here, because they assume the dominance of self-gravity in binding a cloud, and neglect external agents, such as pressure and tidal forces, which are common in galactic nuclei. We formulate a new scheme for judging this stability. We first revisit the conventional Virial theorem, taking into account an external pressure, to identify the correct range of masses that lead to stable clouds. We then extend the theorem to include an external tidal field, crucial for the stability in the region of interest -- in dense star clusters, around SMBHs. We apply our extended Virial theorem to find the correct solutions to practical problems that until now were controversial, namely, the stability of the gas clumps in AGN tori, the circum-nuclear disk in the Galactic Center, and the central molecular zone of the Milky Way. The masses we derive for these structures are orders of magnitude smaller than the commonly-used Virial masses (equivalent to the Jeans mass). Moreover, we prove that these clumps are stable, contrary to what one would naively deduce from the Roche (tidal) limit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Super massive black holes in star forming gaseous circumnuclear discs

Using N-body/SPH simulations we study the evolution of the separation of a pair of SMBHs embedded in a star forming circumnuclear disk (CND). This type of disk is expected to be formed in the central kilo parsec of the remnant of gas-rich galaxy mergers. Our simulations indicate that orbital decay of the SMBHs occurs more quickly when the mean density of the CND is higher, due to increased dynamical friction. However, in simulations where the CND is fragmented in high density gaseous clumps (clumpy CND), the orbits of the SMBHs are erratically perturbed by the gravitational interaction with these clumps, delaying, in some cases, the orbital decay of the SMBHs. The densities of these gaseous clumps in our simulations and in recent studies of clumpy CNDs are significantly higher than the observed density of molecular clouds in isolated galaxies or ULIRGs, thus, we expect that SMBH orbits are perturbed less in real CNDs than in the simulated CNDs of this study and other recent studies. We also find that the migration timescale has a weak dependence on the star formation rate of the CND. Furthermore, the migration timescale of a SMBH pair in a star-forming clumpy CND is at most a factor three longer than the migration timescale of a pair of SMBHs in a CND modeled with more simple gas physics. Therefore, we estimate that the migration timescale of the SMBHs in a clumpy CND is on the order of $10^7$ yrs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The role of feedback in accretion on Low Luminosity AGN: Sgr A* case study

We present numerical models of the gas dynamics in the inner parsec of the Galactic centre. We follow the gas from its origin as stellar winds of several observed young massive stars, until it is either captured by the central black hole, or leaves the system. Unlike our previous models, we include an outflow from the inner accretion flow. Two different kinds of outflows are modelled: (i) an instantaneous-response feedback mode, in which the outflow rate is directly proportional to the current black hole gas capture rate; and (ii) an outburst mode, which is stronger but lasts for a limited time. The latter situation may be particularly relevant to Sgr A*, since there is evidence that Sgr A* was much brighter in the recent past. We find that both types of outflow perturb the gas dynamics near the Bondi radius and the black hole capture rate significantly. The effects persist longer than the outflow itself. We also compare the effects of spherically symmetric and collimated outflows, and find that the latter are far less efficient in transferring its energy to the surrounding gas near the capture radius. Our results imply that accretion feedback is important for non-radiative accretion flows not only within but also outside the capture radius. Steady-state Bondi accretion rate estimates that do not account for feedback outflows over-predict not only the accretion rate onto the black hole but also the capture rate at the Bondi radius itself. Finally, the steady-state assumption under which non-radiative flows have been routinely studied in the literature may have to be abandoned if accretion feedback is bursty in nature.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Tidal disruptions in circumbinary disks. II: Observational signatures in the reverberation spectra

Supermassive Binary Black Holes (SMBBHs) with sub-pc separations form in the course of galaxy mergers, if both galaxies harbour massive black holes. Clear observational evidence for them however still eludes us. We propose a novel method of identifying these systems by means of reverberation mapping their circumbinary disk after a tidal disruption event has ionized it. The tidal disruption of a star at the secondary leads to strong asymmetries in the disk response. We model the shape of the velocity--delay maps for various toy disk models and more realistic gas distributions obtained by SPH simulations. The emissivity of the ionized disk is calculated with {\em Cloudy}. We find peculiar asymmetries in the maps for off center ionizing sources that may help us constrain geometrical parameters of a circumbinary disk such as semi-major axis and orbital phase of the secondary, as well as help strengthen the observational evidence for sub-parsec SMBBHs as such.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the role of initial and boundary conditions in numerical simulations of accretion flows

We study the effects of initial and boundary conditions, taking two-dimensional hydrodynamical numerical simulations of hot accretion flow as an example. The initial conditions considered include a rotating torus, a solution expanded from the one-dimensional global solution of hot accretion flows, injected gas with various angular momentum distributions, and the gas from a large-scale numerical simulation. Special attention is paid to the radial profiles of the mass accretion rate and density. Both can be described by a power-law function, $\dot{M}\propto r^s$ and $ρ\propto r^{-p}$. We find that if the angular momentum is not very low, the value of $s$ is not sensitive to the initial condition and lies within a narrow range, $0.47\la s \la 0.55$. However, the value of $p$ is more sensitive to the initial condition and lies in the range $0.48\la p \la 0.8$. The diversity of the density profile is because different initial conditions give different radial profiles of radial velocity due to the different angular momentum of the initial conditions. When the angular momentum of the accretion flow is very low, the inflow rate is constant with radius. Taking the torus model as an example, we have also investigated the effects of inner and outer boundary conditions by considering the widely adopted "outflow" boundary condition and the "mass flux conservation" condition. We find that the results are not sensitive to these two boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Evolution of binary black holes in self gravitating discs: dissecting the torques

We study the interplay between gas accretion and gravity torques in changing a binary elements and its total angular momentum (L) budget. Especially, we analyse the physical origin of the gravity torques (T_g) and their location within the disc. We analyse 3D SPH simulations of the evolution of initially quasi-circular massive black hole binaries (BHBs) residing in the central hollow of massive self-gravitating circumbinary discs. We use different thermodynamics within the cavity and for the numerical size of the black holes to show that (i) the BHB eccentricity growth found previously is a general result, independent of the accretion and the adopted thermodynamics; (ii) the semi-major axis decay depends both on the T_g and on the interplay with the disc-binary L-transfer due to accretion; (iii) the spectral structure of the T_g is predominately caused by disc edge overdensities and spiral arms developing in the body of the disc and, in general, does not reflect directly the period of the binary; (iv) the net T_g changes sign across the BHB corotation radius. We quantify the relative importance of the two, which appear to depend on the thermodynamical properties of the instreaming gas, and which is crucial in assessing the disc-binary L-transfer; (v) the net torque manifests as a purely kinematic (non-resonant) effect as it stems from the cavity, where the material flows in and out in highly eccentric orbits. Both accretion onto the black holes and the interaction with gas streams inside the cavity must be taken into account to assess the fate of the BHB. Moreover, the total torque exerted by the disc affects L(BHB) by changing all the elements (mass, mass ratio, eccentricity, semimajor axis) of the BHB. Common prescriptions equating tidal torque to semi-major axis shrinking might therefore be poor approximations for real astrophysical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Limiting eccentricity of sub-parsec massive black hole binaries surrounded by self-gravitating gas discs

We study the dynamics of supermassive black hole binaries embedded in circumbinary gaseous discs, with the SPH code Gadget-2. The sub-parsec binary (of total mass M and mass ratio q=1/3) has excavated a gap and transfers its angular momentum to the self--gravitating disc (M_disc=0.2 M). We explore the changes of the binary eccentricity e, by simulating a sequence of binary models that differ in the initial eccentricity e_0, only. In initially low-eccentric binaries, the eccentricity increases with time, while in high-eccentric binaries e declines, indicating the existence of a limiting eccentricity e_crit that is found to fall in the interval [0.6,0.8]. We also present an analytical interpretation for this saturation limit. An important consequence of the existence of e_crit is the detectability of a significant residual eccentricity e_LISA} by the proposed gravitational wave detector LISA. It is found that at the moment of entering the LISA frequency domain e_LISA ~ 10^{-3}-10^{-2}; a signature of its earlier coupling with the massive circumbinary disc. We also observe large periodic inflows across the gap, occurring on the binary and disc dynamical time scales rather than on the viscous time. These periodic changes in the accretion rate (with amplitudes up to ~100%, depending on the binary eccentricity) can be considered a fingerprint of eccentric sub-parsec binaries migrating inside a circumbinary disc.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The butterfly effect in the extreme-mass ratio inspiral problem

Measurements of gravitational waves from the inspiral of a stellar-mass compact object into a massive black hole are unique probes to test General Relativity (GR) and MBH properties, as well as the stellar distribution about these holes in galactic nuclei. Current data analysis techniques can provide us with parameter estimation with very narrow errors. However, an EMRI is not a two-body problem, since other stellar bodies orbiting nearby will influence the capture orbit. Any deviation from the isolated inspiral will induce a small, though observable deviation from the idealised waveform which could be misinterpreted as a failure of GR. Based on conservative analysis of mass segregation in a Milky Way like nucleus, we estimate that the possibility that another star has a semi-major axis comparable to that of the EMRI is non-negligible, although probably very small. This star introduces an observable perturbation in the orbit in the case in which we consider only loss of energy via gravitational radiation. When considering the two first-order non-dissipative post-Newtonian contributions (the periapsis shift of the orbit), the evolution of the orbital elements of the EMRI turns out to be chaotic in nature. The implications of this study are twofold. From the one side, the application to testing GR and measuring MBHs parameters with the detection of EMRIs in galactic nuclei with a millihertz mission will be even more challenging than believed. From the other side, this behaviour could in principle be used as a signature of mass segregation in galactic nuclei.

0 citations0 reviewsNo code linkedOpen access