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.

A fast radio burst 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 $7 \times 10^9$ cm when approaching the magnetar, the accreted material will flow along the magnetic field lines at Alfv\'en radius $\sim 10^7$ cm. After falling onto the magnetar 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 could be interpreted self-consistently in our scenario. Quasi-periodic oscillations in this specific X-ray burst may help to verify our scenario.

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.

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.

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.

Recent observations have suggested that circumstellar disks may commonly form around young stellar objects. Although the formation of circumstellar disks can be a natural result of the conservation of angular momentum in the parent cloud, theoretical studies instead show disk formation to be difficult from dense molecular cores magnetized to a realistic level, owing to efficient magnetic braking that transports a large fraction of the angular momentum away from the circumstellar region. We review recent progress in the formation and early evolution of disks around young stellar objects of both low-mass and high-mass, with an emphasis on mechanisms that may bridge the gap between observation and theory, including non-ideal MHD effects and asymmetric perturbations in the collapsing core (e.g., magnetic field misalignment and turbulence). We also address the associated processes of outflow launching and the formation of multiple systems, and discuss possible implications in properties of protoplanetary disks.

We investigate the turbulence driving mode of ionizing radiation from massive stars on the surrounding interstellar medium (ISM). We run hydrodynamical simulations of a turbulent cloud impinged by a plane-parallel ionization front. We find that the ionizing radiation forms pillars of neutral gas reminiscent of those seen in observations. We quantify the driving mode of the turbulence in the neutral gas by calculating the driving parameter $b$, which is characterised by the relation $\sigma_s^2 = \ln({1+b^2\mathcal{M}^2})$ between the variance of the logarithmic density contrast $\sigma_s^2$ (where $s = \ln({\rho/\rho_0})$ with the gas density $\rho$ and its average $\rho_0$), and the turbulent Mach number $\mathcal{M}$. Previous works have shown that $b\sim1/3$ indicates solenoidal (divergence-free) driving and $b\sim1$ indicates compressive (curl-free) driving, with $b\sim1$ producing up to ten times higher star formation rates than $b\sim1/3$. The time variation of $b$ in our study allows us to infer that ionizing radiation is inherently a compressive turbulence driving source, with a time-averaged $b\sim 0.76 \pm 0.08$. We also investigate the value of $b$ of the pillars, where star formation is expected to occur, and find that the pillars are characterised by a natural mixture of both solenoidal and compressive turbulent modes ($b\sim0.4$) when they form, and later evolve into a more compressive turbulent state with $b\sim0.5$--$0.6$. A virial parameter analysis of the pillar regions supports this conclusion. This indicates that ionizing radiation from massive stars may be able to trigger star formation by producing predominately compressive turbulent gas in the pillars.

Coalescence of intermediate-mass black holes (IMBHs) as a result of the migration toward galactic centers via dynamical friction may contribute to the formation of supermassive BHs. Here we reinvestigate the gaseous dynamical friction, which was claimed to be inefficient with radiative feedback from BHs in literature, by performing 3D radiation-hydrodynamics simulations that solve the flow structure in the vicinity of BHs. We consider a $10^4~M_\odot$ BH moving at the velocity $V_{\rm flow}$ through the homogeneous medium with metallicity $Z$ in the range of $0-0.1~Z_\odot$ and density $n_{\infty}$. We show that, if $n_{\infty} \lesssim 10^{6}~{\rm cm^{-3}}$ and $V_{\rm flow} \lesssim 60~{\rm km~s^{-1}}$, the BH is accelerated forward because of the gravitational pull from a dense shell ahead of an ionized bubble around the BH, regardless of the value of $Z$. If $n_{\infty} \gtrsim 10^{6}~{\rm cm^{-3}}$, however, our simulation shows the opposite result. The ionized bubble and associating shell temporarily appear, but immediately go downstream with significant ram pressure of the flow. They eventually converge into a massive downstream wake, which gravitationally drags the BH backward. The BH decelerates over the timescale of $\sim 0.01$~Myr, much shorter than the dynamical timescale in galactic disks. Our results suggest that IMBHs that encounter the dense clouds rapidly migrate toward galactic centers, where they possibly coalescence with others.

We present high-resolution (30 mas or 130 au at 4.2 kpc) Atacama Large Millimeter/submillimeter Array observations at 1.2 mm of the disc around the forming O-type star AFGL 4176 mm1. The disc (AFGL 4176 mm1-main) has a radius of ~1000 au and contains significant structure, most notably a spiral arm on its redshifted side. We fitted the observed spiral with logarithmic and Archimedean spiral models. We find that both models can describe its structure, but the Archimedean spiral with a varying pitch angle fits its morphology marginally better. As well as signatures of rotation across the disc, we observe gas arcs in CH$_3$CN that connect to other millimetre continuum sources in the field, supporting the picture of interactions within a small cluster around AFGL 4176 mm1-main. Using local thermodynamic equilibrium modelling of the CH$_3$CN K-ladder, we determine the temperature and velocity field across the disc, and thus produce a map of the Toomre stability parameter. Our results indicate that the outer disc is gravitationally unstable and has already fragmented or is likely to fragment in the future, possibly producing further companions. These observations provide evidence that disc fragmentation is one possible pathway towards explaining the high fraction of multiple systems around high-mass stars.

In the disk-mediated accretion scenario for the formation of the most massive stars, gravitational instabilities in the disk can force it to fragment. We investigate the effects of inclination and spatial resolution on observable kinematics and stability of disks in high-mass star formation. We study a high-resolution 3D radiation-hydrodynamic simulation that leads to the fragmentation of a massive disk. Using RADMC-3D we produce 1.3 mm continuum and CH3CN line cubes at different inclinations. The model is set to different distances and synthetic observations are created for ALMA at ~80 mas resolution and NOEMA at ~0.3''. The synthetic ALMA observations resolve all fragments and their kinematics well. The synthetic NOEMA observations at 800 pc (~300 au resolution) are able to resolve the fragments, while at 2000 pc (~800 au resolution) only a single slightly elongated structure is observed. The position-velocity (PV) plots show the differential rotation of material best in the edge-on views. As the observations become less resolved, the inner high-velocity components of the disk become blended with the envelope and the PV plots resemble rigid-body-like rotation. Protostellar mass estimates from PV plots of poorly resolved observations are therefore overestimated. We fit the emission of CH3CN lines and produce maps of gas temperature with values in the range of 100-300 K. Studying the Toomre stability of the disks in the resolved observations, we find Q values below the critical value for stability against gravitational collapse at the positions of the fragments and the arms connecting the fragments. For the poorly resolved observations we find low Q values in the outskirts of the disk. Therefore we are able to predict that the disk is unstable and fragmenting even in poorly resolved observations. This conclusion is true regardless of knowledge about the inclination of the disk.

We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth.

The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions.

The super-Eddington accretion onto intermediate seed BHs is a potential formation mode of supermassive black holes exceeding $10^9~M_\odot$ in the early universe. We here investigate how such rapid accretion may occur with finite amounts of heavy elements contained in the gas and dust. In our 1D radiation-hydrodynamics simulations, the radiative transfer is solved for both the direct UV lights emitted by an accretion disk and the diffuse IR lights thermally emitted by dust grains. Our results show that the radiative force by the IR lights causes a strong feedback to regulate the mass accretion. The resulting mean accretion rate is lower with the higher metallicity, and there is the critical metallicity $Z \sim 10^{-2}~Z_\odot$, above which the super-Eddington accretion is prevented by the radiation pressure of the IR lights. With this taken into account, we examine if the dusty super-Eddington accretion occurs in young galaxies using a simple model. We show that a sufficient number of galaxies at $z \gtrsim 10$ can be such potential sites if BHs accrete the cold dense gas with $T \sim 10^2$ K, approximately the thermal equilibrium value at $Z = 10^{-2}~Z_\odot$. We argue that the efficiency of the BH growth via the rapid accretion depends on the metallicity, and that the metallicity slightly lower than $10^{-2}~Z_\odot$ provides a chance for the most efficient growth.

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

In the core-accretion formation scenario of gas giants, most of the gas accreting onto a planet is processed through an accretion shock. In this series of papers we study this shock since it is key in setting the forming planet's structure and thus its post-formation luminosity, with dramatic observational consequences. We perform one-dimensional grey radiation-hydrodynamical simulations with non-equilibrium (two-temperature) radiation transport and up-to-date opacities. We survey the parameter space of accretion rate, planet mass, and planet radius and obtain post-shock temperatures, pressures, and entropies, as well as global radiation efficiencies. We find that usually, the shock temperature T_shock is given by the "free-streaming" limit. At low temperatures the dust opacity can make the shock hotter but not significantly. We corroborate this with an original semi-analytical derivation of T_shock . We also estimate the change in luminosity between the shock and the nebula. Neither T_shock nor the luminosity profile depend directly on the optical depth between the shock and the nebula. Rather, T_shock depends on the immediate pre-shock opacity, and the luminosity change on the equation of state (EOS). We find quite high immediate post-shock entropies (S ~ 13-20 kB/mH), which makes it seem unlikely that the shock can cool the planet. The global radiation efficiencies are high (eta^phys > 97%) but the remainder of the total incoming energy, which is brought into the planet, exceeds the internal luminosity of classical cold starts by orders of magnitude. Overall, these findings suggest that warm or hot starts are more plausible.

This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars.

We perform a series of relativistic magnetohydrodynamics simulations to investigate how a hot magnetic jet propagates within the dynamical ejecta of a binary neutron star merger, with the focus on how the jet structure depends on the delay time of jet launching with respect to the merger time, $\Delta t_{\rm jet}$. We find that regardless of the jet launching delay time, a structured jet with an angle-dependent luminosity and Lorentz factor is always formed after the jet breaks out the ejecta. On the other hand, the jet launching delay time has an impact on the jet structure. If the jet launching delay time is relatively long, e.g., $\ge 0.5$ s, the line-of-sight material has a dominant contribution from the cocoon. On the other hand, for a relatively short jet launching delay time, the jet penetrates through the ejecta early on and develops an angular structure afterward. The line-of-sight ejecta is dominated by the structured jet itself. We discuss the case of GW170817/GRB 170817A within the framework of both long and short $\Delta t_{\rm jet}$. Future more observations of GW/GRB associations can help to differentiate between these two scenarios.

We study the properties of rotating accretion flows onto supermassive black holes (SMBHs) using axisymmetric two-dimensional hydrodynamical simulations with radiative cooling and BH feedback. The simulations resolve the accretion dynamics of gas outside from the BH influence radius through an inner accretion disk. For lower Bondi accretion rates in units of the Eddington rate ($\dot{M}_{\rm B}\ll 10^{-3}~\dot{M}_{\rm Edd}$), the BH feeding is suppressed due to turbulent motion by several orders of magnitudes from the Bondi rate. Thus, the radiative luminosity results in as low as $\sim 10^{-10}-10^{-7}~L_{\rm Edd}$, where $L_{\rm Edd}$ is the Eddington luminosity. For higher rates of $\dot{M}_{\rm B}> 10^{-3}~\dot{M}_{\rm Edd}$, the optically-thin accreting gas cools via free-free emission and forms a geometrically-thin disk, which feeds the BH efficiently and increases the radiative luminosity to $> 10^{-3}~L_{\rm Edd}$. The transitional behavior of accreting BHs in galactic nuclei from radiatively inefficient phases to cold disk accretion naturally explains (1) the reason for the offset between the observed luminosities and theoretical predictions for nearby quiescent SMBHs, and (2) the conditions to fuel gas into the nuclear SMBH. In addition, the cold disk formed in galactic nuclei tends to be gravitationally unstable and leads to star formation when the Bondi rate is as high as $ \dot{M}_{\rm B} > 10^{-2}~M_\odot~{\rm yr}^{-1}$. This is a plausible explanation of the correlation observed between star formation rates and BH feeding rates in Seyfert galaxies.

It is well established that Solar-mass stars gain mass via disk accretion, until the mass reservoir of the disk is exhausted and dispersed, or condenses into planetesimals. Accretion disks are intimately coupled with mass ejection via polar cavities, in the form of jets and less collimated winds, which allow mass accretion through the disk by removing a substantial fraction of its angular momentum. Whether disk accretion is the mechanism leading to the formation of stars with much higher masses is still unclear. Here, we are able to build a comprehensive picture for the formation of an O-type star, by directly imaging a molecular disk which rotates and undergoes infall around the central star, and drives a molecular jet which arises from the inner disk regions. The accretion disk is truncated between 2000-3000au, it has a mass of about a tenth of the central star mass, and is infalling towards the central star at a high rate (6x10^-4 Msun/yr), as to build up a very massive object. These findings, obtained with the Atacama Large Millimeter/submillimeter Array at 700au resolution, provide observational proof that young massive stars can form via disk accretion much like Solar-mass stars.

Radiative feedback from luminous, massive stars during their formation is a key process in moderating accretion onto the stellar object. In the prior papers in this series, we showed that one form such feedback takes is UV line-driven disc ablation. Extending on this study, we now constrain the strength of this effect in the parameter range of star and disc properties appropriate to forming massive stars. Simulations show that ablation rate depends strongly on stellar parameters, but that this dependence can be parameterized as a nearly constant, fixed enhancement over the wind mass loss rate, allowing us to predict the rate of disc ablation for massive (proto)stars as a function of stellar mass and metallicity. By comparing this to predicted accretion rates, we conclude that ablation is a strong feedback effect for very massive (proto)stars which should be considered in future studies of massive star formation.

Methods: We observed the high-mass hot core region G351.77-0.54 with ALMA and more than 16km baselines. Results: At a spatial resolution of 18/40au (depending on the distance), we identify twelve sub-structures within the inner few thousand au of the region. The brightness temperatures are high, reaching values greater 1000K, signposting high optical depth toward the peak positions. Core separations vary between sub-100au to several 100 and 1000au. The core separations and approximate masses are largely consistent with thermal Jeans fragmentation of a dense gas core. Due to the high continuum optical depth, most spectral lines are seen in absorption. However, a few exceptional emission lines are found that most likely stem from transitions with excitation conditions above1000K. Toward the main continuum source, these emission lines exhibit a velocity gradient across scales of 100-200au aligned with the molecular outflow and perpendicular to the previously inferred disk orientation. While we cannot exclude that these observational features stem from an inner hot accretion disk, the alignment with the outflow rather suggests that it stems from the inner jet and outflow region. The highest-velocity features are found toward the peak position, and no Hubble-like velocity structure can be identified. Therefore, these data are consistent with steady-state turbulent entrainment of the hot molecular gas via Kelvin-Helmholtz instabilities at the interface between the jet and the outflow. Conclusions: Resolving this high-mass star-forming region at sub-50au scales indicates that the hierarchical fragmentation process in the framework of thermal Jeans fragmentation can continue down to the smallest accessible spatial scales. Velocity gradients on these small scales have to be treated cautiously and do not necessarily stem from disks, but may be better explained with outflow emission.

In a comprehensive convergence study, we investigate the computational conditions necessary to resolve disk formation and jet-launching processes, and analyze possible caveats. We explore the magneto-hydrodynamic (MHD) processes of the collapse of massive prestellar cores in detail, including an analysis of the forces involved and their temporal evolution for up to two free-fall times. We conduct MHD simulations, combining nonideal MHD, self-gravity, and very high resolutions as they have never been achieved before. Our setup includes a 100 Msol cloud core that collapses under its own self-gravity to self-consistently form a dense disk structure and launch tightly collimated magneto-centrifugal jets and wide-angle winds. Our high-resolution simulations can resolve a magneto-centrifugal jet and a magnetic pressure-driven outflow, separately. The nature of the outflows depends critically on spatial resolution. Only high-resolution simulations are able to differentiate a magneto-centrifugally launched, highly collimated jet from a slow wide-angle magnetic-pressure-driven tower flow. Of these two outflow components, the tower flow dominates angular-momentum transport. The mass outflow rate is dominated by the entrained material from the interaction of the jet with the stellar environment and only part of the ejected medium is directly launched from the accretion disk. A tower flow can only develop to its full extent when much of the original envelope has already dispersed. Taking into account both the mass launched from the surface of the disk and the entrained material from the envelope, we find an ejection-to-accretion efficiency of 10%. Nonideal MHD is required to form centrifugally supported accretion disks and the disk size is strongly dependent on spatial resolution.

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources.

We study low-density axisymmetric accretion flows onto black holes (BHs) with two-dimensional hydrodynamical simulations, adopting the $\alpha$-viscosity prescription. When the gas angular momentum is low enough to form a rotationally supported disk within the Bondi radius ($R_{\rm B}$), we find a global steady accretion solution. The solution consists of a rotational equilibrium distribution at $r\sim R_{\rm B}$, where the density follows $\rho \propto (1+R_{\rm B}/r)^{3/2}$, surrounding a geometrically thick and optically thin accretion disk at the centrifugal radius, where thermal energy generated by viscosity is transported via strong convection. Physical properties of the inner solution agree with those expected in convection-dominated accretion flows (CDAF; $\rho \propto r^{-1/2}$). In the inner CDAF solution, the gas inflow rate decreases towards the center due to convection ($\dot{M}\propto r$), and the net accretion rate (including both inflows and outflows) is strongly suppressed by several orders of magnitude from the Bondi accretion rate $\dot{M}_{\rm B}$ The net accretion rate depends on the viscous strength, following $\dot{M}/\dot{M}_{\rm B}\propto (\alpha/0.01)^{0.6}$. This solution holds for low accretion rates of $\dot{M}_{\rm B}/\dot{M}_{\rm Edd}< 10^{-3}$ having minimal radiation cooling, where $\dot{M}_{\rm Edd}$ is the Eddington rate. In a hot plasma at the bottom ($r<10^{-3}~R_{\rm B}$), thermal conduction would dominate the convective energy flux. Since suppression of the accretion by convection ceases, the final BH feeding rate is found to be $\dot{M}/\dot{M}_{\rm B} \sim 10^{-3}-10^{-2}$. This rate is as low as $\dot{M}/\dot{M}_{\rm Edd} \sim 10^{-7}-10^{-6}$ inferred for SgrA$^*$ and the nuclear BHs in M31 and M87, and can explain the low luminosities in these sources, without invoking any feedback mechanism.

Orbital mean motion resonances in planetary systems originate from dissipative processes in disk-planet interactions that lead to orbital migration. In multi-planet systems that host giant planets, the perturbation of the protoplanetary disk strongly affects the migration of companion planets. By studying the well-characterized resonant planetary system around GJ 876 we aim to explore which effects shape disk-driven migration in such a multi-planet system to form resonant chains. We modelled the orbital migration of three planets embedded in a protoplanetary disk using two-dimensional locally isothermal hydrodynamical simulations. We performed a parameter study by varying the disk thickness, $\alpha$ viscosity, mass as well as the initial position of the planets. Moreover, we have analysed and compared simulations with various boundary conditions at the disk's inner rim. We find that due to the high masses of the giant planets, substantial eccentricity can be excited in the disk. This results in large variations of the torque acting on the outer lower-mass planet, which we attribute to a shift of Lindblad and corotation resonances due to disk eccentricity. Depending on disk parameters, the migration of the outer planet can be stopped in a non-resonant state. In other models, the outer planet is able to open a partial gap and to circularize the disk again, later entering a 2:1 resonance with the most massive planet in the system to complete the observed Laplace resonance. Disk-mediated interactions between planets due to spiral waves and excitation of disk eccentricity cause deviations from smooth inward migration of exterior lower mass planets. Self-consistent modelling of the disk-driven migration of multi-planet systems is thus mandatory. Our results are compatible with a late migration of the outermost planet into the resonant chain, when the giant planet pair already is in resonance.

We perform a suite of radiation hydrodynamics simulations of photoevaporating disks with varying the metallicity in a wide range of $10^{-3} \, Z_\odot \leq Z \leq 10^{0.5} \, Z_\odot $. We follow the disk evolution for over $\sim 5000$ years by solving hydrodynamics, radiative transfer, and non-equilibrium chemistry. Our chemistry model is updated from the first paper of this series by adding X-ray ionization and heating. We study the metallicity dependence of the disk photoevaporation rate and examine the importance of X-ray radiation. In the fiducial case with solar metallicity, including the X-ray effects does not significantly increase the photoevaporation rate when compared to the case with ultra-violet (UV) radiation only. At sub-solar metallicities in the range of $Z \gtrsim 10^{-1.5} \, Z_\odot $, the photoevaporation rate increases as metallicity decreases owing to the reduced opacity of the disk medium. The result is consistent with the observational trend that disk lifetimes are shorter in low metallicity environments. Contrastingly, the photoevaporation rate decreases at even lower metallicities of $Z \lesssim 10^{-1.5} \, Z_\odot $, because dust-gas collisional cooling remains efficient compared to far UV photoelectric heating whose efficiency depends on metallicity. The net cooling in the interior of the disk suppresses the photoevaporation. However, adding X-ray radiation significantly increases the photoevaporation rate, especially at $Z \sim 10^{-2}\, Z_\odot$. Although the X-ray radiation itself does not drive strong photoevaporative flows, X-rays penetrate deep into the neutral region in the disk, increase the ionization degree there, and reduce positive charges of grains. Consequently, the effect of photoelectric heating by far UV radiation is strengthened by the X-rays and enhances the disk photoevaporation.

The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales.

In this study, the main goal is to understand the molecular cloud core collapse through the stages of first and second hydrostatic core formation. We investigate the properties of Larson's first and second cores following the evolution of the molecular cloud core until formation of Larson's cores. We expand these collapse studies for the first time to span a wide range of initial cloud masses from 0.5 to 100 Msun. Understanding the complexity of the numerous physical processes involved in the very early stages of star formation requires detailed thermodynamical modeling in terms of radiation transport and phase transitions. For this we use a realistic gas equation of state via a density and temperature-dependent adiabatic index and mean molecular weight to model the phase transitions. We use a gray treatment of radiative transfer coupled with hydrodynamics to simulate Larson's collapse in spherical symmetry. We reveal a dependence of a variety of first core properties on the initial cloud mass. The first core radius and mass increase from the low-mass to the intermediate-mass regime and decrease from the intermediate-mass to the high-mass regime. The lifetime of first cores strongly decreases towards the intermediate- and high-mass regime. Our studies show the presence of a transition region in the intermediate-mass regime. Low-mass protostars tend to evolve through two distinct stages of formation which are related to the first and second hydrostatic cores. In contrast, in the high-mass star formation regime, the collapsing cloud cores rapidly evolve through the first collapse phase and essentially immediately form Larson's second cores.

In studying the formation of massive stars, it is essential to consider the strong radiative feedback on the stars' natal environments from their high luminosities ($10^4 \sim 10^6 L_\odot$). Given that massive stars contract to main-sequence-like radii before accretion finishes, one form this feedback takes is UV line-acceleration, resulting in outflows much like those expected from main-sequence massive stars. As shown by the prior papers in this series, in addition to driving stellar winds, such line forces also ablate the surface layers off of circumstellar discs within a few stellar radii of the stellar photosphere. This removal of material from an accretion disc in turn results in a decreased accretion rate onto the forming star. Quantifying this, however, requires accounting for the continuum optical depth of the disc along the non-radial rays required for the three-dimensional line-acceleration prescription used in this paper series. We introduce the "thin disc approximation", allowing these continuum optical depths arising from an optically thick but geometrically thin disc to be dynamically treated in the context of radiation-hydrodynamics simulations. Using this approximation in full dynamical simulations, we show that such continuum optical depth effects only reduce the disc ablation by 30 percent or less relative to previous simulations that ignored continuum absorption.

As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation.

We present the first simulations of the formation and feedback of massive stars which account for radiation forces as well as photoionization feedback (along with protostellar outflows). In two different accretion scenarios modeled, we determine the relative strength of these feedback components and derive the size of the reservoir from which the forming stars gained their masses. We find that protostellar outflows alone limit the stellar mass growth only in an accretion scenario with a finite mass reservoir; when including accretion and ram pressure from large scales (> 0.1 pc), protostellar outflows do not limit stellar mass growth at all. Photoionization and HII regions dominate the feedback ladder only at later times and only on large scales. Specifically, photoionization yields a broadening of the bipolar outflow cavities and a reduction of the gravitational infall momentum by about 50 percent, but does not limit the stellar mass accretion. On the other hand, we find radiation forces restrain the gravitational infall toward the circumstellar disk, impact the gravito-centrifugal equilibrium at the outer edge of the disk, and eventually shut down stellar accretion completely. The most massive star formed in the simulations accreted 95 solar masses before disk destruction; this mass was drawn-in from an accretion reservoir of approximately 240 solar masses and 0.24 pc in radius. In the regime of very massive stars, the final mass of these stars is controlled by their own radiation force feedback.

Protoplanetary disks are thought to have lifetimes of several million years in the solar neighborhood, but recent observations suggest that the disk lifetimes are shorter in a low metallicity environment. We perform a suite of radiation hydrodynamics simulations of photoevaporation of protoplanetary disks to study the disk structure and its long-term evolution of $\sim 10000$ years, and the metallicity dependence of mass-loss rate. Our simulations follow hydrodynamics, extreme and far ultra-violet radiative transfer, and non-equilibrium chemistry in a self-consistent manner. Dust grain temperatures are also calculated consistently by solving the radiative transfer of the stellar irradiation and grain (re-)emission. We vary the disk gas metallicity over a wide range of $10^{-4}~ Z_\odot \leq Z \leq 10 ~Z_\odot$. The photoevaporation rate is lower with higher metallicity in the range of $10^{-1} \,Z_\odot \lesssim Z \lesssim 10 \,Z_\odot$, because dust shielding effectively prevents far-ultra violet (FUV) photons from penetrating into and heating the dense regions of the disk. The photoevaporation rate sharply declines at even lower metallicities in $10^{-2} \,Z_\odot \lesssim Z \lesssim 10^{-1}\,Z_\odot$, because FUV photoelectric heating becomes less effective than dust-gas collisional cooling. The temperature in the neutral region decreases, and photoevaporative flows are excited only in an outer region of the disk. At $10^{-4}\,Z_\odot \leq Z \lesssim 10^{-2}\,Z_\odot$, HI photoionization heating acts as a dominant gas heating process and drives photoevaporative flows with roughly a constant rate. The typical disk lifetime is shorter at $Z=0.3~Z_\odot$ than at $Z = Z_\odot$, being consistent with recent observations of the extreme outer galaxy.

Context: The importance of magnetic fields at the onset of star formation related to the early fragmentation and collapse processes is largely unexplored today. Aims: We want to understand the magnetic field properties at the earliest evolutionary stages of high-mass star formation. Methods: The Atacama Large Millimeter Array is used at 1.3mm wavelength in full polarization mode to study the polarized emission and by that the magnetic field morphologies and strengths of the high-mass starless region IRDC18310-4. Results: The polarized emission is clearly detected in four sub-cores of the region. In general it shows a smooth distribution, also along elongated cores. Estimating the magnetic field strength via the Davis-Chandrasekhar-Fermi method and following a structure function analysis, we find comparably large magnetic field strengths between ~0.6 and 3.7mG. Comparing the data to spectral line observations, the turbulent-to-magnetic energy ratio is low, indicating that turbulence does not significantly contribute to the stability of the gas clump. A mass-to-flux ratio around the critical value 1.0 - depending on column density - indicates that the region starts to collapse which is consistent with the previous spectral line analysis of the region. Conclusions: While this high-mass region is collapsing and thus at the verge of star formation, the high magnetic field values and the smooth spatial structure indicate that the magnetic field is important for the fragmentation and collapse process. This single case study can only be the starting point for larger sample studies of magnetic fields at the onset of star formation.

Highly collimated parsec-scale jets, generally linked to the presence of an accretion disk, are a commonly observed phenomenon from revealed low-mass young stellar objects. In the past two decades, only a very few of these objects have been directly (or indirectly) observed towards high-mass (M > 8 M$_{\odot}$) young stellar objects, adding to the growing evidence that disk-mediated accretion is a phenomenon also occurring in high-mass stars, the formation mechanism of which is still poorly understood. Of the observed jets from massive young stars, none is in the optical regime (due to these being typically highly obscured by their native material), and none are found outside of the Milky Way. Here, we report the detection of HH 1177, the first extragalactic optical ionized jet originating from a massive young stellar object located in the Large Magellanic Cloud. The jet is highly collimated over the entire measured extent of at least 10 pc, and has a bipolar geometry. The presence of a jet indicates ongoing, disk-mediated accretion, and together with the high degree of collimation, this system is therefore likely to be an up-scaled version of low-mass star formation. We conclude that the physics governing jet launching and collimation is independent of stellar mass.

High-mass stars form in much richer environments than those associated with isolated low-mass stars, and once they reach a certain mass, produce ionised (HII) regions. The formation of these pockets of ionised gas are unique to the formation of high-mass stars (M $>8$ M$_\odot$), and present an excellent opportunity to study the final stages of accretion, which could include accretion through the HII region itself. This study of the dynamics of the gas on both sides of these ionisation boundaries in very young HII regions aims to quantify the relationship between the HII regions and their immediate environments.We present high-resolution ($\sim$ 0.5$"$) ALMA observations of nine HII regions selected from the Red MSX Source (RMS) survey with compact radio emission and bolometric luminosities greater than 10$^4$ L$_\odot$. We focus on the initial presentation of the data, including initial results from the radio recombination line H29$\alpha$, some complementary molecules, and the 256 GHz continuum emission. Of the six (out of nine) regions with H29$\alpha$ detections, two appear to have cometary morphologies with velocity gradients across them, and two appear more spherical with velocity gradients suggestive of infalling ionised gas. The remaining two were either observed at low resolution or had signals that were too weak to draw robust conclusions. We also present a description of the interactions between the ionised and molecular gas (as traced by CS (J=5-4)), often (but not always) finding theHII region had cleared its immediate vicinity of molecules. Of our sample of nine, the observations of the two clusters expected to have the youngest HII regions (from previous radio observations) are suggestive of having infalling motions in the H29$\alpha$ emission, which could be indicative of late stage accretion onto the stars despite the presence of an HII region.

We extend our previous study of radiative ablation of circumstellar disks by line-scattering of the star's radiation, accounting now for the effect of multiple line resonances off the stellar limb. For an analytic, three-dimensional model of the velocity structure of an equatorial Keplerian disk bounded at higher latitudes by a radially accelerating stellar wind outflow, we use root-finding methods to identify multiple resonances from a near-disk circumstellar location along starward rays both on and off the stellar core. Compared to our previous study that accounted only for the effect of on-core resonances in reducing the radiative driving through the scattering of radiation away from a near-disk circumstellar location, including off-limb resonances leads to additional radiative driving from scattering toward this location. Instead of the up-to-50% reduction in line-acceleration previously inferred from multiple resonance effects, we now find a more modest 15-20% net reduction.

The star-forming ability of a molecular cloud depends on the fraction of gas it can cycle into the dense-phase. Consequently, one of the crucial questions in reconciling star-formation in clouds is to understand the factors that control this process. While it is widely accepted that the variation in ambient conditions can alter significantly the ability of a cloud to spawn stars, the observed variation in the star-formation rate in nearby clouds that experience similar ambient conditions, presents an interesting question. In this work we attempted to reconcile this variation within the paradigm of colliding flows. To this end we develop self-gravitating, hydrodynamic realisations of identical flows, but allowed to collide off-centre. Typical observational diagnostics such as the gas-velocity dispersion, the fraction of dense-gas, the column density distribution ({\small N-PDF}), the distribution of gas mass as a function of $K$-band extinction and the strength of compressional/solenoidal modes in the post-collision cloud were deduced for different choices of the impact parameter of collision. We find that a strongly sheared cloud is terribly inefficient in cycling gas into the dense phase and that such a cloud can possibly reconcile the sluggish nature of star-formation reported for some clouds. Within the paradigm of cloud-formation via colliding flows this is possible in case of flows colliding with a relatively large impact parameter. We conclude that compressional modes - though probably essential - are insufficient to ensure a relatively higher star-formation efficiency in a cloud.

The origin of super-massive black holes in the early universe remains poorly understood.Gravitational collapse of a massive primordial gas cloud is a promising initial process,but theoretical studies have difficulty growing the black hole fast enough.We report numerical simulations of early black hole formation starting from realistic cosmological conditions.Supersonic gas motions left over from the Big Bang prevent early gas cloud formation until rapid gas condensation is triggered in a proto-galactic halo. A protostar is formed in the dense, turbulent gas cloud, and it grows by sporadic mass accretion until it acquires 34,000 solar masses.The massive star ends its life with a catastrophic collapse to leave a black hole -- a promising seed for the formation of a monstrous black hole.

Using the HPC ressources of the state of Baden-W\"urttemberg, we modelled for the first time the luminous burst from a young massive star by accretion of material from its close environment. We found that the surroundings of young massive stars are shaped as a clumpy disk whose fragments provoke outbursts once they fall onto the protostar and concluded that similar strong luminous events observed in high-mass star forming regions may be a signature of the presence of such disks.

The surroundings of massive protostars constitute an accretion disc which has numerically been shown to be subject to fragmentation and responsible for luminous accretion-driven outbursts. Moreover, it is suspected to produce close binary companions which will later strongly influence the star's future evolution in the Hertzsprung-Russel diagram. We present three-dimensional gravitation-radiation-hydrodynamic numerical simulations of 100 Mo pre-stellar cores. We find that accretion discs of young massive stars violently fragment without preventing the (highly variable) accretion of gaseous clumps onto the protostars. While acquiring the characteristics of a nascent low-mass companion, some disc fragments migrate onto the central massive protostar with dynamical properties showing that its final Keplerian orbit is close enough to constitute a close massive proto-binary system, having a young high-mass and a low-mass component. We conclude on the viability of the disc fragmentation channel for the formation of such short-period binaries, and that both processes -close massive binary formation and accretion bursts- may happen at the same time. FU-Orionis-type bursts, such as observed in the young high-mass star S255IR-NIRS3, may not only indicate ongoing disc fragmentation, but also be considered as a tracer for the formation of close massive binaries - progenitors of the subsequent massive spectroscopic binaries - once the high-mass component of the system will enter the main-sequence phase of its evolution. Finally, we investigate the ALMA-observability of the disc fragments.

The population of close-in super-Earths, with gas mass fractions of up to 10% represents a challenge for planet formation theory: how did they avoid runaway gas accretion and collapsing to hot Jupiters despite their core masses being in the critical range of $M_\mathrm{c} \simeq 10 M_\mathrm{\oplus}$? Previous three-dimensional (3D) hydrodynamical simulations indicate that atmospheres of low-mass planets cannot be considered isolated from the protoplanetary disc, contrary to what is assumed in 1D-evolutionary calculations. This finding is referred to as the recycling hypothesis. In this Paper we investigate the recycling hypothesis for super-Earth planets, accounting for realistic 3D radiation hydrodynamics. Also, we conduct a direct comparison in terms of the evolution of the entropy between 1D and 3D geometries. We clearly see that 3D atmospheres maintain higher entropy: although gas in the atmosphere loses entropy through radiative cooling, the advection of high entropy gas from the disc into the Bondi/Hill sphere slows down Kelvin-Helmholtz contraction, potentially arresting envelope growth at a sub-critical gas mass fraction. Recycling, therefore, operates vigorously, in line with results by previous studies. However, we also identify an "inner core" -- in size $\approx$ 25% of the Bondi radius -- where streamlines are more circular and entropies are much lower than in the outer atmosphere. Future studies at higher resolutions are needed to assess whether this region can become hydrodynamically-isolated on long time-scales.

The radiative transport of photons through arbitrary three-dimensional (3D) structures of dust is a challenging problem due to the anisotropic scattering of dust grains and strong coupling between different spatial regions. The radiative transfer problem in 3D is solved using Monte Carlo or Ray Tracing techniques as no full analytic solution exists for the true 3D structures. We provide the first 3D dust radiative transfer benchmark composed of a slab of dust with uniform density externally illuminated by a star. This simple 3D benchmark is explicitly formulated to provide tests of the different components of the radiative transfer problem including dust absorption, scattering, and emission. This benchmark includes models with a range of dust optical depths fully probing cases that are optically thin at all wavelengths to optically thick at most wavelengths. This benchmark includes solutions for the full dust emission including single photon (stochastic) heating as well as two simplifying approximations: One where all grains are considered in equilibrium with the radiation field and one where the emission is from a single effective grain with size-distribution-averaged properties. A total of six Monte Carlo codes and one Ray Tracing code provide solutions to this benchmark. Comparison of the results revealed that the global SEDs are consistent on average to a few percent for all but the scattered stellar flux at very high optical depths. The image results are consistent within 10%, again except for the stellar scattered flux at very high optical depths. The lack of agreement between different codes of the scattered flux at high optical depths is quantified for the first time. We provide the first 3D dust radiative transfer benchmark and validate the accuracy of this benchmark through comparisons between multiple independent codes and detailed convergence tests.

Aims: We resolve the small-scale structure around the high-mass hot core region G351.77-0.54 to investigate its disk and fragmentation properties. Methods: Using ALMA at 690GHz with baselines exceeding 1.5km, we study the dense gas, dust and outflow emission at an unprecedented spatial resolution of 0.06" (130AU@2.2kpc). Results: Within the inner few 1000AU, G351.77 fragments into at least four cores (brightness temperatures between 58 and 197K). The central structure around the main submm source #1 with a diameter of ~0.5" does not show additional fragmentation. While the CO(6-5) line wing emission shows an outflow lobe in the north-western direction emanating from source #1, the dense gas tracer CH3CN shows a velocity gradient perpendicular to the outflow that is indicative of rotational motions. Absorption profile measurements against the submm source #2 indicate infall rates on the order of 10^{-4} to 10^{-3}M_sun/yr which can be considered as an upper limit of the mean accretion rates. The position-velocity diagrams are consistent with a central rotating disk-like structure embedded in an infalling envelope, but they may also be influenced by the outflow. Using the CH_3CN(37_k-36_k) k-ladder with excitation temperatures up to 1300K, we derive a gas temperature map of source #1 exhibiting temperatures often in excess of 1000K. Brightness temperatures of the submm continuum never exceed 200K. This discrepancy between gas temperatures and submm dust brightness temperatures (in the optically thick limit) indicates that the dust may trace the disk mid-plane whereas the gas could be tracing a hotter gaseous disk surface layer. In addition, we conduct a pixel-by-pixel Toomre gravitational stability analysis of the central rotating structure. The derived high Q values throughout the structure confirm that this central region appears stable against gravitational instability.

The key aspect determining the post-formation luminosity of gas giants has long been considered to be the energetics of the accretion shock at the planetary surface. We use 1D radiation-hydrodynamical simulations to study the radiative loss efficiency and to obtain post-shock temperatures and pressures and thus entropies. The efficiency is defined as the fraction of the total incoming energy flux which escapes the system (roughly the Hill sphere), taking into account the energy recycling which occurs ahead of the shock in a radiative precursor. We focus here on a constant equation of state to isolate the shock physics but use constant and tabulated opacities. While robust quantitative results will require a self-consistent treatment including hydrogen dissocation and ionization, the results show the correct qualitative behavior and can be understood semi-analytically. The shock is found to be isothermal and supercritical for a range of conditions relevant to core accretion (CA), with Mach numbers greater than ca. 3. Across the shock, the entropy decreases significantly, by a few entropy units (k_B/baryon). While nearly 100 percent of the incoming kinetic energy is converted to radiation locally, the efficiencies are found to be as low as roughly 40 percent, implying that a meaningful fraction of the total accretion energy is brought into the planet. For realistic parameter combinations in the CA scenario, a non-zero fraction of the luminosity always escapes the system. This luminosity could explain, at least in part, recent observations in the LkCa 15 and HD 100546 systems.

Recently reported variations in the typical physical properties of Galactic and extra-Galactic molecular clouds (MCs), and in their ability to form stars have been attributed to local variations in the magnitude of interstellar pressure. Inferences from these surveys have called into question two long-standing beliefs that the MCs : 1 are Virialised entities and (2) have approximately constant surface density i.e., the validity of the Larson's third law. In this work we invoke the framework of cloud-formation via collisions between warm gas flows. Post-collision clouds forming in these realisations cool rapidly and evolve primarily via the interplay between the Non-linear Thin Shell Instability (NTSI), and the self-gravity. Over the course of these simulations we traced the temporal evolution of the surface density of the assembled clouds, the fraction of dense gas, the distribution of gas column density (NPDF), and the Virial nature of the assembled clouds. We conclude, these physical properties of MCs not only exhibit temporal variation, but their respective peak-magnitude also increases in proportion with the magnitude of external pressure, $P_{ext}$. The velocity dispersion in assembled clouds appears to follow the power-law, $\sigma_{gas}\propto P_{ext}^{0.23}$. Also, the power-law tail at higher densities becomes shallower with increasing magnitude of external pressure, for magnitudes, $P_{ext}/k_{B}\lesssim 10^{7}$ K cm$^{-3}$, at higher magnitudes such as those typically found in the Galactic CMZ ($P_{ext}/k_{B} > 10^{7}$ K cm$^{-3}$), the power-law shows significant steepening. Thus while our results are broadly consistent with inferences from various recent observational surveys, it appears, MCs hardly exhibit a unique set of properties, but rather a wide variety, that can be reconciled with a range of magnitudes of pressure between 10$^{4}$ K cm$^{-3}$ - 10$^{8}$ K cm$^{-3}$.

The prompt emission of gamma-ray bursts (GRBs) is characterized by rapid variabilities, which may be a direct reflection of the unsteady central engine. We perform a series of axisymmetric 2.5-dimensional simulations to study the propagation of relativistic, hydrodynamic, intermittent jets through the envelope of a GRB progenitor star. A realistic rapidly rotating star is incorporated as the background of jet propagation, and the star is allowed to collapse due to the gravity of the central black hole. By modeling the intermittent jets with constant-luminosity pulses with equal on and off durations, we investigate how the half-period, $T$, affects the jet dynamics. For relatively small $T$ values (e.g. 0.2 s), the jet breakout time $t_{\rm bo}$ depends on the opening angle of the jet, with narrower jets more penetrating and reaching the surface at shorter times. For $T \leq 1$ s, the reverse shock crosses each pulse before the jet penetrates through the stellar envelope. As a result, after the breakout of the first group of pulses at $t_{\rm bo}$, several subsequent pulses vanish before penetrating the star, causing a quiescent gap. For larger half-periods ($T=2.0, 4.0$ s), all the pulses can successfully penetrate through the envelope, since each pulse can propagate through the star before the reverse shock crosses the shell. Our results may interpret the existence of a weak precursor in many long GRBs, given that the GRB central engine injects intermittent pulses with half-period $T \leq 1$ s. The observational data seem to be consistent with such a possibility.

A significant fraction of OB-type, main-sequence massive stars are classified as runaway and move supersonically through the interstellar medium (ISM). Their strong stellar winds interact with their surroundings where the typical strength of the local ISM magnetic field is about 3.5-7 micro-G, which can result in the formation of bow shock nebulae. We investigate the effects of such magnetic fields, aligned with the motion of the flow, on the formation and emission properties of these circumstellar structures. Our axisymmetric, magneto-hydrodynamical simulations with optically-thin radiative cooling, heating and anisotropic thermal conduction show that the presence of the background ISM magnetic field affects the projected optical emission our bow shocks at Ha and [OIII] lambda 5007 which become fainter by about 1-2 orders of magnitude, respectively. Radiative transfer calculations against dust opacity indicate that the magnetic field slightly diminishes their projected infrared emission and that our bow shocks emit brightly at 60 micron. This may explain why the bow shocks generated by ionizing runaway massive stars are often difficult to identify. Finally, we discuss our results in the context of the bow shock of Zeta Ophiuchi and we support the interpretation of its imperfect morphology as an evidence of the presence of an ISM magnetic field not aligned with the motion of its driving star.

We perform two-dimensional axially symmetric radiation-hydrodynamic simulations to assess the impact of outflows and radiative force feedback from massive protostars by varying when the protostellar outflow starts, the ratio of ejection to accretion rates, and the strength of the wide angle disk wind component. The star formation efficiency, i.e. the ratio of final stellar mass to initial core mass, is dominated by radiative forces and the ratio of outflow to accretion rates. Increasing this ratio has three effects: First, the protostar grows slower with a lower luminosity at any given time, lowering radiative feedback. Second, bipolar cavities cleared by the outflow are larger, further diminishing radiative feedback on disk and core scales. Third, the higher momentum outflow sweeps up more material from the collapsing envelope, decreasing the protostar's potential mass reservoir via entrainment. The star formation efficiency varies with the ratio of ejection to accretion rates from 50% in the case of very weak outflows to as low as 20% for very strong outflows. At latitudes between the low density bipolar cavity and the high density accretion disk, wide angle disk winds remove some of the gas, which otherwise would be part of the accretion flow onto the disk; varying the strength of these wide angle disk winds, however, alters the final star formation efficiency by only +/-6%. For all cases, the opening angle of the bipolar outflow cavity remains below 20 degree during early protostellar accretion phases, increasing rapidly up to 65 degree at the onset of radiation pressure feedback.

Accretion-driven luminosity outbursts are a vivid manifestation of variable mass accretion onto protostars. They are known as the so-called FU Orionis phenomenon in the context of low-mass protostars. More recently, this process has been found in models of primordial star formation. Using numerical radiation hydrodynamics simulations, we stress that present-day forming massive stars also experience variable accretion and show that this process is accompanied by luminous outbursts induced by the episodic accretion of gaseous clumps falling from the circumstellar disk onto the protostar. Consequently, the process of accretion-induced luminous flares is also conceivable in the high-mass regime of star formation and we propose to regard this phenomenon as a general mechanism that can affect protostars regardless of their mass and/or the chemical properties of the parent environment in which they form. In addition to the commonness of accretion-driven outbursts in the star formation machinery, we conjecture that luminous flares from regions hosting forming high-mass star may be an observational implication of the fragmentation of their accretion disks.

We present coupled stellar evolution (SE) and 3D radiation-hydrodynamic (RHD) simulations of the evolution of primordial protostars, their immediate environment, and the dynamic accretion history under the influence of stellar ionizing and dissociating UV feedback. Our coupled SE-RHD calculations result in a wide diversity of final stellar masses covering 10 Msun $\lesssim M_* \lesssim$ 1000 Msun. The formation of very massive ($\gtrsim$ 250 Msun) stars is possible under weak UV feedback, whereas ordinary massive (a few x 10 Msun) stars form when UV feedback can efficiently halt the accretion. This may explain the peculiar abundance pattern of a Galactic metal-poor star recently reported by Aoki et al. (2014), possibly the observational signature of very massive precursor primordial stars. Weak UV feedback occurs in cases of variable accretion, in particular when repeated short accretion bursts temporarily exceed 0.01 Msun/yr, causing the protostar to inflate. In the bloated state, the protostar has low surface temperature and UV feedback is suppressed until the star eventually contracts, on a thermal adjustment timescale, to create an HII region. If the delay time between successive accretion bursts is sufficiently short, the protostar remains bloated for extended periods, initiating at most only short periods of UV feedback. Disk fragmentation does not necessarily reduce the final stellar mass. Quite the contrary, we find that disk fragmentation enhances episodic accretion as many fragments migrate inward and are accreted onto the star, thus allowing continued stellar mass growth under conditions of intermittent UV feedback. This trend becomes more prominent as we improve the resolution of our simulations. We argue that simulations with significantly higher resolution than reported previously are needed to derive accurate gas mass accretion rates onto primordial protostars.

We present radiation hydrodynamic simulations of collapsing protostellar cores with initial masses of 30, 100, and 200 M$_{\odot}$. We follow their gravitational collapse and the formation of a massive protostar and protostellar accretion disk. We employ a new hybrid radiative feedback method blending raytracing techniques with flux-limited diffusion for a more accurate treatment of the temperature and radiative force. In each case, the disk that forms becomes Toomre-unstable and develops spiral arms. This occurs between 0.35 and 0.55 freefall times and is accompanied by an increase in the accretion rate by a factor of 2-10. Although the disk becomes unstable, no other stars are formed. In the case of our 100 and 200 M$_{\odot}$ simulation, the star becomes highly super-Eddington and begins to drive bipolar outflow cavities that expand outwards. These radiatively-driven bubbles appear stable, and appear to be channeling gas back onto the protostellar accretion disk. Accretion proceeds strongly through the disk. After 81.4 kyr of evolution, our 30 M$_{\odot}$ simulation shows a star with a mass of 5.48 M$_{\odot}$ and a disk of mass 3.3 M$_{\odot}$, while our 100 M$_{\odot}$ simulation forms a 28.8 M$_{\odot}$ mass star with a 15.8 M$_{\odot}$ disk over the course of 41.6 kyr, and our 200 M$_{\odot}$ simulation forms a 43.7 M$_{\odot}$ star with an 18 M$_{\odot}$ disk in 21.9 kyr. In the absence of magnetic fields or other forms of feedback, the masses of the stars in our simulation do not appear limited by their own luminosities.

Biasing or importance sampling is a powerful technique in Monte Carlo radiative transfer, and can be applied in different forms to increase the accuracy and efficiency of simulations. One of the drawbacks of the use of biasing is the potential introduction of large weight factors. We discuss a general strategy, composite biasing, to suppress the appearance of large weight factors. We use this composite biasing approach for two different problems faced by current state-of-the-art Monte Carlo radiative transfer codes: the generation of photon packages from multiple components, and the penetration of radiation through high optical depth barriers. In both cases, the implementation of the relevant algorithms is trivial and does not interfere with any other optimisation techniques. Through simple test models, we demonstrate the general applicability, accuracy and efficiency of the composite biasing approach. In particular, for the penetration of high optical depths, the gain in efficiency is spectacular for the specific problems that we consider: in simulations with composite path length stretching, high accuracy results are obtained even for simulations with modest numbers of photon packages, while simulations without biasing cannot reach convergence, even with a huge number of photon packages.

Massive stars that have been ejected from their parent cluster and supersonically sailing away through the interstellar medium (ISM) are classified as exiled. They generate circumstellar bow shock nebulae that can be observed. We present two-dimensional, axisymmetric hydrodynamical simulations of a representative sample of stellar wind bow shocks from Galactic OB stars in an ambient medium of densities ranging from n_ISM=0.01 up to 10.0/cm3. Independently of their location in the Galaxy, we confirm that the infrared is the most appropriated waveband to search for bow shocks from massive stars. Their spectral energy distribution is the convenient tool to analyze them since their emission does not depend on the temporary effects which could affect unstable, thin-shelled bow shocks. Our numerical models of Galactic bow shocks generated by high-mass (~40 Mo) runaway stars yield H$\alpha$ fluxes which could be observed by facilities such as the SuperCOSMOS H-Alpha Survey. The brightest bow shock nebulae are produced in the denser regions of the ISM. We predict that bow shocks in the field observed at Ha by means of Rayleigh-sensitive facilities are formed around stars of initial mass larger than about 20 Mo. Our models of bow shocks from OB stars have the emission maximum in the wavelength range 3 <= lambda <= 50 micrometer which can be up to several orders of magnitude brighter than the runaway stars themselves, particularly for stars of initial mass larger than 20 Mo.

The supersonic motion of gravitating objects through a gaseous medium constitutes a classical problem in theoretical astrophysics. Its application covers a broad range of objects and scales from planets up to galaxies. Especially the dynamical friction, caused by the forming wake behind the object, plays an important role for the dynamics of the system. To calculate the dynamical friction, standard formulae, based on linear theory are often used. It is our goal to check the general validity of these formulae and provide suitable expressions for the dynamical friction acting on the moving object, based on the basic physical parameters of the problem. We perform sequences of high resolution numerical studies of rigid bodies moving supersonically through a homogeneous medium, and calculate the total drag acting on the object, which is the sum of gravitational and hydro drag. We study cases without gravity with purely hydrodynamical drag, as well as gravitating objects. From the final equilibrium state of the simulations, we compute for gravitating objects the dynamical friction by direct numerical integration of the gravitational pull acting on the embedded object. The numerical experiments confirm the known scaling laws for the dependence of the dynamical friction on the basic physical parameters as derived in earlier semi-analytical studies. As a new important result we find that the shock's stand-off distance is revealed as the minimum spatial interaction scale of dynamical friction. Below this radius, the gas settles into a hydrostatic state, which causes no net gravitational pull onto the moving body. Finally, we derive an analytic estimate for the stand-off distance that can be used to calculate the dynamical friction force.

We present Atacama Large Millimeter/submillimeter Array (ALMA) line and continuum observations at 1.2mm with ~0.3" resolution that uncover a Keplerian-like disk around the forming O-type star AFGL 4176. The continuum emission from the disk at 1.21 mm (source mm1) has a deconvolved size of 870+/-110 AU x 330+/-300 AU and arises from a structure ~8 M_sun in mass, calculated assuming a dust temperature of 190 K. The first-moment maps, pixel-to-pixel line modeling, assuming local thermodynamic equilibrium (LTE), and position-velocity diagrams of the CH3CN J=13-12 K-line emission all show a velocity gradient along the major axis of the source, coupled with an increase in velocity at small radii, consistent with Keplerian-like rotation. The LTE line modeling shows that where CH3CN J=13-12 is excited, the temperatures in the disk range from ~70 to at least 300 K and that the H2 column density peaks at 2.8x10^24 cm^-2. In addition, we present Atacama Pathfinder Experiment (APEX) 12CO observations which show a large-scale outflow from AFGL 4176 perpendicular to the major axis of mm1, supporting the disk interpretation. Finally, we present a radiative transfer model of a Keplerian disk surrounding an O7 star, with a disk mass and radius of 12 M_sun and 2000 AU, that reproduces the line and continuum data, further supporting our conclusion that our observations have uncovered a Keplerian disk around an O-type star.

Aims: Understanding the fragmentation and collapse properties of the dense gas during the onset of high-mass star formation. Methods: We observed the massive (~800M_sun) starless gas clump IRDC18310-4 with the Plateau de Bure Interferometer (PdBI) at sub-arcsecond resolution in the 1.07mm continuum andN2H+(3-2) line emission. Results: Zooming from a single-dish low-resolution map to previous 3mm PdBI data, and now the new 1.07mm continuum observations, the sub-structures hierarchically fragment on the increasingly smaller spatial scales. While the fragment separations may still be roughly consistent with pure thermal Jeans fragmentation, the derived core masses are almost two orders of magnitude larger than the typical Jeans mass at the given densities and temperatures. However, the data can be reconciled with models using non-homogeneous initial density structures, turbulence and/or magnetic fields. While most sub-cores remain (far-)infrared dark even at 70mum, we identify weak 70mum emission toward one core with a comparably low luminosity of ~16L_sun, re-enforcing the general youth of the region. The spectral line data always exhibit multiple spectral components toward each core with comparably small line widths for the individual components (in the 0.3 to 1.0km/s regime). Based on single-dish C18O(2-1) data we estimate a low virial-to-gas-mass ratio <=0.25. We discuss that the likely origin of these spectral properties may be the global collapse of the original gas clump that results in multiple spectral components along each line of sight. Even within this dynamic picture the individual collapsing gas cores appear to have very low levels of internal turbulence.

StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase of HII region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and SPH- based codes. The first test examines the `early phase' D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the `late phase' D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all twelve participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to $\lesssim2\%$. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a dataset for calibrating the treatment of ionizing radiation hydrodynamics codes.

We carry out radiation hydrodynamical simulations of the formation of massive stars in the super-Eddington regime including both their radiative feedback and protostellar outflows. The calculations start from a prestellar core of dusty gas and continue until the star stops growing. The accretion ends when the remnants of the core are ejected, mostly by the force of the direct stellar radiation in the polar direction and elsewhere by the reradiated thermal infrared radiation. How long the accretion persists depends on whether the protostellar outflows are present. We set the mass outflow rate to 1% of the stellar sink particle's accretion rate. The outflows open a bipolar cavity extending to the core's outer edge, through which the thermal radiation readily escapes. The radiative flux is funneled into the polar directions while the core's collapse proceeds near the equator. The outflow thus extends the "flashlight effect", or anisotropic radiation field, found in previous studies from the few hundred AU scale of the circumstellar disk up to the 0.1 parsec scale of the core. The core's flashlight effect allows core gas to accrete on the disk for longer, in the same way that the disk's flashlight effect allows disk gas to accrete on the star for longer. Thus although the protostellar outflows remove material near the core's poles, causing slower stellar growth over the first few free-fall times, they also enable accretion to go on longer in our calculations. The outflows ultimately lead to stars of somewhat higher mass.

In the core accretion paradigm of planet formation, gas giants only form a massive atmosphere after their progenitors exceeded a threshold mass: the critical core mass. Most (exo)planets, being smaller and rock/ice-dominated, never crossed this line. Nevertheless, they were massive enough to attract substantial amounts of gas from the disc, while their atmospheres remained in pressure-equilibrium with the disc. Our goal is to characterise the hydrodynamical properties of the atmospheres of such embedded planets and their implication for their (long-term) evolution. In this paper -- the first in series -- we start to investigate the properties of an isothermal and inviscid flow past a small, embedded planet by conducting local, 2D hydrodynamical simulations. Using the PLUTO code we confirm that the flow is steady and bound. This steady outcome is most apparent for the log-polar grid (with the grid spacing proportional to the distance from the planet). For low-mass planets, Cartesian grids are somewhat less efficient as they have difficulty to follow the circular, large speeds in the deep atmosphere. Relating the amount of rotation to the gas fraction of the atmosphere, we find that more massive atmospheres rotate faster -- a finding consistent with Kelvin's circulation theorem. Rotation therefore limits the amount of gas that planets can acquire from the nebula. Dependent on the Toomre-Q parameter of the circumstellar disc, the planet's atmosphere will reach Keplerian rotation before self-gravity starts to become important.

Following Paper I we investigate the properties of atmospheres that form around small protoplanets embedded in a protoplanetary disc by conducting hydrodynamical simulations. These are now extended to three dimensions, employing a spherical grid centred on the planet. Compression of gas is shown to reduce rotational motions. Contrasting the 2D case, no clear boundary demarcates bound atmospheric gas from disc material; instead, we find an open system where gas enters the Bondi sphere at high latitudes and leaves through the midplane regions, or, vice versa, when the disc gas rotates sub-Keplerian. The simulations do not converge to a time-independent solution; instead, the atmosphere is characterized by a time-varying velocity field. Of particular interest is the timescale to replenish the atmosphere by nebular gas, $t_\mathrm{replenish}$. It is shown that the replenishment rate, $M_\mathrm{atm}/t_\mathrm{replenish}$, can be understood in terms of a modified Bondi accretion rate, $\sim$$R_\mathrm{Bondi}^2\rho_\mathrm{gas}v_\mathrm{Bondi}$, where $v_\mathrm{Bondi}$ is set by the Keplerian shear or the magnitude of the sub-Keplerian motion of the gas, whichever is larger. In the inner disk, the atmosphere of embedded protoplanets replenishes on a timescale that is shorter than the Kelvin-Helmholtz contraction (or cooling) timescale. As a result, atmospheric gas can no longer contract and the growth of these atmospheres terminates. Future work must confirm whether these findings continue to apply when the (thermodynamical) idealizations employed in this study are relaxed. But if shown to be broadly applicable, replenishment of atmospheric gas provides a natural explanation for the preponderance of gas-rich but rock-dominant planets like super-Earths and mini-Neptunes.

Radiation feedback plays a crucial role in the process of star formation. In order to simulate the thermodynamic evolution of disks, filaments, and the molecular gas surrounding clusters of young stars, we require an efficient and accurate method for solving the radiation transfer problem. We describe the implementation of a hybrid radiation transport scheme in the adaptive grid-based FLASH general magnetohydrodynamics code. The hybrid scheme splits the radiative transport problem into a raytracing step and a diffusion step. The raytracer captures the first absorption event, as stars irradiate their environments, while the evolution of the diffuse component of the radiation field is handled by a flux-limited diffusion (FLD) solver. We demonstrate the accuracy of our method through a variety of benchmark tests including the irradiation of a static disk, subcritical and supercritical radiative shocks, and thermal energy equilibration. We also demonstrate the capability of our method for casting shadows and calculating gas and dust temperatures in the presence of multiple stellar sources. Our method enables radiation-hydrodynamic studies of young stellar objects, protostellar disks, and clustered star formation in magnetized, filamentary environments.

In a molecular cloud dust opacity typically dominates over gas opacity, yet in the vicinities of forming stars dust is depleted, and gas is the sole provider of opacity. In the optically thin circumstellar environments the radiation temperature cannot be assumed to be equal to the gas temperature, hence the two-temperature Planck means are necessary to calculate the radiative equilibrium. By using the two-temperature mean opacity one does obtain the proper equilibrium gas temperature in a circumstellar environment, which is in a chemical equilibrium. A careful consideration of a radiative transfer problem reveals that the equilibrium temperature solution can be degenerate in an optically thin gaseous environment. We compute mean gas opacities based on the publicly available code DFSYNTHE by Kurucz and Castelli. We performed the calculations assuming local thermodynamic equilibrium and an ideal gas equation of state. The values were derived by direct integration of the high-resolution opacity spectrum. We produced two sets of gas opacity tables: Rosseland means and two-temperature Planck means (the tables available via http://cdsweb.u-strasbg.fr/ as well as via http://www.mpia-hd.mpg.de/homes/malygin). For three metallicities [Me/H] = 0.0,+/-0.3 we covered the parameter range 3.48 <= log T_rad[K] <= 4.48 in radiation temperature, 2.8 <= log T_gas[K]} <= 6.0 in gas temperature, and -10 <= log P[dyn cm^-2] <= 6 in gas pressure. We show that in the optically thin circumstellar environment for a given stellar radiation field and local gas density there are several equilibrium gas temperatures possible. We conclude that, in general, equilibrium gas temperature cannot be determined without treating the temperature evolution.

Studies of the evolution of massive protostars and the evolution of their host molecular cloud cores are commonly treated as separate problems. However, interdependencies between the two can be significant. Here, we study the simultaneous evolution of massive protostars and their host molecular cores using a multi-dimensional radiation hydrodynamics code that incorporates the effects of the thermal pressure and radiative acceleration feedback of the centrally forming protostar. The evolution of the massive protostar is computed simultaneously using the stellar evolution code STELLAR, modified to include the effects of variable accretion. The interdependencies are studied in three different collapse scenarios. For comparison, stellar evolutionary tracks at constant accretion rates and the evolution of the host cores using pre-computed stellar evolutionary tracks are computed. The resulting interdependencies of the protostellar evolution and the evolution of the environment are extremely diverse and depend on the order of events, in particular the time of circumstellar accretion disk formation with respect to the onset of the bloating phase of the star. Feedback mechanisms affect the instantaneous accretion rate and the protostar's radius, temperature and luminosity on timescales equal or smaller than 5 kyr, corresponding to the accretion timescale and Kelvin-Helmholtz contraction timescale, respectively. Nevertheless, it is possible to approximate the overall protostellar evolution in many cases by pre-computed stellar evolutionary tracks assuming appropriate constant average accretion rates.

Context: Dynamical studies of irradiated circumstellar disks require an accurate treatment of radiation transport to, for example, properly determine cooling and fragmentation properties. The radiation transport algorithm should be as fast as the (magneto-) hydrodynamics to allow for an efficient usage of computing resources. Methods: We use a setup of a central star and a slightly flared circumstellar disk. We perform simulations for a wide range of optical depths of the disk's midplane from tau(550nm) = 0.1 up to tau(810nm) = 1 million. We check the accuracy of the gray flux-limited diffusion (FLD) approximation and a gray and frequency-dependent ray-tracing plus FLD approximation. Results: 1. For moderate optical depths, a gray approximation of the stellar irradiation yields a slightly hotter inner rim and a slightly cooler midplane of the disk at larger radii, but is otherwise in agreement with the frequency-dependent treatment. 2. The gray FLD approximation fails to compute an appropriate temperature profile in all regimes of optical depth; the maximum deviations to the comparison runs are 50 percent in the optically thin and up to 280 percent in the optically thick limit. For low optical depth, the isotropic assumption within the FLD method yields a too steep decrease of the radial temperature slope. For higher optical depths, the FLD approximation does not reproduce the shadow behind the optically thick inner rim of the circumstellar disk, yielding artificial heating at larger disk radii. 3. The frequency-dependent RT + gray FLD approximation yields remarkable accuracy for the whole range of optical depths. Conclusions: The high accuracy of the frequency-dependent hybrid radiation transport algorithm makes this method ideally suited for (magneto-) hydrodynamical studies of irradiated circumstellar disks.

We review our recent studies demonstrating that the radiation pressure problem in the formation of massive stars can be circumvented via an anisotropy of the thermal radiation field. Such an anisotropy naturally establishes with the formation of a circumstellar disk. The required angular momentum transport within the disk can be provided by developing gravitational torques. Radiative Rayleigh-Taylor instabilities in the cavity regions - as previously suggested in the literature - are not required and are shown to be not occurring in the context of massive star formation.

Numerical simulations have shown that the often cited radiation pressure barrier to accretion onto massive stars can be circumvented, when the radiation field is highly anisotropic in the presence of a circumstellar accretion disk with high optical depth. Here, these studies of the so-called flashlight effect are expanded by including the opacity of the innermost dust-free but potentially optically thick gas regions around forming massive stars. In addition to frequency-dependent opacities for the dust grains, we use temperature- and density-dependent Planck- and Rosseland mean opacities for the gas. The simulations show that the innermost dust-free parts of the accretion disks are optically thick to the stellar radiation over a substantial fraction of the solid angle above and below the disk's midplane. The temperature in the shielded disk region decreases faster with radius than in a comparison simulation with a lower constant gas opacity, and the dust sublimation front is shifted to smaller radii. The shielding by the dust-free gas in the inner disk thus contributes to an enhanced flashlight effect, which ultimately results in a smaller opening angle of the radiation pressure driven outflow and in a much longer timescale of sustained feeding of the circumstellar disk by the molecular cloud core. We conclude that it is necessary to properly account for the opacity of the inner dust-free disk regions around forming massive stars in order to correctly assess the effectiveness of the flashlight effect, the opening angle of radiation pressure driven outflows, and the lifetime and morphological evolution of the accretion disk.

We present a new publicly available tool (DustPol) aimed to model the polarised thermal dust emission. The module DustPol, which is publicly available, is part of the ARTIST (Adaptable Radiative Transfer Innovations for Submillimetre Telescopes) package, which also offers tools for modelling the polarisation of line emission together with a model library and a Python-based user interface. DustPol can easily manage analytical as well as pre-gridded models to generate synthetic maps of the Stokes I, Q, and U parameters. These maps are stored in FITS format which is straightforwardly read by the data reduction software used, e.g., by the Atacama Large Millimeter Array (ALMA). This turns DustPol into a powerful engine for the prediction of the expected polarisation features of a source observed with ALMA or the Planck satellite as well as for the interpretation of existing submillimetre observations obtained with other telescopes. DustPol allows the parameterisation of the maximum degree of polarisation and we find that, in a prestellar core, if there is depolarisation, this effect should happen at densities of 10^6 cm-3 or larger. We compare a model generated by DustPol with the observational polarisation data of the low-mass Class 0 object NGC 1333 IRAS 4A, finding that the total and the polarised emission are consistent.

Context: When massive stars exert a radiation pressure onto their environment that is higher than their gravitational attraction, they launch a radiation-pressure-driven outflow. It has been claimed that a radiative Rayleigh-Taylor instability should lead to the collapse of the outflow cavity and foster the growth of massive stars. Aims: We investigate the stability of idealized radiation-pressure-dominated cavities, focusing on its dependence on the radiation transport approach for the stellar radiation feedback. Methods: We compare two different methods for stellar radiation feedback: gray flux-limited diffusion (FLD) and ray-tracing (RT). We also derive simple analytical models to support our findings. Results: Only the FLD cases lead to prominent instability in the cavity shell. The RT cases do not show such instability. The gray FLD method underestimates the opacity at the location of the cavity shell and leads to extended epochs of marginal Eddington equilibrium in the cavity shell, making them prone to the radiative Rayleigh-Taylor instability. In the RT cases, the radiation pressure exceeds gravity by 1-2 orders of magnitude. The radiative Rayleigh-Taylor instability is then consequently suppressed. Conclusions: Treating the stellar irradiation in the gray FLD approximation underestimates the radiative forces acting on the cavity shell. This can lead artificially to situations that are affected by the radiative Rayleigh-Taylor instability. The proper treatment of direct stellar irradiation by massive stars is crucial for the stability of radiation-pressure-dominated cavities.

The most massive stars can form via standard disk accretion - despite of the radiation pressure generated - due to the fact that the massive accretion disk yields a strong anisotropy in the radiation field, releasing most of the radiation pressure perpendicular to the disk accretion flow. Here, we analyze the self-gravity of the forming circumstellar disk as the potential major driver of the angular momentum transport in such massive disks responsible for the high accretion rates needed for the formation of massive stars. For this purpose, we perform self-gravity radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. The formation and evolution of the resulting circumstellar disk is investigated in 1.) axially symmetric simulations using an alpha-shear-viscosity prescription and 2.) a three-dimensional simulation, in which the angular momentum transport is provided self-consistently by developing gravitational torques in the self-gravitating accretion disk. The simulation series of different strength of the alpha-viscosity shows that the accretion history of the forming star is mostly independent of the alpha-viscosity-parameter. The accretion history of the three-dimensional run driven by self-gravity is more time-dependent than the viscous disk evolution in axial symmetry. The mean accretion rate, i.e. the stellar mass growth, is nearly identical to the alpha-viscosity models. We conclude that the development of gravitational torques in self-gravitating disks around forming massive stars provides a self-consistent mechanism to efficiently transport the angular momentum to outer disk radii. Also the formation of the most massive stars can therefore be understood in the standard accretion disk scenario.

We investigate the radiation pressure feedback in the formation of massive stars in 1, 2, and 3D radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. In contrast to previous research, we consider frequency dependent stellar radiation feedback, resolve the dust sublimation front in the vicinity of the forming star down to 1.27 AU, compute the evolution for several 10^5 yrs covering the whole accretion phase of the forming star, and perform a comprehensive survey of the parameter space. The most fundamental result is that the formation of a massive accretion disk in slowly rotating cores preserves a high anisotropy in the radiation field. The thermal radiation escapes through the optically thin atmosphere, effectively diminishing the radiation pressure feedback onto the accretion flow. Gravitational torques in the self-gravitating disk drive a sufficiently high accretion rate to overcome the residual radiation pressure. Simultaneously, the radiation pressure launches an outflow in the bipolar direction, which grows in angle with time and releases a substantial fraction of the initial core mass from the star-disk system. Summarized, for an initial core mass of 60, 120, 240, and 480 Msol these mechanisms allow the star to grow up to 28.2, 56.5, 92.6, and at least 137.2 Msol respectively.

We present radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. We treat frequency dependent radiative feedback from stellar evolution and accretion luminosity at a numerical resolution down to 1.27 AU. In the 2D approximation of axially symmetric simulations, it is possible for the first time to simulate the whole accretion phase of several 10^5 yr for the forming massive star and to perform a comprehensive scan of the parameter space. Our simulation series show evidently the necessity to incorporate the dust sublimation front to preserve the high shielding property of massive accretion disks. Our disk accretion models show a persistent high anisotropy of the corresponding thermal radiation field, yielding to the growth of the highest-mass stars ever formed in multi-dimensional radiation hydrodynamics simulations. Non-axially symmetric effects are not necessary to sustain accretion. The radiation pressure launches a stable bipolar outflow, which grows in angle with time as presumed from observations. For an initial mass of the pre-stellar host core of 60, 120, 240, and 480 Msol the masses of the final stars formed in our simulations add up to 28.2, 56.5, 92.6, and at least 137.2 Msol respectively.

We present radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. We treat frequency dependent radiative feedback from stellar evolution and accretion luminosity at a numerical resolution down to 1.27 AU. In the 2D approximation of axially symmetric simulations, it is possible for the first time to simulate the whole accretion phase (up to the end of the accretion disk epoch) for the forming massive star and to perform a broad scan of the parameter space. Our simulation series show evidently the necessity to incorporate the dust sublimation front to preserve the high shielding property of massive accretion disks. While confirming the upper mass limit of spherically symmetric accretion, our disk accretion models show a persistent high anisotropy of the corresponding thermal radiation field. This yields to the growth of the highest-mass stars ever formed in multi-dimensional radiation hydrodynamics simulations, far beyond the upper mass limit of spherical accretion. Non-axially symmetric effects are not necessary to sustain accretion. The radiation pressure launches a stable bipolar outflow, which grows in angle with time as presumed from observations. For an initial mass of the pre-stellar host core of 60, 120, 240, and 480 Msun the masses of the final stars formed in our simulations add up to 28.2, 56.5, 92.6, and at least 137.2 Msun respectively.

Context: Radiative feedback plays a crucial role in the formation of massive stars. The implementation of a fast and accurate description of the proceeding thermodynamics in pre-stellar cores and evolving accretion disks is therefore a main effort in current hydrodynamics simulations. Aims: We introduce our newly implemented three-dimensional frequency dependent radiation transport algorithm for hydrodynamics simulations of spatial configurations with a dominant central source. Methods: The module combines the advantage of the speed of an approximate Flux Limited Diffusion (FLD) solver with the high accuracy of a frequency dependent first order ray-tracing routine. Results: We prove the viability of the scheme in a standard radiation benchmark test compared to a full frequency dependent Monte-Carlo based radiative transfer code. The setup includes a central star, a circumstellar flared disk, as well as an envelope. The test is performed for different optical depths. Considering the frequency dependence of the stellar irradiation, the temperature distributions can be described precisely in the optically thin, thick, and irradiated transition regions. Resulting radiative forces onto dust grains are reproduced with high accuracy. The achievable parallel speedup of the method imposes no restriction on further radiative (magneto-) hydrodynamics simulations. Conclusions: The proposed approximate radiation transport method enables frequency dependent radiation hydrodynamics studies of the evolution of pre-stellar cores and circumstellar accretion disks around an evolving massive star in a highly efficient and accurate manner.

We briefly overview our newly developed radiation transport module for MHD simulations and two actual applications. The method combines the advantage of the speed of the Flux-Limited Diffusion approximation and the high accuracy obtained in ray-tracing methods.