Most gas in giant molecular clouds is relatively low-density and forms star inefficiently, converting only a small fraction of its mass to stars per dynamical time. However, star formation models generally predict the existence of a threshold density above which the process is efficient and most mass collapses to stars on a dynamical timescale. A number of authors have proposed observational techniques to search for a threshold density above which star formation is efficient, but it is unclear which of these techniques, if any, are reliable. In this paper we use detailed simulations of turbulent, magnetised star-forming clouds, including stellar radiation and outflow feedback, to investigate whether it is possible to recover star formation thresholds using current observational techniques. Using mock observations of the simulations at realistic resolutions, we show that plots of projected star formation efficiency per free-fall time $\epsilon_{\rm ff}$ can detect the presence of a threshold, but that the resolutions typical of current dust emission or absorption surveys are insufficient to determine its value. In contrast, proposed alternative diagnostics based on a change in the slope of the gas surface density versus star formation rate surface density (Kennicutt-Schmidt relation) or on the correlation between young stellar object counts and gas mass as a function of density are ineffective at detecting thresholds even when they are present. The signatures in these diagnostics sometimes taken as indicative of a threshold in observations, which we generally reproduce in our mock observations, do not prove to correspond to real physical features in the 3D gas distribution.

We present a large suite of simulations of the formation of low-mass star clusters. Our simulations include an extensive set of physical processes -- magnetohydrodynamics, radiative transfer, and protostellar outflows -- and span a wide range of virial parameters and magnetic field strengths. Comparing the outcomes of our simulations to observations, we find that simulations remaining close to virial balance throughout their history produce star formation efficiencies and initial mass function (IMF) peaks that are stable in time and in reasonable agreement with observations. Our results indicate that small-scale dissipation effects near the protostellar surface provide a feedback loop for stabilizing the star formation efficiency. This is true regardless of whether the balance is maintained by input of energy from large scale forcing or by strong magnetic fields that inhibit collapse. In contrast, simulations that leave virial balance and undergo runaway collapse form stars too efficiently and produce an IMF that becomes increasingly top-heavy with time. In all cases we find that the competition between magnetic flux advection toward the protostar and outward advection due to magnetic interchange instabilities, and the competition between turbulent amplification and reconnection close to newly-formed protostars renders the local magnetic field structure insensitive to the strength of the large-scale field, ensuring that radiation is always more important than magnetic support in setting the fragmentation scale and thus the IMF peak mass. The statistics of multiple stellar systems are similarly insensitive to variations in the initial conditions and generally agree with observations within the range of statistical uncertainty.

The turbulent environment from which stars form may lead to misalignment between the stellar spin and the remnant protoplanetary disk. By using hydrodynamic and magnetohydrodynamic simulations, we demonstrate that a wide range of stellar obliquities may be produced as a by-product of forming a star within a turbulent environment. We present a simple semi-analytic model that reveals this connection between the turbulent motions and the orientation of a star and its disk. Our results are consistent with the observed obliquity distribution of hot Jupiters. Migration of misaligned hot Jupiters may, therefore, be due to tidal dissipation in the disk, rather than tidal dissipation of the star-planet interaction.

In regions of star formation, protostars and newborn stars accrete mass from their natal clouds. These clouds are threaded by magnetic fields with a strength characterized by the plasma beta---the ratio of thermal and magnetic pressures. Observations show molecular clouds have beta <= 1, so magnetic fields can play a significant role in the accretion process. We have carried out a numerical study of the effect of large-scale magnetic fields on the rate of accretion onto a uniformly moving point particle from a uniform, non-self-gravitating, isothermal gas. We consider gas moving with sonic Mach numbers of up M ~ 45, magnetic fields that are either parallel, perpendicular, or oriented 45 degrees to the flow, and beta as low as 0.01. Our simulations utilize AMR to obtain high spatial resolution where needed; this also allows the simulation boundaries to be far from the accreting object. Additionally, we show our results are independent of our exact prescription for accreting mass in the sink particle. We give simple expressions for the steady-state accretion rate as a function of beta, M, and field orientation. Using typical molecular clouds values of M ~ 5 and beta ~ 0.04 from the literature, our fits suggest a 0.4 M_Sun star accretes ~ 4*10^{-9} M_Sun/year, almost a factor of two less than accretion rates predicted by hydrodynamic models. This disparity grows to orders of magnitude for stronger fields and lower Mach numbers. We discuss the applicability of these accretion rates versus accretion rates expected from gravitational collapse, and when a steady state is possible. This reduction in Mdot increases the time required to form stars in competitive accretion models, making such models less efficient. In numerical codes, our results should enable accurate subgrid models of sink particles accreting from magnetized media.

We present a set of 3-dimensional, radiation-magnetohydrodynamic calculations of the gravitational collapse of massive (300 Msun), star-forming molecular cloud cores. We show that the combined effects of magnetic fields and radiative feedback strongly suppress core fragmentation, leading to the production of single star systems rather than small clusters. We find that the two processes are efficient at suppressing fragmentation in different regimes, with the feedback most effective in the dense, central region and the magnetic field most effective in more diffuse, outer regions. Thus, the combination of the two is much more effective at suppressing fragmentation than either one considered in isolation. Our work suggests that typical massive cores, which have mass-to-flux ratios of about 2 relative to critical, likely form a single star system, but that cores with weaker fields may form a small star cluster. This result helps us understand why the observed relationship between the core mass function and the stellar initial mass function holds even for ~100 Msun cores with many thermal Jeans masses of material. We also demonstrate that a ~40 AU Keplerian disk is able to form in our simulations, despite the braking effect caused by the strong magnetic field.

We have carried out a numerical study of the effect of large scale magnetic fields on the rate of accretion from a uniform, isothermal gas onto a resistive, stationary point mass. Only mass, not magnetic flux, accretes onto the point mass. The simulations for this study avoid complications arising from boundary conditions by keeping the boundaries far from the accreting object. Our simulations leverage adaptive refinement methodology to attain high spatial fidelity close to the accreting object. Our results are particularly relevant to the problem of star formation from a magnetized molecular cloud in which thermal energy is radiated away on time scales much shorter than the dynamical time scale. Contrary to the adiabatic case, our simulations show convergence toward a finite accretion rate in the limit in which the radius of the accreting object vanishes, regardless of magnetic field strength. For very weak magnetic fields, the accretion rate first approaches the Bondi value and then drops by a factor ~ 2 as magnetic flux builds up near the point mass. For strong magnetic fields, the steady-state accretion rate is reduced by a factor ~ 0.2 \beta^{1/2} compared to the Bondi value, where \beta is the ratio of the gas pressure to the magnetic pressure. We give a simple expression for the accretion rate as a function of the magnetic field strength. Approximate analytic results are given in the Appendixes for both time-dependent accretion in the limit of weak magnetic fields and steady-state accretion for the case of strong magnetic fields.

We report the results of a series of AMR radiation-hydrodynamic simulations of the collapse of massive star forming clouds using the ORION code. These simulations are the first to include the feedback effects protostellar outflows, as well as protostellar radiative heating and radiation pressure exerted on the infalling, dusty gas. We find that that outflows evacuate polar cavities of reduced optical depth through the ambient core. These enhance the radiative flux in the poleward direction so that it is 1.7 to 15 times larger than that in the midplane. As a result the radiative heating and outward radiation force exerted on the protostellar disk and infalling cloud gas in the equatorial direction are greatly diminished. The simultaneously reduces the Eddington radiation pressure barrier to high-mass star formation and increases the minimum threshold surface density for radiative heating to suppress fragmentation compared to models that do not include outflows. The strength of both these effects depends on the initial core surface density. Lower surface density cores have longer free-fall times and thus massive stars formed within them undergo more Kelvin contraction as the core collapses, leading to more powerful outflows. Furthermore, in lower surface density clouds the ratio of the time required for the outflow to break out of the core to the core free-fall time is smaller, so that these clouds are consequently influenced by outflows at earlier stages of collapse. As a result, outflow effects are strongest in low surface density cores and weakest in high surface density one. We also find that radiation focusing in the direction of outflow cavities is sufficient to prevent the formation of radiation pressure-supported circumstellar gas bubbles, in contrast to models which neglect protostellar outflow feedback.

We investigate the role discrete clumps embedded in an astrophysical jet play on the jet's morphology and line emission characteristics. By varying clumps' size, density, position, and velocity, we cover a range of parameter space motivated by observations of objects such as the Herbig Haro object HH~34. We here extend the results presented in Yirak et al. 2009, including how analysis of individual observations may lead to spurious sinusoidal variation whose parameters vary widely over time, owing chiefly to interacts between clumps. The goodness of the fits, while poor in all simulations, are best when clump-clump collisions are minimal. Our results indicate that a large velocity dispersion leads to a clump-clump collision-dominated flow which disrupts the jet beam. Finally, we present synthetic emission images of H-$\alpha$ and [SII] and note an excess of [SII] emission along the jet length as compared to observations. This suggests that observed beams undergo earlier processing, if they are present at all.

The fragmentation of star-forming interstellar clouds, and the resulting stellar initial mass function (IMF), is strongly affected by the temperature structure of the collapsing gas. Since radiation feedback from embedded stars can modify this as collapse proceeds, feedback plays an important role in determining the IMF. However, the effects and importance of radiative heating are likely to depend strongly on the surface density of the collapsing clouds, which determines both their effectiveness at trapping radiation and the accretion luminosities of the stars forming within them. In this paper we report a suite of adaptive mesh refinement radiation-hydrodynamic simulations using the ORION code in which we isolate the effect of column density on fragmentation by following the collapse of clouds of varying column density while holding the mass, initial density and velocity structure, and initial virial ratio fixed. We find that radiation does not significantly modify the overall star formation rate or efficiency, but that it suppresses fragmentation more and more as cloud surface densities increase from those typical of low mass star-forming regions like Taurus, through the typical surface density of massive star-forming clouds in the Galaxy, up to conditions found only in super star clusters. In regions of low surface density, fragmentation during collapse leads to the formation of small clusters rather than individual massive star systems, greatly reducing the fraction of the stellar population with masses >~ 10 Msun. Our simulations have important implications for the formation of massive stars and the universality of the IMF.

Numeric convergence studies demonstrate that the evolution of an adiabatic clump is well-captured by roughly 100 cells per clump radius. The presence of radiative cooling, however, imposes limits on the problem due to the removal of thermal energy. Numerical studies which include radiative cooling typically adopt the 100--200 cells per clump radius resolution. In this paper we present the results of a convergence study for radiatively cooling clumps undertaken over a broad range of resolutions, from 12 to 1,536 cells per clump radius, employing adaptive mesh refinement (AMR) in a 2D axisymmetric geometry ("2.5D"). We also provide a fully 3D simulation, at 192 cells per clump radius, which supports our 2.5D results. We find no appreciable self-convergence at ~100 cells per clump radius as small-scale differences owing to increasingly resolving the "cooling length" have global effects. We therefore conclude that self-convergence is an insufficient criterion to apply on its own when addressing the question of sufficient resolution for radiatively cooled shocked clump simulations. We suggest the adoption of alternate criteria to support a statement of sufficient resolution, such as the demonstration of adequate resolution of the cooling layers behind shocks. We discuss an associated refinement criteria for AMR codes.

A description is given of the algorithms implemented in the AstroBEAR adaptive mesh refinement code for ideal magnetohydrodynamics. The code provides several high resolution, shock capturing schemes which are constructed to maintain conserved quantities of the flow in a finite volume sense. Divergence free magnetic field topologies are maintained to machine precision by collating the components of the magnetic field on a cell-interface staggered grid and utilizing the constrained transport approach for integrating the induction equations. The maintenance of magnetic field topologies on adaptive grids is achieved using prolongation and restriction operators which preserve the divergence and curl of the magnetic field across co-located grids of different resolution. The robustness and correctness of the code is demonstrated by comparing the numerical solution of various tests with analytical solutions or previously published numerical solutions obtained by other codes.

Massive stars produce so much light that the radiation pressure they exert on the gas and dust around them is stronger than their gravitational attraction, a condition that has long been expected to prevent them from growing by accretion. We present three-dimensional radiation-hydrodynamic simulations of the collapse of a massive prestellar core and find that radiation pressure does not halt accretion. Instead, gravitational and Rayleigh-Taylor instabilities channel gas onto the star system through non-axisymmetric disks and filaments that self-shield against radiation, while allowing radiation to escape through optically-thin bubbles. Gravitational instabilities cause the disk to fragment and form a massive companion to the primary star. Radiation pressure does not limit stellar masses, but the instabilities that allow accretion to continue lead to small multiple systems.

Herbig-Haro (HH) jets are commonly thought of as homogeneous beams of plasma traveling at hypersonic velocities. Structure within jet beams is often attributed to periodic or ``pulsed'' variations of conditions at the jet source. Simulations based on this scenario result in knots extending across the jet diameter. Observations and recent high energy density laboratory experiments shed new light on structures below this scale and indicate they may be important for understanding the fundamentals of jet dynamics. In this paper we offer an alternative to ``pulsed'' models of protostellar jets. Using direct numerical simulations we explore the possibility that jets are chains of sub-radial clumps propagating through a moving inter-clump medium. Our models explore an idealization of this scenario by injecting small ($r<r_{jet}$), dense ($\rho>\rho_{jet}$) spheres embedded in an otherwise smooth inter-clump jet flow. The spheres are initialized with velocities differing from the jet velocity by $\sim15$%. We find the consequences of shifting from homogeneous to heterogeneous flows are significant as clumps interact with each other and with the inter-clump medium in a variety of ways. Structures which mimic what is expected from pulsed-jet models can form, as can previously unseen ``sub-radial'' behaviors including backward facing bow shocks and off-axis working surfaces. While these small-scale structures have not been seen before in simulation studies, they are found in high resolution jet observations. We discuss implications of our simulations for the interpretation of protostellar jets with regard to characterization of knots by a ``lifetime'' or ``velocity history'' approach as well as linking observed structures with central engines which produce the jets.

The link between turbulence in star formatting environments and protostellar jets remains controversial. To explore issues of turbulence and fossil cavities driven by young stellar outflows we present a series of numerical simulations tracking the evolution of transient protostellar jets driven into a turbulent medium. Our simulations show both the effect of turbulence on outflow structures and, conversely, the effect of outflows on the ambient turbulence. We demonstrate how turbulence will lead to strong modifications in jet morphology. More importantly, we demonstrate that individual transient outflows have the capacity to re-energize decaying turbulence. Our simulations support a scenario in which the directed energy/momentum associated with cavities is randomized as the cavities are disrupted by dynamical instabilities seeded by the ambient turbulence. Consideration of the energy power spectra of the simulations reveals that the disruption of the cavities powers an energy cascade consistent with Burgers'-type turbulence and produces a driving scale-length associated with the cavity propagation length. We conclude that fossil cavities interacting either with a turbulent medium or with other cavities have the capacity to sustain or create turbulent flows in star forming environments. In the last section we contrast our work and its conclusions with previous studies which claim that jets can not be the source of turbulence.

In this paper we explore the relationship between protostellar outflows and turbulence in molecular clouds. Using 3-D numerical simulations we focus on the hydrodynamics of multiple outflows interacting within a parsec scale volume. We explore the extent to which transient outflows injecting directed energy and momentum into a sub-volume of a molecular cloud can be converted into random turbulent motions. We show that turbulence can readily be sustained by these interactions and show that it is possible to broadly characterize an effective driving scale of the outflows. We compare the velocity spectrum obtained in our studies to that of isotropically forced hydrodynamic turbulence finding that in outflow driven turbulence a power law is indeed achieved. However we find a steeper spectrum (beta ~ 3) is obtained in outflow driven turbulence models than in isotropically forced simulations (beta ~ 2). We discuss possible physical mechanisms responsible for these results as well and their implications for turbulence in molecular clouds where outflows will act in concert with other processes such as gravitational collapse.

Many proto-planetary nebulae (PPN) appear as narrow collimated structures sometimes showing multiple, roughly aligned lobes. In addition, many PPN flows have been shown to have short acceleration times. In this paper we explore whether jet or ``bullet'' (a massive clump) models fit the observations of individual collimated lobes adequately by comparing simulations of both radiatively cooled (stable) jets and bullets. We find that the clump model is somewhat favored over jets because (1) it leads to greater collimation of outflows (2) it accounts better and more naturally for ring-like structures observed in the PPN CRL 618, and (3) it is more successful in reproducing the Hubble-flow character of observed kinematics in some PPN. In addition, bullets naturally account for observed multipolar flows, since the likely MHD launch mechanisms required to drive outflows make multiple non-aligned jets unlikely. Thus we argue that PPN outflows may be driven by explosive MHD launch mechanisms such as those discussed in the context of supernovae (SNe) and gamma-ray bursts(GRB).

We have performed 2D hydrodynamic simulations of a pulsed astrophysical jet propagating through a medium that is populated with spherical inhomogeneities, or "clumps," which are smaller than the jet radius. The clumps are seen to affect the jet in several ways, such as impeding jet propagation and deflecting the jet off-axis. While there has been some debate as to the prevalence of these types of condensations in the ISM or in molecular clouds, the exploration of this region of parameter space nonetheless both shows the potential for these clumps to disrupt astrophysical jets and yields results which recover aspects of recent observations of Herbig-Haro objects. We find that the propagation of the jet and the vorticity induced in the clump/ambient medium correlate well with a "dynamic filling function" $f_d$ across all the simulations.

We investigate the evolution of fossil cavities produced by extinct YSO jets and wide angle outflows. Fossil cavities are ellipsoidal or cylindrical shells of swept-up of ambient (molecular cloud) material moving at low velocities. The cavities form when the momentum in a YSO jet or wide angle outflow decays in time allowing the bowshock or swept-up shell to decelerate to velocities near the turbulent speed in the cloud. It has been suggested in previous studies that cavities provide efficient coupling between the jets/outflows and the cloud and, as such, are the agents by which cloud turbulence is can be re-energized. In this paper we carry forward a series of numerical simulations of jets and outflows whose momentum flux decrease in time. We compare simulations with decaying momentum fluxes to those with constant flux. We show that decaying flux models exhibit deceleration of the outflow head and backfilling via expansion off of the cavity walls. They also have lower density contrast, are longer lived and wider than their continuously driven counterparts. The simulations recover the basic properties of observed fossil cavities. In addition, we provide synthetic observations in terms of P-V diagrams which demonstrate that fossil cavities form both jets and wide angle outflows are characterized by linear "Hubble-law" expansions patterns superimposed on "spur" patterns indicative of the head of a bow shock.

We present a series of numerical studies of the interaction of colliding radiative, hydrodynamic young stellar outflows. We study the effect of the collision impact parameter on the acceleration of ambient material and the degree to which the flow is isotropized by the collision as a mechanism for driving turbulence in the parent molecular cloud. Our results indicate that the high degrees of compression of outflow material, achieved through radiative shocks near the vertex of the interaction, prevents the redirected outflow from spraying over a large spatial region. Furthermore, the collision reduces the redirected outflow's ability to entrain and impart momentum into the ambient cloud. Consideration of the probabilities of outflow collisions leads us to conclude that individual low velocity fossil outflows are the principle coupling between outflows and the cloud.

Outflows from young stellar objects have been identified as a possible source of turbulence in molecular clouds. To investigate the relationship between outflows, cloud dynamics and turbulence, we compare the kinematics of the molecular gas associated with NGC 1333, traced in 13CO(1-0), with the distribution of young stellar objects (YSOs) within. We find a velocity dispersion of ~ 1-1.6 km/s in 13CO that does not significantly vary across the cloud, and is uncorrelated with the number of nearby young stellar outflows identified from optical and submillimeter observations. However, from velocity channel maps we identify about 20 cavities or depressions in the 13CO intensity of scales > 0.1-0.2 pc and velocity widths 1-3 km/s. The cavities exhibit limb brightened rims in both individual velocity channel maps and position velocity diagrams, suggesting that they are slowly expanding. We interpret these cavities to be remnants of past YSO outflow activity: If these cavities are presently empty, they would fill in on time scales of a million years. This can exceed the lifetime of a YSO outflow phase, or the transit time of the central star through the cavity, explaining the the absence of any clear correlation between the cavities and YSO outflows. We find that the momentum and energy deposition associated with the expansion of the cavities is sufficient to power the turbulence in the cloud. In this way we conclude that the cavities are an important intermediary step between the conversion of YSO outflow energy and momentum into cloud turbulent motions.

We carry out high resolution simulations of the inner regions of a wide angle wind driven bipolar outflow using an Adaptive Mesh Refinement code. Our code follows H-He gas with molecular, atomic and ionic components and the associated time dependent molecular chemistry and ionization dynamics with radiative cooling. Our simulations explore the nature of the outflow when a spherical wind expands into a rotating, collapsing envelope. We compare with key observational properties of the outflow system of Source I in the BN/KL region. Our calculations show that the wind evacuates a bipolar outflow cavity in the infalling envelope. We find the head of the outflow to be unstable and that it rapidly fragments into clumps. We resolve the dynamics of the strong shear layer which defines the side walls of the cavity. We conjecture that this layer is the likely site of maser emission and examine its morphology and rotational properties. The shell of swept up ambient gas that delineates the cavity edge retains its angular momentum. This rotation is roughly consistent with that observed in the Source I SiO maser spots. The observed proper motions and line-of-sight velocity are approximately reproduced by the model. The cavity shell at the base of the flow assumes an X-shaped morphology which is also consistent with Source I. We conclude that the wide opening angle of the outflow is evidence that a wide-angle wind drives the Source I outflow and not a collimated jet.

We present two dimensional cylindrically symmetric hydrodynamic simulations and synthetic emission maps of a stellar wind propagating into an infalling, rotating environment. The resulting outflow morphology, collimation and stability observed in these simulations have relevance to the study of young stellar objects, Herbig-Haro jets and molecular outflows. Our code follows hydrogen gas with molecular, atomic and ionic components tracking the associated time dependent molecular chemistry and ionization dynamics with radiative cooling appropriate for a dense molecular gas. We present tests of the code as well as new simulations which indicate the presence of instabilities in the wind-blown bubble's swept-up shell.

We present new data from High Energy Density (HED) laboratory experiments designed to explore the interaction of a heavy hypersonic radiative jet with a cross wind. The jets are generated with the MAGPIE pulsed power machine where converging conical plasma flows are produced from a cylindrically symmetric array of inclined wires. Radiative hypersonic jets emerge from the convergence point. The cross wind is generated by ablation of a plastic foil via soft-X-rays from the plasma convergence region. Our experiments show that the jets are deflected by the action of the cross wind with the angle of deflection dependent on the proximity of the foil. Shocks within the jet beam are apparent in the data. Analysis of the data shows that the interaction of the jet and cross wind is collisional and therefore in the hydro-dynamic regime. MHD plasma code simulations of the experiments are able to recover the deflection behaviour seen in the experiments. We consider the astrophysical relevance of these experiments applying published models of jet deflection developed for AGN and YSOs. Fitting the observed jet deflections to quadratic trajectories predicted by these models allows us to recover a set of plasma parameters consistent with the data. We also present results of 3-D numerical simulations of jet deflection using a new astrophysical Adaptive Mesh Refinement code. These simulations show highly structured shocks occurring within the beam similar to what was observed in the experiments