Improving the thin-disk models of circumstellar disk evolution. The 2+1-dimensional model

Circumstellar disks of gas and dust are naturally formed from contracting pre-stellar molecular cores during the star formation process. To study various dynamical and chemical processes that take place in circumstellar disks prior to their dissipation and transition to debris disks, the appropriate numerical models capable of studying the long-term disk chemodynamical evolution are required. We present a new 2+1-dimensional numerical hydrodynamics model of circumstellar disk evolution, in which the thin-disk model is complemented with the procedure for calculating the vertical distributions of gas volume density and temperature in the disk. The reconstruction of the disk vertical structure is performed at every time step via the solution of the time-dependent radiative transfer equations coupled to the equation of the vertical hydrostatic equilibrium. We perform a detailed comparison between circumstellar disks produced with our previous 2D model and with the improved 2+1D approach. The structure and evolution of resulting disks, including the differences in temperatures, densities, disk masses and protostellar accretion rates, are discussed in detail. The new 2+1D model yields systematically colder disks, while the in-falling parental clouds are warmer. Both effects act to increase the strength of disk gravitational instability and, as a result, the number of gravitationally bound fragments that form in the disk via gravitational fragmentation as compared to the purely 2D thin-disk simulations with a simplified thermal balance calculation.Comment: Accepted for publication in Astronomy & Astrophysic

Stellar hydrodynamical modeling of dwarf galaxies: simulation methodology, tests, and first results

Cosmological simulations still lack numerical resolution or physical processes to simulate dwarf galaxies in sufficient details. Accurate numerical simulations of individual dwarf galaxies are thus still in demand. We aim at (i) studying in detail the coupling between stars and gas in a galaxy, exploiting the so-called stellar hydrodynamical approach, and (ii) studying the chemo-dynamical evolution of individual galaxies starting from self-consistently calculated initial gas distributions. We present a novel chemo-dynamical code in which the dynamics of gas is computed using the usual hydrodynamics equations, while the dynamics of stars is described by the stellar hydrodynamics approach, which solves for the first three moments of the collisionless Boltzmann equation. The feedback from stellar winds and dying stars is followed in detail. In particular, a novel and detailed approach has been developed to trace the aging of various stellar populations, which enables an accurate calculation of the stellar feedback depending on the stellar age. We build initial equilibrium models of dwarf galaxies that take gas self-gravity into account and present different levels of rotational support. Models with high rotational support develop prominent bipolar outflows; a newly-born stellar population in these models is preferentially concentrated to the galactic midplane. Models with little rotational support blow away a large fraction of the gas and the resulting stellar distribution is extended and diffuse. The stellar dynamics turns out to be a crucial aspect of galaxy evolution. If we artificially suppress stellar dynamics, supernova explosions occur in a medium heated and diluted by the previous activity of stellar winds, thus artificially enhancing the stellar feedback (abridged).Comment: 22 pages, 19 figures, accepted for publication in Astronomy & Astrophysic

Lifetime of the embedded phase of low-mass star formation and the envelope depletion rates

Motivated by a considerable scatter in the observationally inferred lifetimes of the embedded phase of star formation, we study the duration of the Class 0 and Class I phases in upper-mass brown dwarfs and low-mass stars using numerical hydrodynamics simulations of the gravitational collapse of a large sample of cloud cores. We resolve the formation of a star/disk/envelope system and extend our numerical simulations to the late accretion phase when the envelope is nearly totally depleted of matter. We adopted a classification scheme of Andre et al. and calculate the lifetimes of the Class 0 and Class I phases (\tau_C0 and \tau_CI, respectively) based on the mass remaining in the envelope. When cloud cores with various rotation rates, masses, and sizes (but identical otherwise) are considered, our modeling reveals a sub-linear correlation between the Class 0 lifetimes and stellar masses in the Class 0 phase with the least-squares fit exponent m=0.8 \pm 0.05. The corresponding correlation between the Class I lifetimes and stellar masses in the Class I is super-linear with m=1.2 \pm 0.05. If a wider sample of cloud cores is considered, which includes possible variations in the initial gas temperature, cloud core truncation radii, density enhancement amplitudes, initial gas density and angular velocity profiles, and magnetic fields, then the corresponding exponents may decrease by as much as 0.3. The duration of the Class I phase is found to be longer than that of the Class~0 phase in most models, with a mean ratio \tau_CI / \tau_C0 \approx 1.5--2. A notable exception are YSOs that form from cloud cores with large initial density enhancements, in which case \tau_C0 may be greater than \tau_CI. Moreover, the upper-mass (>= 1.0 Msun) cloud cores with frozen-in magnetic fields and high cloud core rotation rates may have the \tau_CI / \tau_C0 ratios as large as 3.0--4.0. (Abdridged).Comment: Accepted for publication by The Astrophysical Journa

The Ejection of Low Mass Clumps During Star Formation

Modeling of the self-consistent formation and evolution of disks as a result of prestellar core collapse reveals an intense early phase of recurrent gravitational instability and clump formation. These clumps generally migrate inward due to gravitational interaction with trailing spiral arms, and can be absorbed into the central object. However, in situations of multiple clump formation, gravitational scattering of clumps can result in the ejection of a low mass clump. These clumps can then give rise to free-floating low mass stars, brown dwarfs, or even giant planets. Detailed modeling of this process in the context of present-day star formation reveals that these clumps start out essentially as Larson first cores and grow subsequently by accretion. In the context of Pop III star formation, preliminary indications are that the disk clumps may also be of low mass. This mechanism of clump formation and possible ejection provides a channel for the formation of low mass objects in the first generation of stars.Comment: 4 pages, 2 figures, to appear in proceedings of First Stars IV meeting (Kyoto, Japan; 2012

The effect of episodic accretion on the phase transition of CO and CO_2 in low-mass star formation

We study the evaporation and condensation of CO and CO_2 during the embedded stages of low-mass star formation by using numerical simulations. We focus on the effect of luminosity bursts, similar in magnitude to FUors and EXors, on the gas-phase abundance of CO and CO_2 in the protostellar disk and infalling envelope. The evolution of a young protostar and its environment is followed based on hydrodynamical models using the thin-disk approximation, coupled with a stellar evolution code and phase transformations of CO and CO_2. The accretion and associated luminosity bursts in our model are caused by disk gravitational fragmentation followed by quick migration of the fragments onto the forming protostar. We found that bursts with luminosity on the order of 100-200 L_sun can evaporate CO ices in part of the envelope. The typical freeze-out time of the gas-phase CO onto dust grains in the envelope (a few kyr) is much longer than the burst duration (100-200 yr). This results in an increased abundance of the gas-phase CO in the envelope long after the system has returned into a quiescent stage. In contrast, luminosity bursts can evaporate CO_2 ices only in the disk, where the freeze-out time of the gas-phase CO_2 is comparable to the burst duration. We thus confirm that luminosity bursts can leave long-lasting traces in the abundance of gas-phase CO in the infalling envelope, enabling the detection of recent bursts as suggested by previous semi-analytical studies.Comment: 12 pages, 6 figures, accepted for publication in Astronomy & Astrophysic

A Hybrid Scenario for the Formation of Brown Dwarfs and Very Low Mass Stars

We present a calculation of protostellar disk formation and evolution in which gaseous clumps (essentially, the first Larson cores formed via disk fragmentation) are ejected from the disk during the early stage of evolution. This is a universal process related to the phenomenon of ejection in multiple systems of point masses. However, it occurs in our model entirely due to the interaction of compact, gravitationally-bound gaseous clumps and is free from the smoothing-length uncertainty that is characteristic of models using sink particles. Clumps that survive ejection span a mass range of 0.08--0.35 $M_\odot$, and have ejection velocities $0.8 \pm 0.35$ km s$^{-1}$, which are several times greater than the escape speed. We suggest that, upon contraction, these clumps can form substellar or low-mass stellar objects with notable disks, or even close-separation very-low-mass binaries. In this hybrid scenario, allowing for ejection of clumps rather than finished protostars/proto--brown-dwarfs, disk formation and the low velocity dispersion of low-mass objects are naturally explained, while it is also consistent with the observation of isolated low-mass clumps that are ejection products. We conclude that clump ejection and the formation of isolated low mass stellar and substellar objects is a common occurrence, with important implications for understanding the initial mass function, the brown dwarf desert, and the formation of stars in all environments and epochs.Comment: 20 pages, 6 figures, to appear in The Astrophysical Journa

Formation of freely floating sub-stellar objects via close encounters

We numerically studied close encounters between a young stellar system hosting a massive, gravitationally fragmenting disk and an intruder diskless star with the purpose to determine the evolution of fragments that have formed in the disk prior to the encounter. Numerical hydrodynamics simulations in the non-inertial frame of reference of the host star were employed to simulate the prograde and retrograde co-planar encounters. The initial configuration of the target system (star plus disk) was obtained via a separate numerical simulation featuring the gravitational collapse of a solar-mass pre-stellar core. We found that close encounters can lead to the ejection of fragments that have formed in the disk of the target prior to collision. In particular, prograde encounters are more efficient in ejecting the fragments than the retrograde encounters. The masses of ejected fragments are in the brown-dwarf mass regime. They also carry away an appreciable amount of gas in their gravitational radius of influence, implying that these objects may possess extended disks or envelopes, as also suggested by Thies et al. (2015). Close encounters can also lead to the ejection of entire spiral arms, followed by fragmentation and formation of freely-floating objects straddling the planetary mass limit. However, numerical simulations with a higher resolution are needed to confirm this finding.Comment: 12 pages, 7 figures, accepted for publication by Astronomy & Astrophysic