Researcher profile

Rolf Kuiper

Rolf Kuiper contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

17 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.

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$.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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&#39;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.

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.

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.

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.

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$.

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.