Showing 1–50 of 90 results for author: Cuadra, J

Multistructured accretion flow of Sgr A* II: Signatures of a Cool Accretion Disk in Hydrodynamic Simulations of Stellar Winds

Authors: Mayura Balakrishnan, Christopher M. P. Russell, Lia Corrales, Diego Calderón, Jorge Cuadra, Daryl Haggard, Sera Markoff, Joey Neilsen, Michael Nowak, Q. Daniel Wang, Fred Baganoff

Abstract: Hydrodynamic simulations of the stellar winds from Wolf-Rayet stars within the Galactic Center can provide predictions for the X-ray spectrum of supermassive black hole Sgr A*. Herein, we present results from updated smooth particle hydrodynamics simulations, building on the architecture of Cuadra et al. (2015); Russell et al. (2017), finding that a cold gas disk forms around Sgr A* with a simulat… ▽ More Hydrodynamic simulations of the stellar winds from Wolf-Rayet stars within the Galactic Center can provide predictions for the X-ray spectrum of supermassive black hole Sgr A*. Herein, we present results from updated smooth particle hydrodynamics simulations, building on the architecture of Cuadra et al. (2015); Russell et al. (2017), finding that a cold gas disk forms around Sgr A* with a simulation runtime of 3500 years. This result is consistent with previous grid-based simulations, demonstrating that a cold disk can form regardless of numerical method. We examine the plasma scenarios arising from an environment with and without this cold disk, by generating synthetic spectra for comparison to the quiescent Fe K alpha Sgr A* spectrum from Chandra HETG-S, taken through the Chandra X-ray Visionary Program. We find that current and future X-ray missions are unlikely to distinguish between the kinematic signatures in the plasma in these two scenarios. Nonetheless, the stellar wind plasma model presents a good fit to the dispersed Chandra spectra within 1.5" of Sgr A*. We compare our results to the Radiatively Inefficient Accretion Flow (RIAF) model fit to the HETG-S spectrum presented in Paper I and find that the Bayesian model evidence does not strongly favor either model. With 9" angular resolution and high spectral resolution of the X-IFU, NewAthena will offer a clearer differentiation between the RIAF plasma model and hydrodynamic simulations, but only a future X-ray mission with arcsecond resolution will significantly advance our understanding of Sgr A*'s accretion flow in X-rays. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 14 pages, 6 figures. Accepted for publication in ApJ

Multistructured accretion flow of Sgr A* I: Examination of a RIAF model

Authors: Mayura Balakrishnan, Lia Corrales, Sera Markoff, Michael Nowak, Daryl Haggard, Q. Daniel Wang, Joey Neilsen, Christopher M. P. Russell, Diego Calderón, Jorge Cuadra, Fred Baganoff

Abstract: The extreme low-luminosity supermassive black hole Sgr A* provides a unique laboratory in which to test radiatively inefficient accretion flow (RIAF) models. Previous fits to the quiescent Chandra ACIS-S spectrum found a RIAF model with an equal inflow-outflow balance works well. In this work, we apply the RIAF model to the Chandra HETG-S spectrum obtained through the Chandra X-ray Visionary Progr… ▽ More The extreme low-luminosity supermassive black hole Sgr A* provides a unique laboratory in which to test radiatively inefficient accretion flow (RIAF) models. Previous fits to the quiescent Chandra ACIS-S spectrum found a RIAF model with an equal inflow-outflow balance works well. In this work, we apply the RIAF model to the Chandra HETG-S spectrum obtained through the Chandra X-ray Visionary Program, which displays features suggestive of temperature and velocity structures within the plasma. A comprehensive forward model analysis accounting for the accretion flow geometry and HETG-S instrumental effects is required for a full interpretation of the quiescent Chandra HETG-S spectrum. We present a RIAF model that takes these effects into account. Our fits to the high-resolution gratings spectrum indicate an inflow balanced by an outflow ($s \sim 1$) alongside a temperature profile that appears shallower than what would be expected from a gravitational potential following $1/r$. The data require that the abundance of Iron relative to solar is $Z_{Fe} < 0.32 Z_\odot$ (90\% credible interval), much lower than the $2~Z_\odot$ metallicity measured in nearby late-type giants. While future missions like NewAthena will provide higher spectral resolution, source separation will continue to be a problem. Leveraging Chandra's unparalleled spatial resolution, which is not expected to be surpassed for decades, remains essential for detailed investigations of the densely populated Galactic Center in X-rays. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 11 pages, 5 figures, 1 table. Accepted for publication in ApJ

Evolution of rotating massive stars adopting a newer, self-consistent wind prescription at SMC metallicity

Authors: Alex Camilo Gormaz-Matamala, Jorge Cuadra, Sylvia Ekström, Georges Meynet, Michel Curé, Krzysztof Belczynski

Abstract: We use Geneva-evolution-code to run evolutionary tracks for stellar masses ranging from $20$ to $85$ $M_\odot$ at SMC metallicity ($Z=0.002$). We upgrade the recipe for stellar winds by adopting our self-consistent m-CAK prescription, which reduces the value of mass-loss rate by a factor between 2 and 6 depending on the mass range. The impact of our new winds is wide, and it can be divided betwe… ▽ More We use Geneva-evolution-code to run evolutionary tracks for stellar masses ranging from $20$ to $85$ $M_\odot$ at SMC metallicity ($Z=0.002$). We upgrade the recipe for stellar winds by adopting our self-consistent m-CAK prescription, which reduces the value of mass-loss rate by a factor between 2 and 6 depending on the mass range. The impact of our new winds is wide, and it can be divided between direct and indirect impact. For the most massive models ($60$ and $85$ $M_\odot$) with $\dot M\gtrsim2\times10^{-7}$ $M_\odot$ yr$^{-1}$, the impact is direct because lower mass loss make stars remove less envelope and therefore remain more massive and less chemically enriched at their surface at the end of their MS phase. For the less massive models ($20$ and $25$ $M_\odot$) with $\dot M\lesssim2\times10^{-8}$ $M_\odot$ yr$^{-1}$, the impact is indirect because lower mass loss make the stars keep high rotational velocities for a longer period of time, then extending the H-core burning lifetime and reaching the end of the MS with higher surface enrichment. Given that the conditions at the H-depletion change, the stars will lose more mass during their He-core burning stages anyways. For $M_\text{zams}=20$ to $40$ $M_\odot$, our models predict stars will evolve through the Hertzsprung gap, from O-type supergiants to BSG and finally RSG, with larger mass fractions of helium compared to old evolution models. New models also set down to $M_\text{zams}=85\,M_\odot$ the minimal initial mass required for a single star to become WR at metallicity $Z=0.002$. New values for $\dot M$ need to be complemented with upgrades in additional features such as convective core overshooting and distribution of rotational velocities, besides more detailed observations from projects such as XShootU, in order to provide a robust framework for the study of massive stars at low metallicity environments. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

Long-term Evolution of Binary Orbits Induced by Circumbinary Disks

Authors: Ruggero Valli, Christopher Tiede, Alejandro Vigna-Gómez, Jorge Cuadra, Magdalena Siwek, Jing-Ze Ma, Daniel J. D'Orazio, Jonathan Zrake, Selma E. de Mink

Abstract: Circumbinary disks are found in a variety of astrophysical scenarios, spanning binary star formation to accreting supermassive black hole binaries. The interaction with a circumbinary disk can yield opposite effects on the binary orbit leading to circularization, or exciting the eccentricity, widening the orbit or shrinking it and facilitating mergers. We present a new formalism for the long-term… ▽ More Circumbinary disks are found in a variety of astrophysical scenarios, spanning binary star formation to accreting supermassive black hole binaries. The interaction with a circumbinary disk can yield opposite effects on the binary orbit leading to circularization, or exciting the eccentricity, widening the orbit or shrinking it and facilitating mergers. We present a new formalism for the long-term evolution of the disk-binary interaction based on the results of recent suites of hydrodynamic simulations, which resolve the complex geometry of the gas in the vicinity of the binary and fully account for the gravitational and accretion forces. We release a python package, \texttt{spindler}, that implements our model. We show that, unless the mass reservoir feeding the disk is comparable to the mass of the binary, accretion onto the binary depletes the disk mass before inducing a significant change in orbital separation or mass ratio. This finding implies that, in most scenarios, interaction with a circumbinary disk is not an efficient mechanism to shrink the orbit of the binary. However, as long as the mass of the disk is at least a few percent of the mass of the binary, the interaction can excite the eccentricity up to an equilibrium value, and induce a statistical correlation between mass ratio and eccentricity. We consider the applicability of our model to a variety of astrophysical scenarios: during star formation, in evolved stellar binaries, triples and in supermassive black hole binaries. We discuss the theoretical and observational implications of our predictions. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: 18 pages, 12 figures. Submitted. Comments welcome. Spindler available as a Python package on Zenodo: https://doi.org/10.5281/zenodo.10529060

Modeling the Galactic center gamma-ray emission with more realistic cosmic-ray dynamics

Authors: Andrés Scherer, Jorge Cuadra, Franz E. Bauer

Abstract: Very-high-energy gamma-ray observations of the Galactic center (GC) show extended emission that is strongly correlated with the morphology of the central molecular zone (CMZ). The best explanation for that emission is a hadronic interaction between cosmic rays (CRs) and ambient gas, where a CR central and continuous source accelerates protons up to 1 PeV ("PeVatron"). However, current models assum… ▽ More Very-high-energy gamma-ray observations of the Galactic center (GC) show extended emission that is strongly correlated with the morphology of the central molecular zone (CMZ). The best explanation for that emission is a hadronic interaction between cosmic rays (CRs) and ambient gas, where a CR central and continuous source accelerates protons up to 1 PeV ("PeVatron"). However, current models assume very simplistic CR dynamics. Our goal is to verify if more realistic CR dynamics for the GC environment are consistent with current gamma-ray observations, and whether they could be constrained by upcoming observations with the Cherenkov Telescope Array (CTA). We generated synthetic gamma-ray maps using a CR transport model with spherical injection, different diffusion regimes (in and out of the CMZ), polar advection, and mono-energetic particles of 1 PeV, and including different CR populations injected from the Arches, Quintuplet, and nuclear clusters of young massive stars, plus supernova Sgr A East. We adopted two different 3D gas distributions consistent with the observed gas column density, either with or without an inner cavity. In order to reproduce the existing observations detected by the High Energy Stereoscopic System (HESS), a ring-like gas distribution, with its mass set by the standard Galactic CO-to-H$_2$ conversion factor, and CR acceleration from all relevant sources are required. For a conversion factor one order of magnitude lower, injection rates that are ten times higher are needed. We show that CTA will be able to differentiate between models with different CR dynamics, proton sources, and CMZ morphology, owing to its unprecedented sensitivity and angular resolution. More realistic CR dynamics suggest that the CMZ has a large inner cavity and that the GC PeVatron is a composite CR population accelerated by the Arches, Quintuplet, and nuclear star clusters, and Sgr A East. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Comments: Accepted for publication in A&A. 15 pages, 11 figures

Journal ref: A&A 679, A114 (2023)

Resolving the Binary Components of the Outbursting Protostar HBC 494 with ALMA

Authors: Pedro Henrique Nogueira, Alice Zurlo, Sebastián Pérez, Camilo González-Ruilova, Lucas A. Cieza, Antonio Hales, Trisha Bhowmik, Dary A. Ruíz-Rodríguez, David A. Principe, Gregory J. Herczeg, Jonathan P. Williams, Jorge Cuadra, Matías Montesinos, Nicolás Cuello, Prachi Chavan, Simon Casassus, Zhaohuan Zhu, Felipe G. Goicovic

Abstract: Episodic accretion is a low-mass pre-main sequence phenomenon characterized by sudden outbursts of enhanced accretion. These objects are classified into two: protostars with elevated levels of accretion that lasts for decades or more, called FUors, and protostars with shorter and repetitive bursts, called EXors. HBC 494 is a FUor object embedded in the Orion Molecular Cloud. Earlier Atacama Large… ▽ More Episodic accretion is a low-mass pre-main sequence phenomenon characterized by sudden outbursts of enhanced accretion. These objects are classified into two: protostars with elevated levels of accretion that lasts for decades or more, called FUors, and protostars with shorter and repetitive bursts, called EXors. HBC 494 is a FUor object embedded in the Orion Molecular Cloud. Earlier Atacama Large (sub-)Millimeter Array (ALMA) continuum observations showed an asymmetry in the disk at 0.''2 resolution. Here, we present follow-up observations at ~0.''03, resolving the system into two components: HBC 494 N (primary) and HBC 494 S (secondary). No circumbinary disk was detected. Both disks are resolved with a projected separation of ~0.''18 (75 au). Their projected dimensions are 84+/-1.8 x 66.9+/-1.5 mas for HBC 494 N and 64.6+/-2.5 x 46.0+/-1.9 mas for HBC 494 S. The disks are almost aligned and with similar inclinations. The observations show that the primary is ~5 times brighter/more massive and ~2 times bigger than the secondary. We notice that the northern component has a similar mass to the FUors, while the southern has to EXors. The HBC 494 disks show individual sizes that are smaller than single eruptive YSOs. In this work, we also report 12CO, 13CO, and C18O molecular line observations. At large scale, the 12CO emission shows bipolar outflows, while the 13CO and C18O maps show a rotating and infalling envelope around the system. At a smaller scale, the 12CO and 13CO moment zero maps show cavities within the continuum disks' area, which may indicate continuum over-subtraction or slow-moving jets and chemical destruction along the line-of-sight. △ Less

Submitted 26 June, 2023; v1 submitted 24 May, 2023; originally announced May 2023.

Comments: Published in MNRAS. Main text (11 pages, 8 figures). Appendix (11 pages, 9 figures). Total (22 pages, 17 figures)

Protoplanet Express, a video game based on numerical simulations

Authors: Jorge Cuadra, Miguel Vergara, Bruno Escárate, Josué Sandoval, Benjamín Cordero

Abstract: Astronomical images can be fascinating to the general public, but the interaction is typically limited to contemplation. Numerical simulations of astronomical systems do permit a closer interaction, but are generally unknown outside the research community. We are developing "Protoplanet Express", a video game based on hydrodynamical simulations of protoplanetary discs. In the game, the player visi… ▽ More Astronomical images can be fascinating to the general public, but the interaction is typically limited to contemplation. Numerical simulations of astronomical systems do permit a closer interaction, but are generally unknown outside the research community. We are developing "Protoplanet Express", a video game based on hydrodynamical simulations of protoplanetary discs. In the game, the player visits several discs, finds its relevant features and learns about them. Here we present the current version of the game, discuss its reception, and consider its further development. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Comments: 6 pages, 4 figures, submitted to the Proceedings of ASP 2022. Beta version of game available at https://questionablegames.itch.io/protoplanet-express

Evolution of rotating massive stars with new hydrodynamic wind models

Authors: Alex Camilo Gormaz-Matamala, Jorge Cuadra, Georges Meynet, Michel Curé

Abstract: Mass loss due to line-driven winds is central to our understanding of the evolution of massive stars. We extend the evolution models introduced in Paper I, where the mass loss recipe is based on the simultaneous calculation of the wind hydrodynamics and the line-acceleration, by incorporating the effects of stellar rotation. We introduce a grid of self-consistent line-force parameters for a set… ▽ More Mass loss due to line-driven winds is central to our understanding of the evolution of massive stars. We extend the evolution models introduced in Paper I, where the mass loss recipe is based on the simultaneous calculation of the wind hydrodynamics and the line-acceleration, by incorporating the effects of stellar rotation. We introduce a grid of self-consistent line-force parameters for a set of standard evolutionary tracks. With that, we generate a new set of evolutionary tracks with rotation for $M_\text{ZAMS}=25,40,70,$ and $120\,M_\odot$, and metallicities $Z=0.014$ and $0.006$. The self-consistent approach gives lower mass loss rates than the standard values adopted in previous evolution models. This decrease impacts strongly on the tracks of the most massive models. Weaker winds allow the star to retain more mass, but also more angular momentum. As a consequence, weaker wind models rotate faster and show a less efficient mixing in their inner stellar structure. The new tracks predict an evolution of the rotational velocities through the MS in close agreement with the range of $\varv\sin i$ values found by recent surveys of Galactic O-type stars. As subsequent implications, the weaker winds from self-consistent models suggest a reduction of the contribution of the isotope $^{26}$Al to the ISM due to stellar winds of massive stars during the MS phase. Moreover, the higher luminosities found for the self-consistent evolutionary models suggest that some populations of massive stars might be less massive than previously thought, as in the case of Ofpe stars at the Galactic Centre. Therefore, this study opens a wide range of consequences for further research based on the evolution of massive stars. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: Accepted for publication in Astronomy & Astrophysics journal

Journal ref: A&A 673, A109 (2023)

Infrared variability of young solar analogs in the Lagoon Nebula

Authors: C. Ordenes-Huanca, M. Zoccali, A. Bayo, J. Cuadra, R. Contreras Ramos, L. A. Hillenbrand, I. Lacerna, S. Abarzua, C. Avendaño, P. Diaz, I. Fernandez, G. Lara

Abstract: T Tauri stars are low-mass pre-main sequence stars that are intrinsically variable. Due to the intense magnetic fields they possess, they develop dark spots on their surface that, because of rotation, introduce a periodic variation of brightness.In addition, the presence of surrounding disks could generate flux variations by variable extinction or accretion. Both can lead to a brightness decrease… ▽ More T Tauri stars are low-mass pre-main sequence stars that are intrinsically variable. Due to the intense magnetic fields they possess, they develop dark spots on their surface that, because of rotation, introduce a periodic variation of brightness.In addition, the presence of surrounding disks could generate flux variations by variable extinction or accretion. Both can lead to a brightness decrease or increase, respectively. Here, we have compiled a catalog of light curves for 379 T Tauri stars in the Lagoon Nebula (M8) region, using VVVX survey data in the Ks-band. All these stars were already classified as pre-MS stars based on other indicators. The data presented here are spread over a period of about eight years, which gives us a unique follow-up time for these sources at this wavelength. The light curves were classified according to their degree of periodicity and asymmetry, to constrain the physical processes responsible for their variation. Periods were compared with the ones found in literature, on a much shorter baseline. This allowed us to prove that for 126 stars, the magnetically active regions remain stable for several years. Besides, our near-IR data were compared with the optical Kepler/K2 light curves, when available, giving us a better understanding of the mechanisms responsible for the brightness variations observed and how they manifest at different bands. We found that the periodicity in both bands is in fairly good agreement, but the asymmetry will depend on the amplitude of the bursts or dips events and the observation cadence. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Comments: 15 pages, 18 figures, accepted for publication in MNRAS

Evolution of massive stars with new hydrodynamic wind models

Authors: A. C. Gormaz-Matamala, M. Curé, G. Meynet, J. Cuadra, J. H. Groh, L. J. Murphy

Abstract: Here we present evolutionary models for a set of massive stars, introducing a new prescription for the mass-loss rate obtained from hydrodynamical calculations in which the wind velocity profile, $v(r)$, and the line-acceleration, $g_\text{line}$, are obtained in a self consistently way. Replacing mass-loss rates at the Main Sequence stage from the standard Vink's formula by our new recipe, we gen… ▽ More Here we present evolutionary models for a set of massive stars, introducing a new prescription for the mass-loss rate obtained from hydrodynamical calculations in which the wind velocity profile, $v(r)$, and the line-acceleration, $g_\text{line}$, are obtained in a self consistently way. Replacing mass-loss rates at the Main Sequence stage from the standard Vink's formula by our new recipe, we generate a new set of evolutionary tracks for $M_\text{ZAMS}=25,40,70$ and $120\,M_\odot$ and metallicities $Z=0.014$ (Galactic), $Z=0.006$ (LMC), and $Z=0.002$ (SMC). Our new derived formula for mass-loss rate predicts a dependence $\dot M\propto Z^a$, where $a$ is not longer constant but dependent on the stellar mass: ranging from $a\sim0.53$ when $M_*\sim120\;M_\odot$, to $a\sim1.02$ when $M_*\sim25\;M_\odot$. We found that models adopting the new recipe for $\dot M$ retain more mass during their evolution, which is expressed in larger radii and consequently more luminous tracks over the Hertzsprung-Russell diagram. These differences are more prominent for the cases of $M_\text{ZAMS}=70$ and 120 $M_\odot$ at solar metallicity, where we found self-consistent tracks are $\sim0.1$ dex brighter and keep extra mass up to 20 $M_\odot$, compared with the classical models using the previous formulation for mass-loss rate. Moreover, we observed remarkable differences for the evolution of the radionuclide isotope $^{26}$Al in the core and the surface of the star. Since $\dot M_\text{sc}$ are weaker than the commonly adopted values for evolutionary tracks, self-consistent tracks predict a later modification in the abundance number of $^{26}$Al in the stellar winds. This new behaviour could provide useful information about the real contribution of this isotope from massive stars to the Galactic interstellar medium. △ Less

Submitted 11 July, 2022; originally announced July 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 665, A133 (2022)

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

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

Authors: Alex C. Gormaz-Matamala, Michel Curé, Alex Lobel, Jorge A. Panei, Jorge Cuadra, Ignacio Araya, Catalina Arcos, Felipe Figueroa-Tapia

Abstract: 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,de… ▽ More 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. △ Less

Submitted 3 March, 2022; v1 submitted 17 February, 2022; originally announced February 2022.

Comments: Accepted for publication in Astronomy & Astrophysics, in 17 February 2022 In this updated version, we corrected information about funding data in Acknowledgements section

Journal ref: A&A 661, A51 (2022)

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

Authors: Andrés Scherer, Jorge Cuadra, Franz E. Bauer

Abstract: 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 cosmi… ▽ More 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. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: 9 pages, 5 figures

Journal ref: A&A 659, A105 (2022)

Cronomoons: origin, dynamics, and light-curve features of ringed exomoons

Authors: Mario Sucerquia, Jaime A. Alvarado-Montes, Amelia Bayo, Jorge Cuadra, Nicolás Cuello, Cristian A. Giuppone, Matías Montesinos, J. Olofsson, Christian Schwab, Lee Spitler, Jorge I. Zuluaga

Abstract: In recent years, technical and theoretical work to detect moons and rings around exoplanets has been attempted. The small mass/size ratios between moons and planets means this is very challenging, having only one exoplanetary system where spotting an exomoon might be feasible (i.e. Kepler-1625b i). In this work, we study the dynamical evolution of ringed exomoons, dubbed "cronomoons" after their s… ▽ More In recent years, technical and theoretical work to detect moons and rings around exoplanets has been attempted. The small mass/size ratios between moons and planets means this is very challenging, having only one exoplanetary system where spotting an exomoon might be feasible (i.e. Kepler-1625b i). In this work, we study the dynamical evolution of ringed exomoons, dubbed "cronomoons" after their similarity with Cronus (Greek for Saturn), and after Chronos (the epitome of time), following the Transit Timing Variations (TTV) and Transit Duration Variation (TDV) that they produce on their host planet. Cronomoons have extended systems of rings that make them appear bigger than they actually are when transiting in front of their host star. We explore different possible scenarios that could lead to the formation of such circumsatellital rings, and through the study of the dynamical/thermodynamic stability and lifespan of their dust and ice ring particles, we found that an isolated cronomoon can survive for time-scales long enough to be detected and followed up. If these objects exist, cronomoons' rings will exhibit gaps similar to Saturn's Cassini Division and analogous to the asteroid belt's Kirkwood gaps, but instead raised due to resonances induced by the host planet. Finally, we analyse the case of Kepler-1625b i under the scope of this work, finding that the controversial giant moon could instead be an Earth-mass cronomoon. From a theoretical perspective, this scenario can contribute to a better interpretation of the underlying phenomenology in current and future observations. △ Less

Submitted 5 December, 2021; originally announced December 2021.

Comments: Accepted to MNRAS. See a Twitter thread at @MarioSucerquia. Press release: http://www.npf.cl/ 14 pages, 11 figures

The protoplanetary disc around HD 169142: circumstellar or circumbinary?

Authors: P. P. Poblete, N. Cuello, S. Pérez, S. Marino, J. Calcino, E. Macías, Á. Ribas, A. Zurlo, J. Cuadra, M. Montesinos, S. Zúñiga-Fernández, A. Bayo, C. Pinte, F. Ménard, D. J. Price

Abstract: Stellar binaries represent a substantial fraction of stellar systems, especially among young stellar objects. Accordingly, binaries play an important role in setting the architecture of a large number of protoplanetary discs. Binaries in coplanar and polar orientations with respect to the circumbinary disc are stable configurations and could induce non-axisymmetric structures in the dust and gas d… ▽ More Stellar binaries represent a substantial fraction of stellar systems, especially among young stellar objects. Accordingly, binaries play an important role in setting the architecture of a large number of protoplanetary discs. Binaries in coplanar and polar orientations with respect to the circumbinary disc are stable configurations and could induce non-axisymmetric structures in the dust and gas distributions. In this work, we suggest that the structures shown in the central region of the protoplanetary disc HD 169142 are produced by the presence of an inner stellar binary and a circumbinary (P-type) planet. We find that a companion with a mass-ratio of 0.1, semi-major axis of 9.9 au, eccentricity of 0.2, and inclination of 90°, together with a 2 Jupiter Mass coplanar planet on a circular orbit at 45 au reproduce the structures at the innermost ring observed at 1.3 mm and the shape of spiral features in scattered light observations. The model predicts changes in the disc's dust structure, and star's astrometric parameters, which would allow testing its veracity by monitoring this system over the next 20 years. △ Less

Submitted 14 December, 2021; v1 submitted 26 November, 2021; originally announced November 2021.

Comments: 11 pages, 7 figures, accepted for publication in MNRAS

The importance of thermal torques on the migration of planets growing by pebble accretion

Authors: O. M. Guilera, M. M. Miller Bertolami, F. Masset, J. Cuadra, J. Venturini, M. P. Ronco

Abstract: A key process in planet formation is the exchange of angular momentum between a growing planet and the protoplanetary disc, which makes the planet migrate through the disc. Several works show that in general low-mass and intermediate-mass planets migrate towards the central star, unless corotation torques become dominant. Recently, a new kind of torque, called the thermal torque, was proposed as a… ▽ More A key process in planet formation is the exchange of angular momentum between a growing planet and the protoplanetary disc, which makes the planet migrate through the disc. Several works show that in general low-mass and intermediate-mass planets migrate towards the central star, unless corotation torques become dominant. Recently, a new kind of torque, called the thermal torque, was proposed as a new source that can generate outward migration of low-mass planets. While the Lindblad and corotation torques depend mostly on the properties of the protoplanetary disc and on the planet mass, the thermal torque depends also on the luminosity of the planet, arising mainly from the accretion of solids. Thus, the accretion of solids plays an important role not only in the formation of the planet but also in its migration process. In a previous work, we evaluated the thermal torque effects on planetary growth and migration mainly in the planetesimal accretion paradigm. In this new work, we study the role of the thermal torque within the pebble accretion paradigm. Computations are carried out consistently in the framework of a global model of planet formation that includes disc evolution, dust growth and evolution, and pebble formation. We also incorporate updated prescriptions of the thermal torque derived from high resolution hydrodynamical simulations. Our simulations show that the thermal torque generates extended regions of outward migration in low viscosity discs. This has a significant impact in the formation of the planets. △ Less

Submitted 10 August, 2021; originally announced August 2021.

Comments: Accepted for publication in MNRAS

Long live the disk: lifetimes of protoplanetary disks in hierarchical triple star systems and a possible explanation for HD 98800 B

Authors: María Paula Ronco, Octavio M. Guilera, Jorge Cuadra, Marcelo M. Miller Bertolami, Nicolás Cuello, Camilo Fontecilla, Pedro Poblete, Amelia Bayo

Abstract: The gas dissipation from a protoplanetary disk is one of the key processes affecting planet formation, and it is widely accepted that it happens on timescales of a few million years for disks around single stars. Over the last years, several protoplanetary disks have been discovered in multiple star systems, and despite the complex environment in which they find themselves, some of them seem to be… ▽ More The gas dissipation from a protoplanetary disk is one of the key processes affecting planet formation, and it is widely accepted that it happens on timescales of a few million years for disks around single stars. Over the last years, several protoplanetary disks have been discovered in multiple star systems, and despite the complex environment in which they find themselves, some of them seem to be quite old, a situation that may favor planet formation. A clear example of this is the disk around HD 98800 B, a binary in a hierarchical quadruple stellar system, which at a $\sim$10 Myr age seems to still be holding significant amounts of gas. Here we present a 1D+1D model to compute the vertical structure and gas evolution of circumbinary disks in hierarchical triple star systems considering different stellar and disk parameters. We show that tidal torques due to the inner binary together with the truncation of the disk due to the external companion strongly reduce the viscous accretion and expansion of the disk. Even allowing viscous accretion by tidal streams, disks in these kind of environments can survive for more than 10 Myr, depending on their properties, with photoevaporation being the main gas dissipation mechanism. We particularly apply our model to the circumbinary disk around HD 98800 B and confirm that its longevity, along with the current non-existence of a disk around the companion binary HD 98800 A, can be explained with our model and by this mechanism. △ Less

Submitted 19 May, 2021; originally announced May 2021.

Comments: 20 pages, 9 figures. Accepted for publication in ApJ

A detailed characterization of HR 8799's debris disk with ALMA in Band 7

Authors: Virginie Faramaz, Sebastian Marino, Mark Booth, Luca Matrà, Eric E. Mamajek, Geoffrey Bryden, Karl R. Stapelfeldt, Simon Casassus, Jorge Cuadra, Antonio S. Hales, Alice Zurlo

Abstract: The exoplanetary system of HR 8799 is one of the rare systems in which multiple planets have been directly imaged. Its architecture is strikingly similar to that of the Solar System, with the four imaged giant planets surrounding a warm dust belt analogous to the Asteroid Belt, and themselves being surrounded by a cold dust belt analogue to the Kuiper Belt. Previous observations of this cold belt… ▽ More The exoplanetary system of HR 8799 is one of the rare systems in which multiple planets have been directly imaged. Its architecture is strikingly similar to that of the Solar System, with the four imaged giant planets surrounding a warm dust belt analogous to the Asteroid Belt, and themselves being surrounded by a cold dust belt analogue to the Kuiper Belt. Previous observations of this cold belt with ALMA in Band 6 (1.3 mm) revealed its inner edge, but analyses of the data differ on its precise location. It was therefore unclear whether the outermost planet HR 8799 b was dynamically sculpting it or not. We present here new ALMA observations of this debris disk in Band 7 (340 GHz, 880 micron). These are the most detailed observations of this disk obtained so far, with a resolution of 1" (40 au) and sensitivity of $9.8\,μ\mathrm{Jy\,beam^{-1}}$, which allowed us to recover the disk structure with high confidence. In order to constrain the disk morphology, we fit its emission using radiative transfer models combined with a MCMC procedure. We find that this disk cannot be adequately represented by a single power law with sharp edges. It exhibits a smoothly rising inner edge and smoothly falling outer edge, with a peak in between, as expected from a disk that contains a high eccentricity component, hence confirming previous findings. Whether this excited population and inner edge shape stem from the presence of an additional planet remains, however, an open question. △ Less

Submitted 5 April, 2021; originally announced April 2021.

Comments: 30 pages, 12 figures, 7 tables. Accepted in AJ

Radiative scale-height and shadows in protoplanetary disks

Authors: Matías Montesinos, Nicolás Cuello, Johan Olofsson, Jorge Cuadra, Amelia Bayo, Gesa H. -M. Bertrang, Clément Perrot

Abstract: Planets form in young circumstellar disks called protoplanetary disks. However, it is still difficult to catch planet formation in-situ. Nevertheless, from recent ALMA/SPHERE data, encouraging evidence of the direct and indirect presence of embedded planets has been identified in disks around young stars: co-moving point sources, gravitational perturbations, rings, cavities, and emission dips or s… ▽ More Planets form in young circumstellar disks called protoplanetary disks. However, it is still difficult to catch planet formation in-situ. Nevertheless, from recent ALMA/SPHERE data, encouraging evidence of the direct and indirect presence of embedded planets has been identified in disks around young stars: co-moving point sources, gravitational perturbations, rings, cavities, and emission dips or shadows cast on disks. The interpretation of these observations needs a robust physical framework to deduce the complex disk geometry. In particular, protoplanetary disk models usually assume the gas pressure scale-height given by the ratio of the sound speed over the azimuthal velocity $H/r = c_{s\rm }/v_{\rm k}$. By doing so, \textit{radiative} pressure fields are often ignored, which could lead to a misinterpretation of the real vertical structure of such disks. We follow the evolution of a gaseous disk with an embedded Jupiter mass planet through hydrodynamical simulations, computing the disk scale-height including radiative pressure, which was derived from a generalization of the stellar atmosphere theory. We focus on the vertical impact of the radiative pressure in the vicinity of circumplanetary disks, where temperatures can reach $\gtrsim 1000$ K for an accreting planet, and radiative forces can overcome gravitational forces from the planet. The radiation-pressure effects create a vertical optically thick column of gas and dust at the proto-planet location, casting a shadow in scattered light. This mechanism could explain the peculiar illumination patterns observed in some disks around young stars such as HD 169142 where a moving shadow has been detected, or the extremely high aspect-ratio $H/r \sim 0.2$ observed in systems like AB Aur and CT Cha. △ Less

Submitted 4 February, 2021; originally announced February 2021.

Comments: 12 pages, 7 figures. Accepted for publication by ApJ

The effect of cooling on the accretion of circumprimary discs inmerging supermassive black hole binaries

Authors: Camilo Fontecilla, Giuseppe Lodato, Jorge Cuadra

Abstract: At the final stages of a supermassive black hole coalescence, the emission of gravitational waves will efficiently remove energy and angular momentum from the binary orbit, allowing the separation between the compact objects to shrink. In the scenario where a circumprimary disc is present, a squeezing phase will develop, in which the tidal interaction between the disc and the secondary black hole… ▽ More At the final stages of a supermassive black hole coalescence, the emission of gravitational waves will efficiently remove energy and angular momentum from the binary orbit, allowing the separation between the compact objects to shrink. In the scenario where a circumprimary disc is present, a squeezing phase will develop, in which the tidal interaction between the disc and the secondary black hole could push the gas inwards, enhancing the accretion rate on to the primary and producing what is known as an electromagnetic precursor. In this context, using 3D hydrodynamic simulations, we study how an adiabatic circumprimary accretion disc responds to the varying gravitational potential as the secondary falls onto the more massive object. We included a cooling prescription controlled by the parameter β= Ωt_{cool}, which will determine how strong the final accretion rate is: a hotter disc is thicker, and the tidal interaction is suppressed for the gas outside the binary plane. Our main results are that for scenarios where the gas cannot cool fast enough (β>30) the disc becomes thick and renders the system invisible, while for β<10 the strong cooling blocks any leakage on to the secondary's orbit, allowing an enhancement in the accretion rate two orders of magnitude stronger than the average through the rest of the simulation. △ Less

Submitted 1 October, 2020; originally announced October 2020.

Comments: to be published in MNRAS

Dust trapping around Lagrangian points in protoplanetary disks

Authors: Matias Montesinos, Juan Garrido-Deutelmoser, Johan Olofsson, Cristian A. Giuppone, Jorge Cuadra, Amelia Bayo, Mario Sucerquia, Nicolas Cuello

Abstract: Trojans are defined as objects that share the orbit of a planet at the stable Lagrangian points $L_4$ and $L_5$. In the Solar System, these bodies show a broad size distribution ranging from micrometer($μ$m) to centimeter(cm) particles (Trojan dust) and up to kilometer (km) rocks (Trojan asteroids). It has also been theorized that earth-like Trojans may be formed in extra-solar systems. The Trojan… ▽ More Trojans are defined as objects that share the orbit of a planet at the stable Lagrangian points $L_4$ and $L_5$. In the Solar System, these bodies show a broad size distribution ranging from micrometer($μ$m) to centimeter(cm) particles (Trojan dust) and up to kilometer (km) rocks (Trojan asteroids). It has also been theorized that earth-like Trojans may be formed in extra-solar systems. The Trojan formation mechanism is still under debate, especially theories involving the effects of dissipative forces from a viscous gaseous environment. We perform hydro-simulations to follow the evolution of a protoplanetary disk with an embedded 1--10 Jupiter-mass planet. On top of the gaseous disk, we set a distribution of $μ$m--cm dust particles interacting with the gas. This allows us to follow dust dynamics as solids get trapped around the Lagrangian points of the planet. We show that large vortices generated at the Lagrangian points are responsible for dust accumulation, where the leading Lagrangian point $L_4$ traps a larger amount of submillimeter (submm) particles than the trailing $L_5$, which traps mostly mm--cm particles. However, the total bulk mass, with typical values of $\sim M_{\rm moon}$, is more significant in $L_5$ than in $L_4$, in contrast to what is observed in the current Solar System a few gigayears later. Furthermore, the migration of the planet does not seem to affect the reported asymmetry between $L_4$ and $L_5$. The main initial mass reservoir for Trojan dust lies in the same co-orbital path of the planet, while dust migrating from the outer region (due to drag) contributes very little to its final mass, imposing strong mass constraints for the in situ formation scenario of Trojan planets. △ Less

Submitted 12 October, 2020; v1 submitted 22 September, 2020; originally announced September 2020.

Comments: 18 pages, 18 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 642, A224 (2020)

How Jupiters save or destroy inner Neptunes around evolved stars

Authors: María Paula Ronco, Matthias R. Schreiber, Cristian A. Giuppone, Dimitri Veras, Jorge Cuadra, Octavio M. Guilera

Abstract: 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,… ▽ More 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. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Comments: Accepted for publication in ApJL. 10 pages, 5 figures

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

Authors: Pedro P. Poblete, Josh Calcino, Nicolás Cuello, Enrique Macías, Álvaro Ribas, Daniel J. Price, Jorge Cuadra, Christophe Pinte

Abstract: 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 circumbina… ▽ More 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. △ Less

Submitted 8 June, 2020; v1 submitted 21 May, 2020; originally announced May 2020.

Comments: Accepted for publication in MNRAS. 11 Pages, 7 Figures

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

Authors: Q. Daniel Wang, Jun Li, Christopher M. P. Russell, Jorge Cuadra

Abstract: 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 opticall… ▽ More 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. △ Less

Submitted 27 December, 2019; originally announced December 2019.

Comments: 17 pages, accepted for publication in MNRAS

Stellar winds pump the heart of the Milky Way

Authors: Diego Calderón, Jorge Cuadra, Marc Schartmann, Andreas Burkert, Christopher M. P. Russell

Abstract: 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 ga… ▽ More 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. △ Less

Submitted 4 December, 2019; v1 submitted 15 October, 2019; originally announced October 2019.

Comments: 6 pages and 5 figures. Accepted for publication in ApJL

Flybys in protoplanetary discs: II. Observational signatures

Authors: Nicolás Cuello, Fabien Louvet, Daniel Mentiplay, Christophe Pinte, Daniel J. Price, Andrew J. Winter, Rebecca Nealon, François Ménard, Giuseppe Lodato, Giovanni Dipierro, Valentin Christiaens, Matías Montesinos, Jorge Cuadra, Guillaume Laibe, Lucas Cieza, Ruobing Dong, Richard Alexander

Abstract: Tidal encounters in star clusters perturb discs around young protostars. In Cuello et al. (2019a, Paper I) we detailed the dynamical signatures of a stellar flyby in both gas and dust. Flybys produce warped discs, spirals with evolving pitch angles, increasing accretion rates, and disc truncation. Here we present the corresponding observational signatures of these features in optical/near-infrared… ▽ More Tidal encounters in star clusters perturb discs around young protostars. In Cuello et al. (2019a, Paper I) we detailed the dynamical signatures of a stellar flyby in both gas and dust. Flybys produce warped discs, spirals with evolving pitch angles, increasing accretion rates, and disc truncation. Here we present the corresponding observational signatures of these features in optical/near-infrared scattered light and (sub-) millimeter continuum and CO line emission. Using representative prograde and retrograde encounters for direct comparison, we post-process hydrodynamical simulations with radiative transfer methods to generate a catalogue of multi-wavelength observations. This provides a reference to identify flybys in recent near-infrared and sub-millimetre observations (e.g., RW Aur, AS 205, HV Tau & DO Tau, FU Ori, V2775 Ori, and Z CMa). △ Less

Submitted 15 October, 2019; originally announced October 2019.

Comments: 11 pages, 5 figures, accepted for publication in MNRAS

From scattered-light to millimeter emission: A comprehensive view of the Gyr-old system of HD 202628 and its eccentric debris ring

Authors: Virginie Faramaz, John Krist, Karl R. Stapelfeldt, Geoffrey Bryden, Eric E. Mamajek, Luca Matrà, Mark Booth, Kevin Flaherty, Antonio S. Hales, A. Meredith Hughes, Amelia Bayo, Simon Casassus, Jorge Cuadra, Johan Olofsson, Kate Y. L. Su, David J. Wilner

Abstract: We present here new observations of the eccentric debris ring surrounding the Gyr-old solar-type star HD 202628: at millimeter wavelengths with ALMA, at far-infrared wavelengths with \textit{Herschel}, and in scattered light with \textit{HST}. The ring inner edge is found to be consistent between ALMA and \textit{HST} data. As radiation pressure affects small grains seen in scattered-light, the ri… ▽ More We present here new observations of the eccentric debris ring surrounding the Gyr-old solar-type star HD 202628: at millimeter wavelengths with ALMA, at far-infrared wavelengths with \textit{Herschel}, and in scattered light with \textit{HST}. The ring inner edge is found to be consistent between ALMA and \textit{HST} data. As radiation pressure affects small grains seen in scattered-light, the ring appears broader at optical than at millimeter wavelengths. The best fit to the ring seen with ALMA has inner and outer edges at $143.1 \pm 1.7$ AU and $165.5 \pm 1.4$, respectively, and an inclination of $57.4^\circ \pm 0.4$ from face-on. The offset of the ring centre of symmetry from the star allows us to quantify its eccentricity to be $e=0.09_{-0.01}^{+0.02}$. This eccentric feature is also detected in low resolution \textit{Herschel}/PACS observations, under the form of a pericenter-glow. Combining the infrared and millimeter photometry, we retrieve a disk grain size distribution index of $\sim -3.4$, and therefore exclude in-situ formation of the inferred belt-shaping perturber, for which we provide new dynamical constraints. Finally, ALMA images show four point-like sources that exceed 100$\,μ$Jy, one of them being just interior to the ring. Although the presence of a background object cannot be excluded, we cannot exclude either that this source is circumplanetary material surrounding the belt-shaper, in which case degeneracies between its mass and orbital parameters could be lifted, allowing us to fully characterize such a distant planet in this mass and age regime for the very first time. △ Less

Submitted 9 September, 2019; originally announced September 2019.

Comments: 29 pages, 15 figures

Dusty clumps in circumbinary discs

Authors: Pedro P. Poblete, Nicolás Cuello, Jorge Cuadra

Abstract: Recent observations have revealed that protoplanetary discs often exhibit cavities and azimuthal asymmetries such as dust traps and clumps. The presence of a stellar binary system in the inner disc regions has been proposed to explain the formation of these structures. Here, we study the dust and gas dynamics in circumbinary discs around eccentric and inclined binaries. This is done through two-fl… ▽ More Recent observations have revealed that protoplanetary discs often exhibit cavities and azimuthal asymmetries such as dust traps and clumps. The presence of a stellar binary system in the inner disc regions has been proposed to explain the formation of these structures. Here, we study the dust and gas dynamics in circumbinary discs around eccentric and inclined binaries. This is done through two-fluid simulations of circumbinary discs, considering different values of the binary eccentricity and inclination. We find that two kinds of dust structures can form in the disc: a single horseshoe-shaped clump, on top of a similar gaseous over-density; or numerous clumps, distributed along the inner disc rim. The latter features form through the complex interplay between the dust particles and the gaseous spirals caused by the binary. All these clumps survive between one and several tens of orbital periods at the feature location. We show that their evolution strongly depends on the gas-dust coupling and the binary parameters. Interestingly, these asymmetric features could in principle be used to infer or constrain the orbital parameters of a stellar companion - potentially unseen - inside the inner disc cavity. Finally, we apply our findings to the disc around AB Aurigae. △ Less

Submitted 20 August, 2019; v1 submitted 15 August, 2019; originally announced August 2019.

Comments: 12 pages, 13 figures, accepted for publication in MNRAS

3D simulations of clump formation in stellar wind collisions

Authors: Diego Calderón, Jorge Cuadra, Marc Schartmann, Andreas Burkert, Joaquín Prieto, Christopher M. P. Russell

Abstract: 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 p… ▽ More 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. △ Less

Submitted 9 January, 2020; v1 submitted 10 June, 2019; originally announced June 2019.

Comments: 21 pages, 20 Figures, 2 Tables. Accepted for publication in MNRAS

Thermal torque effects on the migration of growing low-mass planets

Authors: O. M. Guilera, N. Cuello, M. Montesinos, M. M. Miller Bertolami, M. P. Ronco, J. Cuadra, F. S. Masset

Abstract: As planets grow the exchange of angular momentum with the gaseous component of the protoplanetary disc produces a net torque resulting in a variation of the semi-major axis of the planet. For low-mass planets not able to open a gap in the gaseous disc this regime is known as type I migration. Pioneer works studied this mechanism in isothermal discs finding fast inward type I migration rates that w… ▽ More As planets grow the exchange of angular momentum with the gaseous component of the protoplanetary disc produces a net torque resulting in a variation of the semi-major axis of the planet. For low-mass planets not able to open a gap in the gaseous disc this regime is known as type I migration. Pioneer works studied this mechanism in isothermal discs finding fast inward type I migration rates that were unable to reproduce the observed properties of extrasolar planets. In the last years, several improvements have been made in order to extend the study of type I migration rates to non-isothermal discs. Moreover, it was recently shown that if the planet's luminosity due to solid accretion is taken into account, inward migration could be slowed down and even reversed. In this work, we study the planet formation process incorporating, and comparing, updated type I migration rates for non-isothermal discs and the role of planet's luminosity over such rates. We find that the latter can have important effects on planetary evolution, producing a significant outward migration for the growing planets. △ Less

Submitted 24 April, 2019; originally announced April 2019.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Flybys in protoplanetary discs: I. Gas and dust dynamics

Authors: Nicolás Cuello, Giovanni Dipierro, Daniel Mentiplay, Daniel J. Price, Christophe Pinte, Jorge Cuadra, Guillaume Laibe, François Ménard, Pedro P. Poblete, Matías Montesinos

Abstract: We present 3D smoothed particle hydrodynamics simulations of protoplanetary discs undergoing a flyby by a stellar perturber on a parabolic orbit lying in a plane inclined relative to the disc mid-plane. We model the disc as a mixture of gas and dust, with grains ranging from 1 μm to 10 cm in size. Exploring different orbital inclinations, periastron distances and mass ratios, we investigate the di… ▽ More We present 3D smoothed particle hydrodynamics simulations of protoplanetary discs undergoing a flyby by a stellar perturber on a parabolic orbit lying in a plane inclined relative to the disc mid-plane. We model the disc as a mixture of gas and dust, with grains ranging from 1 μm to 10 cm in size. Exploring different orbital inclinations, periastron distances and mass ratios, we investigate the disc dynamical response during and after the flyby. We find that flybys induce evolving spiral structure in both gas and dust which can persist for thousands of years after periastron. Gas and dust structures induced by the flyby differ because of drag-induced effects on the dust grains. Variations in the accretion rate by up to an order of magnitude occur over a time-scale of order 10 years or less, inducing FU Orionis-like outbursts. The remnant discs are truncated and warped. The dust disc is left more compact than the gas disc, both because of disc truncation and accelerated radial drift of grains induced by the flyby. △ Less

Submitted 3 December, 2018; originally announced December 2018.

Comments: 27 pages, 24 figures, accepted for publication in MNRAS

Dusty spirals triggered by shadows in transition discs

Authors: Nicolás Cuello, Matías Montesinos, Sebastian M. Stammler, Fabien Louvet, Jorge Cuadra

Abstract: Context. Despite the recent discovery of spiral-shaped features in protoplanetary discs in the near-infrared and millimetric wavelengths, there is still an active discussion to understand how they formed. In fact, the spiral waves observed in discs around young stars can be due to different physical mechanisms: planet/companion torques, gravitational perturbations or illumination effects. Aims. We… ▽ More Context. Despite the recent discovery of spiral-shaped features in protoplanetary discs in the near-infrared and millimetric wavelengths, there is still an active discussion to understand how they formed. In fact, the spiral waves observed in discs around young stars can be due to different physical mechanisms: planet/companion torques, gravitational perturbations or illumination effects. Aims. We study the spirals formed in the gaseous phase due to two diametrically opposed shadows cast at fixed disc locations. The shadows are created by an inclined non-precessing disc inside the cavity, which is assumed to be optically thick. In particular, we analyse the effect of these spirals on the dynamics of the dust particles and discuss their detectability in transition discs. Methods. We perform gaseous hydrodynamical simulations with shadows, then we compute the dust evolution on top of the gaseous distribution, and finally we produce synthetic ALMA observations of the dust emission based on radiative transfer calculations. Results. Our main finding is that mm- to cm-sized dust particles are efficiently trapped inside the shadow-triggered spirals. We also observe that particles of various sizes starting at different stellocentric distances are well mixed inside these pressure maxima. This dynamical effect would favour grain growth and affect the resulting composition of planetesimals in the disc. In addition, our radiative transfer calculations show spiral patterns in the disc at 1.6 μm and 1.3 mm. Due to their faint thermal emission (compared to the bright inner regions of the disc) the spirals cannot be detected with ALMA. Our synthetic observations prove however that shadows are observable as dips in the thermal emission. △ Less

Submitted 20 November, 2018; originally announced November 2018.

Comments: 15 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 622, A43 (2019)

F.S.I.: Flyby Scene Investigation

Authors: Nicolás Cuello, Giovanni Dipierro, Daniel Mentiplay, Daniel J. Price, Christophe Pinte, Jorge Cuadra, Guillaume Laibe, François Ménard, Pedro P. Poblete, Matías Montesinos

Abstract: We present a methodology to interpret observations of protoplanetary discs where a flyby, also called a tidal encounter, is suspected. In case of a flyby, protoplanetary discs can be significantly disturbed. The resulting dynamical and kinematical signatures can last for several thousands of years after the flyby and hence deeply affect the evolution of the disc. These effects are stronger for clo… ▽ More We present a methodology to interpret observations of protoplanetary discs where a flyby, also called a tidal encounter, is suspected. In case of a flyby, protoplanetary discs can be significantly disturbed. The resulting dynamical and kinematical signatures can last for several thousands of years after the flyby and hence deeply affect the evolution of the disc. These effects are stronger for closer encounters and more massive perturbers. For the very same flyby parameters, varying the inclination of the perturber's orbit produces a broad range of disc structures: spirals, bridges, warps and cavities. We study this kind of features both in the gas and in the dust for grains ranging from 1 μm to 10 cm in size. Interestingly, the dust exhibits a different dynamical behaviour compared to the gas because of gas-drag effects. Finally, flybys can also trigger high accretion events in the disc-hosting star, readily similar to FU Orionis-type outbursts. All this information can be used to infer the flyby parameters from an incomplete set of observations at different wavelengths. Therefore, the main scope of our flyby scene investigation (FSI) methodology is to help to interpret recent "puzzling" disc observations. △ Less

Submitted 24 October, 2018; originally announced October 2018.

Comments: Proceedings of the SF2A 2018 Conference; 7 pages, 5 figures

Non-steady-state long-term evolution of supermassive black hole binaries surrounded by accretion discs

Authors: Camilo Fontecilla, Zoltán Haiman, Jorge Cuadra

Abstract: Supermassive black holes (SMBHs) pair and form bound binaries after their host galaxies merge. In a gas-rich merger, accretion discs are expected to form around the binary and its components. These discs control the binary orbital evolution until the system is compact enough for gravitational waves to drive the SMBHs to coalescence. In this work, we implemented a time-dependent 1D model to follow… ▽ More Supermassive black holes (SMBHs) pair and form bound binaries after their host galaxies merge. In a gas-rich merger, accretion discs are expected to form around the binary and its components. These discs control the binary orbital evolution until the system is compact enough for gravitational waves to drive the SMBHs to coalescence. In this work, we implemented a time-dependent 1D model to follow the long-term evolution of the coupled binary + discs system, from a separation of 10^5 down to 20 Schwarzschild radii. We run different models changing the system parameters, including the binary mass ratio q \leq 0.3 and a factor γ that controls the inflow across the gap created by the secondary. We find that our implementation yields higher residual masses and longer binary residence times than previous studies. Our main conclusion is the non-steady-state nature of the evolution of the system: the properties the disc had when the binary was still at large separations influence its whole evolution. To recover steady state, the binary residence time would have to be much longer than the inflow time-scale of the disc throughout their entire history, which in general is not satisfied. △ Less

Submitted 5 October, 2018; originally announced October 2018.

Comments: submitted

Is there really a debris disc around $ζ^2\,\mathrm{Reticuli}$ ?

Authors: V. Faramaz, G. Bryden, K. R. Stapelfeldt, M. Booth, A. Bayo, H. Beust, S. Casassus, J. Cuadra, A. Hales, A. M. Hughes, J. Olofsson, K. Y. L. Su, D. J. Wilner

Abstract: The presence of a debris disc around the Gyr-old solar-type star $ζ^2\,\mathrm{Reticuli}$ was suggested by the $\mathit{Spitzer}$ infrared excess detection. Follow-up observations with $\mathit{Herschel}$/PACS revealed a double-lobed feature, that displayed asymmetries both in brightness and position. Therefore, the disc was thought to be edge-on and significantly eccentric. Here we present ALMA/A… ▽ More The presence of a debris disc around the Gyr-old solar-type star $ζ^2\,\mathrm{Reticuli}$ was suggested by the $\mathit{Spitzer}$ infrared excess detection. Follow-up observations with $\mathit{Herschel}$/PACS revealed a double-lobed feature, that displayed asymmetries both in brightness and position. Therefore, the disc was thought to be edge-on and significantly eccentric. Here we present ALMA/ACA observations in Band 6 and 7 which unambiguously reveal that these lobes show no common proper motion with $ζ^2\,\mathrm{Reticuli}$. In these observations, no flux has been detected around $ζ^2\,\mathrm{Reticuli}$ that exceeds the $3σ$ levels. We conclude that surface brightness upper limits of a debris disc around $ζ^2\,\mathrm{Reticuli}$ are $5.7\,\mathrm{μJy/arcsec^2}$ at 1.3 mm, and $26\,\mathrm{μJy/arcsec^2}$ at 870 microns. Our results overall demonstrate the capability of the ALMA/ACA to follow-up $\mathit{Herschel}$ observations of debris discs and clarify the effects of background confusion. △ Less

Submitted 3 September, 2018; originally announced September 2018.

Comments: 6 pages, 2 figures, 2 tables

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 481, Issue 1, 21 November 2018, Pages 44-48

NPF: mirror development in Chile

Authors: Sebastián Zúñiga-Fernández, Amelia Bayo, Johan Olofsson, Leslie Pedrero, Claudio Lobos, Elias Rozas, Nicolás Soto, Matthias Schreiber, Pedro Escárate, Christian Romero, Hayk Hakobyan, Jorge Cuadra, Cristopher Rozas, John D. Monnier, Stefan Kraus, Mike J. Ireland, Pedro Mardones

Abstract: In the era of high-angular resolution astronomical instrumentation, where long and very long baseline interferometers (constituted by many, $\sim$ 20 or more, telescopes) are expected to work not only in the millimeter and submillimeter domain, but also at near and mid infrared wavelengths (experiments such as the Planet Formation Imager, PFI, see Monnier et al. 2018 for an update on its design);… ▽ More In the era of high-angular resolution astronomical instrumentation, where long and very long baseline interferometers (constituted by many, $\sim$ 20 or more, telescopes) are expected to work not only in the millimeter and submillimeter domain, but also at near and mid infrared wavelengths (experiments such as the Planet Formation Imager, PFI, see Monnier et al. 2018 for an update on its design); any promising strategy to alleviate the costs of the individual telescopes involved needs to be explored. In a recent collaboration between engineers, experimental physicists and astronomers in Valparaiso, Chile, we are gaining expertise in the production of light carbon fiber polymer reinforced mirrors. The working principle consists in replicating a glass, or other substrate, mandrel surface with the mirrored adequate curvature, surface characteristics and general shape. Once the carbon fiber base has hardened, previous studies have shown that it can be coated (aluminum) using standard coating processes/techniques designed for glass-based mirrors. The resulting surface quality is highly dependent on the temperature and humidity control among other variables. Current efforts are focused on improving the smoothness of the resulting surfaces to meet near/mid infrared specifications, overcoming, among others, possible deteriorations derived from the replication process. In a second step, at the validation and quality control stage, the mirrors are characterized using simple/traditional tools like spherometers (down to micron precision), but also an optical bench with a Shack-Hartman wavefront sensor. This research line is developed in parallel with a more classical glass-based approach, and in both cases we are prototyping at the small scale of few tens of cms. We here present our progress on these two approaches. △ Less

Submitted 19 July, 2018; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: 13 pages, SPIE Proceedings Volume 10700, Ground-based and Airborne Telescopes VII; 107003X (2018). Event: SPIE Astronomical Telescopes + Instrumentation, 2018, Austin, Texas, United States

The Galactic Centre source G2 was unlikely born in any of the known massive binaries

Authors: Diego Calderón, Jorge Cuadra, Marc Schartmann, Andreas Burkert, Philipp Plewa, Frank Eisenhauer, Maryam Habibi

Abstract: The source G2 has already completed its pericentre passage around Sgr A*, the super-massive black hole in the centre of our Galaxy. Although it has been monitored for 15 years, its astrophysical nature and origin still remain unknown. In this work, we aim to test the hypothesis of G2 being the result of a stellar wind collision. To do so, we study the motion and final fate of gas clumps formed as… ▽ More The source G2 has already completed its pericentre passage around Sgr A*, the super-massive black hole in the centre of our Galaxy. Although it has been monitored for 15 years, its astrophysical nature and origin still remain unknown. In this work, we aim to test the hypothesis of G2 being the result of a stellar wind collision. To do so, we study the motion and final fate of gas clumps formed as a result of collisions of stellar winds in massive binaries. Our approach is based on a test-particle model in order to describe the trajectories of such clumps. The model takes into account the gravitational field of Sgr A*, the interaction of the clumps with the interstellar medium as well as their finite lifetimes. Our analysis allows us to reject the hypothesis based on four arguments: i) if G2 has followed a purely Keplerian orbit since its formation, it cannot have been produced in any of the known massive binaries since their motions are not consistent; ii) in general, gas clumps are evaporated through thermal conduction on very short timescale (< 100yr) before getting close enough to Sgr A*; iii) IRS 16SW, the best candidate for the origin of G2, cannot generate clumps as massive as G2; and iv) clumps ejected from IRS 16SW describe trajectories significantly different to the observed motion of G2. △ Less

Submitted 17 May, 2018; v1 submitted 14 May, 2018; originally announced May 2018.

Comments: 14 pages, 10 figures. Accepted for publication in MNRAS

Circumbinary, not transitional: On the spiral arms, cavity, shadows, fast radial flows, streamers and horseshoe in the HD142527 disc

Authors: Daniel J. Price, Nicolas Cuello, Christophe Pinte, Daniel Mentiplay, Simon Casassus, Valentin Christiaens, Grant M. Kennedy, Jorge Cuadra, Sebastian Perez M., Sebastian Marino, Philip J. Armitage, Alice Zurlo, Attila Juhasz, Enrico Ragusa, Guillaume Laibe, Giuseppe Lodato

Abstract: We present 3D hydrodynamical models of the HD142527 protoplanetary disc, a bright and well studied disc that shows spirals and shadows in scattered light around a 100 au gas cavity, a large horseshoe dust structure in mm continuum emission, together with mysterious fast radial flows and streamers seen in gas kinematics. By considering several possible orbits consistent with the observed arc, we sh… ▽ More We present 3D hydrodynamical models of the HD142527 protoplanetary disc, a bright and well studied disc that shows spirals and shadows in scattered light around a 100 au gas cavity, a large horseshoe dust structure in mm continuum emission, together with mysterious fast radial flows and streamers seen in gas kinematics. By considering several possible orbits consistent with the observed arc, we show that all of the main observational features can be explained by one mechanism - the interaction between the disc and the observed binary companion. We find that the spirals, shadows and horseshoe are only produced in the correct position angles by a companion on an inclined and eccentric orbit approaching periastron - the 'red' family from Lacour et al. (2016). Dust-gas simulations show radial and azimuthal concentration of dust around the cavity, consistent with the observed horseshoe. The success of this model in the HD142527 disc suggests other mm-bright transition discs showing cavities, spirals and dust asymmetries may also be explained by the interaction with central companions. △ Less

Submitted 6 March, 2018; originally announced March 2018.

Comments: 16 pages, 12 figures, accepted to MNRAS. Movies at http://users.monash.edu.au/~dprice/pubs/HD142527/

Accretion of clumpy cold gas onto massive black hole binaries: a possible fast route to binary coalescence

Authors: Felipe G. Goicovic, Cristian Maureira-Fredes, Alberto Sesana, Pau Amaro-Seoane, Jorge Cuadra

Abstract: In currently favoured hierarchical cosmologies, the formation of massive black hole binaries (MBHBs) following galaxy mergers is unavoidable. Still, due the complex physics governing the (hydro)dynamics of the post-merger dense environment of stars and gas in galactic nuclei, the final fate of those MBHBs is still unclear. In gas-rich environments, it is plausible that turbulence and gravitational… ▽ More In currently favoured hierarchical cosmologies, the formation of massive black hole binaries (MBHBs) following galaxy mergers is unavoidable. Still, due the complex physics governing the (hydro)dynamics of the post-merger dense environment of stars and gas in galactic nuclei, the final fate of those MBHBs is still unclear. In gas-rich environments, it is plausible that turbulence and gravitational instabilities feed gas to the nucleus in the form of a series of cold incoherent clumps, thus providing a way to exchange energy and angular momentum between the MBHB and its surroundings. Within this context, we present a suite of smoothed-particle-hydrodynamical models to study the evolution of a sequence of near-radial turbulent gas clouds as they infall towards equal-mass, circular MBHBs. We focus on the dynamical response of the binary orbit to different levels of anisotropy of the incoherent accretion events. Compared to a model extrapolated from a set of individual cloud-MBHB interactions, we find that accretion increases considerably and the binary evolution is faster. This occurs because the continuous infall of clouds drags inwards circumbinary gas left behind by previous accretion events, thus promoting a more effective exchange of angular momentum between the MBHB and the gas. These results suggest that sub-parsec MBHBs efficiently evolve towards coalescence during the interaction with a sequence of individual gas pockets. △ Less

Submitted 27 June, 2018; v1 submitted 15 January, 2018; originally announced January 2018.

Comments: 18 pages, 17 figures. Accepted for publication by MNRAS. Find companion paper at arXiv:1801.06179 Animations available at http://multipleclouds.xyz/movies/

Feedback-limited Accretion: Variable Luminosity from Growing Planets

Authors: Matías Gárate, Jorge Cuadra, Matias Montesinos, Patricia Arévalo

Abstract: 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… ▽ More 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. △ Less

Submitted 14 December, 2020; v1 submitted 3 November, 2017; originally announced November 2017.

Comments: 11 Pages, 10 Figures, accepted to MNRAS. Manuscript updated following the referee report, changing the focus towards the accretion variability

Science with the Cherenkov Telescope Array

Authors: The Cherenkov Telescope Array Consortium, :, B. S. Acharya, I. Agudo, I. Al Samarai, R. Alfaro, J. Alfaro, C. Alispach, R. Alves Batista, J. -P. Amans, E. Amato, G. Ambrosi, E. Antolini, L. A. Antonelli, C. Aramo, M. Araya, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, M. Ashley, M. Backes, C. Balazs, M. Balbo, O. Ballester , et al. (558 additional authors not shown)

Abstract: The Cherenkov Telescope Array, CTA, will be the major global observatory for very high energy gamma-ray astronomy over the next decade and beyond. The scientific potential of CTA is extremely broad: from understanding the role of relativistic cosmic particles to the search for dark matter. CTA is an explorer of the extreme universe, probing environments from the immediate neighbourhood of black ho… ▽ More The Cherenkov Telescope Array, CTA, will be the major global observatory for very high energy gamma-ray astronomy over the next decade and beyond. The scientific potential of CTA is extremely broad: from understanding the role of relativistic cosmic particles to the search for dark matter. CTA is an explorer of the extreme universe, probing environments from the immediate neighbourhood of black holes to cosmic voids on the largest scales. Covering a huge range in photon energy from 20 GeV to 300 TeV, CTA will improve on all aspects of performance with respect to current instruments. The observatory will operate arrays on sites in both hemispheres to provide full sky coverage and will hence maximize the potential for the rarest phenomena such as very nearby supernovae, gamma-ray bursts or gravitational wave transients. With 99 telescopes on the southern site and 19 telescopes on the northern site, flexible operation will be possible, with sub-arrays available for specific tasks. CTA will have important synergies with many of the new generation of major astronomical and astroparticle observatories. Multi-wavelength and multi-messenger approaches combining CTA data with those from other instruments will lead to a deeper understanding of the broad-band non-thermal properties of target sources. The CTA Observatory will be operated as an open, proposal-driven observatory, with all data available on a public archive after a pre-defined proprietary period. Scientists from institutions worldwide have combined together to form the CTA Consortium. This Consortium has prepared a proposal for a Core Programme of highly motivated observations. The programme, encompassing approximately 40% of the available observing time over the first ten years of CTA operation, is made up of individual Key Science Projects (KSPs), which are presented in this document. △ Less

Submitted 21 January, 2018; v1 submitted 22 September, 2017; originally announced September 2017.

Comments: 213 pages, including references and glossary. Version 2: credits and references updated, some figures updated, and author list updated

Cherenkov Telescope Array Contributions to the 35th International Cosmic Ray Conference (ICRC2017)

Authors: F. Acero, B. S. Acharya, V. Acín Portella, C. Adams, I. Agudo, F. Aharonian, I. Al Samarai, A. Alberdi, M. Alcubierre, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, J. -P. Amans, E. Amato, L. Ambrogi, G. Ambrosi, M. Ambrosio, J. Anderson, M. Anduze, E. O. Angüner, E. Antolini, L. A. Antonelli, V. Antonuccio , et al. (1117 additional authors not shown)

Abstract: List of contributions from the Cherenkov Telescope Array Consortium presented at the 35th International Cosmic Ray Conference, July 12-20 2017, Busan, Korea. List of contributions from the Cherenkov Telescope Array Consortium presented at the 35th International Cosmic Ray Conference, July 12-20 2017, Busan, Korea. △ Less

Submitted 24 October, 2017; v1 submitted 11 September, 2017; originally announced September 2017.

Comments: Index of Cherenkov Telescope Array conference proceedings at the ICRC2017, Busan, Korea

Modelling the thermal X-ray emission around the Galactic centre from colliding Wolf-Rayet winds

Authors: Christopher M. P. Russell, Q. Daniel Wang, Jorge Cuadra

Abstract: We compute the thermal X-ray emission from hydrodynamic simulations of the 30 Wolf-Rayet (WR) stars orbiting within a parsec of Sgr A$^*$, with the aim of interpreting the Chandra X-ray observations of this region. The model well reproduces the spectral shape of the observations, indicating that the shocked WR winds are the dominant source of this thermal emission. The model X-ray flux is tied to… ▽ More We compute the thermal X-ray emission from hydrodynamic simulations of the 30 Wolf-Rayet (WR) stars orbiting within a parsec of Sgr A$^*$, with the aim of interpreting the Chandra X-ray observations of this region. The model well reproduces the spectral shape of the observations, indicating that the shocked WR winds are the dominant source of this thermal emission. The model X-ray flux is tied to the strength of the Sgr A$^*$ outflow, which clears out hot gas from the vicinity of Sgr A$^*$. A moderate outflow best fits the present-day observations, even though this supermassive black hole (SMBH) outflow ended $\sim$100 yr ago. △ Less

Submitted 21 July, 2017; originally announced July 2017.

Comments: 1 page, 1 figure; to appear in Proceedings of IAU Symposium 329, "The Lives and Death-Throes of Massive Stars"

Satellite Infall and Mass Deposition on the Galactic Centre

Authors: Sofia G. Gallego, Jorge Cuadra

Abstract: We model the infall of a 2e5 Msun satellite galaxy on to the inner 200 parsec of our Galaxy, to test whether the satellite could perturb the gas previously on stable orbits in the central molecular zone (CMZ), as proposed by Lang et al. (2013). This process would have driven a large gas inflow around 10 Myr ago, necessary to explain the past high accretion rate onto the super-massive black hole, a… ▽ More We model the infall of a 2e5 Msun satellite galaxy on to the inner 200 parsec of our Galaxy, to test whether the satellite could perturb the gas previously on stable orbits in the central molecular zone (CMZ), as proposed by Lang et al. (2013). This process would have driven a large gas inflow around 10 Myr ago, necessary to explain the past high accretion rate onto the super-massive black hole, and the presence of young stars in the inner parsecs of the Galaxy. Our hydrodynamical simulations show a much smaller inflow of gas, not sufficient to produce the aforementioned effects. △ Less

Submitted 16 November, 2016; originally announced November 2016.

Comments: 6 pages, 2 figures, 1 table

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

Authors: Virginie Faramaz, Steve Ertel, Mark Booth, Jorge Cuadra, Charlotte Simmonds

Abstract: 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 system… ▽ More 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. △ Less

Submitted 7 November, 2016; originally announced November 2016.

Comments: 15 pages, 12 figures; Accepted for publication in MNRAS

A Second Decoupling Between Merging Binary Black Holes and the Inner Disc--Impact on the Electromagnetic Counterpart

Authors: Camilo Fontecilla, Xian Chen, Jorge Cuadra

Abstract: The coalescence of two supermassive black holes (SMBHs) produces powerful gravitational-wave (GW) radiation and, if gas is present in the vicinity, also an electromagnetic (EM) counterpart. In the standard picture, an EM outburst will be produced when the binary "decouples" from the circum-binary disc and starts "squeezing" the disc inside the secondary orbit, resulting in its quick accretion on t… ▽ More The coalescence of two supermassive black holes (SMBHs) produces powerful gravitational-wave (GW) radiation and, if gas is present in the vicinity, also an electromagnetic (EM) counterpart. In the standard picture, an EM outburst will be produced when the binary "decouples" from the circum-binary disc and starts "squeezing" the disc inside the secondary orbit, resulting in its quick accretion on to the primary black hole. Here we use analytical arguments and numerical simulations to show that the disc within about $20~R_S$ of a SMBH survives the merger without being depleted. The reason is a "second decoupling": the inner disc thickens due to tidal heating and inefficient cooling, effectively decoupling from the interaction of the binary. We show that this second decoupling quenches the heating sources in the disc ${\cal O}(10^2)$ days before coalescence. This will render the peak UV/X-ray luminosity significantly weaker than previously thought. After the merger, the residual disc cools down and expands, merging with the outer disc rather than being completely accreted. This results in continuous EM emission, hindering the detection of the cut-off and re-brightening proposed in earlier studies. △ Less

Submitted 28 October, 2016; originally announced October 2016.

Comments: Submitted. Comments welcome!

Journal ref: MNRAS, 468, L50 (2017)

Contributions of the Cherenkov Telescope Array (CTA) to the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016)

Authors: The CTA Consortium, :, A. Abchiche, U. Abeysekara, Ó. Abril, F. Acero, B. S. Acharya, C. Adams, G. Agnetta, F. Aharonian, A. Akhperjanian, A. Albert, M. Alcubierre, J. Alfaro, R. Alfaro, A. J. Allafort, R. Aloisio, J. -P. Amans, E. Amato, L. Ambrogi, G. Ambrosi, M. Ambrosio, J. Anderson, M. Anduze, E. O. Angüner , et al. (1387 additional authors not shown)

Abstract: List of contributions from the Cherenkov Telescope Array (CTA) Consortium presented at the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016), July 11-15, 2016, in Heidelberg, Germany. List of contributions from the Cherenkov Telescope Array (CTA) Consortium presented at the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016), July 11-15, 2016, in Heidelberg, Germany. △ Less

Submitted 17 October, 2016; originally announced October 2016.

Comments: Index of CTA conference proceedings for the Gamma 2016, Heidelberg, Germany

Modelling the thermal X-ray emission around the Galactic centre from colliding Wolf-Rayet winds

Authors: Christopher M. P. Russell, Q. Daniel Wang, Jorge Cuadra

Abstract: The Galactic centre is a hotbed of astrophysical activity, with the injection of wind material from $\sim$30 massive Wolf-Rayet (WR) stars orbiting within 12" of the super-massive black hole (SMBH) playing an important role. Hydrodynamic simulations of such colliding and accreting winds produce a complex density and temperature structure of cold wind material shocking with the ambient medium, crea… ▽ More The Galactic centre is a hotbed of astrophysical activity, with the injection of wind material from $\sim$30 massive Wolf-Rayet (WR) stars orbiting within 12" of the super-massive black hole (SMBH) playing an important role. Hydrodynamic simulations of such colliding and accreting winds produce a complex density and temperature structure of cold wind material shocking with the ambient medium, creating a large reservoir of hot, X-ray-emitting gas. This work aims to confront the 3Ms of Chandra X-ray Visionary Program (XVP) observations of this diffuse emission by computing the X-ray emission from these hydrodynamic simulations of the colliding WR winds, amid exploring a variety of SMBH feedback mechanisms. The major success of the model is that it reproduces the spectral shape from the 2"--5" ring around the SMBH, where most of the stellar wind material that is ultimately captured by Sgr A* is shock-heated and thermalised. This naturally explains that the hot gas comes from colliding WR winds, and that the wind speeds of these stars are in general well constrained. The flux level of these spectra, as well as 12"$\times$12" images of 4--9 keV, show the X-ray flux is tied to the SMBH feedback strength; stronger feedback clears out more hot gas, thereby decreasing the thermal X-ray emission. The model in which Sgr A* produced an intermediate-strength outflow during the last few centuries best matches the observations to within about 10%, showing SMBH feedback is required to interpret the X-ray emission in this region. △ Less

Submitted 7 October, 2016; originally announced October 2016.

Comments: 4 pages, 3 figures, to appear in "Proceedings of IAU Symposium 322: The Multi-Messenger Astrophysics of the Galactic Centre," Steve Longmore, Geoff Bicknell and Roland Crocker, eds

Modelling the thermal X-ray emission around the Galactic Centre from colliding Wolf-Rayet winds

Authors: Christopher M. P. Russell, Q. Daniel Wang, Jorge Cuadra

Abstract: The Galactic Centre is a hotbed of astrophysical activity, with the injection of wind material from $\sim$30 massive Wolf-Rayet (WR) stars orbiting within 12 arcsec of the supermassive black hole (SMBH) playing an important role. Hydrodynamic simulations of such colliding and accreting winds produce a complex density and temperature structure of cold wind material shocking with the ambient medium,… ▽ More The Galactic Centre is a hotbed of astrophysical activity, with the injection of wind material from $\sim$30 massive Wolf-Rayet (WR) stars orbiting within 12 arcsec of the supermassive black hole (SMBH) playing an important role. Hydrodynamic simulations of such colliding and accreting winds produce a complex density and temperature structure of cold wind material shocking with the ambient medium, creating a large reservoir of hot, X-ray-emitting gas. This work aims to confront the 3 Ms of Chandra X-ray Visionary Program observations of this diffuse emission by computing the X-ray emission from these hydrodynamic simulations of the colliding WR winds, amid exploring a variety of SMBH feedback mechanisms. The major success of the model is that it reproduces the spectral shape from the 2-5 arcsec ring around the SMBH, where most of the stellar wind material that is ultimately captured by Sgr A* is shock-heated and thermalised. This naturally explains that the hot gas comes from colliding WR winds, and that the wind speeds of these stars are, in general, well constrained. The flux level of these spectra, as well as 12 $\times$ 12-arcsec$^2$ images of 4-9 keV, shows that the X-ray flux is tied to the SMBH feedback strength; stronger feedback clears out more hot gas, thereby decreasing the thermal X-ray emission. The model in which Sgr A* produced an intermediate-strength outflow during the last few centuries best matches the observations to within about 10 per cent, showing that SMBH feedback is required to interpret the X-ray emission in this region. △ Less

Submitted 1 August, 2017; v1 submitted 6 July, 2016; originally announced July 2016.

Comments: 9 pages, 8 figures, published in MNRAS

Journal ref: 2017, MNRAS, 464, 4958-4965

Resolving the Planetesimal Belt of HR 8799 with ALMA

Authors: Mark Booth, Andrés Jordán, Simon Casassus, Antonio S. Hales, William R. F. Dent, Virginie Faramaz, Luca Matrà, Denis Barkats, Rafael Brahm, Jorge Cuadra

Abstract: 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… ▽ More 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. △ Less

Submitted 15 March, 2016; originally announced March 2016.

Comments: 5 pages, accepted by MNRAS Letters

Journal ref: MNRAS Letters, 460, 1, L10-14 (2016)