Self-gravitating equilibrium models of dwarf galaxies and the minimum mass for star formation

We construct a series of model galaxies in rotational equilibrium consisting of gas, stars, and a fixed dark matter (DM) halo and study how these equilibrium systems depend on the mass and form of the DM halo, gas temperature, non-thermal and rotation support against gravity, and also on the redshift of galaxy formation. For every model galaxy we find the minimum gas mass M_g^min required to achieve a state in which star formation (SF) is allowed according to contemporary SF criteria. The obtained M_g^min--M_DM relations are compared against the baryon-to-DM mass relation M_b--M_DM inferred from the LambdaCDM theory and WMAP4 data. Our aim is to construct realistic initial models of dwarf galaxies (DGs), which take into account the gas self-gravity and can be used as a basis to study the dynamical and chemical evolution of DGs. Rotating equilibria are found by solving numerically the steady-state momentum equation for the gas component in the combined gravitational potential of gas, stars, and DM halo using a forward substitution procedure. We find that for a given M_DM the value of M_g^min depends crucially on the gas temperature T_g, gas spin parameter alpha, degree of non-thermal support sigma_eff, and somewhat on the redshift for galaxy formation z_gf. Depending on the actual values of T_g, alpha, sigma_eff, and z_gf, model galaxies may have M_g^min that are either greater or smaller than M_b. Galaxies with M_DM ga 10^9 M_sun are usually characterized by M_g^min la M_b, implying that SF in such objects is a natural outcome as the required gas mass is consistent with what is available according to the LambdaCDM theory. On the other hand, models with M_DM la 10^9 M_sun are often characterized by M_g^min >> M_b, implying that they need much more gas than available to achieve a state in which SF is allowed. Abridged.

Destruction of massive fragments in protostellar disks and crystalline silicate production

We present a mechanism for the crystalline silicate production associated with the formation and subsequent destruction of massive fragments in young protostellar disks. The fragments form in the embedded phase of star formation via disk fragmentation at radial distances ga 50-100 AU and anneal small amorphous grains in their interior when the gas temperature exceeds the crystallization threshold of ~ 800 K. We demonstrate that fragments that form in the early embedded phase can be destroyed before they either form solid cores or vaporize dust grains, thus releasing the processed crystalline dust into various radial distances from sub-AU to hundred-AU scales. Two possible mechanisms for the destruction of fragments are the tidal disruption and photoevaporation as fragments migrate radially inward and approach the central star and also dispersal by tidal torques exerted by spiral arms. As a result, most of the crystalline dust concentrates to the disk inner regions and spiral arms, which are the likely sites of fragment destruction.

Embedded protostellar disks around (sub-)solar stars. II. Disk masses, sizes, densities, temperatures and the planet formation perspective

We present basic properties of protostellar disks in the embedded phase of star formation (EPSF), which is difficult to probe observationally using available observational facilities. We use numerical hydrodynamics simulations of cloud core collapse and focus on disks formed around stars in the 0.03-1.0 Msun mass range. Our obtained disk masses scale near-linearly with the stellar mass. The mean and median disk masses in the Class 0 and I phases (M_{d,C0}^{mean}=0.12 Msun, M_{d,C0}^{mdn}=0.09 Msun and M_{d,CI}^{mean}=0.18 Msun, M_{d,CI}^{mdn}=0.15 Msun, respectively) are greater than those inferred from observations by (at least) a factor of 2--3. We demonstrate that this disagreement may (in part) be caused by the optically thick inner regions of protostellar disks, which do not contribute to millimeter dust flux. We find that disk masses and surface densities start to systematically exceed that of the minimum mass solar nebular for objects with stellar mass as low as M_st=0.05-0.1 Msun. Concurrently, disk radii start to grow beyond 100 AU, making gravitational fragmentation in the disk outer regions possible. Large disk masses, surface densities, and sizes suggest that giant planets may start forming as early as in the EPSF, either by means of core accretion (inner disk regions) or direct gravitational instability (outer disk regions), thus breaking a longstanding stereotype that the planet formation process begins in the Class II phase.

Embedded protostellar disks around (sub-)solar protostars. I. Disk structure and evolution

We perform a comparative numerical hydrodynamics study of embedded protostellar disks formed as a result of the gravitational collapse of cloud cores of distinct mass (M_cl=0.2--1.7 M_sun) and ratio of rotational to gravitational energy (beta=0.0028--0.023). An increase in M_cl and/or beta leads to the formation of protostellar disks that are more susceptible to gravitational instability. Disk fragmentation occurs in most models but its effect is often limited to the very early stage, with the fragments being either dispersed or driven onto the forming star during tens of orbital periods. Only cloud cores with high enough M_cl or beta may eventually form wide-separation binary/multiple systems with low mass ratios and brown dwarf or sub-solar mass companions. It is feasible that such systems may eventually break up, giving birth to rogue brown dwarfs. Protostellar disks of {it equal} age formed from cloud cores of greater mass (but equal beta) are generally denser, hotter, larger, and more massive. On the other hand, protostellar disks formed from cloud cores of higher beta (but equal M_cl) are generally thinner and colder but larger and more massive. In all models, the difference between the irradiation temperature and midplane temperature triangle T is small, except for the innermost regions of young disks, dense fragments, and disks outer edge where triangle T is negative and may reach a factor of two or even more. Gravitationally unstable, embedded disks show radial pulsations, the amplitude of which increases along the line of increasing M_cl and beta but tends to diminish as the envelope clears. We find that single stars with a disk-to-star mass ratio of order unity can be formed only from high-beta cloud cores, but such massive disks are unstable and quickly fragment into binary/multiple systems.

The age of blue LSB galaxies

(Abridged). Low metallicities, large gas-to-star mass ratios, and blue colors of most low surface brightness (LSB) galaxies imply that these systems may be younger than their high surface brightness counterparts. We seek to find observational signatures that can help to constrain the age of blue LSB galaxies. We use numerical hydrodynamic modelling to study the long-term (~13 Gyr) dynamical and chemical evolution of blue LSB galaxies adopting a sporadic scenario for star formation. Our models utilize various rates of star formation and different shapes of the initial mass function (IMF). We complement hydrodynamic modelling with population synthesis modelling to produce the integrated B-V colors and Halpha equivalent widths (EW(Ha)). We find that the mean oxygen abundances, B-V colors, EW(Ha), and the radial fluctuations in the oxygen abundance, when considered altogether, can be used to constrain the age of blue LSB galaxies if some independent knowledge of the IMF is available. Our modelling strongly suggests the existence of a minimum age for blue LSB galaxies. Model B-V colors and mean oxygen abundances set a tentative minimum age at 1.5-3.0 Gyr, whereas model EW(Ha) suggest a larger value of order 5-6 Gyr. The latter value may decrease somewhat, if blue LSB galaxies host IMFs with a truncated upper mass limit. We found no firm evidence that the age of blue LSB galaxies is significantly smaller than 13 Gyr.

Peculiarities in the stellar velocity distribution of galaxies with a two-armed spiral structure

We expand our pervious numerical study of the properties of the stellar velocity distribution within the disk of a two-armed spiral galaxy by considering spiral stellar density waves with different global Fourier amplitudes, C_2. We confirm our previous conclusion that the ratio sigma_1:sigma_2 of smallest versus largest principal axes of the stellar velocity ellipsoid becomes abnormally small near the outer edges of the stellar spiral arms. The extent to which the stellar velocity ellipsoid is elongated (as compared to the unperturbed value typical for the axisymmetric disk) increases with the strength of the spiral density wave. In particular, the C_2=0.06 spiral can decrease the unperturbed value of sigma_1:sigma_2 by 20%, while the C_2=0.13 spiral can decrease the unperturbed sigma_1:sigma_2 by a factor of 3. The abnormally small values of the sigma_1:sigma_2 ratio can potentially be used to track the position of {it stellar} spiral density waves. The sigma_{phiphi}:sigma_{rr} ratio is characterized by a more complex behaviour and exhibits less definite minima near the outer edges of the spiral arms. We find that the epicycle approximation is violated near the spiral arms and cannot be used in spiral galaxies with C_2 >= 0.05-0.06 or in galaxies with the amplitude of the spiral stellar density wave (relative to the unperturbed background) of order 0.1 or greater.

Mass accretion rates in self-regulated disks of T Tauri stars

We have studied numerically the evolution of protostellar disks around intermediate and upper mass T Tauri stars (0.25 M_sun < M_st < 3.0 M_sun) that have formed self-consistently from the collapse of molecular cloud cores. In the T Tauri phase, disks settle into a self-regulated state, with low-amplitude nonaxisymmetric density perturbations persisting for at least several million years. Our main finding is that the global effect of gravitational torques due to these perturbations is to produce disk accretion rates that are of the correct magnitude to explain observed accretion onto T Tauri stars. Our models yield a correlation between accretion rate M_dot and stellar mass M_st that has a best fit M_dot propto M_st^{1.7}, in good agreement with recent observations. We also predict a near-linear correlation between the disk accretion rate and the disk mass.