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

Variable protostellar accretion with episodic bursts

We present the latest development of the disk gravitational instability and fragmentation model, originally introduced by us to explain episodic accretion bursts in the early stages of star formation. Using our numerical hydrodynamics model with improved disk thermal balance and star-disk interaction, we computed the evolution of protostellar disks formed from the gravitational collapse of prestellar cores. In agreement with our previous studies, we find that cores of higher initial mass and angular momentum produce disks that are more favorable to gravitational instability and fragmentation, while a higher background irradiation and magnetic fields moderate the disk tendency to fragment. The protostellar accretion in our models is time-variable, thanks to the nonlinear interaction between different spiral modes in the gravitationally unstable disk, and can undergo episodic bursts when fragments migrate onto the star owing to the gravitational interaction with other fragments or spiral arms. Most bursts occur in the partly embedded Class I phase, with a smaller fraction taking place in the deeply embedded Class 0 phase and a few possible bursts in the optically visible Class II phase. The average burst duration and mean luminosity are found to be in good agreement with those inferred from observations of FU-Orionis-type eruptions. The model predicts the existence of two types of bursts: the isolated ones, showing well-defined luminosity peaks separated with prolonged periods (~ 10^4 yr) of quiescent accretion, and clustered ones, demonstrating several bursts occurring one after another during just a few hundred years. Finally, we estimate that 40%--70% of the star-forming cores can display bursts after forming a star-disk system.

The effect of external environment on the evolution of protostellar disks

Using numerical hydrodynamics simulations we studied the gravitational collapse of pre-stellar cores of sub-solar mass embedded into a low-density external environment. Four models with different magnitude and direction of rotation of the external environment with respect to the central core were studied and compared with an isolated model. We found that the infall of matter from the external environment can significantly alter the disk properties as compared to those seen in the isolated model. Depending on the magnitude and direction of rotation of the external environment, a variety of disks can form including compact (<= 200 AU) ones shrinking in size due to infall of external matter with low angular momentum, as well as extended disks forming due to infall of external matter with high angular momentum. The former are usually stable against gravitational fragmentation, while the latter are prone to fragmentation and formation of stellar systems with sub-stellar/very-low-mass companions. In the case of counterrotating external environment, very compact (< 5 AU) and short-lived (<= a few * 10^5 yr) disks can form when infalling material has low angular momentum. The most interesting case is found for the infall of counterrotating external material with high angular momentum, leading to the formation of counterrotating inner and outer disks separated by a deep gap at a few tens AU. The gap migrates inward due to accretion of the inner disk onto the protostar, turns into a central hole, and finally disappears giving way to the outer strongly gravitationally unstable disk. This model may lead to the emergence of a transient stellar system with sub-stellar/very-low-mass components counterrotating with respect to that of the star.

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.

Fragmenting protostellar disks: properties and observational signatures

Using numerical hydrodynamic simulations, we study the gravitational fragmentation of an unstable protostellar disc formed during the collapse of a pre-stellar core with a mass of 1.2 M_sun. The forming fragments span a mass range from about a Jupiter mass to very-low-mass protostars and are located at distances from a few tens to a thousand AU, with a dearth of objects at < 100 AU. We explore the possibility of observational detection of the fragments in discs viewed through the outflow cavity at a distance of 250 pc. We demonstrate that one hour of integration time with the Atacama Large Millimeter/sub-millimeter Array (ALMA) is sufficient to detect the fragments with masses as low as 1.5 M_Jup at orbital distances up to 800 AU from the protostar. The ALMA resolution sets the limit on the minimum orbital distance of detectable fragments. For the adopted resolution of our simulated ALMA images of 0.1, the fragments can be detected at distances down to 50 AU. At smaller distances, the fragments usually merge with the central density peak. The likelihood for detecting the fragments reduces significantly for a lower resolution of 0.5. Some of the most massive fragments, regardless of their orbital distance, can produce characteristic peaks at approximately 5 micron and hence their presence can be indirectly inferred from the observed spectral energy distributions of protostars.

Formation of giant planets and brown dwarfs on wide orbits

Context: We studied numerically the formation of giant planet (GP) and brown dwarf (BD) embryos in gravitationally unstable protostellar disks and compared our findings with directly-imaged, wide-orbit (>= 50 AU) companions known to-date. The viability of the disk fragmentation scenario for the formation of wide-orbit companions in protostellar disks around (sub-)solar mass stars was investigated. Methods: We used numerical hydrodynamics simulations of disk formation and evolution with an accurate treatment of disk thermodynamics. The use of the thin-disk limit allowed us to probe the long-term evolution of protostellar disks. We focused on models that produced wide-orbit GP/BD embryos, which opened a gap in the disk and showed radial migration timescales similar to or longer than the typical disk lifetime. Results: While disk fragmentation was seen in the majority of our models, only 6 models out of 60 revealed the formation of quasi-stable, wide-orbit GP/BD embryos. Disk fragmentation produced GP/BD embryos with masses in the 3.5-43 M_J range, covering the whole mass spectrum of directly-imaged, wide-orbit companions to (sub-)solar mass stars. On the other hand, our modelling failed to produce embryos on orbital distances <= 170 AU, whereas several directly-imaged companions were found at smaller orbits down to a few AU. Disk fragmentation also failed to produce wide-orbit companions around stars with mass <= 0.7 Msun, in disagreement with observations. Conclusions: Disk fragmentation is unlikely to explain the whole observed spectrum of wide-orbit companions to (sub-)solar-mass stars and other formation mechanisms, e.g., dynamical scattering of closely-packed companions onto wide orbits, should be invoked to account for companions at orbital distance from a few tens to approx 150 AU and wide-orbit companions with masses of the host star <= 0.7 Msun. (abridged)