Stars form in clusters, while planets form in gaseous disks around young stars. Cluster dissolution occurs on longer time scales than disk dispersal. Planet formation thus typically takes place while the host star is still inside the cluster. We explore how the presence of other stars affects the evolution of circumstellar disks. Our numerical approach requires multi-scale and multi-physics simulations where the relevant components and their interactions are resolved. The simulations start with the collapse of a turbulent cloud, from which stars with disks form, which are able to influence each other. We focus on the effect of extinction due to residual cloud gas on the early evolution of circumstellar disks. We find that this extinction protects circumstellar disks against external photoevaporation, but these disks then become vulnerable to dynamic truncation by passing stars. We conclude that circumstellar disk evolution is heavily affected by the early evolution of the cluster.

The Atacama Large (Sub-)millimeter Array (ALMA) has provided glimpse of the interstellar medium (ISM) properties of galaxies at the Epoch of Reionization (EoR); however, detailed understanding of their internal structure is still lacking. We present properties of molecular cloud complexes (MCCs) in a prototypical galaxy at this epoch studied in cosmological zoom-in simulations (Leung et al. 2019c). Typical MCC mass and size are comparable to nearby spirals and starburst galaxies (Mgas∼106.5Mȯ and R≃45–100 pc). MCCs are highly supersonic, with velocity dispersion of σgas≃20–100 km s−1 and pressure of P/kB ≃107.6Kcm−3, which are comparable to gas-rich starburst galaxies. In addition, we perform stability analysis to understand the origin and dynamical properties of MCCs. We find that MCCs are globally stable in the main disk of Althæa. Densest regions where star formation is expected to take place in clumps and cores on even smaller scales instead have lower virial parameter and Toomre-Q values. Detailed studies of the star-forming gas dynamics at the EoR thus require a spatial resolution of < 40 pc ( ≃ 0.01″), which is within reach of ALMA, to complement studies of stellar populations at EoR using the James Webb Space Telescope (JWST).

This review covers four current questions in the behavior of the atomic and molecular interstellar medium. These include whether the atomic gas originates primarily in cold streams or hot flows onto galaxies; what the filling factor of cold gas actually is in galactic regions observationally determined to be completely molecular; whether molecular hydrogen determines or merely traces star formation; and whether gravity or turbulence drives the dynamical motions observed in interstellar clouds, with implications on their star formation properties.

The formation of stars from gas drives the evolution of galaxies. Yet, it remains one of the hardest processes to understand when trying to connect observations of modern and high-redshift stellar and galaxy populations to models of large scale structure formation. It has become clear that the star formation rate at redshifts z > 2 drops off rather more quickly than was thought even five years ago. Theoretical models have tended to overpredict the star formation rate at these high redshifts substantially, primarily due to overcooling. Overcooling in galaxies typically occurs because of unphysical radiative cooling. As a result, insufficient turbulence is driven by stellar feedback in galaxies. I show that such turbulence has the net effect of strongly inhibiting star formation, despite its ability to locally promote star formation by compression. Radiation pressure appears less likely to be a dominant driver of the turbulence than has been argued, but supernova and magnetorotational instabilities remain viable mechanisms. Gravity alone cannot be the main driver, as otherwise well-resolved models without feedback would accurately predict star formation rates. Star formation rate surface density correlates well with observed molecular gas surface density, as well as with other tracers of high density material. Correlation does not, however, necessarily imply causation. In this case, it appears that both molecule formation and star formation occur as a consequence of gravitational collapse, with molecules typically playing an important but not an essential role in cooling. The basic concept that gravitational instability drives star formation remains a true guide through the thickets of complexity surrounding this topic. I finally briefly note that understanding ionization heating and radiation pressure from the most massive stars will likely require much higher resolution models (sub-parsec scale) than resolving supernova feedback.

Due to the gravitational influence of density fluctuations driven by magneto-rotational instability in the gas disk, planetesimals and protoplanets undergo diffusive radial migration as well as changes of other orbital properties. The magnitude of the effect on particle orbits has important consequences for planet formation scenarios. We use the local-shearing-box approximation to simulate an ideal, isothermal, magnetized gas disk with vertical density stratification and simultaneously evolve numerous massless particles moving under the gravity of the gas and the host star. Although the results converge with resolution for fixed box dimensions, we find there exists no convergence of the response of the particles to the gravity of the gas against the horizontal box size, up to 16 disk scale heights. This lack of convergence indicate that caution should be exercised when interpreting local-shearing-box models involving gravitational physics of magneto-rotational turbulence.

Star formation in galaxies has been suggested to depend on large-scale gravitational instability or on the pressure required to form molecular hydrogen. I present numerical models and analysis of observations in support of the gravitational instability hypothesis. I also consider whether the correlation between the surface densities of molecular hydrogen and star formation implies causation, and if so in which direction.

Massive stars influence the surrounding universe far out of proportion to their numbers through ionizing radiation, supernova explosions, and heavy element production. Their formation requires the collapse of massive interstellar gas clouds with very high accretion rates. We discuss results from the first three-dimensional simulations of the gravitational collapse of a massive, rotating molecular cloud core that include heating by both non-ionizing and ionizing radiation. Local gravitational instabilities in the accretion flow lead to the build-up of a small cluster of stars. These lower-mass companions subsequently compete with the high-mass star for the same common gas reservoir and limit its overall mass growth. This process is called fragmentation-induced starvation, and explains why massive stars are usually found as members of high-order stellar systems. These simulations also show that the H ii regions forming around massive stars are initially trapped by the infalling gas, but soon begin to fluctuate rapidly. Over time, the same ultracompact H ii region can expand anisotropically, contract again, and take on any of the observed morphological classes. The total lifetime of H ii regions is given by the global accretion timescale, rather than their short internal sound-crossing time. This solves the so-called lifetime problem of ultracompact H ii region. We conclude that the the most significant differences between the formation of low-mass and high-mass stars are all explained as the result of rapid accretion within a dense, gravitationally unstable flow.

The formation of the first stars out of metal-free gas appears to result in stars at least an order of magnitude more massive than in the present-day case. We here consider what controls the transition from a primordial to a modern initial mass function. We study the influence of low levels of metal enrichment and different initial conditions on the cooling and collapse of initially ionized gas in small protogalactic halos using three-dimensional, smoothed particle hydrodynamics simulations. We argue that fragmentation at moderate density depends on the initial conditions for star formation more than on the metal abundances present.

We report on a study of interstellar turbulence driven by both correlated and isolated supernova explosions. We use three-dimensional hydrodynamic models of a vertically stratified interstellar medium run with the adaptive mesh refinement code Flash at a maximum resolution of 2 pc, with a grid size of 0.5 × 0.5 × 10 kpc. Cold dense clouds form even in the absence of self-gravity due to the collective action of thermal instability and supersonic turbulence. Studying these clouds, we show that it can be misleading to predict physical properties such as the star formation rate or the stellar initial mass function using numerical simulations that do not include self-gravity of the gas. Even if all the gas in turbulently Jeans unstable regions in our simulation is assumed to collapse and form stars in local freefall times, the resulting total collapse rate is significantly lower than the value consistent with the input supernova rate. The amount of mass available for collapse depends on scale, suggesting a simple translation from the density PDF to the stellar IMF may be questionable. Even though the supernova-driven turbulence does produce compressed clouds, it also opposes global collapse. The net effect of supernova-driven turbulence is to inhibit star formation globally by decreasing the amount of mass unstable to gravitational collapse.

We present numerical experiments that demonstrate that the nonlinear development of the Toomre instability in disks of isothermal gas, stars, and dark matter reproduces the observed Schmidt Law for star formation. The rate of gas collapse depends exponentially on the (minimum) value of the Toomre parameter in the disk. We demonstrate that spurious fragmentation occurs in the absence of sufficient resolution in our SPH model. Our models also reproduce observed star formation thresholds in disk galaxies. We finally briefly discuss the application of our models to the study of globular cluster formation in merging galaxies.

In cold dark matter cosmological models, the first stars to form are believed to do so within small protogalaxies. We study the influence of low levels of metal enrichment on the cooling and collapse of ionized gas in these protogalactic halos using three-dimensional, smoothed particle hydrodynamics simulations.

We present preliminary results of a high resolution simulation of globular cluster formation in a galaxy merger using GADGET (Springel et al. 2001). A barotropic equation of state (Li et al 2003) is implemented to include effects of cooling and heating. After one orbital period, a dozen proto-globular clusters are identified in the tidal tails.

The Hubble Space Telescope observations showed that some nova outbursts may not produce smooth shells of gas, but thousands of gaseous blobs. In the case of T Pyx the blobs are distributed in eight concentric circles around the exploding star, producing a pattern similar to tree rings. We explore here the formation of the tree rings by hydrodynamical simulations.

The formation of ring nebulae and their shapes is studied by numerical simulations. The old problem of the origin of the nebular shape is addressed. Two proposed mechanisms are discussed: stellar rotation and stellar magnetic fields, or a combination of both. Two and three-dimensional hydrodynamical and magneto-hydrodynamical simulations indicate that these processes are equally efficient in driving elliptical nebulae such as those surrounding Wolf-Rayet stars, like NGC 6888 around WR136. However, bipolar nebulae are better reproduced by rotation.

The physical basis for interpreting observations of nebular morphology around massive stars in terms of the evolution of the central stars is reviewed, and examples are discussed, including NGC 6888, OMC-1, and η Carinae.

We numerically model the impact of superbubbles produced by starbursts on the gas in dwarf galaxies, using realistic gravitational potentials including the contributions from dark matter haloes, and galaxy radii based on empirical laws. We find that the metal-rich gas is very efficiently ejected in many cases, contradicting previous results; the ambient gas is far more difficult to eject. X-ray emitting plumes can last more than 108 yr after the end of the active starburst.

This review attempts to give a coherent explanation of the main observations of the entry Comet Shoemaker-Levy 9 and the aftermath of the resulting explosions by using models of the tidal breakup of the comet, the entry of individual fragments into the jovian atmosphere, and the resulting fireballs and plumes. A critical review shows that the models appear reasonably well understood. The biggest theoretical uncertainties currently concern how to best tie models of the entry to models of the resulting fireballs. The key unknown before the impact was the size and kinetic energy of the comet fragments. The evidence now available includes the behavior of the chain of fragments, the luminosity of the observed visible fireballs and later infrared emission, the chemistry of the spots, and the lack of seismic waves or perturbations at the water cloud pressure level. These observations point to the fragments having diameters under a kilometer, densities of order 0.5 gem“3, and kinetic energies of order 1027 erg.

We follow the interaction of massive stars with their circumstellar media over their entire lifetimes by combining several numerical models. We use a stellar evolution model for 35 and 60 M⊙ stars (Langer et al. in preparation) as inner boundary conditions for one- and two-dimensional gas dynamical models of the circumstellar medium using ZEUS (Stone & Norman 1992). In this poster, we focus on circumstellar evolution during the Luminous Blue Variable (LBV) and Wolf-Rayet (WR) stages.