Showing 1–50 of 120 results for author: Kuiper, R

X-Shooting ULLYSES: Massive Stars at low metallicity -- IV. Spectral analysis methods and exemplary results for O stars

Authors: A. A. C. Sander, J. -C. Bouret, M. Bernini-Peron, J. Puls, F. Backs, S. R. Berlanas, J. M. Bestenlehner, S. A. Brands, A. Herrero, F. Martins, O. Maryeva, D. Pauli, V. Ramachandran, P. A. Crowther, V. M. A. Gómez-González, A. C. Gormaz-Matamala, W. -R. Hamann, D. J. Hillier, R. Kuiper, C. J. K. Larkin, R. R. Lefever, A. Mehner, F. Najarro, L. M. Oskinova, E. C. Schösser , et al. (4 additional authors not shown)

Abstract: CONTEXT: The spectral analysis of hot, massive stars is a fundamental astrophysical method to obtain their intrinsic properties and their feedback. Quantitative spectroscopy for hot, massive stars requires detailed numerical modeling of the atmosphere and an iterative treatment to obtain the best solution within a given framework. AIMS: We present an overview of different techniques for the quanti… ▽ More CONTEXT: The spectral analysis of hot, massive stars is a fundamental astrophysical method to obtain their intrinsic properties and their feedback. Quantitative spectroscopy for hot, massive stars requires detailed numerical modeling of the atmosphere and an iterative treatment to obtain the best solution within a given framework. AIMS: We present an overview of different techniques for the quantitative spectroscopy of hot stars employed within the X-Shooting ULLYSES collaboration, from grid-based approaches to tailored fits. By performing a blind test, we gain an overview about the similarities and differences of the resulting parameters. Our study aims to provide an overview of the parameter spread caused by different approaches. METHODS: For three different stars from the sample (SMC O5 star AzV 377, LMC O7 star Sk -69 50, and LMC O9 star Sk -66 171), we employ different atmosphere codes (CMFGEN, Fastwind, PoWR) and strategies to determine their best-fitting model. For our analyses, UV and optical spectra are used to derive the properties with some methods relying purely on optical data for comparison. To determine the overall spectral energy distribution, we further employ additional photometry from the literature. RESULTS: Effective temperatures for each of three sample stars agree within 3 kK while the differences in log g can be up to 0.2 dex. Luminosity differences of up to 0.1 dex result from different reddening assumptions, which seem to be larger for the methods employing a genetic algorithm. All sample stars are nitrogen-enriched. CONCLUSIONS: We find a reasonable agreement between the different methods. Tailored fitting tends to be able to minimize discrepancies obtained with more course or automatized treatments. UV spectral data is essential for the determination of realistic wind parameters. For one target (Sk -69 50), we find clear indications of an evolved status. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 18+15 pages, 21+4 figures, under review at A&A, condensed abstract

X-Shooting ULLYSES: Massive Stars at low metallicity VIII. Stellar and wind parameters of newly revealed stripped stars in Be binaries

Authors: V. Ramachandran, A. A. C. Sander, D. Pauli, J. Klencki, F. Backs, F. Tramper, M. Bernini-Peron, P. Crowther, W. -R. Hamann, R. Ignace, R. Kuiper, S. Oey, L. M. Oskinova, T. Shenar, H. Todt, J. S. Vink, L. Wang, A. Wofford, the XShootU collaboration

Abstract: On the route towards merging neutron stars and stripped-envelope supernovae, binary population synthesis predicts a large number of post-interaction systems with massive stars that have stripped off their outer layers. Yet, observations of such stars in the intermediate-mass regime below the Wolf-Rayet masses are rare. Using X-Shooting ULLYSES (XShootU) data, we discovered three partially stripped… ▽ More On the route towards merging neutron stars and stripped-envelope supernovae, binary population synthesis predicts a large number of post-interaction systems with massive stars that have stripped off their outer layers. Yet, observations of such stars in the intermediate-mass regime below the Wolf-Rayet masses are rare. Using X-Shooting ULLYSES (XShootU) data, we discovered three partially stripped star + Be/Oe binaries in the Magellanic Clouds. We analyzed the UV and optical spectra using the PoWR model atmosphere code by superimposing model spectra corresponding to each component. The estimated current masses of the partially stripped stars fall within the intermediate mass range of 4-8 $M_{\odot}$. These objects are overluminous for their stellar masses, matching core He-burning luminosities. Their Be/Oe secondaries have much higher masses than their stripped primaries (mass ratio > 2). All three partially stripped stars show significant nitrogen enrichment and carbon and oxygen depletion on their surfaces. Additionally, one of our sample stars exhibits significant helium enrichment. Our study provides the first comprehensive determination of the wind parameters of partially stripped stars in the intermediate mass range. The wind mass-loss rates of these stars are found to be on the order of $10^{-7} M_\odot$ yr$^{-1}$, which is over ten times higher than that of OB stars of the same luminosity. Current evolutionary models characterizing this phase typically employ OB or WR mass-loss rates, which underestimate or overestimate stripped stars' mass-loss rates by an order of magnitude. Binary evolution models indicate that the observed primaries had initial masses of 12-17 $M_{\odot}$, making them potential candidates for stripped-envelope supernovae that form neutron stars. If they survive the explosion, these systems may become Be X-ray binaries and later double neutron stars. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

Radiation hydrodynamical simulations of super-Eddington mass transfer and black hole growth in close binaries

Authors: Daisuke Toyouchi, Kenta Hotokezaka, Kohei Inayoshi, Rolf Kuiper

Abstract: Radiation-driven outflows play a crucial role in extracting mass and angular momentum from binary systems undergoing rapid mass transfer at super-Eddington rates. To study the mass transfer process from a massive donor star to a stellar-mass black hole (BH), we perform multi-dimensional radiation-hydrodynamical simulations that follow accretion flows from the first Lagrange point down to about a h… ▽ More Radiation-driven outflows play a crucial role in extracting mass and angular momentum from binary systems undergoing rapid mass transfer at super-Eddington rates. To study the mass transfer process from a massive donor star to a stellar-mass black hole (BH), we perform multi-dimensional radiation-hydrodynamical simulations that follow accretion flows from the first Lagrange point down to about a hundred times the Schwarzschild radius of the accreting BH. Our simulations reveal that rapid mass transfer occurring at over a thousand times the Eddington rate leads to significant mass loss from the accretion disk via radiation-driven outflows. Consequently, the inflow rates at the innermost radius are regulated by two orders of magnitude smaller than the transfer rates. We find that convective motions within the accretion disk drive outward energy and momentum transport, enhancing the radiation pressure in the outskirts of the disk and ultimately generating large-scale outflows with sufficient energy to leave the binary. Furthermore, we observe strong anisotropy in the outflows, which occur preferentially toward both the closest and furthest points from the donor star. However, when averaged over all directions, the specific angular momentum of the outflows is nearly comparable to the value predicted in the isotropic emission case. Based on our simulation results, we propose a formula that quantifies the mass growth rates on BHs and the mass loss rates from binaries due to radiation-driven outflows. This formula provides important implications for the binary evolution and the formation of merging binary BHs. △ Less

Submitted 12 May, 2024; originally announced May 2024.

Comments: 16 pages, 10 figures, 2 tables, submitted to MNRAS, comments welcome

Mixing is easy: New insights for cosmochemical evolution from pre-stellar core collapse

Authors: Asmita Bhandare, Benoît Commerçon, Guillaume Laibe, Mario Flock, Rolf Kuiper, Thomas Henning, Andrea Mignone, Gabriel-Dominique Marleau

Abstract: Signposts of early planet formation are ubiquitous in substructured young discs. Dense, hot and high-pressure regions formed during gravitational collapse process, integral to star formation, facilitate dynamical mixing of dust within the protostellar disc. This provides an incentive to constrain the role of gas-dust interaction and resolve zones of dust concentration during star-disc formation. W… ▽ More Signposts of early planet formation are ubiquitous in substructured young discs. Dense, hot and high-pressure regions formed during gravitational collapse process, integral to star formation, facilitate dynamical mixing of dust within the protostellar disc. This provides an incentive to constrain the role of gas-dust interaction and resolve zones of dust concentration during star-disc formation. We explore if thermal and dynamical conditions developed during disc formation can generate gas flows that efficiently mix and transport well-coupled gas and dust components. We simulated the collapse of dusty molecular cloud cores with the hydrodynamics code PLUTO augmented with radiation transport and self-gravity. We used a 2D axisymmetric geometry and follow the azimuthal component of velocity. Dust was treated as Lagrangian particles that are subject to drag from the gas, whose motion is computed on a Eulerian grid. We considered 1, 10 and 100 micron-sized neutral spherical dust. Importantly, the equation of state accurately includes molecular hydrogen dissociation. We focus on molecular cloud core masses of 1 and 3 Msun and explore effects of initial rotation rates and cloud core sizes. Our study underlines mechanisms for early transport of dust from inner hot disc regions via the occurrence of meridional flows and outflow. The vortical flow fosters dynamical mixing and retention of dust while thermal pressure driven outflow replenishes dust in the outer disc. Young dynamical precursors to planet-forming discs exhibit regions with complex hydrodynamical gas features and high-temperature structures. These can play a crucial role in concentrating dust for subsequent growth into protoplanets. Dust transport, especially, from sub-au scales surrounding the protostar to outer relatively cooler parts, offers an efficient pathway for thermal reprocessing during pre-stellar core collapse. [Abridged] △ Less

Submitted 12 May, 2024; v1 submitted 14 April, 2024; originally announced April 2024.

Comments: Language edited version

X-Shooting ULLYSES: Massive Stars at low metallicity II. DR1: Advanced optical data products for the Magellanic Clouds

Authors: H. Sana, F. Tramper, M. Abdul-Masih, R. Blomme, K. Dsilva, G. Maravelias, L. Martins, A. Mehner, A. Wofford, G. Banyard, C. L. Barbosa, J. Bestenlehner, C. Hawcroft, D. John Hillier, H. Todt, C. J. K. Larkin, L. Mahy, F. Najarro, V. Ramachandran, M. C. Ramirez-Tannus, M. M. Rubio-Diez, A. A. C. Sander, T. Shenar, J. S. Vink, F. Backs , et al. (12 additional authors not shown)

Abstract: Using the medium resolution spectrograph X-shooter, spectra of 235 OB and Wolf-Rayet (WR) stars in sub-solar metallicity environments have been secured. [...]This second paper focuses on the optical observations of 232 Magellanic Clouds targets. It describes the uniform reduction of the UVB (300 - 560 nm) and VIS (550 - 1020 nm) XShootU data as well as the preparation of advanced data products [..… ▽ More Using the medium resolution spectrograph X-shooter, spectra of 235 OB and Wolf-Rayet (WR) stars in sub-solar metallicity environments have been secured. [...]This second paper focuses on the optical observations of 232 Magellanic Clouds targets. It describes the uniform reduction of the UVB (300 - 560 nm) and VIS (550 - 1020 nm) XShootU data as well as the preparation of advanced data products [...] . The data reduction of the raw data is based on the ESO CPL X-shooter pipeline. We paid particular attention to the determination of the response curves [...] We implemented slit-loss correction, absolute flux calibration, (semi-)automatic rectification to the continuum, and a correction for telluric lines. The spectra of individual epochs were corrected for the barycentric motion, re-sampled and co-added, and the spectra from the two arms were merged into a single flux calibrated spectrum covering the entire optical range with maximum signal-to-noise ratio. [...] We provide three types of data products: (i) two-dimensional spectra for each UVB and VIS exposure; (ii) one-dimensional UVB and VIS spectra before and after response-correction, as well as after applying various processing, including absolute flux calibration, telluric removal, normalisation and barycentric correction; and (iii) co-added flux-calibrated and rectified spectra over the full optical range, for which all available XShootU exposures were combined. For many of the targets, the final signal-to-noise ratio per resolution element is above 200 in both the UVB and the VIS co-added spectra. The reduced data and advanced scientific data products will be made available to the community upon publication of this paper. [...] △ Less

Submitted 26 February, 2024; originally announced February 2024.

Comments: 22 pages, 19 figures; accepted for publication in Astronomy & Astrophysics. Links to online tables and databases will be included upon publication

Observations of high-order multiplicity in a high-mass stellar protocluster

Authors: Shanghuo Li, Patricio Sanhueza, Henrik Beuther, Huei-Ru Vivien Chen, Rolf Kuiper, Fernando A. Olguin, Ralph E. Pudritz, Ian W. Stephens, Qizhou Zhang, Fumitaka Nakamura, Xing Lu, Rajika L. Kuruwita, Takeshi Sakai, Thomas Henning, Kotomi Taniguchi, Fei Li

Abstract: The dominant mechanism forming multiple stellar systems in the high-mass regime (M$_\ast \gtrsim $ 8 $M_{\odot}$) remained unknown because direct imaging of multiple protostellar systems at early phases of high-mass star formation is very challenging. High-mass stars are expected to form in clustered environments containing binaries and higher-order multiplicity systems. So far only a few high-mas… ▽ More The dominant mechanism forming multiple stellar systems in the high-mass regime (M$_\ast \gtrsim $ 8 $M_{\odot}$) remained unknown because direct imaging of multiple protostellar systems at early phases of high-mass star formation is very challenging. High-mass stars are expected to form in clustered environments containing binaries and higher-order multiplicity systems. So far only a few high-mass protobinary systems, and no definitive higher-order multiples, have been detected. Here we report the discovery of one quintuple, one quadruple, one triple and four binary protostellar systems simultaneously forming in a single high-mass protocluster, G333.23--0.06, using Atacama Large Millimeter/submillimeter Array high-resolution observations. We present a new example of a group of gravitationally bound binary and higher-order multiples during their early formation phases in a protocluster. This provides the clearest direct measurement of the initial configuration of primordial high-order multiple systems, with implications for the in situ multiplicity and its origin. We find that the binary and higher-order multiple systems, and their parent cores, show no obvious sign of disk-like kinematic structure. We conclude that the observed fragmentation into binary and higher-order multiple systems can be explained by core fragmentation, indicating its crucial role in establishing the multiplicity during high-mass star cluster formation. △ Less

Submitted 12 January, 2024; originally announced January 2024.

Comments: Published in Nature Astronomy at https://www.nature.com/articles/s41550-023-02181-9. 27 pages, 12 figures, 1 table

A probable Keplerian disk feeding an optically revealed massive young star

Authors: Anna F. McLeod, Pamela D. Klaassen, Megan Reiter, Jonathan Henshaw, Rolf Kuiper, Adam Ginsburg

Abstract: The canonical picture of star formation involves disk-mediated accretion, with Keplerian accretion disks and associated bipolar jets primarily observed in nearby, low-mass young stellar objects (YSOs). Recently, rotating gaseous structures and Keplerian disks have been detected around a number of massive (M > 8 solar masses) YSOs (MYSOs) including several disk-jet systems. All of the known MYSO sy… ▽ More The canonical picture of star formation involves disk-mediated accretion, with Keplerian accretion disks and associated bipolar jets primarily observed in nearby, low-mass young stellar objects (YSOs). Recently, rotating gaseous structures and Keplerian disks have been detected around a number of massive (M > 8 solar masses) YSOs (MYSOs) including several disk-jet systems. All of the known MYSO systems are located in the Milky Way, and all are embedded in their natal material. Here we report the detection of a rotating gaseous structure around an extragalactic MYSO in the Large Magellanic Cloud. The gas motions show radial flow of material falling from larger scales onto a central disk-like structure, the latter exhibiting signs of Keplerian rotation, i.e., a rotating toroid feeding an accretion disk and thus the growth of the central star. The system is in almost all aspects comparable to Milky Way high-mass young stellar objects accreting gas via a Keplerian disk. The key difference between this source and its Galactic counterparts is that it is optically revealed, rather than being deeply embedded in its natal material as is expected of such a young massive star. We suggest that this is the consequence of the star having formed in a low-metallicity and low-dust content environment, thus providing important constraints for models of the formation and evolution of massive stars and their circumstellar disks. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: 20 pages, 9 pages

Density distributions, magnetic field structures and fragmentation in high-mass star formation

Authors: H. Beuther, C. Gieser, J. D. Soler, Q. Zhang, R. Rao, D. Semenov, Th. Henning, R. Pudritz, T. Peters, P. Klaassen, M. T. Beltran, A. Palau, T. Moeller, K. G. Johnston, H. Zinnecker, J. Urquhart, R. Kuiper, A. Ahmadi, A. Sanchez-Monge, S. Feng, S. Leurini, S. E. Ragan

Abstract: Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures. Results… ▽ More Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures. Results: We infer density distributions n~r^{-p} of the regions with typical power-law slopes p around ~1.5. There is no obvious correlation between the power-law slopes of the density structures on larger clump scales (~1pc) and the number of fragments on smaller core scales (<0.1pc). Comparing the large-scale single-dish density profiles to those derived earlier from interferometric observations at smaller spatial scales, we find that the smaller-scale power-law slopes are steeper, typically around ~2.0. The flattening toward larger scales is consistent with the star-forming regions being embedded in larger cloud structures that do not decrease in density away from a particular core. Regarding the magnetic field, for several regions it appears aligned with filamentary structures leading toward the densest central cores. Furthermore, we find different polarization structures with some regions exhibiting central polarization holes whereas other regions show polarized emission also toward the central peak positions. Nevertheless, the polarized intensities are inversely related to the Stokes I intensities. We estimate magnetic field strengths between ~0.2 and ~4.5mG, and we find no clear correlation between magnetic field strength and the fragmentation level of the regions. Comparison of the turbulent to magnetic energies shows that they are of roughly equal importance in this sample. The mass-to-flux ratios range between ~2 and ~7, consistent with collapsing star-forming regions. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: Accepted for Astronomy & Astrophysics, 14 pages, 14 figures plus appendices, also download option at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

The magnetic field of a magnetohydrodynamic disk wind: Water maser observations and simulations

Authors: Luca Moscadelli, André Oliva, Gabriele Surcis, Alberto Sanna, Maria Teresa Beltrán, Rolf Kuiper

Abstract: Our goal is to measure and model the magnetic field distribution in the disk wind of the young stellar object (YSO) IRAS 21078+5211. We performed sensitive global Very Long Baseline Interferometry observations of the polarized emission of the 22 GHz water masers tracing individual streamlines of the magnetohydrodynamic (MHD) disk wind in IRAS 21078+5211. Our resistive-radiative-gravito-MHD simulat… ▽ More Our goal is to measure and model the magnetic field distribution in the disk wind of the young stellar object (YSO) IRAS 21078+5211. We performed sensitive global Very Long Baseline Interferometry observations of the polarized emission of the 22 GHz water masers tracing individual streamlines of the magnetohydrodynamic (MHD) disk wind in IRAS 21078+5211. Our resistive-radiative-gravito-MHD simulations of a jet around a forming massive star are able to reproduce the maser kinematics in the inner jet cavity. We measure a weak level of 0.3%-3.2% of linear and circular polarization in 24 and 8 water masers, respectively. The detected polarized masers sample the direction and the strength of the magnetic field along five distinct streamlines within the inner 100 au region of the disk wind. Along the four streamlines at smaller radii from the jet axis (< 25 au), the sky-projected direction of the magnetic field forms, in most cases, a small offset angle of < 30$^{\circ}$ with the tangent to the streamline. Along the stream at larger radii (50-100 au), the magnetic field is sampled at only three separated positions, and it is found to be approximately perpendicular to the streamline tangent at heights of approximately 10 and 40 au, and parallel to the tangent at about 70 au. According to our simulations, the magnetic field lines should coincide with the flow streamlines in the inner jet cavity. The small tilt in the magnetic field direction observed along the inner streams can be well explained by Faraday rotation, assuming a realistic low level of ionization for the molecular shell of the jet of namely 10$^{-2}$. The magnetic field amplitudes measured from maser circular polarization are all within a relatively small range of 100-700 mG, which is in good agreement with the simulation results and consistent with reduced magnetic diffusivity in the jet cavity owing to efficient shock ionization. △ Less

Submitted 21 October, 2023; originally announced October 2023.

Comments: 14 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 680, A107 (2023)

Atmospheric Recyling of Volatiles by Pebble-Accreting Planets

Authors: Yu Wang, Chris W. Ormel, Pinghui Huang, Rolf Kuiper

Abstract: Planets, embedded in their natal discs, harbour hot envelopes. When pebbles are accreted by these planets, the contained volatile components may sublimate, enriching the envelope and potentially changing its thermodynamical properties. However, the envelopes of embedded planets actively exchange material with the disc, which would limit the buildup of a vapour-rich atmosphere. To properly investig… ▽ More Planets, embedded in their natal discs, harbour hot envelopes. When pebbles are accreted by these planets, the contained volatile components may sublimate, enriching the envelope and potentially changing its thermodynamical properties. However, the envelopes of embedded planets actively exchange material with the disc, which would limit the buildup of a vapour-rich atmosphere. To properly investigate these processes, we have developed a new phase change module to treat the sublimation process with hydrodynamical simultions. Combined with the recently developed multi-dust fluid approach, we conduct 2D self-consistent hydrodynamic simulations to study how pebble sublimation influences the water content of super-Earths and sub-Neptunes. We find the extent and the amount of vapour that a planet is able to hold on to is determined by the relative size of the sublimation front and the atmosphere. When the sublimation front lies far inside the atmosphere, vapour tends to be locked deep in the atmosphere and keeps accumulating through a positive feedback mechanism. On the other hand, when the sublimation front exceeds the (bound) atmosphere, the ice component of incoming pebbles can be fully recycled and the vapour content reaches a low, steady value. Low disc temperature, small planet mass and high pebble flux (omitting accretion heating by pebbles) render the planet atmosphere vapour-rich while the reverse changes render it vapour-poor. The phase change module introduced here can in future studies also be employed to model the chemical composition of the gas in the vicinity of accreting planets and around snowlines. △ Less

Submitted 9 June, 2023; originally announced June 2023.

Comments: 21 pages, 16 figures. Accepted for publication in MNRAS on June 5th 2023

X-Shooting ULLYSES: massive stars at low metallicity. I. Project Description

Authors: Jorick S. Vink, A. Mehner, P. A. Crowther, A. Fullerton, M. Garcia, F. Martins, N. Morrell, L. M. Oskinova, N. St-Louis, A. ud-Doula, A. A. C. Sander, H. Sana, J. -C. Bouret, B. Kubatova, P. Marchant, L. P. Martins, A. Wofford, J. Th. van Loon, O. Grace Telford, Y. Gotberg, D. M. Bowman, C. Erba, V. M. Kalari, M. Abdul-Masih, T. Alkousa , et al. (56 additional authors not shown)

Abstract: Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity… ▽ More Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity (Z). The Hubble Space Telescope has devoted 500 orbits to observe 250 massive stars at low Z in the ultraviolet (UV) with the COS and STIS spectrographs under the ULLYSES program. The complementary ``X-Shooting ULLYSES'' (XShootU) project provides enhanced legacy value with high-quality optical and near-infrared spectra obtained with the wide-wavelength coverage X-shooter spectrograph at ESO's Very Large Telescope. We present an overview of the XShootU project, showing that combining ULLYSES UV and XShootU optical spectra is critical for the uniform determination of stellar parameters such as effective temperature, surface gravity, luminosity, and abundances, as well as wind properties such as mass-loss rates in function of Z. As uncertainties in stellar and wind parameters percolate into many adjacent areas of Astrophysics, the data and modelling of the XShootU project is expected to be a game-changer for our physical understanding of massive stars at low Z. To be able to confidently interpret James Webb Space Telescope (JWST) spectra of the first stellar generations, the individual spectra of low Z stars need to be understood, which is exactly where XShootU can deliver. △ Less

Submitted 1 June, 2023; v1 submitted 10 May, 2023; originally announced May 2023.

Comments: Accepted in A&A - 35 Pages, 12 Figures, 4 Tables, 2 Large Tables

Journal ref: A&A 675, A154 (2023)

The Planetary Accretion Shock. III. Smoothing-free 2.5D simulations and calculation of H alpha emission

Authors: Gabriel-Dominique Marleau, Rolf Kuiper, William Béthune, Christoph Mordasini

Abstract: Surveys have looked for H alpha emission from accreting gas giants but found very few objects. Analyses of the detections and non-detections have assumed that the entire gas flow feeding the planet is in radial free-fall. However, hydrodynamical simulations suggest that this is far from reality. We calculate the H alpha emission from multidimensional accretion onto a gas giant, following the gas f… ▽ More Surveys have looked for H alpha emission from accreting gas giants but found very few objects. Analyses of the detections and non-detections have assumed that the entire gas flow feeding the planet is in radial free-fall. However, hydrodynamical simulations suggest that this is far from reality. We calculate the H alpha emission from multidimensional accretion onto a gas giant, following the gas flow from Hill-sphere scales down to the circumplanetary disc (CPD) and the planetary surface. We perform azimuthally-symmetric radiation-hydrodynamics simulations around the planet and use modern tabulated gas and dust opacities. Crucially, contrasting with most previous simulations, we do not smooth the gravitational potential and do follow the flow down to the planetary surface, where grid cells are 0.01 Jupiter radii small radially. We find that only roughly one percent of the net gas inflow into the Hill sphere reaches directly the planet. As expected for ballistic infall trajectories, most of the gas falls at too large a distance on the CPD to generate H alpha. Including radiation transport removes the high-velocity sub-surface flow previously seen in hydrodynamics-only simulations, so that only the free planet surface and the inner regions of the CPD emit substantially H alpha. Unless magnetospheric accretion, which we neglect here, additionally produces H alpha, the corresponding H alpha production efficiency is much smaller than usually assumed, which needs to be taken into account when analysing (non-)detection statistics. △ Less

Submitted 19 July, 2023; v1 submitted 2 May, 2023; originally announced May 2023.

Comments: Main text: 16 pages, 7 figures; appendices: 4 pages, 4 figures; references: 1.7 pages. Printed at ApJ on 2023 July 20th. Version 4: This information and correcting R_Hill -> R_Jup (obvious) in the axis labels of Fig. 9

Journal ref: The Astrophysical Journal, Volume 952, Issue 1, article id. 89, 18 pp. (2023)

Kinematics and stability of high-mass protostellar disk candidates at sub-arcsecond resolution -- Insights from the IRAM NOEMA large program CORE

Authors: Aida Ahmadi, H. Beuther, F. Bosco, C. Gieser, S. Suri, J. C. Mottram, R. Kuiper, Th. Henning, Á. Sánchez-Monge, H. Linz, R. E. Pudritz, D. Semenov, J. M. Winters, T. Möller, M. T. Beltrán, T. Csengeri, R. Galván-Madrid, K. G. Johnston, E. Keto, P. D. Klaassen, S. Leurini, S. N. Longmore, S. L. Lumsden, L. T. Maud, L. Moscadelli , et al. (6 additional authors not shown)

Abstract: The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in… ▽ More The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars. △ Less

Submitted 3 May, 2023; v1 submitted 28 April, 2023; originally announced May 2023.

Comments: 27 pages, 12 figures, 6 appendices - accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 677, A171 (2023)

Spectroscopic and evolutionary analyses of the binary system AzV 14 outline paths toward the WR stage at low metallicity

Authors: D. Pauli, L. M. Oskinova, W. -R. Hamann, D. M. Bowman, H. Todt, T. Shenar, A. A. C. Sander, C. Erba, V. M. A. Gómez-González, C. Kehrig, J. Klencki, R. Kuiper, A. Mehner, S. E. de Mink, M. S. Oey, V. Ramachandran, A. Schootemeijer, S. Reyero Serantes, A. Wofford

Abstract: The origin of the observed population of Wolf-Rayet (WR) stars in low-metallicity (low-Z) galaxies, such as the Small Magellanic Cloud (SMC), is not yet understood. Standard, single-star evolutionary models predict that WR stars should stem from very massive O-type star progenitors, but these are very rare. On the other hand, binary evolutionary models predict that WR stars could originate from pr… ▽ More The origin of the observed population of Wolf-Rayet (WR) stars in low-metallicity (low-Z) galaxies, such as the Small Magellanic Cloud (SMC), is not yet understood. Standard, single-star evolutionary models predict that WR stars should stem from very massive O-type star progenitors, but these are very rare. On the other hand, binary evolutionary models predict that WR stars could originate from primary stars in close binaries. We conduct an analysis of the massive O star, AzV 14, to spectroscopically determine its fundamental and stellar wind parameters, which are then used to investigate evolutionary paths from the O-type to the WR stage with stellar evolutionary models. Multi-epoch UV and optical spectra of AzV 14 are analyzed using the non-LTE stellar atmosphere code PoWR. An optical TESS light curve was extracted and analyzed using the PHOEBE code. The obtained parameters are put into an evolutionary context, using the MESA code. AzV 14 is a close binary system consisting of two similar main sequence stars with masses of 32 Msol. Both stars have weak stellar winds with mass-loss rates of log $\dot{M}$ = -7.7. Binary evolutionary models can explain the empirically derived stellar and orbital parameters. The model predicts that the primary will evolve into a WR star with T = 100 kK, while the secondary, which will accrete significant amounts of mass during the first mass transfer phase, will become a cooler WR star with T = 50 kK and are predicted to have compared to other WR stars increased oxygen abundances. This model prediction is supported by a spectroscopic analysis of a WR star in the SMC. We hypothesize that the populations of WR stars in low-Z galaxies may have bimodal temperature distributions. Hotter WR stars might originate from primary stars, while cooler WR stars are the evolutionary descendants of the secondary stars if they accreted a significant amount of mass. △ Less

Submitted 7 March, 2023; originally announced March 2023.

Comments: 21 pages (13 main body + 8 appendix), 16 figures, 9 tables

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

Disk fragmentation around a massive protostar: a comparison of two three-dimensional codes

Authors: Raphaël Mignon-Risse, André Oliva, Matthias González, Rolf Kuiper, Benoît Commerçon

Abstract: (Abridged) Most massive stars are located in multiple systems. The modeling of disk fragmentation, a possible mechanism leading to stellar multiplicity, relies on parallel 3D simulation codes whose agreement remains to be evaluated. Using the Cartesian AMR code RAMSES, we compare disk fragmentation in a centrally-condensed protostellar system to the study of Oliva & Kuiper (2020) performed on a gr… ▽ More (Abridged) Most massive stars are located in multiple systems. The modeling of disk fragmentation, a possible mechanism leading to stellar multiplicity, relies on parallel 3D simulation codes whose agreement remains to be evaluated. Using the Cartesian AMR code RAMSES, we compare disk fragmentation in a centrally-condensed protostellar system to the study of Oliva & Kuiper (2020) performed on a grid in spherical coordinates using PLUTO. Two RAMSES runs are considered and give qualitatively distinct pictures. When allowing for unlimited sink particle creation, gas fragmentation leads to a multiple stellar system whose multiplicity is affected by the grid when triggering fragmentation and by numerically-assisted mergers. On the other hand, using a unique, central, fixed sink particle, a centrally-condensed system forms, similar to that reported in PLUTO. The RAMSES-PLUTO comparison is performed with the latter: agreement between the two codes is found regarding the first rotationally-supported disk formation, the presence of an accretion shock onto it, the first fragmentation phase. Gaseous fragments form and their properties are in agreement between the two codes. As a minor difference, fragments dynamics causes the disk structure to be sub-Keplerian in RAMSES whereas it is found to be Keplerian and reaches quiescence in PLUTO. We attribute this discrepancy to the central star being twice less massive in RAMSES because of the different stellar accretion subgrid models. In a centrally-condensed system, the agreement between RAMSES and PLUTO regarding many of the collapse properties and fragmentation process is good. Fragmentation occurring in the innermost region and numerical choices (use of sink particles, grid) have a crucial impact when similar but smooth initial conditions are employed - more crucial than the code's choice - on the system's outcome, multiple or centrally-condensed. △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: 17 pages, 16 figures; accepted for publication in A&A

Journal ref: A&A 672, A88 (2023)

Modeling disks and magnetic outflows around a forming massive star: I. Investigating the two layer-structure of the accretion disk

Authors: André Oliva, Rolf Kuiper

Abstract: Like their lower mass siblings, massive protostars can be expected to: a) be surrounded by circumstellar disks and b) launch magnetically-driven jets and outflows. The disk formation and global evolution is thereby controlled by advection of angular momentum from large scales, the efficiency of magnetic braking and the resistivity of the medium, and the internal thermal and magnetic pressures of t… ▽ More Like their lower mass siblings, massive protostars can be expected to: a) be surrounded by circumstellar disks and b) launch magnetically-driven jets and outflows. The disk formation and global evolution is thereby controlled by advection of angular momentum from large scales, the efficiency of magnetic braking and the resistivity of the medium, and the internal thermal and magnetic pressures of the disk. We perform a series of 30 simulations of a massive star forming from the gravitational collapse of a molecular cloud threaded by an initially-uniform magnetic field, starting from different values for the mass of the cloud, its initial density and rotation profiles, its rotational energy content, the magnetic field strength, and the resistivity of the material. The gas and dust is modeled with the methods of resistive magnetohydrodynamics, also considering radiation transport of thermal emission and self-gravity. After the initial infall phase dominated by the gravitational collapse, an accretion disk is formed, shortly followed by the launching of magnetically-driven outflows. Two layers can be distinguished in the accretion disk: a thin layer, vertically supported by thermal pressure, and a thick layer, vertically supported by magnetic pressure. We observe the effects of magnetic braking in the inner ~50 au of the disk at late times in our fiducial case. The parameter study reveals that the size of the disk is mostly determined by the density and rotation profiles of the initial mass reservoir and not by the magnetic field strength. Magnetic pressure can slightly increase the size of the accretion disk, while magnetic braking is more relevant in the innermost parts of the disk as opposed to the outer disk. From the parameter study, we infer that multiple initial conditions for the onset of gravitational collapse are able to produce a given disk size and protostellar mass. △ Less

Submitted 31 October, 2022; originally announced October 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 669, A80 (2023)

Modeling disks and magnetic outflows around a forming massive star: II. Dynamics of jets from massive protostars

Authors: André Oliva, Rolf Kuiper

Abstract: Forming massive stars launch outflows of magnetic origin, which in fact serve as a marker for finding sites of massive star formation. However, both the theoretical and observational study of the mechanisms that intervene in the formation and propagation of such outflows has been possible only until recent years. With this work, we aim to study in detail the mechanisms that drive highly collimated… ▽ More Forming massive stars launch outflows of magnetic origin, which in fact serve as a marker for finding sites of massive star formation. However, both the theoretical and observational study of the mechanisms that intervene in the formation and propagation of such outflows has been possible only until recent years. With this work, we aim to study in detail the mechanisms that drive highly collimated outflows from early stages of the formation of a massive star, and how those processes are impacted by the properties of the natal environment of the forming massive star. We perform a series of 31 simulations with the aim of building a unified theoretical picture of these mechanisms, and see how the impact of different environments alter their morphology and momentum output. The magnetohydrodynamical simulations consider also Ohmic dissipation as a nonideal effect, self-gravity, and diffusive radiation transport for thermal absorption and emission by the dust and gas. We start from a collapsing cloud core that is threaded by an initially-uniform magnetic field and which is slowly rotating. We utilize a two-dimensional axisymmetric grid in spherical coordinates. In the simulations, we can clearly distinguish a fast, magneto-centrifugally launched and collimated jet (of speeds > 100 km/s), from a wider magnetic tower flow driven by magnetic pressure which broadens in time. We analyze in detail the acceleration of the flow, and its re-collimation by magnetic forces happening at distances of several hundreds of astronomical units. We quantify the impact of magnetic braking in the outflows, which narrows the outflow cavity for the late evolution of the system. We observe the presence of the same jet-driving mechanisms for a wide range of assumptions on the natal environment of the massive protostar, but with changes to their morphology and mechanical feedback into larger scales over time. △ Less

Submitted 4 October, 2022; originally announced October 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 669, A81 (2023)

Snapshot of a magnetohydrodynamic disk wind traced by water maser observations

Authors: Luca Moscadelli, Alberto Sanna, Henrik Beuther, André Oliva, Rolf Kuiper

Abstract: The formation of astrophysical objects of different nature and size, from black holes to gaseous giant planets, involves a disk-jet system, where the disk drives the mass accretion onto a central compact object and the jet is a fast collimated ejection along the disk rotation axis. Magnetohydrodynamic disk winds can provide the link between mass accretion and ejection, which is essential to ensure… ▽ More The formation of astrophysical objects of different nature and size, from black holes to gaseous giant planets, involves a disk-jet system, where the disk drives the mass accretion onto a central compact object and the jet is a fast collimated ejection along the disk rotation axis. Magnetohydrodynamic disk winds can provide the link between mass accretion and ejection, which is essential to ensure that the excess angular momentum is removed from the system and accretion onto the central object can proceed. However, up to now, we have been lacking direct observational proof of disk winds. This work presents a direct view of the velocity field of a disk wind around a forming massive star. Achieving a very high spatial resolution of ~0.05 au, our water maser observations trace the velocities of individual streamlines emerging from the disk orbiting the forming star. We find that, at low elevation above the disk midplane, the flow co-rotates with its launch point in the disk, in agreement with magneto-centrifugal acceleration where the gas is flung away along the magnetic field line anchored to the disk. Beyond the co-rotation point, the flow rises spiraling around the disk rotation axis along a helical magnetic field. We have performed (resistive-radiative-gravito-) magnetohydrodynamic simulations of the formation of a massive star and record the development of a magneto-centrifugally launched jet presenting many properties in agreement with our observations. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Comments: Main paper: 15 pages, 6 figures, 2 tables; Appendix: 21 pages, 7 figures, 1 table

Journal ref: published on-line in Nature Astronomy on August 11th 2022; view-only at: https://rdcu.be/cTukQ ; You Tube video: https://www.youtube.com/watch?v=EcmmvuVWa3Q

Radiative feedback on supermassive star formation: the massive end of the Population III initial mass function

Authors: Daisuke Toyouchi, Kohei Inayoshi, Wenxiu Li, Zoltán Haiman, Rolf Kuiper

Abstract: Supermassive stars (SMSs) with masses of $M_\ast \simeq 10^4$--$10^5~{\rm M_\odot}$ are invoked as possible seeds of high-redshift supermassive black holes, but it remains under debate whether their protostar indeed acquires sufficient mass via gas accretion overcoming radiative feedback. We investigate protostellar growth in dynamically heated atomic-cooling haloes (ACHs) found in recent cosmolog… ▽ More Supermassive stars (SMSs) with masses of $M_\ast \simeq 10^4$--$10^5~{\rm M_\odot}$ are invoked as possible seeds of high-redshift supermassive black holes, but it remains under debate whether their protostar indeed acquires sufficient mass via gas accretion overcoming radiative feedback. We investigate protostellar growth in dynamically heated atomic-cooling haloes (ACHs) found in recent cosmological simulations, performing three-dimensional radiation hydrodynamical (RHD) simulations that consider stellar evolution under variable mass accretion. We find that one of the ACHs feeds the central protostar at rates exceeding a critical value, above which the star evolves in a cool bloating phase and hardly produces ionizing photons. Consequently, the stellar mass reaches $M_\ast \gtrsim 10^4~{\rm M_\odot}$ unimpeded by radiative feedback. In the other ACH, where the mass supply rate is lower, the star spends most of its life as a hot main-sequence star, emitting intense ionizing radiation. Then, the stellar mass growth is terminated around $500~{\rm M_\odot}$ by photoevaporation of the circumstellar disk. A series of our RHD simulations provide a formula of the final stellar mass determined either by stellar feedback or their lifetime as a function of the mass supply rate from the parent cloud in the absence of stellar radiation. Combining the results with the statistical properties of SMS-forming clouds in high-redshift quasar progenitor haloes, we construct a top-heavy mass distribution of primordial stars over $M_\ast \simeq 100$--$10^5~{\rm M_\odot}$, approximately following a power-law spectrum of $\propto M_\ast^{-1.3}$ with a steeper decline at $M_\ast \gtrsim 2 \times 10^4~{\rm M_\odot}$. Their massive BH remnants would be further fed via the dense debris disk, powering "milli-quasars" with a bolometric luminosity of $L_{\rm bol}~\gtrsim~10^{43}~{\rm erg~s^{-1}}$. △ Less

Submitted 29 June, 2022; originally announced June 2022.

Comments: 16 pages, 12 figures, and 3 tables, submitted to MNRAS

Supercritical growth pathway to overmassive black holes at cosmic dawn: coevolution with massive quasar hosts

Authors: Haojie Hu, Kohei Inayoshi, Zoltán Haiman, Wenxiu Li, Eliot Quataert, Rolf Kuiper

Abstract: Observations of the most luminous quasars at high redshifts ($z > 6$) have revealed that the largest supermassive black holes (SMBHs) at those epochs tend to be substantially overmassive relative to their host galaxies compared to the local relations, suggesting they experienced rapid early growth phases. We propose an assembly model for the SMBHs that end up in rare massive… ▽ More Observations of the most luminous quasars at high redshifts ($z > 6$) have revealed that the largest supermassive black holes (SMBHs) at those epochs tend to be substantially overmassive relative to their host galaxies compared to the local relations, suggesting they experienced rapid early growth phases. We propose an assembly model for the SMBHs that end up in rare massive $\sim10^{12}~M_{\odot}$ host halos at $z \sim 6-7$, applying a kinetic feedback prescription for BHs accreting above the Eddington rate, provided by radiation hydrodynamic simulations for the long-term evolution of the accretion-flow structure. The large inflow rates into these halos during their assembly enable the formation of $>10^9~M_{\odot}$ SMBHs by $z \sim 6$, even starting from stellar-mass seeds at $z \sim 30$, and even in the presence of outflows that reduce the BH feeding rate, especially at early times. This mechanism also naturally yields a high BH-to-galaxy mass ratio of $> 0.01$ before the SMBH mass reaches $M_{\rm BH} > 10^9~M_{\odot}$ by $z \sim 6$. These fast-growing SMBH progenitors are bright enough to be detected by upcoming observations with the James Webb Space Telescope over a wide range of redshift ($7 < z < 15$), regardless of how they were seeded. △ Less

Submitted 30 June, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: 4 figures, accepted for publication in ApJ

Long-term evolution of supercritical black hole accretion with outflows: a subgrid feedback model for cosmological simulations

Authors: Haojie Hu, Kohei Inayoshi, Zoltán Haiman, Eliot Quataert, Rolf Kuiper

Abstract: We study the long-term evolution of the global structure of axisymmetric accretion flows onto a black hole (BH) at rates substantially higher than the Eddington value ($\dot{M}_{\rm Edd}$), performing two-dimensional hydrodynamical simulations with and without radiative diffusion. In the high-accretion optically-thick limit, where the radiation energy is efficiently trapped within the inflow, the… ▽ More We study the long-term evolution of the global structure of axisymmetric accretion flows onto a black hole (BH) at rates substantially higher than the Eddington value ($\dot{M}_{\rm Edd}$), performing two-dimensional hydrodynamical simulations with and without radiative diffusion. In the high-accretion optically-thick limit, where the radiation energy is efficiently trapped within the inflow, the accretion flow becomes adiabatic and comprises of turbulent gas in the equatorial region and strong bipolar outflows. As a result, the mass inflow rate decreases toward the center as $\dot{M}_{\rm in}\propto r^{p}$ with $p\sim 0.5-0.7$ and a small fraction of the inflowing gas feeds the nuclear BH. Thus, super-Eddington accretion is sustained only when a larger amount of gas is supplied from larger radii at $> 100-1000~\dot{M}_{\rm Edd}$. The global structure of the flow settles down to a quasi-steady state in millions of the orbital timescale at the BH event horizon, which is $> 10-100$ times longer than that addressed in previous (magneto-)RHD simulation studies. Energy transport via radiative diffusion accelerates the outflow near the poles in the inner region but does not change the overall properties of the accretion flow compared to the cases without diffusion. Based on our simulation results, we provide a mechanical feedback model for super-Eddington accreting BHs. This can be applied as a sub-grid model in large-scale cosmological simulations that do not sufficiently resolve galactic nuclei, and to the formation of the heaviest gravitational-wave sources via accretion in dense environments. △ Less

Submitted 1 June, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: 17 pages, 13 figures, 1 table (re-submitted to APJ, after revision)

A radiation hydrodynamics scheme on adaptive meshes using the Variable Eddington Tensor (VET) closure

Authors: Shyam H. Menon, Christoph Federrath, Mark R. Krumholz, Rolf Kuiper, Benjamin D. Wibking, Manuel Jung

Abstract: We present a new algorithm to solve the equations of radiation hydrodynamics (RHD) in a frequency-integrated, two-moment formulation. Novel features of the algorithm include i) the adoption of a non-local Variable Eddington Tensor (VET) closure for the radiation moment equations, computed with a ray-tracing method, ii) support for adaptive mesh refinement (AMR), iii) use of a time-implicit Godunov… ▽ More We present a new algorithm to solve the equations of radiation hydrodynamics (RHD) in a frequency-integrated, two-moment formulation. Novel features of the algorithm include i) the adoption of a non-local Variable Eddington Tensor (VET) closure for the radiation moment equations, computed with a ray-tracing method, ii) support for adaptive mesh refinement (AMR), iii) use of a time-implicit Godunov method for the hyperbolic transport of radiation, and iv) a fixed-point Picard iteration scheme to accurately handle the stiff nonlinear gas-radiation energy exchange. Tests demonstrate that our scheme works correctly, yields accurate rates of energy and momentum transfer between gas and radiation, and obtains the correct radiation field distribution even in situations where more commonly used -- but less accurate -- closure relations like the Flux-limited Diffusion and Moment-1 approximations fail. Our scheme presents an important step towards performing RHD simulations with increasing spatial and directional accuracy, effectively improving their predictive capabilities. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 7 pages, 2 figures, to appear in the proceedings of IAUS362: The predictive power of computational astrophysics as discovery tool (D. Bisikalo, T. Hanawa, C. Boily, J. Stone, eds.). arXiv admin note: substantial text overlap with arXiv:2202.08778

Recycling of the first atmospheres of embedded planets: Dependence on core mass and optical depth

Authors: T. W. Moldenhauer, R. Kuiper, W. Kley, C. W. Ormel

Abstract: Recent observations found close-in planets with significant atmospheres of hydrogen and helium in great abundance. These are the so-called super-Earths and mini-Neptunes. Their atmospheric composition suggests that they formed early during the gas-rich phase of the circumstellar disk and were able to avoid becoming hot Jupiters. As a possible explanation, recent studies explored the recycling hypo… ▽ More Recent observations found close-in planets with significant atmospheres of hydrogen and helium in great abundance. These are the so-called super-Earths and mini-Neptunes. Their atmospheric composition suggests that they formed early during the gas-rich phase of the circumstellar disk and were able to avoid becoming hot Jupiters. As a possible explanation, recent studies explored the recycling hypothesis and showed that atmosphere-disk recycling is able to fully compensate for radiative cooling and thereby halt Kelvin-Helmholtz contraction to prevent runaway gas accretion. To understand the parameters that determine the efficiency of atmospheric recycling, we extend our earlier studies by exploring the effects of the core mass, the effect of circumstellar gas on sub-Keplerian orbits (headwind), and the optical depth of the surrounding gas on the recycling timescale. Additionally, we analyze their effects on the size and mass of the forming atmosphere. For the explored parameter space, all simulations eventually reach an equilibrium where heating due to hydrodynamic recycling fully compensates radiative cooling. In this equilibrium, the atmosphere-to-core mass ratio stays well below $10 \, \%$, preventing the atmosphere from becoming self-gravitating and entering runaway gas accretion. Higher core masses cause the atmosphere to become turbulent, which further enhances recycling. Even for our highest core mass of $10 \, M_\mathrm{Earth}$, atmosphere-disk recycling is efficient enough to fully compensate for radiative cooling and prevent the atmosphere from becoming self-gravitating. Hence, in-situ formation of hot Jupiters is very unlikely, and migration of gas giants is a key process required to explain their existence. Our findings imply that atmosphere-disk recycling is the most natural explanation for the prevalence of close-in super-Earths and mini-Neptunes. △ Less

Submitted 24 February, 2022; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: 18 pages, 13 figures, Accepted for publication by A&A on 2/7/2022

VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method

Authors: Shyam H. Menon, Christoph Federrath, Mark R. Krumholz, Rolf Kuiper, Benjamin D. Wibking, Manuel Jung

Abstract: We present Variable Eddington Tensor-closed Transport on Adaptive Meshes (\texttt{VETTAM}), a new algorithm to solve the equations of radiation hydrodynamics (RHD) with support for adaptive mesh refinement (AMR) in a frequency-integrated, two-moment formulation. The method is based on a non-local Variable Eddington Tensor (VET) closure computed with a hybrid characteristics scheme for ray tracing.… ▽ More We present Variable Eddington Tensor-closed Transport on Adaptive Meshes (\texttt{VETTAM}), a new algorithm to solve the equations of radiation hydrodynamics (RHD) with support for adaptive mesh refinement (AMR) in a frequency-integrated, two-moment formulation. The method is based on a non-local Variable Eddington Tensor (VET) closure computed with a hybrid characteristics scheme for ray tracing. We use a Godunov method for the hyperbolic transport of radiation with an implicit backwards-Euler temporal update to avoid the explicit timestep constraint imposed by the light-crossing time, and a fixed-point Picard iteration scheme to handle the nonlinear gas-radiation exchange term, with the two implicit update stages jointly iterated to convergence. We also develop a modified wave-speed correction method for AMR, which we find to be crucial for obtaining accurate results in the diffusion regime. We demonstrate the robustness of our scheme with a suite of pure radiation and RHD tests, and show that it successfully captures the streaming, static diffusion, and dynamic diffusion regimes and the spatial transitions between them, casts sharp shadows, and yields accurate results for rates of momentum and energy exchange between radiation and gas. A comparison between different closures for the radiation moment equations, with the Eddington approximation (0th-moment closure) and the $M_1$ approximation (1st-moment closure), demonstrates the advantages of the VET method (2nd-moment closure) over the simpler closure schemes. \texttt{VETTAM} has been coupled to the AMR \texttt{FLASH} (magneto-)hydrodynamics code and we summarize by reporting performance features and bottlenecks of our implementation. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: 23 pages, 15 figures. Accepted for publication in MNRAS

The effect of ultra-violet photon pumping of H$_2$ in dust-deficient protoplanetary disks

Authors: Ayano Komaki, Rolf Kuiper, Naoki Yoshida

Abstract: We perform radiation hydrodynamics simulations to study the structure and evolution of a photoevaporating protoplanetary disk. Ultraviolet and X-ray radiation from the host star heats the disk surface, where H$_2$ pumping also operates efficiently. We run a set of simulations with varying the amount of dust grains, or the dust-to-gas mass ratio, which determines the relative importance between pho… ▽ More We perform radiation hydrodynamics simulations to study the structure and evolution of a photoevaporating protoplanetary disk. Ultraviolet and X-ray radiation from the host star heats the disk surface, where H$_2$ pumping also operates efficiently. We run a set of simulations with varying the amount of dust grains, or the dust-to-gas mass ratio, which determines the relative importance between photoelectric heating and H$_2$ pumping. We show that H$_2$ pumping and X-ray heating contribute stronger to the mass-loss of the disk if the dust-to-gas mass ratio is $\mathcal{D}\leq10^{-3}$.The disk mass loss rate decreases with a lower dust amount, but remains around $10^{-10-11} M_{\odot} {\rm yr}^{-1}$. In such dust-deficient disks, H$_2$ pumping enhances photoevaporation from the inner disk region and shapes the disk mass-loss profile. We thus argue that the late-stage disk evolution is affected by the ultra-violet H$_2$ pumping effect. The mass-loss rates derived from our simulations can be used in the study of long-term disk evolution. △ Less

Submitted 20 March, 2024; v1 submitted 6 February, 2022; originally announced February 2022.

Comments: 12 pages, 3 figures, 1 table

The sharp ALMA view of infall and outflow in the massive protocluster G31.41+0.31

Authors: M. T. Beltrán, V. M. Rivilla, R. Cesaroni, D. Galli, L. Moscadelli, A. Ahmadi, H. Beuther, S. Etoka, C. Goddi, P. D. Klaassen, R. Kuiper, M. S. N. Kumar, A. Lorenzani, T. Peters, Á. Sánchez-Monge, P. Schilke, F. van der Tak, S. Vig

Abstract: Context. To better understand the formation of high-mass stars, it is fundamental to investigate how matter accretes onto young massive stars, how it is ejected, and how all this differs from the low-mass case. The massive protocluster G31.41+0.31 is the ideal target to study all these processes because observations at millimeter and centimeter wavelengths have resolved the emission of the Main co… ▽ More Context. To better understand the formation of high-mass stars, it is fundamental to investigate how matter accretes onto young massive stars, how it is ejected, and how all this differs from the low-mass case. The massive protocluster G31.41+0.31 is the ideal target to study all these processes because observations at millimeter and centimeter wavelengths have resolved the emission of the Main core into at least four massive dust continuum sources, named A, B, C, and D, within 1" or 0.018 pc, and have identified signatures of infall and several outflows associated with the core. Aims. We study the interplay between infall and outflow in G31.41+0.31 by investigating at a spatial resolution of a few 100 au their properties and their possible impact on the core. Methods. We carried out molecular line observations of typical high-density tracers, such as CH3CN or H2CO, and shock and outflow tracers, such as SiO, with ALMA at 1.4 mm that achieved an angular resolution of 0.09" (340 au). Results. The observations have revealed inverse P-Cygni profiles in CH3CN and H2CO toward the four sources embedded in the Main core, suggesting that all of them are undergoing collapse. The infall rates, estimated from the red-shifted absorption are on the order of 1E-2 Msun/yr. The individual infall rates imply that the accretion timescale of the Main core is an order of magnitude smaller than its rotation timescale. This confirms that rotating toroids such as the G31 Main core are non-equilibrium, transient collapsing structures that need to be constantly replenished with fresh material from a large-scale reservoir. For sources B, C, and D, the infall could be accelerating inside the sources, while for source A, the presence of a second emission component complicates the interpretation. The SiO observations have revealed the presence of at least six outflows in the G31.41+0.31 star-forming region, ... △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 17 pages, 9 figures, 2 tables. Accepted by A&A

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

Accreting protoplanets: Spectral signatures and magnitude of gas and dust extinction at H alpha

Authors: G. -D. Marleau, Y. Aoyama, R. Kuiper, K. Follette, N. J. Turner, G. Cugno, C. F. Manara, S. Y. Haffert, D. Kitzmann, S. C. Ringqvist, K. R. Wagner, R. van Boekel, S. Sallum, M. Janson. T. O. B. Schmidt, L. Venuti, Ch. Lovis, C. Mordasini

Abstract: Accreting planets have been seen at Ha (H alpha), but targeted searches have not been fruitful. For planets, accretion tracers should come from the shock itself, exposing them to extinction by the accreting material. High-resolution (R>5e4) spectrographs at Ha should soon allow studying how the incoming material shapes the line profile. We calculate how much the gas and dust accreting onto a plane… ▽ More Accreting planets have been seen at Ha (H alpha), but targeted searches have not been fruitful. For planets, accretion tracers should come from the shock itself, exposing them to extinction by the accreting material. High-resolution (R>5e4) spectrographs at Ha should soon allow studying how the incoming material shapes the line profile. We calculate how much the gas and dust accreting onto a planet reduce the Ha flux from the shock at the planetary surface and how they affect line shapes. We also study the absorption-modified relationship between Ha luminosity and Mdot. We compute the high-resolution radiative transfer of the Ha line using a 1D velocity-density-temperature structure for the inflowing matter in three representative accretion geometries: spherical symmetry, polar inflow, and magnetospheric accretion. For each, we explore wide ranges of Mdot and planet mass M. We use detailed gas opacities and estimate dust opacities. At Mdot<3e-6 MJ/yr, gas extinction is negligible for spherical or polar inflow and at most A_Ha<0.5 mag for magnetospheric accretion. Up to Mdot~3e-4 MJ/yr, the gas has A_Ha<4 mag. This decreases with M. We estimate realistic dust opacities at Ha as ~0.01-10 cm^2/g, i.e., 10-1e4 times lower than in the ISM. Extinction flattens the L_Ha-Mdot relationship, which becomes non-monotonic with a maximum L_Ha~1e-4 LSun near Mdot~1e-4 MJ/yr for M~10 MJ. In magnetospheric accretion, the gas can introduce features in line profiles, but the velocity gradient smears them out in other geometries. For most of parameter space, extinction by the accreting matter should be negligible, simplifying interpretation of observations, especially for planets in gaps. At high Mdot, strong absorption reduces the Ha flux, and some measurements can be interpreted as two Mdot values. Line profiles at R~1e5 can provide complex constraints on the accretion flow's thermal-dynamical structure. △ Less

Submitted 23 November, 2021; originally announced November 2021.

Comments: 25 pages, 11 figures (main text) + 3 pages, 4 figures (three appendices). In press at A&A. Comments welcome

Journal ref: A&A 657, A38 (2022)

Rapid growth of seed black holes during early bulge formation

Authors: Kohei Inayoshi, Riouhei Nakatani, Daisuke Toyouchi, Takashi Hosokawa, Rolf Kuiper, Masafusa Onoue

Abstract: We study the early growth of massive seed black holes (BHs) via accretion in protogalactic nuclei where the stellar bulge component is assembled, performing axisymmetric two-dimensional radiation hydrodynamical simulations. We find that when a seed BH with $M_\bullet \sim 10^5~M_\odot$ is embedded in dense metal-poor gas ($Z=0.01~Z_\odot$) with a density of $\gtrsim 100~{\rm cm}^{-3}$ and bulge st… ▽ More We study the early growth of massive seed black holes (BHs) via accretion in protogalactic nuclei where the stellar bulge component is assembled, performing axisymmetric two-dimensional radiation hydrodynamical simulations. We find that when a seed BH with $M_\bullet \sim 10^5~M_\odot$ is embedded in dense metal-poor gas ($Z=0.01~Z_\odot$) with a density of $\gtrsim 100~{\rm cm}^{-3}$ and bulge stars with a total mass of $M_\star \gtrsim 100~M_\bullet$, a massive gaseous disk feeds the BH efficiently at rates of $\gtrsim 0.3-1~M_\odot~{\rm yr}^{-1}$ and the BH mass increases nearly tenfold within $\sim 2$ Myr. This rapid accretion phase lasts until a good fraction of the gas bounded within the bulge accretes onto the BH, although the feeding rate is regulated owing to strong outflows driven by ionizing radiation emitted from the accreting BH. The transient growing mode can be triggered for seed BHs formed in massive dark-matter halos with masses of $\gtrsim 10^9~M_\odot$ at $z\sim 15-20$ (the virial temperature is $T_{\rm vir}\simeq 10^5~{\rm K}$). The host halos are heavier and rarer than those of typical first galaxies, but are more likely to end up in quasar hosts by $z\simeq 6$. This mechanism naturally yields a mass ratio of $M_\bullet/M_\star >0.01$ higher than the value seen in the local universe and the existence of such overmassive BHs provides us a unique opportunity of detecting highly accreting seed BHs at $z\sim 15$ with AB magnitude of $m_{\rm AB} \sim26 - 29$ mag at $2~μ{\rm m}$ (rest-frame 10 eV) by the upcoming observations by the James Webb Space Telescope and Nancy Grace Roman Space Telescope. △ Less

Submitted 20 January, 2022; v1 submitted 20 October, 2021; originally announced October 2021.

Comments: 28 pages, 11 figures, accepted for publication to ApJ. Movies of the simulation are available at https://www.youtube.com/playlist?list=PLqDrKvxmu0l4fyJ74u5ZYyQ8PmwuNZA32

Stability of an Ionization Front in Bondi Accretion

Authors: Eric Keto, Rolf Kuiper

Abstract: Spherical Bondi accretion is used in astrophysics as an approximation to investigate many types of accretion processes. Two-phase accretion flows that transition from neutral to ionized have observational support in high-mass star formation, and have application to accretion flows around any ionizing source, but the hydrodynamic stability of two-phase Bondi accretion is not understood. With both s… ▽ More Spherical Bondi accretion is used in astrophysics as an approximation to investigate many types of accretion processes. Two-phase accretion flows that transition from neutral to ionized have observational support in high-mass star formation, and have application to accretion flows around any ionizing source, but the hydrodynamic stability of two-phase Bondi accretion is not understood. With both semi-analytic and fully numerical methods we find that these flows may be stable, conditionally stable or unstable depending on the initial conditions. The transition from an R-type to a D-type ionization front plays a key role in conditionally stable and unstable flows. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Comments: Accepted MNRAS; 8 pages, 5 figures

Clustered star formation at early evolutionary stages. Physical and chemical analysis of the young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953

Authors: C. Gieser, H. Beuther, D. Semenov, S. Suri, J. D. Soler, H. Linz, J. Syed, Th. Henning, S. Feng, T. Möller, A. Palau, J. M. Winters, M. T. Beltrán, R. Kuiper, L. Moscadelli, P. Klaassen, J. S. Urquhart, T. Peters, S. N. Longmore, Á. Sánchez-Monge, R. Galván-Madrid, R. E. Pudritz, K. G. Johnston

Abstract: We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from t… ▽ More We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from the 1.3 mm continuum visibility profiles and the radial temperature profiles are estimated from H2CO rotation temperature maps. Molecular column densities are derived with the line fitting tool XCLASS. The physical and chemical properties are combined by applying the physical-chemical model MUSCLE in order to constrain the chemical timescales of a few line-rich cores. The morphology and spatial correlations of the molecular emission are analyzed using the HOG method. The mid-infrared data show that both regions contain a cluster of young stellar objects. Bipolar molecular outflows are observed in the CO 2-1 transition toward the strong mm cores indicating protostellar activity. We find strong molecular emission of SO, SiO, H2CO, and CH3OH in locations which are not associated with the mm cores. These shocked knots can be either associated with the bipolar outflows or, in the case of ISOSS J23053+5953, with a colliding flow that creates a large shocked region between the mm cores. The mean chemical timescale of the cores is lower (20 000 yr) compared to that of the sources of the more evolved CORE sample (60 000 yr). With the HOG method, we find that the spatial emission of species tracing the extended emission and of shock-tracing molecules are well correlated within transitions of these groups. △ Less

Submitted 14 October, 2021; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: 34 pages, 17 figures, accepted for publication in A&A

Journal ref: A&A 657, A3 (2022)

Discovery of Molecular Line Polarization in the Disk of TW Hya

Authors: Richard Teague, Chat L. H. Hull, Stéphane Guilloteau, Edwin A. Bergin, Anne Dutrey, Thomas Henning, Rolf Kuiper, Dmitry Semenov, Ian W. Stephens, Wouter H. T. Vlemmings

Abstract: We report observations of polarized line and continuum emission from the disk of TW~Hya using the Atacama Large Millimeter/submillimeter Array. We target three emission lines, $^{12}$CO (3-2), $^{13}$CO (3-2) and CS (7-6), to search for linear polarization due to the Goldreich-Kylafis effect, while simultaneously tracing the continuum polarization morphology at 332\,GHz (900\,\micron{}), achieving… ▽ More We report observations of polarized line and continuum emission from the disk of TW~Hya using the Atacama Large Millimeter/submillimeter Array. We target three emission lines, $^{12}$CO (3-2), $^{13}$CO (3-2) and CS (7-6), to search for linear polarization due to the Goldreich-Kylafis effect, while simultaneously tracing the continuum polarization morphology at 332\,GHz (900\,\micron{}), achieving a spatial resolution of 0.5\arcsec{} (30~au). We detect linear polarization in the dust continuum emission; the polarization position angles show an azimuthal morphology, and the median polarization fraction is $\sim$\,0.2\%, comparable to previous, lower frequency observations. Adopting a `shift-and-stack' technique to boost the sensitivity of the data, combined with a linear combination of the $Q$ and $U$ components to account for their azimuthal dependence, we detect weak linear polarization of $^{12}$CO and $^{13}$CO line emission at a $\sim 10σ$ and $\sim 5σ$ significance, respectively. The polarization was detected in the line wings, reaching a peak polarization fraction of $\sim 5\%$ and $\sim 3\%$ for the two molecules between disk radii of 0.5" and 1". The sign of the polarization was found to flip from the blue-shifted side of the emission to the red-shifted side, suggesting a complex, asymmetric polarization morphology. Polarization is not robustly detected for the CS emission; however, a tentative signal, comparable in morphology to that found for the $^{12}$CO and $^{13}$CO emission, is found at a $\lesssim 3σ$ significance. We are able to reconstruct a polarization morphology, consistent with the azimuthally averaged profiles, under the assumption that this is also azimuthally symmetric, which can be compared with future higher-sensitivity observations. △ Less

Submitted 22 September, 2021; v1 submitted 19 September, 2021; originally announced September 2021.

Comments: Accepted by ApJ

Disk fragmentation in high-mass star formation. High-resolution observations towards AFGL 2591-VLA 3

Authors: S. Suri, H. Beuther, C. Gieser, A. Ahmadi, Á. Sánchez-Monge, J. M. Winters, H. Linz, Th. Henning, M. T. Beltrán, F. Bosco, R. Cesaroni, T. Csengeri, S. Feng, M. G. Hoare, K. G. Johnston, P. Klaasen, R. Kuiper, S. Leurini, S. Longmore, S. Lumsden, L. Maud, L. Moscadelli, T. Möller, A. Palau, T. Peters , et al. (7 additional authors not shown)

Abstract: Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk ar… ▽ More Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk around the high-mass star forming region AFGL 2591-VLA 3 which was hypothesized to be fragmenting based on the observations that show multiple outflow directions. We use a new set of high-resolution (0.19 arcsec) IRAM/NOEMA observations at 843 micron towards VLA 3 which allow us to resolve its disk, characterize the fragmentation, and study its kinematics. In addition to the 843 micron continuum emission, our spectral setup targets warm dense gas and outflow tracers such as HCN, HC$_3$N and SO$_2$, as well as vibrationally excited HCN lines. The high resolution continuum and line emission maps reveal multiple fragments with subsolar masses within the inner 1000 AU of VLA 3. Furthermore, the velocity field of the inner disk observed at 843 micron shows a similar behavior to that of the larger scale velocity field studied in the CORE project at 1.37 mm. We present the first observational evidence for disk fragmentation towards AFGL 2591-VLA 3, a source that was thought to be a single high-mass core. While the fragments themselves are low-mass, the rotation of the disk is dominated by the protostar with a mass of 10.3$\pm 1.8~M_{\odot}$. These data also show that NOEMA Band 4 can obtain the highest currently achievable spatial resolution at (sub-)mm wavelengths in observations of strong northern sources. △ Less

Submitted 10 September, 2021; originally announced September 2021.

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

Journal ref: A&A 655, A84 (2021)

Physical conditions in the warped accretion disk of a massive star. 349 GHz ALMA observations of G023.01$-$00.41

Authors: A. Sanna, A. Giannetti, M. Bonfand, L. Moscadelli, R. Kuiper, J. Brand, R. Cesaroni, A. Caratti o Garatti, T. Pillai, K. M. Menten

Abstract: Young massive stars warm up the large amount of gas and dust which condenses in their vicinity, exciting a forest of lines from different molecular species. Their line brightness is a diagnostic tool of the gas physical conditions locally, which we use to set constraints on the environment where massive stars form. We made use of the Atacama Large Millimeter/submillimeter Array at frequencies near… ▽ More Young massive stars warm up the large amount of gas and dust which condenses in their vicinity, exciting a forest of lines from different molecular species. Their line brightness is a diagnostic tool of the gas physical conditions locally, which we use to set constraints on the environment where massive stars form. We made use of the Atacama Large Millimeter/submillimeter Array at frequencies near 349 GHz, with an angular resolution of $0.1''$, to observe the methyl cyanide (CH$_3$CN) emission which arises from the accretion disk of a young massive star. We sample the disk midplane with twelve distinct beams, where we get an independent measure of the gas (and dust) physical conditions. The accretion disk extends above the midplane showing a double-armed spiral morphology projected onto the plane of the sky, which we sample with ten additional beams: along these apparent spiral features, gas undergoes velocity gradients of about $\rm 1 km s^{-1}$ per 2000 au. The gas temperature (T) rises symmetrically along each side of the disk, from about 98 K at 3000 au to 289 K at 250 au, following a power law with radius, R$^{-0.43}$. The CH$_3$CN column density (N) increases from $\rm 9.2\times10^{15} cm^{-2}$ to $\rm 8.7\times10^{17} cm^{-2}$ at the same radii, following a power law with radius, R$^{-1.8}$. In the framework of a circular gaseous disk observed approximately edge-on, we infer an H$_2$ volume density in excess of $\rm 4.8\times10^9 cm^{-3}$ at a distance of 250 au from the star. We study the disk stability against fragmentation following the methodology by Kratter et al. (2010), appropriate under rapid accretion, and we show that the disk is marginally prone to fragmentation along its whole extent. △ Less

Submitted 12 July, 2021; originally announced July 2021.

Comments: 17 pages, 10 figures, 3 tables, accepted by Astronomy & Astrophysics

Journal ref: A&A 655, A72 (2021)

Fragmentation and kinematics in high-mass star formation: CORE-extension targeting two very young high-mass star-forming regions

Authors: H. Beuther, C. Gieser, S. Suri, H. Linz, P. Klaassen, D. Semenov, J. M. Winters, Th. Henning, J. D. Soler, J. S. Urquhart, J. Syed, S . Feng, T. Moeller, M. T. Beltran, A. Sanchez-Monge, S. N. Longmore, T. Peters, J. Ballesteros-Paredes, P. Schilke, L. Moscadelli, A. Palau, R. Cesaroni, S. Lumsden, R. Pudritz, F. Wyrowski , et al. (2 additional authors not shown)

Abstract: Context: The formation of high-mass star-forming regions from their parental gas cloud and the subsequent fragmentation processes lie at the heart of star formation research. Aims: We aim to study the dynamical and fragmentation properties at very early evolutionary stages of high-mass star formation. Methods: Employing the NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30m telescope, we… ▽ More Context: The formation of high-mass star-forming regions from their parental gas cloud and the subsequent fragmentation processes lie at the heart of star formation research. Aims: We aim to study the dynamical and fragmentation properties at very early evolutionary stages of high-mass star formation. Methods: Employing the NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30m telescope, we observed two young high-mass star-forming regions, ISOSS22478 and ISOSS23053, in the 1.3mm continuum and spectral line emission at a high angular resolution (~0.8''). Results: We resolved 29 cores that are mostly located along filament-like structures. Depending on the temperature assumption, these cores follow a mass-size relation of approximately M~r^2.0, corresponding to constant mean column densities. However, with different temperature assumptions, a steeper mass-size relation up to M~r^3.0, which would be more likely to correspond to constant mean volume densities, cannot be ruled out. The correlation of the core masses with their nearest neighbor separations is consistent with thermal Jeans fragmentation. We found hardly any core separations at the spatial resolution limit, indicating that the data resolve the large-scale fragmentation well. Although the kinematics of the two regions appear very different at first sight - multiple velocity components along filaments in ISOSS22478 versus a steep velocity gradient of more than 50km/s/pc in ISOSS23053 - the findings can be explained within the framework of a dynamical cloud collapse scenario. Conclusions: While our data are consistent with a dynamical cloud collapse scenario and subsequent thermal Jeans fragmentation, the importance of additional environmental properties, such as the magnetization of the gas or external shocks triggering converging gas flows, is nonetheless not as well constrained and would require future investigation. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 16 pages, 16 figures, accepted for Astronomy and Astrophysics, a higher-resolution version can also be found at https://www.mpia.de/homes/beuther/papers.html

Journal ref: A&A 649, A113 (2021)

Fragmentation in the massive G31.41+0.31 protocluster

Authors: M. T. Beltrán, V. M. Rivilla, R. Cesaroni, L. T. Maud, D. Galli, L. Moscadelli, A. Lorenzani, A. Ahmadi, H. Beuther, T. Csengeri, S. Etoka, C. Goddi, P. D. Klaassen, R. Kuiper, M. S. N. Kumar, T. Peters, Á. Sánchez-Monge, P. Schilke, F. van der Tak, S. Vig, H. Zinnecker

Abstract: Context. ALMA observations at 1.4 mm and 0.2'' (750au) angular resolution of the Main core in the high-mass star forming region G31.41+0.31 have revealed a puzzling scenario: on the one hand, the continuum emission looks very homogeneous and the core appears to undergo solid-body rotation, suggesting a monolithic core stabilized by the magnetic field; on the other hand, rotation and infall speed u… ▽ More Context. ALMA observations at 1.4 mm and 0.2'' (750au) angular resolution of the Main core in the high-mass star forming region G31.41+0.31 have revealed a puzzling scenario: on the one hand, the continuum emission looks very homogeneous and the core appears to undergo solid-body rotation, suggesting a monolithic core stabilized by the magnetic field; on the other hand, rotation and infall speed up toward the core center, where two massive embedded free-free continuum sources have been detected, pointing to an unstable core having undergone fragmentation. Aims. To establish whether the Main core is indeed monolithic or its homogeneous appearance is due to a combination of large dust opacity and low angular resolution, we carried out millimeter observations at higher angular resolution and different wavelengths. Methods. We carried out ALMA observations at 1.4 mm and 3.5 mm that achieved angular resolutions of 0.1''(375 au) and 0.075'' (280 au), respectively. VLA observations at 7 mm and 1.3 cm at even higher angular resolution, 0.05'' (190 au) and 0.07'' (260 au), respectively, were also carried out to better study the nature of the free-free continuum sources detected in the core. Results. The millimeter continuum emission of the Main core has been clearly resolved into at least four sources, A, B, C, and D, within 1'', indicating that the core is not monolithic. The deconvolved radii of the dust emission of the sources, estimated at 3.5 mm, are 400-500au, their masses range from 15 to 26 Msun, and their number densities are several 1E9 cm-3. Sources A and B, located closer to the center of the core and separated by 750 au, are clearly associated with two free-free continuum sources, likely thermal radio jets, and are the brightest in the core. The spectral energy distribution of these two sources and their masses and sizes are similar and suggest a common origin. △ Less

Submitted 8 March, 2021; originally announced March 2021.

Comments: 18 pages,9 figures, 5 tables. Accepted by A&A

Journal ref: A&A 648, A100 (2021)

The physical and chemical structure of high-mass star-forming regions. Unraveling chemical complexity with the NOEMA large program "CORE"

Authors: C. Gieser, H. Beuther, D. Semenov, A. Ahmadi, S. Suri, T. Möller, M. T. Beltran, P. Klaassen, Q. Zhang, J. S. Urquhart, Th. Henning, S. Feng, R. Galván-Madrid, V. de Souza Magalhães, L. Moscadelli, S. Longmore, S. Leurini, R. Kuiper, T. Peters, K. M. Menten, T. Csengeri, G. Fuller, F. Wyrowski, S. Lumsden, Á. Sánchez-Monge , et al. (8 additional authors not shown)

Abstract: We use sub-arcsecond resolution ($\sim$0.4$''$) observations with NOEMA at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. We combine the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structure of these regions are determined by fitting H2CO and CH3CN line emission. The density profil… ▽ More We use sub-arcsecond resolution ($\sim$0.4$''$) observations with NOEMA at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. We combine the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structure of these regions are determined by fitting H2CO and CH3CN line emission. The density profiles are inferred from the 1.37 mm continuum visibilities. The column densities of 11 different species are determined by fitting the emission lines with XCLASS. Within the 18 observed regions, we identify 22 individual cores with associated 1.37 mm continuum emission and with a radially decreasing temperature profile. We find an average temperature power-law index of q = 0.4$\pm$0.1 and an average density power-law index of p = 2.0$\pm$0.2 on scales on the order of several 1 000 au. Comparing these results with values of p derived in the literature suggest that the density profiles remain unchanged from clump to core scales. The column densities relative to N(C18O) between pairs of dense gas tracers show tight correlations. We apply the physical-chemical model MUSCLE to the derived column densities of each core and find a mean chemical age of $\sim$60 000 yrs and an age spread of 20 000-100 000 yrs. With this paper we release all data products of the CORE project available at https://www.mpia.de/core. The CORE sample reveals well constrained density and temperature power-law distributions. Furthermore, we characterize a large variety in molecular richness that can be explained by an age spread confirmed by our physical-chemical modeling. The hot molecular cores show the most emission lines, but we also find evolved cores at an evolutionary stage, in which most molecules are destroyed and thus the spectra appear line-poor again. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Comments: 29 pages, 14 figures, accepted for publication in A&A

Journal ref: A&A 648, A66 (2021)

Multi-scale view of star formation in IRAS 21078+5211: From clump fragmentation to disk wind

Authors: L. Moscadelli, H. Beuther, A. Ahmadi, C. Gieser, F. Massi, R. Cesaroni, Á. Sánchez-Monge, F. Bacciotti, M. T. Beltrán, T. Csengeri, R. Galván-Madrid, Th. Henning, P. D. Klaassen, R. Kuiper, S. Leurini, S. N. Longmore, L. T. Maud, T. Möller, A. Palau, T. Peters, R. E. Pudritz, A. Sanna, D. Semenov, J. S. Urquhart, J. M. Winters , et al. (1 additional authors not shown)

Abstract: In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we deri… ▽ More In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of about 10^{-4}~M_{sun}~yr^{-1}. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. We reveal an extended (size about 0.1~pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by a (previously discovered) compact (size <= 1000~au) radio jet, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (14~km/s over 500~au) Vlsr gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6+/-2.0~M_{sun}. The water masers (previously observed with VLBI) emerge within 100-300~au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be about 16~au and the corresponding launching radius and terminal velocity to be about 2.2~au and 200~km/s, respectively. △ Less

Submitted 9 February, 2021; originally announced February 2021.

Comments: 25 pages, 16 figures

Journal ref: A&A 647, A114 (2021)

Steady State by Recycling prevents Premature Collapse of Protoplanetary Atmospheres

Authors: T. W. Moldenhauer, R. Kuiper, W. Kley, C. W. Ormel

Abstract: In recent years, space missions such as Kepler and TESS have discovered many close-in planets with significant atmospheres consisting of hydrogen and helium: mini-Neptunes. This indicates that these planets formed early in gas-rich disks while avoiding the runaway gas accretion that would otherwise have turned them into hot-Jupiters. A solution is to invoke a long Kelvin-Helmholtz contraction (or… ▽ More In recent years, space missions such as Kepler and TESS have discovered many close-in planets with significant atmospheres consisting of hydrogen and helium: mini-Neptunes. This indicates that these planets formed early in gas-rich disks while avoiding the runaway gas accretion that would otherwise have turned them into hot-Jupiters. A solution is to invoke a long Kelvin-Helmholtz contraction (or cooling) timescale, but it has also been suggested that thermodynamical cooling can be prevented by hydrodynamical planet atmosphere-disk recycling. We investigate the efficacy of the recycling hypothesis in preventing the collapse of the atmosphere, check for the existence of a steady state configuration, and determine the final atmospheric mass to core mass ratio. We use three-dimensional radiation-hydrodynamic simulations to model the formation of planetary proto-atmospheres. Equations are solved in a local frame centered on the planet. Ignoring small oscillations that average to zero over time, the simulations converge to a steady state where the velocity field of the gas becomes constant in time. In a steady state, the energy loss by radiative cooling is fully compensated by the recycling of the low entropy gas in the planetary atmosphere with high entropy gas from the circumstellar disk. For close-in planets, recycling naturally halts the cooling of planetary proto-atmospheres, preventing them from contracting toward the runaway regime and collapsing into gas giants. △ Less

Submitted 1 February, 2021; originally announced February 2021.

Comments: 5 pages, 4 figures, accepted to A&A

Journal ref: A&A 646, L11 (2021)

On the compressive nature of turbulence driven by ionising feedback in the pillars of the Carina Nebula

Authors: Shyam H. Menon, Christoph Federrath, Pamela Klaassen, Rolf Kuiper, Megan Reiter

Abstract: The ionizing radiation of massive stars sculpts the surrounding neutral gas into pillar-like structures. Direct signatures of star formation through outflows and jets are observed in these structures, typically at their tips. Recent numerical simulations have suggested that this star formation could potentially be triggered by photoionising radiation, driving compressive modes of turbulence in the… ▽ More The ionizing radiation of massive stars sculpts the surrounding neutral gas into pillar-like structures. Direct signatures of star formation through outflows and jets are observed in these structures, typically at their tips. Recent numerical simulations have suggested that this star formation could potentially be triggered by photoionising radiation, driving compressive modes of turbulence in the pillars. In this study we use recent high-resolution ALMA observations of $^{12}\mathrm{CO}$, $^{13}\mathrm{CO}$, and $\mathrm{C}^{18} \mathrm{O}, \; J = 2-1$ emission to test this hypothesis for pillars in the Carina Nebula. We analyse column density and intensity-weighted velocity maps, and subtract any large-scale bulk motions in the plane of the sky to isolate the turbulent motions. We then reconstruct the dominant turbulence driving mode in the pillars, by computing the turbulence driving parameter $b$, characterised by the relation $σ_{ρ/ρ_0} = b \mathcal{M}$ between the standard deviation of the density contrast $σ_{ρ/ρ_0}$ (with gas density $ρ$ and its average $ρ_0$) and the turbulent Mach number $\mathcal{M}$. We find values of \mbox{$b\sim0.7$--$1.0$} for most of the pillars, suggesting that predominantly compressive modes of turbulence are driven in the pillars by the ionising radiation from nearby massive stars. We find that this range of $b$ values can produce star formation rates in the pillars that are a factor $\sim 3$ greater than with $b \sim 0.5$, a typical average value of $b$ for spiral-arm molecular clouds. Our results provide further evidence for the potential triggering of star formation in pillars through compressive turbulent motions. △ Less

Submitted 19 October, 2020; originally announced October 2020.

Comments: 18 pages, 8 figures. Accepted for publication in MNRAS

Gravity and rotation drag the magnetic field in high-mass star formation

Authors: H. Beuther, J. D. Soler, H. Linz, Th. Henning, C. Gieser, R. Kuiper, W. Vlemmings, P. Hennebelle, S. Feng, R. Smith, A. Ahmadi

Abstract: The formation of hot stars out of the cold interstellar medium lies at the heart of astrophysical research. Understanding the importance of magnetic fields during star formation remains a major challenge. With the advent of the Atacama Large Millimeter Array, the potential to study magnetic fields by polarization observations has tremendously progressed. However, the major question remains how muc… ▽ More The formation of hot stars out of the cold interstellar medium lies at the heart of astrophysical research. Understanding the importance of magnetic fields during star formation remains a major challenge. With the advent of the Atacama Large Millimeter Array, the potential to study magnetic fields by polarization observations has tremendously progressed. However, the major question remains how much magnetic fields shape the star formation process or whether gravity is largely dominating. Here, we show that for the high-mass star-forming region G327.3 the magnetic field morphology appears to be dominantly shaped by the gravitational contraction of the central massive gas core where the star formation proceeds. We find that in the outer parts of the region, the magnetic field is directed toward the gravitational center of the region. Filamentary structures feeding the central core exhibit U-shaped magnetic field morphologies directed toward the gravitational center as well, again showing the gravitational drag toward the center. The inner part then shows rotational signatures, potentially associated with an embedded disk, and there the magnetic field morphology appears to be rotationally dominated. Hence, our results demonstrate that for this region gravity and rotation are dominating the dynamics and shaping the magnetic field morphology. △ Less

Submitted 12 October, 2020; originally announced October 2020.

Comments: 10 pages, 4 figures, accepted for the Astrophysical Journal, also available at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

Super-Eddington mass growth of intermediate-mass black holes embedded in dusty circumnuclear disks

Authors: Daisuke Toyouchi, Kohei Inayoshi, Takashi Hosokawa, Rolf Kuiper

Abstract: We perform the first three-dimensional radiation hydrodynamical simulations that investigate the growth of intermediate-mass BHs (IMBHs) embedded in massive self-gravitating, dusty nuclear accretion disks. We explore the dependence of mass accretion efficiency on the gas metallicity $Z$ and mass injection at super-Eddington accretion rates from the outer galactic disk $\dot{M}_{\rm in}$, and find… ▽ More We perform the first three-dimensional radiation hydrodynamical simulations that investigate the growth of intermediate-mass BHs (IMBHs) embedded in massive self-gravitating, dusty nuclear accretion disks. We explore the dependence of mass accretion efficiency on the gas metallicity $Z$ and mass injection at super-Eddington accretion rates from the outer galactic disk $\dot{M}_{\rm in}$, and find that the central BH can be fed at rates exceeding the Eddington rate only when the dusty disk becomes sufficiently optically thick to ionizing radiation. In this case, mass outflows from the disk owing to photoevaporation is suppressed and thus a large fraction ($\gtrsim 40\%$) of the mass injection rate can feed the central BH. The conditions are expressed as $\dot{M}_{\rm in} > 2.2\times 10^{-1}~M_\odot ~{\rm yr}^{-1} (1+Z/10^{-2}~Z_\odot)^{-1}(c_{\rm s}/10~{\rm km~s}^{-1})$, where $c_{\rm s}$ is the sound speed in the gaseous disk. With increasing numerical resolution, vigorous disk fragmentation reduces the disk surface density and dynamical heating by formed clumps makes the disk thickness higher. As a result, the photoevaorative mass-loss rate rises and thus the critical injection rate increases for fixed metallicity. This process enables super-Eddington growth of BHs until the BH mass reaches $M_{\rm BH} \sim 10^{7-8}~M_\odot$, depending on the properties of the host dark-matter halo and metal-enrichment history. In the assembly of protogalaxies, seed BHs that form in overdense regions with a mass variance of 3-4$σ$ at $z\sim 15-20$ are able to undergo short periods of their rapid growth and transits into the Eddington-limited growth phase afterwards to be supermassive BHs observed at $z>6-7$. △ Less

Submitted 30 September, 2020; originally announced September 2020.

Comments: 23 pages, 10 figures, submitted to ApJ

Makemake + Sedna: A Continuum Radiation Transport and Photoionization Framework for Astrophysical Newtonian Fluid Dynamics

Authors: Rolf Kuiper, Harold W. Yorke, Andrea Mignone

Abstract: Astrophysical fluid flow studies often encompass a wide range of physical processes to account for the complexity of the system under consideration. In addition to gravity, a proper treatment of thermodynamic processes via continuum radiation transport and/or photoionization is becoming the state of the art. We present a major update of our continuum radiation transport module, MAKEMAKE, and a new… ▽ More Astrophysical fluid flow studies often encompass a wide range of physical processes to account for the complexity of the system under consideration. In addition to gravity, a proper treatment of thermodynamic processes via continuum radiation transport and/or photoionization is becoming the state of the art. We present a major update of our continuum radiation transport module, MAKEMAKE, and a newly developed module for photoionization, SEDNA, coupled to the magnetohydrodynamics code PLUTO. These extensions are currently not publicly available; access can be granted on a case-by-case basis. We explain the theoretical background of the equations solved, elaborate on the numerical layout, and present a comprehensive test suite for radiation-ionization hydrodynamics. The grid-based radiation and ionization modules support static one-dimensional, two-dimensional, and three-dimensional grids in Cartesian, cylindrical, and spherical coordinates. Each module splits the radiation field into two components, one originating directly from a point source - solved using a ray-tracing scheme - and a diffuse component - solved with a three-dimensional flux-limited diffusion (FLD) solver. The FLD solver for the continuum radiation transport makes use of either the equilibrium one-temperature approach or the linearization two-temperature approach. The FLD solver for the photoionization module enables accounting for the temporal evolution of the radiation field from direct recombination of free electrons into hydrogen's ground state as an alternative to the on-the-spot approximation. A brief overview of completed and ongoing scientific studies is given to explicitly illustrate the multipurpose nature of the numerical framework presented. △ Less

Submitted 25 September, 2020; originally announced September 2020.

Journal ref: The Astrophysical Journal Supplement Series, 2020, Volume 250, Issue 1, id.13, 28 pp

Photoevaporation of Grain-Depleted Protoplanetary Disks around Intermediate-Mass Stars: Investigating Possibility of Gas-Rich Debris Disks as Protoplanetary Remnants

Authors: Riouhei Nakatani, Hiroshi Kobayashi, Rolf Kuiper, Hideko Nomura, Yuri Aikawa

Abstract: Debris disks are classically considered to be gas-less systems, but recent (sub)millimeter observations have detected tens of those with rich gas content. The origin of the gas component remains unclear; namely, it can be protoplanetary remnants and/or secondary products deriving from large bodies. In order to be protoplanetary in origin, the gas component of the parental protoplanetary disk is re… ▽ More Debris disks are classically considered to be gas-less systems, but recent (sub)millimeter observations have detected tens of those with rich gas content. The origin of the gas component remains unclear; namely, it can be protoplanetary remnants and/or secondary products deriving from large bodies. In order to be protoplanetary in origin, the gas component of the parental protoplanetary disk is required to survive for $\gtrsim10{\,\rm Myr}$. However, previous models predict $\lesssim 10{\,\rm Myr}$ lifetimes because of efficient photoevaporation at the late stage of disk evolution. In the present study, we investigate photoevaporation of gas-rich, optically-thin disks around intermediate-mass stars at a late stage of the disk evolution. The evolved system is modeled as those where radiation force is sufficiently strong to continuously blow out small grains ($\lesssim 4 {\,\rm μm}$), which are an essential component for driving photoevaporation via photoelectric heating induced by stellar far-ultraviolet (FUV). We find that the grain depletion reduces photoelectric heating, so that FUV photoevaporation is not excited. Extreme-ultraviolet (EUV) photoevaporation is dominant and yields a mass-loss rate of $2$--$5\times10^{-10}(Φ_{\rm EUV}/10^{41}{\,\rm s}^{-1})^{1/2}\,M_\odot\,{\rm yr}^{-1}$, where $Φ_{\rm EUV}$ is the EUV emission rate. The estimated lifetimes of the gas component are $\sim 50 (M_{\rm disk}/10^{-2}\,M_\odot)(Φ_{\rm EUV}/10^{41}\,{\rm s}^{-1})^{1/2}\,{\rm Myr}$ and depend on the ``initial'' disk mass at the point small grains have been depleted in the system. With an order estimation, we show that the gas component can survive for a much longer time around A-type stars than lower-mass stars. This trend is consistent with the higher frequency of gas-rich debris disks around A-type stars, implying the possibility of the gas component being protoplanetary remnants. △ Less

Submitted 24 May, 2021; v1 submitted 14 September, 2020; originally announced September 2020.

Comments: Accepted for publication in ApJ

Modeling disk fragmentation and multiplicity in massive star formation

Authors: G. André Oliva, Rolf Kuiper

Abstract: We investigate the formation and early evolution and fragmentation of an accretion disk around a forming massive protostar. We use a grid-based self-gravity-radiation-hydrodynamics code including a sub-grid module for stellar and dust evolution. On purpose, we do not use sink particles to allow for all paths of fragment formation and destruction, but instead keeping the spatial grid resolution hig… ▽ More We investigate the formation and early evolution and fragmentation of an accretion disk around a forming massive protostar. We use a grid-based self-gravity-radiation-hydrodynamics code including a sub-grid module for stellar and dust evolution. On purpose, we do not use sink particles to allow for all paths of fragment formation and destruction, but instead keeping the spatial grid resolution high enough to properly resolve the physical length scales of the problem. We use a 3D grid in spherical coordinates with a logarithmic scaling in the radial direction and cosine scaling in the polar direction. Because of that, roughly 25% of the total number of grid cells, corresponding to $\sim$ 26 million grid cells, are used to model the disk physics. They constitute the highest resolution simulations performed up to now on disk fragmentation around a forming massive star with the physics considered here. We study the convergence of our results by performing the same simulation for 5 different resolutions. We start from the collapse of a molecular cloud; a massive (proto)star is formed in its center, surrounded by a fragmenting Keplerian-like accretion disk with spiral arms. The fragments have masses of $\sim 1 M_\odot$, and their continuous interactions with the disk, spiral arms and other fragments results in eccentric orbits. Fragments form hydrostatic cores, surrounded by secondary disks with spiral arms that also produce new fragments. We identified several mechanisms of fragment formation, interaction and destruction. Central temperatures of the fragments can reach the hydrogen dissociation limit, form second Larson cores and evolve into companion stars. Based on this, we study the multiplicity predicted by the simulations and find $\sim 6$ companions at different distances from the primary: from possible spectroscopic multiples, to companions at distances between 1000 and 2000 au. △ Less

Submitted 31 August, 2020; originally announced August 2020.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 644, A41 (2020)

Multi-wavelength modelling of the circumstellar environment of the massive proto-star AFGL 2591 VLA 3

Authors: F. A. Olguin, M. G. Hoare, K. G. Johnston, F. Motte, H. -R. V. Chen, H. Beuther, J. C. Mottram, A. Ahmadi, C. Gieser, D. Semenov, T. Peters, A. Palau, P. D. Klaassen, R. Kuiper, Á. Sánchez-Monge, Th. Henning

Abstract: We have studied the dust density, temperature and velocity distributions of the archetypal massive young stellar object (MYSO) AFGL 2591. Given its high luminosity ($L=2 \times 10^5$ L$_\odot$) and distance ($d=3.3$ kpc), AFGL 2591 has one of the highest $\sqrt{L}/d$ ratio, giving better resolved dust emission than any other MYSO. As such, this paper provides a template on how to use resolved mult… ▽ More We have studied the dust density, temperature and velocity distributions of the archetypal massive young stellar object (MYSO) AFGL 2591. Given its high luminosity ($L=2 \times 10^5$ L$_\odot$) and distance ($d=3.3$ kpc), AFGL 2591 has one of the highest $\sqrt{L}/d$ ratio, giving better resolved dust emission than any other MYSO. As such, this paper provides a template on how to use resolved multi-wavelength data and radiative transfer to obtain a well-constrained 2-D axi-symmetric analytic rotating infall model. We show for the first time that the resolved dust continuum emission from Herschel 70 $μ$m observations is extended along the outflow direction, whose origin is explained in part from warm dust in the outflow cavity walls. However, the model can only explain the kinematic features from CH$_3$CN observations with unrealistically low stellar masses ($<15$ M$_\odot$), indicating that additional physical processes may be playing a role in slowing down the envelope rotation. As part of our 3-step continuum and line fitting, we have identified model parameters that can be further constrained by specific observations. High-resolution mm visibilities were fitted to obtain the disc mass (6 M$_\odot$) and radius (2200 au). A combination of SED and near-IR observations were used to estimate the luminosity and envelope mass together with the outflow cavity inclination and opening angles. △ Less

Submitted 12 August, 2020; originally announced August 2020.

Comments: 34 pages, 22 figures, 15 tables, 4 appendices. Accepted for publication in MNRAS

The Dynamic Proto-atmospheres around Low-Mass Planets with Eccentric Orbits

Authors: Chuhong Mai, Steven J. Desch, Rolf Kuiper, Gabriel-Dominique Marleau, Cornelis Dullemond

Abstract: Protoplanets are able to accrete primordial atmospheres when embedded in the gaseous protoplanetary disk. The formation and structure of the proto-atmosphere are subject to the planet--disk environment and orbital effects. Especially, when planets are on eccentric orbits, their relative velocities to the gas can exceed the sound speed. The planets generate atmosphere-stripping bow shocks. We inves… ▽ More Protoplanets are able to accrete primordial atmospheres when embedded in the gaseous protoplanetary disk. The formation and structure of the proto-atmosphere are subject to the planet--disk environment and orbital effects. Especially, when planets are on eccentric orbits, their relative velocities to the gas can exceed the sound speed. The planets generate atmosphere-stripping bow shocks. We investigate the proto-atmospheres on low-mass planets with eccentric orbits with radiation-hydrodynamics simulations. A 2D radiative model of the proto-atmosphere is established with tabulated opacities for the gas and dust. The solutions reveal large-scale gas recycling inside a bow shock structure. The atmospheres on eccentric planets are typically three to four orders of magnitude less massive than those of planets with circular orbits. Overall, however, a supersonic environment is favorable for planets to keep an early stable atmosphere, rather than harmful, due to the steady gas supply through the recycling flow. We also quantitatively explore how such atmospheres are affected by the relative velocity of the planet to the gas, the planet mass, and the background gas density. Our time-dependent simulations track the orbital evolution of the proto-atmosphere with the planet--disk parameters changing throughout the orbit. Atmospheric properties show oscillatory patterns as the planet travels on an eccentric orbit, with a lag in phase. To sum up, low-mass eccentric planets can retain small proto-atmospheres despite the stripping effects of bow shocks. The atmospheres are always connected to and interacting with the disk gas. These findings provide important insights into the impacts of migration and scattering on planetary proto-atmospheres. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Comments: Accepted for publication in ApJ (23 pages, 15 figures, 2 tables)

FRB 200428: an Impact between an Asteroid and a Magnetar

Authors: Jin-Jun Geng, Bing Li, Long-Biao Li, Shao-Lin Xiong, Rolf Kuiper, Yong-Feng Huang

Abstract: A fast radio burst (FRB) was recently detected to be associated with a hard X-ray burst from the Galactic magnetar SGR 1935+2154. Scenarios involving magnetars for FRBs are hence highly favored. In this work, we suggest that the impact between an asteroid and a magnetar could explain such a detection. According to our calculations, an asteroid of mass $10^{20}$ g will be disrupted at a distance of… ▽ More A fast radio burst (FRB) was recently detected to be associated with a hard X-ray burst from the Galactic magnetar SGR 1935+2154. Scenarios involving magnetars for FRBs are hence highly favored. In this work, we suggest that the impact between an asteroid and a magnetar could explain such a detection. According to our calculations, an asteroid of mass $10^{20}$ g will be disrupted at a distance of $7 \times 10^9$ cm when approaching the magnetar. The accreted material will flow along the magnetic field lines from the Alfvén radius $\sim 10^7$ cm. After falling onto the magnetar's surface, an instant accretion column will be formed, producing a Comptonized X-ray burst and an FRB in the magnetosphere. We show that all the observational features of FRB 200428 could be interpreted self-consistently in this scenario. We predict quasi-periodic oscillations in this specific X-ray burst, which can serve as an independent observational test. △ Less

Submitted 22 July, 2020; v1 submitted 8 June, 2020; originally announced June 2020.

Comments: 5 pages, 1 figure. References updated, accepted by ApJL

A Detailed View of the Circumstellar Environment and Disk of the Forming O-star AFGL 4176

Authors: K. G. Johnston, M. G. Hoare, H. Beuther, H. Linz, P. Boley, R. Kuiper, N. D. Kee, T. P. Robitaille

Abstract: We present a detailed analysis of the disk and circumstellar environment of the forming O-type star AFGL 4176 mm1, placing results from the Atacama Large Millimeter/submillimeter Array (ALMA) into context with multiwavelength data. With ALMA, we detect seventeen 1.2 mm continuum sources within 5$''$ (21,000 au) of AFGL 4176 mm1. We find that mm1 has a spectral index of 3.4$\pm$0.2 across the ALMA… ▽ More We present a detailed analysis of the disk and circumstellar environment of the forming O-type star AFGL 4176 mm1, placing results from the Atacama Large Millimeter/submillimeter Array (ALMA) into context with multiwavelength data. With ALMA, we detect seventeen 1.2 mm continuum sources within 5$''$ (21,000 au) of AFGL 4176 mm1. We find that mm1 has a spectral index of 3.4$\pm$0.2 across the ALMA band, with $>$87% of its 1.2 mm continuum emission from dust. The source mm2, projected 4200 au from mm1, may be a companion or a blueshifted knot in a jet. We also explore the morphological differences between the molecular lines detected with ALMA, finding 203 lines from 25 molecules, which we categorize into several morphological types. Our results show that AFGL 4176 mm1 provides an example of a forming O-star with a large and chemically complex disk, which is mainly traced by nitrogen-bearing molecules. Lines that show strong emission on the blueshifted side of the disk are predominantly oxygen-bearing, which we suggest are tracing a disk accretion shock. The molecules C$^{34}$S, H$_2$CS and CH$_{3}$CCN trace a slow wide-angle wind or dense structures in the outflow cavity walls. With the Australia Telescope Compact Array (ATCA), we detect a compact continuum source ($<$2000 $\times$ 760 au) at 1.2 cm, associated with mm1, of which $>$96% is from ionized gas. The ATCA NH$_3$(1,1) and (2,2) emission traces a large-scale (r$\sim$0.5 pc) rotating toroid with the disk source mm1 in the blueshifted part of this structure offset to the northwest. △ Less

Submitted 23 May, 2020; v1 submitted 28 April, 2020; originally announced April 2020.

Comments: Accepted for publication in The Astrophysical Journal, 50 pages, 26 figures. Data available online at http://doi.org/10.5281/zenodo.3369188

Birth of convective low-mass to high-mass second Larson cores

Authors: Asmita Bhandare, Rolf Kuiper, Thomas Henning, Christian Fendt, Mario Flock, Gabriel-Dominique Marleau

Abstract: Stars form as an end product of the gravitational collapse of cold, dense gas in magnetized molecular clouds. This multi-scale scenario occurs via the formation of two quasi-hydrostatic cores and involves complex physical processes, which require a robust, self-consistent numerical treatment. The aim of this study is to understand the formation and evolution of the second Larson core and the depen… ▽ More Stars form as an end product of the gravitational collapse of cold, dense gas in magnetized molecular clouds. This multi-scale scenario occurs via the formation of two quasi-hydrostatic cores and involves complex physical processes, which require a robust, self-consistent numerical treatment. The aim of this study is to understand the formation and evolution of the second Larson core and the dependence of its properties on the initial cloud core mass. We used the PLUTO code to perform high resolution, 1D and 2D RHD collapse simulations. We include self-gravity and use a grey FLD approximation for the radiative transfer. Additionally, we use for the gas EOS density- and temperature-dependent thermodynamic quantities to account for the effects such as dissociation, ionisation, and molecular vibrations and rotations. Properties of the second core are investigated using 1D studies spanning a wide range of initial cloud core masses from 0.5 to 100 $M_{\odot}$. Furthermore, we expand to 2D collapse simulations for a few cases of 1, 5, 10, and 20 $M_{\odot}$. We follow the evolution of the second core for $\geq$ 100 years after its formation, for each of these non-rotating cases. Our results indicate a dependence of several second core properties on the initial cloud core mass. For the first time, due to an unprecedented resolution, our 2D non-rotating collapse studies indicate that convection is generated in the outer layers of the second core, which is formed due to the gravitational collapse of a 1 $M_{\odot}$ cloud core. Additionally, we find large-scale oscillations of the second accretion shock front triggered by the standing accretion shock instability, which has not been seen before in early evolutionary stages of stars. We predict that the physics within the second core would not be significantly influenced by the effects of magnetic fields or an initial cloud rotation. △ Less

Submitted 26 April, 2020; originally announced April 2020.

Comments: 22 pages, 20 figures, and 2 tables. Accepted for publication in A&A. Collapse movie: https://keeper.mpdl.mpg.de/f/f04abdeabdf3472fb56d/

Journal ref: A&A 638, A86 (2020)

Formation of massive stars under protostellar radiation feedback: Very metal-poor stars

Authors: Hajime Fukushima, Takashi Hosokawa, Gen Chiaki, Kazuyuki Omukai, Naoki Yoshida, Rolf Kuiper

Abstract: We study the formation of very metal-poor stars under protostellar radiative feedback effect. We use cosmological simulations to identify low-mass dark matter halos and star-forming gas clouds within them. We then follow protostar formation and the subsequent long-term mass accretion phase of over one million years using two-dimensional radiation-hydrodynamics simulations. We show that the critica… ▽ More We study the formation of very metal-poor stars under protostellar radiative feedback effect. We use cosmological simulations to identify low-mass dark matter halos and star-forming gas clouds within them. We then follow protostar formation and the subsequent long-term mass accretion phase of over one million years using two-dimensional radiation-hydrodynamics simulations. We show that the critical physical process that sets the final mass is formation and expansion of a bipolar HII region. The process is similar to the formation of massive primordial stars, but radiation pressure exerted on dust grains also contributes to halting the accretion flow in the low-metallicity case. We find that the net feedback effect in the case with metallicity $Z = 10^{-2}~Z_{\odot}$ is stronger than in the case with $Z \sim 1~Z_{\odot}$. With decreasing metallicity, the radiation pressure effect becomes weaker, but photoionization heating of the circumstellar gas is more efficient owing to the reduced dust attenuation. In the case with $Z = 10^{-2}~Z_{\odot}$, the central star grows as massive as 200 solar-masses, similarly to the case of primordial star formation. We conclude that metal-poor stars with a few hundred solar masses can be formed by gas accretion despite the strong radiative feedback. △ Less

Submitted 22 April, 2021; v1 submitted 5 April, 2020; originally announced April 2020.

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