We investigate how each aspect of a multi-channel stellar feedback model drives the chemodynamical evolution of a low-mass, isolated dwarf galaxy using a suite of high-resolution simulations. Our model follows individual star particles sampled randomly from an adopted initial mass function, considering independently feedback from: supernovae; stellar radiation causing photoelectric heating of dust grains, ionization and associated heating, Lyman-Werner (LW) dissociation of H$_2$, and radiation pressure; and winds from massive main sequence (neglecting their energy input) and asymptotic giant branch (AGB) stars. Radiative transfer is done by ray tracing. We consider the effects each of these processes have on regulating the star formation rate, global properties, multi-phase interstellar medium (ISM), and driving of galactic winds. We follow individual metal species from distinct nucleosynthetic enrichment channels (AGB winds, massive star stellar winds, core collapse and Type Ia supernovae) and pay particular attention to how these feedback processes regulate metal mixing in the ISM, the metal content of outflows, and the stellar abundance patterns in our galaxy. We find that---for a low-metallicity, low-mass dwarf galaxy ---stellar radiation, particularly ionizing radiation and LW radiation, are important sources of stellar feedback whose effects dominate over photoelectric heating and HI radiation pressure. However, feedback is coupled non-linearly, and the inclusion or exclusion of each process produces non-negligible effects. We find strong variations with: the star formation history; the ejection fractions of metals, mass, and energy; and the distribution of elements from different nucleosynthetic sources in both the gas and stars.

The mass, rate, and spins of merging stellar-mass binary black holes (BBHs) detected by aLIGO and Advanced Virgo provide challenges to traditional BBH formation and merger scenarios. An active galactic nucleus (AGN) disk provides a promising additional merger channel, because of the powerful influence of the gas that drives orbital evolution, makes encounters dissipative, and leads to migration. Previous work showed that stellar mass black holes (sBHs) in an AGN disk migrate to regions of the disk, known as migration traps, where positive and negative gas torques cancel out, leading to frequent BBH formation. Here we build on that work by simulating the evolution of additional sBHs that enter the inner disk by either migration or inclination reduction. We also examine whether the BBHs formed in our models have retrograde or prograde orbits around their centers of mass with respect to the disk, determining the orientation, relative to the disk, of the spin of the merged BBHs. We find that orbiters entering the inner disk form BBHs, that are commonly asymmetric in mass, with sBHs on resonant orbits near the migration trap. When these black holes reach 80 solar masses, they form BBHs with sBHs in the migration trap, which over 10 Myr reach about 1000 solar masses. We also find that 68% of the BBHs in our simulation orbit in the retrograde direction, which implies our merger channel could easily produce the low values of the dimensionless aligned spin, chi_eff, observed in a majority of BBH mergers detected thus far.

We study the properties of molecular-forming gas clumps (MGCs) at the epoch of reionization using cosmological zoom-in simulations. We identify MGCs in a z=6 prototypical galaxy ("Althaea") using an H2 density-based clump finder. We compare their mass, size, velocity dispersion, gas surface density, and virial parameter (alpha_vir) to observations. In Althaea, the typical MGC mass and size are Mgas=10^6.5 Msun and R=45-100 pc, which are comparable to those found in nearby spirals and starburst galaxies. MGCs are highly supersonic and supported by turbulence, with rms velocity dispersions of sigma_gas=20-100 km s^-1 and pressure of P/k_B=10^7.6 K cm^-3 (i.e., >1000x with respect to the Milky Way), similar to those found in nearby and z~2 gas-rich starburst galaxies. In addition, we perform stability analysis to understand the origin and dynamical properties of MGCs. We find that MGCs are globally stable in the main disk of Althaea. Densest regions where star formation is expected to take place in clouds and cores on even smaller scales instead have lower alpha_vir and Toomre-Q values. Detailed studies of the star-forming gas dynamics at the epoch of reionization thus require a spatial resolution of <40 pc (=0.01"), which is within reach with the Atacama Large (sub-)Millimeter Array and the Next Generation Very Large Array.

Understanding the formation of stellar clusters requires following the interplay between gas and newly formed stars accurately. We therefore couple the magnetohydrodynamics code FLASH to the N-body code ph4 and the stellar evolution code SeBa using the Astrophysical Multipurpose Software Environment (AMUSE) to model stellar dynamics, evolution, and collisional N-body dynamics and the formation of binary and higher-order multiple systems, while implementing stellar feedback in the form of radiation, stellar winds and supernovae in FLASH. We use this novel numerical method to simulate the formation and early evolution of open clusters of $\sim 1000$ stars formed from clouds with a mass range of $10^3$-$10^5$ M$_\odot$. Analyzing the effects of stellar feedback on the gas and stars of the natal cluster, we find that our clusters are resilient to disruption, even in the presence of intense feedback. This can even slightly increase the amount of dense, Jeans unstable gas by sweeping up shells; thus, a stellar wind strong enough to trap its own H II region shows modest triggering of star formation. Our clusters are born moderately mass segregated, an effect enhanced by feedback, and retained after the ejection of their natal gas, in agreement with observations.

One-zone models constructed to match observed stellar abundance patterns have been used extensively to constrain the sites of nucleosynthesis with sophisticated libraries of stellar evolution and stellar yields. The metal mixing included in these models is usually highly simplified, although it is likely to be a significant driver of abundance evolution. In this work we use high-resolution hydrodynamics simulations to investigate how metals from individual enrichment events with varying source energies $E_{\rm ej}$ mix throughout the multi-phase interstellar medium (ISM) of a low-mass ($M_{\rm gas}=2\times 10^{6}$~M$_{\odot}$), low-metallicity, isolated dwarf galaxy. These events correspond to the characteristic energies of both common and exotic astrophysical sites of nucleosynthesis, including: asymptotic giant branch winds ($E_{\rm ej}\sim$10$^{46}$~erg), neutron star-neutron star mergers ($E_{\rm ej}\sim$10$^{49}$~erg), supernovae ($E_{\rm ej}\sim$10$^{51}$~erg), and hypernovae ($E_{\rm ej}\sim$10$^{52}$~erg). We find the mixing timescales for individual enrichment sources in our dwarf galaxy to be long (100~Myr--1~Gyr), with a clear trend of increasing homogeneity for the more energetic events. Given these timescales, we conclude that the spatial distribution and frequency of events are important drivers of abundance homogeneity on large scales; rare, low $E_{\rm ej}$ events should be characterized by particularly broad abundance distributions. The source energy $E_{\rm ej}$ also correlates with the fraction of metals ejected in galactic winds, ranging anywhere from 60\% at the lowest energy to 95\% for hypernovae. We conclude by examining how the radial position, local ISM density, and global star formation rate influence these results.

The merger rate of stellar-mass black hole binaries (sBHBs) inferred by the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) suggests the need for an efficient source of sBHB formation. Active galactic nucleus (AGN) disks are a promising location for the formation of these sBHBs, as well as binaries of other compact objects, because of powerful torques exerted by the gas disk. These gas torques cause orbiting compact objects to migrate towards regions in the disk where inward and outward torques cancel, known as migration traps. We simulate the migration of stellar mass black holes in an example of a model AGN disk, using an augmented N-body code that includes analytic approximations to migration torques, stochastic gravitational forces exerted by turbulent density fluctuations in the disk, and inclination and eccentricity dampening produced by passages through the gas disk, in addition to the standard gravitational forces between objects. We find that sBHBs form rapidly in our model disk as stellar-mass black holes migrate towards the migration trap. These sBHBs are likely to subsequently merge on short time-scales. The process continues, leading to the build-up of a population of over-massive stellar-mass black holes. The formation of sBHBs in AGN disks could contribute significantly to the sBHB merger rate inferred by LIGO.

We present here the first of a series of papers aimed at better understanding the evolution and properties of giant molecular clouds (GMCs) in a galactic context. We perform high resolution, three-dimensional {\sc arepo} simulations of an interacting galaxy inspired by the well-observed M51 galaxy. Our fiducial simulations include a non-equilibrium, time-dependent, chemical network that follows the evolution of atomic and molecular hydrogen as well as carbon and oxygen self-consistently. Our calculations also treat gas self-gravity and subsequent star formation (described by sink particles), and coupled supernova feedback. In the densest parts of the simulated interstellar medium (ISM) we reach sub-parsec resolution, granting us the ability to resolve individual GMCs and their formation and destruction self-consistently throughout the galaxy. In this initial work we focus on the general properties of the ISM with a particular focus on the cold star-forming gas. We discuss the role of the interaction with the companion galaxy in generating cold molecular gas and controlling stellar birth. We find that while the interaction drives large-scale gas flows and induces spiral arms in the galaxy, it is of secondary importance in determining gas fractions in the different ISM phases and the overall star-formation rate. The behaviour of the gas on small GMC scales instead is mostly controlled by the self-regulating property of the ISM driven by coupled feedback.

Orbital properties of stars, computed from their six-dimensional phase space measurements and an assumed Galactic potential, are used to understand the structure and evolution of the Galaxy. Stellar actions, computed from orbits, have the attractive quality of being invariant under certain assumptions and are therefore used as quantitative labels of a star's orbit. We report a subtle but important systematic error that is induced in the actions as a consequence of local midplane variations expected for the Milky Way. This error is difficult to model because it is non-Gaussian and bimodal, with neither mode peaking on the null value. An offset in the vertical position of the Galactic midplane of $\sim15\,\text{pc}$ for a thin disk-like orbit or $\sim 120\,\text{pc}$ for a thick disk-like orbit induces a $25\%$ systematic error in the vertical action $J_z$. In FIRE simulations of Milky Way-mass galaxies, these variations are on the order of $\sim100\,\text{pc}$ at the solar circle. From observations of the mean vertical velocity variation of $\sim5\text{--}10\,\text{km}\,\text{s}^{-1}$ with radius, we estimate that the Milky Way midplane variations are $\sim60\text{--}170\,\text{pc}$, consistent with three-dimensional dust maps. Action calculations and orbit integrations, which assume the global and local midplanes are identical, are likely to include this induced error, depending on the volume considered. Variation in the local standard of rest or distance to the Galactic center causes similar issues. The variation of the midplane must be taken into account when performing dynamical analysis across the large regions of the disk accessible to Gaia and future missions.

Context. Supersonic disordered flows accompany the formation and evolution of MCs. It has been argued that turbulence can support against gravitational collapse and form hierarchical sub-structures. Aims. We study the time evolution of simulated MCs to investigate: What physical process dominates the driving of turbulent flows? How can these flows be characterised? Do the simulated flows agree with observations? Methods. We analyse three MCs that have formed self-consistently within kpc-scale numerical simulations of the ISM. The simulated ISM evolves under the influence of physical processes including self-gravity, stratification, magnetic fields, supernova-driven turbulence, and radiative heating and cooling. We characterise the flows using VSFs with/out density weighting or a density cutoff, and computed in one or three dimensions. Results. In regions with sufficient resolution, the density-weighted VSFs initially appear to follow the expectations for uniform turbulence consistent with Larson's size-velocity relationship. SN blast wave impacts produce short-lived coherent motions at large scales, increasing the scaling exponents for a crossing time. Gravitational contraction drives small-scale motions, producing scaling coefficients that drop or even turn negative as small scales become dominant. Conclusions. We conclude that two different effects coincidentally reproduce Larson's size velocity relationship. Initially, uniform turbulence dominates, so the energy cascade produces VSFs consistent with Larson's relationship. Later, contraction dominates, the density-weighted VSFs become much shallower or even inverted, but the relationship of the global average velocity dispersion of the MCs to their radius follows Larson's relationship, reflecting virial equilibrium or free-fall collapse. The injection of energy by shocks is visible in the VSFs, but decays within a crossing time.

High Mach number shocks are ubiquitous in interstellar turbulence. The Pencil Code is particularly well suited to the study of magnetohydrodynamics in weakly compressible turbulence and the numerical investigation of dynamos because of its high-order advection and time evolution algorithms. However, the high-order algorithms and lack of Riemann solver to follow shocks make it less well suited to handling high Mach number shocks, such as those produced by supernovae (SNe). Here, we outline methods required to enable the code to efficiently and accurately model SNe, using parameters that allow stable simulation of SN-driven turbulence, in order to construct a physically realistic galactic dynamo model. These include the resolution of shocks with artificial viscosity, thermal conductivity, and mass diffusion; the correction of the mass diffusion terms; and a novel generalization of the Courant condition to include all source terms in the momentum and energy equations. We test our methods with the numerical solution of the one-dimensional (1D) Riemann shock tube (Sod, J. Comput. Phys. 1978, 27), also extended to a 1D adiabatic shock with parameters and Mach number relevant to SN shock evolution, including shocks with radiative losses. We extend our test with the three-dimensional (3D) numerical simulation of individual SN remnant evolution for a range of ambient gas densities typical of the interstellar medium and compare these to the analytical solutions of Sedov-Taylor (adiabatic) and the snowplough and Cioffi, McKee and Bertschinger (Astrophys. J. 1988, 334) results incorporating cooling and heating processes. We show that our new timestep algorithm leads to linear rather than quadratic resolution dependence as the strength of the artificial viscosity varies, because of the corresponding change in the strength of interzone gradients.

The stellar mass binary black hole (sBBH) mergers presently detected by LIGO may originate wholly or in part from binary black hole mergers embedded in disks of gas around supermassive black holes. Determining the contribution of these active galactic nucleus (AGN) disks to the sBBH merger rate enables us to uniquely measure important parameters of AGN disks, including their typical density, aspect ratio, and lifetime, thereby putting unique limits on an important element of galaxy formation. For the first time, gravitational waves will allow us to reveal the properties of the hidden interior of AGN disks, while electromagnetic radiation (EM) probes the disk photosphere. The localization of sBBH merger events from LIGO is generally insufficient for association with a single EM counterpart. However, the contribution to the LIGO event rate from rare source types (such as AGNs) can be determined on a statistical basis. To determine the contribution to the sBBH rate from AGNs in the next decade requires: {\em 1) a complete galaxy catalog for the LIGO search volume, 2) strategic multi-wavelength EM follow-up of LIGO events and 3) significant advances in theoretical understanding of AGN disks and the behavior of objects embedded within them.}

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

We describe a star cluster formation model that includes individual star formation from self-gravitating, magnetized gas, coupled to collisional stellar dynamics. The model uses the Astrophysical Multi-purpose Software Environment (AMUSE) to integrate an adaptive-mesh magnetohydrodynamics code (FLASH) with a fourth order Hermite N-body code (ph4), a stellar evolution code (SeBa), and a method for resolving binary evolution (multiples). This combination yields unique star formation simulations that allow us to study binaries formed dynamically from interactions with both other stars and dense, magnetized gas subject to stellar feedback during the birth and early evolution of stellar clusters. We find that for massive stars, our simulations are consistent with the observed dynamical binary fractions and mass ratios. However, our binary fraction drops well below observed values for lower mass stars, presumably due to unincluded binary formation during initial star formation. Further, we observe a build up of binaries near the hard-soft boundary that may be an important mechanism driving early cluster contraction.

In order to understand the origin of observed molecular cloud properties, it is critical to understand how clouds interact with their environments during their formation, growth, and collapse. It has been suggested that accretion-driven turbulence can maintain clouds in a highly turbulent state, preventing runaway collapse, and explaining the observed non-thermal velocity dispersions. We present 3D, AMR, MHD simulations of a kiloparsec-scale, stratified, supernova-driven, self-gravitating, interstellar medium, including diffuse heating and radiative cooling. These simulations model the formation and evolution of a molecular cloud population in the turbulent interstellar medium. We use zoom-in techniques to focus on the dynamics of the mass accretion and its history for individual molecular clouds. We find that mass accretion onto molecular clouds proceeds as a combination of turbulent and near free-fall accretion of a gravitationally bound envelope. Nearby supernova explosions have a dual role, compressing the envelope, boosting accreted mass, but also disrupting parts of the envelope and eroding mass from the cloud's surface. It appears that the inflow rate of kinetic energy onto clouds from supernova explosions is insufficient to explain the net rate of charge of the cloud kinetic energy. In the absence of self-consistent star formation, conversion of gravitational potential into kinetic energy during contraction seems to be the main driver of non-thermal motions within clouds. We conclude that although clouds interact strongly with their environments, bound clouds are always in a state of gravitational contraction, close to runaway, and their properties are a natural result of this collapse.

Using a high resolution simulation of an isolated dwarf galaxy, accounting for multi-channel stellar feedback and chemical evolution on a star-by-star basis, we investigate how each of 15 metal species are distributed within our multi-phase interstellar medium (ISM) and ejected from our galaxy by galactic winds. For the first time, we demonstrate that the mass fraction probability distribution functions (PDFs) of individual metal species in the ISM are well described by a piecewise log-normal and power-law distribution. The PDF properties vary within each ISM phase. Hot gas is dominated by recent enrichment, with a significant power-law tail to high metal fractions, while cold gas is predominately log-normal. In addition, elements dominated by asymptotic giant branch (AGB) wind enrichment (e.g. N and Ba) mix less efficiently than elements dominated by supernova enrichment (e.g. $\alpha$ elements and Fe). This result is driven by the differences in source energetics and source locations, particularly the higher chance compared to massive stars for AGB stars to eject material into cold gas. Nearly all of the produced metals are ejected from the galaxy (only 4% are retained), but over 20% of metals dominated by AGB enrichment are retained. In dwarf galaxies, therefore, elements synthesized predominately through AGB winds should be both overabundant and have a larger spread compared to elements synthesized in either core collapse or Type Ia supernovae. We discuss the observational implications of these results, their potential use in developing improved models of galactic chemical evolution, and their generalization to more massive galaxies.

The streaming instability is a promising mechanism to drive the formation of planetesimals in protoplanetary disks. To trigger this process, it has been argued that sedimentation of solids onto the mid-plane needs to be efficient and therefore that a quiescent gaseous environment is required. It is often suggested that dead-zone or disk-wind structure created by non-ideal magnetohydrodynamical (MHD) effects meets this requirement. However, simulations have shown that the midplane of a dead zone is not completely quiescent. In order to examine the concentration of solids in such an environment, we use the local-shearing-box approximation to simulate a particle-gas system with an Ohmic dead zone including mutual drag force between the gas and the solids. We systematically compare the evolution of the system with ideal or non-ideal MHD, with or without back-reaction drag force from particles on gas, and with varying solid abundances. Similar to previous investigations of dead zone dynamics, we find that particles of dimensionless stopping time $\tau_s=0.1$ do not sediment appreciably more than those in ideal magneto-rotational turbulence, resulting in a vertical scale height an order of magnitude larger than in a laminar disk. Contrary to the expectation that this should curb the formation of planetesimals, we nevertheless find that strong clumping of solids still occurs in the dead zone when solid abundances are similar to the critical value for a laminar environment. This can be explained by the weak radial diffusion of particles near the mid-plane. The results imply that the sedimentation of particles to the mid-plane is not a necessary criterion for the formation of planetesimals by the streaming instability.

In order to better understand the relationship between feedback and galactic chemical evolution, we have developed a new model for stellar feedback at grid resolutions of only a few parsecs in global disk simulations, using the adaptive mesh refinement hydrodynamics code Enzo. For the first time in galaxy scale simulations, we simulate detailed stellar feedback from individual stars including asymptotic giant branch winds, photoelectric heating, Lyman-Werner radiation, ionizing radiation tracked through an adaptive ray-tracing radiative transfer method, and core collapse and Type Ia supernovae. We furthermore follow the star-by-star chemical yields using tracer fields for 15 metal species: C, N, O, Na, Mg, Si, S, Ca, Mn, Fe, Ni, As, Sr, Y, and Ba. We include the yields ejected in massive stellar winds, but greatly reduce the winds' velocities due to computational constraints. We describe these methods in detail in this work and present the first results from 500~Myr of evolution of an isolated dwarf galaxy with properties similar to a Local Group, low-mass dwarf galaxy. We demonstrate that our physics and feedback model is capable of producing a dwarf galaxy whose evolution is consistent with observations in both the Kennicutt-Schmidt relationship and extended Schmidt relationship. Effective feedback drives outflows with a greater metallicity than the ISM, leading to low metal retention fractions consistent with observations. Finally, we demonstrate that these simulations yield valuable information on the variation in mixing behavior of individual metal species within the multi-phase interstellar medium.

Effective stellar feedback is used in models of galaxy formation to drive realistic galaxy evolution. Models typically include energy injection from supernovae as the dominant form of stellar feedback, often in some form of sub-grid recipe. However, it has been recently suggested that pre-SN feedback (stellar winds or radiation) is necessary in high-resolution simulations of galaxy evolution to properly regulate star formation and properties of the interstellar medium (ISM). Following these processes is computationally challenging, so many prescriptions model this feedback approximately, accounting for the local destruction of dense gas clouds around newly formed stars in lieu of a full radiative transfer calculation. In this work we examine high resolution simulations (1.8~pc) of an isolated dwarf galaxy with detailed stellar feedback tracked on a star-by-star basis. By following stellar ionizing radiation with an adaptive ray-tracing radiative transfer method, we test its importance in regulating star formation and driving outflows in this galaxy. We find that including ionizing radiation reduces the star formation rate (SFR) by over a factor of 5, and is necessary to produce the ISM conditions needed for supernovae to drive significant outflows. We find that a localized approximation for radiation feedback is sufficient to regulate the SFR on short timescales, but does not allow significant outflows. Short and long range radiation effects are both important in driving the evolution of our low-metallicity, low-mass dwarf galaxy. Generalizing these results to more massive galaxies would be a valuable avenue of future research.

A number of ultracompact H II regions in Galactic star forming environments have been observed to vary significantly in radio flux density on timescales of 10-20 years. Theory predicted that such variations should occur when the accretion flow that feeds a young massive star becomes unstable and clumpy. We have targeted the massive star-forming region W49A with the Karl G. Jansky Very Large Array (VLA) for observations at 3.6 cm with the B-configuration at 0.8'' resolution, to compare to nearly identical observations taken almost 21 years earlier (February 2015 and August 1994). Most of the sources in the crowded field of ultracompact and hypercompact H II regions exhibit no significant changes over this time period. However, one source, W49A/G2, decreased by 20% in peak intensity (from 71+/-4 mJy/beam to 57+/-3 mJy/beam), and 40% in integrated flux (from 0.109+/-0.011 Jy to 0.067+/-0.007 Jy), where we cite 5 sigma errors in peak intensity, and 10% errors in integrated flux. We present the radio images of the W49A region at the two epochs, the difference image that indicates the location of the flux density decrease, and discuss explanations for the flux density decrease near the position of W49A/G2.

We investigate the impact of the far ultraviolet (FUV) heating rate on the stability of the three-phase interstellar medium using three-dimensional simulations of a $1$ kpc$^2$, vertically-extended domain. The FUV heating rate sets the range of thermal pressures across which the cold ($\sim10^2$ K) and warm ($\sim10^4$ K) neutral media (CNM and WNM) can coexist in equilibrium. Even absent a variable star formation rate regulating the FUV heating rate, the gas physics keeps the pressure in the two-phase regime: because radiative heating and cooling processes happen on shorter timescales than sound wave propagation, turbulent compressions tend to keep the interstellar medium within the CNM-WNM pressure regime over a wide range of heating rates. The thermal pressure is set primarily by the heating rate with little influence from the hydrostatics. The vertical velocity dispersion adjusts as needed to provide hydrostatic support given the thermal pressure: when the turbulent pressure $\langle\rho\rangle\sigma_z^2$ is calculated over scales $\gtrsim500$ pc, the thermal plus turbulent pressure approximately equals the weight of the gas. The warm gas volume filling fraction is $0.2<f_w<0.8$ over a factor of less than three in heating rate, with $f_w$ near unity at higher heating rates and near zero at lower heating rates. We suggest that cosmological simulations that do not resolve the CNM should maintain an interstellar thermal pressure within the two-phase regime.

In models of galaxy formation, feedback driven both by supernova (SN) and active galactic nucleus (AGN) is not efficient enough to quench star formation in massive galaxies. Models of smaller galaxies have suggested that cosmic rays (CRs) play a major role in expelling material from the star forming regions by diffusing SN energy to the lower density outskirts. We therefore run gas dynamical simulations of galactic outflows from a galaxy contained in a halo with $5\times10^{12}$~M$_{\odot}$ that resembles a local ultraluminous galaxy, including both SN thermal energy and a treatment of CRs using the same diffusion approximation as Salem & Bryan. We find that CR pressure drives a low-density bubble beyond the edge of the shell swept up by thermal pressure, but the main bubble driven by SN thermal pressure overtakes it later, which creates a large-scale biconical outflow. CRs diffusing into the disk are unable to entrain its gas in the outflows, yielding a mass-loading rate of only $\sim0.1~\%$ with varied CR diffusion coefficients. We find no significant difference in mass-loading rates in SN driven outflows with or without CR pressure. Our simulations strongly suggest that it is hard to drive a heavily mass-loaded outflow with CRs from a massive halo potential, although more distributed star formation could lead to a different result.

Meteoritical and astrophysical models of planet formation make contradictory predictions for dust concentration factors in chondrule forming regions of the solar nebula. Meteoritical and cosmochemical models strongly suggest that chondrules, a key component of the meteoritical record, formed in regions with solids-to-gas mass ratios orders of magnitude above background. However, models of dust grain dynamics in protoplanetary disks struggle to surpass factors of a few outside of very brief windows in the lifetime of the dust grains. Worse, those models do not predict significant concentration factors for dust grains the size of chondrule precursors. We briefly develop the difficulty in concentrating dust particles in the context of nebular chondrule formation and show that the disagreement is sufficiently stark that cosmochemists should explore ideas that might revise the concentration factor requirements downwards.

${\bf Aim:}$ To test the hypothesis that radiation pressure from star clusters acting on dust is the dominant feedback agent disrupting the largest star-forming molecular clouds and thus regulating the star-formation process. ${\bf Methods:}$ We perform multi-frequency, 3D, RT calculations including scattering, absorption, and re-emission to longer wavelengths for clouds with masses of $10^4$-$10^7\,$M$_{\odot}$, with embedded clusters and a star formation efficiencies of 0.009%-91%, and varying maximum grain sizes up to 200$\,\mu$m. We calculate the ratio between radiative force and gravity to determine whether radiation pressure can disrupt clouds. ${\bf Results:}$ We find that radiation acting on dust almost never disrupts star-forming clouds. UV and optical photons to which the cloud is optically thick do not scatter much. Instead, they quickly get absorbed and re-emitted by at thermal wavelengths. As the cloud is typically optically thin to far-IR radiation, it promptly escapes, depositing little momentum. The resulting spectrum is more narrowly peaked than the corresponding Planck function with an extended tail at longer wavelengths. As the opacity drops significantly across the sub-mm and mm, the resulting radiative force is even smaller than for the corresponding single-temperature black body. The force from radiation pressure falls below the strength of gravitational attraction by an order of magnitude or more for either Milky Way or starbust conditions. For unrealistically large maximum grain sizes, and star formation efficiencies far exceeding 50% do we find that the strength of radiation pressure can exceed gravity. ${\bf Conclusions:}$ We conclude that radiation pressure acting on dust does not disrupt star-forming molecular clouds in any Local Group galaxies. Radiation pressure thus appears unlikely to regulate the star-formation process on either local or global scales.

We assess the contribution of dynamical hardening by direct three-body scattering interactions to the rate of stellar-mass black hole binary (BHB) mergers in galactic nuclei. We derive an analytic model for the single-binary encounter rate in a nucleus with spherical and disk components hosting a super-massive black hole (SMBH). We determine the total number of encounters $N_{\rm GW}$ needed to harden a BHB to the point that inspiral due to gravitational wave emission occurs before the next three-body scattering event. This is done independently for both the spherical and disk components. Using a Monte Carlo approach, we refine our calculations for $N_{\rm GW}$ to include gravitational wave emission between scattering events. For astrophysically plausible models we find that typically $N_{\rm GW} \lesssim$ 10. We find two separate regimes for the efficient dynamical hardening of BHBs: (1) spherical star clusters with high central densities, low velocity dispersions and no significant Keplerian component; and (2) migration traps in disks around SMBHs lacking any significant spherical stellar component in the vicinity of the migration trap, which is expected due to effective orbital inclination reduction of any spherical population by the disk. We also find a weak correlation between the ratio of the second-order velocity moment to velocity dispersion in galactic nuclei and the rate of BHB mergers, where this ratio is a proxy for the ratio between the rotation- and dispersion-supported components. Because disks enforce planar interactions that are efficient in hardening BHBs, particularly in migration traps, they have high merger rates that can contribute significantly to the rate of BHB mergers detected by the advanced Laser Interferometer Gravitational-Wave Observatory.

The fragmentation of filaments in molecular clouds has attracted a lot of attention as there seems to be a relation between the evolution of filaments and star formation. The study of the fragmentation process has been motivated by simple analytical models. However, only a few comprehensive studies have analysed the evolution of filaments using numerical simulations where the filaments form self-consistently as part of molecular clouds. We address the early evolution of pc-scale filaments that form within individual clouds. We focus on three questions: How do the line masses of filaments evolve? How and when do the filaments fragment? How does the fragmentation relate to the line masses of the filaments? We examine three simulated molecular clouds formed in kpc-scale numerical simulations performed with the FLASH code. We compare the properties of the identified filaments with the predictions of analytic filament stability models. The line masses and mass fraction enclosed in the identified filaments increase continuously after the onset of self-gravity. The first fragments appear early when the line masses lie well below the critical line mass of Ostriker's hydrostatic equilibrium solution. The average line masses of filaments identified in 3D density cubes increases far more quickly than those identified in 2D column density maps. Our results suggest that hydrostatic or dynamic compression from the surrounding cloud has a significant impact on the early dynamical evolution of filaments. A simple model of an isolated, isothermal cylinder may not provide a good approach for fragmentation analysis. Caution must be exercised in interpreting distributions of properties of filaments identified in column density maps, especially in the case of low-mass filaments. Comparing or combining results from studies that use different filament finding techniques is strongly discouraged.

Recent observations of the protoplanetary disk around the Herbig Be star HD 100546 show two bright features in infrared (H and L' bands) at about 50 AU, with one so far unexplained. We explore the observational signatures of a high mass planet causing shock heating in order to determine if it could be the source of the unexplained infrared feature in HD 100546. More fundamentally, we identify and characterize planetary shocks as an extra, hitherto ignored, source of luminosity in transition disks. The RADMC-3D code is used to perform dust radiative transfer calculations on the hydrodynamical disk models, including volumetric heating. A stronger shock heating rate by a factor 20 would be necessary to qualitatively reproduce the morphology of the second infrared source. Instead, we find that the outer edge of the gap carved by the planet heats up by about 50% relative to the initial reference temperature, which leads to an increase in the scale height. The bulge is illuminated by the central star, producing a lopsided feature in scattered light, as the outer gap edge shows an asymmetry in density and temperature attributable to a secondary spiral arm launched not from the Lindblad resonances but from the 2:1 resonance. We conclude that high-mass planets lead to shocks in disks that may be directly observed, particularly at wavelengths of 10 mm or longer, but that they are more likely to reveal their presence in scattered light by puffing up their outer gap edges and exciting multiple spiral arms.

A large fraction of the gas in the Galaxy is cold, dense, and molecular. If all this gas collapsed under the influence of gravity and formed stars in a local free-fall time, the star formation rate in the Galaxy would exceed that observed by more than an order of magnitude. Other star-forming galaxies behave similarly. Yet observations and simulations both suggest that the molecular gas is indeed gravitationally collapsing, albeit hierarchically. Prompt stellar feedback offers a potential solution to the low observed star formation rate if it quickly disrupts star-forming clouds during gravitational collapse. However, this requires that molecular clouds must be short-lived objects, raising the question of how so much gas can be observed in the molecular phase. This can occur only if molecular clouds form as quickly as they are destroyed, maintaining a global equilibrium fraction of dense gas. We therefore examine cloud formation timescales. We first demonstrate that supernova and superbubble sweeping cannot produce dense gas at the rate required to match the cloud destruction rate. On the other hand, Toomre gravitational instability can reach the required production rate. We thus argue that, although dense, star-forming gas may last only around a single global free-fall time, the dense gas in star-forming galaxies can globally exist in a state of dynamic equilibrium between formation by gravitational instability, and disruption by stellar feedback. At redshift z >~ 2, the Toomre instability timescale decreases, resulting in a prediction of higher molecular gas fractions at early times, in agreement with observations.

The evolution of dwarf satellites of the Milky Way is affected by the combination of ram pressure and tidal stripping, and internal feedback from massive stars. We investigate gas loss processes in the smallest satellites of the Milky Way using three-dimensional, high resolution, idealized wind tunnel simulations, accounting for gas loss through both ram pressure stripping and expulsion by supernova feedback. Using initial conditions appropriate for a dwarf galaxy like Leo T, we investigate whether or not environmental gas stripping and internal feedback can quench these low mass galaxies on the expected timescales, shorter than 2 Gyr. We find that supernova feedback contributes negligibly to the stripping rate for these low star formation rate galaxies. However, we also find that ram pressure stripping is less efficient than expected in the stripping scenarios we consider. Our work suggests that, although ram pressure stripping can eventually completely strip these galaxies, other physics is likely at play to reconcile our computed stripping times with the rapid quenching timescales deduced from observations of low mass Milky Way dwarf galaxies. We discuss the roles additional physics may play in this scenario, including host-satellite tidal interactions, cored vs. cuspy dark matter profiles, reionization, and satellite pre-processing. We conclude that a proper accounting of these physics together is necessary to understand the quenching of low mass Milky Way satellites.

Molecular cloud observations show that clouds have non-thermal velocity dispersions that scale with the cloud size as $\sigma\propto R^{1/2}$ at constant surface density, and for varying surface density scale with both the cloud`s size and surface density, $\sigma^2 \propto R \Sigma$. The energy source driving these chaotic motions remains poorly understood. We describe the velocity dispersions observed in a cloud population formed in a kiloparsec-scale numerical simulation of a magnetized, supernova-driven, self-gravitating, interstellar medium, including diffuse heating and radiative cooling. We compare the relationships between velocity dispersion, size, and surface density measured in the simulated cloud population to those found in observations of Galactic molecular clouds. We find that external supernova explosions can not drive turbulent motions of the observed magnitudes within dense clouds. On the other hand, self-gravity also induces non-thermal motions as gravitationally bound clouds begin to collapse in our model, and by doing so their internal velocity dispersions recover the observed relations. Energy conservation suggests that the observed behavior is consistent with the kinetic energy being proportional to the gravitational energy. However, the clouds in our model show no sign of reaching a stable equilibrium state at any time, even for strongly magnetized clouds. We conclude that gravitationally bound molecular clouds are always in a state of gravitational collapse and their properties are a natural result of this chaotic collapse. In order to agree with observed star formation efficiencies, this process must be terminated by the early destruction of the clouds, presumably from internal stellar feedback.

Accretion disks around supermassive black holes (SMBHs) in active galactic nuclei contain stars, stellar mass black holes, and other stellar remnants, which perturb the disk gas gravitationally. The resulting density perturbations in turn exert torques on the embedded masses causing them to migrate through the disk in a manner analogous to the behavior of planets in protoplanetary disks. We determine the strength and direction of these torques using an empirical analytic description dependent on local disk gradients, applied to two different analytic, steady-state disk models of SMBH accretion disks. We find that there are radii in such disks where the gas torque changes sign, trapping migrating objects. Our analysis shows that major migration traps generally occur where the disk surface density gradient changes sign from positive to negative, around 20--300$R_{\rm g}$, where $R_{\rm g}=2GM/c^{2}$ is the Schwarzschild radius. At these traps, massive objects in the AGN disk can accumulate, collide, scatter, and accrete. Intermediate mass black hole formation is likely in these disk locations, which may lead to preferential gap and cavity creation at these radii. Our model thus has significant implications for SMBH growth as well as gravitational wave source populations.

We present a hydrodynamical simulation of the turbulent, magnetized, supernova (SN)-driven interstellar medium (ISM) in a stratified box that dynamically couples the injection and evolution of cosmic rays (CRs) and a self-consistent evolution of the chemical composition. CRs are treated as a relativistic fluid in the advection-diffusion approximation. The thermodynamic evolution of the gas is computed using a chemical network that follows the abundances of H+, H, H2, CO, C+, and free electrons and includes (self-)shielding of the gas and dust. We find that CRs perceptibly thicken the disk with the heights of 90% (70%) enclosed mass reaching ~1.5 kpc (~0.2 kpc). The simulations indicate that CRs alone can launch and sustain strong outflows of atomic and ionized gas with mass loading factors of order unity, even in solar neighborhood conditions and with a CR energy injection per SN of 10^50 erg, 10% of the fiducial thermal energy of an SN. The CR-driven outflows have moderate launching velocities close to the midplane (~100 km/s) and are denser (\rho~1e-24 - 1e-26 g/cm^3), smoother, and colder than the (thermal) SN-driven winds. The simulations support the importance of CRs for setting the vertical structure of the disk as well as the driving of winds.

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.

We have recently published observations of significant flux density variations at 1.3 cm in HII regions in the star forming regions Sgr B2 Main and North (De Pree et al. 2014). To further study these variations, we have made new 7 mm continuum and recombination line observations of Sgr B2 at the highest possible angular resolution of the Karl G. Jansky Very Large Array (VLA). We have observed Sgr B2 Main and North at 42.9 GHz and at 45.4 GHz in the BnA configuration (Main) and the A configuration (North). We compare these new data to archival VLA 7 mm continuum data of Sgr B2 Main observed in 2003 and Sgr B2 North observed in 2001. We find that one of the 41 known ultracompact and hypercompact HII regions in Sgr B2 (K2-North) has decreased $\sim$27% in flux density from 142$\pm$14 mJy to 103$\pm$10 mJy (2.3$\sigma$) between 2001 and 2012. A second source, F3c-Main has increased $\sim$30% in flux density from 82$\pm$8 mJy to 107 $\pm$11 mJy (1.8$\sigma$) between 2003 and 2012. F3c-Main was previously observed to increase in flux density at 1.3 cm over a longer time period between 1989 and 2012 (De Pree et al. 2014). An observation of decreasing flux density, such as that observed in K2-North, is particularly significant since such a change is not predicted by the classical hypothesis of steady expansion of HII regions during massive star accretion. Our new observations at 7 mm, along with others in the literature, suggest that the formation of massive stars occurs through time-variable and violent accretion.

Recent high-resolution, near-infrared images of protoplanetary disks have shown that these disks often present spiral features. Spiral arms are among the structures predicted decades ago by numerical simulations of disk-planet interaction and thus it is tempting to suspect that planetary perturbers are responsible for the observed signatures. However, such interpretation is not free of problems. The spirals are found to have large pitch angles, and in at least one case (HD 100546) the spiral feature appears effectively unpolarized, which implies thermal emission of the order of 1000K (465$\pm$40K at closer inspection). We have recently shown in two-dimensional models that shock dissipation in the supersonic wake of high-mass planets can lead to significant heating if the disk is sufficiently adiabatic. In this paper we extend this analysis to three dimensions in thermodynamically evolving disks. We use the Pencil Code in spherical coordinates for our models, with a prescription for thermal cooling based on the optical depth of the local vertical gas column. We use a 5$M_J$ planet, and show that shocks in the region around the planet where the Lindblad resonances occur heat the gas to substantially higher temperatures than the ambient disk gas at that radius. The gas is accelerated vertically away from the midplane by the shocks to form shock bores, and the gas falling back toward the midplane breaks up into a turbulent surf near the Lindblad resonances. This turbulence, although localized, has high $\alpha$ values, reaching 0.05 in the inner Lindblad resonance, and 0.1 in the outer one. We also find evidence that the disk regions heated up by the planetary shocks eventually becomes superadiabatic, generating convection far from the planet's orbit.

The Rotten Egg Nebula has at its core a binary composed of a Mira star and an A-type companion at a separation >10 au. It has been hypothesized to have formed by strong binary interactions between the Mira and a companion in an eccentric orbit during periastron passage ~800 years ago. We have performed hydrodynamic simulations of an asymptotic giant branch star interacting with companions with a range of masses in orbits with a range of initial eccentricities and periastron separations. For reasonable values of the eccentricity, we find that Roche lobe overflow can take place only if the periods are <<100 years. Moreover, mass transfer causes the system to enter a common envelope phase within several orbits. Since the central star of the Rotten Egg nebula is an AGB star, we conclude that such a common envelope phase must have lead to a merger, so the observed companion must have been a tertiary companion of a binary that merged at the time of nebula ejection. Based on the mass and timescale of the simulated disc formed around the companion before the common envelope phase, we analytically estimate the properties of jets that could be launched. Allowing for super-Eddington accretion rates, we find that jets similar to those observed are plausible, provided that the putative lost companion was relatively massive.

Luminous compact blue galaxies (LCBGs) have enhanced star formation rates and compact morphologies. We combine Sloan Digital Sky Survey data with HI data of 29 LCBGs at redshift z~0 to understand their nature. We find that local LCBGs have high atomic gas fractions (~50%) and star formation rates per stellar mass consistent with some high redshift star forming galaxies. Many local LCBGs also have clumpy morphologies, with clumps distributed across their disks. Although rare, these galaxies appear to be similar to the clumpy star forming galaxies commonly observed at z~1-3. Local LCBGs separate into three groups: 1. Interacting galaxies (~20%); 2. Clumpy spirals (~40%); 3. Non-clumpy, non-spirals with regular shapes and smaller effective radii and stellar masses (~40%). It seems that the method of building up a high gas fraction, which then triggers star formation, is not the same for all local LCBGs. This may lead to a dichotomy in galaxy characteristics. We consider possible gas delivery scenarios and suggest that clumpy spirals, preferentially located in clusters and with companions, are smoothly accreting gas from tidally disrupted companions and/or intracluster gas enriched by stripped satellites. Conversely, as non-clumpy galaxies are preferentially located in the field and tend to be isolated, we suggest clumpy, cold streams, which destroy galaxy disks and prevent clump formation, as a likely gas delivery mechanism for these systems. Other possibilities include smooth cold streams, a series of minor mergers, or major interactions.

Classical Sweet-Parker models of reconnection predict that reconnection rates depend inversely on the resistivity, usually parameterized using the dimensionless Lundquist number ($\Lund$). We describe magnetohydrodynamic (MHD) simulations using a static, nested grid that show the development of a three-dimensional instability in the plane of a current sheet between reversing field lines without a guide field. The instability leads to rapid reconnection of magnetic field lines at a rate independent of $\Lund$ over at least the range $3.2\times 10^3 \lesssim \Lund \lesssim 3.2 \times 10^5$ resolved by the simulations. We find that this instability occurs even for cases with $\Lund \lesssim 10^4$ that in our models appear stable to the recently described, two-dimensional, plasmoid instability. Our results suggest that three-dimensional, MHD processes alone produce fast (resistivity independent) reconnection without recourse to kinetic effects or external turbulence. The unstable reconnection layers provide a self-consistent environment in which the extensively studied turbulent reconnection process can occur.

Recent observations of gaps and non-axisymmetric features in the dust distributions of transition disks have been interpreted as evidence of embedded massive protoplanets. However, comparing the predictions of planet-disk interaction models to the observed features has shown far from perfect agreement. This may be due to the strong approximations used for the predictions. For example, spiral arm fitting typically uses results that are based on low-mass planets in an isothermal gas. In this work, we describe two-dimensional, global, hydrodynamical simulations of disks with embedded protoplanets, with and without the assumption of local isothermality, for a range of planet-to-star mass ratios 1-10 M_jup for a 1 M_sun star. We use the Pencil Code in polar coordinates for our models. We find that the inner and outer spiral wakes of massive protoplanets (M>5 M_jup) produce significant shock heating that can trigger buoyant instabilities. These drive sustained turbulence throughout the disk when they occur. The strength of this effect depends strongly on the mass of the planet and the thermal relaxation timescale; for a 10 M_jup planet embedded in a thin, purely adiabatic disk, the spirals, gaps, and vortices typically associated with planet-disk interactions are disrupted. We find that the effect is only weakly dependent on the initial radial temperature profile. The spirals that form in disks heated by the effects we have described may fit the spiral structures observed in transition disks better than the spirals predicted by linear isothermal theory.

T Pyxidis is the only recurrent nova surrounded by knots of material ejected in previous outbursts. Following the eruption that began on 2011 April 14.29, we obtained seven epochs (from 4 to 383 days after eruption) of Hubble Space Telescope narrowband Ha images of T Pyx . The flash of radiation from the nova event had no effect on the ejecta until at least 55 days after the eruption began. Photoionization of hydrogen located north and south of the central star was seen 132 days after the beginning of the eruption. That hydrogen recombined in the following 51 days, allowing us to determine a hydrogen atom density of at least 7e5 cm^-3 - at least an order of magnitude denser than the previously detected, unresolved [NII] knots surrounding T Pyx. Material to the northwest and southeast was photoionized between 132 and 183 days after the eruption began. 99 days later that hydrogen had recombined. Both then (282 days after outburst) and 101 days later, we detected almost no trace of hydrogen emission around T Pyx. There is a large reservoir of previously unseen, cold diffuse hydrogen overlapping the previously detected, [NII] - emitting knots of T Pyx ejecta. The mass of this newly detected hydrogen is probably an order of magnitude larger than that of the [NII] knots. We also determine that there is no significant reservoir of undetected ejecta from the outer boundaries of the previously detected ejecta out to about twice that distance, near the plane of the sky. The lack of distant ejecta is consistent with the Schaefer et al (2010) scenario for T Pyx, in which the star underwent its first eruption within five years of 1866 after many millennia of quiescence, followed by the six observed recurrent nova eruptions since 1890. This lack of distant ejecta is not consistent with scenarios in which T Pyx has been erupting continuously as a recurrent nova for many centuries or millennia.

Chondrules are millimeter-sized spherules that dominate primitive meteorites (chondrites) originating from the asteroid belt. The incorporation of chondrules into asteroidal bodies must be an important step in planet formation, but the mechanism is not understood. We show that the main growth of asteroids can result from gas-drag-assisted accretion of chondrules. The largest planetesimals of a population with a characteristic radius of 100 km undergo run-away accretion of chondrules within ~3 Myr, forming planetary embryos up to Mars sizes along with smaller asteroids whose size distribution matches that of main belt asteroids. The aerodynamical accretion leads to size-sorting of chondrules consistent with chondrites. Accretion of mm-sized chondrules and ice particles drives the growth of planetesimals beyond the ice line as well, but the growth time increases above the disk life time outside of 25 AU. The contribution of direct planetesimal accretion to the growth of both asteroids and Kuiper belt objects is minor. In contrast, planetesimal accretion and chondrule accretion play more equal roles for the formation of Moon-sized embryos in the terrestrial planet formation region. These embryos are isolated from each other and accrete planetesimals only at a low rate. However, the continued accretion of chondrules destabilizes the oligarchic configuration and leads to the formation of Mars-sized embryos and terrestrial planets by a combination of direct chondrule accretion and giant impacts.

We use hydrodynamical simulations in a $(256\;{\rm pc})^3$ periodic box to model the impact of supernova (SN) explosions on the multi-phase interstellar medium (ISM) for initial densities $n = 0.5-30$ cm$^{-3}$ and SN rates $1-720$ Myr$^{-1}$. We include radiative cooling, diffuse heating, and the formation of molecular gas using a chemical network. The SNe explode either at random positions, at density peaks, or both. We further present a model combining thermal energy for resolved and momentum input for unresolved SNe. Random driving at high SN rates results in hot gas ($T\gtrsim 10^6$ K) filling $> 90$% of the volume. This gas reaches high pressures ($10^4 < P/k_\mathrm{B} < 10^7$ K cm$^{-3}$) due to the combination of SN explosions in the hot, low density medium and confinement in the periodic box. These pressures move the gas from a two-phase equilibrium to the single-phase, cold branch of the cooling curve. The molecular hydrogen dominates the mass ($>50$%), residing in small, dense clumps. Such a model might resemble the dense ISM in high-redshift galaxies. Peak driving results in huge radiative losses, producing a filamentary ISM with virtually no hot gas, and a small molecular hydrogen mass fraction ($\ll 1$%). Varying the ratio of peak to random SNe yields ISM properties in between the two extremes, with a sharp transition for equal contributions. The velocity dispersion in HI remains $\lesssim 10$ km s$^{-1}$ in all cases. For peak driving the velocity dispersion in H$_\alpha$ can be as high as $70$ km s$^{-1}$ due to the contribution from young, embedded SN remnants.

We test the X-ray emission predictions of galactic fountain models against XMM-Newton measurements of the emission from the Milky Way's hot halo. These measurements are from 110 sight lines, spanning the full range of Galactic longitudes. We find that a magnetohydrodynamical simulation of a supernova-driven interstellar medium, which features a flow of hot gas from the disk to the halo, reproduces the temperature but significantly underpredicts the 0.5-2.0 keV surface brightness of the halo (by two orders of magnitude, if we compare the median predicted and observed values). This is true for versions of the model with and without an interstellar magnetic field. We consider different reasons for the discrepancy between the model predictions and the observations. We find taking into account overionization in cooled halo plasma, which could in principle boost the predicted X-ray emission, is unlikely in practice to bring the predictions in line with the observations. We also find that including thermal conduction, which would tend to increase the surface brightnesses of interfaces between hot and cold gas, would not overcome the surface brightness shortfall. However, charge exchange emission from such interfaces, not included in the current model, may be significant. The faintness of the model may also be due to the lack of cosmic ray driving, meaning that the model may underestimate the amount of material transported from the disk to halo. In addition, an extended hot halo of accreted material may be important, by supplying hot electrons that could boost the emission of the material driven out from the disk. Additional model predictions are needed to test the relative importance of these processes in explaining the observed halo emission.

Early in the study of viscous accretion disks it was realized that energy transfers from distant sources must be important, not least because the flow at the disk midplane in the bulk of the disk is likely outwards, out of the gravitational potential well. If the source of the viscosity is powered by accretion, such as in the case of the magneto-rotational instability, such distant energy sources must lie in the innermost regions of the disk, where accretion occurs even at the midplane. We argue here that modulations in this energy supply can alter the accretion rate on dynamical, rather than far longer viscous, time scales. This means that both the steady state value of and fluctuations in the inner disk's accretion rate, depending on the details of the inner boundary condition and occurring on the inner disk's rapid evolution time, can affect the outer disk. This is particularly interesting because observations have shown that disk accretion is not steady (e.g.~EX Lupi type objects). We also note that the power supplied to shearing boxes is set by the boxes themselves rather than the physical energy fluxes in a global disk. That is, their saturated magnetic field is not subject to the full set of energy constraints present in an actual disk. Our analysis suggests that large scale radial transport of energy has a critical impact on the evolution and variability of accretion disks.

The magnetorotational instability (MRI) drives magnetized turbulence in sufficiently ionized regions of protoplanetary disks, leading to mass accretion. The dissipation of the potential energy associated with this accretion determines the thermal structure of accreting regions. Until recently, the heating from the turbulence has only been treated in an azimuthally averaged sense, neglecting local fluctuations. However, magnetized turbulence dissipates its energy intermittently in current sheet structures. We study this intermittent energy dissipation using high resolution numerical models including a treatment of radiative thermal diffusion in an optically thick regime. Our models predict that these turbulent current sheets drive order unity temperature variations even where the MRI is damped strongly by Ohmic resistivity. This implies that the current sheet structures where energy dissipation occurs must be well resolved to correctly capture the flow structure in numerical models. Higher resolutions are required to resolve energy dissipation than to resolve the magnetic field strength or accretion stresses. The temperature variations are large enough to have major consequences for mineral formation in disks, including melting chondrules, remelting calcium-aluminum rich inclusions, and annealing silicates; and may drive hysteresis: current sheets in MRI active regions could be significantly more conductive than the remainder of the disk.

The formation of high-mass stars is usually accompanied by powerful protostellar outflows. Such high-mass outflows are not simply scaled-up versions of their lower-mass counterparts, since observations suggest that the collimation degree degrades with stellar mass. Theoretically, the origins of massive outflows remain open to question because radiative feedback and fragmentation of the accretion flow around the most massive stars, with M > 15 M_Sun, may impede the driving of magnetic disk winds. We here present a three-dimensional simulation of the early stages of core fragmentation and massive star formation that includes a subgrid-scale model for protostellar outflows. We find that stars that form in a common accretion flow tend to have aligned outflow axes, so that the individual jets of multiple stars can combine to form a collective outflow. We compare our simulation to observations with synthetic H_2 and CO observations and find that the morphology and kinematics of such a collective outflow resembles some observed massive outflows, such as Cepheus A and DR 21. We finally compare physical quantities derived from simulated observations of our models to the actual values in the models to examine the reliability of standard methods for deriving physical quantities, demonstrating that those methods indeed recover the actual values to within a factor of 2-3.

Accretion flows onto massive stars must transfer mass so quickly that they are themselves gravitationally unstable, forming dense clumps and filaments. These density perturbations interact with young massive stars, emitting ionizing radiation, alternately exposing and confining their HII regions. As a result, the HII regions are predicted to flicker in flux density over periods of decades to centuries rather than increasing monotonically in size as predicted by simple Spitzer solutions. We have recently observed the Sgr B2 region at 1.3 cm with the VLA in its three hybrid configurations (DnC, CnB and BnA) at a resolution of 0.25''. These observations were made to compare in detail with matched continuum observations from 1989. At 0.25'' resolution, Sgr B2 contains 41 UC HII regions, 6 of which are hypercompact. The new observations of Sgr B2 allow comparison of relative peak flux densites for the HII regions in Sgr B2 over a 23 year time baseline (1989-2012) in one of the most source-rich massive star forming regions in the Milky Way. The new 1.3 cm continuum images indicate that four of the 41 UC HII regions exhibit significant changes in their peak flux density, with one source (K3) dropping in peak flux density, and the other 3 sources (F10.303, F1 and F3) increasing in peak flux density. The results are consistent with statistical predictions from simulations of high mass star formation, suggesting that they offer a solution to the lifetime problem for ultracompact HII regions.

Giant molecular clouds (GMCs) are the primary reservoirs of cold, star-forming molecular gas in the Milky Way and similar galaxies, and thus any understanding of star formation must encompass a model for GMC formation, evolution, and destruction. These models are necessarily constrained by measurements of interstellar molecular and atomic gas, and the emergent, newborn stars. Both observations and theory have undergone great advances in recent years, the latter driven largely by improved numerical simulations, and the former by the advent of large-scale surveys with new telescopes and instruments. This chapter offers a thorough review of the current state of the field.

Hubble Space Telescope} images of the ejecta surrounding the nova T Pyxidis resolve the emission into more than two thousand bright knots. We simulate the dynamical evolution of the ejecta from T Pyxidis during its multiple eruptions over the last 150 years using the adaptive mesh refinement capability of the gas dynamics code Ramses. We demonstrate that the observed knots are the result of Richtmeyer-Meshkov gas dynamical instabilities (the equivalent of Rayleigh-Taylor instabilities in an accelerated medium). These instabilities are caused by the overrunning of the ejecta from the classical nova of 1866 by fast moving ejecta from the subsequent six recurrent nova outbursts. The model correctly predicts the observed expansion and dimming of the T Pyx ejecta as well as the knotty morphology. The model also predicts that deeper, high resolution imagery will show filamentary structure connecting the knots. We show reprocessed Hubble Space Telescope imagery that shows the first hints of such structure.

From the time the first stars formed over 13 billion years ago to the present, star formation has had an unexpectedly dynamic history. At first, the star formation rate density increased dramatically, reaching a peak 10 billion years ago more than ten times the present day value. Observations of the initial rise in star formation remain difficult, poorly constraining it. Theoretical modeling has trouble predicting this history because of the difficulty in following the feedback of energy from stellar radiation and supernova explosions into the gas from which further stars form. Observations from the ground and space with the next generation of instruments should reveal the full history of star formation in the universe, while simulations appear poised to accurately predict the observed history.

Meteoritic chondrules were formed in the early solar system by brief heating of silicate dust to melting temperatures. Some highly refractory grains (Type B calcium-aluminum-rich inclusions, CAIs) also show signs of transient heating. A similar process may occur in other protoplanetary disks, as evidenced by observations of spectra characteristic of crystalline silicates. One possible environment for this process is the turbulent magnetohydrodynamic flow thought to drive accretion in these disks. Such flows generally form thin current sheets, which are sites of magnetic reconnection, and dissipate the magnetic fields amplified by a disk dynamo. We suggest that it is possible to heat precursor grains for chondrules and other high-temperature minerals in current sheets that have been concentrated by our recently described short-circuit instability. We extend our work on this process by including the effects of radiative cooling, taking into account the temperature dependence of the opacity; and by examining current sheet geometry in three-dimensional, global models of magnetorotational instability. We find that temperatures above 1600 K can be reached for favorable parameters that match the ideal global models. This mechanism could provide an efficient means of tapping the gravitational potential energy of the protoplanetary disk to heat grains strongly enough to form high-temperature minerals. The volume-filling nature of turbulent magnetic reconnection is compatible with constraints from chondrule-matrix complementarity, chondrule-chondrule complementarity, the occurrence of igneous rims, and compound chondrules. The same short-circuit mechanism may perform other high-temperature mineral processing in protoplanetary disks such as the production of crystalline silicates and CAIs.

We simulate the evolution of a 10^9 Msun dark matter halo in a cosmological setting with an adaptive-mesh refinement code as an analogue to local low luminosity dwarf irregular and dwarf spheroidal galaxies. The primary goal of our study is to investigate the roles of reionization and supernova feedback in determining the star formation histories of low mass dwarf galaxies. We include a wide range of physical effects, including metal cooling, molecular hydrogen formation and cooling, photoionization and photodissociation from a metagalactic background, a simple prescription for self-shielding, star formation, and a simple model for supernova driven energetic feedback. We carry out simulations excluding each major effect in turn. We find that reionization is primarily responsible for expelling most of the gas in our simulations, but that supernova feedback is required to disperse the dense, cold gas in the core of the halo. Moreover, we show that the timing of reionization can produce an order of magnitude difference in the final stellar mass of the system. For our full physics run with reionization at z=9, we find a stellar mass of about 10^5 Msun at z=0, and a mass-to-light ratio within the half-light radius of approximately 130 Msun/Lsun, consistent with observed low-luminosity dwarfs. However, the resulting median stellar metallicity is 0.06 Zsun, considerably larger than observed systems. In addition, we find star formation is truncated between redshifts 4 and 7, at odds with the observed late time star formation in isolated dwarf systems but in agreement with Milky Way ultrafaint dwarf spheroidals. We investigate the efficacy of energetic feedback in our simple thermal-energy driven feedback scheme, and suggest that it may still suffer from excessive radiative losses, despite reaching stellar particle masses of about 100 Msun, and a comoving spatial resolution of 11 pc.

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 stellar and galaxy populations to models of large scale structure formation. 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. 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. 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, imply causation. It appears that both molecule formation and star formation occur as a consequence of gravitational collapse, with molecules not essential 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. (lightly abridged)

Many astrophysical systems of interest, including protoplanetary accretion disks, are made of turbu- lent magnetized gas with near solar metallicity. Thermal ionization of alkali metals in such gas exceeds non-thermal ionization when temperatures climb above roughly 1000 K. As a result, the conductiv- ity, proportional to the ionization fraction, gains a strong, positive dependence on temperature. In this paper, we demonstrate that this relation between the temperature and the conductivity triggers an exponential instability that acts similarly to an electrical short, where the increased conductivity concentrates the current and locally increases the Ohmic heating. This contrasts with the resistiv- ity increase expected in an ideal magnetic reconnection region. The instability acts to focus narrow current sheets into even narrower sheets with far higher currents and temparatures. We lay out the basic principles of this behavior in this paper using protoplanetary disks as our example host system, motivated by observations of chondritic meteorites and their ancestors, dust grains in protoplanetary disks, that reveal the existence of strong, frequent heating events that this instability could explain.

The recent discovery of two Earth-mass planets in close orbits around an evolved star has raised questions as to whether substellar companions can survive encounters with their host stars. We consider whether these companions could have been stripped of significant amounts of mass during the phase when they orbited through the dense inner envelopes of the giant. We apply the criterion derived by Murray et al. for disruption of gravitationally bound objects by ram pressure, to determine whether mass loss may have played a role in the histories of these and other recently discovered low-mass companions to evolved stars. We find that the brown dwarf and Jovian mass objects circling WD 0137-349, SDSS J08205+0008, and HIP 13044 are most unlikely to have lost significant mass during the common envelope phase. However, the Earth-mass planets found around KIC 05807616 could well be the remnant of one or two Jovian mass planets that lost extensive mass during the common envelope phase.

The formation of massive stars exceeding 10 solar masses usually results in large-scale molecular outflows. Numerical simulations, including ionization, of the formation of such stars show evidence for ionization-driven molecular outflows. We here examine whether the outflows seen in these models reproduce the observations. We compute synthetic ALMA and CARMA maps of CO emission lines of the outflows, and compare their signatures to existing single-dish and interferometric data. We find that the ionization-driven models can only reproduce weak outflows around high-mass star-forming regions. We argue that expanding H II regions probably do not represent the dominant mechanism for driving observed outflows. We suggest instead that observed outflows are driven by the collective action of the outflows from the many lower-mass stars that inevitably form around young massive stars in a cluster.

It has been suggested that the transition between magnetorotationally active and dead zones in protoplanetary disks should be prone to the excitation of vortices via Rossby wave instability (RWI). However, the only numerical evidence for this has come from alpha disk models, where the magnetic field evolution is not followed, and the effect of turbulence is parametrized by Laplacian viscosity. We aim to establish the phenomenology of the flow in the transition in 3D resistive-magnetohydrodynamical models. We model the transition by a sharp jump in resistivity, as expected in the inner dead zone boundary, using the Pencil Code to simulate the flow. We find that vortices are readily excited in the dead side of the transition. We measure the mass accretion rate finding similar levels of Reynolds stress at the dead and active zones, at the $\alpha\approx 10^{-2}$ level. The vortex sits in a pressure maximum and does not migrate, surviving until the end of the simulation. A pressure maximum in the active zone also triggers the RWI. The magnetized vortex that results should be disrupted by parasitical magneto-elliptic instabilities, yet it subsists in high resolution. This suggests that either the parasitic modes are still numerically damped, or that the RWI supplies vorticity faster than they can destroy it. We conclude that the resistive transition between the active and dead zones in the inner regions of protoplanetary disks, if sharp enough, can indeed excite vortices via RWI. Our results lend credence to previous works that relied on the alpha-disk approximation, and caution against the use of overly reduced azimuthal coverage on modeling this transition.

We present an algorithm for simulating the equations of ideal magnetohydrodynamics and other systems of differential equations on an unstructured set of points represented by sample particles. Local, third-order, least-squares, polynomial interpolations (Moving Least Squares interpolations) are calculated from the field values of neighboring particles to obtain field values and spatial derivatives at the particle position. Field values and particle positions are advanced in time with a second order predictor-corrector scheme. The particles move with the fluid, so the time step is not limited by the Eulerian Courant-Friedrichs-Lewy condition. Full spatial adaptivity is implemented to ensure the particles fill the computational volume, which gives the algorithm substantial flexibility and power. A target resolution is specified for each point in space, with particles being added and deleted as needed to meet this target. Particle addition and deletion is based on a local void and clump detection algorithm. Dynamic artificial viscosity fields provide stability to the integration. The resulting algorithm provides a robust solution for modeling flows that require Lagrangian or adaptive discretizations to resolve. This paper derives and documents the Phurbas algorithm as implemented in Phurbas version 1.1. A following paper presents the implementation and test problem results.

We present an algorithm for simulating the equations of ideal magnetohydrodynamics and other systems of differential equations on an unstructured set of points represented by sample particles. The particles move with the fluid, so the time step is not limited by the Eulerian Courant-Friedrichs-Lewy condition. Full spatial adaptivity is required to ensure the particles fill the computational volume, and gives the algorithm substantial flexibility and power. A target resolution is specified for each point in space, with particles being added and deleted as needed to meet this target. We have parallelized the code by adapting the framework provided by GADGET-2. A set of standard test problems, including 1e-6 amplitude linear MHD waves, magnetized shock tubes, and Kelvin-Helmholtz instabilities is presented. Finally we demonstrate good agreement with analytic predictions of linear growth rates for magnetorotational instability in a cylindrical geometry. This paper documents the Phurbas algorithm as implemented in Phurbas version 1.1.

The torques exerted by a locally isothermal disk on an embedded planet lead to rapid inward migration. Recent work has shown that modeling the thermodynamics without the assumption of local isothermality reveals regions where the net torque on an embedded planet is positive, leading to outward migration of the planet. When a region with negative torque lies directly exterior to this, planets in the inner region migrate outwards and planets in the outer region migrate inwards, converging where the torque is zero. We incorporate the torques from an evolving non-isothermal disk into an N-body simulation to examine the behavior of planets or planetary embryos interacting in the convergence zone. We find that mutual interactions do not eject objects from the convergence zone. Small numbers of objects in a laminar disk settle into near resonant orbits that remain stable over the 10 Myr periods that we examine. However, either or both increasing the number of planets or including a correlated, stochastic force to represent turbulence drives orbit crossings and mergers in the convergence zone. These processes can build gas giant cores with masses of order ten Earth masses from sub-Earth mass embryos in 2-3 Myr.

Stellar feedback drives the circulation of matter from the disk to the halo of galaxies. We perform three-dimensional magnetohydrodynamic simulations of a vertical column of the interstellar medium with initial conditions typical of the solar circle in which supernovae drive turbulence and determine the vertical stratification of the medium. The simulations were run using a stable, positivity-preserving scheme for ideal MHD implemented in the FLASH code. We find that the majority (\approx 90 %) of the mass is contained in thermally-stable temperature regimes of cold molecular and atomic gas at T < 200 K or warm atomic and ionized gas at 5000 K < T < 10^{4.2} K, with strong peaks in probability distribution functions of temperature in both the cold and warm regimes. The 200 - 10^{4.2} K gas fills 50-60 % of the volume near the plane, with hotter gas associated with supernova remnants (30-40 %) and cold clouds (< 10 %) embedded within. At |z| ~ 1-2 kpc, transition-temperature (10^5 K) gas accounts for most of the mass and volume, while hot gas dominates at |z| > 3 kpc. The magnetic field in our models has no significant impact on the scale heights of gas in each temperature regime; the magnetic tension force is approximately equal to and opposite the magnetic pressure, so the addition of the field does not significantly affect the vertical support of the gas. The addition of a magnetic field does reduce the fraction of gas in the cold (< 200 K) regime with a corresponding increase in the fraction of warm (~ 10^4 K) gas. However, our models lack rotational shear and thus have no large-scale dynamo, which reduces the role of the field in the models compared to reality. The supernovae drive oscillations in the vertical distribution of halo gas, with the period of the oscillations ranging from ~ 30 Myr in the T < 200 K gas to ~ 100 Myr in the 10^6 K gas, in line with predictions by Walters & Cox.

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 in other orbital properties. The magnitude of the effect on particle orbits can have 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 gravitational field of the gas and the host star. We measure the evolution of the particle orbital properties, including mean radius, eccentricity, inclination, and velocity dispersion, and its dependence on the disk properties and the particle initial conditions. Although the results converge with resolution for fixed box dimensions, we find the response of the particles to the gravity of the turbulent gas correlates with the horizontal box size, up to 16 disk scale heights. This correlation indicates that caution should be exercised when interpreting local-shearing-box models involving gravitational physics of magneto-rotational turbulence. Based on heuristic arguments, nevertheless, the criterion L_h / R ~ O(1), where L_h is the horizontal box size and R is the distance to the host star, is proposed to possibly circumvent this conundrum. If this criterion holds, we can still conclude that magneto-rotational turbulence seems likely to be ineffective at driving either diffusive migration or collisional erosion under most circumstances.

Observations of spiral galaxies show a strong linear correlation between the ratio of molecular to atomic hydrogen surface density R_mol and midplane pressure. To explain this, we simulate three-dimensional, magnetized turbulence, including simplified treatments of non-equilibrium chemistry and the propagation of dissociating radiation, to follow the formation of H_2 from cold atomic gas. The formation time scale for H_2 is sufficiently long that equilibrium is not reached within the 20-30 Myr lifetimes of molecular clouds. The equilibrium balance between radiative dissociation and H_2 formation on dust grains fails to predict the time-dependent molecular fractions we find. A simple, time-dependent model of H_2 formation can reproduce the gross behavior, although turbulent density perturbations increase molecular fractions by a factor of few above it. In contradiction to equilibrium models, radiative dissociation of molecules plays little role in our model for diffuse radiation fields with strengths less than ten times that of the solar neighborhood, because of the effective self-shielding of H_2. The observed correlation of R_mol with pressure corresponds to a correlation with local gas density if the effective temperature in the cold neutral medium of galactic disks is roughly constant. We indeed find such a correlation of R_mol with density. If we examine the value of R_mol in our local models after a free-fall time at their average density, as expected for models of molecular cloud formation by large-scale gravitational instability, our models reproduce the observed correlation over more than an order of magnitude range in density.

We present new 3D magnetohydrodynamic (MHD) simulations of a supernova-driven, stratified interstellar medium. These simulations were run using the Waagan (2009) positivity preserving scheme for ideal MHD implemented in the Flash code. The scheme is stable even for the Mach numbers approaching 100 found in this problem. We have previously shown that the density distribution arising from hydrodynamical versions of these simulations creates low-density pathways through which Lyman continuum photons can travel to heights |z| > 1 kpc. This naturally produces the warm ionized medium through photoionization due primarily to O stars near the plane. However, our earlier models reproduce the peak but not the width of the observed emission measure distribution. Here, we examine whether inclusion of magnetic fields and a greater vertical extent to the simulation domain produce a gas distribution that better matches the observations. We further study the change of magnetic energy over time in our models, showing that it appears to reach a steady state after a few hundred megayears, presumably supported by a turbulent dynamo driven by the supernova explosions.

Understanding the origin of high-mass stars is central to modern astrophysics. We shed light on this problem with simulations using a novel, adaptive-mesh, ray-tracing algorithm. These simulations consistently follow the gravitational collapse of a massive molecular cloud core, the subsequent build-up and fragmentation of the accretion disk surrounding the nascent star, and, for the first time, the interaction between its intense UV radiation field and the infalling material. We show that ionization feedback can neither stop protostellar mass growth nor suppress fragmentation. We discuss the effects of feedback by ionizing and non-ionizing radiation on the evolution of the stellar cluster. The accretion is not limited by radiative feedback but by the formation of low-mass companions in a process we call "fragmentation-induced starvation". This behavior consistently reproduces the observed relation between the most massive star and the total mass of stars in a cluster. We show that magnetic fields reduce the star formation rate and lead to the formation of more massive stars.

We use three-dimensional hydrodynamical simulations to study the rapid infall phase of the common envelope interaction of a red giant branch star of mass equal to 0.88 \msun and a companion star of mass ranging from 0.9 down to 0.1 \msun. We first compare the results obtained using two different numerical techniques with different resolutions, and find overall very good agreement. We then compare the outcomes of those simulations with observed systems thought to have gone through a common envelope. The simulations fail to reproduce those systems in the sense that most of the envelope of the donor remains bound at the end of the simulations and the final orbital separations between the donor's remnant and the companion, ranging from 26.8 down to 5.9 \rsun, are larger than the ones observed. We suggest that this discrepancy vouches for recombination playing an essential role in the ejection of the envelope and/or significant shrinkage of the orbit happening in the subsequent phase.

The magnetorotational instability (MRI) may dominate outward transport of angular momentum in accretion disks, allowing material to fall onto the central object. Previous work has established that the MRI can drive a mean-field dynamo, possibly leading to a self-sustaining accretion system. Recently, however, simulations of the scaling of the angular momentum transport parameter $\alphaSS$ with the magnetic Prandtl number $\Prandtl$ have cast doubt on the ability of the MRI to transport astrophysically relevant amounts of angular momentum in real disk systems. Here, we use simulations including explicit physical viscosity and resistivity to show that when vertical stratification is included, mean field dynamo action operates, driving the system to a configuration in which the magnetic field is not fully helical. This relaxes the constraints on the generated field provided by magnetic helicity conservation, allowing the generation of a mean field on timescales independent of the resistivity. Our models demonstrate the existence of a critical magnetic Reynolds number $\Rmagc$, below which transport becomes strongly $\Prandtl$-dependent and chaotic, but above which the transport is steady and $\Prandtl$-independent. Prior simulations showing $\Prandtl$-dependence had $\Rmag < \Rmagc$. We conjecture that this steady regime is possible because the mean field dynamo is not helicity-limited and thus does not depend on the details of the helicity ejection process. Scaling to realistic astrophysical parameters suggests that disks around both protostars and stellar mass black holes have $\Rmag >> \Rmagc$. Thus, we suggest that the strong $\Prandtl$ dependence seen in recent simulations does not occur in real systems.

Ultracompact and hypercompact HII regions appear when a star with a mass larger than about 15 solar masses starts to ionize its own environment. Recent observations of time variability in these objects are one of the pieces of evidence that suggest that at least some of them harbor stars that are still accreting from an infalling neutral accretion flow that becomes ionized in its innermost part. We present an analysis of the properties of the HII regions formed in the 3D radiation-hydrodynamic simulations presented by Peters et al. as a function of time. Flickering of the HII regions is a natural outcome of this model. The radio-continuum fluxes of the simulated HII regions, as well as their flux and size variations are in agreement with the available observations. From the simulations, we estimate that a small but non-negligible fraction (~ 10 %) of observed HII regions should have detectable flux variations (larger than 10 %) on timescales of ~ 10 years, with positive variations being more likely to happen than negative variations. A novel result of these simulations is that negative flux changes do happen, in contrast to the simple expectation of ever growing HII regions. We also explore the temporal correlations between properties that are directly observed (flux and size) and other quantities like density and ionization rates.

We present a 1.42 square degree mosaic of diffuse X-ray emission in the Great Nebula in Carina from the Chandra X-ray Observatory Advanced CCD Imaging Spectrometer camera. After removing >14,000 X-ray point sources from the field, we smooth the remaining unresolved emission, tessellate it into segments of similar apparent surface brightness, and perform X-ray spectral fitting on those tessellates to infer the intrinsic properties of the X-ray-emitting plasma. By modeling faint resolved point sources, we estimate the contribution to the extended X-ray emission from unresolved point sources and show that the vast majority of Carina's unresolved X-ray emission is truly diffuse. Line-like correlated residuals in the X-ray spectral fits suggest that substantial X-ray emission is generated by charge exchange at the interfaces between Carina's hot, rarefied plasma and its many cold neutral pillars, ridges, and clumps.

The Great Nebula in Carina provides an exceptional view into the violent massive star formation and feedback that typifies giant HII regions and starburst galaxies. We have mapped the Carina star-forming complex in X-rays, using archival Chandra data and a mosaic of 20 new 60ks pointings using the Chandra X-ray Observatory's Advanced CCD Imaging Spectrometer, as a testbed for understanding recent and ongoing star formation and to probe Carina's regions of bright diffuse X-ray emission. This study has yielded a catalog of properties of >14,000 X-ray point sources; >9800 of them have multiwavelength counterparts. Using Chandra's unsurpassed X-ray spatial resolution, we have separated these point sources from the extensive, spatially-complex diffuse emission that pervades the region; X-ray properties of this diffuse emission suggest that it traces feedback from Carina's massive stars. In this introductory paper, we motivate the survey design, describe the Chandra observations, and present some simple results, providing a foundation for the 15 papers that follow in this Special Issue and that present detailed catalogs, methods, and science results.

Massive stars disproportionately influence their surroundings. How they form has only started to become clear recently through radiation gas dynamical simulations. However, until now, no simulation has simultaneously included both magnetic fields and ionizing radiation. Here we present the results from the first radiation-magnetohydrodynamical (RMHD) simulation including ionization feedback, comparing an RMHD model of a 1000 M_sol rotating cloud to earlier radiation gas dynamical models with the same initial density and velocity distributions. We find that despite starting with a strongly supercritical mass to flux ratio, the magnetic field has three effects. First, the field offers locally support against gravitational collapse in the accretion flow, substantially reducing the amount of secondary fragmentation in comparison to the gas dynamical case. Second, the field drains angular momentum from the collapsing gas, further increasing the amount of material available for accretion by the central, massive, protostar, and thus increasing its final mass by about 50% from the purely gas dynamical case. Third, the field is wound up by the rotation of the flow, driving a tower flow. However, this flow never achieves the strength seen in low-mass star formation simulations for two reasons: gravitational fragmentation disrupts the circular flow in the central regions where the protostars form, and the expanding H II regions tend to further disrupt the field geometry. Therefore, outflows driven by ionization heating look likely to be more dynamically important in regions of massive star formation.

We study the relationship between the H2 and CO abundances in simulated molecular clouds using a fully dynamical model of magnetized turbulence coupled to a detailed chemical network. We find that the CO-to-H2 conversion factor for a given molecular cloud, the so-called X-factor, is determined primarily by the mean extinction of the cloud, rather than by its metallicity. Our results explain the discrepancy observed in low metallicity systems between cloud masses derived from CO observations and other techniques such as infrared emission, and predict that CO-bright clouds in low metallicity systems should be systematically larger and/or denser than Milky Way clouds.

The {\alpha}-formalism is a common way to parametrize the common envelope interaction between a giant star and a more compact companion. The {\alpha} parameter describes the fraction of orbital energy released by the companion that is available to eject the giant star's envelope. By using new, detailed stellar evolutionary calculations we derive a user-friendly prescription for the {\lambda} parameter and an improved approximation for the envelope binding energy, thus revising the {\alpha} equation. We then determine {\alpha} both from simulations and observations in a self consistent manner. By using our own stellar structure models as well as population considerations to reconstruct the primary's parameters at the time of the common envelope interaction, we gain a deeper understanding of the uncertainties. We find that systems with very low values of q (the ratio of the companion's mass to the mass of the primary at the time of the common envelope interaction) have higher values of {\alpha}. A fit to the data suggests that lower mass companions are left at comparable or larger orbital separations to more massive companions. We conjecture that lower mass companions take longer than a stellar dynamical time to spiral in to the giant's core, and that this is key to allowing the giant to use its own thermal energy to help unbind its envelope. As a result, although systems with light companions might not have enough orbital energy to unbind the common envelope, they might stimulate a stellar reaction that results in the common envelope ejection.

The most usual tracer of molecular gas is line emission from CO. However, the reliability of that tracer has long been questioned in environments different from the Milky Way. We study the relationship between H2 and CO abundances using a fully dynamical model of magnetized turbulence coupled to a chemical network simplified to follow only the dominant pathways for H2 and CO formation and destruction, and including photodissociation using a six-ray approximation. We find that the abundance of H2 is primarily determined by the amount of time available for its formation, which is proportional to the product of the density and the metallicity, but insensitive to photodissociation. Photodissociation only becomes important at extinctions under a few tenths of a visual magnitude, in agreement with both observational and prior theoretical work. On the other hand, CO forms quickly, within a dynamical time, but its abundance depends primarily on photodissociation, with only a weak secondary dependence on H2 abundance. As a result, there is a sharp cutoff in CO abundance at mean visual extinctions A_V < 3. At lower values of A_V we find that the ratio of H2 column density to CO emissivity X_CO is proportional to A_V^(-3.5). This explains the discrepancy observed in low metallicity systems between cloud masses derived from CO observations and other techniques such as infrared emission. Our work predicts that CO-bright clouds in low metallicity systems should be systematically larger or denser than Milky Way clouds, or both. Our results further explain the narrow range of observed molecular cloud column densities as a threshold effect, without requiring the assumption of virial equilibrium.

Massive stars influence their surroundings through radiation, winds, and supernova explosions far out of proportion to their small numbers. However, the physical processes that initiate and govern the birth of massive stars remain poorly understood. Two widely discussed models are monolithic collapse of molecular cloud cores and competitive accretion. To learn more about massive star formation, we perform simulations of the collapse of rotating, massive, cloud cores including radiative heating by both non-ionizing and ionizing radiation using the FLASH adaptive mesh refinement code. These simulations show fragmentation from gravitational instability in the enormously dense accretion flows required to build up massive stars. Secondary stars form rapidly in these flows and accrete mass that would have otherwise been consumed by the massive star in the center, in a process that we term fragmentation-induced starvation. This explains why massive stars are usually found as members of high-order stellar systems that themselves belong to large clusters containing stars of all masses. The radiative heating does not prevent fragmentation, but does lead to a higher Jeans mass, resulting in fewer and more massive stars than would form without the heating. This mechanism reproduces the observed relation between the total stellar mass in the cluster and the mass of the largest star. It predicts strong clumping and filamentary structure in the center of collapsing cores, as has recently been observed. We speculate that a similar mechanism will act during primordial star formation.

We compare the predictions of three physical models for the origin of the hot halo gas with the observed halo X-ray emission, derived from 26 high-latitude XMM-Newton observations of the soft X-ray background between $l=120\degr$ and $l=240\degr$. These observations were chosen from a much larger set of observations as they are expected to be the least contaminated by solar wind charge exchange emission. We characterize the halo emission in the XMM-Newton band with a single-temperature plasma model. We find that the observed halo temperature is fairly constant across the sky (~1.8e6-2.3e6 K), whereas the halo emission measure varies by an order of magnitude (~0.0005-0.006 cm^-6 pc). When we compare our observations with the model predictions, we find that most of the hot gas observed with XMM-Newton does not reside in isolated extraplanar supernova remnants -- this model predicts emission an order of magnitude too faint. A model of a supernova-driven interstellar medium, including the flow of hot gas from the disk into the halo in a galactic fountain, gives good agreement with the observed 0.4-2.0 keV surface brightness. This model overpredicts the halo X-ray temperature by a factor of ~2, but there are a several possible explanations for this discrepancy. We therefore conclude that a major (possibly dominant) contributor to the halo X-ray emission observed with XMM-Newton is a fountain of hot gas driven into the halo by disk supernovae. However, we cannot rule out the possibility that the extended hot halo of accreted material predicted by disk galaxy formation models also contributes to the emission.

We investigate models for the photoionization of the widespread diffuse ionized gas in galaxies. In particular we address the long standing question of the penetration of Lyman continuum photons from sources close to the galactic midplane to large heights in the galactic halo. We find that recent hydrodynamical simulations of a supernova-driven interstellar medium have low density paths and voids that allow for ionizing photons from midplane OB stars to reach and ionize gas many kiloparsecs above the midplane. We find ionizing fluxes throughout our simulation grids are larger than predicted by one dimensional slab models, thus allowing for photoionization by O stars of low altitude neutral clouds in the Galaxy that are also detected in Halpha. In previous studies of such clouds the photoionization scenario had been rejected and the Halpha had been attributed to enhanced cosmic ray ionization or scattered light from midplane H II regions. We do find that the emission measure distributions in our simulations are wider than those derived from Halpha observations in the Milky Way. In addition, the horizontally averaged height dependence of the gas density in the hydrodynamical models is lower than inferred in the Galaxy. These discrepancies are likely due to the absence of magnetic fields in the hydrodynamic simulations and we discuss how magnetohydrodynamic effects may reconcile models and observations. Nevertheless, we anticipate that the inclusion of magnetic fields in the dynamical simulations will not alter our primary finding that midplane OB stars are capable of producing high altitude diffuse ionized gas in a realistic three-dimensional interstellar medium.

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 HII regions forming around massive stars are initially trapped by the infalling gas, but soon begin to fluctuate rapidly. Over time, the same ultracompact HII region can expand anisotropically, contract again, and take on any of the observed morphological classes. The total lifetime of HII 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 HII 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 spatial morphology, spectral characteristics, and time variability of ultracompact H II regions provide strong constraints on the process of massive star formation. We have performed simulations of the gravitational collapse of rotating molecular cloud cores, including treatments of the propagation of ionizing and non-ionizing radiation. We here present synthetic radio continuum observations of H II regions from our collapse simulations, to investigate how well they agree with observation, and what we can learn about how massive star formation proceeds. We find that intermittent shielding by dense filaments in the gravitationally unstable accretion flow around the massive star leads to highly variable H II regions that do not grow monotonically, but rather flicker, growing and shrinking repeatedly. This behavior appears able to resolve the well-known lifetime problem. We find that multiple ionizing sources generally form, resulting in groups of ultracompact H II regions, consistent with observations. We confirm that our model reproduces the qualitative H II region morphologies found in surveys, with generally consistent relative frequencies. We also find that simulated spectral energy distributions (SEDs) from our model are consistent with the range of observed H II region SEDs, including both regions showing a normal transition from optically thick to optically thin emission, and those with intermediate spectral slopes. In our models, anomalous slopes are solely produced by inhomogeneities in the H II region, with no contribution from dust emission at millimeter or submillimeter wavelengths. We conclude that many observed characteristics of ultracompact H II regions appear consistent with massive star formation in fast, gravitationally unstable, accretion flows.

Outward migration of low-mass planets has recently been shown to be a possibility in non-barotropic disks. We examine the consequences of this result in evolutionary models of protoplanetary disks. Planet migration occurs towards equilibrium radii with zero torque. These radii themselves migrate inwards because of viscous accretion and photoevaporation. We show that as the surface density and temperature fall, the planet orbital migration and disk depletion timescales eventually become comparable, with the precise timing depending on the mass of the planet. When this occurs, the planet decouples from the equilibrium radius. At this time, however, the gas surface density is already too low to drive substantial further migration. A higher mass planet, of 10 Earth masses, can open a gap during the late evolution of the disk, and stops migrating. Low mass planets, with 1 or 0.1 Earth masses, released beyond 1 AU in our models, avoid migrating into the star. Our results provide support for the reduced migration rates adopted in recent planet population synthesis models.

We study two limiting cases of turbulence forcing in numerical experiments: solenoidal (divergence-free) forcing, and compressive (curl-free) forcing, and compare our results to observations reported in the literature. We solve the equations of hydrodynamics on grids with up to 1024^3 cells for purely solenoidal and purely compressive forcing. Eleven lower-resolution models with mixtures of both forcings are also analysed. We find velocity dispersion--size relations consistent with observations and independent numerical simulations, irrespective of the type of forcing. However, compressive forcing yields stronger turbulent compression at the same RMS Mach number than solenoidal forcing, resulting in a three times larger standard deviation of volumetric and column density probability distributions (PDFs). We conclude that the strong dependence of the density PDF on the type of forcing must be taken into account in any theory using the PDF to predict properties of star formation. We supply a quantitative description of this dependence. We find that different observed regions show evidence of different mixtures of compressive and solenoidal forcing, with more compressive forcing occurring primarily in swept-up shells.

We describe the first three-dimensional simulation of the gravitational collapse of a massive, rotating molecular cloud that includes heating by both non-ionizing and ionizing radiation. We find that as the first protostars gain sufficient mass to ionize the accretion flow, their H II regions are initially gravitationally trapped, but soon begin to rapidly fluctuate between trapped and extended states, in agreement with observations. Over time, the same ultracompact H II region can expand anisotropically, contract again, and take on any of the observed morphological classes. In their extended phases, expanding H II regions drive bipolar neutral outflows characteristic of high-mass star formation. The total lifetime of H II regions is given by the global accretion timescale, rather than their short internal sound-crossing time. The pressure of the hot, ionized gas does not terminate accretion. Instead the final stellar mass is set by fragmentation-induced starvation. Local gravitational instabilities in the accretion flow lead to the build-up of a small cluster of stars, all with relatively high masses due to heating from accretion radiation. These companions subsequently compete with the initial high-mass star for the same common gas reservoir and limit its mass growth. Our findings show that 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, ionized flow.

The dynamics of planetesimals and planetary cores may be strongly influenced by density perturbations driven by magneto-rotational turbulence in their natal protoplanetary gas disks. Using the local shearing box approximation, we perform numerical simulations of planetesimals moving as massless particles in a turbulent, magnetized, unstratified gas disk. Our fiducial disk model shows turbulent accretion characterized by a Shakura-Sunyaev viscosity parameter of $\alpha \sim 10^{-2}$, with root-mean-square density perturbations of $\sim$10%. We measure the statistical evolution of particle orbital properties in our simulations including mean radius, eccentricity, and velocity dispersion. We confirm random walk growth in time of all three properties, the first time that this has been done with direct orbital integration in a local model. We find that the growth rate increases with the box size used at least up to boxes of eight scale heights in horizontal size. However, even our largest boxes show velocity dispersions sufficiently low that collisional destruction of planetesimals should be unimportant in the inner disk throughout its lifetime. Our direct integrations agree with earlier torque measurements showing that type I migration dominates over diffusive migration by stochastic torques for most objects in the planetary core and terrestrial planet mass range. Diffusive migration remains important for objects in the mass range of kilometer-sized planetesimals. Discrepancies in the derived magnitude of turbulence between local and global simulations of magneto-rotationally unstable disks remains an open issue, with important consequences for planet formation scenarios.

We analyze the statistics and star formation rate obtained in high-resolution numerical experiments of forced supersonic turbulence, and compare with observations. We concentrate on a systematic comparison of solenoidal (divergence-free) and compressive (curl-free) forcing, which are two limiting cases of turbulence driving. Our results show that for the same RMS Mach number, compressive forcing produces a three times larger standard deviation of the density probability distribution. When self-gravity is included in the models, the star formation rate is more than one order of magnitude higher for compressive forcing than for solenoidal forcing.

Supersonic turbulence in molecular clouds is a dominant agent that strongly affects the clouds' evolution and star formation activity. Turbulence may be initiated and maintained by a number of processes, acting at a wide range of physical scales. By examining the dynamical state of molecular clouds, it is possible to assess the primary candidates for how the turbulent energy is injected. The aim of this paper is to constrain the scales at which turbulence is driven in the molecular interstellar medium, by comparing simulated molecular spectral line observations of numerical magnetohydrodynamic (MHD) models and molecular spectral line observations of real molecular clouds. We use principal component analysis, applied to both models and observational data, to extract a quantitative measure of the driving scale of turbulence. We find that only models driven at large scales (comparable to, or exceeding, the size of the cloud) are consistent with observations. This result applies also to clouds with little or no internal star formation activity. Astrophysical processes acting on large scales, including supernova-driven turbulence, magnetorotational instability, or spiral shock forcing, are viable candidates for the generation and maintenance of molecular cloud turbulence. Small scale driving by sources internal to molecular clouds, such as outflows, can be important on small scales, but cannot replicate the observed large-scale velocity fluctuations in the molecular interstellar medium.

We present results from high-resolution three-dimensional simulations of turbulent interstellar gas that self-consistently follow its coupled thermal, chemical and dynamical evolution, with a particular focus on the formation and destruction of H2 and CO. We quantify the formation timescales for H2 and CO in physical conditions corresponding to those found in nearby giant molecular clouds, and show that both species form rapidly, with chemical timescales that are comparable to the dynamical timescale of the gas. We also investigate the spatial distributions of H2 and CO, and how they relate to the underlying gas distribution. We show that H2 is a good tracer of the gas distribution, but that the relationship between CO abundance and gas density is more complex. The CO abundance is not well-correlated with either the gas number density n or the visual extinction A_V: both have a large influence on the CO abundance, but the inhomogeneous nature of the density field produced by the turbulence means that n and A_V are only poorly correlated. There is a large scatter in A_V, and hence CO abundance, for gas with any particular density, and similarly a large scatter in density and CO abundance for gas with any particular visual extinction. This will have important consequences for the interpretation of the CO emission observed from real molecular clouds. Finally, we also examine the temperature structure of the simulated gas. We show that the molecular gas is not isothermal. Most of it has a temperature in the range of 10--20 K, but there is also a significant fraction of warmer gas, located in low-extinction regions where photoelectric heating remains effective.

Contradicting results have been reported in the literature with respect to the performance of the numerical techniques employed for the study of supersonic turbulence. We aim at characterising the performance of different particle-based and grid-based techniques on the modelling of decaying supersonic turbulence. Four different grid codes (ENZO, FLASH, TVD, ZEUS) and three different SPH codes (GADGET, PHANTOM, VINE) are compared. We additionally analysed two calculations denoted as PHANTOM A and PHANTOM B using two different implementations of artificial viscosity. Our analysis indicates that grid codes tend to be less dissipative than SPH codes, though details of the techniques used can make large differences in both cases. For example, the Morris & Monaghan viscosity implementation for SPH results in less dissipation (PHANTOM B and VINE versus GADGET and PHANTOM A). For grid codes, using a smaller diffusion parameter leads to less dissipation, but results in a larger bottleneck effect (our ENZO versus FLASH runs). As a general result, we find that by using a similar number of resolution elements N for each spatial direction means that all codes (both grid-based and particle-based) show encouraging similarity of all statistical quantities for isotropic supersonic turbulence on spatial scales k<N/32 (all scales resolved by more than 32 grid cells), while scales smaller than that are significantly affected by the specific implementation of the algorithm for solving the equations of hydrodynamics. At comparable numerical resolution, the SPH runs were on average about ten times more computationally intensive than the grid runs, although with variations of up to a factor of ten between the different SPH runs and between the different grid runs. (abridged)

We present three-dimensional numerical simulations of particle clumping and planetesimal formation in protoplanetary disks with varying amounts of solid material. As centimeter-size pebbles settle to the mid-plane, turbulence develops through vertical shearing and streaming instabilities. We find that when the pebble-to-gas column density ratio is 0.01, corresponding roughly to solar metallicity, clumping is weak, so the pebble density rarely exceeds the gas density. Doubling the column density ratio leads to a dramatic increase in clumping, with characteristic particle densities more than ten times the gas density and maximum densities reaching several thousand times the gas density. This is consistent with unstratified simulations of the streaming instability that show strong clumping in particle dominated flows. The clumps readily contract gravitationally into interacting planetesimals of order 100 km in radius. Our results suggest that the correlation between host star metallicity and exoplanets may reflect the early stages of planet formation. We further speculate that initially low metallicity disks can be particle enriched during the gas dispersal phase, leading to a late burst of planetesimal formation.

We investigate fluid motions near the midplane of vertically stratified accretion disks with highly resistive midplanes. In such disks, the magnetorotational instability drives turbulence in thin layers surrounding a resistive, stable dead zone. The turbulent layers in turn drive motions in the dead zone. We examine the properties of these motions using three-dimensional, stratified, local, shearing-box, non-ideal, magnetohydrodynamical simulations. Although the turbulence in the active zones provides a source of vorticity to the midplane, no evidence for coherent vortices is found in our simulations. It appears that this is because of strong vertical oscillations in the dead zone. By analyzing time series of azimuthally-averaged flow quantities, we identify an axisymmetric wave mode particular to models with dead zones. This mode is reduced in amplitude, but not suppressed entirely, by changing the equation of state from isothermal to ideal. These waves are too low-frequency to affect sedimentation of dust to the midplane, but may have significance for the gravitational stability of the resulting midplane dust layers.

Feedback from massive stars is one of the least understood aspects of galaxy formation. We perform a suite of vertically stratified local interstellar medium (ISM) simulations in which supernova rates and vertical gas column densities are systematically varied based on the Schmidt-Kennicutt law. Our simulations have a sufficiently high spatial resolution (1.95 pc) to follow the hydrodynamic interactions among multiple supernovae that structure the ISM. At a given supernova rate, we find that the mean mass-weighted sound speed and velocity dispersion decrease as the inverse square root of gas density, indicating that both thermal and turbulent pressures are nearly constant in the midplane, so the effective equation of state is isobaric. In contrast, across our four models having supernova rates that range from one to 512 times the Galactic supernova rate, the mass-weighted velocity dispersion remains in the range 4-6 km/s. Hence, gas averaged over ~100 pc regions follows an effective equation of state that is close to isothermal. Simulated H I emission lines have widths of 10-18 km/s, comparable to observed values. In our highest supernova rate model, superbubble blow-outs occur, and the turbulent pressure on large scales is >~4 times higher than the thermal pressure. We find a tight correlation between the thermal and turbulent pressures averaged over ~100 pc regions in the midplane of each model, as well as across the four ISM models. We construct a subgrid model for turbulent pressure based on analytic arguments and explicitly calibrate it against our stratified ISM simulations. The subgrid model provides a simple yet physically motivated way to include supernova feedback in cosmological simulations.

We study the origin of Na I absorbing gas in ultraluminous infrared galaxies motivated by the recent observations by Martin of extremely superthermal linewidths in this cool gas. We model the effects of repeated supernova explosions driving supershells in the central regions of molecular disks with M_d=10^10 M_\sun, using cylindrically symmetric gas dynamical simulations run with ZEUS-3D. The shocked swept-up shells quickly cool and fragment by Rayleigh-Taylor instability as they accelerate out of the dense, stratified disks. The numerical resolution of the cooling and compression at the shock fronts determines the peak shell density, and so the speed of Rayleigh-Taylor fragmentation. We identify cooled shells and shell fragments as Na I absorbing gas and study its kinematics. We find that simulations with a numerical resolution of \le 0.2 pc produce multiple Rayleigh-Taylor fragmented shells in a given line of sight. We suggest that the observed wide Na I absorption lines, <v> = 320 \pm 120 km s^-1 are produced by these multiple fragmented shells traveling at different velocities. We also suggest that some shell fragments can be accelerated above the observed average terminal velocity of 750 km s^-1 by the same energy-driven wind with an instantaneous starburst of \sim 10^9 M_\sun. The bulk of mass is traveling with the observed average shell velocity 330 \pm 100 km s^-1. Our results show that an energy-driven bubble causing Rayleigh-Taylor instabilities can explain the kinematics of cool gas seen in the Na I observations without invoking additional physics relying primarily on momentum conservation, such as entrainment of gas by Kelvin-Helmholtz instabilities, ram pressure driving of cold clouds by a hot wind, or radiation pressure acting on dust. (abridged)

Most stars are born in rich young stellar clusters (YSCs) embedded in giant molecular clouds. The most massive stars live out their short lives there, profoundly influencing their natal environments by ionizing HII regions, inflating wind-blown bubbles, and soon exploding as supernovae. Thousands of lower-mass pre-main sequence stars accompany the massive stars, and the expanding HII regions paradoxically trigger new star formation as they destroy their natal clouds. While this schematic picture is established, our understanding of the complex astrophysical processes involved in clustered star formation have only just begun to be elucidated. The technologies are challenging, requiring both high spatial resolution and wide fields at wavelengths that penetrate obscuring molecular material and remove contaminating Galactic field stars. We outline several important projects for the coming decade: the IMFs and structures of YSCs; triggered star formation around YSC; the fate of OB winds; the stellar populations of Infrared Dark Clouds; the most massive star clusters in the Galaxy; tracing star formation throughout the Galactic Disk; the Galactic Center region and YSCs in the Magellanic Clouds. Programmatic recommendations include: developing a 30m-class adaptive optics infrared telescope; support for high-resolution and wide field X-ray telescopes; large-aperture sub-millimeter and far-infrared telescopes; multi-object infrared spectrographs; and both numerical and analytical theory.

We explore what dominant physical mechanism sets the kinetic energy contained in neutral, atomic (HI) gas. We compare the HI line widths predicted from turbulence driven by supernova (SN) explosions and magneto-rotational instability (MRI) to direct observations in 11 disk galaxies. We use high-quality maps of the HI mass surface density and line width, obtained by the THINGS survey. We show that all sample galaxies exhibit a systematic radial decline in the HI line width, which appears to be a generic property of HI disks and also implies a radial decline in kinetic energy density of HI. At a galactocentric radius of r25 there is a characteristic value of the HI velocity dispersion of $10\pm2$ \kms. Inside this radius, galaxies show HI line widths above the thermal value expected from a warm HI component, implying that turbulence drivers must be responsible for maintaining this line width. Therefore, we compare maps of HI kinetic energy to maps of the star formation rate (SFR) and to predictions for energy generated by MRI. We find a positive correlation between kinetic energy of HI and SFR. For a given turbulence dissipation timescale we can estimate the energy input required to maintain the observed kinetic energy. The SN rate implied by the observed recent SFR is sufficient to maintain the observed velocity dispersion, if the SN feedback efficiency is at least \epsilon_SN\simeq0.1. Beyond r25, this efficiency would have to increase to unrealistic values, $\epsilon>1$, suggesting that mechanical energy from young stars does not supply most energy in outer disks. On the other hand, both thermal broadening and turbulence driven by MRI can produce the velocity dispersions and kinetic energies that we observe in this regime.

Great strides have been made in the last two decades in determining how galaxies evolve from their initial dark matter seeds to the complex structures we observe at z=0. The role of mergers has been documented through both observations and simulations, numerous satellites that may represent these initial dark matter seeds have been discovered in the Local Group, high redshift galaxies have been revealed with monstrous star formation rates, and the gaseous cosmic web has been mapped through absorption line experiments. Despite these efforts, the dark matter simulations that include baryons are still unable to accurately reproduce galaxies. One of the major problems is our incomplete understanding of how a galaxy accretes its baryons and subsequently forms stars. Galaxy formation simulations have been unable to accurately represent the required gas physics on cosmological timescales, and observations have only just begun to detect the star formation fuel over a range of redshifts and environments. How galaxies obtain gas and subsequently form stars is a major unsolved, yet tractable problem in contemporary extragalactic astrophysics. In this paper we outline how progress can be made in this area in the next decade.

Stellar feedback in galactic bulges plays an essential role in shaping the evolution of galaxies. To quantify this role and facilitate comparisons with X-ray observations, we conduct 3D hydrodynamical simulations with the adaptive mesh refinement code, FLASH, to investigate the physical properties of hot gas inside a galactic bulge, similar to that of our Galaxy or M31. We assume that the dynamical and thermal properties of the hot gas are dominated by mechanical energy input from SNe, primarily Type Ia, and mass injection from evolved stars as well as iron enrichment from SNe. We study the bulge-wide outflow as well as the SN heating on scales down to ~4 pc. An embedding scheme that is devised to plant individual SNR seeds, allows to examine, for the first time, the effect of sporadic SNe on the density, temperature, and iron ejecta distribution of the hot gas as well as the resultant X-ray morphology and spectrum. We find that the SNe produce a bulge wind with highly filamentary density structures and patchy ejecta. Compared with a 1D spherical wind model, the non-uniformity of simulated gas density, temperature, and metallicity substantially alters the spectral shape and increases the diffuse X-ray luminosity. The differential emission measure as a function of temperature of the simulated gas exhibits a log-normal distribution, with a peak value much lower than that of the corresponding 1D model. The bulk of the X-ray emission comes from the relatively low temperature and low abundance gas shells associated with SN blastwaves. SN ejecta are not well mixed with the ambient medium, at least in the bulge region. These results, at least partly, account for the apparent lack of evidence for iron enrichment in the soft X-ray-emitting gas in galactic bulges and intermediate-mass elliptical galaxies.[...]

Primordial star formation appears to result in stars at least an order of magnitude more massive than modern star formation. It has been proposed that the transition from primordial to modern initial mass functions occurs due to the onset of effective metal line cooling at a metallicity Z/Z_sun = 10^{-3.5}. However, these simulations neglected molecular hydrogen cooling. We perform simulations using the same initial conditions, but including molecular cooling, using a complex network that follows molecular hydrogen formation and also directly follows carbon monoxide and water. We find that molecular hydrogen cooling allows roughly equivalent fragmentation to proceed even at zero metallicity for these initial conditions. The apparent transition just represents the point where metal line cooling becomes more important than molecular cooling. In all cases, the fragments are massive enough to be consistent with models of primordial stellar masses, suggesting that the transition to the modern initial mass function may be determined by other physics such as dust formation. We conclude that such additional cooling mechanisms, combined with the exact initial conditions produced by cosmological collapse are likely more important than metal line cooling in determining the initial mass function, and thus that there is unlikely to be a sharp transition in the initial mass function at Z/Z_sun = 10^{-3.5}.

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. It has been proposed that this occurs when effective metal line cooling occurs at a metallicity threshold of Z/Z_sun > 10^{-3.5}. We study the influence of low levels of metal enrichment on the cooling and collapse of initially ionized gas in small protogalactic halos using three-dimensional, smoothed particle hydrodynamics simulations with particle splitting. Our initial conditions represent protogalaxies forming within a previously ionized H ii region that has not yet had time to cool and recombine. These differ considerably from those used in simulations predicting a metallicity threshold, where the gas was initially cold and only partially ionized. In the centrally condensed potential that we study here, a wide variety of initial conditions for the gas yield a monolithic central collapse. Our models show no fragmentation during collapse to number densities as high as 10^5 cm^{-3}, for metallicities reaching as high as 10^{-1} Z_sun in one rotating case, far above the threshold suggested by previous work. Rotation allows for the formation of gravitationally stable gas disks over large fractions of the local Hubble time. Turbulence slows the growth of the central density slightly, but both spherically symmetric and turbulent initial conditions collapse and form a single sink particle. We therefore argue that fragmentation at moderate density depends on the initial conditions for star formation more than on the metal abundances present.

Finite-difference simulations of fluid dynamics and magnetohydrodynamics generally require an explicit diffusion operator, either to maintain stability by attenuating grid-scale structure, or to implement physical diffusivities such as viscosity or resistivity. If the goal is stability only, the diffusion must act at the grid scale, but should affect structure at larger scales as little as possible. For physical diffusivities the diffusion scale depends on the problem, and diffusion may act at larger scales as well. Diffusivity undesirably limits the computational timestep in both cases. We construct tuned finite-difference diffusion operators that minimally limit the timestep while acting as desired near the diffusion scale. Such operators reach peak values at the diffusion scale rather than at the grid scale, but behave as standard operators at larger scales. We focus on the specific applications of hyperdiffusivity for numerical stabilization, and high Schmidt and high Prandtl number simulations where the diffusion scale greatly exceeds the grid scale.

We estimate a characteristic timescale for star formation in the spiral arms of disk galaxies, going from atomic hydrogen (HI) to dust-enshrouded massive stars. Drawing on high-resolution HI data from The HI Nearby Galaxy Survey and 24$\mu$m images from the Spitzer Infrared Nearby Galaxies Survey we measure the average angular offset between the HI and 24$\mu$m emissivity peaks as a function of radius, for a sample of 14 nearby disk galaxies. We model these offsets assuming an instantaneous kinematic pattern speed, $\Omega_p$, and a timescale, t(HI-->24$\mu$m), for the characteristic time span between the dense \hi phase and the formation of massive stars that heat the surrounding dust. Fitting for $\Omega_p$ and t(HI-->24$\mu$m), we find that the radial dependence of the observed angular offset (of the \hi and 24$\mu$m emission) is consistent with this simple prescription; the resulting corotation radii of the spiral patterns are typically $R_{cor}\simeq 2.7 R_{s}$, consistent with independent estimates. The resulting values of t(HI-->24$\mu$m) for the sample are in the range 1--4 Myr. We have explored the possible impact of non-circular gas motions on the estimate of t(HI-->24$\mu$m) and have found it to be substantially less than a factor of 2. This implies that a short timescale for the most intense phase of the ensuing star formation in spiral arms, and implies a considerable fraction of molecular clouds exist only for a few Myr before forming stars. However, our analysis does not preclude that some molecular clouds persist considerably longer. If much of the star formation in spiral arms occurs within this short interval t(HI-->24$\mu$m), then star formation must be inefficient, in order to avoid the short-term depletion of the gas reservoir.

I review the effects of magnetic fields on star formation in galaxies. This includes the effects of the magnetorotational instability (MRI) at galactic scales, magneto-Jeans and swing instabilities, Parker instabilities, and the effects of magnetic fields on the evolution of supernova-driven turbulence. I argue that currently turbulent support by the MRI appears likely to be the most important of these processes to regulating star formation.

(Abridged) The formation of molecular clouds (MCs) from the diffuse interstellar gas appears to be a necessary step for star formation, as young stars invariably occur within them. However, the mechanisms controlling the formation of MCs remain controversial. In this contribution, we focus on their formation in compressive flows driven by interstellar turbulence and large-scale gravitational instability. Turbulent compression driven by supernovae appears insufficient to explain the bulk of cloud and star formation. Rather, gravity must be important at all scales, driving the compressive flows that form both clouds and cores. Cooling and thermal instability allow the formation of dense gas out of moderate, transonic compressions in the warm diffuse gas, and drive turbulence into the dense clouds. MCs may be produced by an overshoot beyond the thermal-pressure equilibrium between the cold and warm phases of atomic gas, caused by some combination of the ram pressure of compression and the self-gravity of the compressed gas. In this case, properties of the clouds such as their mass, mass-to-magnetic flux ratio, and total kinetic and gravitational energies are in general time-variable quantities. MCs may never enter a quasi-equilibrium or virial equilibrium state but rather continuously collapse to stars.

I review the evidence for the importance of feedback from massive stars at small and large scales. The feedback mechanisms include accretion luminosity, ionizing radiation, collimated outflows, and stellar winds. The good news is that feedback doesn't entirely prevent the formation of massive stars, while the bad news is that we don't know what does limit their masses. Feedback from massive stars also influences their surroundings. I argue that this does not produce a triggering efficiency above unity, nor does it prevent lots of prompt star formation in GMCs, though it may preserve massive remnants of the clouds for many dynamical times.

Numerical simulations including magnetic fields have become important in many fields of astrophysics. Evolution of magnetic fields by the constrained transport algorithm preserves magnetic divergence to machine precision, and thus represents one preferred method for the inclusion of magnetic fields in simulations. We show that constrained transport can be implemented with volume-centered fields and hyperresistivity on a high-order finite difference stencil. Additionally, the finite-difference coefficients can be tuned to enhance high-wavenumber resolution. Similar techniques can be used for the interpolations required for dealiasing corrections at high wavenumber. Together, these measures yield an algorithm with a wavenumber resolution that approaches the theoretical maximum achieved by spectral algorithms. Because this algorithm uses finite differences instead of fast Fourier transforms, it runs faster and isn't restricted to periodic boundary conditions. Also, since the finite differences are spatially local, this algorithm is easily scalable to thousands of processors. We demonstrate that, for low-Mach-number turbulence, the results agree well with a high-order, non-constrained-transport scheme with Poisson divergence cleaning.

In cold dark matter cosmological models, the first stars to form are believed to do so within small protogalaxies. We wish to understand how the evolution of these early protogalaxies changes once the gas forming them has been enriched with small quantities of heavy elements, which are produced and dispersed into the intergalactic medium by the first supernovae. Our initial conditions represent protogalaxies forming within a fossil H II region, a previously ionized region that has not yet had time to cool and recombine. We study the influence of low levels of metal enrichment on the cooling and collapse of ionized gas in small protogalactic halos using three-dimensional, smoothed particle hydrodynamics (SPH) simulations that incorporate the effects of the appropriate chemical and thermal processes. Our previous simulations demonstrated that for metallicities Z < 0.001 Z_sun, metal line cooling alters the density and temperature evolution of the gas by less than 1% compared to the metal-free case at densities below 1 cm-3) and temperatures above 2000 K. Here, we present the results of high-resolution simulations using particle splitting to improve resolution in regions of interest. These simulations allow us to address the question of whether there is a critical metallicity above which fine structure cooling from metals allows efficient fragmentation to occur, producing an initial mass function (IMF) resembling the local Salpeter IMF, rather than only high-mass stars.

The initial stages of planet formation in circumstellar gas discs proceed via dust grains that collide and build up larger and larger bodies (Safronov 1969). How this process continues from metre-sized boulders to kilometre-scale planetesimals is a major unsolved problem (Dominik et al. 2007): boulders stick together poorly (Benz 2000), and spiral into the protostar in a few hundred orbits due to a head wind from the slower rotating gas (Weidenschilling 1977). Gravitational collapse of the solid component has been suggested to overcome this barrier (Safronov 1969, Goldreich & Ward 1973, Youdin & Shu 2002). Even low levels of turbulence, however, inhibit sedimentation of solids to a sufficiently dense midplane layer (Weidenschilling & Cuzzi 1993, Dominik et al. 2007), but turbulence must be present to explain observed gas accretion in protostellar discs (Hartmann 1998). Here we report the discovery of efficient gravitational collapse of boulders in locally overdense regions in the midplane. The boulders concentrate initially in transient high pressures in the turbulent gas (Johansen, Klahr, & Henning 2006), and these concentrations are augmented a further order of magnitude by a streaming instability (Youdin & Goodman 2005, Johansen, Henning, & Klahr 2006, Johansen & Youdin 2007) driven by the relative flow of gas and solids. We find that gravitationally bound clusters form with masses comparable to dwarf planets and containing a distribution of boulder sizes. Gravitational collapse happens much faster than radial drift, offering a possible path to planetesimal formation in accreting circumstellar discs.

This document contains refereed supplementary information for the paper ``Rapid planetesimal formation in turbulent circumstellar discs''. It contains 15 sections (\S1.1 -- \S1.15) that address a number of subjects related to the main paper. We describe in detail the Poisson solver used to find the self-potential of the solid particles, including a linear and a non-linear test problem (\S1.3). Dissipative collisions remove energy from the motion of the particles by collisional cooling (\S1.4), an effect that allows gravitational collapse to occur in somewhat less massive discs (\S1.7). A resolution study of the gravitational collapse of the boulders is presented in \S1.6. We find that gravitational collapse can occur in progressively less massive discs as the grid resolution is increased, likely due to the decreased smoothing of the particle-mesh self-gravity solver with increasing resolution. In \S1.10 we show that it is in good agreement with the Goldreich & Ward (1973) stability analysis to form several-hundred-km-sized bodies, when the analysis is applied to 5 AU and to regions of increased boulder column density. \S11 is devoted to the measurement of random speeds and collision speeds between boulders. We find good agreement between our measurements and analytical theory for the random speeds, but the measured collision speeds are 3 times lower than expected from analytical theory. Higher resolution studies, and an improved analytical theory of collision speeds that takes into account epicyclic motion, will be needed to determine whether collision speeds have converged. In \S1.12 we present models with no magnetic fields. The boulder layer still exhibits strong clumping, due to the streaming instability, if the global solids-to-gas ratio is increased by a factor 3. Gravitational collapse occurs as readily as in magnetised discs.

Large-scale star formation in disk galaxies is hypothesized to be driven by global gravitational instability. The observed gas surface density is commonly used to compute the strength of gravitational instability, but according to this criterion star formation often appears to occur in gravitationally stable regions. One possible reason is that the stellar contribution to the instability has been neglected. We have examined the gravitational instability of the Large Magellanic Cloud (LMC) considering the gas alone, and considering the combination of collisional gas and collisionless stars. We compare the gravitationally unstable regions with the on-going star formation revealed by Spitzer observations of young stellar objects. Although only 62% of the massive young stellar object candidates are in regions where the gas alone is unstable, some 85% lie in regions unstable due to the combination of gas and stars. The combined stability analysis better describes where star formation occurs. In agreement with other observations and numerical models, a small fraction of the star formation occurs in regions with gravitational stability parameter Q > 1. We further measure the dependence of the star formation timescale on the strength of gravitational instability, and quantitatively compare it to the exponential dependence expected from numerical simulations.

Migration of protoplanets in their gaseous host disks may be largely responsible for the observed orbital distribution of extrasolar planets. Recent simulations have shown that the magnetorotational turbulence thought to drive accretion in protoplanetary disks can affect migration by turning it into an orbital random walk. However, these simulations neglected the disk's ionization structure. Low ionization fraction near the midplane of the disk can decouple the magnetic field from the gas, forming a dead zone with reduced or no turbulence. Here, to understand the effect of dead zones on protoplanetary migration, we perform numerical simulations of a small region of a stratified disk with magnetorotational turbulence confined to thin active layers above and below the midplane. Turbulence in the active layers exerts decreased, but still measurable, gravitational torques on a protoplanet located at the disk midplane. We find a decrease of two orders of magnitude in the diffusion coefficient for dead zones with dead-to-active surface density ratios approaching realistic values in protoplanetary disks. This torque arises primarily from density fluctuations within a distance of one scale height of the protoplanet. Turbulent torques have correlation times of only $\sim 0.3$ orbital periods and apparently time-stationary distributions. These properties are encouraging signs that stochastic methods can be used to determine the orbital evolution of populations of protoplanets under turbulent migration. Our results indicate that dead zones may be dynamically distinct regions for protoplanetary migration.

The nature of ultracompact H II regions (UCHRs) remains poorly determined. In particular, they are about an order of magnitude more common than would be expected if they formed around young massive stars and lasted for one dynamical time, around 10^4 yr. We here perform three-dimensional numerical simulations of the expansion of an H II region into self-gravitating, radiatively cooled gas, both with and without supersonic turbulent flows. In the laminar case, we find that H II region expansion in a collapsing core produces nearly spherical shells, even if the ionizing source is off-center in the core. This agrees with analytic models of blast waves in power-law media. In the turbulent case, we find that the H II region does not disrupt the central collapsing region, but rather sweeps up a shell of gas in which further collapse occurs. Although this does not constitute triggering, as the swept-up gas would eventually have collapsed anyway, it does expose the collapsing regions to ionizing radiation. We suggest that these regions of secondary collapse, which will not all themselves form massive stars, may form the bulk of observed UCHRs. As the larger shell will take over 10^5 years to complete its evolution, this could solve the timescale problem. Our suggestion is supported by the ubiquitous observation of more diffuse emission surrounding UCHRs.

Recent observations by Ferrarese et al. (2006) and Wehner et al. (2006) reveal that a majority of galaxies contain a central massive object (CMO), either a supermassive black hole (SMBH) or a compact stellar nucleus, regardless of the galaxy mass or morphological type, and that there is a tight relation between the masses of CMOs and those of the host galaxies. Several recent studies show that feedback from black holes can successfully explain the $\msigma$ correlation in massive elliptical galaxies that contain SMBHs. However, puzzles remain in spirals or dwarf spheroids that do not appear to have black holes but instead harbor a compact central stellar cluster. Here we use three-dimensional, smoothed particle hydrodynamics simulations of isolated galaxies to study the formation and evolution of CMOs in bulgeless disk galaxies, and simulations of merging galaxies to study the transition of the CMO--host mass relation from late-type bulgeless spirals to early-type ellipticals. Our results suggest that the observed correlations may be established primarily by the depletion of gas in the central region by accretion and star-formation, and may hold for all galaxy types. A systematic search for CMOs in the nuclei of bulgeless disk galaxies would offer a test of this conclusion. (Abridged)

We study the influence of low levels of metal enrichment on the cooling and collapse of ionized gas in small protogalactic halos using three-dimensional, smoothed particle hydrodynamics simulations. Our initial conditions represent protogalaxies forming within a fossil HII region -- a previously ionized HII region which has not yet had time to cool and recombine. Prior to cosmological reionization, such regions should be relatively common, since the characteristic lifetimes of the likely ionizing sources are significantly shorter than a Hubble time. We show that in these regions, H_2 is the dominant and most effective coolant, and that it is the amount of H_2 formed that determines whether or not the gas can collapse and form stars. At the low metallicities (Z < 10^{-3} Z_sun) thought to be associated with the transition from population III to early population II star formation, metal line cooling has an almost negligible effect on the evolution of low density gas, altering the density and temperature evolution of the gas by less than 1% compared to the metal-free case at densities below 1 cm^{-3} and temperatures above 2000 K. Although there is evidence that metal line cooling becomes more effective at higher density, we find no significant differences in behaviour from the metal-free case at any density below our sink particle creation threshold at n = 500 cm^{-3}. Increasing the metallicity also increases the importance of metal line cooling, but it does not significantly affect the dynamical evolution of the low density gas until Z = 0.1 Z_sun. This result holds regardless of whether or not an ultraviolet background is present.

We study the formation of H2 in the ISM, using a modified version of the astrophysical magnetohydrodynamical code ZEUS-MP that includes a non-equilibrium treatment of the formation and destruction of H2. We examine two different approximations to treat the shielding of H2 against photodissociation: a local approximation, which gives us a solid lower bound on the amount of shielding, and a method based on ray-tracing that is considerably more accurate in some circumstances but that produces results that are harder to clearly interpret. Either approximation allows one to perform three-dimensional high-resolution simulations of cloud formation with only modest computational resources. We also include a detailed treatment of the thermal behaviour of the gas. In this paper, we focus on the problem of molecular cloud formation in gravitationally unstable, initially static gas. We show that in these conditions, and for initial densities consistent with those observed in the cold, neutral atomic phase of the interstellar medium, H2 formation occurs on a timescale t > 10 Myr, comparable to or longer than the gravitational free-fall timescale of the cloud. We also show that the collapsing gas very quickly reaches thermal equilibrium and that the equation of state of the gas is generally softer than isothermal. Finally, we demonstrate that although these results show little sensitivity to variations in most of our simulation parameters, they are highly sensitive to the assumed initial density n_i. Reducing n_i significantly increases the cloud formation timescale and decreases the amount of hydrogen ultimately converted to H2. (Abridged).

(Abridged). In this paper, we present results from a large set of numerical simulations that demonstrate that H2 formation occurs rapidly in turbulent gas. Starting with purely atomic hydrogen, large quantities of molecular hydrogen can be produced on timescales of 1 -- 2 Myr, given turbulent velocity dispersions and magnetic field strengths consistent with observations. Moreover, as our simulations underestimate the effectiveness of H2 self-shielding and dust absorption, we can be confident that the molecular fractions that we compute are strong lower limits on the true values. The formation of large quantities of H2 on the timescale required by rapid cloud formation models therefore appears to be entirely plausible. We also investigate the density and temperature distributions of gas in our model clouds. We show that the density probability distribution function is approximately log-normal, with a dispersion that agrees well with the prediction of Padoan, Nordlund & Jones (1997). The temperature distribution is similar to that of a polytrope, with an effective polytropic index gamma_eff \simeq 0.8, although at low gas densities, the scatter of the actual gas temperature around this mean value is considerable, and the polytropic approximation does not capture the full range of behaviour of the gas.

This paper describes ZEUS-MP, a multi-physics, massively parallel, message- passing implementation of the ZEUS code. ZEUS-MP differs significantly from the ZEUS-2D code, the ZEUS-3D code, and an early "version 1" of ZEUS-MP distributed publicly in 1999. ZEUS-MP offers an MHD algorithm better suited for multidimensional flows than the ZEUS-2D module by virtue of modifications to the Method of Characteristics scheme first suggested by Hawley and Stone (1995), and is shown to compare quite favorably to the TVD scheme described by Ryu et. al (1998). ZEUS-MP is the first publicly-available ZEUS code to allow the advection of multiple chemical (or nuclear) species. Radiation hydrodynamic simulations are enabled via an implicit flux-limited radiation diffusion (FLD) module. The hydrodynamic, MHD, and FLD modules may be used in one, two, or three space dimensions. Self gravity may be included either through the assumption of a GM/r potential or a solution of Poisson's equation using one of three linear solver packages (conjugate-gradient, multigrid, and FFT) provided for that purpose. Point-mass potentials are also supported. Because ZEUS-MP is designed for simulations on parallel computing platforms, considerable attention is paid to the parallel performance characteristics of each module. Strong-scaling tests involving pure hydrodynamics (with and without self-gravity), MHD, and RHD are performed in which large problems (256^3 zones) are distributed among as many as 1024 processors of an IBM SP3. Parallel efficiency is a strong function of the amount of communication required between processors in a given algorithm, but all modules are shown to scale well on up to 1024 processors for the chosen fixed problem size.

We review the properties of turbulent molecular clouds (MCs), focusing on the physical processes that influence star formation (SF). MC formation appears to occur during large-scale compression of the diffuse ISM driven by supernovae, magnetorotational instability, or gravitational instability in galactic disks of stars and gas. The compressions generate turbulence that can accelerate molecule production and produce the observed morphology. We then review the properties of MC turbulence, including density enhancements observed as clumps and cores, magnetic field structure, driving scales, the relation to observed scaling relations, and the interaction with gas thermodynamics. We argue that MC cores are dynamical, not quasistatic, objects with relatively short lifetimes not exceeding a few megayears. We review their morphology, magnetic fields, density and velocity profiles, and virial budget. Next, we discuss how MC turbulence controls SF. On global scales turbulence prevents monolithic collapse of the clouds; on small scales it promotes local collapse. We discuss its effects on the SF efficiency, and critically examine the possible relation between the clump mass distribution and the initial mass function, and then turn to the redistribution of angular momentum during collapse and how it determines the multiplicity of stellar systems. Finally, we discuss the importance of dynamical interactions between protostars in dense clusters, and the effect of the ionization and winds from those protostars on the surrounding cloud. We conclude that the interaction of self-gravity and turbulence controls MC formation and behavior, as well as the core and star formation processes within them.

(Abridged). We model gravitational instability in a wide range of isolated disk galaxies, using GADGET, a three-dimensional, smoothed particle hydrodynamics code. The model galaxies include a dark matter halo and a disk of stars and isothermal gas. The global Schmidt law observed in disk galaxies is quantitatively reproduced by our models. We find that the surface density of star formation rate directly correlates with the strength of local gravitational instability. The local Schmidt laws of individual galaxies in our models show clear evidence of star formation thresholds. Our results suggest that the non-linear development of gravitational instability determines the local and global Schmidt laws, and the star formation thresholds.

A roughly constant temperature over a wide range of densities is maintained in molecular clouds through radiative heating and cooling. An isothermal equation of state is therefore frequently employed in molecular cloud simulations. However, the dynamical processes in molecular clouds include shock waves, expansion waves, cooling induced collapse and baroclinic vorticity, all incompatible with the assumption of a purely isothermal flow. Here, we incorporate an energy equation including all the important heating and cooling rates and a simple chemical network into simulations of three-dimensional, hydrodynamic, decaying turbulence. This allows us to test the accuracy of the isothermal assumption by directly comparing a model run with the modified energy equation to an isothermal model. We compute an extreme case in which the initial turbulence is sufficiently strong to dissociate much of the gas and alter the specific heat ratio. The molecules then reform as the turbulence weakens. We track the true specific heat ratio as well as its effective value. We analyse power spectra, vorticity and shock structures, and discuss scaling laws for decaying turbulence. We derive some limitations to the isothermal approximation for simulations of the interstellar medium using simple projection techniques. Overall, even given the extreme conditions, we find that an isothermal flow provides an adequate physical and observational description of many properties. The main exceptions revealed here concern behaviour directly related to the high temperature zones behind the shock waves.

To study how supernova feedback structures the turbulent interstellar medium, we construct 3D models of vertically stratified gas stirred by discrete supernova explosions, including vertical gravitational field and parametrized heating and cooling. The models reproduce many observed characteristics of the Galaxy such as global circulation of gas (i.e., galactic fountain) and the existence of cold dense clouds in the galactic disk. Global quantities of the model such as warm and hot gas filling factors in the midplane, mass fraction of thermally unstable gas, and the averaged vertical density profile are compared directly with existing observations, and shown to be broadly consistent. We find that energy injection occurs over a broad range of scales. There is no single effective driving scale, unlike the usual assumption for idealized models of incompressible turbulence. However, >90% of the total kinetic energy is contained in wavelengths shortward of 200 pc. The shape of the kinetic energy spectrum differs substantially from that of the velocity power spectrum, which implies that the velocity structure varies with the gas density. Velocity structure functions demonstrate that the phenomenological theory proposed by Boldyrev is applicable to the medium. We show that it can be misleading to predict physical properties such as the stellar initial mass function based on 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. Supernova-driven turbulence inhibits star formation globally rather than triggering it.

We investigate the question of whether ambipolar diffusion (ion-neutral drift) determines the smallest length and mass scale on which structure forms in a turbulent molecular cloud. We simulate magnetized turbulence in a mostly neutral, uniformly driven, turbulent medium, using a three-dimensional, two-fluid, magnetohydrodynamics (MHD) code modified from Zeus-MP. We find that substantial structure persists below the ambipolar diffusion scale because of the propagation of compressive slow MHD waves at smaller scales. Contrary to simple scaling arguments, ambipolar diffusion thus does not suppress structure below its characteristic dissipation scale as would be expected for a classical diffusive process. We have found this to be true for the magnetic energy, velocity, and density. Correspondingly, ambipolar diffusion leaves the clump mass spectrum unchanged. Ambipolar diffusion appears unable to set a characteristic scale for gravitational collapse and star formation in turbulent molecular clouds.

We model gravitational collapse leading to star formation in a wide range of isolated disk galaxies using a three-dimensional, smoothed particle hydrodynamics code. The model galaxies include a dark matter halo and a disk of stars and isothermal gas. Absorbing sink particles are used to directly measure the mass of gravitationally collapsing gas. They reach masses characteristic of stellar clusters. In this paper, we describe our galaxy models and numerical methods, followed by an investigation of the gravitational instability in these galaxies. Gravitational collapse forms star clusters with correlated positions and ages, as observed, for example, in the Large Magellanic Cloud. Gravitational instability alone acting in unperturbed galaxies appears sufficient to produce flocculent spiral arms, though not more organized patterns. Unstable galaxies show collapse in thin layers in the galactic plane; associated dust will form thin dust lanes in those galaxies, in agreement with observations. (abridged)

We study gravitational instability and consequent star formation in a wide range of isolated disk galaxies, using three-dimensional, smoothed particle hydrodynamics simulations at resolution sufficient to fully resolve gravitational collapse. Stellar feedback is represented by an isothermal equation of state. Absorbing sink particles are inserted in dynamically bound, converging regions with number density $n > 10^3$ cm$^{-3}$ to directly measure the mass of gravitationally collapsing gas available for star formation. Our models quantitatively reproduce not only the observed Schmidt law, but also the observed star formation threshold in disk galaxies. Our results suggest that the dominant physical mechanism determining the star formation rate is just the strength of gravitational instability, with feedback primarily functioning to maintain a roughly constant effective sound speed.

We present a high-resolution simulation of globular cluster formation in a galaxy merger. For the first time in such a simulation, individual star clusters are directly identified and followed on their orbits. We quantitatively compare star formation in the merger to that in the unperturbed galaxies. The merging galaxies show a strong starburst, in sharp contrast to their isolated progenitors. Most star clusters form in the tidal features. With a mass range of $5\times10^{5}$--$5\times 10^{6} M_{\odot}$, they are identified as globular clusters. The merger remnant is an elliptical galaxy. Clusters with different mass or age have different radial distributions in the galaxy. Our results show that the high specific frequency and bimodal distribution of metallicity observed in elliptical galaxies are natural products of gas-rich mergers, supporting a merger origin for the ellipticals and their globular cluster systems.

Observations have suggested substantial departures from pressure equilibrium in the interstellar medium (ISM) in the plane of the Galaxy, even on scales under 50 pc. Nevertheless, multi-phase models of the ISM assume at least locally isobaric gas. The pressure then determines the density reached by gas cooling to stable thermal equilibrium. We use numerical models of the magnetized ISM to examine the consequences of supernova driving for interstellar pressures. In this paper we examine a (200 pc)^3 periodic domain threaded by magnetic fields. Individual parcels of gas at different pressures reach widely varying points on the thermal equilibrium curve: no unique set of phases is found, but rather a dynamically-determined continuum of densities and temperatures. A substantial fraction of the gas remains entirely out of thermal equilibrium. Our results appear consistent with observations of interstellar pressures. They also suggest that the high pressures observed in molecular clouds may be due to ram pressures in addition to gravitational forces. Much of the gas in our model lies far from equipartition between thermal and magnetic pressures, with ratios ranging from 0.1 to $10^4$ and ratios of uniform to fluctuating magnetic field of 0.5--1. Our models show broad pressure probability distribution functions with log-normal functional forms produced by both shocks and rarefaction waves, rather than power-law distributions produced by isolated supernova remnants. The width of the distribution can be described quantitatively by a formula derived from the work of Padoan, Nordlund, & Jones (1997).

Identifying the processes that determine the initial mass function of stars (IMF) is a fundamental problem in star formation theory. One of the major uncertainties is the exact chemical state of the star forming gas and its influence on the dynamical evolution. Most simulations of star forming clusters use an isothermal equation of state (EOS). However, theoretical predictions and observations suggest that the effective polytropic exponent gamma in the EOS varies with density. We address these issues and study the effect of a piecewise polytropic EOS on the formation of stellar clusters in turbulent, self-gravitating molecular clouds using three-dimensional, smoothed particle hydrodynamics simulations. To approximate the results of published predictions of the thermal behavior of collapsing clouds, we increase the polytropic exponent gamma from 0.7 to 1.1 at some chosen density n_c, which we vary. The change of thermodynamic state at n_c selects a characteristic mass scale for fragmentation M_ch, which we relate to the peak of the observed IMF. Our investigation generally supports the idea that the distribution of stellar masses depends mainly on the thermodynamic state of the star-forming gas. The thermodynamic state of interstellar gas is a result of the balance between heating and cooling processes, which in turn are determined by fundamental atomic and molecular physics and by chemical abundances. Given the abundances, the derivation of a characteristic stellar mass can thus be based on universal quantities and constants.

Several topics in the theory of magnetized turbulence are reviewed with application to star formation and the interstellar medium. The density, pressure, and temperature distribution in a turbulent interstellar medium is described in comparison to a medium dominated by thermal instability. The derivation of the empirical theory for the energy spectrum of hydrodynamic and MHD turbulence is outlined, and comparisons are given to numerical models for the magnetized case. Next discussed is how density fluctuations modify observations of velocity fluctuations in line centroids, coincidentally cancelling the effects of projection smoothing. The analytic description of velocity structure functions as a function of the dimensionality of the dissipative structures is then covered. Finally the question of whether turbulence can support against gravitational collapse is addressed. Turbulence both prevents and promotes collapse, but its net effect is to inhibit collapse, so it does not trigger star formation in any significant sense.

(abridged version) Identifying the processes that determine the initial mass function of stars (IMF) is a fundamental problem in star formation theory. One of the major uncertainties is the exact chemical state of the star forming gas and its influence on the dynamical evolution. Most simulations of star forming clusters use an isothermal equation of state (EOS). We address these issues and study the effect of a piecewise polytropic EOS on the formation of stellar clusters in turbulent, self-gravitating molecular clouds. We increase the polytropic exponent gamma from 0.7 to 1.1 at some chosen density n_c, which we vary from from 4.3x10^4 cm^-3 to 4.3x10^7 cm^-3. The change of thermodynamic state at n_c selects a characteristic mass scale for fragmentation M_ch, which we relate to the peak of the observed IMF. We find a relation M_ch ~ n_c^-0.5, supporting the idea that the distribution of stellar masses largely depends on the thermodynamic state of the star-forming gas.

Several lines of evidence suggest that protogalactic and galactic environments manage to grow magnetic field very rapidly. This makes the theoretical problem of achieving rapid growth of magnetic fields in such environments difficult. Following Balsara, Benjamin and Cox (2001), we show that a multiphase ISM that is dominated by realistic energy inputs via supernova explosions will naturally become a strongly turbulent medium with large positive and negative values of the fluid helicity. The resultant interstellar turbulence seems to have some of the essential characteristics of a fast dynamo. Magnetic field growth in such a turbulent dynamo is studied in this paper. It is shown that the field undergoes a linear phase of evolution where it has a time of growth that is comparable to the eddy turn-over time. Energetics, power spectra, statistics and structures of the turbulent dynamo are studied here. It is shown that despite the fact that the ISM is driven by compressive motions, the kinetic energy in this high Mach number flow is mainly concentrated in solenoidal motions rather than compressive motions.

We study the effects of the initial conditions of turbulent molecular clouds on the ionization structure in newly formed H_{ii} regions, using three-dimensional, photon-conserving radiative transfer in a pre-computed density field from three-dimensional compressible turbulence. Our results show that the initial density structure of the gas cloud can play an important role in the resulting structure of the H_{ii} region. The propagation of the ionization fronts, the shape of the resulting H_{ii} region, and the total mass ionized depend on the properties of the turbulent density field. Cuts through the ionized regions generally show ``butterfly'' shapes rather than spherical ones, while emission measure maps are more spherical if the turbulence is driven on scales small compared to the size of the H_{ii} region. The ionization structure can be described by an effective clumping factor $\zeta=< n > \cdot < n^2>/< n>^2$, where $n$ is number density of the gas. The larger the value of $\zeta$, the less mass is ionized, and the more irregular the H_{ii} region shapes. Because we do not follow dynamics, our results apply only to the early stage of ionization when the speed of the ionization fronts remains much larger than the sound speed of the ionized gas, or Alfv\'en speed in magnetized clouds if it is larger, so that the dynamical effects can be negligible.

We model how repeated supernova explosions in high-redshift dwarf starburst galaxies drive superbubbles and winds out of the galaxies. We compute the efficiencies of metal and mass ejection and energy transport from the galactic potentials, including the effect of cosmological infall of external gas. The starburst bubbles quickly blow out of small, high-redshift, galactic disks, but must compete with the ram pressure of the infalling gas to escape into intergalactic space. We show that the assumed efficiency of the star formation rate dominates the bubble evolution and the metal, mass, and energy feedback efficiencies. With star formation efficiency f*=0.01, the ram pressure of infall can confine the bubbles around high-redshift dwarf galaxies with circular velocities v_c>52 km/s. We can expect high metal and mass ejection efficiencies, and moderate energy transport efficiencies in halos with v_c~30-50 km/s and f*~0.01 as well as in halos with v_c~100 km/s and f*>>0.01. Such haloes collapse successively from 1-2 sigma peaks in LambdaCDM Gaussian density perturbations as time progresses. These dwarf galaxies can probably enrich low and high-density regions of intergalactic space with metals to 10^-3-10^-2 Zsun as they collapse at z~8 and z<5 respectively. They also may be able to provide adequate turbulent energy to prevent the collapse of other nearby halos, as well as to significantly broaden Lyman-alpha absorption lines to v_rms~20-40 km/s. We compute the timescales for the next starbursts if gas freely falls back after a starburst, and find that, for star formation efficiencies as low as f*<0.01, the next starburst should occur in less than half the Hubble time at the collapse redshift. This suggests that episodic star formation may be ubiquitous in dwarf galaxies.

We study the escape of Lyman continuum radiation from the disks of dwarf starburst galaxies, with and without supershells, by solving the radiation transfer problem of stellar radiation through them. We model disks with Md=10^8-10Msun, with exponential surface density profiles as a function of redshift, and model the of repeated supernova explosions driving supershells out of the disks, using the hydrodynamic simulation code, ZEUS-3D. The amount of star formation is assumed proportional to mass above some density threshold. We vary the threshold to explore the range of star formation efficiencies, f*=0.006, 0.06, and 0.6. We find that the interstellar gas swept up in dense supershells can effectively trap the ionizing photons, before the supershells blow out of the disks. The blow-outs then create galactic outflows, chimneys which allow the photons to escape directly to the intergalactic medium. Our results are consistent with escape fractions of less than 0.1 measured in local dwarf starburst galaxies, because they are likely observed while the starbursts are young, before blow-out. We suggest that high-redshift dwarf starburst galaxies may make a substantial contribution to the UV background radiation with total escape fractions >0.2, as expected if star formation efficiencies >0.06.

We study the effect of varying the equation of state on the formation of stellar clusters in turbulent molecular clouds, using three-dimensional, smoothed particle hydrodynamics simulations. Our results show that the equation of state helps determine how strongly self-gravitating gas fragments. The degree of fragmentation decreases with increasing \new{polytropic exponent} $\gamma$ in the range $0.2 < \gamma < 1.4$, although the total amount of \new{mass accreted onto collapsed fragments} appears to remain roughly constant through that range. Low values of $\gamma$ are expected to lead to the formation of dense clusters of low-mass stars, while $\gamma>1$ probably results in the formation of isolated and massive stars. Fragmentation and collapse ceases entirely for $\gamma > 1.4$ as expected from analytic arguments. The mass spectrum of overdense gas clumps is roughly log-normal for {\em non}-self-gravitating turbulent gas, but changes to a power-law under the action of gravity. The spectrum of collapsed cores, on the other hand, remains log-normal for $\gamma\le 1$, but flattens markedly for $\gamma >1$. The density PDFs approach log-normal, with widths that decrease with increasing $\gamma$. Primordial gas may have effective $\gamma > 1$, in which case these results could help explain why models of the formation of the first stars tend to produce isolated, massive objects.

Stars form by gravoturbulent fragmentation of interstellar gas clouds. The supersonic turbulence ubiquitously observed in Galactic molecular gas generates strong density fluctuations with gravity taking over in the densest and most massive regions. Collapse sets in to build up stars and star clusters. Turbulence plays a dual role. On global scales it provides support, while at the same time it can promote local collapse. Stellar birth is thus intimately linked to the dynamical behavior of parental gas cloud, which determines when and where protostellar cores form, and how they contract and grow in mass via accretion from the surrounding cloud material to build up stars. Slow, inefficient, isolated star formation is a hallmark of turbulent support, whereas fast, efficient, clustered star formation occurs in its absence. The fact that Galactic molecular clouds are highly filamentary can be explained by a combination of compressional flows and shear. The dynamical evolution of nascent star clusters is very complex. This strongly influences the stellar mass spectrum. The equation of state (EOS) plays a pivotal role in the fragmentation process. Under typical cloud conditions, massive stars form as part of dense clusters. However, for gas with effective polytropic index greater than unity star formation becomes biased towards isolated massive stars, which may be of relevance for understanding Pop III stars.

We use ZEUS-MP to perform high resolution, three-dimensional, super-Alfvenic turbulent simulations in order to investigate the role of magnetic fields in self-gravitating core formation within turbulent molecular clouds. Statistical properties of our super-Alfvenic model without gravity agree with previous similar studies. Including self-gravity, our models give the following results. They are consistent with the turbulent fragmentation prediction of the core mass distribution of Padoan & Nordlund. They also confirm that local gravitational collapse is not prevented by magnetohydrodynamic waves driven by turbulent flows, even when the turbulent Jeans mass exceeds the mass in the simulation volume. Comparison of results between 256^3 and 512^3 zone simulations reveals convergence in the collapse rate. Analysis of self-gravitating cores formed in the simulation shows that: (1) All cores formed are magnetically supercritical by at least an order of magnitude. (2) A power law relation between central magnetic field strength and density B_c propto rho_c^{1/2} is observed despite the cores being strongly supercritical. (3) Specific angular momentum j propto R^{3/2} for cores with radius R. (4) Most cores are prolate and triaxial in shape, in agreement with the results of Gammie et al.

(Modified) The scaling of velocity fluctuation, dv, as a function of spatial scale L in molecular clouds can be measured from size-linewidth relations, principal component analysis, or line centroid variation. Differing values of the power law index of the scaling relation dv = L^(g3D) in 3D are given by these different methods: the first two give g3D=0.5, while line centroid analysis gives g3D=0. This discrepancy has previously not been fully appreciated, as the variation of projected velocity line centroid fluctuations (dv_{lc} = L^(g2D)) is indeed described, in 2D, by g2D=0.5. However, if projection smoothing is accounted for, this implies that g3D=0. We suggest that a resolution of this discrepancy can be achieved by accounting for the effect of density inhomogeneity on the observed g2D obtained from velocity line centroid analysis. Numerical simulations of compressible turbulence are used to show that the effect of density inhomogeneity statistically reverses the effect of projection smoothing in the case of driven turbulence so that velocity line centroid analysis does indeed predict that g2D=g3D=0.5. Using our numerical results we can restore consistency between line centroid analysis, principal component analysis and size-linewidth relations, and we derive g3D=0.5, corresponding to shock-dominated (Burgers) turbulence. We find that this consistency requires that molecular clouds are continually driven on large scales or are only recently formed.

How does turbulence contribute to the formation and structure of the dense interstellar medium (ISM)? Molecular clouds are dense, high-pressure objects. It is usually argued that gravitational confinement causes the high pressures, and that the clouds are (magneto)hydrostaticobjects supported by a balance between magnetic and turbulent pressures and gravity. However, magnetic pressures appear too weak, and MHD turbulent support not only requires driving, but also results in continuing gravitational collapse, as has now been demonstrated in simulations reaching 512^3 zones. Models of supernova-driven, magnetized turbulence readily form transient, high-pressure, dense regions that may form molecular clouds. They are contained not by self-gravity, but by turbulent ram pressures from the larger flow. Apparent virialization may actually be a geometrical effect. Turbulent clouds are unlikely to be in hydrostatic equilibrium, instead either collapsing or expanding, although they may appear well-fit by projected equilibrium Bonnor-Ebert spheres. Collapsing clouds probably form stars efficiently, while expanding ones can still form stars by turbulent compression, but rather inefficiently.

"The purpose of numerical models is not numbers but insight." (Hamming) In the spirit of this adage, and of Don Cox's approach to scientific speaking, we discuss the questions that the latest generation of numerical models of the interstellar medium raise, at least for us. The energy source for the interstellar turbulence is still under discussion. We review the argument for supernovae dominating in star forming regions. Magnetorotational instability has been suggested as a way of coupling disk shear to the turbulent flow. Models make evident that the unstable wavelengths are very long compared to thermally unstable wavelengths, with implications for star formation in the outer galaxy and low surface brightness disks. The perennial question of the factors determining the hot gas filling factor in a SN-driven medium remains open, in particular because of the unexpectedly strong turbulent mixing at the boundaries of hot cavities seen in the models. The formation of molecular clouds in the turbulent flow is also poorly understood. Dense regions suitable for cloud formation clearly form even in the absence of self-gravity, although their ultimate evolution remains to be computed.

Understanding the formation of stars in galaxies is central to much of modern astrophysics. For several decades it has been thought that stellar birth is primarily controlled by the interplay between gravity and magnetostatic support, modulated by ambipolar diffusion. Recently, however, both observational and numerical work has begun to suggest that support by supersonic turbulence rather than magnetic fields controls star formation. In this review we outline a new theory of star formation relying on the control by turbulence. We demonstrate that although supersonic turbulence can provide global support, it nevertheless produces density enhancements that allow local collapse. Inefficient, isolated star formation is a hallmark of turbulent support, while efficient, clustered star formation occurs in its absence. The consequences of this theory are then explored for both local star formation and galactic scale star formation. (ABSTRACT ABBREVIATED)

We study magnetic field steepening due to ambipolar diffusion (Brandenburg & Zweibel 1994) in protoplanetary disk environments and draw the following conclusions. Current sheets are generated in magnetically active regions of the disk where the ionization fraction is high enough for the magnetorotational instability to operate. In late stages of solar nebula evolution, the surface density is expected to have lowered and dust grains to have gravitationally settled to the midplane. If the local dust-to-gas mass ratio near the midplane is increased above cosmic abundances by factors > 10^3, current sheets reach high enough temperatures to melt millimeter-sized dust grains, and hence may provide the mechanism to form meteoritic chondrules. In addition, these current sheets possibly explain the near-infrared excesses observed in spectral energy distributions (SEDs) of young stellar objects. Direct imaging of protoplanetary disks via a nulling interferometer or, in the future, a multi-band, adaptive optics coronagraph can test this hypothesis.

Interstellar turbulence is expected to dissipate quickly in the absence of continuous energy input. I examine the energy available for driving the turbulence from six likely mechanisms: magnetorotational instability, gravitational instability, protostellar outflows, H II region expansion, stellar winds, and supernovae. I conclude that supernovae contribute far more energy than the other mechanisms, and so form the most likely drving mechanism.

We study numerically the effect of varying the equation of state (EOS) on the formation of stellar clusters in turbulent molecular clouds. Our results show that the EOS plays an important role in the fragmentation of the clouds, and the determination of the initial mass function (IMF) of the protostellar cores. Fragmentation decreases with increasing adiabatic index $\gamma$ in the range $0.2 < \gamma < 1.4$, although the total amount of mass accreted appears to remain roughly constant through that range, resulting in more massive cores at higher $\gamma$. Fragmentation and collapse ceases entirely for $\gamma > 1.4$ as expected from analytic arguments. Primordial gas may have effective $\gamma > 1$, in which case these results may help explain why models of the formation of the first stars produce isolated, massive objects.

MHD turbulence plays a central role in the physics of star-forming molecular clouds and the interstellar medium. I here show that MHD turbulence in molecular clouds must be driven to account for the observed supersonic motions in the clouds, as even strongly magnetized turbulence decays quickly. I then show that driven MHD turbulence can globally support gravitationally unstable regions, but local collapse inevitably occurs. Differences in the strength of driving and the gas density may explain the very different rates of star formation observed in different galaxies. The driving mechanism for the turbulence is likely a combination of field supernovae in star-forming sections of galactic disks, and magnetorotational instabilities in outer disks and low surface brightness galaxies. Supernova-driven turbulence has a broad range of pressures with a roughly log-normal distribution. High-pressure, cold regions can be formed even in the absence of self-gravity.

We study the mixing of chemical species in the interstellar medium (ISM). Recent observations suggest that the distribution of species such as deuterium in the ISM may be far from homogeneous. This raises the question of how long it takes for inhomogeneities to be erased in the ISM, and how this depends on the length scale of the inhomogeneities. We added a tracer field to the three-dimensional, supernova-driven ISM model of Avillez (2000) to study mixing and dispersal in kiloparsec-scale simulations of the ISM with different supernova (SN) rates and different inhomogeneity length scales. We find several surprising results. Classical mixing length theory fails to predict the very weak dependence of mixing time on length scale that we find on scales of 25--500 pc. Derived diffusion coefficients increase exponentially with time, rather than remaining constant. The variance of composition declines exponentially, with a time constant of tens of Myr, so that large differences fade faster than small ones. The time constant depends on the inverse square root of the supernova rate. One major reason for these results is that even with numerical diffusion exceeding physical values, gas does not mix quickly between hot and cold regions.

We present the results from three dimensional hydrodynamical simulations of decaying high-speed turbulence in dense molecular clouds. We compare our results, which include a detailed cooling function, molecular hydrogen chemistry and a limited C and O chemistry, to those previously obtained for decaying isothermal turbulence. After an initial phase of shock formation, power-law decay regimes are uncovered, as in the isothermal case. We find that the turbulence decays faster than in the isothermal case because the average Mach number remains higher, due to the radiative cooling. The total thermal energy, initially raised by the introduction of turbulence, decays only a little slower than the kinetic energy. We discover that molecule reformation, as the fast turbulence decays, is several times faster than that predicted for a non-turbulent medium. This is caused by moderate speed shocks which sweep through a large fraction of the volume, compressing the gas and dust. Through reformation, the molecular density and molecular column appear as complex patterns of filaments, clumps and some diffuse structure. In contrast, the molecular fraction has a wider distribution of highly distorted clumps and copious diffuse structure, so that density and molecular density are almost identically distributed during the reformation phase. We conclude that molecules form in swept-up clumps but effectively mix throughout via subsequent expansions and compressions.

We compare velocity structure in the Polaris Flare molecular cloud at scales ranging from 0.015 pc to 20 pc to simulations of supersonic hydrodynamic and MHD turbulence computed with the ZEUS MHD code. We use several different statistical methods to compare models and observations. The Delta-variance wavelet transform is most sensitive to characteristic scales and scaling laws, but is limited by a lack of intensity weighting. The scanning-beam size-linewidth relation is more robust with respect to noisy data. Obtaining the global velocity scaling behaviour requires that large-scale trends in the maps not be removed but treated as part of the turbulent cascade. We compare the true velocity PDF in our models to velocity centroids and average line profiles in optically thin lines, and find that the line profiles reflect the true PDF better unless the map size is comparable to the total line-of-sight thickness of the cloud. Comparison of line profiles to velocity centroid PDFs can thus be used to measure the line-of-sight depth of a cloud. The observed density and velocity structure is consistent with supersonic turbulence with a driving scale at or above the size of the molecular cloud and dissipative processes below 0.05 pc. Ambipolar diffusion could explain the dissipation. The velocity PDFs exclude small-scale driving such as that from stellar outflows as a dominant process in the observed region. In the models, large-scale driving is the only process that produces deviations from a Gaussian PDF shape consistent with observations. Strong magnetic fields impose a clear anisotropy on the velocity field, reducing the velocity variance in directions perpendicular to the field. (abridged)

We examine how well the physical properties of clumps in turbulent molecular clouds can be determined by measurements of observed clump structures. We compare simulated observations of three-dimensional numerical models of isothermal, magnetized, supersonic turbulence to the actual physical structure of the models. We determine how changing the parameters of turbulence changes the structure of the simulations. Stronger driving produces greater density fluctuations, and longer wavelength driving produces larger structures. Magnetic fields have a less pronounced effect on structure, and one that is not monotonic with field strength. Aligned structures are seen only with low-density tracers, and when the intensity of the field is large. Comparing different regions with the same tracers (or the same region with different tracers) can give information about physical conditions: different density tracers can help determine the size of the density fluctuations and thus the strength of the driving. Nevertheless, velocity superposition of multiple physical clumps can fully obscure the physical properties of those clumps, and short wavelength driving worsens this effect. Regarding Larson's relationships, we confirm previous claims that the mean density-size relationship is an artifact of the observations; and the velocity dispersion-size relationship, is consistent with observations. Regarding the mass spectra, we show that, when we look for clumps with high enough resolution, they converge to a log-normal function, rather than the power-law obtained in the literature.

In order to understand the time scales for complete mixing in a supernova-driven ISM we used the model of Avillez (2000) to follow the mixing and dispersal of inhomogeneities with different length scales. We included tracer fields to follow the inhomogeneities in kpc-scale simulations with different supernova rates. We find that an initially inhomogeneous ISM still shows local variations and a clumpy distribution on time scales of several hundred Myr. This clumpy structure does not disapear over a time of a few tens of Myr, even if the rate of SNe is increased to ten times the Galactic rate. These local variations occur because, even in the presence of numerical diffusion, gas does not mix quickly between hot and cold regions. The lower limit that we are able to place on the mixing time at the kpc scale is already longer than the time for chemistry to reach local equilibrium in the gas. The simulations also show that the time scale to erase inhomogeneities in the ISM is nearly independent of their length scale over the range 25-500 pc, contrary to what is expected from classical mixing length theory.

The one-point statistics of column density distributions of turbulent molecular cloud models are investigated and compared with observations. In agreement with the observations, the number N of pixels with surface density S is distributed exponentially N(S)=exp(-S/S0) in models of driven compressible supersonic turbulence. However, in contrast to the observations, the exponential slope defined by S0 is not universal but instead depends strongly on the adopted rms Mach number and on the smoothing of the data cube. We demonstrate that this problem can be solved if one restricts the analysis of the surface density distribution to subregions with sizes equal to the correlation length of the flow which turns out to be given by the driving scale. In this case, the column density distributions are universal with a slope that is in excellent agreement with the observations and independent of the Mach number or smoothing. The observed molecular clouds therefore are coherent structures with sizes of order their correlation lengths. Turbulence inside these clouds must be driven on the largest scales, if at all. Numerical models of turbulent molecular clouds have to be restricted to cubes with sizes similar to the correlation lengthscale in order to be compared with observations. Our results imply that turbulence is generated on scales that are much larger that the Jeans length. In this case, gravitational collapse cannot be suppressed and star formation should start in molecular clouds within a dynamical timescale after their formation.

Observations have suggested substantial departures from pressure equilibrium in the interstellar medium (ISM) in the plane of the Galaxy, even on scales under 50 pc. Nevertheless, multi-phase models of the ISM assume at least locally isobaric gas. The pressure then determines the density reached by gas cooling to stable thermal equilibrium. We use two different sets of numerical models of the ISM to examine the consequences of supernova driving for interstellar pressures. The first set of models is hydrodynamical, and uses adaptive mesh refinement to allow computation of a 1 x 1 x 20 kpc section of a stratified galactic disk. The second set of models is magnetohydrodynamical, using an independent code framework, and examines a 200 pc cubed periodic domain threaded by magnetic fields. Both of these models show broad pressure distributions with roughly log-normal functional forms produced by both shocks and rarefaction waves, rather than the power-law distributions predicted by previous work, with rather sharp thermal pressure gradients. The width of the distribution of the logs of pressure in gas with log T < 3.9 is proportional to the rms Mach number in that gas, while the distribution in hotter gas is broader, but not so broad as would be predicted by the Mach numbers in that gas. Individual parcels of gas reach widely varying points on the thermal equilibrium curve: no unique set of phases is found, but rather a dynamically-determined continuum of densities and temperatures. Furthermore, a substantial fraction of the gas remains entirely out of thermal equilibrium. Our results appear consistent with observations of interstellar pressures, and suggest that the pressures observed in molecular clouds may be due to ram pressure rather than gravitational confinement.

The dynamical state of star-forming molecular clouds cannot be understood without determining the structure and strength of their magnetic fields. Measurements of polarized far-infrared radiation from thermally aligned dust grains are used to map the orientation of the field and estimate its strength, but the accuracy of the results has remained in doubt. In order to assess the reliability of this method, we apply it to simulated far-infrared polarization maps derived from three-dimensional simulations of supersonic magnetohydrodynamical turbulence, and compare the estimated values to the known magnetic field strengths in the simulations. We investigate the effects of limited telescope resolution and self-gravity on the structure of the maps. Limited observational resolution affects the field structure such that small scale variations can be completely suppressed, thus giving the impression of a very homogeneous field. The Chandrasekhar-Fermi method of estimating the mean magnetic field in a turbulent medium is tested, and we suggest an extension to measure the rms field. Both methods yield results within a factor of 2 for field strengths typical of molecular clouds, with the modified version returning more reliable estimates for slightly weaker fields. However, neither method alone works well for very weak fields, missing them by a factor of up to 150. Taking the geometric mean of both methods estimates even the weakest fields accurately within a factor of 2.5. Limited telescope resolution leads to a systematic overestimation of the field strengths for all methods. We discuss the effects responsible for this overestimation and show how to extract information on the underlying (turbulent) power spectrum.

Recent HI emission observations of the Southern Galactic hemisphere have revealed a mushroom-like structure extending from z=-70 to -450 pc, composed of a stem and a cap. Similar structures occur in three-dimensional simulations of a dynamic galactic disk driven by isolated and clustered supernovae. Using these simulations, we show that hot gas in the Galactic disk that is not evacuated through chimneys expands into the cooler gas of the thick disk, forming mushroom-shaped structures. This new class of objects traces buoyant flow of hot gas into the thick disk.

I argue that star formation is controlled by supersonic turbulence, drawing for support on a number of 3D hydrodynamical and MHD simulations as well as theoretical arguments. Clustered star formation appears to be a natural result of a lack of turbulent support, while isolated star formation is a signpost of global turbulent support.

Hydrodynamic supersonic turbulence can only prevent local gravitational collapse if the turbulence is driven on scales smaller than the local Jeans lengths in the densest regions, a very severe requirement (Paper I). Magnetic fields have been suggested to support molecular clouds either magnetostatically or via magnetohydrodynamic (MHD) waves. Whereas the first mechanism would form sheet-like clouds, the second mechanism not only could exert a pressure onto the gas counteracting the gravitational forces, but could lead to a transfer of turbulent kinetic energy down to smaller spatial scales via MHD wave interactions. This turbulent magnetic cascade might provide sufficient energy at small scales to halt local collapse. We test this hypothesis with MHD simulations at resolutions up to 256^3 zones, done with ZEUS-3D. We first derive a resolution criterion for self-gravitating, magnetized gas: in order to prevent collapse of magnetostatically supported regions due to numerical diffusion, the minimum Jeans length must be resolved by four zones. Resolution of MHD waves increases this requirement to roughly six zones. We then find that magnetic fields cannot prevent local collapse unless they provide magnetostatic support. Weaker magnetic fields do somewhat delay collapse and cause it to occur more uniformly across the supported region in comparison to the hydrodynamical case. However, they still cannot prevent local collapse for much longer than a global free-fall time.

Supersonic turbulence generates distributions of shock waves. Here, we analyse the shock waves in three-dimensional numerical simulations of uniformly driven supersonic turbulence, with and without magnetohydrodynamics and self-gravity. We can identify the nature of the turbulence by measuring the distribution of the shock strengths. We find that uniformly driven turbulence possesses a power law distribution of fast shocks with the number of shocks inversely proportional to the square root of the shock jump speed. A tail of high speed shocks steeper than Gaussian results from the random superposition of driving waves which decay rapidly. The energy is dissipated by a small range of fast shocks. These results contrast with the exponential distribution and slow shock dissipation associated with decaying turbulence. A strong magnetic field enhances the shock number transverse to the field direction at the expense of parallel shocks. A simulation with self-gravity demonstrates the development of a number of highly dissipative accretion shocks. Finally, we examine the dynamics to demonstrate how the power-law behaviour arises.

Numerical models of the dynamical interstellar medium show that interactions between structures such as supernova remnants and superbubbles are more important than the structures themselves in determining the behavior of the ISM. I review the techniques and conclusions of recent global models, focussing on what they tell us about the formation of star-forming giant molecular clouds.

We here analyse numerical simulations of supersonic, hypersonic and magnetohydrodynamic turbulence that is free to decay. Our goals are to understand the dynamics of the decay and the characteristic properties of the shock waves produced. This will be useful for interpretation of observations of both motions in molecular clouds and sources of non-thermal radiation. We find that decaying hypersonic turbulence possesses an exponential tail of fast shocks and an exponential decay in time, i.e. the number of shocks is proportional to t exp (-ktv) for shock velocity jump v and mean initial wavenumber k. In contrast to the velocity gradients, the velocity Probability Distribution Function remains Gaussian with a more complex decay law. The energy is dissipated not by fast shocks but by a large number of low Mach number shocks. The power loss peaks near a low-speed turn-over in an exponential distribution. An analytical extension of the mapping closure technique is able to predict the basic decay features. Our analytic description of the distribution of shock strengths should prove useful for direct modeling of observable emission. We note that an exponential distribution of shocks such as we find will, in general, generate very low excitation shock signatures.

Observed molecular clouds often appear to have very low star formation efficiencies and lifetimes an order of magnitude longer than their free-fall times. Their support is attributed to the random supersonic motions observed in them. We study the support of molecular clouds against gravitational collapse by supersonic, gas dynamical turbulence using direct numerical simulation. Computations with two different algorithms are compared: a particle-based, Lagrangian method (SPH), and a grid-based, Eulerian, second-order method (ZEUS). The effects of both algorithm and resolution can be studied with this method. We find that, under typical molecular cloud conditions, global collapse can indeed be prevented, but density enhancements caused by strong shocks nevertheless become gravitationally unstable and collapse into dense cores and, presumably, stars. The occurance and efficiency of local collapse decreases as the driving wave length decreases and the driving strength increases. It appears that local collapse can only be prevented entirely with unrealistically short wave length driving, but observed core formation rates can be reproduced with more realistic driving. At high collapse rates, cores are formed on short time scales in coherent structures with high efficiency, while at low collapse rates they are scattered randomly throughout the region and exhibit considerable age spread. We suggest that this naturally explains the observed distinction between isolated and clustered star formation.

I review what is known about the temperature of the plasma within stellar wind bubbles and superbubbles. Classical theory suggests that it should be hot, with characteristic temperatures of order a million degrees. This temperature should be set by the balance between heating by the internal termination shocks of the central stellar winds and supernovae, which expand at thousands of km/s, and cooling by conductive evaporation of cold gas off the shell walls. However, if the hot interior gas becomes dense enough due to evaporation or ablation off of interior clouds, it will cool in less than a dynamical time, leading to a cold interior. The observational evidence appears mixed. On the one hand, X-ray emission has been observed from both stellar wind bubbles and superbubbles. On the other hand, no stellar wind bubble or superbubble has yet been observed emitting at the rate predicted by the classical theory: they are either too faint or too bright, by up to an order of magnitude. Alternate explanations have been proposed for the observed emission, including off-center supernova remnants hitting the shell walls of superbubbles, and residual emission from highly-ionized gas out of coronal equilibrium. Furthermore, the structures of post-main sequence stellar wind bubbles, expanding into what are presumably old stellar wind bubbles, appear in at least some cases to show that the bubble interior is cold, not hot. (The classical example of this is NGC 6888.) What is the actual state of bubble and superbubble interiors?

We review how supersonic turbulence can both prevent and promote the collapse of molecular clouds into stars. First we show that decaying turbulence cannot significantly delay collapse under conditions typical of molecular clouds, regardless of magnetic field strength so long as the fields are not supporting the cloud magnetohydrostatically. Then we review possible drivers and examine simulations of driven supersonic and trans Alfv\'enic turbulence, finally including the effects of self-gravity. Our preliminary results show that, although turbulence can support regions against gravitational collapse, the strong compressions associated with the required velocities will tend to promote collapse of local condensations.

Modeling the structure of molecular clouds depends on good methods to statistically compare simulations with observations in order to constrain the models. Here we characterize a suite of hydrodynamical and magnetohydrodynamical (MHD) simulations of supersonic turbulence using an averaged wavelet transform, the Delta-variance, that has been successfully used to characterize observations. We find that, independent of numerical resolution and dissipation, the only models that produce scale-free, power-law Delta-variance spectra are those with hypersonic Mach numbers above M ~ 4, while slower supersonic turbulence tends to show characteristic scales and produce non-power-law spectra. Magnetic fields have only a minor influence on this tendency, though they tend to reduce the scale-free nature of the turbulence, and increase the transfer of energy from large to small scales. The evolution of the characteristic length scale seen in supersonic turbulence follows exactly the t^{1/2} power-law predicted from recent studies of the kinetic energy decay rate.

The observed giant outbursts of Luminous Blue Variables (LBVs) may occur when these massive stars approach their Eddington limits. When this happens, they must reach a point where the centrifugal force and the radiative acceleration cancel out gravity at the equator. We call this the Omega-limit. When stars are close to the Omega-limit, strong non-spherical mass loss should occur. This suggests a scenario where a slow and very dense wind, strongly confined to the equatorial plane, is followed by a fast and almost spherical wind. We compute two-dimensional hydrodynamic models of the evolution of the nebula formed from such interacting winds, using parameters consistent with the outburst of eta Carinae in the last century. This outburst gave birth to the Homunculus, the hourglass-shaped inner part of a highly structured circumstellar nebula. Assuming the star was very close to the Omega-limit during outburst, our models produce gas distributions that strongly resemble the Homunculus on large and small scale. This supports the general conjecture that giant outbursts in LBVs occur when they approach the Eddington limit. Our models constrains the average mass loss rate since the outburst to values smaller than the present-day mass loss rate and suggest that eta Car is approaching another outburst. Our models imply that the occurrence of giant LBV outbursts depends on the initial stellar rotation rate, and that the initial angular momentum is as important to the evolution of very massive stars as their initial mass or metallicity.

In this review I focus on how massive stars in the disks of galaxies determine the properties of the ionized, gaseous haloes of star-forming galaxies. I first examine the ionization state of the halo, and conclude that diluted radiation from massive stars, along with a possible secondary heating source such as shocks, appears sufficient to explain the observations. I then look at how single superbubbles evolve, as they must provide one of the main mechanisms for transfer of mass and energy to the halo, as well as opening up lines of sight for ionizing radiation to escape from the disk. Finally, I look at the interactions of superbubbles, as they will rarely evolve in isolation, and describe recent hydrodynamical and MHD simulations of the interstellar medium shaped by multiple superbubbles.

Molecular clouds have broad linewidths suggesting turbulent supersonic motions in the clouds. These motions are usually invoked to explain why molecular clouds take much longer than a free-fall time to form stars. It has classically been thought that supersonic hydrodynamical turbulence would dissipate its energy quickly, but that the introduction of strong magnetic fields could maintain these motions. In a previous paper it has been shown, however, that isothermal, compressible, MHD and hydrodynamical turbulence decay at virtually the same rate, requiring that constant driving occur to maintain the observed turbulence. In this paper direct numerical computations of uniformly driven turbulence with the ZEUS astrophysical MHD code are used to derive the absolute value of energy dissipation as a function of the driving wavelength and amplitude. The ratio of the formal decay time of turbulence E_{kin}/\dot{E}_{kin} to the free-fall time of the gas can then be derived as a function of the ratio of driving wavelength to Jeans wavelength and rms Mach number, and shown to be most likely far less than unity, again showing that turbulence in molecular clouds must be constantly and strongly driven. (abridged)

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 eta Carinae.

We model the effects of repeated supernova explosions from starbursts in dwarf galaxies on the interstellar medium of these galaxies, taking into account the gravitational potential of their dominant dark matter haloes. We explore supernova rates from one every 30,000 yr to one every 3 million yr, equivalent to steady mechanical luminosities of L=0.1-10 x 10^38 ergs/s, occurring in dwarf galaxies with gas masses M_g=10^6 - 10^9 solar masses. We address in detail, both analytically and numerically, the following three questions: 1. When do the supernova ejecta blow out of the disk of the galaxy? 2. When blowout occurs, what fraction of the interstellar gas is blown away, escaping the potential of the galactic halo? 3. What happens to the metals ejected from the massive stars of the starburst? Are they retained or blown away? We give quantitative results for when blowout will or will not occur in galaxies with 10^6 \leq M_g \leq 10^9 solar masses. Surprisingly, we find that the mass ejection efficiency is very low for galaxies with mass M_g \geq 10^7 solar masses. Only galaxies with M_g \leq 10^6 solar masses have their interstellar gas blown away, and then virtually independently of L. On the other hand, metals from the supernova ejecta are accelerated to velocities larger than the escape speed from the galaxy far more easily than the gas. We find that for L_38=1, only about 30% of the metals are retained by a 10^9 solar mass galaxy, and virtually none by smaller galaxies. We discuss the implications of our results for the evolution, metallicity and observational properties of dwarf galaxies.

We present numerical studies of compressible, decaying turbulence, with and without magnetic fields, with initial rms Alfven and Mach numbers ranging up to five, and apply the results to the question of the support of star-forming interstellar clouds of molecular gas. We find that, in 1D, magnetized turbulence actually decays faster than unmagnetized turbulence. In all the regimes that we have studied 3D turbulence-super-Alfvenic, supersonic, sub-Alfvenic, and subsonic-the kinetic energy decays as (t-t0)^(-x), with 0.85 < x < 1.2. We compared results from two entirely different algorithms in the unmagnetized case, and have performed extensive resolution studies in all cases, reaching resolutions of 256^3 zones or 350,000 particles. We conclude that the observed long lifetimes and supersonic motions in molecular clouds must be due to external driving, as undriven turbulence decays far too fast to explain the observations.

This review summarizes the argument for molecular clouds being dominated by turbulence, most likely super-Alfvenic turbulence. Five lines of observational evidence are given: molecular linewidths and line shapes, nonequilibrium chemical abundances, fractal cloud shapes, measurements of the dispersion of dust extinction through clouds, and Zeeman measurements of the field strength. Recent models by Padoan & Nordlund are summarized, that show that super-Alfvenic turbulence appears consistent with these observations. I then present recent computations of my own, with numerical resolution as high as 256^3 zones, confirming the basic picture proposed by Padoan & Nordlund, but showing that the decay timescales they quote are actually too long. My computations show that decaying turbulence loses its kinetic energy with a ~1/t dependence on time regardless of whether the turbulence is subsonic, supersonic, isothermal, adiabatic, unmagnetized or magnetized. Finally, the implications for star formation and turbulence theory are discussed.

The large H II complex N11 in the Large Magellanic Cloud contains OB associations at several different stages in their life histories. We have obtained ROSAT PSPC and HRI X-ray observations, Curtis Schmidt CCD images, echelle spectra in H-alpha and [N II] lines, and IUE interstellar absorption line observations of this region. The central bubble of N11 has an X-ray luminosity a factor of only 3-7 brighter than predicted for an energy-conserving superbubble, making this the first detection of X-ray emission from a superbubble without a strong X-ray excess. The region N11B contains an extremely young OB association analogous to the central association of the Carina nebula, apparently still embedded in its natal molecular cloud. We find that N11B emits diffuse X-ray emission, probably powered by stellar winds. Finally, we compare the tight cluster HD32228 in N11 to R136 in 30 Dor. The latter is a strong X-ray source, while the former is not detected, showing that strong X-ray emission from compact objects is not a universal property of such tight clusters.

We compute the nonlinear development of the instabilities in C-shocks first described by Wardle, using a version of the ZEUS code modified to include a semi-implicit treatment of ambipolar diffusion. We find that, in three dimensions, thin sheets parallel to the shock velocity and perpendicular to the magnetic field lines form. High resolution, two-dimensional models show that the sheets are confined by the Brandenburg & Zweibel ambipolar diffusion singularity, forcing them to numerically unresolvable thinness. Hot and cold regions form around these filaments, disrupting the uniform temperature structure characteristic of a steady-state C-shock. This filamentary region steadily grows as the shock progresses. We compare steady-state to unstable C-shocks, showing excitation diagrams, line ratios, and line profiles for molecular hydrogen lines visible in the K-band, with the Infrared Space Observatory, and with NICMOS on the Hubble Space Telescope.

Shock waves in molecular clouds should evolve into continuous or C-type structures due to the magnetic field and ion-neutral friction. We here determine whether and how this is achieved through plane-parallel numerical simulations using an extended version of ZEUS. We first describe and test the adapted code against analytical results, laying the necessary foundations for subsequent works on supersonic ambipolar diffusion, including C-type jets and shock instability. The evolution away from jump shocks toward the numerous steady C-shock sub-types is then investigated. The evolution passes through four stages, which possess distinctive observational properties. The time scales and length scales cover broad ranges. Specific results are included for shock types including switch, absorber, neutralised, oblique, transverse and intermediate. Only intermediate Type II shocks and `slow shocks', including switch-off shocks, remain as J-type under the low ion levels assumed. Other shocks transform via a steadily growing neutral precursor to a diminishing jump. For neutralised shocks, this takes the form of an extended long-lived ramp. Molecular hydrogen emission signatures are presented. After the jump speed has dropped to under 25 km/s, a non-dissociative jump section can dominate the spectra for a long period. This produces a high-excitation spectrum. Once the jump has further weakened, to < 8 km/s, the fully developed ion front is responsible for brisk progress towards a constant C-type excitation. The time scale for emission-line variations is ~(6 /n_i) yr, where n_i is the pre-shock ion number density.

We present new images of the OMC-1 molecular hydrogen outflow, made using long-slit spectroscopy in order to accurately subtract the underlying continuum emission. These images reveal an extremely clumpy, quasi-spherical inner shell that breaks up at larger radii into bow-shocks and trailing wakes in the north-west, as originally described by Allen & Burton (1993); a fainter counter-finger to the south-east is newly discovered in the present data. While the outflow appears to be broadly bipolar, this is probably due to an interaction between an initially spherical wind from the source and a large-scale density enhancement surrounding it, rather than direct collimation imposed close to the source. The clumpy appearance of the inner shell confirms the prediction of the recent model of Stone, Xu, & Mundy (1995), in which a spherical and time-varying wind fragments a swept-up shell, producing high-velocity shrapnel, which in turn drives bow-shocks into the surrounding gas, resulting in the observed "fingers". As an alternative to the single varying source proposed by Stone et al., we speculate that several young sources in the BN-KL cluster may have been responsible for first sweeping up the shell and then fragmenting it.

With the High Resolution Spectrograph aboard the Hubble Space Telescope we obtained high signal-to-noise (S/N > 200 - 600 per 17 km/s resolution element) spectra of narrow absorption lines toward the Wolf-Rayet star HD 50896. The ring nebula S308 that surrounds this star is thought to be caused by a pressure-driven bubble bounded by circumstellar gas (most likely from a red supergiant or luminous blue variable progenitor) pushed aside by a strong stellar wind. Our observation has shown for the first time that blueshifted (approximately 70 km/s relative to the star) absorption components of C IV and N V arise in a conduction front between the hot interior of the bubble and the cold shell of swept-up material. These lines set limits on models of the conduction front. Nitrogen in the shell appears to be overabundant by a factor ~10. The P Cygni profiles of N V and C IV are variable, possibly due to a suspected binary companion to HD 50896.