BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Rolf Kuiper

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.GAastro-ph.EPastro-ph.HEastro-ph.IMastro-ph.COphysics.comp-ph

Catalog footprint

What is connected

38works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

38 published item(s)

preprint2026arXiv

The Structure and Kinematics of Three Class 0 Protostellar Jets from JWST

We present observations of jets within 2000 au of three deeply embedded protostars using 2.9-27 micron observations with JWST. These observations show the morphologies and kinematics of the collimated jets from three protostars, the low-mass Class 0 protostars B335 and HOPS 153, and the intermediate-mass protostar HOPS 370. These jets are traced by shock-ionized fine-structure line emission observed with the JWST NIRSpec and MIRI IFUs. We find that [Fe II] emission traces the full extent of the inner 1000 to 2000 au of the jets, depending on distance to the protostar, while other ions mostly trace isolated shocked knots. The jets show evidence of wiggling motion in the plane of the sky as well as asymmetries between blue and red-shifted lobes. The widths of the jets increase non-monotonically with distance from the central protostar, with opening angles ranging from 2.1 degrees to < 10.1 degrees for the three protostars in the sample. The jets have total velocities ranging from 147 to 184 km/s after correcting for disk inclination. For B335, an 8-month gap between NIRSpec and MIRI MRS observations enabled measurement of the tangential velocity of a shocked knot; in combination with the radial velocity, this shows that the jet has a different inclination than the outflow cavity. We find multiple knots before and during a recent outburst in B335, although the knots were more frequent during the burst. The asymmetries between blue- and red-shifted lobes strongly suggest complex interactions between the circumstellar disks and magnetic fields.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Unveiling the First O-Type Bloated Star Candidate through ALMA and EVLA Observations

We investigate the circumstellar environment of the O-type bloated star candidate IRAS 19520+2759 (I19520) using high-resolution observations from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Expanded Very Large Array (EVLA). Radio continuum emission traced by the EVLA (C, K, and Q bands) exhibits a spectral index of 0.5, consistent with a thermal jet. ALMA 1.3 mm continuum map reveals a compact source coincident with the optical counterpart of I19520, likely tracing the dense core hosting the central massive young stellar object. A prominent molecular outflow in the east-west direction, along with a possible secondary outflow oriented northeast-southwest, is identified in the $^{13}\mathrm{CO}$ emission. A hot molecular core and a Keplerian disk are detected in several $\mathrm{SO}_2$ transitions. Assuming an edge-on disk geometry, the dynamical mass of the central object is estimated to be in the range of $10$-$15~M_\odot$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Rapid growth of seed black holes during early bulge formation

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Snapshot of a magnetohydrodynamic disk wind traced by water maser observations

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Formation and Evolution of Disks around Young Stellar Objects

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

FRB 200428: an Impact between an Asteroid and a Magnetar

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gaseous dynamical friction under radiative feedback: do intermediate-mass black holes speed up or down?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modeling disk fragmentation and multiplicity in massive star formation

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the turbulence driving mode of expanding HII regions

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 $σ_s^2 = \ln({1+b^2\mathcal{M}^2})$ between the variance of the logarithmic density contrast $σ_s^2$ (where $s = \ln({ρ/ρ_0})$ with the gas density $ρ$ and its average $ρ_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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Composite biasing in Monte Carlo radiative transfer

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical friction for supersonic motion in a homogeneous gaseous medium

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Formation of Massive Primordial Stars: Intermittent UV Feedback with Episodic Mass Accretion

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Protostellar Outflows and Radiative Feedback from Massive Stars. II. Feedback, Star Formation Efficiency, and Outflow Broadening

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Simulating the Formation of Massive Protostars: I. Radiative Feedback and Accretion Disks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Keplerian-like disk around the forming O-type star AFGL 4176

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A general hybrid radiation transport scheme for star formation simulations on an adaptive grid

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hydrodynamics of Embedded Planets' First Atmospheres. I. A Centrifugal Growth Barrier for 2D Flows

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hydrodynamics of Embedded Planets' First Atmospheres. II. A Rapid Recycling of Atmospheric Gas

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ρ_\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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Protostellar Outflows and Radiative Feedback from Massive Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the simultaneous evolution of massive protostars and their host cores

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The reliability of approximate radiation transport methods for irradiated disk studies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A solution to the radiation pressure problem in the formation of massive stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the effects of optically thick gas (disks) around massive stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the stability of radiation-pressure-dominated cavities

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radiation pressure feedback in the formation of massive stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Three-dimensional simulation of massive star formation in the disk accretion scenario

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Circumventing the radiation pressure barrier in the formation of massive stars via disk accretion

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fast and accurate frequency-dependent radiation transport for hydrodynamics simulations in massive star formation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The role of accretion disks in the formation of massive stars

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.

0 citations0 reviewsNo code linkedOpen access