Building on our previous hydrodynamic study of the angular momenta of cloud cores formed during gravitational collapse of star-forming molecular gas in our previous work, we now examine core properties assuming ideal magnetohydrodynamics (MHD). Using the same sink-patch implementation for the \emph{Athena} MHD code, we characterize the statistical properties of cores, including the mass accretion rates, specific angular momenta, and alignments between the magnetic field and the spin axis of the core on the $0.1 \ \mathrm{pc}$ scale. Our simulations, which reproduce the observed relation between magnetic field strength and gas density, show that magnetic fields can help collimate low density flows and help seed the locations of filamentary structures. Consistent with our previous purely hydrodynamic simulations, stars (sinks) form within the heterogeneous environments of filaments, such that accretion onto cores is highly episodic leading to short-term variability but no long-term monotonic growth of the specific angular momenta. With statistical characterization of protostellar cores properties and behaviors, we aim to provide a starting point for building more realistic and self-consistent disk formation models, helping to address whether magnetic fields can prevent the development of (large) circumstellar disks in the ideal MHD limit.

Numerical models of gas inflow towards a supermassive black hole (SMBH) show that star formation may occur in such an environment through the growth of a gravitationally unstable gas disc. We consider the effect of nuclear activity on such a scenario. We present the first three-dimensional grid-based radiative hydrodynamic simulations of direct collisions between infalling gas streams and a $4 \times 10^6~\text{M}_\odot$ SMBH, using ray-tracing to incorporate radiation consistent with an active galactic nucleus (AGN). We assume inflow masses of $ \approx 10^5~\text{M}_\odot$ and explore radiation fields of 10% and 100% of the Eddington luminosity ($L_\text{edd}$). We follow our models to the point of central gas disc formation preceding star formation and use the Toomre Q parameter ($Q_T$) to test for gravitational instability. We find that radiation pressure from UV photons inhibits inflow. Yet, for weak radiation fields, a central disc forms on timescales similar to that of models without feedback. Average densities of $> 10^{8}~\text{cm}^{-3}$ limit photo-heating to the disc surface allowing for $Q_T\approx1$. For strong radiation fields, the disc forms more gradually resulting in lower surface densities and larger $Q_T$ values. Mass accretion rates in our models are consistent with 1%--60% of the Eddington limit, thus we conclude that it is unlikely that radiative feedback from AGN activity would inhibit circumnuclear star formation arising from a massive inflow event.

We present the results of a suite of numerical simulations designed to explore the origin of the angular momenta of protostellar cores. Using the hydrodynamic grid code \emph{Athena} with a sink implementation, we follow the formation of protostellar cores and protostars (sinks) from the subvirial collapse of molecular clouds on larger scales to investigate the range and relative distribution of core properties. We find that the core angular momenta are relatively unaffected by large-scale rotation of the parent cloud; instead, we infer that angular momenta are mainly imparted by torques between neighboring mass concentrations and exhibit a log-normal distribution. Our current simulation results are limited to size scales $\sim 0.05$~pc ($\sim 10^4 \rm AU$), but serve as first steps toward the ultimate goal of providing initial conditions for higher-resolution studies of core collapse to form protoplanetary disks.

Numerical simulations of star formation have found that a power-law mass function can develop at high masses. In a previous paper, we employed isothermal simulations which created large numbers of sinks over a large range in masses to show that the power law exponent of the mass function, $dN/d\log M \propto M^{\Gamma}$, asymptotically and accurately approaches $\Gamma = -1.$ Simple analytic models show that such a power law can develop if the mass accretion rate $\dot{M} \propto M^2$, as in Bondi-Hoyle accretion; however, the sink mass accretion rates in the simulations show significant departures from this relation. In this paper we show that the expected accretion rate dependence is more closely realized provided the gravitating mass is taken to be the sum of the sink mass and the mass in the near environment. This reconciles the observed mass functions with the accretion rate dependencies, and demonstrates that power-law upper mass functions are essentially the result of gravitational focusing, a mechanism present in, for example, the competitive accretion model.

We present the results from an analysis of deep Herschel Space Observatory observations of six nearby dwarf galaxies known to host galactic-scale winds. The superior far-infrared sensitivity and angular resolution of Herschel have allowed detection of cold circumgalactic dust features beyond the stellar components of the host galaxies traced by Spitzer 4.5 $\mu$m images. Comparisons of these cold dust features with ancillary data reveal an imperfect spatial correlation with the ionized gas and warm dust wind components. We find that typically $\sim$10-20\% of the total dust mass in these galaxies resides outside of their stellar disks, but this fraction reaches $\sim$60\% in the case of NGC 1569. This galaxy also has the largest metallicity (O/H) deficit in our sample for its stellar mass. Overall, the small number of objects in our sample precludes drawing strong conclusions on the origin of the circumgalactic dust. We detect no statistically significant trends with star formation properties of the host galaxies, as might be expected if the dust were lifted above the disk by energy inputs from on-going star formation activity. Although a case for dust entrained in a galactic wind is seen in NGC 1569, in all cases, we cannot rule out the possibility that some of the circumgalactic dust might be associated instead with gas accreted or removed from the disk by recent galaxy interaction events, or that it is part of the outer gas-rich portion of the disk that lies below the sensitivity limit of the Spitzer 4.5 $\mu$m data.

We present numerical simulations of reorienting oblique shocks that form in the collision layer between magnetized colliding flows. Reorientation aligns parsec-scale post-shock filaments normal to the background magnetic field. We find that reorientation begins with pressure gradients between the collision region and the ambient medium. This drives a lateral expansion of post-shock gas, which reorients the growing filament from the outside-in (i.e. from the flow-ambient boundary, toward the colliding flows axis). The final structures of our simulations resemble polarization observations of filaments in Taurus and Serpens South, as well as the integral-shaped filament in Orion A. Given the ubiquity of colliding flows in the interstellar medium, shock reorientation may be relevant to the formation of filaments normal to magnetic fields.

Turbulence models attempt to account for unresolved dynamics and diffusion in hydrodynamical simulations. We develop a common framework for two-equation Reynolds-Averaged Navier-Stokes (RANS) turbulence models, and we implement six models in the Athena code. We verify each implementation with the standard subsonic mixing layer, although the level of agreement depends on the definition of the mixing layer width. We then test the validity of each model into the supersonic regime, showing that compressibility corrections can improve agreement with experiment. For models with buoyancy effects, we also verify our implementation via the growth of the Rayleigh-Taylor instability in a stratified medium. The models are then applied to the ubiquitous astrophysical shock-cloud interaction in three dimensions. We focus on the mixing of shock and cloud material, comparing results from turbulence models to high-resolution simulations (up to 200 cells per cloud radius) and ensemble-averaged simulations. We find that the turbulence models lead to increased spreading and mixing of the cloud, although no two models predict the same result. Increased mixing is also observed in inviscid simulations at resolutions greater than 100 cells per radius, which suggests that the turbulent mixing begins to be resolved.

We present three-dimensional magneto-hydrodynamical simulations of the self-gravitating interstellar medium (ISM) in a periodic (256 pc)$^3$ box with a mean number density of 0.5 cm$^{-3}$. At a fixed supernova rate we investigate the multi-phase ISM structure, H$_{2}$ molecule formation and density-magnetic field scaling for varying initial magnetic field strengths (0, $6\times 10^{-3}$, 0.3, 3 $\mu$G). All magnetic runs saturate at mass weighted field strengths of $\sim$ 1 $-$ 3 $\mu$G but the ISM structure is notably different. With increasing initial field strengths (from $6\times 10^{-3}$ to 3 $\mu$G) the simulations develop an ISM with a more homogeneous density and temperature structure, with increasing mass (from 5% to 85%) and volume filling fractions (from 4% to 85%) of warm (300 K $<$ T $<$ 8000 K) gas, with decreasing volume filling fractions (VFF) from $\sim$ 35% to $\sim$ 12% of hot gas (T $> 10^5$ K) and with a decreasing H$_{2}$ mass fraction (from 70% to $<$ 1%). Meanwhile the mass fraction of gas in which the magnetic pressure dominates over the thermal pressure increases by a factor of 10, from 0.07 for an initial field of $6\times 10^{-3}$ $\mu$G to 0.7 for a 3 $\mu$G initial field. In all but the simulations with the highest initial field strength self-gravity promotes the formation of dense gas and H$_{2}$, but does not change any other trends. We conclude that magnetic fields have a significant impact on the multi-phase, chemical and thermal structure of the ISM and discuss potential implications and limitations of the model.

High levels of deuterium fraction in N$_2$H$^+$ are observed in some pre-stellar cores. Single-zone chemical models find that the timescale required to reach observed values ($D_{\rm frac}^{{\rm N}_2{\rm H}^+} \equiv {\rm N}_2{\rm D}^+/{\rm N}_2{\rm H}^+ \gtrsim 0.1$) is longer than the free-fall time, possibly ten times longer. Here, we explore the deuteration of turbulent, magnetized cores with 3D magnetohydrodynamics simulations. We use an approximate chemical model to follow the growth in abundances of N$_2$H$^+$ and N$_2$D$^+$. We then examine the dynamics of the core using each tracer for comparison to observations. We find that the velocity dispersion of the core as traced by N$_2$D$^+$ appears slightly sub-virial compared to predictions of the Turbulent Core Model of McKee & Tan, except at late times just before the onset of protostar formation. By varying the initial mass surface density, the magnetic energy, the chemical age, and the ortho-to-para ratio of H$_2$, we also determine the physical and temporal properties required for high deuteration. We find that low initial ortho-to-para ratios ($\lesssim 0.01$) and/or multiple free-fall times ($\gtrsim 3$) of prior chemical evolution are necessary to reach the observed values of deuterium fraction in pre-stellar cores.

The presence of short-lived radioisotopes (SLRs) in solar system meteorites has been interpreted as evidence that the solar system was exposed to a supernova shortly before or during its formation. Yet results from hydrodynamical models of SLR injection into the proto-solar cloud or disc suggest that gas-phase mixing may not be efficient enough to reproduce the observed abundances. As an alternative, we explore the injection of SLRs via dust grains as a way to overcome the mixing barrier. We numerically model the interaction of a supernova remnant containing SLR-rich dust grains with a nearby molecular cloud. The dust grains are subject to drag forces and both thermal and non-thermal sputtering. We confirm that the expanding gas shell stalls upon impact with the dense cloud and that gas-phase SLR injection occurs slowly due to hydrodynamical instabilities at the cloud surface. In contrast, dust grains of sufficient size (> 1 micron) decouple from the gas and penetrate into the cloud within 0.1 Myr. Once inside the cloud, the dust grains are destroyed by sputtering, releasing SLRs and rapidly enriching the dense (potentially star-forming) regions. Our results suggest that SLR transport on dust grains is a viable mechanism to explain SLR enrichment.

Our three-dimensional hydro-dynamical simulations of starbursts examine the formation of superbubbles over a range of driving luminosities and mass loadings that determine superbubble growth and wind velocity. From this we determine the relationship between the velocity of a galactic wind and the power of the starburst. We find a threshold for the formation of a wind, above which the speed of the wind is not affected by grid resolution or the temperature floor of our radiative cooling. We investigate the effect two different temperature floors in our radiative cooling prescription have on wind kinematics and content. We find that cooling to $10$ K instead of to $10^4$ K increases the mass fraction of cold neutral and hot X-ray gas in the galactic wind while halving that in warm H$\alpha$. Our simulations show the mass of cold gas transported into the lower halo does not depend on the starburst strength. Optically bright filaments form at the edge of merging superbubbles, or where a cold dense cloud has been disrupted by the wind. Filaments formed by merging superbubbles will persist and grow to $>400$ pc in length if anchored to a star forming complex. Filaments embedded in the hot galactic wind contain warm and cold gas that moves $300-1200$ km s$^{-1}$ slower than the surrounding wind, with the coldest gas hardly moving with respect to the galaxy. Warm and cold matter in the galactic wind show asymmetric absorption profiles consistent with observations, with a thin tail up to the wind velocity.

Using synthetic absorption lines generated from 3D hydro-dynamical simulations we explore how the velocity of a starburst-driven galactic wind correlates with the star formation rate (SFR) and SFR density. We find strong correlations until the scaling relations flatten abruptly at a point set by the mass loading of the starburst. Below this point the scaling relation depends on the temperature regime being probed by the absorption line, not on the mass loading. The exact scaling relation depends on whether the maximum or mean velocity of the absorption line is used. We find that the outflow velocity of neutral gas is four to five times lower than the average velocity of the hottest gas, with the difference in velocity between the neutral and ionized gas increasing with gas ionization. Thus, absorption lines of neutral or low ionized gas will underestimate the outflow velocity of hot gas, severely underestimating outflow energetics.

We present a proof-of-concept study of a method to estimate the inclination angle of compact high velocity clouds (CHVCs), i.e. the angle between a CHVC's trajectory and the line-of-sight. The inclination angle is derived from the CHVC's morphology and kinematics. We calibrate the method with numerical simulations, and we apply it to a sample of CHVCs drawn from HIPASS. Implications for CHVC distances are discussed.

The colliding flows (CF) model is a well-supported mechanism for generating molecular clouds. However, to-date most CF simulations have focused on the formation of clouds in the normal-shock layer between head-on colliding flows. We performed simulations of magnetized colliding flows that instead meet at an oblique-shock layer. Oblique shocks generate shear in the post-shock environment, and this shear creates inhospitable environments for star formation. As the degree of shear increases (i.e. the obliquity of the shock increases), we find that it takes longer for sink particles to form, they form in lower numbers, and they tend to be less massive. With regard to magnetic fields, we find that even a weak field stalls gravitational collapse within forming clouds. Additionally, an initially oblique collision interface tends to reorient over time in the presence of a magnetic field, so that it becomes normal to the oncoming flows. This was demonstrated by our most oblique shock interface, which became fully normal by the end of the simulation.

The Smith Cloud is a gaseous high-velocity cloud (HVC) in an advanced state of accretion, only 2.9 kpc below the Galactic plane and due to impact the disk in 27 Myr. It is unique among HVCs in having a known distance (12.4+/-1.3 kpc) and a well-constrained 3D velocity (296 km/s), but its origin has long remained a mystery. Here we present the first absorption-line measurements of its metallicity, using HST/COS UV spectra of three AGN lying behind the Cloud together with Green Bank Telescope 21 cm spectra of the same directions. Using Voigt-profile fitting of the S II 1250, 1253, 1259 triplet together with ionization corrections derived from photoionization modeling, we derive the sulfur abundance in each direction; a weighted average of the three measurements gives [S/H]=-0.28+/-0.14, or 0.53+0.21-0.15 solar metallicity. The finding that the Smith Cloud is metal-enriched lends support to scenarios where it represents recycled Galactic material rather than the remnant of a dwarf galaxy or accreting intergalactic gas. The metallicity and trajectory of the Cloud are both indicative of an origin in the outer disk. However, its large mass and prograde kinematics remain to be fully explained. If the cloud has accreted cooling gas from the corona during its fountain trajectory, as predicted in recent theoretical work, its current mass would be higher than its launch mass, alleviating the mass concern.

To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models.

We present a detailed analysis of deep far-infrared observations of the nearby edge-on star-forming galaxy NGC 4631 obtained with the Herschel Space Observatory. Our PACS images at 70 and 160 um show a rich complex of filaments and chimney-like features that extends up to a projected distance of 6 kpc above the plane of the galaxy. The PACS features often match extraplanar Halpha, radio-continuum, and soft X-ray features observed in this galaxy, pointing to a tight disk-halo connection regulated by star formation. On the other hand, the morphology of the colder dust component detected on larger scale in the SPIRE 250, 350, and 500 um data matches the extraplanar H~I streams previously reported in NGC 4631 and suggests a tidal origin. The PACS 70/160 ratios are elevated in the central ~3.0 kpc region above the nucleus of this galaxy (the "superbubble"). A pixel-by-pixel analysis shows that dust in this region has a higher temperature and/or an emissivity with a steeper spectral index (beta > 2) than the dust in the disk, possibly the result of the harsher environment in the superbubble. Star formation in the disk seems energetically insufficient to lift the material out of the disk, unless it was more active in the past or the dust-to-gas ratio in the superbubble region is higher than the Galactic value. Some of the dust in the halo may also have been tidally stripped from nearby companions or lifted from the disk by galaxy interactions.

A comparison of accretion and (turbulent) magnetic diffusion timescales for sheets and filaments demonstrates that dense star-forming clouds generally will -- under realistic conditions -- become supercritical due to mass accretion on timescales at least an order of magnitude shorter than ambipolar and/or turbulent diffusion timescales. Thus, ambipolar or turbulent diffusion -- while present -- is unlikely to control the formation of cores and stars.

Numerical simulations have explored the possibility to form molecular clouds through either a quasi-static, self-gravitating mechanism or the collision of gas streams or lower-density clouds. They also quantitatively predict the distribution of matter at the transition from atomic to molecular gases. We aim to observationally test these models by studying the environment of W43, a molecular cloud complex near the tip of the Galactic long bar. Using Galaxy-wide HI and 12CO surveys we searched for gas flowing toward the W43 molecular cloud complex. We also estimated the HI and H2 mass surface densities to constrain the transition from atomic to molecular gas around and within W43. We found 3 cloud ensembles within the position-velocity diagrams of 12CO and HI gases. They are separated by 20km/s along the line of sight and extend into the 13CO velocity structure of W43. Since their velocity gradients are consistent with free-fall, they could be nearby clouds attracted by, and streaming toward, the W43 10^7Msun potential well. We show that the HI surface density, Sigma_HI=45-85Msun/pc2, does not reach any threshold level but increases when entering the 130pc-wide molecular complex previously defined. This suggests that an equilibrium between H2 formation and photodissociation has not yet been reached. The H2-to-HI ratio measured over the W43 region and its surroundings, R_H2~3.5, is high, indicating that most of the gas is already in molecular form in W43 and in structures several hundreds of parsecs downstream along the Scutum-Centaurus arm. The W43 molecular cloud complex may have formed, and in fact may still be accreting mass from the agglomeration of clouds. Already in the molecular-dominated regime, most of these clouds are streaming from the Scutum-Centaurus arm. This is in clear disagreement with quasi-static and steady-state models of molecular cloud formation.

We aim to give a full description of the distribution and location of dense molecular clouds in the giant molecular cloud complex W43. It has previously been identified as one of the most massive star-forming regions in our Galaxy. To trace the moderately dense molecular clouds in the W43 region, we initiated an IRAM 30m large program, named W43-HERO, covering a large dynamic range of scales (from 0.3 to 140 pc). We obtained on-the-fly-maps in 13CO (2-1) and C18O (2-1) with a high spectral resolution of 0.1 km/s and a spatial resolution of 12". These maps cover an area of ~1.5 square degrees and include the two main clouds of W43, as well as the lower density gas surrounding them. A comparison with Galactic models and previous distance calculations confirms the location of W43 near the tangential point of the Scutum arm at a distance from the Sun of approximately 6 kpc. The resulting intensity cubes of the observed region are separated into sub-cubes, centered on single clouds which are then analyzed in detail. The optical depth, excitation temperature, and H2 column density maps are derived out of the 13CO and C18O data. These results are then compared with those derived from Herschel dust maps. The mass of a typical cloud is several 10^4 solar masses while the total mass in the dense molecular gas (>100 cm^-3) in W43 is found to be about 1.9e6 solar masses. Probability distribution functions obtained from column density maps derived from molecular line data and Herschel imaging show a log-normal distribution for low column densities and a power-law tail for high densities. A flatter slope for the molecular line data PDF may imply that those selectively show the gravitationally collapsing gas.

In an extension of Fischera & Martin (2012a) and Heitsch (2013), two aspects of the evolution of externally pressurized, hydrostatic filaments are discussed. (a) The free-fall accretion of gas onto such a filament will lead to filament parameters (specifically, FWHM--column density relations) inconsistent with the observations of Arzoumanian et al. (2011), except for two cases: For low-mass, isothermal filaments, agreement is found as in the analysis by Fischera & Martin (2012b). Magnetized cases, for which the field scales weakly with the density as $B\propto n^{1/2}$, also reproduce observed parameters. (b) Realistically, the filaments will be embedded not only in gas of non-zero pressure, but also of non-zero density. Thus, the appearance of sheet-embedded filaments is explored. Generating a grid of filament models and comparing the resulting column density ratios and profile shapes with observations suggests that the three-dimensional filament profiles are intrinsically flatter than isothermal, beyond projection and evolution effects.

The formation of high-mass stars is tightly linked to that of their parental clouds. We here focus on the high-density parts of W43, a molecular cloud undergoing an efficient event of formation. The cloud structure is studied with a column density image derived from Herschel continuum maps obtained at 70, 160, 250, 350, and 500 micron. We identify two high-column density filamentary clouds, quoted as the W43-MM1 and W43-MM2 ridges, which both account for 1.5x10^4 Msun gas mass above 10^23 cm-2 and within areas of 5 and 14pc^2, respectively. We used the N_2H^+ 1--0 line to confirm that the W43-MM1 and W43-MM2 ridges are structures coherent in velocity and gravitationally bound, despite their large velocity dispersion and ~5 kms line widths. The most intriguing result of the W43 large program is the bright wide-spread SiO 2--1 emission: 1--11 K kms$ stretching an area of ~28 pc^2. Concentrated toward the W43-MM1 and W43-MM2 ridges and their immediate surroundings, it leads to a total luminosity of L_SiO 2-1 ~4 10^4 K kms kpc^2pc^2. We measured a steep relation between the luminosity and velocity extent of the SiO~2--1 lines and propose to use it to distinguish the low-velocity shocks observed here from the more classical high-velocity ones associated with outflows of high-mass young stellar objects. We used state-of-the-art shock models to demonstrate that low-velocity (<10 kms^-1) shocks with a small amount (10%) of Si atoms initially in gas phase or in grain mantles can explain the observed SiO column density in W43. The spatial and velocity overlaps between the ridges high-density gas (n_H2>10^4-10^5 cm^-3) and the shocked SiO gas suggests that ridges could be forming via colliding flows driven by gravity and accompanied by low-velocity shocks. This mechanism may be the initial conditions for the formation of young massive clusters in these ridges.

Two aspects of filamentary molecular cloud evolution are addressed: (1) Exploring analytically the role of the environment for the evolution of filaments demonstrates that considering them in isolation (i.e. just addressing the fragmentation stability) will result in unphysical conclusions about the filament's properties. Accretion can also explain the observed decorrelation between FWHM and peak column density. (2) Free-fall accretion onto finite filaments can lead to the characteristic "fans" of infrared-dark clouds around star-forming regions. The fans may form due to tidal forces mostly arising at the ends of the filaments, consistent with numerical models and earlier analytical studies.

Observational evidence from local star-forming regions mandates that star formation occurs shortly after, or even during, molecular cloud formation. Models of the formation of molecular clouds in large-scale converging flows have identified the physical mechanisms driving the necessary rapid fragmentation. They also point to global gravitational collapse driving supersonic turbulence in molecular clouds. Previous cloud formation models have focused on turbulence generation, gravitational collapse, magnetic fields, and feedback. Here, we explore the effect of structure in the flow on the resulting clouds and the ensuing gravitational collapse. We compare two extreme cases, one with a collision between two smooth streams, and one with streams containing small clumps. We find that structured converging flows lead to a delay of local gravitational collapse ("star formation"). Thus, more gas has time to accumulate, eventually leading to a strong global collapse, and thus to a high star formation rate. Uniform converging flows fragment hydrodynamically early on, leading to the rapid onset of local gravitational collapse and an overall low sink formation rate.

We investigate the collapse of non-spherical substructures, such as sheets and filaments, which are ubiquitous in molecular clouds. Such non-spherical substructures collapse homologously in their interiors but are influenced by an edge effect that causes their edges to be preferentially accelerated. We analytically compute the homologous collapse timescales of the interiors of uniform-density, self-gravitating filaments and find that the homologous collapse timescale scales linearly with the aspect ratio. The characteristic timescale for an edge driven collapse mode in a filament, however, is shown to have a square root dependence on the aspect ratio. For both filaments and circular sheets, we find that selective edge acceleration becomes more important with increasing aspect ratio. In general, we find that lower dimensional objects and objects with larger aspect ratios have longer collapse timescales. We show that estimates for star formation rates, based upon gas densities, can be overestimated by an order of magnitude if the geometry of a cloud is not taken into account.

We present an interferometric kinematic study of morphologically complex protostellar envelopes based on observations of the dense gas tracers N2H+ and NH3. The strong asymmetric nature of most envelopes in our sample leads us to question the common interpretation of velocity gradients as rotation, given the possibility of projection effects in the observed velocities. Several "idealized" sources with well-ordered velocity fields and envelope structures are now analyzed in more detail. We compare the interferometric data to position-velocity diagrams of kinematic models for spherical rotating collapse and filamentary rotating collapse. For this purpose, we developed a filamentary parametrization of the rotating collapse model to explore the effects of geometric projection on the observed velocity structures. We find that most envelopes in our sample have PV structures that can be reproduced by an infalling filamentary envelope projected at different angles within the plane of the sky. The infalling filament produces velocity shifts across the envelope that can mimic rotation, especially when viewed at single-dish resolutions and the axisymmetric rotating collapse model does not uniquely describe any dataset. Furthermore, if the velocities are assumed to reflect rotation, then the inferred centrifugal radii are quite large in most cases, indicating significant fragmentation potential or more likely another component to the line-center velocity. We conclude that ordered velocity gradients cannot be interpreted as rotation alone when envelopes are non-axisymmetric and that projected infall velocities likely dominate the velocity field on scales larger than 1000 AU.

Infrared-dark clouds (IRDCs) are believed to be the birthplaces of rich clusters and thus contain the earliest phases of high-mass star formation. We use the Green Bank Telescope (GBT) and Very Large Array (VLA) maps of ammonia (NH3) in six IRDCs to measure their column density and temperature structure (Paper 1), and here, we investigate the kinematic structure and energy content. We find that IRDCs overall display organized velocity fields, with only localized disruptions due to embedded star formation. The local effects seen in NH3 emission are not high velocity outflows but rather moderate (few km/s) increases in the line width that exhibit maxima near or coincident with the mid-infrared emission tracing protostars. These line width enhancements could be the result of infall or (hidden in NH3 emission) outflow. Not only is the kinetic energy content insufficient to support the IRDCs against collapse, but also the spatial energy distribution is inconsistent with a scenario of turbulent cloud support. We conclude that the velocity signatures of the IRDCs in our sample are due to active collapse and fragmentation, in some cases augmented by local feedback from stars.

Most young stars in nearby molecular clouds have estimated ages of 1-2 Myr, suggesting that star formation is rapid. However, small numbers of stars in these regions with inferred ages of >= 5-10 Myr have been cited to argue that star formation is instead a slow, quasi-static process. When considering these alternative pictures it is important to recognize that the age spread in a given star-forming cloud is necessarily an upper limit to the timescales of local collapse, as not all spatially-distinct regions will start contracting at precisely the same instant. Moreover, star-forming clouds may dynamically evolve on timescales of a few Myr; in particular, global gravitational contraction will tend to yield increasing star formation rates with time due to generally increasing local gas densities. We show that two different numerical simulations of dynamic, flow-driven molecular cloud formation and evolution 1) predict age spreads for the main stellar population roughly consistent with observations, and 2) raise the possibility of forming small numbers of stars early in cloud evolution, before global contraction concentrates the gas and the bulk of the stellar population is produced. In general, the existence of a small number of older stars among a generally much-younger population is consistent with the picture of dynamic star formation, and may even provide clues to the time evolution of star-forming clouds.

We analytically investigate the modes of gravity-induced star formation possible in idealized finite molecular clouds where global collapse competes against both local Jeans instabilities and discontinuous edge instabilities. We examine these timescales for collapse in spheres, discs, and cylinders, with emphasis on the structure, size, and degree of internal perturbations required in order for local collapse to occur before global collapse. We find that internal, local collapse is more effective for the lower dimensional objects. Spheres and discs, if unsupported against global collapse, must either contain strong perturbations or must be unrealistically large in order for small density perturbations to collapse significantly faster than the entire cloud. We find, on the other hand, that filamentary geometry is the most favorable situation for the smallest perturbations to grow before global collapse overwhelms them and that filaments containing only a few Jeans masses and weak density perturbations can readily fragment. These idealized solutions are compared with simulations of star-forming regions in an attempt to delineate the role of global, local, and edge instabilities in determining the fragmentation properties of molecular clouds. The combined results are also discussed in the context of recent observations of Galactic molecular clouds.

We present a study of dense molecular gas kinematics in seventeen nearby protostellar systems using single-dish and interferometric molecular line observations. The non-axisymmetric envelopes around a sample of Class 0/I protostars were mapped in the N2H+ (J=1-0) tracer with the IRAM 30m, CARMA and PdBI as well as NH3 (1,1) with the VLA. The molecular line emission is used to construct line-center velocity and linewidth maps for all sources to examine the kinematic structure in the envelopes on spatial scales from 0.1 pc to ~1000 AU. The direction of the large-scale velocity gradients from single-dish mapping is within 45 degrees of normal to the outflow axis in more than half the sample. Furthermore, the velocity gradients are often quite substantial, the average being ~2.3 km\s\pc. The interferometric data often reveal small-scale velocity structure, departing from the more gradual large-scale velocity gradients. In some cases, this likely indicates accelerating infall and/or rotational spin-up in the inner envelope; the median velocity gradient from the interferometric data is ~10.7 km/s/pc. In two systems, we detect high-velocity HCO+ (J=1-0) emission inside the highest-velocity \nthp\ emission. This enables us to study the infall and rotation close to the disk and estimate the central object masses. The velocity fields observed on large and small-scales are more complex than would be expected from rotation alone, suggesting that complex envelope structure enables other dynamical processes (i.e. infall) to affect the velocity field.

Context: The large Coalsack dark cloud is one of the rare molecular clouds without clear signs of star formation. Aims: We investigate the dynamical properties of the gas within the Coalsack. Methods: The two highest extinction regions were mapped with the APEX telescope in 13CO(2-1) comprising a region of ~1 square degree. Results: In addition to the well-known, nearby gas component around -4km/s, we identified additional molecular gas components - in particular a second extended molecular cloud at a velocity of ~-30km/s and an estimated distance of ~3.1kpc - that dominate the column density and visual extinction distributions in the northeastern part of the Coalsack. Although comprising ~2600Msun, the mass of this distant cloud is distributed over an extent of ~73pc, much larger than typical high-mass infrared dark clouds. Its filamentary structure is consistent with a compressible gaseous self-gravitating cylinder, and its low mass per length indicates that it may be stable against gravitational collapse. We find barely any mid-infrared emission in archival MSX data, which is indicative of almost no star-formation activity in the near and far cloud complexes. The nearby clouds have narrow, almost thermal velocity dispersions with median values between 0.2 and 0.4km/s, which is also consistent with low star-formation activity. Only Tapia's Globule 2 exhibits a velocity dispersion increase toward the extinction peak and peak-velocity gradients over the core, which is indicative of a state of elevated dynamical properties. Conclusions: The Coalsack is not one single coherent structure, but consists of several cloud complexes nearby as well as at several kpc distance. All studied clouds appear as starless low-turbulence regions that may not even collapse in the future. Only one globule exhibits more dynamical signatures and is a good candidate for present/future star formation.

It has been recently shown that molecular clouds do not exhibit a unique shape for the column density probability distribution function (Npdf). Instead, clouds without star formation seem to possess a lognormal distribution, while clouds with active star formation develope a power-law tail at high column densities. The lognormal behavior of the Npdf has been interpreted in terms of turbulent motions dominating the dynamics of the clouds, while the power-law behavior occurs when the cloud is dominated by gravity. In the present contribution we use thermally bi-stable numerical simulations of cloud formation and evolution to show that, indeed, these two regimes can be understood in terms of the formation and evolution of molecular clouds: a very narrow lognormal regime appears when the cloud is being assembled. However, as the global gravitational contraction occurs, the initial density fluctuations are enhanced, resulting, first, in a wider lognormal Npdf, and later, in a power-law Npdf. We thus suggest that the observed Npdf of molecular clouds are a manifestation of their global gravitationally contracting state. We also show that, contrary to recent suggestions, the exact value of the power-law slope is not unique, as it depends on the projection in which the cloud is being observed.

The fragmentation of shocked flows in a thermally bistable medium provides a natural mechanism to form turbulent cold clouds as precursors to molecular clouds. Yet because of the large density and temperature differences and the range of dynamical scales involved, following this process with numerical simulations is challenging. We compare two-dimensional simulations of flow-driven cloud formation without self-gravity, using the Lagrangian Smoothed Particle Hydrodynamics (SPH) code VINE and the Eulerian grid code Proteus. Results are qualitatively similar for both methods, yet the variable spatial resolution of the SPH method leads to smaller fragments and thinner filaments, rendering the overall morphologies different. Thermal and hydro-dynamical instabilities lead to rapid cooling and fragmentation into cold clumps with temperatures below 300K. For clumps more massive than 1 Msun/pc, the clump mass function has an average slope of -0.8. The internal velocity dispersion of the clumps is nearly an order of magnitude smaller than their relative motion, rendering it subsonic with respect to the internal sound speed of the clumps, but supersonic as seen by an external observer. For the SPH simulations most of the cold gas resides at temperatures below 100K, while the grid-based models show an additional, substantial component between 100 and 300K. Independently of the numerical method our models confirm that converging flows of warm neutral gas fragment rapidly and form high-density, low-temperature clumps as possible seeds for star formation.

We used a large database extracted from Galaxy-wide surveys of H {\scriptsize I}, $^{13}$CO 1-0, $8 \micron$ and $870 \micron$ continuum to trace diffuse atomic gas, low- to medium-density molecular gas, high-density molecular gas, and star formation activity which we complemented by dedicated $^{12}$CO 2--1, 3--2 observations of the W43 region. From the detailed 3D (space-space-velocity) analysis of the molecular and atomic cloud tracers through the region and despite its wide velocity range (\emph{FWHM}$\sim22.3 \kms$ around $\vlsr\sim 95.9 \kms$), we identified W43 as a large (equivalent diameter $\sim 140$ pc) and coherent complex of molecular clouds which is surrounded by an atomic gas envelope (equivalent diameter $\sim 290$ pc). We measured the total mass of this newly-identified molecular complex ($M_{{\tiny total}}\sim 7.1 \times 10^6 \msun$), the mass contained in dense $870 \micron$ clumps ($<5$ pc dense cloud structures, $M_{{\tiny clumps}}\sim 8.4 \times 10^5 \msun$) and conclude that W43 is particularly massive and concentrated. The distance we assume for the W43 complex is 6 kpc from the Sun, which may place it at the meeting point of the Scutum-Centaurus (or Scutum-Crux) Galactic arm and the Bar, a dynamically complex region where high-velocity streams could easily collide. The star formation rate of W43 is suggested not to be steady but it is increasing from $\, \sim 0.01 \msun\, {yr}^{-1}$ (measured from its $8\,\mu$m luminosity) to $\sim 0.1 \msun\, {yr}^{-1}$ (measured from its molecular content). From the global properties of W43, we claim it is an extreme molecular complex in the Milky Way and it could even be forming starburst clusters in the near future.

The expansion and collision of two wind-blown superbubbles is investigated numerically. Our models go beyond previous simulations of molecular cloud formation from converging gas flows by exploring this process with realistic flow parameters, sizes and timescales. The superbubbles are blown by time-dependent winds and supernova explosions, calculated from population synthesis models. They expand into a uniform or turbulent diffuse medium. We find that dense, cold gas clumps and filaments form naturally in the compressed collision zone of the two superbubbles. Their shapes resemble the elongated, irregular structure of observed cold, molecular gas filaments and clumps. At the end of the simulations, between 65 and 80 percent of the total gas mass in our simulation box is contained in these structures. The clumps are found in a variety of physical states, ranging from pressure equilibrium with the surrounding medium to highly under-pressured clumps with large irregular internal motions and structures which are rotationally supported.

Observations from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey of the tail of Complex C are presented and the halo clouds associated with this complex cataloged. The properties of the Complex C clouds are compared to clouds cataloged at the tail of the Magellanic Stream to provide insight into the origin and destruction mechanism of Complex C. Magellanic Stream and Complex C clouds show similarities in their mass distributions (slope = -0.7 and -0.6, respectively) and have a common linewidth of 20 - 30 km/s (indicative of a warm component), which may indicate a common origin and/or physical process breaking down the clouds. The clouds cataloged at the tail of Complex C extend over a mass range of 10^1.1 to 10^4.8 solar masses, sizes of 10^1.2 to 10^2.6 pc, and have a median volume density of 0.065 cm^(-3) and median pressure of (P/k) = 580 K cm^{-3}. We do not see a prominent two-phase structure in Complex C, possibly due to its low metallicity and inefficient cooling compared to other halo clouds. From assuming the Complex C clouds are in pressure equilibrium with a hot halo medium, we find a median halo density of 5.8 x 10^(-4) cm^(-3), which given a constant distance of 10 kpc, is at a z-height of ~3 kpc. Using the same argument for the Stream results in a median halo density of 8.4 x 10^(-5) x (60kpc/d) cm^(-3). These densities are consistent with previous observational constraints and cosmological simulations. We also assess the derived cloud and halo properties with three dimensional grid simulations of halo HI clouds and find the temperature is generally consistent within a factor of 1.5 and the volume densities, pressures and halo densities are consistent within a factor of 3.

We discuss the nature of the velocity dispersion vs. size relation for molecular clouds. In particular, we add to previous observational results showing that the velocity dispersions in molecular clouds and cores are not purely functions of spatial scale but involve surface gas densities as well. We emphasize that hydrodynamic turbulence is required to produce the first condensations in the progenitor medium. However, as the cloud is forming, it also becomes bound, and gravitational accelerations dominate the motions. Energy conservation in this case implies $|E_g| \sim E_k$, in agreement with observational data, and providing an interpretation for two recent observational results: the scatter in the $\delta v-R$ plane, and the dependence of the velocity dispersion on the surface density ${\delta v^2/ R} \propto \Sigma$. We argue that the observational data are consistent with molecular clouds in a state of hierarchical gravitational collapse, i.e., developing local centers of collapse throughout the whole cloud while the cloud itself is collapsing, and making equilibrium unnecessary at all stages prior to the formation of actual stars. Finally, we discuss how this mechanism need not be in conflict with the observed star formation rate.

Extinction maps at 8 micron from the Spitzer Space Telescope show that many Class 0 protostars exhibit complex, irregular, and on-axisymmetric structure within the densest regions of their dusty envelopes. Many of the systems have highly irregular and on-axisymmetric morphologies on scales $\sim$1000 AU, with a quarter of the sample exhibiting filamentary or flattened dense structures. Complex envelope structure is observed in regions spatially distinct from outflow cavities, and the densest structures often show no systematic alignment perpendicular to the cavities. We suggest that the observed envelope complexity is the result of collapse from protostellar cores with initially non-equilibrium structures. The striking non-axisymmetry in many envelopes could provide favorable conditions for the formation of binary systems. We then show that the kinematics around L1165 as probed with N2H+ are indicative of asymmetric infall; the velocity gradient is not perpendicular to the outflow.

We report a set of numerical experiments aimed at addressing the applicability of competitive accretion to explain the high-mass end of the stellar initial mass function in a sheet geometry with shallow gravitational potential, in contrast to most previous simulations which have assumed formation in a cluster gravitational potential. Our flat cloud geometry is motivated by models of molecular cloud formation due to large-scale flows in the interstellar medium. The experiments consisted of SPH simulations of gas accretion onto sink particles formed rapidly from Jeans-unstable dense clumps placed randomly in the finite sheet. These simplifications allow us to study accretion with a minimum of free parameters, and to develop better statistics on the resulting mass spectra. We considered both clumps of equal mass and gaussian distributions of masses, and either uniform or spatially-varying gas densities. In all cases, the sink mass function develops a power law tail at high masses, with $dN/dlog M \propto M^{-\Gamma}$. The accretion rates of individual sinks follow $\dot{M} \propto M^2$ at high masses; this results in a continual flattening of the slope of the mass function towards an asymptotic form $\Gamma \sim 1$ (where the Salpeter slope is $\Gamma = 1.35$). The asymptotic limit is most rapidly reached when starting from a relatively broad distribution of initial sink masses. In general the resulting upper mass slope is correlated with the maximum sink mass; higher sink masses are found in simulations with flatter upper mass slopes. Although these simulations are of a highly idealized situation, the results suggest that competitive accretion may be relevant in a wider variety of environments than previously considered, and in particular that the upper mass distribution may generally evolve towards a limiting value of $\Gamma \sim 1$.

The HII complex N159 in the Large Magellanic Cloud (LMC) is used to study massive star formation in different environments, as it contains three giant molecular clouds (GMCs) that have similar sizes and masses but exhibit different intensities of star formation. We identify candidate massive young stellar objects (YSOs) using infrared photometry, and model their SEDs to constrain mass and evolutionary state. Good fits are obtained for less evolved Type I, I/II, and II sources. Our analysis suggests that there are massive embedded YSOs in N159B, a maser source, and several ultracompact HII regions. Massive O-type YSOs are found in GMCs N159-E and N159-W, which are associated with ionized gas, i.e., where massive stars formed a few Myr ago. The third GMC, N159-S, has neither O-type YSOs nor evidence of previous massive star formation. This correlation between current and antecedent formation of massive stars suggests that energy feedback is relevant. We present evidence that N159-W is forming YSOs spontaneously, while collapse in N159-E may be triggered. Finally, we compare star formation rates determined from YSO counts with those from integrated H-alpha and 24 micron luminosities and expected from gas surface densities. Detailed dissection of extragalactic GMCs like the one presented here is key to revealing the physics underlying commonly used star formation scaling laws.

Given the importance of shear flows for astrophysical gas dynamics, we study the evolution of the Kelvin-Helmholtz instability (KHI) analytically and numerically. We derive the dispersion relation for the two-dimensional KHI including viscous dissipation. The resulting expression for the growth rate is then used to estimate the intrinsic viscosity of four numerical schemes depending on code-specific as well as on physical parameters. Our set of numerical schemes includes the Tree-SPH code VINE, an alternative SPH formulation developed by Price (2008), and the finite-volume grid codes FLASH and PLUTO. In the first part, we explicitly demonstrate the effect of dissipation-inhibiting mechanisms such as the Balsara viscosity on the evolution of the KHI. With VINE, increasing density contrasts lead to a continuously increasing suppression of the KHI (with complete suppression from a contrast of 6:1 or higher). The alternative SPH formulation including an artificial thermal conductivity reproduces the analytically expected growth rates up to a density contrast of 10:1. The second part addresses the shear flow evolution with FLASH and PLUTO. Both codes result in a consistent non-viscous evolution (in the equal as well as in the different density case) in agreement with the analytical prediction. The viscous evolution studied with FLASH shows minor deviations from the analytical prediction.

We present Spitzer observations of a sample of 12 starless cores selected to have prominent 24 micron shadows. The Spitzer images show 8 and 24 micron shadows and in some cases 70 micron shadows; these spatially resolved absorption features trace the densest regions of the cores. We have carried out a 12CO (2-1) and 13CO (2-1) mapping survey of these cores with the Heinrich Hertz Telescope (HHT). We use the shadow features to derive optical depth maps. We derive molecular masses for the cores and the surrounding environment; we find that the 24 micron shadow masses are always greater than or equal to the molecular masses derived in the same region, a discrepancy likely caused by CO freeze--out onto dust grains. We combine this sample with two additional cores that we studied previously to bring the total sample to 14 cores. Using a simple Jeans mass criterion we find that ~ 2/3 of the cores selected to have prominent 24 micron shadows are collapsing or near collapse, a result that is supported by millimeter line observations. Of this subset at least half have indications of 70 micron shadows. All cores observed to produce absorption features at 70 micron are close to collapse. We conclude that 24 micron shadows, and even more so the 70 micron ones, are useful markers of cloud cores that are approaching collapse.

Recent models of molecular cloud formation and evolution suggest that such clouds are dynamic and generally exhibit gravitational collapse. We present a simple analytic model of global collapse onto a filament and compare this with our numerical simulations of the flow-driven formation of an isolated molecular cloud to illustrate the supersonic motions and infall ram pressures expected in models of gravity-driven cloud evolution. We apply our results to observations of the Pipe Nebula, an especially suitable object for our purposes as its low star formation activity implies insignifcant perturbations from stellar feedback. We show that our collapsing cloud model can explain the magnitude of the velocity dispersions seen in the $^{13}$CO filamentary structure by Onishi et al. and the ram pressures required by Lada et al. to confine the lower-mass cores in the Pipe nebula. We further conjecture that higher-resolution simulations will show small velocity dispersions in the densest core gas, as observed, but which are infall motions and not supporting turbulence. Our results point out the inevitability of ram pressures as boundary conditions for molecular cloud filaments, and the possibility that especially lower-mass cores still can be accreting mass at significant rates, as suggested by observations.

We discuss the support for the cooling of gas directly at the disk-halo interface as a mechanism to continually fuel a galaxy. This may be an important reservoir as there is not enough cold gas observed in galaxy halos and simulations indicate the existing cold clouds will be rapidly destroyed as they move through the surrounding halo medium. We show possible evidence for a net infall of the WIM layer in the Milky Way, simulation results showing the recooling of warm clouds at the disk-halo interface, and GALFA HI data of small, cold HI clouds that could represent this recooling.

The formation of stars is a key process in astrophysics. Detailed knowledge of the physical mechanisms that govern stellar birth is a prerequisite for understanding the formation and evolution of our galactic home, the Milky Way. A theory of star formation is an essential part of any model for the origin of our solar system and of planets around other stars. Despite this pivotal importance, and despite many decades of research, our understanding of the processes that initiate and regulate star formation is still limited. Stars are born in cold interstellar clouds of molecular hydrogen gas. Star formation in these clouds is governed by the complex interplay between the gravitational attraction in the gas and agents such as turbulence, magnetic fields, radiation and thermal pressure that resist compression. The competition between these processes determines both the locations at which young stars form and how much mass they ultimately accrete. It plays out over many orders of magnitude in space and time, ranging from galactic to stellar scales. In addition, star formation is a highly stochastic process in which rare and hard-to-predict events, such as the formation of very massive stars and the resulting feedback, can play a dominant role in determining the evolution of a star-forming cloud. As a consequence of the wide range of scales and processes that control star formation, analytic models are usually restricted to highly idealized cases. These can yield insight, but the complexity of the problem means that they must be used in concert with large-scale numerical simulations. Here we summarize the state of modern star formation theory and review the recent advances in numerical simulation techniques.

We present two sets of grid-based hydrodynamical simulations of high-velocity clouds (HVCs) traveling through the diffuse, hot Galactic halo. These HI clouds have been suggested to provide fuel for ongoing star formation in the Galactic disk. The first set of models is best described as a wind-tunnel experiment in which the HVC is exposed to a wind of constant density and velocity. In the second set of models we follow the trajectory of the HVC on its way through an isothermal hydrostatic halo towards the disk. Thus, we cover the two extremes of possible HVC trajectories. The resulting cloud morphologies exhibit a pronounced head-tail structure, with a leading dense cold core and a warm diffuse tail. Morphologies and velocity differences between head and tail are consistent with observations. For typical cloud velocities and halo densities, clouds with H{\small{I}} masses $< 10^{4.5}$ M$_\odot$ will lose their H{\small{I}} content within 10 kpc or less. Their remnants may contribute to a population of warm ionized gas clouds in the hot coronal gas, and they may eventually be integrated in the warm ionized Galactic disk. Some of the (still over-dense, but now slow) material might recool, forming intermediate or low velocity clouds close to the Galactic disk. Given our simulation parameters and the limitation set by numerical resolution, we argue that the derived disruption distances are strong upper limits.

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.

We present high resolution simulations on the impact of ionizing radiation of massive O-stars on the surrounding turbulent interstellar medium (ISM). The simulations are performed with the newly developed software iVINE which combines ionization with smoothed particle hydrodynamics (SPH) and gravitational forces. We show that radiation from hot stars penetrates the ISM, efficiently heats cold low density gas and amplifies over-densities seeded by the initial turbulence. The formation of observed pillar-like structures in star forming regions (e.g. in M16) can be explained by this scenario. At the tip of the pillars gravitational collapse can be induced, eventually leading to the formation of low mass stars. Detailed analysis of the evolution of the turbulent spectra shows that UV-radiation of O-stars indeed provides an excellent mechanism to sustain and even drive turbulence in the parental molecular cloud.

The HII complex N44 in the Large Magellanic Cloud (LMC) provides an excellent site to perform a detailed study of star formation in a mild starburst, as it hosts three regions of star formation at different evolutionary stages and it is not as complicated and confusing as the 30 Doradus giant HII region. We have obtained Spitzer Space Telescope observations and complementary ground-based 4m uBVIJK observations of N44 to identify candidate massive young stellar objects (YSOs). We further classify the YSOs into Types I, II, and III, according to their spectral energy distributions (SEDs). In our sample of 60 YSO candidates, ~65% of them are resolved into multiple components or extended sources in high-resolution ground-based images. We have modeled the SEDs of 36 YSOs that appear single or dominant within a group. We find good fits for Types I and I/II YSOs,but Types II and II/III YSOs show deviations between their observed SEDs and models that do not include PAH emission. We have also found that some Type III YSOs have central holes in their disk components. YSO counterparts are found in four ultracompact HII regions and their stellar masses determined from SED model fits agree well with those estimated from the ionization requirements of the HII regions. The distribution of YSOs is compared with those of the underlying stellar population and interstellar gas conditions to illustrate a correlation between the current formation of O-type stars and previous formation of massive stars. Evidence of triggered star formation is also presented.

We present a set of numerical simulations addressing the effects of magnetic field strength and orientation on the flow-driven formation of molecular clouds. Fields perpendicular to the flows sweeping up the cloud can efficiently prevent the formation of massive clouds but permit the build-up of cold, diffuse filaments. Fields aligned with the flows lead to substantial clouds, whose degree of fragmentation and turbulence strongly depends on the background field strength. Adding a random field component leads to a "selection effect" for molecular cloud formation: high column densities are only reached at locations where the field component perpendicular to the flows is vanishing. Searching for signatures of colliding flows should focus on the diffuse, warm gas, since the cold gas phase making up the cloud will have lost the information about the original flow direction because the magnetic fields redistribute the kinetic energy of the inflows.

We present a three-dimensional, fully parallelized, efficient implementation of ionizing UV radiation for smoothed particle hydrodynamics (SPH) including self-gravity. Our method is based on the SPH/tree code VINE. We therefore call it iVINE (for Ionization + VINE). This approach allows detailed high-resolution studies of the effects of ionizing radiation from e.g. young massive stars on their turbulent parental molecular clouds. In this paper we describe the concept and the numerical implementation of the radiative transfer for a plain-parallel geometry and we discuss several test cases demonstrating the efficiency and accuracy of the new method. As a first application, we study the radiatively driven implosion of marginally stable molecular clouds at various distances of a strong UV source and show that they are driven into gravitational collapse. The resulting cores are very compact and dense exactly as it is observed in clustered environments. Our simulations indicate that the time of triggered collapse depends on the distance of the core from the UV source. Clouds closer to the source collapse several $10^5$ years earlier than more distant clouds. This effect can explain the observed age spread in OB associations where stars closer to the source are found to be younger. We discuss possible uncertainties in the observational derivation of shock front velocities due to early stripping of proto-stellar envelopes by ionizing radiation.

We demonstrate that the observationally inferred rapid onset of star formation after parental molecular clouds have assembled can be achieved by flow-driven cloud formation of atomic gas, using our previous three-dimensional numerical simulations. We post-process these simulations to approximate CO formation, which allows us to investigate the times at which CO becomes abundant relative to the onset of cloud collapse. We find that global gravity in a finite cloud has two crucial effects on cloud evolution. (a) Lateral collapse (perpendicular to the flows sweeping up the cloud) leads to rapidly increasing column densities above the accumulation from the one-dimensional flow. This in turn allows fast formation of CO, allowing the molecular cloud to ``appear'' rapidly. (b) Global gravity is required to drive the dense gas to the high pressures necessary to form solar-mass cores, in support of recent analytical models of cloud fragmentation. While the clouds still appear ``supersonically turbulent'', this turbulence is relegated to playing a secondary role, in that it is to some extent a consequence of gravitational forces.

The observed rapid onset of star formation in molecular clouds requires rapid formation of dense fragments which can collapse individually before being overtaken by global gravitationally-driven flows. Many previous investigations have suggested that supersonic turbulence produces the necessary fragmentation, without addressing however the source of this turbulence. Motivated by our previous (numerical) work on the flow-driven formation of molecular clouds, we investigate the expected timescales of the dynamical and thermal instabilities leading to the rapid fragmentation of gas swept up by large-scale flows, and compare them with global gravitational collapse timescales. We identify parameter regimes in gas density, temperature and spatial scale within which a given instability will dominate. We find that dynamical instabilities disrupt large-scale coherent flows via generation of turbulence, while strong thermal fragmentation amplifies the resulting low-amplitude density perturbations, thus leading to small-scale, high-density fragments as seeds for {\em local} gravity to act upon. Global gravity dominates only on the largest scales; large-scale gravitationally-driven flows promote the formation of groups and clusters of stars formed by turbulence, thermal fragmentation, and rapid cooling.

We present numerical simulations of the growth and saturation of the Kelvin-Helmholtz instability in a compressible fluid layer with and without a weak magnetic field. In the absence of a magnetic field, the instability generates a single eddy which flattens the velocity profile, stabilizing it against further perturbations. Adding a weak magnetic field - weak in the sense that it has almost no effect on the linear instability - leads to a complex flow morphology driven by MHD forces and to enhanced broadening of the layer, due to Maxwell stresses. We corroborate earlier studies which showed that magnetic fields destroy the large scale eddy structure through periodic cycles of windup and resistive decay, but we show that the rate of decay decreases with decreasing plasma resistivity, at least within the range of resistivity accessible to our simulations. Magnetization increases the efficiency of momentum transport, and the transport increases with decreasing resistivity.

Computational gas dynamics has become a prominent research field both in astrophysics and cosmology. In the first part of this review we intend to briefly describe several of the numerical methods used in this field, discuss their range of application and present strategies for converting conditionally-stable numerical methods into unconditionally-stable solution procedures. The underlying aim of the conversion is to enhance the robustness and unification of numerical methods and subsequently enlarge their range of applications considerably. In the second part Fabian Heitsch presents and discusses the implementation of a time-explicit MHD Boltzmann solver.

We study numerically the formation of molecular clouds in large-scale colliding flows including self-gravity. The models emphasize the competition between the effects of gravity on global and local scales in an isolated cloud. Global gravity builds up large-scale filaments, while local gravity -- triggered by a combination of strong thermal and dynamical instabilities -- causes cores to form. The dynamical instabilities give rise to a local focusing of the colliding flows, facilitating the rapid formation of massive protostellar cores of a few 100 M$_\odot$. The forming clouds do not reach an equilibrium state, though the motions within the clouds appear comparable to ``virial''. The self-similar core mass distributions derived from models with and without self-gravity indicate that the core mass distribution is set very early on during the cloud formation process, predominantly by a combination of thermal and dynamical instabilities rather than by self-gravity.

We have identified 230 Tycho-2 Spectral Catalog stars that exhibit 8 micron mid-infrared extraphotospheric excesses in the MidCourse Space Experiment (MSX) and Spitzer Space Telescope Galactic Legacy MidPlane Survey Extraordinaire (GLIMPSE) surveys. Of these, 183 are either OB stars earlier than B8 in which the excess plausibly arises from a thermal bremsstrahlung component or evolved stars in which the excess may be explained by an atmospheric dust component. The remaining 47 stars have spectral classifications B8 or later and appear to be main sequence or late pre-main-sequence objects harboring circumstellar disks. Six of the 47 stars exhibit multiple signatures characteristic of pre-main-sequence circumstellar disks, including emission lines, near-infrared K-band excesses, and X-ray emission. Approximately one-third of the remaining 41 sources have emission lines suggesting relative youth. Of the 25 GLIMPSE stars with SST data at >24 microns, 20 also show an excess at 24 microns. Three additional objects have 24 micron upper limits consistent with possible excesses, and two objects have photospheric measurements at 24 microns. Six MSX sources had a measurement at wavelengths >8 microns. We modeled the excesses in 26 stars having two or more measurements in excess of the expected photospheres as single-component blackbodies. We determine probable disk temperatures and fractional infrared luminosities in the range 191 < T < 787 and 3.9x10^-4 < L_IR/L_* < 2.7x10^-1. We estimate a lower limit on the fraction of Tycho-2 Spectral Catalog main-sequence stars having mid-IR, but not near-IR, excesses to be 1.0+-0.3%.

We revisit the analysis of the Non-linear Thin Shell Instability (NTSI) numerically, including magnetic fields. The magnetic tension force is expected to work against the main driver of the NTSI -- namely transverse momentum transport. However, depending on the field strength and orientation, the instability may grow. For fields aligned with the inflow, we find that the NTSI is suppressed only when the Alfv\'en speed surpasses the (supersonic) velocities generated along the collision interface. Even for fields perpendicular to the inflow, which are the most effective at preventing the NTSI from developing, internal structures form within the expanding slab interface, probably leading to fragmentation in the presence of self-gravity or thermal instabilities. High Reynolds numbers result in local turbulence within the perturbed slab, which in turn triggers reconnection and dissipation of the excess magnetic flux. We find that when the magnetic field is initially aligned with the flow, there exists a (weak) correlation between field strength and gas density. However, for transverse fields, this correlation essentially vanishes. In light of these results, our general conclusion is that instabilities are unlikely to be erased unless the magnetic energy in clouds is much larger than the turbulent energy. Finally, while our study is motivated by the scenario of molecular cloud formation in colliding flows, our results span a larger range of applicability, from supernovae shells to colliding stellar winds.

Cold clouds embedded in warm media are very common objects in astrophysics. Their disruption timescale depends strongly on the dynamical configuration. We discuss the evolution of an initially homogeneous cold cloud embedded in warm turbulent gas. Within a couple of dynamical timescales, the filling factor of the cold gas within the original cloud radius drops below 50%. Turbulent diffusivities estimated from the time evolution of radial filling factor profiles are not constant with time. Cold and warm gas are bodily transported by turbulence and mixed. This is only mildly indicated by column density maps. The radiation field within the cloud, however, increases by several orders of magnitudes due to the mixing, with possible consequences for cloud chemistry and evolution within a few dynamical timescales.

Recent observations with the Spitzer Space Telescope show clear evidence that star formation takes place in the surrounding of young massive O-type stars, which are shaping their environment due to their powerful radiation and stellar winds. In this work we investigate the effect of ionising radiation of massive stars on the ambient interstellar medium (ISM): In particular we want to examine whether the UV-radiation of O-type stars can lead to the observed pillar-like structures and can trigger star formation. We developed a new implementation, based on a parallel Smooth Particle Hydrodynamics code (called IVINE), that allows an efficient treatment of the effect of ionising radiation from massive stars on their turbulent gaseous environment. Here we present first results at very high resolution. We show that ionising radiation can trigger the collapse of an otherwise stable molecular cloud. The arising structures resemble observed structures (e.g. the pillars of creation in the Eagle Nebula (M16) or the Horsehead Nebula B33). Including the effect of gravitation we find small regions that can be identified as formation places of individual stars. We conclude that ionising radiation from massive stars alone can trigger substantial star formation in molecular clouds.

Diffuse emission in the mid-infrared shows a wealth of structure, that lends itself to high-resolution structure analysis of the interstellar gas. A large part of the emission comes from polycyclic aromatic hydrocarbons, excited by nearby ultra-violet sources. Can the observed diffuse emission structure be interpreted as column density structure? We discuss this question with the help of a set of model molecular clouds bathed in the radiation field of a nearby O-star. The correlation strength between column density and ``observed'' flux density strongly depends on the absolute volume density range in the region. Shadowing and irradiation effects may completely alter the appearance of an object. Irradiation introduces additional small-scale structure and it can generate structures resembling shells around HII-regions in objects that do not possess any shell-like structures whatsoever. Nevertheless, structural information about the underlying interstellar medium can be retrieved. In the more diffuse regime ($n({HI})\lesssim 100$cm$^{-3}$), flux density maps may be used to trace the 3D density structure of the cloud via density gradients. Thus, while caution definitely is in order, mid-infrared surveys such as GLIMPSE will provide quantitative insight into the turbulent structure of the interstellar medium.

Molecular Cloud Complexes (MCCs) are highly structured and ``turbulent''. Observational evidence suggests that MCCs are dynamically dominated systems, rather than quasi-equilibrium entities. The observed structure is more likely a consequence of the formation process rather than something that is imprinted after the formation of the MCC. Converging flows provide a natural mechanism to generate MCC structure. We present a detailed numerical analysis of this scenario. Our study addresses the evolution of a MCC from its birth in colliding atomic hydrogen flows up until the point when H$_2$ may begin to form. A combination of dynamical and thermal instabilities breaks up coherent flows efficiently, seeding the small-scale non-linear density perturbations necessary for local gravitational collapse and thus allowing (close to) instantaneous star formation. Many observed properties of MCCs come as a natural consequence of this formation scenario. Since converging flows are omnipresent in the ISM, we discuss the general applicability of this mechanism, from local star formation regions to galaxy mergers.

Molecular clouds (MCs) are highly structured and ``turbulent''. Colliding gas streams of atomic hydrogen have been suggested as a possible source of MCs, imprinting the filamentary structure as a consequence of dynamical and thermal instabilities. We present a 2D numerical analysis of MC formation via converging HI flows. Even with modest flow speeds and completely uniform inflows, non-linear density perturbations as possible precursors of MCs arise. Thus, we suggest that MCs are inevitably formed with substantial structure, e.g., strong density and velocity fluctuations, which provide the initial conditions for subsequent gravitational collapse and star formation in a variety of galactic and extragalactic environments.

Spitzer IRAC 3.6-8 micron photometry obtained as part of the GLIMPSE survey has revealed mid-infrared excesses for 33 field stars with known spectral types in a 1.2 sq. degree field centered on the southern Galactic HII region RCW49. These stars comprise a subset of 184 stars with known spectral classification, most of which were pre-selected to have unusually red IR colors. We propose that the mid-IR excesses are caused by circumstellar dust disks that are either very late remnants of stellar formation or debris disks generated by planet formation. Of these 33 stars, 29 appear to be main-sequence stars based on optical spectral classifications. Five of the 29 main-sequence stars are O or B stars with excesses that can be plausibly explained by thermal bremsstrahlung emission, and four are post main-sequence stars. The lone O star is an O4V((f)) at a spectrophotometric distance of 3233+ 540- 535 pc and may be the earliest member of the Westerlund 2 cluster. Of the remaining 24 main-sequence stars, 18 have SEDs that are consistent with hot dusty debris disks, a possible signature of planet formation. Modeling the excesses as blackbodies demonstrates that the blackbody components have fractional bolometric disk-to-star luminosity ratios, LIR/L*, ranging from 10^-3 to 10^-2 with temperatures ranging from 220 to 820 K. The inferred temperatures are more consistent with asteroid belts rather than the cooler temperatures expected for Kuiper belts. Mid-IR excesses are found in all spectral types from late B to early K.

Spitzer Space Telescope imaging as part of the Galactic Legacy Mid-Plane Survey Extraordinaire (GLIMPSE) reveals a previously unidentified low-latitude rich star cluster near l=31.3 degrees, b=-0.1 degrees. Near-infrared JHK' photometry from the Wyoming Infrared Observatory indicates an extinction of A_V ~ 15+/-3 mag for cluster members. Analysis of 13CO features along the same sightline suggests a probable kinematic distance of 3.1 -- 5.2 kpc. The new cluster has an angular diameter of ~1-2 pc, a total magnitude m_{K_0}=2.1 corrected for extinction, and a luminosity of M_K ~ -10.3 at 3.1 kpc. In contrast to young massive Galactic clusters with ages <100 Myr, the new cluster has no significant radio emission. Comparison to theoretical K-band luminosity functions indicates an age of at least several Gyr and a mass of at least 10^5 solar masses. Unlike known old open clusters, this new cluster lies in the inner Galaxy at R_{GC} ~ 6.1 kpc. We designate this object ``GLIMPSE-C01'' and classify it as a Milky Way globular cluster passing through the Galactic disk. We also identify a region of star formation and fan-shaped outflows from young stellar objects in the same field as the cluster. The cluster's passage through the Galactic molecular layer may have triggered this star formation activity.

The currently most viable methods to estimate magnetic field strengths in molecular cloud cores are Zeeman measurements and the Chandrasekhar-Fermi (CF) method. The CF-method estimates magnetic field strengths from polarimetry and relies on equipartition between the turbulent kinetic and turbulent magnetic energy in the observed region. Thus, its application to objects not dominated by hydromagnetic turbulence is questionable. We calibrate the CF-method and Zeeman measurements against numerical models of self-gravitating molecular cloud cores. We find -- in agreement with previous results -- that the CF-method is on the average accurate up to a factor of 2, while a single estimate can be off by a factor of 7. The CF-method is surprisingly robust even when applied to regions not dominated by hydromagnetic turbulence and in the face of small-number statistics. Zeeman measurements can systematically overestimate the field strength in high-density regions, if the field is predominantly oriented along the line of sight, while they tend to underestimate the field in other cases.

GLIMPSE imaging using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope indicates that star formation is ongoing in the RCW 49 giant H II region. A photometric comparison of the sources in RCW 49 to a similar area to its north finds that at least 300 stars brighter than 13th magnitude in band [3.6] have infrared excesses inconsistent with reddening due to foreground extinction. These are likely young stellar objects (YSOs) more massive than 2.5 M_sun suggesting that thousands more low-mass stars are forming in this cloud. Some of the YSOs are massive (B stars) and therefore very young, suggesting that a new generation of star formation is occurring, possibly triggered by stellar winds and shocks generated by the older (2-3 Myr) central massive cluster. The Spitzer IRAC camera has proven to be ideally suited for distinguishing young stars from field stars, and the GLIMPSE survey of the Galactic Plane will likely find thousands of new star formation regions.

We have calculated the radiation field, dust grain temperatures, and far infrared emissivity of numerical models of turbulent molecular clouds. When compared to a uniform cloud of the same mean optical depth, most of the volume inside the turbulent cloud is brighter, but most of the mass is darker. There is little mean attenuation from center to edge, and clumping causes the radiation field to be somewhat bluer. There is also a large dispersion, typically by a few orders of magnitude, of all quantities relative to their means. However, despite the scatter, the 850 micron emission maps are well correlated with surface density. The fraction of mass as a function of intensity can be reproduced by a simple hierarchical model of density structure.

Multi-frequency radio polarimetric observations of the diffuse Galactic synchrotron background enable us to study the structure of the diffuse ionized gas via rotation measure maps. However, depolarization will introduce artifacts in the resulting rotation measure, most notably in the form of narrow, elongated ``depolarization canals''. We use numerical models of a non-emitting Faraday rotating medium to study the RM distribution needed to create depolarization canals by depolarization due to a finite beam width, and to estimate the influence of this depolarization mechanism on the determination of RM. We argue that the depolarization canals indeed can be caused by beam depolarization, which in turn is a natural consequence when observing a turbulent medium with limited resolution. Furthermore, we estimate that beam depolarization can induce an additional error of about 20% in RM determinations, and considerably less in regions that are not affected by depolarization canals.

Ambient interstellar material may become entrained in outflows from massive stars as a result of shear flow instabilities. We study the linear theory of the Kelvin - Helmholtz instability, the simplest example of shear flow instability, in a partially ionized medium. We model the interaction as a two fluid system (charged and neutral) in a planar geometry. Our principal result is that for much of the relevant parameter space, neutrals and ions are sufficiently decoupled that the neutrals are unstable while the ions are held in place by the magnetic field. Thus, we predict that there should be a detectably narrower line profile in ionized species tracing the outflow compared with neutral species since ionized species are not participating in the turbulent interface with the ambient ISM. Since the magnetic field is frozen to the plasma, it is not tangled by the turbulence in the boundary layer.

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.

Multi-frequency radio polarimetry of the diffuse Galactic synchrotron background gives new viewpoints on the Galactic magnetic field. Rotation measure maps reveal magnetic structures on arcminute to degree scales, such as a ring in polarization that we interpret as a magnetic tunnel. A complication using this technique is depolarization across the beam and along the line of sight. The influence of beam depolarization has been estimated using numerical models of the magneto-ionic ISM, through which polarized radiation propagates. The models show that depolarization canals similar to those observed can be caused by beam depolarization, and that the one-dimensional gradients in RM needed to produce these canals are ubiquitous in the medium.

Under ideal MHD conditions the magnetic field strength should be correlated with density in the interstellar medium (ISM). However, observations indicate that this correlation is weak. Ambipolar diffusion can decrease the flux-to-mass ratio in weakly ionized media; however, it is generally thought to be too slow to play a significant role in the ISM except in the densest molecular clouds. Turbulence is often invoked in astrophysical problems to increase transport rates above the (very slow) laminar values predicted by kinetic theory. We describe a series of numerical experiments addressing the problem of turbulent transport of magnetic fields in weakly ionized gases. We show, subject to various geometrical and physical restrictions, that turbulence in a weakly ionized medium rapidly diffuses the magnetic flux to mass ratio through the buildup of appreciable ion-neutral drifts on small scales. These results are applicable to the fieldstrength - density correlation in the ISM, as well as the merging of flux systems such as protostar and accretion disk fields or protostellar jets with ambient matter, and the vertical transport of galactic magnetic fields.

Magnetic reconnection plays an essential role in the generation and evolution of astrophysical magnetic fields. The best tested and most robust reconnection theory is that of Parker and Sweet. According to this theory, the reconnection rate scales with magnetic diffusivity lambda as lambda^0.5. In the interstellar medium, the Parker-Sweet reconnection rate is far too slow to be of interest. Thus, a mechanism for fast reconnection seems to be required. We have studied the magnetic merging of two oppositely directed flux systems in weakly ionized, but highly conducting, compressible gas. In such systems, ambipolar diffusion steepens the magnetic profile, leading to a thin current sheet. If the ion pressure is small enough, and the recombination of ions is fast enough, the resulting rate of magnetic merging is fast, and independent of lambda. Slow recombination or sufficiently large ion pressure leads to slower merging which scales with lambda as lambda^0.5. We derive a criterion for distinguishing these two regimes, and discuss applications to the weakly ionized ISM and to protoplanetary accretion disks.

The generation and evolution of astrophysical magnetic fields occurs largely through the action of turbulence. In many situations, the magnetic field is strong enough to influence many important properties of turbulence itself. Numerical simulation of magnetized turbulence is especially challenging in the astrophysical regime because of the high magnetic Reynolds numbers involved, but some aspects of this difficulty can be avoided in weakly ionized systems.

To study the interaction of star-formation and turbulent molecular cloud structuring, we analyse numerical models and observations of self-gravitating clouds using the Delta-variance as statistical measure for structural characteristics. In the models we resolve the transition from purely hydrodynamic turbulence to gravitational collapse associated with the formation and mass growth of protostellar cores. We compare models of driven and freely decaying turbulence with and without magnetic fields. Self-gravitating supersonic turbulence always produces a density structure that contains most power on the smallest scales provided by collapsed cores as soon as local collapse sets in. This is in contrast to non-self-gravitating hydrodynamic turbulence where the Delta-variance is dominated by large scale structures. To detect this effect in star-forming regions observations have to resolve the high density contrast of protostellar cores with respect to their ambient molecular cloud. Using the 3mm continuum map of a star-forming cluster in Serpens we show that the dust emission traces the full density evolution. On the contrary, the density range accessible by molecular line observations is insufficient for this analysis. Only dust emission and dust extinction observations are able to to determine the structural parameters of star-forming clouds following the density evolution during the gravitational collapse.

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.

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.

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.

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.

We present new (V,V-I)-photometry of the metal-rich globular clusters NGC 5927, 6316, 6342, 6441 and 6760. The clusters show differential reddening up to dE(V-I)=0.32 mag, for which the CMDs are corrected via extinction maps. There are hints of a variation in the extinction law. Two different ways to determine the parameters metallicity, reddening and distance lead to consistent results. The metallicities of the clusters range between -0.7 <= [M/H] <= 0.0 dex and the absolute reddening between 0.43 <= E(V-I) <= 0.76 mag. Taking the differential reddening into account leads to slightly increased distances. From the resulting parameters we conclude that the usual halo-disk-distinction in the system of globular clusters seems questionable.