The early evolution of protostellar disks with metallicities in the $Z=1.0-0.01~Z_\odot$ range was studied with a particular emphasis on the strength of gravitational instability and the nature of protostellar accretion in low-metallicity systems. Numerical hydrodynamics simulations in the thin-disk limit were employed that feature separate gas and dust temperatures, and disk mass-loading from the infalling parental cloud cores. Models with cloud cores of similar initial mass and rotation pattern, but distinct metallicity were considered to distinguish the effect of metallicity from that of initial conditions. The early stages of disk evolution in low-metallicity models are characterized by vigorous gravitational instability and fragmentation. Disk instability is sustained by continual mass-loading from the collapsing core. The time period that is covered by this unstable stage is much shorter in the $Z=0.01~Z_\odot$ models as compared to their higher metallicity counterparts thanks to the higher mass infall rates caused by higher gas temperatures (that decouple from lower dust temperatures) in the inner parts of collapsing cores. Protostellar accretion rates are highly variable in the low-metallicity models reflecting a highly dynamical nature of the corresponding protostellar disks. The low-metallicity systems feature short, but energetic episodes of mass accretion caused by infall of inward-migrating gaseous clumps that form via gravitational fragmentation of protostellar disks. These bursts seem to be more numerous and last longer in the $Z=0.1~Z_\odot$ models in comparison to the $Z=0.01~Z_\odot$ case. Variable protostellar accretion with episodic bursts is not a particular feature of solar metallicity disks. It is also inherent to gravitationally unstable disks with metallicities up to 100 times lower than solar.

While accreting through a circumstellar disk, young stellar objects are observed to undergo sudden and powerful accretion events known as FUor or EXor outbursts. Although such episodic accretion is considered to be an integral part of the star formation process, the triggers and mechanisms behind them remain uncertain. We conducted global numerical hydrodynamics simulations of protoplanetary disk formation and evolution in the thin-disk limit, assuming both magnetically layered and fully magnetorotational instability (MRI)-active disk structure. In this paper, we characterize the nature of the outbursts occurring in these simulations. The instability in the dead zone of a typical layered disk results in "MRI outbursts". We explore their progression and their dependence on the layered disk parameters as well as cloud core mass. The simulations of fully MRI-active disks showed an instability analogous to the classical thermal instability. This instability manifested at two temperatures--above approximately 1400 K and 3500 K--due to the steep dependence of Rosseland opacity on the temperature. The origin of these thermally unstable regions is related to the bump in opacity resulting from molecular absorption by water vapor and may be viewed as a novel mechanism behind some of the shorter duration accretion events. Although we demonstrated local thermal instability in the disk, more investigations are needed to confirm that a large-scale global instability will ensue. We conclude that the magnetic structure of a disk, its composition, as well as the stellar mass, can significantly affect the nature of episodic accretion in young stellar objects.

Aims: We explore the long-term evolution of young protoplanetary disks with different approaches to computing the thermal structure determined by various cooling and heating processes in the disk and its surroundings. Methods: Numerical hydrodynamics simulations in the thin-disk limit were complemented with three thermal evolution schemes: a simplified $\beta$-cooling approach with and without irradiation, in which the rate of disk cooling is proportional to the local dynamical time, a fiducial model with equal dust and gas temperatures calculated taking viscous heating, irradiation, and radiative cooling into account, and also a more sophisticated approach allowing decoupled dust and gas temperatures. Results: We found that the gas temperature may significantly exceed that of dust in the outer regions of young disks thanks to additional compressional heating caused by the infalling envelope material in the early stages of disk evolution and slow collisional exchange of energy between gas and dust in low-density disk regions. The outer envelope however shows an inverse trend with the gas temperatures dropping below that of dust. The global disk evolution is only weakly sensitive to temperature decoupling. Nevertheless, separate dust and gas temperatures may affect the chemical composition, dust evolution, and disk mass estimates. Constant-$\beta$ models without stellar and background irradiation fail to reproduce the disk evolution with more sophisticated thermal schemes because of intrinsically variable nature of the $\beta$-parameter. Constant-$\beta$ models with irradiation can better match the dynamical and thermal evolution, but the agreement is still incomplete. Conclusions: Models allowing separate dust and gas temperatures are needed when emphasis is placed on the chemical or dust evolution in protoplanetary disks, particularly in sub-solar metallicity environments.

We study the dynamics and growth of dust particles in circumstellar disks of different masses that are prone to gravitational instability during the critical first Myr of their evolution. The dust component is made up of two different components: micron-sized dust and grown dust of evolving size. For the dust component, we considered the dust coagulation, fragmentation, momentum exchange with the gas, and dust self-gravity. We found that the micron-sized dust particles grow rapidly in the circumstellar disk, reaching a few cm in size in the inner 100 au of the disk during less than 100 kyr after the disk formation, provided that fragmentation velocity is $30\rm~ms^{-1}$. Due to the accretion of micron dust particles from the surrounding envelope, which serves as a micron dust reservoir, the approximately cm-sized dust particles continue to be present in the disk for more than 900 kyr after the disk formation and maintain a dust-to-gas ratio close to 0.01. We show that a strong correlation exists between the gas and pebble fluxes in the disk. We find that radial surface density distribution of pebbles in the disk shows power-law distribution with an index similar to that of the Minimum-mass solar nebula (MMSN) regardless the disk mass. We also show that the gas surface density in our models agrees well with measurements of dust in protoplanetary disks of AS 209, HD 163296, and DoAr 25 systems. Pebbles are formed during the very early stages of protoplanetary disk evolution. They play a crucial role in the planet formation process. Our disc simulations reveal the early onset ($<10^5$ yr) of an inwards-drifting flux of pebble-sized particles that makes up approximately between one hundredth and one tenth of the gas mass flux, which appears consistent with mm-observations of discs. (Abridged)

We study the origin of tail-like structures recently detected around the disk of SU Aurigae and several FU~Orionis-type stars. Dynamic protostellar disks featuring ejections of gaseous clumps and quiescent protoplanetary disks experiencing a close encounter with an intruder star were modeled using the numerical hydrodynamics code FEOSAD. Both the gas and dust dynamics were taken into account, including dust growth and mutual friction between the gas and dust components. Only plane-of-the-disk encounters were considered. Ejected clumps produce a unique type of tail that is characterized by a bow-shock shape. Such tails originate from the supersonic motion of ejected clumps through the dense envelope that often surrounds young gravitationally unstable protostellar disks. The ejected clumps either sit at the head of the tail-like structure or disperse if their mass is insufficient to withstand the head wind of the envelope. On the other hand, close encounters with quiescent protoplanetary disks produce three types of the tail-like structure; we define these as pre-collisional, post-collisional, and spiral tails. These tails can in principle be distinguished from one another by particular features of the gas and dust flow in and around them. We find that the brown-dwarf-mass intruders do not capture circumintruder disks during the encounter, while the subsolar-mass intruders can acquire appreciable circumintruder disks with elevated dust-to-gas ratios, which can ease their observational detection. However, this is true only for prograde collisions; the retrograde intruders fail to collect appreciable amounts of gas or dust from the disk of the target. The predicted mass of dust in the model tail-like structures is higher than what was inferred for similar structures in SU~Aur, FU~Ori, and Z~CMa, making their observational detection feasible. Abridged.

Luminosity outbursts of the FUOri type stars, which have a magnitude of $\sim100 L_{\odot}$ and last for decades, may affect chemical composition of the surrounding protoplanetary disk. Using astrochemical modeling we analyze the changes induced by the outburst and search for species sensitive to the luminosity rise. Some changes in the disk molecular composition appear not only during the outburst itself but can also retain for decades after the end of the outburst. We analyze main chemical processes responsible for these effects and assess timescales at which chemically inert species return to the pre-outburst abundances.

It is believed that large-scale horseshoe-like brightness asymmetries found in dozens of transitional protoplanetary discs are caused by anticyclonic vortices. These vortices can play a key role in planet formation, as mm-sized dust -- the building blocks of planets -- can be accumulated inside them. Anticyclonic vortices are formed by the Rossby wave instability, which can be excited at the gap edges opened by a giant planet or at sharp viscosity transitions of accretionally inactive regions. It is known that vortices are prone to stretching and subsequent dissolution due to disc self-gravity for canonical disc masses in the isothermal approximation. To improve the hydrodynamic model of protoplanetary discs, we include the disc thermodynamics in our model. In this paper, we present our results on the evolution of the vortices formed at the outer edge of an accretionally inactive region (dead zone) assuming an ideal equation of state and taking $PdV$ work, disc cooling in the $\beta$-approximation, and disc self-gravity into account. Thermodynamics affects the offset and the mode number (referring to the number of small vortices at the early phase) of the RWI excitation, as well as the strength, shape, and lifetime of the large-scale vortex formed through merging of the initial small vortices. We found that the inclusion of gas thermodynamics results in stronger, however decreased lifetime vortices. Our results suggest that a hypothetical vortex-aided planet formation scenario favours effectively cooling discs.

The role of convection in the gas-dust accretion disk around a young star is studied. The evolution of a Keplerian disk is modeled using the Pringle equation, which describes the time variations of the surface density under the action of turbulent viscosity. The distributions of the density and temperature in the polar directions are computed simultaneously in the approximation that the disk is hydrostatically stable. The computations of the vertical structure of the disk take into account heating by stellar radiation, interstellar radiation, and viscous heating. The main factor governing evolution of the disk in this model is the dependence of the viscosity coefficient on the radius of the disk. The computations of this coefficient take into account the background viscosity providing the continuous accretion of the gas and the convective viscosity, which depends on the parameters of the convection at a given radius. The results of computations of the global evolution and morphology of the disk obtained in this approach are presented. It is shown that, in the adopted model, the accretion has burst-like character: after the inner part of the disk ( R < 3 AU) is filled with matter, this material is relatively fast discharged onto the star, after which the process is repeated. Our results may be useful for explaining the activity of young FU Ori and EX Lup objects. It is concluded that convection may be one of the mechanisms responsible for the non-steady pattern of accretion in protostellar disks.

We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $\mu$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/.

Massive young stellar object (MYSOs) form during the collapse of high-mass pre-stellar cores, where infalling molecular material is accreted through a centrifugally-balanced accretion disc that is subject to efficient gravitational instabilities. In the resulting fragmented accretion disc of the MYSO, gaseous clumps and low-mass stellar companions can form, which will influence the future evolution of massive protostars in the Hertzsprung-Russell diagram. We perform dust continuum radiative transfer calculations and compute synthetic images of disc structures modelled by the gravito-radiation-hydrodynamics simulation of a forming MYSO, in order to investigate the Atacama Large Millimeter/submillimeter Array (ALMA) observability of circumstellar gaseous clumps and forming multiple systems. Both spiral arms and gaseous clumps located at ~a few 100 au from the protostar can be resolved by interferometric ALMA Cycle 7 C43-8 and C43-10 observations at band 6 (1.2 mm), using a maximal 0.015" beam angular resolution and at least 10-30 min exposure time for sources at distances of 1-2 kpc. Our study shows that substructures are observable regardless of their viewing geometry or can be inferred in the case of an edge-viewed disc. The observation probability of the clumps increases with the gradually increasing efficiency of gravitational instability at work as the disc evolves. As a consequence, large discs around MYSOs close to the zero-age-main-sequence line exhibit more substructures than at the end of the gravitational collapse. Our results motivate further observational campaigns devoted to the close surroundings of the massive protostars S255IR-NIRS3 and NGC 6334I-MM1, whose recent outbursts are a probable signature of disc fragmentation and accretion variability.

We present ALMA observations of 12CO, 13CO, and C18O J=2--1 lines and the 230 GHz continuum for the FU Ori-type object (FUor) V900 Mon (d~1.5 kpc), for which the accretion burst was triggered between 1953 and 2009. We identified CO emission associated with a molecular bipolar outflow extending up to a ~10^4 au scale and a rotating molecular envelope extending over >10^4 au. The interaction with the hot energetic FUor wind, which was observed using optical spectroscopy, appears limited to a region within ~400 au of the star. The envelope mass and the collimation of the extended CO outflow suggest that the progenitor of this FUor is a low-mass Class I young stellar object (YSO). These parameters for V900 Mon, another FUor, and a few FUor-like stars are consistent with the idea that FUor outbursts are associated with normal YSOs. The continuum emission is marginally resolved in our observations with a 0."2x0."15 (~300x225 au) beam, and a Gaussian model provides a deconvolved FWHM of ~90 au. The emission is presumably associated with a dusty circumstellar disk, plus a possible contribution from a wind or a wind cavity close to the star. The warm compact nature of the disk continuum emission could be explained with viscous heating of the disk, while gravitational fragmentation in the outer disk and/or a combination of grain growth and their inward drift may also contribute to its compact nature.

Spatial distribution and growth of dust in a clumpy protoplanetary disk subject to vigorous gravitational instability and fragmentation is studied numerically with sub-au resolution using the FEOSAD code. Hydrodynamics equations describing the evolution of self-gravitating and viscous protoplanetary disks in the thin-disk limit were modified to include a dust component consisting of two parts: sub-micron-sized dust and grown dust with a variable maximum radius. The conversion of small to grown dust, dust growth, friction of dust with gas, and dust self-gravity were also considered. We found that the disk appearance is notably time-variable with spiral arms, dusty rings, and clumps, constantly forming, evolving, and decaying. As a consequence, the total dust-to-gas mass ratio is highly non-homogeneous throughout the disk extent, showing order-of-magnitude local deviations from the canonical 1:100 value. Gravitationally bound clumps formed through gravitational fragmentation have a velocity pattern that deviates notably from the Keplerian rotation. Small dust is efficiently converted into grown dust in the clump interiors, reaching a maximum radius of several decimeters. Concurrently, grown dust drifts towards the clump center forming a massive compact central condensation (70-100 $M_\oplus$). We argue that protoplanets may form in the interiors of inward migrating clumps before they disperse through the action of tidal torques. We foresee the formation of protoplanets at orbital distances of several tens of au with initial masses of gas and dust in the protoplanetary seed in the (0.25-1.6) $M_{\rm Jup}$ and (1.0-5.5) $M_\oplus$ limits, respectively. The final masses of gas and dust in the protoplanets may however be much higher due to accretion from surrounding massive metal-rich disks/envelopes.

Episodic accretion plays an important role in the evolution of young stars. Although it has been under investigation for a long time, the origin of such episodic accretion events is not yet understood. We investigate the dust and gas emission of a sample of young outbursting sources in the infrared to get a better understanding of their properties and circumstellar material, and we use the results in a further work to model the objects. We used Herschel data, from our PI program of 12 objects and complemented with archival observations to obtain the spectral energy distributions (SEDs) and spectra of our targets. We report here the main characteristics of our sample, focussing on the SED properties and on the gas emission lines detected in the PACS and SPIRE spectra. The SEDs of our sample show the diversity of the outbursting sources, with several targets showing strong emission in the far-infrared from the embedded objects. Most of our targets reside in a complex environment, which we discuss in detail. We detected several atomic and molecular lines, in particular rotational CO emission from several transitions from J=38-37 to J=4-3. We constructed rotational diagrams for the CO lines, and derived in three domains of assumed local thermodynamic equilibrium (LTE) temperatures and column densities, ranging mainly between 0-100 K and 400-500K. We confirm correlation in our sample between intense CO $J=16-15$ emission and the column density of the warm domain of CO, N(warm). We notice a strong increase in luminosity of HH 381 IRS and a weaker increase for PP 13 S, which shows the beginning of an outburst.

The central region of a circumstellar disk is difficult to resolve in global numerical simulations of collapsing cloud cores, but its effect on the evolution of the entire disk can be significant. We use numerical hydrodynamics simulations to model the long-term evolution of self-gravitating and viscous circumstellar disks in the thin-disk limit. Simulations start from the gravitational collapse of prestellar cores of 0.5--1.0~$M_\odot$ and both gaseous and dusty subsystems were considered, including a model for dust growth. The inner unresolved 1.0 au of the disk is replaced with a central "smart" cell (CSC) -- a simplified model that simulates physical processes that may occur in this region. We found that the mass transport rate through the CSC has an appreciable effect on the evolution of the entire disk. Models with slow mass transport form more massive and warmer disks and they are more susceptible to gravitational instability and fragmentation, including a newly identified episodic mode of disk fragmentation in the T Tauri phase of disk evolution. Models with slow mass transport through the CSC feature episodic accretion and luminosity bursts in the early evolution, while models with fast transport are characterized by a steadily declining accretion rate with low-amplitude flickering. Dust grows to a larger, decimeter size in the slow transport models and efficiently drifts in the CSC, where it accumulates reaching the limit when streaming instability becomes operational. We argue that gravitational instability, together with streaming instability likely operating in the inner disk regions, constitute two concurrent planet-forming mechanisms, which may explain the observed diversity of exoplanetary orbits (Abridged).

We report new Atacama Large Millimeter/submillimeter Array Band 3 (86-100 GHz; $\sim$80 mas angular resolution) and Band 4 (146-160 GHz; $\sim$50 mas angular resolution) observations of the dust continuum emission towards the archetypal and ongoing accretion burst young stellar object FU Ori, which simultaneously covered its companion, FU Ori S. In addition, we present near-infrared (2-2.45 $\mu$m) observations of FU Ori taken with the General Relativity Analysis via VLT InTerferometrY (GRAVITY; $\sim$1 mas angular resolution) instrument on the Very Large Telescope Interferometer (VLTI). We find that the emission in both FU Ori and FU Ori S at (sub)millimeter and near infrared bands is dominated by structures inward of $\sim$10 au radii. We detected closure phases close to zero from FU Ori with VLTI/GRAVITY, which indicate the source is approximately centrally symmetric and therefore is likely viewed nearly face-on. Our simple model to fit the GRAVITY data shows that the inner 0.4 au radii of the FU Ori disk has a triangular spectral shape at 2-2.45 $\mu$m, which is consistent with the H$_{2}$O and CO absorption features in a $\dot{M}\sim$10$^{-4}$ $M_{\odot}\,yr^{-1}$, viscously heated accretion disk. At larger ($\sim$0.4-10 au) radii, our analysis shows that viscous heating may also explain the observed (sub)millimeter and centimeter spectral energy distribution when we assume a constant, $\sim$10$^{-4}$ $M_{\odot}\,yr^{-1}$ mass inflow rate in this region. This explains how the inner 0.4 au disk is replenished with mass at a modest rate, such that it neither depletes nor accumulates significant masses over its short dynamic timescale. Finally, we tentatively detect evidence of vertical dust settling in the inner 10 au of the FU Ori disk, but confirmation requires more complete spectral sampling in the centimeter bands.

Global numerical simulations of protoplanetary disk formation and evolution were conducted in thin-disk limit, where the model included magnetically layered disk structure, a self-consistent treatment for the infall from cloud core as well as the smallest possible inner computational boundary. We compared the evolution of a layered disk with a fully magnetically active disk. We also studied how the evolution depends on the parameters of the layered disk model - the MRI triggering temperature and active layer thickness - as well as the mass of the prestellar cloud core. With the canonical values of parameters a dead zone formed within the inner $\approx$ 15 au region of the magnetically layered disk. The dead zone was not a uniform structure and long-lived, axisymmetric, gaseous rings ubiquitously formed within this region due to the action of viscous torques. The rings showed a remarkable contrast in the disk environment as compared to a fully magnetically active disk and were characterized by high surface density and low effective viscosity. Multiple gaseous rings could form simultaneously in the dead zone region which were highly dynamical and showed complex, time-dependent behavior such as inward migration, vortices, gravitational instability and large-scale spiral waves. An increase in MRI triggering temperature had only marginal effects, while changes in active layer thickness as well as the initial cloud core mass had significant effects on the structure and evolution of the inner disk. Dust with large fragmentation barrier could be trapped in the rings, which may play a key role in planet formation.

The accretion-driven outbursts of young FU Orionis-type stars may be a common stage of pre-main sequence evolution and can have a significant impact on the circumstellar environment as it pertains to the growth of solids and eventually planets. This episodic accretion is thought to be sustained by additional gas in-falling from the circumstellar envelope and disk. We present APEX observations of the CO gas in the envelope around V883 Orionis, a young outbursting star. The observations mapped the $^{12}$CO(4-3), $^{12}$CO(3-2), and $^{13}$CO(3-2) lines with the FLASH$^{+}$ instrument and the $^{12}$CO(6-5) line with the SEPIA instrument. We detected high signal-to-noise emission extending out to radii $>10000$ au and calculated integrated fluxes of $1100~\rm Jy~km~s^{-1}$ for $^{12}$CO(6-5), $2400~\rm Jy~km~s^{-1}$ for $^{12}$CO(4-3), $1600~\rm Jy~km~s^{-1}$ for $^{12}$CO(3-2), and $450~\rm Jy~km~s^{-1}$ for $^{13}$CO(3-2). We used the thermo-chemical code P{\small RO}D{\small I}M{\small O} to test several models and find the data are best described by an envelope structure with $\rm M_{env}\approx 0.2-0.4\,M_{\odot}$ and a mass-infall rate of $\rm \dot{M}_{inf}=1-2\times10^{-6}\,M_{\odot}\,yr^{-1}$. We infer that the observed envelope and outflow structure around V883 Ori could be caused by multiple outbursts, consistent with episodic accretion.

Dust is known to drift and grow in protoplanetary discs, which results in dust segregation over the disc extent. Maps of the spectral index $\alpha$ are a common tool for studying the dust content in protoplanetary discs. The analysis of observationally derived maps reveals significant gradients of the spectral index, confirming that dust evolves in the disc, but a more detailed information about the dust redistribution is required to make inferences about the early stages of dust growth. We calculated the spectral index maps based on the results of numerical hydrodynamical simulations using the FEOSAD code, which allows studying a long-term dynamics of a self-gravitating viscous disc populated with coagulating, drifting, and fragmenting dust. Here we demonstrate that values of the spectral index estimated for different wavelength intervals within the far-infrared and radio bands reveal the presence of dust grains of various sizes. Specifically, we show that the disc regions with the maximal spectral index in a specific wavelength interval are the regions with the prevalence of dust grains of a specific size. Thus, a set of spectral index maps derived using different wavelength intervals can be used to recover the dust size-distribution over the disc extent.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the CO ($J$=2--1) line emission from the protoplanetary disk around T-Tauri star SU Aurigae (hereafter SU Aur). Previous observations in optical and near infrared wavelengths find a unique structure in SU Aur. One of the highlights of the observational results is that an extended tail-like structure is associated with the disk, indicating mass transfer from or into the disk. Here we report the discovery of the counterpart of the tail-like structure in CO gas extending more than 1000 au long. Based on geometric and kinematic perspectives, both of the disk and the tail-like structure components physically connect to each other. Several theoretical studies predicted the observed tail-like structure via the following possible scenarios, 1) a gaseous stream from the molecular cloud remnant, 2) collision with a (sub)stellar intruder or a gaseous blob from the ambient cloud, and 3) ejection of a planetary or brown dwarf mass object due to gravitational instability via multi-body gravitational interaction. Since the tail-like structures associated with the SU Aur disk is a new example following RW Aurigae, some disks may experience the internal or external interaction and drastically lose mass during disk evolution.

The chemical influence of luminosity outbursts on the environments of young solar-type stars is explored. Species are categorised into several types according to their response to the outburst. The first and second types imply chemical changes only during the outburst (with slightly different behaviours). These response types are mostly observed close to the star and are caused by icy mantle evaporation. However, mantles recover after the outburst almost immediately. A notable exception is benzene ice, which is accumulated on dust surfaces during and after the outburst, so that its abundance exceeds the pre-outburst level by orders of magnitude. The third type of response is mostly seen at the disc periphery and implies alteration of abundances during the outburst and preservation of these `abnormal' abundances for centuries. This behaviour is typical of organic compounds, like HCOOCH$_3$, CH$_3$CN, CH$_2$CO. Their presence in the dark disc regions can be a manifestation of the past outburst. CO and CO$_2$ only trace past outbursts at the remote disc regions. The outburst changes the C/O ratio, but it quickly returns to the pre-outburst value almost everywhere in the disc. An important factor determining the sensitivity of molecular composition to the outburst is the dust size distribution. The duration of the pre-outburst stage and of the outburst itself influence the chemical effects, if the burst duration is shorter than 50 yr and the duration of the quiescent phase between the bursts is shorter than 100 kyr.

Following our recent work devoted to the effect of accretion on the pre-main-sequence evolution of low-mass stars, we perform a detailed analysis of episodic excursions of low-mass protostars in the Hertzsprung-Russell (H-R) diagram triggered by strong mass accretion bursts typical of FU Orionis-type objects (FUors). These excursions reveal themselves as sharp increases in the stellar total luminosity and/or effective temperature of the protostar and can last from hundreds to a few thousands of years, depending on the burst strength and characteristics of the protostar. During the excursions, low-mass protostars occupy the same part of the H-R diagram as young intermediate-mass protostars in the quiescent phase of accretion. Moreover, the time spent by low-mass protostars in these regions is on average a factor of several longer than that spent by the intermediate-mass stars in quiescence. During the excursions, low-mass protostars pass close to the position of most known FUors in the H-R diagram, but owing to intrinsic ambiguity the model stellar evolutionary tracks are unreliable in determining the FUor properties. We find that the photospheric luminosity in the outburst state may dominate the accretion luminosity already after a few years after the onset of the outburst, meaning that the mass accretion rates of known FUors inferred from the bolometric luminosity may be systematically overestimated, especially in the fading phase.

Massive protostars grow and evolve under the effect of rapid accretion of circumstellar gas and dust, falling at high rates ($\ge 10^{-4}$-$10^{-3}\, \rm M_{\odot}\, \rm yr^{-1}$). This mass infall has been shown, both numerically and observationally, to be episodically interspersed by accretion of dense gaseous clumps migrating through the circumstellar disc to the protostellar surface, causing sudden accretion and luminous bursts. Using numerical gravito-radiation-hydrodynamics and stellar evolution calculations, we demonstrate that, in addition to the known bloating of massive protostars, variable episodic accretion further influences their evolutionary tracks of massive young stellar objects (MYSOs). For each accretion-driven flare, they experience rapid excursions toward more luminous, but colder regions of the Hertzsprung-Russell diagram. During these excursions, which can occur up to the end of the pre-main-sequence evolution, the photosphere of massive protostars can episodically release much less energetic photons and MYSOs surreptitiously adopt the same spectral type as evolved massive (supergiants) stars. Each of these evolutionary loop brings the young high-mass stars close to the forbidden Hayashi region and might make their surrounding H II regions occasionally fainter, before they recover their quiescent, pre-burst surface properties. We interpret such cold, intermittent pre-main-sequence stellar evolutionary excursions and the dipping variability of HII regions as the signature of the presence of a fragmenting circumstellar accretion disc surrounding the MYSOs. We conjecture that this mechanism might equivalently affect young stars in the intermediate-mass regime.

Circumstellar discs, from which planetary systems are formed, consist of gas, dust and solids. Simulations of self-consistent dynamics of gas, dust and solids in circumstellar discs is a challenging problem. In the paper we present fast algorithms for computing the drag force (momentum transfer) between solid phase and gas. These algorithms (a) are universal and applicable to dust and solids with any sizes smaller than the mean free path of gas molecules, (b) can be used to calculate the momentum transfer between dust and gas instead of one-way effect, as it is done in many models, (c) can perform simulations, without a loss in accuracy, with the time step determined by gas-dynamic parameters rather than by drag force, and (d) are compatible with the widely used parallel algorithms for solving 3D equations of gas dynamics, hydrodynamic equations for dust, and the collisionless Boltzmann equation for large bodies. Preliminary results of supercomputer simulation of the gas-dust disc dynamics within the developed approach are reported.

Aims: We aim to use the high spatial resolution of the Atacama Large Millimeter/submillimeter Array (ALMA) to map the flow pattern of molecular gas near DG Tau and its disk, a young stellar object driving a jet and a molecular outflow. Methods: We use observations from ALMA in the J = 2 - 1 transition of 12CO, 13CO, and C18O to study the Keplerian disk of DG Tau and outflows that may be related to the disk and the jet. Results: We find a new wind component flowing radially at a steep angle (~25 deg from the vertical) above the disk with a velocity of ~ 3.1 km/s. It continues the trend of decreasing velocity for increasing distance from the jet axis ("onion-like velocity structure"). Conclusions: The new component is located close to the protostellar disk surface and may be related to photoevaporative winds.

Episodic accretion may be a common occurrence in the evolution of young pre-main sequence stars and has important implications for our understanding of star and planet formation. Many fundamental aspects of what drives the accretion physics, however, are still unknown. The ngVLA will be a key tool in understanding the nature of these events. The high spatial resolution, broad spectral coverage, and unprecedented sensitivity will allow for the detailed analysis of outburst systems. The proposed frequency range of the ngVLA allows for observations of the gas, dust, and non-thermal emission from the star and disk.

The FU Ori-type young stellar objects are characterized by a sudden increase in luminosity by 1$-$2 orders of magnitude, followed by slow return to the pre-outburst state on timescales of $\sim$10$-$100 yr. The outburst strongly affects the entire disk, changing its thermal structure and radiation field. In this paper, using a detailed physical-chemical model we study the impact of the FU Ori outburst on the disk chemical inventory. Our main goal is to identify gas-phase molecular tracers of the outburst activity that could be observed after the outburst with modern telescopes such as ALMA and NOEMA. We find that the majority of molecules experience a considerable increase in the total disk gas-phase abundances due to the outburst, mainly due to the sublimation of their ices. Their return to the pre-outburst chemical state takes different amounts of time, from nearly instantaneous to very long. Among the former ones we identify CO, NH$_3$, C$_2$H$_6$, C$_3$H$_4$, etc. Their abundance evolution tightly follows the luminosity curve. For CO the abundance increase does not exceed an order of magnitude, while for other tracers the abundances increase by 2$-$5 orders of magnitude. Other molecules like H$_2$CO and NH$_2$OH have longer retention timescales, remaining in the gas phase for $\sim10-10^3$ yr after the end of the outburst. Thus H$_2$CO could be used as an indicator of the previous outbursts in the post-outburst FU Ori systems. We investigate the corresponding time-dependent chemistry in detail and present the most favorable transitions and ALMA configurations for future observations.

A systematic analysis of methods for computing the trajectories of solid-phase particles applied in modern astrophysics codes designed for modeling gas-dust circumstellar disks has been carried out for the first time. The motion of grains whose velocities are determined mainly by the gas drag, that is, for which the stopping time or relaxation time for the velocity of the dust to the velocity of the gas t_{stop} is less than or comparable to the rotation period, are considered. The methods are analyzed from the point of view of their stability for computing the motions of small bodies, including dust grains less than 1 micron, which are strongly coupled to the gas. Two test problems are with analytical solutions. Fast first order accurate methods that make it possible to avoid additional restriction on the time step due to gas drag in computations of the motion of grains of any size are presented. For the conditions of a circumstellar disks, the error in the velocity computations obtained when using some stable methods become unacceptably large when the time step size \tau>t_{stop}. For the radial migration of bodies that exhibit drifts along nearly Keplerian orbits, an asymptotic approximation, sometimes called the short friction time approximation or drift flux model, gives a relative error for the radial velocity equals to St^2 where St is the Stokes number, the ratio of the stopping time of the body to some fraction of the rotation period (dynamical time scale) in the disk.

We present near-IR imaging polarimetry of five classical FU Ori-type objects (FU Ori, V1057 Cyg, V1515 Cyg, V1735 Cyg, Z CMa) with a $\sim$0\farcs1 resolution observed using HiCIAO+AO188 at Subaru Telescope. We observed scattered light associated with circumstellar dust around four of them (i.e., all but V1515 Cyg). Their polarized intensity distribution shows a variety of morphologies with arms, tails or streams, spikes and fragmented distributions, many of which were reported in our previous paper. The morphologies of these reflection nebulae significantly differ from many other normal young stellar objects (Class I-II objects). These structures are attributed to gravitationally unstable disks, trails of clump ejections, dust blown by a wind or a jet, and a stellar companion. We can consistently explain our results with the scenario that their accretion outbursts (FUor outbursts) are triggered by gravitationally fragmenting disks, and with the hypothesis that many low-mass young stellar objects experience such outbursts.

Approaches used in modern numerical simulations of the dynamics of dust and gas in circumstellar disks are tested. The gas and dust are treated like interpenetrating continuous media that can exchange momentum. A stiff coupling between the gas and dust phases is typical for such disks, with the dust stopping time much less than the characteristic dynamical time scale. This imposes high demands on the methods used to simulate the dust dynamics. A grid, piecewise-parabolic method is used as the basic algorithm for solving the gas-dynamical equations. Numerical solutions obtained using various methods to compute the momentum exchanges are presented for the case of monodisperse dust. Numerical solutions are obtained for shock tube problem and the propagation of sound waves in a gas-dust medium. The studied methods are compared in terms of their ability to model media with (a) an arbitrary (short or long) dust stopping time, and (b) an arbitrary dust concentration in the gas (varying the dust to gas mass ratio from 0.01 to 1). A method for computing the momentum exchange with infinite-order accuracy in time is identified, which makes it possible to satisfy the conditions (a) and (b) with minimal computational costs. A first-order method that shows similar results in the test computations is also presented. It is shown that the proposed first-order method for monodisperse dust can be extended to a regime when the dust is polydisperse; i.e., a regime represented by several fractions with different stopping times. Formulas for computing the gas and dust velocities for polydisperse dust with each fraction exchanging momentum with the gas are presented.

Migration of dense gaseous clumps that form in young protostellar disks via gravitational fragmentation is investigated to determine the likelihood of giant planet formation. High-resolution numerical hydrodynamics simulations in the thin-disk limit are employed to compute the formation and long-term evolution of a gravitationally unstable protostellar disk around a solar-mass star. We show that gaseous clumps that form in the outer regions of the disk (>100 AU) through disk fragmentation are often perturbed by other clumps or disk structures, such as spiral arms, and migrate toward the central star on timescales from a few 10^3 to few 10^4 yr. The migration timescale is slowest when stellar motion in response to the disk gravity is considered. When approaching the star, the clumps first gain mass (up to several tens of M_Jup), but then quickly lose most of their diffuse envelopes through tidal torques. Part of the clump envelope can be accreted on the central star causing an FU-Ori-type accretion and luminosity outburst. The tidal mass loss helps the clumps to significantly slow down or even halt their inward migration at a distance of a few tens of AU from the protostar. The resulting clumps are heavily truncated both in mass and size compared to their wider-orbit counterparts, keeping only a dense and hot nucleus. During the inward migration, the temperature in the clump interiors may exceed the molecular hydrogen dissociation limit (2000 K) and the central region of the clump can collapse into a gas giant protoplanet. We argue that FU-Orionis-type luminosity outbursts may be the end product of disk fragmentation and clump inward migration, ushering the formation of giant protoplanets on tens-of-au orbits in systems such as HR~8799.

We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we can characterize the evolution of protostars. This paper releases the MASSES 1.3 mm data from the subcompact configuration ($\sim$4$^{\prime\prime}$ or $\sim$1000 au resolution), which is the SMA's most compact array configuration. We release both $uv$ visibility data and imaged data for the spectral lines CO(2-1),$^{13}$CO(2-1), C$^{18}$O(2-1), and N$_2$D$^+$(3-2), as well as for the 1.3 mm continuum. We identify the tracers that are detected toward each source. We also show example images of continuum and CO(2-1) outflows, analyze C$^{18}$O(2-1) spectra, and present data from the SVS 13 star-forming region. The calculated envelope masses from the continuum show a decreasing trend with bolometric temperature (a proxy for age). Typical C$^{18}$O(2-1) linewidths are 1.45 km s$^{-1}$, which is higher than the C$^{18}$O linewidths detected toward Perseus filaments and cores. We find that N$_2$D$^+$(3-2) is significantly more likely to be detected toward younger protostars. We show that the protostars in SVS 13 are contained within filamentary structures as traced by C$^{18}$O(2-1) and N$_2$D$^+$(3-2). We also present the locations of SVS 13A's high velocity (absolute line-of-sight velocities $>$150 km s$^{-1}$) red and blue outflow components. Data can be downloaded from https://dataverse.harvard.edu/dataverse/MASSES .

The long-term evolution of a circumstellar disk starting from its formation and ending in the T Tauri phase was simulated numerically with the purpose of studying the evolution of dust in the disk with distinct values of viscous \alpha-parameter and dust fragmentation velocity v_frag. We solved numerical hydrodynamics equations in the thin-disk limit, which are modified to include a dust component consisting of two parts: sub-micron-sized dust and grown dust with a maximum radius a_r. The former is strictly coupled to the gas, while the latter interacts with the gas via friction. The conversion of small to grown dust, dust growth, and dust self-gravity are also considered. We found that the process of dust growth known for the older protoplanetary phase also holds for the embedded phase of disk evolution. The dust growth efficiency depends on the radial distance from the star - a_r is largest in the inner disk and gradually declines with radial distance. In the inner disk, a_r is limited by the dust fragmentation barrier. The process of small-to-grown dust conversion is very fast once the disk is formed. The total mass of grown dust in the disk (beyond 1 AU) reaches tens or even hundreds of Earth masses already in the embedded phase of star formation and even a greater amount of grown dust drifts in the inner, unresolved 1 AU of the disk. Dust does not usually grow to radii greater than a few cm. A notable exception are models with \alpha <= 10^{-3}, in which case a zone with reduced mass transport develops in the inner disk and dust can grow to meter-sized boulders in the inner 10 AU. Grown dust drifts inward and accumulates in the inner disk regions. This effect is most pronounced in the \alpha <= 10^{-3} models where several hundreds of Earth masses can be accumulated in a narrow region of several AU from the star by the end of embedded phase. (abridged).

We aim at studying the causal link between the knotty jet structure in CARMA 7, a young Class 0 protostar in the Serpens South cluster, and episodic accretion in young protostellar disks. We used numerical hydrodynamics simulations to derive the protostellar accretion history in gravitationally unstable disks around solar-mass protostars. We compared the time spacing between luminosity bursts \Delta\tau_mod, caused by dense clumps spiralling on the protostar, with the differences of dynamical timescales between the knots \Delta\tau_obs in CARMA 7. We found that the time spacing between the bursts have a bi-modal distribution caused by isolated and clustered luminosity bursts. The former are characterized by long quiescent periods between the bursts with \Delta\tau_mod = a few * (10^3-10^4) yr, whereas the latter occur in small groups with time spacing between the bursts \Delta\tau_mod= a few * (10-10^2) yr. For the clustered bursts, the distribution of \Delta\tau_mod in our models can be fit reasonably well to the distribution of \Delta\tau_obs in the protostellar jet of CARMA 7, if a certain correction for the (yet unknown) inclination angle with respect to the line of sight is applied. The K-S test on the model and observational data sets suggests the best-fit values for the inclination angles of 55-80 deg., which become narrower (75-80 deg.) if only strong luminosity bursts are considered. The dynamical timescales of the knots in the jet of CARMA 7 are too short for a meaningful comparison with the long time spacings between isolated bursts in our models. The exact sequences of time spacings between the luminosity bursts in our models and knots in the jet of CARMA 7 were found difficult to match. (abridged)

Young stellar objects (YSOs) may undergo periods of active accretion (outbursts), during which the protostellar accretion rate is temporarily enhanced by a few orders of magnitude. Whether or not these accretion outburst YSOs possess similar dust/gas reservoirs to each other, and whether or not their dust/gas reservoirs are similar as quiescent YSOs, are issues not yet clarified. The aim of this work is to characterize the millimeter thermal dust emission properties of a statistically significant sample of long and short duration accretion outburst YSOs (i.e., FUors and EXors) and the spectroscopically identified candidates of accretion outbursting YSOs (i.e., FUor-like objects). We have carried out extensive Submillimeter Array (SMA) observations mostly at $\sim$225 GHz (1.33 mm) and $\sim$272 GHz (1.10 mm), from 2008 to 2017. We covered accretion outburst YSOs located at $<$1 kpc distances from the solar system. We analyze all the existing SMA data of such objects, both published and unpublished, in a coherent way to present a millimeter interferometric database of 29 objects. We obtained 21 detections at $>$3-$\sigma$ significance. Detected sources except for the two cases of V883 Ori and NGC 2071 MM3 were observed with $\sim$1$"$ angular resolution. Overall our observed targets show a systematically higher millimeter luminosity distribution than those of the $M_{*}>$0.3 $M_{\odot}$ Class II YSOs in the nearby ($\lesssim$400 pc) low-mass star-forming molecular clouds (e.g., Taurus, Lupus, Upp Scorpio, and Chameleon I). In addition, at 1 mm our observed confirmed binaries or triple-system sources are systematically fainter than the rest of the sources even though their 1 mm fluxes are broadly distributed. We may have detected $\sim$30-60\% millimeter flux variability from V2494 Cyg and V2495 Cyg, from the observations separated by $\sim$1 year.

Using the HPC ressources of the state of Baden-W\"urttemberg, we modelled for the first time the luminous burst from a young massive star by accretion of material from its close environment. We found that the surroundings of young massive stars are shaped as a clumpy disk whose fragments provoke outbursts once they fall onto the protostar and concluded that similar strong luminous events observed in high-mass star forming regions may be a signature of the presence of such disks.

The surroundings of massive protostars constitute an accretion disc which has numerically been shown to be subject to fragmentation and responsible for luminous accretion-driven outbursts. Moreover, it is suspected to produce close binary companions which will later strongly influence the star's future evolution in the Hertzsprung-Russel diagram. We present three-dimensional gravitation-radiation-hydrodynamic numerical simulations of 100 Mo pre-stellar cores. We find that accretion discs of young massive stars violently fragment without preventing the (highly variable) accretion of gaseous clumps onto the protostars. While acquiring the characteristics of a nascent low-mass companion, some disc fragments migrate onto the central massive protostar with dynamical properties showing that its final Keplerian orbit is close enough to constitute a close massive proto-binary system, having a young high-mass and a low-mass component. We conclude on the viability of the disc fragmentation channel for the formation of such short-period binaries, and that both processes -close massive binary formation and accretion bursts- may happen at the same time. FU-Orionis-type bursts, such as observed in the young high-mass star S255IR-NIRS3, may not only indicate ongoing disc fragmentation, but also be considered as a tracer for the formation of close massive binaries - progenitors of the subsequent massive spectroscopic binaries - once the high-mass component of the system will enter the main-sequence phase of its evolution. Finally, we investigate the ALMA-observability of the disc fragments.

Horseshoe-shaped brightness asymmetries of several transitional discs are thought to be caused by large-scale vortices. Anticyclonic vortices are efficiently collect dust particles, therefore they can play a major role in planet formation. Former studies suggest that the disc self-gravity weakens vortices formed at the edge of the gap opened by a massive planet in discs whose masses are in the range of 0.01<=M_disc/M_*<=0.1. Here we present an investigation on the long-term evolution of the large-scale vortices formed at the viscosity transition of the discs' dead zone outer edge by means of two-dimensional hydrodynamic simulations taking disc self-gravity into account. We perform a numerical study of low mass, 0.001<=M_disc/M_*<=0.01, discs, for which cases disc self-gravity was previously neglected. The large-scale vortices are found to be stretched due to disc self-gravity even for low-mass discs with M_disc/M_*>=0.005 where initially the Toomre Q-parameter was <=50 at the vortex distance. As a result of stretching, the vortex aspect ratio increases and a weaker azimuthal density contrast develops. The strength of the vortex stretching is proportional to the disc mass. The vortex stretching can be explained by a combined action of a non-vanishing gravitational torque caused by the vortex, and the Keplerian shear of the disc. Self-gravitating vortices are subject to significantly faster decay than non-self-gravitating ones. We found that vortices developed at sharp viscosity transitions of self-gravitating discs can be described by a GNG model as long as the disc viscosity is low, i.e. alpha_dz<=10^-5.

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.

A long-standing open issue of the paradigm of low-mass star formation is the luminosity problem: most protostars are less luminous than theoretically predicted. One possible solution is that the accretion process is episodic. FU Ori-type stars (FUors) are thought to be the visible examples for objects in the high accretion state. FUors are often surrounded by massive envelopes, which replenish the disk material and enable the disk to produce accretion outbursts. However, we have insufficient information on the envelope dynamics in FUors, about where and how mass transfer from the envelope to the disk happens. Here we present ALMA observations of the FUor-type star V346 Nor at 1.3 mm continuum and in different CO rotational lines. We mapped the density and velocity structure of its envelope and analyze the results using channel maps, position-velocity diagrams, and spectro-astrometric methods. We found that V346 Nor is surrounded by gaseous material on 10000 au scale in which a prominent outflow cavity is carved. Within the central $\sim$700 au, the circumstellar matter forms a flattened pseudo-disk where material is infalling with conserved angular momentum. Within $\sim$350 au, the velocity profile is more consistent with a disk in Keplerian rotation around a central star of 0.1 $M_{\odot}$. We determined an infall rate from the envelope onto the disk of 6$\times$10$^{-6}\,M_{\odot}$yr$^{-1}$, a factor of few higher than the quiescent accretion rate from the disk onto the star, hinting for a mismatch between the infall and accretion rates as the cause of the eruption.

We present new Submillimeter Array (SMA) observations of CO(2-1) outflows toward young, embedded protostars in the Perseus molecular cloud as part of the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. For 57 Perseus protostars, we characterize the orientation of the outflow angles and compare them with the orientation of the local filaments as derived from $Herschel$ observations. We find that the relative angles between outflows and filaments are inconsistent with purely parallel or purely perpendicular distributions. Instead, the observed distribution of outflow-filament angles are more consistent with either randomly aligned angles or a mix of projected parallel and perpendicular angles. A mix of parallel and perpendicular angles requires perpendicular alignment to be more common by a factor of $\sim$3. Our results show that the observed distributions probably hold regardless of the protostar's multiplicity, age, or the host core's opacity. These observations indicate that the angular momentum axis of a protostar may be independent of the large-scale structure. We discuss the significance of independent protostellar rotation axes in the general picture of filament-based star formation.

The pre-main-sequence evolution of low-mass stars and brown dwarfs is studied numerically starting from the formation of a protostellar/proto-brown dwarf seed and taking into account the mass accretion onto the central object during the initial several Myr of evolution. The stellar evolution was computed using the STELLAR evolution code developed by Yorke & Bodenheimer with recent modifications by Hosokawa et al. The mass accretion rates were taken from numerical hydrodynamics models of Vorobyov & Basu computing the circumstellar disk evolution starting from the gravitational collapse of pre-stellar cloud cores of various mass and angular momentum. The resulting stellar evolution tracks were compared with the isochrones and isomasses calculated using non-accreting models. We find that mass accretion in the initial several Myr of protostellar evolution can have a strong effect on the subsequent evolution of young stars and brown dwarfs. The disagreement between accreting and non-accreting models in terms of the total stellar luminosity L_st, stellar radius R_st and effective temperature T_eff depends on the thermal efficiency of accretion, i.e., on the fraction of accretion energy absorbed by the central object. The largest mismatch is found for the cold accretion case, in which essentially all accretion energy is radiated away. The relative deviations in L_st and R_st in this case can reach 50% for 1.0-Myr-old objects and remain notable even for 10-Myr-old objects. In the hot and hybrid accretion cases, in which a constant fraction of accretion energy is absorbed, the disagreement between accreting and non-accreting models becomes less pronounced, but still remains notable for 1.0-Myr-old objects. These disagreements may lead to the wrong age estimate for objects of (sub-)solar mass when using the isochrones based on non-accreting models (abridged).

Episodic accretion is an important process in the evolution of young stars and their environment. The observed strong luminosity bursts of young stellar objects likely have a long lasting impact on the chemical evolution of the disk and envelope structure. We want to investigate observational signatures of the chemical evolution in the post-burst phase for embedded sources. With such signatures it is possible to identify targets that experienced a recent luminosity burst. We present a new model for episodic accretion chemistry based on the 2D, radiation thermo-chemical disk code ProDiMo. We have extended ProDiMo with a proper treatment for envelope structures. For a representative Class I model, we calculated the chemical abundances in the post-burst phase and produced synthetic observables like intensity maps and radial profiles. During a burst many chemical species, like CO, sublimate from the dust surfaces. As the burst ends they freeze out again (post-burst phase). This freeze-out happens from inside-out due to the radial density gradient in the disk and envelope structure. This inside-out freeze-out produces clear observational signatures in spectral line emission, like rings and distinct features in the slope of radial intensity profiles. We fitted synthetic C18O J=2-1 observations with single and two component fits and find that post-burst images are much better matched by the latter. Comparing the quality of such fits allows identification of post-burst targets in a model-independent way. Our models confirm that it is possible to identify post-burst objects from spatially resolved CO observations. However, to derive proper statistics, like frequencies of bursts, from observations it is important to consider aspects like the inclination and structure of the target and also dust properties as those have a significant impact on the freeze-out timescale.

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.

Context: Understanding how accretion proceeds is a key question of star formation, with important implications for both the physical and chemical evolution of young stellar objects. In particular, very little is known about the accretion variability in the earliest stages of star formation. Aims: To characterise protostellar accretion histories towards individual sources by utilising sublimation and freeze-out chemistry of CO. Methods: A sample of 24 embedded protostars are observed with the Submillimeter Array (SMA) in context of the large program "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). The size of the C$^{18}$O emitting region, where CO has sublimated into the gas-phase, is measured towards each source and compared to the expected size of the region given the current luminosity. The SMA observations also include 1.3 mm continuum data, which are used to investigate whether a link can be established between accretion bursts and massive circumstellar disks. Results: Depending on the adopted sublimation temperature of the CO ice, between 20% and 50% of the sources in the sample show extended C$^{18}$O emission indicating that the gas was warm enough in the past that CO sublimated and is currently in the process of refreezing; something which we attribute to a recent accretion burst. Given the fraction of sources with extended C$^{18}$O emission, we estimate an average interval between bursts of 20000-50000 yr, which is consistent with previous estimates. No clear link can be established between the presence of circumstellar disks and accretion bursts, however the three closest known binaries in the sample (projected separations <20 AU) all show evidence of a past accretion burst, indicating that close binary interactions may also play a role in inducing accretion variability.

We resolved FU Ori at 29-37 GHz using the JVLA with $\sim$0$''$.07 resolution, and performed the complementary JVLA 8-10 GHz observations, the SMA 224 GHz and 272 GHz observations, and compared with archival ALMA 346 GHz observations to obtain the SEDs. Our 8-10 GHz observations do not find evidence for the presence of thermal radio jets, and constrain the radio jet/wind flux to at least 90 times lower than the expected value from the previously reported bolometric luminosity-radio luminosity correlation. The emission at $>$29 GHz may be dominated by the two spatially unresolved sources, which are located immediately around FU Ori and its companion FU Ori S, respectively. Their deconvolved radii at 33 GHz are only a few au. The 8-346 GHz SEDs of FU Ori and FU Ori S cannot be fit by constant spectral indices (over frequency). The more sophisticated models for SEDs suggest that the $>$29 GHz emission is contributed by a combination of free-free emission from ionized gas, and thermal emission from optically thick and optically thin dust components. We hypothesize that dust in the innermost parts of the disks ($\lesssim$0.1 au) has been sublimated, and thus the disks are no more well shielded against the ionizing photons. The estimated overall gas and dust mass based on SED modeling, can be as high as a fraction of a solar mass, which is adequate for developing disk gravitational instability. Our present explanation for the observational data is that the massive inflow of gas and dust due to disk gravitational instability or interaction with a companion/intruder, was piled up at the few au scale due to the development of a deadzone with negligible ionization. The piled up material subsequently triggered the thermal and the MRI instabilities when the ionization fraction in the inner sub-au scale region exceeded a threshold value, leading to the high protostellar accretion rate.

Grid-based hydrodynamics simulations of circumstellar disks are often performed in the curvilinear coordinate system, in which the center of the computational domain coincides with the motionless star. However, the center of mass may be shifted from the star due to the presence of any non-axisymmetric mass distribution. As a result, the system exerts a gravity force on the star, causing the star to move in response, which can affect the evolution of the circumstellar disk. We aim at studying the effects of stellar motion on the evolution of protostellar and protoplanetary disks. In protostellar disks, a non-axisymmetric distribution of matter in the form of spiral arms or massive clumps can form due to gravitational instability. Protoplanetary disks can also feature non-axisymmetric structures caused by a high-mass planet or a large-scale vortex. We use 2D grid-based hydrodynamic simulations to explore the effect of stellar motion. We adopt a non-inertial polar coordinate system centered on the star, in which the stellar motion is taken into account by calculating the indirect potential caused by the non-axisymmetric disk, a high-mass planet, or a large-scale vortex. We found that the stellar motion has a moderate effect on the evolution history in protostellar disks, reducing somewhat the disk size and mass, while having a profound effect on the collapsing envelope, changing its inner shape from an initially axisymmetric to a non-axisymmetric configuration. Protoplanetary disk simulations show that the stellar motion slightly reduces the width of the gap opened by a high-mass planet, decreases the planet migration rate, and strengthens the large-scale vortices formed at the viscosity transition. We conclude that the inclusion of the indirect potential is recommended in grid-based hydrodynamics simulations of circumstellar disks which use the curvilinear coordinate system.

The nature of very low luminosity objects with the internal luminosity L_obj <= 0.1 Lsun is investigated by means of numerical modeling coupling the core collapse simulations with the stellar evolution calculations. The gravitational collapse of a large sample of model cores in the mass range 0.1-2.0 Msun is investigated. Numerical simulations were started at the pre-stellar phase and terminated at the end of the embedded phase when 90% of the initial core mass had been accreted onto the forming protostar plus disk system. The disk formation and evolution was studied using numerical hydrodynamics simulations, while the formation and evolution of the central star was calculated using a stellar evolution code. Three scenarios for mass accretion from the disk onto the star were considered: hybrid accretion in which a fraction of accreted energy absorbed by the protostar depends on the accretion rate, hot accretion wherein a fraction of accreted energy is constant, and cold accretion wherein all accretion energy is radiated away. Our conclusions on the nature of VeLLOs depend crucially on the character of protostellar accretion. In the hybrid accretion scenario, most VeLLOs (90.6%) are expected to be the first hydrostatic cores (FHSCs) and only a small fraction (9.4%) are true protostars. In the hot accretion scenario, all VeLLOs are FHSCs due to overly high photospheric luminosity of protostars. In the cold accretion scenario, on the contrary, the majority of VeLLOs belong to the Class I phase of stellar evolution. The reason is that the stellar photospheric luminosity, which sets the floor for the total internal luminosity of a young star, is lower in cold accretion, thus enabling more VeLLOs in the protostellar stage. VeLLOs are relatively rare objects occupying 7%-11% of the total duration of the embedded phase and their masses do not exceed 0.3 Msun. (abridged).

Protoplanetary disks are the birthplaces of planetary systems. The evolution of the star-disk system and the disk chemical composition determines the initial conditions for planet formation. Therefore a comprehensive understanding of the main physical and chemical processes in disks is crucial for our understanding of planet formation. We give an overview of the early evolution of disks, discuss the importance of the stellar high-energy radiation for disk evolution and describe the general thermal and chemical structure of disks. Finally we provide an overview of observational tracers of the gas component and disk winds.

We present results from a deep sub-millimeter survey in the Serpens Main and Serpens/G3-G6 clusters, conducted with the Submillimetre Common-User Bolometer Array (SCUBA-2) at the James Clerk Maxwell Telescope. We have combined Herschel PACS far-infrared photometry, sub-millimeter continuum and molecular gas line observations, with the aim to conduct a detailed multi-wavelength characterization of `proto-brown dwarf' candidates in Serpens. We have performed continuum and line radiative transfer modeling, and have considered various classification schemes to understand the structure and the evolutionary stage of the system. We have identified four proto-brown dwarf candidates, of which the lowest luminosity source has an Lbol ~0.05 Lsun. Two of these candidates show characteristics consistent with Stage 0/I systems, while the other two are Stage I-T/Class Flat systems with tenuous envelopes. Our work has also revealed a ~20% fraction of mis-identified Class 0/I/Flat sources that show characteristics consistent with Class II edge-on disk systems. We have set constraints on the mass of the central object using the measured bolometric luminosities and numerical simulations of stellar evolution. Considering the available gas+dust mass reservoir and the current mass of the central source, three of these candidates are likely to evolve into brown dwarfs.

Accretion-driven luminosity outbursts are a vivid manifestation of variable mass accretion onto protostars. They are known as the so-called FU Orionis phenomenon in the context of low-mass protostars. More recently, this process has been found in models of primordial star formation. Using numerical radiation hydrodynamics simulations, we stress that present-day forming massive stars also experience variable accretion and show that this process is accompanied by luminous outbursts induced by the episodic accretion of gaseous clumps falling from the circumstellar disk onto the protostar. Consequently, the process of accretion-induced luminous flares is also conceivable in the high-mass regime of star formation and we propose to regard this phenomenon as a general mechanism that can affect protostars regardless of their mass and/or the chemical properties of the parent environment in which they form. In addition to the commonness of accretion-driven outbursts in the star formation machinery, we conjecture that luminous flares from regions hosting forming high-mass star may be an observational implication of the fragmentation of their accretion disks.

We present self-consistent calculations coupling numerical hydrodynamics simulations of collapsing pre-stellar cores and stellar evolution models of accreting objects. We analyse the main impact of consistent accretion history on the evolution and lithium depletion of young low-mass stars and brown dwarfs. These consistent models confirm the generation of a luminosity spread in Herzsprung-Russell diagrams at ages $\sim$ 1-10 Myr. They also confirm that early accretion can produce objects with abnormal Li depletion, as found in a previous study that was based on arbitrary accretion rates. The results strengthen that objects with anomalously high level of Li depletion in young clusters should be extremely rare. We also find that early phases of burst accretion can produce coeval models of similar mass with a range of different Li surface abundances, and in particular with Li-excess compared to the predictions of non-accreting counterparts. This result is due to a subtle competition between the effect of burst accretion and its impact on the central stellar temperature, the growth of the stellar radiative core and the accretion of fresh Li from the accretion disk. Only consistent models could reveal such a subtle combination of effects. This new result could explain the recent, puzzling observations of Li-excess of fast rotators in the young cluster NGC 2264. Present self-consistent accreting models are available in electronic form.

An approach for constructing a low-dissipation numerical method is described. The method is based on a combination of the operator-splitting method, Godunov method, and piecewise-parabolic method on the local stencil. Numerical method was tested on a standard suite of hydrodynamic test problems. In addition, the performance of the method is demonstrated on a global test problem showing the development of a spiral structure in a gravitationally unstable gaseous galactic disk.

Protostellar (class 0/I) disks, having masses comparable to those of their nascent host stars, and fed continuously from their natal infalling envelopes, are prone to gravitational instability (GI). Motivated by advances in near-infrared (NIR) adaptive optics imaging and mm-wave interferometry, we explore the observational signatures of GI in disks, using hydrodynamical and Monte Carlo radiative transfer simulations to synthesize NIR scattered light images and mm dust continuum maps. Spiral arms induced by GI, located at disk radii of hundreds of AUs, are local overdensities and have their photospheres displaced to higher altitudes above the disk midplane, arms therefore scatter more NIR light from their central stars than inter-arm regions, and are detectable at distances up to 1 kpc by Gemini/GPI, VLT/SPHERE, and Subaru/HiCIAO/SCExAO. By contrast, collapsed clumps formed by disk fragmentation have such strong local gravitational fields that their scattering photospheres are at lower altitudes, such fragments appear fainter than their surroundings in NIR scattered light. Spiral arms and streamers recently imaged in four FU Ori systems at NIR wavelengths resemble GI-induced structures and support the interpretation that FUors are gravitationally unstable protostellar disks. At mm wavelengths, both spirals and clumps appear brighter in thermal emission than the ambient disk and can be detected by ALMA at distances up to 0.4 kpc with one-hour integration times at ~0.1 arcsec resolution. Collapsed fragments having masses >=1 MJ can be detected by ALMA within ~10 minutes.

We investigate the dynamics of gaseous clumps formed via gravitational fragmentation in young protostellar disks, focusing on the fragments that are ejected from the disk via many-body gravitational interaction. Numerical hydrodynamics simulations were employed to study the evolution of young protostellar disks formed from the collapse of rotating pre-stellar cores with mass in the 1.1-1.6 M_sun range. Protostellar disks formed in our models undergo gravitational fragmentation driven by continuing mass loading from parental collapsing cores. A few fragments can be ejected from the disk during the early evolution, but the low-mass fragments (< 15~M_Jup) disperse creating spectacular bow-type structures while passing through the disk and collapsing core. The least massive fragment that survived the ejection (21 M_Jup) straddles the planetary-mass limit, while the most massive ejected fragments (145 M_Jup) can break up into several pieces, leading to the ejection of wide separation binary clumps in the brown-dwarf mass range. About half of the ejected fragments are gravitationally bound, the majority is supported by rotation against gravity, and all fragments have the specific angular momentum that is much higher than that expected for brown dwarfs. We found that the internal structure of the ejected fragments is distinct from what would be expected for gravitationally contracting clumps formed via cloud core fragmentation, which can help to differentiate their origin. The ejection of fragments is an important process inherent to massive protostellar disks, which produces freely-floating pre-brown dwarf cores, regulates the disk and stellar masses, and potentially enriches the intracluster medium with processed dust and complex organics.

We investigate the alignment between outflow axes in nine of the youngest binary/multiple systems in the Perseus Molecular Cloud. These systems have typical member spacing larger than 1000 AU. For outflow identification, we use 12CO(2-1) and 12CO(3-2) data from a large survey with the Sub- millimeter Array: Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES). The distribution of outflow orientations in the binary pairs is consistent with random or preferentially anti-aligned distributions, demonstrating that these outflows are misaligned. This result suggests that these systems are possibly formed in environments where the distribution of angular momentum is complex and disordered, and these systems do not come from the same co-rotating structures or from an initial cloud with aligned vectors of angular momentum.

Young stellar objects (YSOs) may not accumulate their mass steadily, as was previously thought, but in a series of violent events manifesting themselves as sharp stellar brightening. These events can be caused by fragmentation due to gravitational instabilities in massive gaseous disks surrounding young stars, followed by migration of dense gaseous clumps onto the star. We report our high angular resolution, coronagraphic near-infrared polarization imaging observations using the High Contrast Instrument for the Subaru Next Generation Adaptive Optics (HiCIAO) of the Subaru 8.2 m Telescope, towards four YSOs which are undergoing luminous accretion outbursts. The obtained infrared images have verified the presence of several hundred AUs scale arms and arcs surrounding these YSOs. In addition, our hydrodynamics simulations and radiative transfer models further demonstrate that these observed structures can indeed be explained by strong gravitational instabilities occurring at the beginning of the disk formation phase. The effect of those tempestuous episodes of disk evolution on star and planet formation remains to be understood.

Motivated by recent observational and numerical studies suggesting that collapsing protostellar cores may be replenished from the local environment, we explore the evolution of protostellar cores submerged in the external counter-rotating environment. These models predict the formation of counter-rotating disks with a deep gap in the gas surface density separating the inner disk (corotating with the star) and the outer counter-rotating disk. The properties of these gaps are compared to those of planet-bearing gaps that form in disks hosting giant planets. We employ numerical hydrodynamics simulations of collapsing cores that are replenished from the local counter-rotating environment, as well as numerical hydrodynamic simulations of isolated disks hosting giant planets, to derive the properties of the gaps that form in both cases. Our numerical simulations demonstrate that counter-rotating disks can form for a wide range of mass and angular momentum available in the local environment. The gap that separates both disks has a depletion factor smaller than 1%, can be located at a distance from ten to over a hundred AU from the star, and can propagate inward with velocity ranging from 1 AU/Myr to >100 AU/Myr. Unlike our previous conclusion, the gap can therefore be a long-lived phenomenon, comparable in some cases to the lifetime of the disk itself. For a proper choice of the planetary mass, the viscous \alpha-parameter and the disk mass, the planet-bearing gaps and the gaps in counter-rotating disks may show a remarkable similarity in the gas density profile and depletion factor, which may complicate their observational differentiation.

We report Submillimeter Array (SMA) 1.3 mm high angular resolution observations towards the four EXor type outbursting young stellar objects (YSOs) VY Tau, V1118 Ori, V1143 Ori, and NY Ori. The data mostly show low dust masses $M_{dust}$ in the associated circumstellar disks. Among the sources, NY Ori possesses a relatively massive disk with $M_{dust} \sim 9 \times 10^{-4}$ $M_{\odot}$. V1118 Ori has a marginal detection equivalent to $M_{dust} \sim 6 \times 10^{-5}$ $M_{\odot}$. V1143 Ori has a non-detection also equivalent to $M_{dust} < 6 \times 10^{-5}$ $M_{\odot}$. For the nearest source VY Tau, we get a surprising non-detection which provides a stringent upper limit $M_{dust} < 6 \times 10^{-6}$ $M_{\odot}$. We interpret our findings as suggesting that the gas and dust reservoirs that feed the short duration, repetitive optical outbursts seen in some EXors may be limited to the small scale, innermost region of their circumstellar disks. This hot dust may have escaped our detection limits. Follow-up, more sensitive millimeter observations are needed to improve our understanding of the triggering mechanisms of EXor type outbursts.

The formation of SMSs is a potential pathway to seed SMBHs in the early universe. A critical issue for forming SMSs is stellar UV feedback, which may limit the stellar mass growth via accretion. In this paper we study the evolution of an accreting SMS and its UV emissivity under conditions of realistic variable accretion from a self-gravitating circumstellar disc. First we conduct a 2D hydrodynamical simulation to follow the long-term protostellar accretion until the stellar mass exceeds $10^4~M_\odot$. The disc fragments due to gravitational instability, creating a number of small clumps that rapidly migrate inward to fall onto the star. The resulting accretion history is thus highly time-dependent: short episodic accretion bursts are followed by longer, relative quiescent phases. We show that the circumstellar disc for the so-called direct collapse model is more unstable and generates greater variability over shorter timescales than normal Pop III cases. We conduct a post-process stellar evolution calculation using the obtained accretion history. Our results show that, regardless of the strong variability of the accretion rates, the stellar radius monotonically increases with almost constant effective temperature at $T_{\rm eff} \simeq 5000$ K as the stellar mass increases. The resulting UV feedback is too weak to hinder mass accretion due to the low flux of stellar UV photons, thus verifying our implicit assumption of no stellar feedback during the hydrodynamic simulations. The insensitivity of stellar evolution to variable accretion is attributed to the fact that typical timescales of variability, $\lesssim 10^3$ years, are too short to affect the stellar structure. We argue that this evolution will continue until the SMS eventually collapses to produce a massive black hole by the general relativistic instability after the stellar mass reaches $\gtrsim 10^5~M_\odot$.

We present continuum and molecular line observations at 230 GHz and 345 GHz from the Sub-millimeter Array (SMA) toward three protostars in the Perseus L1448N region. The data are from the large project "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C18O (J=2-1), 13CO (J=2-1), and HCO+ (J=4-3), and have offsets with N2D+ (J=3-2) peaks. High angular resolution data from a complimentary continuum survey with the Karl G. Jansky Very Large Array show that Source B is associated with three 8 mm continuum objects, Source NW with two, and Source A remains single. These results suggest that multiplicity in L1448N exists at different spatial scales from a few thousand AU to < 100 AU. Velocity gradients in each source obtained from two-dimensional fits to the SMA C18O emission are found to be perpendicular to within 20 degrees of the outflow directions as revealed by 12CO (J=2-1). We have observed that Sources B and NW with multiplicity have higher densities than Source A without multiplicity. This suggests that thermal Jeans fragmentation can be relevant in the fragmentation process. However, we have not observed a difference in the ratio between rotational and gravitational energy between sources with and without multiplicity. We also have not observed a trend between non-thermal velocity dispersions and the level of fragmentation. Our study has provided the first direct and comprehensive comparison between multiplicity and core properties in low-mass protostars, although based on small number statistics.

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

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.

A numerical hydrodynamical model for the evolution of spherically symmetric collapsing clouds, designed for the calculation of the thermal structure of these objects in both the prestellar and protostellar stages of their evolution, is presented. Distinctive features of the model include the possibility of independently describing the temperatures of the gas and dust, which is extremely important when calculating the thermal structure of prestellar and protostellar clouds, and the account of the radiation flux from the central protostar. This model is used to compare the theoretical density and temperature distributions with observations for nearby sites of star formation obtained with the Herschel Space Observatory. Application of the diffusion approximation with a flux limiter describes well the radial density and temperature distributions in protostellar clouds. However, significant differences between the model and observational density profiles were found for prestellar stages, suggesting the presence of appreciable deviations from equilibrium in the prestellar clouds. An approximate method for calculating the thermal structure of a cloud based on the adaptive $\tau$-approximation is presented. Application of the $\tau$-approximation yields good agreement with the diffusion approximation for the prestellar phase, but produces appreciable discrepancies for the protostellar phase, when the thermal structure of the accreting envelope is determined by the radiation of the protostar.

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.

A physical mechanism that drives FU Orionis-type outbursts is reconsidered. We study the effect of inner part of a circumstellar disk covering a region from near the central star to the radius of approximately $5$ AU (hereafter, the inner disk). Using the fluctuated mass accretion rate onto the inner disk $\dot{M}_{\rm out}$, we consider the viscous evolution of the inner disk and the time variability of the mass accretion rate onto the central star $\dot{M}_{\rm in}$ by means of numerical calculation of an unsteady viscous accretion disk in a one-dimensional axisymmetric model. First, we calculate the evolution of the inner disk assuming an oscillating $\dot{M}_{\rm out}$. It is shown that the time variability of $\dot{M}_{\rm in}$ does not coincide with $\dot{M}_{\rm out}$ due to viscous diffusion. Second, we investigate the properties of spontaneous outbursts with temporally constant $\dot{M}_{\rm out}$. Outburst occur only in a limited range of mass accretion rates onto the inner disk $10^{-10}<\dot{M}_{\rm out}< 3\times 10^{-6}~{\rm M}_{\odot} {\rm yr}^{-1}$ due to gravo-magneto limit cycle (GML). Finally, we discuss the case with a combination of episodic $\dot{M}_{\rm out}$ and accretion outbursts cause by the GML in the inner disk. The GML can drive accretion outbursts onto the star even for the case of fluctuating $\dot{M}_{\rm out}$, although fluctuations of $\dot{M}$ decay during transmitting the inner disk inwards. We newly identified two modes of outburst which are spontaneous one and stimulated one. In a stimulated mode of outburst, $\dot{M}_{\rm out}$ does appear directly in $\dot{M}_{\rm in}$ (the latter defining the stellar accretion luminosity). In a spontaneous mode of outburst, $\dot{M}_{\rm out}$ appears as the interval between outbursts.

We couple non-magnetic, hydrodynamical simulations of collapsing protostellar cores with radiative transfer evolutionary models to generate synthetic observations. We then use these synthetic observations to investigate the extent to which a simple method for measuring protostellar disc masses used in the literature recovers the intrinsic masses of the discs formed in the simulations. We evaluate the effects of contamination from the surrounding core, partially resolving out the disc, optical depth, fixed assumed dust temperatures, inclination, and the dust opacity law. We show that the combination of these effects can lead to disc mass underestimates by up to factors of 2-3 at millimeter wavelengths and up to an order of magnitude or larger at submillimeter wavelengths. The optically thin portions of protostellar discs are generally cooler in the Class I stage than the Class 0 stage since Class I discs are typically larger and more optically thick, and thus more shielded. The observed disc mass distribution closely resembles the intrinsic distribution if this effect is taken into account, especially at millimeter wavelengths where optical depth effects are minimized. Approximately 50%-70% of protostellar discs observed to date with this method are consistent with the masses of the gravitationally unstable discs formed in the simulations, suggesting that at least some protostellar discs are likely sufficiently massive to fragment. We emphasize key future work needed to confirm these results, including assembling larger, less biased samples, and using molecular line observations to distinguish between rotationally supported, Keplerian discs and magnetically supported pseudodiscs.

In the last twenty years, the topic of episodic accretion has gained significant interest in the star formation community. It is now viewed as a common, though still poorly understood, phenomenon in low-mass star formation. The FU Orionis objects (FUors) are long-studied examples of this phenomenon. FUors are believed to undergo accretion outbursts during which the accretion rate rapidly increases from typically $10^{-7}$ to a few $10^{-4}$ $M_\odot$ yr$^{-1}$, and remains elevated over several decades or more. EXors, a loosely defined class of pre-main sequence stars, exhibit shorter and repetitive outbursts, associated with lower accretion rates. The relationship between the two classes, and their connection to the standard pre-main sequence evolutionary sequence, is an open question: do they represent two distinct classes, are they triggered by the same physical mechanism, and do they occur in the same evolutionary phases? Over the past couple of decades, many theoretical and numerical models have been developed to explain the origin of FUor and EXor outbursts. In parallel, such accretion bursts have been detected at an increasing rate, and as observing techniques improve each individual outburst is studied in increasing detail. We summarize key observations of pre-main sequence star outbursts, and review the latest thinking on outburst triggering mechanisms, the propagation of outbursts from star/disk to disk/jet systems, the relation between classical EXors and FUors, and newly discovered outbursting sources -- all of which shed new light on episodic accretion. We finally highlight some of the most promising directions for this field in the near- and long-term.

Stars form from the gravitational collapse of dense molecular cloud cores. In the protostellar phase, mass accretes from the core onto a protostar, likely through an accretion disk, and it is during this phase that the initial masses of stars and the initial conditions for planet formation are set. Over the past decade, new observational capabilities provided by the Spitzer Space Telescope and Herschel Space Observatory have enabled wide-field surveys of entire star-forming clouds with unprecedented sensitivity, resolution, and infrared wavelength coverage. We review resulting advances in the field, focusing both on the observations themselves and the constraints they place on theoretical models of star formation and protostellar evolution. We also emphasize open questions and outline new directions needed to further advance the field.

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.

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.

We study the formation and long-term evolution of primordial protostellar disks harbored by first stars using numerical hydrodynamics simulations in the thin-disk limit. The initial conditions are specified by pre-stellar cores with distinct mass, angular momentum, and temperature. This allows us to probe several tens of thousand years of the disk's initial evolution, during which we observe multiple episodes of fragmentation leading to the formation of gravitationally bound gaseous clumps within spiral arms. These fragments are torqued inward due to gravitational interaction with the spiral arms on timescales of 10^3 - 10^4 yr and accreted onto the growing protostar, giving rise to accretion and luminosity bursts. The burst phenomenon is fueled by continuing accretion of material falling onto the disk from the collapsing parent core, which replenishes the mass lost by the disk due to accretion, and triggers repetitive episodes of disk fragmentation. We show that the burst phenomenon is expected to occur for a wide spectrum of initial conditions in primordial pre-stellar cores and speculate on how the intense luminosities (~10^7 solar luminosities) produced by this mechanism may have important consequences for the disk evolution and subsequent growth of the protostar.

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)

A new view of disk evolution is emerging from self-consistent numerical simulation modeling of the formation of circumstellar disks from the direct collapse of prestellar cloud cores. This has implications for many aspects of star and planet formation, including the growth of dust and high-temperature processing of materials. A defining result is that the early evolution of a disk is crucially affected by the continuing mass loading from the core envelope, and is driven into recurrent phases of gravitational instability. Nonlinear spiral arms formed during these episodes fragment to form gaseous clumps in the disk. These clumps generally migrate inward due to gravitational torques arising from their interaction with a trailing spiral arm. Occasionally, a clump can open up a gap in the disk and settle into a stable orbit, revealing a direct pathway to the formation of companion stars, brown dwarfs, or giant planets. At other times, when multiple clumps are present, a low mass clump may even be ejected from the system, providing a pathway to the formation of free-floating brown dwarfs and giant planets in addition to low mass stars. Finally, it has been suggested that the inward migration of gaseous clumps can provide the proper conditions for the transport of high-temperature processed solids from the outer disk to the inner disk, and even possibly accelerate the formation of terrestrial planets in the inner disk. All of these features arising from clump formation and migration can be tied together conceptually in a Migrating Embryo model for disk evolution that can complement the well-known Core Accretion model for planet formation.

The paper considers the evolution of the supernova envelopes produced by Population III stars with masses of $M_*\sim 25-200 M_\odot$ located in non-rotating protogalaxies with masses of $M\sim 10^7 M_\odot$ at redshifts $z=12$, with dark-matter density profiles in the form of modified isothermal spheres. The supernova explosion occurs in the ionization zone formed by a single parent star. The properties of the distribution of heavy elements (metals) produced by the parent star are investigated, as well as the efficiency with which they are mixed with the primordial gas in the supernova envelope. In supernovae with high energies ($E\simgt 5\times 10^{52}$ erg), an appreciable fraction of the gas can be ejected from the protogalaxy, but nearly all the heavy elements remain in the protogalaxy. In explosions with lower energies ($E\simlt 3\times 10^{52}$ erg), essentially no gas and heavy elements are lost from the protogalaxy: during the first one to threemillion years, the gas and heavy elements are actively carried from the central region of the protogalaxy ($r\sim 0.1 r_{vir}$, where $r_{vir}$ is the virial radius of the protogalaxy), but an appreciable fraction of the mass of metals subsequently returns when the hot cavity cools and the envelope collapses. Supernovae with high energies ($E\simgt 5\times 10^{52}$ erg) are characterized by a very low efficiency of mixing of metals; their heavy elements are located in the small volume occupied by the disrupted envelope (in a volume comparable with that of the entire envelope), with most of the metals remaining inside the hot, rarified cavity of the envelope. (abridged)

For simulations that deal only with dark matter or stellar systems, the conventional N-body technique is fast, memory efficient, and relatively simple to implement. However when including the effects of gas physics, mesh codes are at a distinct disadvantage compared to SPH. Whilst implementing the N-body approach into SPH codes is fairly trivial, the particle-mesh technique used in mesh codes to couple collisionless stars and dark matter to the gas on the mesh, has a series of significant scientific and technical limitations. These include spurious entropy generation resulting from discreteness effects, poor load balancing and increased communication overhead which spoil the excellent scaling in massively parallel grid codes. We propose the use of the collisionless Boltzmann moment equations as a means to model collisionless material as a fluid on the mesh, implementing it into the massively parallel FLASH AMR code. This approach, which we term "collisionless stellar hydrodynamics" enables us to do away with the particle-mesh approach. Since the parallelisation scheme is identical to that used for the hydrodynamics, it preserves the excellent scaling of the FLASH code already demonstrated on peta-flop machines. We find the classic hydrodynamic equations and Boltzmann moment equations can be reconciled under specific conditions, allowing us to generate analytic solutions for collisionless systems using conventional test problems. We confirm the validity of our approach using a suite of demanding test problems, including the use of a modified Sod shock test. We conclude by demonstrating the ability of our code to model complex phenomena by simulating the evolution of a spiral galaxy whose properties agree with those predicted by swing amplification theory. (Abridged)

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.

We present simulation results for the formation and long-term evolution of a primordial protostellar disk harbored by a first star. Using a 2+1D nonaxisymmetric thin disk numerical simulation, together with a barotropic relation for the gas, we are able to probe ~20 kyr of the disk's evolution. During this time period we observe fragmentation leading to loosely bound gaseous clumps within the disk. These are then torqued inward and accreted onto the growing protostar, giving rise to a burst phenomenon. The luminous feedback produced by this mechanism may have important consequences for the subsequent growth of the protostar.

We consider the evolution of the ionization zone around Population III stars with $M_*\sim 25-200 M_\odot$ in protogalaxies with $M\sim 10^7 M_\odot$ at redshifts $z = 12$, assuming that the dark matter profile is a modified isothermal sphere. We study the conditions for the growth of instabilities in the ionization zones. The Rayleigh-Taylor and thermal instabilities develop efficiently in the ionization zones around 25-40 $M_\odot$ stars, while this efficiency is lower for stars with $\sim 120 M_\odot$. For more massive stars ($\sim 200 M_\odot$), the flux of ionizing photons is strong enough to considerably reduce the gas density in the ionization zone, and the typical lifetimes of stars ($\sim 2$ Myr) are insufficient for the growth of instabilities. The gas in a protogalaxy with $M\sim 10^7 M_\odot$ with a 200 $M_\odot$ central star is completely ionized by the end of the star's lifetime; in the case of a 120 $M_\odot$ central star, only one-third of the total mass of gas is ionized. Thus, ionizing photons from stars with $M_*\simlt 120 M_\odot$ cannot leave protogalaxies with $M\simgt 10^7 M_\odot$. If the masses of the central stars are 25 and 40 $M_\odot$, the gas in protogalaxies of this mass remains essentially neutral. We discuss the consequences of the evolution of the ionization zones for the propagation of the envelope after the supernova explosions of the stars and the efficiency of enrichment of the intergalactic medium in heavy elements.

The idea that non steady accretion during the embedded phase of protostar evolution can produce the observed luminosity spread in the Herzsprung-Russell diagram (HRD) of young clusters has recently been called into question. Observations of Fu Ori, for instance, suggest an expansion of the star during strong accretion events whereas the luminosity spread implies a contraction of the accreting objects, decreasing their radiating surface. In this paper, we present a global scenario based on calculations coupling episodic accretion histories derived from numerical simulations of collapsing cloud prestellar cores of various masses and subsequent protostar evolution. Our calculations show that, assuming an initial protostar mass $\mi \sim 1\,\mjup$, typical of the second Larson's core, both the luminosity spread in the HRD and the inferred properties of Fu Ori events (mass, radius, accretion rate) can be explained by this scenario, providing two conditions. First, there must be some variation within the fraction of accretion energy absorbed by the protostar during the accretion process. Second the range of this variation should increase with increasing accretion burst intensity, and thus with the initial core mass and final star mass. The numerical hydrodynamics simulations of collapsing cloud prestellar cores indeed show that the intensity of the accretion bursts correlates with the mass and initial angular momentum of the prestellar core. Massive prestellar cores with high initial angular momentum are found to produce intense bursts characteristic of Fu Ori like events. Our results thus suggest a link between the burst intensities and the fraction of accretion energy absorbed by the protostar, with some threshold in the accretion rate, of the order of $10^{-5}\msolyr$, delimitating the transition from "cold" to "hot" accretion. [Abridged]

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.

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.

We determine the observational signatures of protostellar cores by coupling two-dimensional radiative transfer calculations with numerical hydrodynamical simulations that predict accretion rates that both decline with time and feature short-term variability and episodic bursts caused by disk gravitational instability and fragmentation. We calculate the radiative transfer of the collapsing cores throughout the full duration of the collapse, using as inputs the core, disk, and protostellar masses, radii, and mass accretion rates predicted by the hydrodynamical simulations. From the resulting spectral energy distributions, we calculate standard observational signatures (bolometric luminosity, bolometric temperature, ratio of bolometric to submillimeter luminosity) to directly compare to observations. We show that the accretion process predicted by these models reproduces the full spread of observed protostars in both Lbol - Tbol and Lbol - core mass space, including very low luminosity objects, provides a reasonable match to the observed protostellar luminosity distribution, and resolves the long-standing luminosity problem. These models predict an embedded phase duration shorter than recent observationally determined estimates (0.12 Myr vs. 0.44 Myr), and a fraction of total time spent in Stage 0 of 23%, consistent with the range of values determined by observations. On average, the models spend 1.3% of their total time in accretion bursts, during which 5.3% of the final stellar mass accretes, with maximum values being 11.8% and 35.5% for the total time and accreted stellar mass, respectively. Time-averaged models that filter out the accretion variability and bursts do not provide as good of a match to the observed luminosity problem, suggesting that the bursts are required.

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.

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.

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 disk's 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.

We revisit our original papers on the burst mode of accretion by incorporating a detailed energy balance equation into a thin-disk model for the formation and evolution of circumstellar disks around low-mass protostars.Our model includes the effect of radiative cooling, viscous and shock heating, and heating due to stellar and background irradiation. Following the collapse from the prestellar phase allows us to model the early embedded phase of disk formation and evolution. During this time, the disk is susceptible to fragmentation, depending upon the properties of the initial prestellar core. Globally, we find that higher initial core angular momentum and mass content favors more fragmentation, but higher levels of background radiation can moderate the tendency to fragment. A higher rate of mass infall onto the disk than that onto the star is a necessary but not sufficient condition for disk fragmentation. More locally, both the Toomre Q-parameter needs to be below a critical value _and_ the local cooling time needs to be shorter than a few times the local dynamical time. Fragments that form during the early embedded phase tend to be driven into the inner disk regions, and likely trigger mass accretion and luminosity bursts that are similar in magnitude to FU-Orionis-type or EX-Lupi-like events. Disk accretion is shown to be an intrinsically variable process, thanks to disk fragmentation, nonaxisymmetric structure, and the effect of gravitational torques. The additional effect of a generic \alpha-type viscosity acts to reduce burst frequency and accretion variability, and is likely to not be viable for values of \alpha significantly greater than 0.01.

Motivated by the recent discovery of massive planets on wide orbits, we present a mechanism for the formation of such planets via disk fragmentation in the embedded phase of star formation. In this phase, the forming disk intensively accretes matter from the natal cloud core and undergoes several fragmentation episodes. However, most fragments are either destroyed or driven into the innermost regions (and probably onto the star) due to angular momentum exchange with spiral arms, leading to multiple FU-Ori-like bursts and disk expansion. Fragments that are sufficiently massive and form in the late embedded phase (when the disk conditions are less extreme) may open a gap and evolve into giant planets on typical orbits of several tens to several hundreds of AU. For this mechanism to work, the natal cloud core must have sufficient mass and angular momentum to trigger the burst mode and also form extended disks of the order of several hundreds of AU. When mass loading from the natal cloud core diminishes and the main fragmentation phase ends, such extended disks undergo a transient episode of contraction and density increase, during which they may give birth to a last and survivable set of giant planets on wide and relatively stable orbits.

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

Motivated by the recent detection of a large number of embedded young stellar objects (YSOs) with mass accretion rates that are inconsistent with the predictions of the standard model of inside-out collapse (Shu 1977), we perform a series on numerical hydrodynamic simulations of the gravitational collapse of molecular cloud cores with various initial masses, rotation rates, and sizes. We focus on the early Class I stage of stellar evolution when circumstellar disks are exposed to high rates of mass deposition from infalling envelopes. Our numerical modeling reproduces the large observed spread in accretion rates inferred for embedded YSOs in Perseus, Serpens, and Ophiuchus star forming regions by Enoch et al. (2009), yielding 37%--75% of objects with "sub-Shu" accretion rates \dot{M} \la 10^{-6} Msun/yr and 1%--2% of objects with "super-Shu" accretion rates \dot{M}>10^{-5} Msun/yr. Mass accretion rates in the Class I stage have a log-normal distribution, with its shape controlled by disk viscosity and disk temperature. The spread in $\dot{M}$ is greater in models with lower viscosity and smaller in models with higher viscosity and higher disk temperature, suggesting that gravitational instability may be a dominant cause of the observed diversity in $\dot{M}$ in embedded YSOs. Our modeling predicts a weak dependence between the mean mass accretion rates and stellar masses in the Class I stage, in sharp contrast to the corresponding steep dependence for evolved T Tauri stars and brown dwarfs.

We present numerical solutions of the collapse of prestellar cores that lead to the formation and evolution of circumstellar disks. The disk evolution is then followed for up to three million years. A variety of models of different initial masses and rotation rates allows us to study disk accretion around brown dwarfs and low-mass T Tauri stars, with central object mass $M_st < 0.2 Msun$, as well as intermediate and upper-mass T Tauri stars (0.2 Msun < M_st < 3.0 Msun). Our models include self-gravity and allow for nonaxisymmetric motions. In addition to the self-consistently generated gravitational torques, we introduce an effective turbulent \alpha-viscosity with \alpha = 0.01, which allows us particularly to model accretion in the low-mass regime where disk self-gravity is diminishing. A range of models with observationally-motivated values of the initial ratio of rotational to gravitational energy yields a correlation between mass accretion rate \dot{M} and M_st that is relatively steep, as observed. Additionally, our modeling reveals evidence for a bimodality in the \dot{M}--M_st correlation, with a steeper slope at lower masses and a shallower slope at intermediate and upper masses, as also implied by observations. We show that the neglect of disk self-gravity leads to a much steeper \dot{M}--M_st relation for intermediate and upper-mass T Tauri stars. This demonstrates that an accurate treatment of global self-gravity is essential to understanding observations of circumstellar disks.

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

(Abridged) We study numerically the applicability of the effective-viscosity approach for simulating the effect of gravitational instability (GI) in disks of young stellar objects with different disk-to-star mass ratios \xi. We adopt two \alpha$-parameterizations for the effective viscosity based on Lin & Pringle (1990) and Kratter et al (2008) and compare the resultant disk structure, disk and stellar masses, and mass accretion rates with those obtained directly from numerical simulations of self-gravitating disks around low-mass (M_\ast ~ 1.0 M_sun) protostars. We find that the effective viscosity can, in principle, simulate the effect of GI in stellar systems with \xi <= 0.2-0.3, thus corroborating a similar conclusion by Lodato & Rice (2004) that was based on a different \alpha-parameterization. In particular, the Kratter et al's \alpha-parameterization has proven superior to that of Lin & Pringle's, because the success of the latter depends crucially on the proper choice of the \alpha-parameter. However, the \alpha-parameterization generally fails in stellar systems with \xi >= 0.3, particularly in the Class 0 and Class I phases of stellar evolution, yielding too small stellar masses and too large disk-to-star mass ratios. The failure of the \alpha-parameterization in the case of large \xi is caused by a growing strength of low-order spiral modes in massive disks.

We add the effect of turbulent viscosity via the \alpha-prescription to models of the self-consistent formation and evolution of protostellar discs. Our models are non-axisymmetric and carried out using the thin-disc approximation. Self-gravity plays an important role in the early evolution of a disc, and the later evolution is determined by the relative importance of gravitational and viscous torques. In the absence of viscous torques, a protostellar disc evolves into a self-regulated state with disk-averaged Toomre parameter Q \sim 1.5-2.0, non-axisymmetric structure diminishing with time, and maximum disc-to-star mass ratio \xi = 0.14. We estimate an effective viscosity parameter \alpha_eff associated with gravitational torques at the inner boundary of our simulation to be in the range 10^{-4}-10^{-3} during the late evolution. Addition of viscous torques with a low value \alpha = 10^{-4} has little effect on the evolution, structure, and accretion properties of the disc, and the self-regulated state is largely preserved. A sequence of increasing values of \alpha results in the discs becoming more axisymmetric in structure, being more gravitationally stable, having greater accretion rates, larger sizes, shorter lifetimes, and lower disc-to-star mass ratios. For \alpha=10^{-2}, the model is viscous-dominated and the self-regulated state largely disappears by late times. (Abridged)

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_{\phi\phi}:\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.

(Abridged). Motivated by a growing concern that masses of circumstellar disks may have been systematically underestimated by conventional observational methods, we present a numerical hydrodynamics study of time-averaged disk masses (<M_d>) around low-mass Class 0, Class I, and Class II objects. Mean disk masses (\overline{M}_d}) are then calculated by weighting the time-averaged disk masses according to the corresponding stellar masses using a power-law weight function with a slope typical for the Kroupa initial mass function of stars. Two distinct types of disks are considered: self-gravitating disks, in which mass and angular momentum are redistributed exclusively by gravitational torques, and viscous disks, in which both the gravitational and viscous torques are at work. We find that self-gravitating disks have mean masses that are slowly increasing along the sequence of stellar evolution phases. More specifically, Class 0/I/II self-gravitating disks have mean masses \overline{M}_d=0.09, 0.10, and 0.12 M_sun, respectively. Viscous disks have similar mean masses (\overline{M}_d=0.10-0.11 M_sun) in the Class 0/I phases but almost a factor of 2 lower mean mass in the Class II phase (\overline{M}_d,CII=0.06 M_sun). In each evolution phase, time-averaged disk masses show a large scatter around the mean value. Our obtained mean disk masses are larger than those recently derived by Andrews & Williams and Brown et al., regardless of the physical mechanisms of mass transport in the disk.

We present numerical studies of the properties of the stellar velocity distribution in galactic disks which have developed a saturated, two-armed spiral structure. In previous papers we used the Boltzmann moment equations (BME) up to second order for our studies of the velocity structure in self-gravitating stellar disks. A key assumption of our BME approach is the zero-heat flux approximation, i.e. the neglection of third order velocity terms. We tested this assumption by performing test particle simulations for stars in a disk galaxy subject to a rotating spiral perturbation. As a result we corroborated qualitatively the complex velocity structure found in the BME approach. It turned out that an equilibrium configuration in velocity space is only slowly established on a typical timescale of 5 Gyrs or more. Since many dynamical processes in galaxies (like the growth of spirals or bars)act on shorter timescales, pure equilibrium models might not be fully appropriate for a detailed comparison with observations like the local Galactic velocity distribution. Third order velocity moments were typically small and uncorrelated over almost all of the disk with the exception of the 4:1 resonance region (UHR). Near the UHR (normalized) fourth and fifth order velocity moments are still of the same order as the second and third order terms. Thus, at the UHR higher order terms are not negligible.

(Abridged) We explore numerically the chemical, thermal, and dynamical evolution of a shell formed by a high-energy supernova explosion ($10^{53}$ erg) in dwarf protogalaxies with total mass $10^7 M_\odot$ at a redshift $z=12$. We consider two initial configurations for the baryonic matter, one without rotation and the other having the ratio of rotational to gravitational energy $\beta=0.17$. The (non-rotating) dark matter halo is described by a quasi-isothermal sphere. We find that the dynamics of the shell is different in protogalaxies with and without rotation. For instance, the Rayleigh-Taylor instability in the shell develops faster in protogalaxies without rotation. The fraction of a blown-away baryonic mass is approximately twice as large in models with rotation ($\sim 20%$) than in models without rotation. On the other hand, the chemical evolution of gas in protogalaxies with and without rotation is found to be similar. The relative number densities of molecular hydrogen and HD molecules in the cold gas ($T \le 10^3$ K) saturate at typical values of $10^{-3}$ and $10^{-7}$, respectively. The clumps formed in the fragmented shell move with velocities that are at least twice as large as the escape velocity. The mass of the clumps is $\sim 0.1-10 \msun$, which is lower than the Jeans mass. We conclude that the clumps are pressure supported. A supernova explosion with energy $10^{53}$ ergs destructs our model protogalaxy. The clumps formed in the fragmented shell are pressure supported. We conclude that protogalaxies with total mass $\sim 10^{7} M_\odot$ are unlikely to form stars due to high-energy supernova explosions of the first stars.

We present new multi-band imaging data in the optical (BVRI and Halpha) and near infrared bands (JHK) of 15 candidate ring galaxies from the sample of Appleton & Marston (1997). We use these data to obtain color composite images, global magnitudes and colors of both the ring galaxy and its companion(s), and radial profiles of intensity and colors. We find that only nine of the observed galaxies have multi-band morphologies expected for the classical collisional scenario of ring formation, indicating the high degree of contamination of the ring galaxy sample by galaxies without a clear ring morphology. The radial intensity profiles, obtained by masking the off-centered nucleus, peak at the position of the ring, with the profiles in the continuum bands broader than that in the Halpha line. The images as well as the radial intensity and color profiles clearly demonstrate the existence of the pre-collisional stellar disk outside the star-forming ring, which is in general bluer than the disk internal to the ring. The stellar disk seems to have retained its size, with the disk outside the ring having a shorter exponential scale length as compared to the values expected in normal spiral galaxies of comparable masses. The rings in our sample of galaxies are found to be located preferentially at around half-way through the stellar disk. The most likely reason for this preference is bias against detecting rings when they are close to the center (they would be confused with the resonant rings), and at the edge of the disk the gas surface density may be below the critical density required for star formation. Most of the observed characteristics point to relatively recent collisions (<80 Myr ago) according to the N-body simulations of Gerber et al. (1996).

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.

We present a numerical study of the properties of the stellar velocity distribution in stellar discs which have developed a saturated, two-armed spiral structure. We follow the growth of the spiral structure deeply into the non-linear regime by solving the Boltzmann moment equations up to second order. By adopting the thin-disc approximation, we restrict our study of the stellar velocity distribution to the plane of the stellar disc. We find that the outer (convex) edges of stellar spiral arms are characterized by peculiar properties of the stellar velocity ellipsoids, which make them distinct from most other galactic regions. In particular, the ratio \sigma_1:\sigma_2 of the smallest versus largest principal axes of the stellar velocity ellipsoid can become abnormally small (as compared to the rest of the disc) near the outer edges of spiral arms. Moreover, the epicycle approximation fails to reproduce the ratio of the tangential versus radial velocity dispersions in these regions. These peculiar properties of the stellar velocity distribution are caused by large-scale non-circular motions of stars, which in turn are triggered by the non-axisymmetric gravitational field of stellar spiral arms. The magnitude of the vertex deviation appears to correlate globally with the amplitude of the spiral stellar density perturbations. However, locally there is no simple correlation between the vertex deviation and the density perturbations. (Abstract abridged).

We present a numerical model for the evolution of a protostellar disc that has formed self-consistently from the collapse of a molecular cloud core. The global evolution of the disc is followed for several million years after its formation. The capture of a wide range of spatial and temporal scales is made possible by use of the thin-disc approximation. We focus on the role of gravitational torques in transporting mass inward and angular momentum outward during different evolutionary phases of a protostellar disc with disc-to-star mass ratio of order 0.1. In the early phase, when the infall of matter from the surrounding envelope is substantial, mass is transported inward by the gravitational torques from spiral arms that are a manifestation of the envelope-induced gravitational instability in the disc. In the late phase, when the gas reservoir of the envelope is depleted, the distinct spiral structure is replaced by ongoing irregular nonaxisymmetric density perturbations. The amplitude of these density perturbations decreases with time, though this process is moderated by swing amplification aided by the existence of the disc's sharp outer edge. Our global modelling of the protostellar disc reveals that there is typically a residual nonzero gravitational torque from these density perturbations, i.e. their effects do not exactly cancel out in each region. In particular, the net gravitational torque in the inner disc tends to be negative during first several million years of the evolution, while the outer disc has a net positive gravitational torque. Our global model of a self-consistently formed disc shows that it is also self-regulated in the late phase, so that it is near the Toomre stability limit, with a near-uniform Toomre parameter Q\approx 1.5-2.0. (Abstract abridged).

Propagating star formation in a collisional ring galaxy Arp10 is investigated by a complex approach, which includes the broad- and narrow-band photometry, long- slit spectroscopy, and scanning Fabry-Perot spectroscopy. The ionized gas velocity field obtained with best spatial resolution to date indicates a non- isotropic expansion of the outer ring with a maximum velocity 110km/s. Strong vertical and non-circular motions are also seen in the vicinity of the inner ring. Our kinematic data suggest that Arp10 has a small inclination i=22\degr and high total (luminous plus dark matter) mass of about $10^{12} M_{\odot}$ within a 50 kpc radius.The abundance of oxygen 12 + log(O/H) in both star- forming rings is about 8.6.The analysis of spectral indices provides an estimate on the propagation velocities of both rings and metallicity of the pre-collision stellar population.A small "knot" near the nucleus of Arp10,which was previously suspected as a possible candidate for collision, is now unambiguously identified as the "intruder" of at least 1/4 of the total mass of Arp~10.We use a simplified two-dimensional hydrodynamic modeling of galaxy collisions to test a collisional origin of Arp10. We confirm that the sizes of the inner and outer rings, maximum expansion velocity of the outer ring, and radial profile of the gas circular velocity can be reproduced by a near-central collision with the intruder galaxy, which occurred approximately 85Myr ago. We acknowledge that an apparent crescent- shaped distribution of H$\alpha$ emission in the outer ring is caused by a star formation threshold in the gas disk of Arp10.

We present the Boltzmann moment equation approach for the dynamics of stars (BEADS-2D), which is a finite-difference Eulerian numerical code designed for the modelling of anisotropic and non-axisymmetric flat stellar disks. The BEADS-2D code solves the Boltzmann moment equations up to second order in the thin-disk approximation. This allows us to obtain the anisotropy of the velocity ellipsoid and the vertex deviation in the plane of the disk. We study the time-dependent evolution of exponential stellar disks in the linear regime and beyond. The disks are initially characterized by different values of the Toomre parameter Q_s and are embedded in a dark matter halo, yielding a rotation curve composed of a rigid central part and a flat outer region. Starting from a near equilibrium state, several unstable modes develop in the disk. In the early linear phase, the very centre and the large scales are characterized by growing one-armed and bisymmetric positive density perturbations, respectively. This is in agreement with expectations from the swing amplification mechanism of short-wavelength trailing disturbances, propagating through the disk centre. In the late linear phase, the overall appearance is dominated by a two-armed spiral structure localized within the outer Lindblad resonance (OLR). During the non-linear evolutionary phase, radial mass redistribution due to the gravitational torques of spiral arms produces an outflow of mass, which forms a ring at the OLR, and an inflow of mass, which forms a transient central bar. This process of mass redistribution is self-regulatory and it terminates when spiral arms diminish due to a shortage of matter. Finally, a compact central disk and a diffuse ring at the OLR are formed (the abstract is abridged).

We present new numerical simulations in the thin-disk approximation which characterize the burst mode of protostellar accretion. The burst mode begins upon the formation of a centrifugally balanced disk around a newly formed protostar. It is comprised of prolonged quiescent periods of low accretion rate (typically $\la 10^{-7} \Msun$ yr$^{-1}$) which are punctuated by intense bursts of accretion (typically $\ga 10^{-4} \Msun$ yr$^{-1}$, with duration $\la 100$ yr) during which most of the protostellar mass is accumulated. The accretion bursts are associated with the formation of dense protostellar/protoplanetary embryos, which are later driven onto the protostar by the gravitational torques that develop in the disk. Gravitational instability in the disk, driven by continuing infall from the envelope, is shown to be an effective means of transporting angular momentum outward, and mass inward to the protostar. We show that the disk mass always remains significantly less than the central protostar mass throughout this process. The burst phenomenon is robust enough to occur for a variety of initial values of rotation rate, frozen-in (supercritical) magnetic field, and density-temperature relations. Even in cases where the bursts are nearly entirely suppressed, a moderate increase in cloud size or rotation rate can lead to vigorous burst activity. We conclude that most (if not all) protostars undergo a burst mode of evolution during their early accretion history, as inferred empirically from observations of FU Orionis variables.

Motivated by recent observations of plateaus and minima in the radial abundance distributions of heavy elements in the Milky Way and some other spiral galaxies, we propose a dynamical mechanism for the formation of such features around corotation. Our numerical simulations show that the non-axisymmetric gravitational field of spiral density waves generates cyclone and anticylone gas flows in the vicinity of corotation. The anticyclones flatten the pre-existing negative abundance gradients by exporting many more atoms of heavy elements outside corotation than importing inside it. This process is very efficient and forms plateaus of several kiloparsec in size around corotation after two revolution periods of a galaxy. The strength of anticyclones and, consequently, the sizes of plateaus depend on the pitch angle of spiral arms and are expected to increase along the Hubble sequence.

We use long-slit spectra obtained with the 6-m telescope BTA (Russian Academy of Sciences) to investigate the history of star formation in the peculiar ring galaxy Arp 10. The radial distributions of quasi-Lick spectral indices are calculated from the observed spectra and are compared with the model spectra, the latter are generated using the Starburst99 population synthesis code. Our model includes an outward propagating density wave which triggers star formation in the gas disk, and an old stellar disk. The metallicity of both old and young stellar populations is a function of radius. We show that a mix of the young and old populations is required to explain the radial distribution of the spectral indices. The density wave propagates outward with a moderate speed (of order 40 km/s), and the metallicities of both young and old populations decrease with the radius. The model indices corresponding to the alpha-elements require somewhat higher metallicities as compared to the Fe-peak elements. We acknowledge partial support from grant RFBR 04-02-16518.

We study numerically the evolution of rotating cloud cores, from the collapse of a magnetically supercritical core to the formation of a protostar and the development of a protostellar disk during the main accretion phase. We find that the disk quickly becomes unstable to the development of a spiral structure similar to that observed recently in AB Aurigae. A continuous infall of matter from the protostellar envelope makes the protostellar disk unstable, leading to spiral arms and the formation of dense protostellar/protoplanetary clumps within them. The growing strength of spiral arms and ensuing redistribution of mass and angular momentum creates a strong centrifugal disbalance in the disk and triggers bursts of mass accretion during which the dense protostellar/protoplanetary clumps fall onto the central protostar. These episodes of clump infall may manifest themselves as episodes of vigorous accretion rate (\ge 10^{-4} M_sun/yr) as is observed in FU Orionis variables. Between these accretion bursts, the protostar is characterized by a low accretion rate (< 10^{-6} M_sun/yr). During the phase of episodic accretion, the mass of the protostellar disk remains less than or comparable to the mass of the protostar.

We investigate the role of non-isothermality in gravitational collapse and protostellar accretion by explicitly including the effects of molecular radiative cooling, gas-dust energy transfer, and cosmic ray heating in models of spherical hydrodynamic collapse. Isothermal models have previously shown an initial decline in the mass accretion rate \dot{M}, due to a gradient of infall speed that develops in the prestellar phase. Our results show that: (1) in the idealized limit of optically thin cooling, a positive temperature gradient is present in the prestellar phase which effectively cancels out the effect of the velocity gradient, producing a near constant \dot{M} in the early accretion phase; (2) in the more realistic case including cooling saturation at higher densities, \dot{M} may initially be either weakly increasing or weakly decreasing with time, for low (T_d ~ 6 K) and high dust temperature (T_d ~ 10 K) cases, respectively. Hence, our results show that the initial decline in \dot{M} seen in isothermal models is definitely not enhanced by non-isothermal effects, and is often suppressed by them. In all our models, \dot{M} does eventually decline rapidly due to the finite mass condition on our cores and a resulting inward propagating rarefaction wave. Thus, any explanation for a rapid decline of $\dot{M}$ in the accretion phase likely needs to appeal to the global molecular cloud structure and possible envelope support, which results in a finite mass reservoir for cores.

Motivated by recent observations which detect dust at large galactocentric distances in the disks of spiral galaxies, we propose a mechanism of outward radial transport of dust by spiral stellar density waves. We consider spiral galaxies in which most of dust formation is localized inside the corotation radius. We show that in the disks of such spiral galaxies, the dust grains can travel over radial distances that exceed the corotation radius by roughly 25%. A fraction of the dust grains can be trapped on kidney-shaped stable orbits between the stellar spiral arms and thus can escape the destructive effect of supernova explosions. These grains form diffuse dusty spiral arms, which stretch 4-5 kpc from the sites of active star formation. About 10% of dust by mass injected inside corotation, can be transported over radial distances 3-4 kpc during approximately 1.0 Gyr. This is roughly an order of magnitude more efficient than can be provided by the turbulent motions.

Motivated by recent observations which detect an outer boundary for starless cores, and evidence for time-dependent mass accretion in the Class 0 and Class I protostellar phases, we reexamine the case of spherical isothermal collapse in the case of a finite mass reservoir. The presence of a core boundary results in the generation of an inward propagating rarefaction wave. This steepens the gas density profile from r^{-2} to r^{-3} or steeper. After a protostar forms, the mass accretion rate \dot{M} evolves through three distinct phases: (1) an early phase of decline in \dot{M}, which is a non-self-similar effect due to spatially nonuniform infall in the prestellar phase; (2) for large cores, an intermediate phase of near-constant \dot{M} from the infall of the outer part of the self-similar density profile; (3) a late phase of rapid decline in \dot{M} when accretion occurs from the region affected by the inward propagating rarefaction wave. Our model clouds of small to intermediate size make a direct transition from phase (1) to phase (3) above. Both the first and second phase are characterized by a temporally increasing bolometric luminosity L_bol, while L_bol is decreasing in the third (final) phase. We identify the period of temporally increasing L_bol with the Class 0 phase, and the later period of terminal accretion and decreasing L_bol with the Class I phase. The peak in L_bol corresponds to the evolutionary time when 50% \pm 10% of the cloud mass has been accreted by the protostar. This is in agreement with the classification scheme proposed by Andre et al. (1993). We show how our results can be used to explain tracks of envelope mass M_env versus L_bol for protostars in Taurus and Ophiuchus. We also develop an analytic formalism which reproduces the protostellar accretion rate.

We present results of Fabry-Perot and long-slit spectroscopy of the peculiar galaxy Arp 10. The ionized gas velocity field shows evidence for significant radial motions in both outer and inner galactic rings. Long-slit spectroscopy reveals gradients of age and metallicity of stellar population in agreement with the propagating nature of star formation in the galaxy. We present strong evidence that a small ``knot'' at 5 arcsec from the center of Arp 10 is its dwarf elliptic satellite, the most probable ``intruder'' responsible for triggering the expanding rings in Arp 10.

New sensitive Very Large Array 20 cm continuum observations of the Cartwheel, the prototypical collisional ring galaxy, were carried out with the principal aim of tracing supernova remnants that are expected to lie in the wake of the expanding ring and in the ring itself. We detect predominantly non-thermal radio continuum emission from regions associated with 13 ring HII complexes. The emission interior to the ring is confined to structures that resemble spokes of the wheel. The spokes start near bright HII complexes, and extend to around 6 arcsec (4 kpc) inward in the direction of the geometrical center of the ring. There is no apparent positional coincidence between the radio continuum and optical spokes. Radial distribution of intensity along the spokes suggests that the past star formation rate (SFR) in the Cartwheel was much lower than the current SFR. New Halpha observations were used to revise the current SFR in the Cartwheel. The revised value is 18 Msun/yr, which is a factor of 4 lower than the value reported previously, but is in good agreement with the SFR estimated from far infrared luminosity. About 30% of the observed 20 cm continuum non-thermal emission seems to originate in processes that are not related to star formation. Revised SFR in the Cartwheel is comparable to that in the rest of the ring galaxies.

We perform numerical hydrodynamic modeling of various physical processes that can form an HI ring as is observed in Holmberg I. Three energetic mechanisms are considered: multiple supernova explosions (SNe), a hypernova explosion associated with a gamma ray burst (GRB), and the vertical impact of a high velocity cloud (HVC). The total released energy has an upper limit of 10^54 ergs. We find that multiple SNe are in general more effective in producing shells that break out of the disk than a hypernova explosion of the same total energy. As a consequence, multiple SNe form rings with a high ring-to-center contrast K<100 in the HI column density, whereas single hypernova explosions form rings with K<10. Only multiple SNe can reproduce both the size (diameter \~1.7 kpc) and the ring-to-center contrast (K ~ 15-20) of the HI ring in Hoolmberg I. High velocity clouds create HI rings that are much smaller in size (< 0.8 kpc) and contrast (K < 4.5) than seen in Holmberg I. We construct model position-velocity (pV) diagrams and find that they can be used to distinguish among different HI ring formation mechanisms. The observed pV-diagrams of Holmberg I are best reproduced by multiple SNe. We conclude that the giant HI ring in Holmberg I is most probably formed by multiple SNe. We also find that the appearance of the SNe-driven shell in the integrated HI image depends on the inclination angle of the galaxy. In nearly face-on galaxies, the integrated HI image shows a ring of roughly constant HI column density surrounding a deep central depression, whereas in considerably inclined galaxies (i > 45 deg) the HI image is characterized by two kidney-shaped density enhancements and a mild central depression.

Multiple SN explosions in disk galaxies efficiently evacuate gas and form cavities with the sizes and shapes of the surrounding envelopes determined by the total amount of injected energy and by the initial gas distribution. Such cavities are seen as HI holes when observed in face-on galaxies. Gas hydrodynamics simulations are performed to obtain the quantitative characteristics of HI holes that could serve for the determination of the gas vertical scale height and the corresponding dark matter content and its distribution. Among these characteristics is the ratio of the maximum column density in the HI ring surrounding the hole to the background HI column density and the vertical expansion velocity of gas in the HI ring. We show that in some cases the extragalactic background ionizing radiation may produce HI holes in the outer regions of galaxies, and can account for the existence of HI holes in nearby face-on galaxies with the apparent lack of an underlying stellar population.

Numerical hydrodynamical modelling of supernova-driven shell formation is performed with a purpose to reproduce a giant HI ring (diameter 1.7 kpc) in the dwarf irregular galaxy Holmberg I (Ho I). We find that the contrast in HI surface density between the central HI depression and the ring is sensitive to the shape of the gravitational potential. This circumstance can be used to constrain the total mass (including the dark matter halo) of nearly face-on dwarf irregulars. We consider two models of Ho I, which differ by an assumed mass of the dark matter halo M_h. The contrast in HI surface density between the central HI depression and the ring, as well as the lack of gas expansion in the central hole, are better reproduced by the model with a massive halo of M_h=6.0*10^9 M_sun than by that with a small halo of M_h=4.0*10^8 M_sun, implying that Ho I is halo-dominated. Assuming the halo mass of 6.0*10^9 M_sum, we determine the mechanical energy required to form the observed ring equal to (3.0 +- 0.5)*10^53 ergs, equivalent 300+-50 Type II supernovae. The inclination of Ho I is constrained to 15-20 degrees by comparing the modelled HI spectrum and channel maps with those observed.

The propagation velocity of the first gas ring in collisional ring galaxies, i.e. the velocity at which the maximum in the radial gas density profile propagates radially in the galactic disk, is usually inferred from the radial expansion velocity of gas in the first ring. Our numerical hydrodynamics modeling of ring galaxy formation however shows that the maximum radial expansion velocity of gas in the first ring ($v_{gas}$) is invariably below the propagation velocity of the first gas ring itself ($v_{ring}$). Modeling of the Cartwheel galaxy indicates that the outer ring is currently propagating at $v_{ring} \approx$ 100 km/s, while the maximum radial expansion velocity of gas in the outer ring is currently $v_{gas} \approx$ 65 km/s. Modeling of the radial B-V/V-K color gradients of the Cartwheel ring galaxy also indicates that the outer ring is propagating at $v_{ring} \ge $ 90 km/s. We show that a combined effect of inclination, finite thickness, and warping of the Cartwheel's disk might be responsible for the lack of angular difference in the peak positions found for the azimuthally averaged $H\alpha$, K and B surface brightness profiles of the Cartwheel's outer ring. Indeed, the radial $H\alpha$ surface brightness profiles obtained along the Cartwheel's major axis, where effects of inclination and finite thickness are minimized, do peak exterior to those at K- and B-bands. The angular difference in peak positions implies $v_{ring}$ = 110 km/s, which is in agreement with the model predictions. We briefly discuss the utility of radio continuum emission and spectral line equivalent widths for determining the propagation velocity of gas rings in collisional ring galaxies.

In this paper we model and analyse the B-V/V-K radial color gradients observed in the Cartwheel ring galaxy. Along with the color-color diagrams, we use the Q_BVK combined color indices, which minimise the uncertainties in the observed B-V and V-K colors introduced by dust extinction. To model the optical and near-infrared color properties of the Cartwheel galaxy, we assume that an intruder galaxy generates an expanding ring density wave in the Cartwheel's disk, which in its turn triggers massive star formation along the wave's perimeter according to the Schmidt law. We use the population synthesis to calculate the color properties of stellar populations formed in the expanding density wave. The results of color modelling suggest that the pre-collision Cartwheel was a late-type spiral, embedded in an extensive gaseous disk of sub-critical surface density. The properties of the old stellar disk are typical for the late-type Freeman disks, with the central surface brightness in V-band and the scale length being 21.0 mag arcsec^-2 and 2 kpc respectively. The pre-collision gaseous disk has a metallicity gradient ranging from z=z_sun/5 at the outer regions to z=z_sun in the central regions. At present, the wave of star formation has passed the initial extent of the pre-collision, old stellar disk and is currently moving in the predominantly gaseous, low-metallicity disk at the radius of 16 kpc. Neither young stellar populations formed in an expanding density wave, nor their mixture with the old, pre-collision stellar populations can reproduce the B-V and V-K colors of the Cartwheel's nucleus+inner ring. We find that an additional 10-Myr-old burst of star formation in the nuclear regions, along with the visual extinction of 1.3 mag, might be responsible for the peculiar colors of the Cartwheel's nucleus.

We compute radial color gradients produced by an outwardly propagating circular wave of star formation and compare our results with color gradients observed in the classical ring galaxy, the ``Cartwheel''. We invoke two independent models of star formation in the ring galaxies. The first one is the conventional density wave scenario, in which an intruder galaxy creates a radially propagating density wave accompanied by an enhanced star formation following the Schmidt law. The second scenario is a pure self-propagating star formation model, in which the intruder only sets off the first burst of stars at the point of impact. Both models give essentially the same results. Systematic reddening of B-V, V-K colors towards the center, such as that observed in the Cartwheel, can be obtained only if the abundance of heavy elements in the star-forming gas is a few times below solar. The B-V and V-K color gradients observed in the Cartwheel can be explained as a result of mixing of stellar populations born in a star-forming wave propagating through a low-metallicity gaseous disk, and a pre-existing stellar disk of the size of the gaseous disk with color properties typical to those observed in nearby disk galaxies.

Ring waves of star formation, propagating outwardly in the galactic disks, leave chemical abundance gradients in their wakes. We show that the relative [Fe/O] abundance gradients in ring galaxies can be used as a tool for determining the role of the SNIa explosions in their chemical enrichment. We consider two mechanisms which can create outwardly propagating star forming rings in a purely gaseous disk -- a self-induced wave and a density wave, and demonstrate that the radial distribution of the relative [Fe/O] abundance gradients does not depend on the particular mechanism of the wave formation or on the parameters of the star-forming process. We show that the [Fe/O] profile is determined by the velocity of the wave, initial mass function, and the initial chemical composition of the star-forming gas. If the role of SNIa explosions is negligible in the chemical enrichment, the ratio [Fe/O] remains constant throughout the galactic disk with a steep gradient at the wave front. If SNIa stars are important in the production of cosmic iron, the [Fe/O] ratio has gradient in the wake of the star-forming wave with the value depending on the frequency of SNIa explosions.

We compute surface brightness profiles of galactic disks for outwardly propagating waves of star formation with a view to investigate the stellar populations in ring galaxies. We consider two mechanisms which can create outwardly propagating star forming rings in a purely gaseous disk --- a self-induced wave and a density wave. We show that the surface brightness profiles produced by both scenarios of ring formation are similar and are strongly sensitive to the velocity of the wave. The results of our computations are compared with the observational quantities sensitive to the young and old stellar populations in the ring galaxies A0035-335 (the Cartwheel galaxy) and VIIZw466. The best fit to the observed radial H_alpha surface brightness distribution in the Cartwheel galaxy is obtained for a wave velocity of about 90 km/s. The red continuum brightness of the ring can be fully explained by the evolving stars present in the trailing part of the wave. However the red continuum brightness in regions internal to the ring indicates that the wave of star formation propagates in a pre-existing stellar disk in the Cartwheel. The H_alpha and K-band surface brightness profiles in VIIZw466 match the values expected from stellar populations produced by a wave of star formation propagating in a purely gaseous disk very well. We conclude that VIIZw466 is probably experiencing the first event of star formation in the disk.