Recent ALMA observations indicate that while a range of disk sizes exist, typical disk radii are small, and that radial dust drift affects the distribution of solids in disks. Here we explore the consequences of these features in planet population synthesis models. A key feature of our model is planet traps - barriers to otherwise rapid type-I migration of forming planets - for which we include the ice line, heat transition, and outer edge of the dead zone. We find that the ice line plays a fundamental role in the formation of warm Jupiters. In particular, the ratio of super Earths to warm Jupiters formed at the ice line depend sensitively on the initial disk radius. Initial gas disk radii of $\sim$50 AU results in the largest super Earth populations, while both larger and smaller disk sizes result in the ice line producing more gas giants near 1 AU. This transition between typical planet class formed at the ice line at various disk radii confirms that planet formation is fundamentally linked to disk properties (in this case, disk size), and is a result that is only seen when dust evolution effects are included in our models. Additionally, we find that including radial dust drift results in the formation of more super Earths between 0.1 - 1 AU, having shorter orbital radii than those produced in models where dust evolution effects are not included.

The role of outflows in the formation of stars and the protostellar disks that generate them is a central question in astrophysics. Outflows are associated with star formation across the entire stellar mass spectrum. In this review, we describe the observational, theoretical, and computational advances on magnetized outflows, and their role in the formation of disks and stars of all masses in turbulent, magnetized clouds. The ability of torques exerted on disks by magnetized winds to efficiently extract and transport disk angular momentum was developed in early theoretical models and confirmed by a variety of numerical simulations. The recent high resolution ALMA observations of disks and outflows now confirm several key aspects of these ideas, e.g. that jets rotate and originate from large regions of their underlying disks. New insights on accretion disk physics show that magneto-rotational instability (MRI) turbulence is strongly damped, leaving magnetized disk winds as the dominant mechanism for transporting disk angular momentum. This has major consequences for star formation, as well as planet formation. Outflows also play an important role in feedback processes particularly in the birth of low mass stars and cluster formation. Despite being almost certainly fundamental to their production and focusing, magnetic fields in outflows in protostellar systems, and even in the disks, are notoriously difficult to measure. Most methods are indirect and lack precision, as for example, when using optical/near-infrared line ratios. Moreover, in those rare cases where direct measurements are possible - where synchrotron radiation is observed, one has to be very careful in interpreting derived values. Here we also explore what is known about magnetic fields from observations, and take a forward look to the time when facilities such as SPIRou and the SKA are in routine operation.

Star formation involves the flow of gas and dust within molecular clouds into protostars and young stellar objects (YSOs) due to gravity. Along the way, these flows are shaped significantly by many other mechanisms, including pressure, turbulent motions, magnetic fields, stellar feedback, jets, and angular momentum. How all these mechanisms interact nonlinearly with each other on various length scales leads to the formation and evolution of substructures within clouds, including filaments, clumps, cores, disks, outflows, the protostars/YSOs themselves, and planets. In this white paper, prepared for the 2020 Long Range Plan panel which will recommend Canada's future directions for astronomy, we describe the observational and theoretical leadership in the star formation field that Canada's vibrant community has demonstrated over the past decade. Drawing from this extensive background, we identify five key questions that must be addressed for further progress to be made in understanding star formation in the next decade. Addressing these questions will improve our understanding of the dynamics of the dense gas and the role of the magnetic field in star formation, the optical properties of the dust used to trace mass and magnetic fields, the sources of variability in star-forming objects on short timescales, and the physical processes that specifically promote the clustering of stars. We further highlight key facilities in which Canada should become involved to continue making progress in this field. Single-dish facilities we recommend include LSST, trans-atmospheric far-infrared telescopes like BLAST-TNG and SPICA, and ground-based telescopes like JCMT, GBT, and CCAT-p. Interferometric facilities we recommend include ALMA, ngVLA, and SKA1.

To understand the role that planet formation history has on the observable atmospheric carbon-to-oxygen ratio (C/O) we have produced a population of astrochemically evolving protoplanetary disks. Based on the parameters used in a pre-computed population of growing planets their combination allows us to trace the molecular abundances of the gas that is being collected into planetary atmospheres. We include atmospheric pollution of incoming (icy) planetesimals as well as the effect of refractory carbon erosion noted to exist in our own solar system. We find that the carbon and oxygen content of Neptune-mass planets are determined primarily through solid accretion and result in more oxygen-rich (by roughly two orders of magnitude) atmospheres than Hot-Jupiters, whose C/O are primarily determined by gas accretion. Generally we find a `main-sequence' between the fraction of planetary mass accreted through solid accretion and the resulting atmospheric C/O - with planets of higher solid accretion fraction having lower C/O. Hot-Jupiters whose atmospheres have been chemically characterized agree well with our population of planets, and our results suggest that Hot Jupiter formation typically begins near the water ice line. Lower mass hot-Neptunes are observed to be much more carbon-rich (with 0.33 $\lesssim$ C/O $\lesssim$ 1) than is found in our models (C/O $\sim 10^{-2}$), and suggest that some form of chemical processing may affect their observed C/O over the few Gyrs between formation and observation. Our population reproduces the general mass-metallicity trend of the solar system and qualitatively reproduces the C/O-metallicity anti-correlation that has been inferred for the population of characterized exoplanetary atmospheres.

We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth.

Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary disks around young stars, the birth cradles of planets. The high spatial resolution has revealed astonishing structures in disks, such as rings, gaps, asymmetries and spiral arms, and the enormous jump in sensitivity has provided the tools for both large, statistically relevant surveys and deep, sensitive molecular line studies. These observations have revolutionized our view of planet formation, disk formation and disk evolution, bringing model simulations and observations closer to the same level of detail, with many contributions from Canadian researchers on theoretical, observational and technological sides. The new results have inevitably led to a range of new questions, which require next generation instruments such as the Next Generation Very Large Array (ngVLA) and large scale optical infrared facilities. In this white paper we will discuss the current transformation in our understanding of planet formation and the next steps and challenges in connecting theory with exoplanet demographics and protoplanetary disk observations for Canadian research.

Magnetic fields are an elemental part of the interstellar medium in galaxies. However, their impact on gas dynamics and star formation in galaxies remains controversial. We use a suite of global magnetohydrodynamical simulations of isolated disc galaxies to study the influence of magnetic fields on the diffuse and dense gas in the discs. We find that the magnetic field acts in multiple ways. Stronger magnetised discs fragment earlier due to the shorter growth time of the Parker instability. Due to the Parker instability in the magnetised discs we also find cold ($T<50\,\mathrm{K}$) and dense ($n\sim10^3-10^4\,\mathrm{cm}^{-3}$) gas several hundred pc above/below the midplane without any form of stellar feedback. In addition, magnetic fields change the fragmentation pattern. While in the hydrodynamical case, the disc breaks up into ring-like structures, magnetised discs show the formation of filamentary entities that extent both in the azimuthal and radial direction. These kpc scale filaments become magnetically (super-)critical very quickly and allow for the rapid formation of massive giant molecular clouds. Our simulations suggest that major differences in the behaviour of star formation - due to a varying magnetisation - in galaxies could arise.

The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions.

This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars.

We investigate the possibility that multiple populations in globular clusters arise as a natural by-product of massive star-cluster formation. We use 3D radiative hydrodynamics simulations for the formation of young massive clusters to track their chemical self-enrichment during their first 5 Myr. These clusters form embedded within filamentary Giant Molecular Clouds by a combination of gas accretion and rapid merging of protoclusters. Chemical enrichment is a dynamic process happening as the young cluster assembles, so that the original (1P) and enriched (2P) subpopulations of stars form almost simultaneously. Here we test two simple and opposite extremes for the injection of enriched material into the intracluster gas: we assume either continuous injection in a way that tracks the star formation rate; or sudden injection by a single instantaneous event. Using helium abundance as a proxy for the enrichment, we find that realistic multiple population features can be reproduced by injecting a total helium mass amounting to a few percent of the cluster's total mass. The differences in individual growth histories can lead to widely differing 1P/2P outcomes. These models suggest that dual or multiple populations can emerge rapidly in massive star clusters undergoing the typical mode of star cluster formation.

HCN is a key ingredient for synthesizing biomolecules such as nucleobases and amino acids. We calculate 42 reaction rate coefficients directly involved with or in competition with the production of HCN in the early Earth or Titan atmospheres. These reactions are driven by methane and nitrogen radicals produced via UV photodissociation or lightning. For every reaction in this network, we calculate rate coefficients at 298 K using canonical variational transition state theory (CVT) paired with computational quantum chemistry simulations at the BHandHLYP/augcc-pVDZ level of theory. We also calculate the temperature dependence of the rate coefficients for the reactions that have barriers from 50 to 400 K. We present 15 new reaction rate coefficients with no previously known value; 93% of our calculated coefficients are within an order of magnitude of the nearest experimental or recommended values. Above 320 K, the rate coefficient for the new reaction H2CN -> HCN + H dominates. Contrary to experiments, we find the HCN reaction pathway, N + CH3 -> HCN + H2, to be inefficient and suggest that the experimental rate coefficient actually corresponds to an indirect pathway, through the H2CN intermediate. We present CVT using energies computed with density functional theory as a feasible and accurate method for calculating a large network of rate coefficients of small-molecule reactions.

We explore planet formation in the HL Tau disk and possible origins of the prominent gaps and rings observed by ALMA. We investigate whether dust gaps are caused by dynamically trapped planetary embryos at the ice lines of abundant volatiles. The global properties of the HL Tau disk (total mass, size) at its current age are used to constrain an evolving analytic disk model describing its temperature and density profiles. By performing a detailed analysis of the planet-disk interaction for a planet near the water ice line including a rigorous treatment of the dust opacity, we confirm that water is sufficiently abundant (1.5x10^-4 molecules per H) to trap planets at its ice line due to an opacity transition. When the abundance of water is reduced by 50% planet trapping disappears. We extend our analysis to other planet traps: the heat transition, dead zone edge, and the CO_2 ice line and find similar trapping. The formation of planets via planetesimal accretion is computed for dynamically trapped embryos at the water ice line, dead zone, and heat transition. The end products orbit in the inner disk (R < 3 AU), unresolved by ALMA, with masses that range between sub-Earth to 5 Jupiter masses. While we find that the dust gaps correspond well with the radial positions of the CO_2 , CH_4 , and CO ice lines, the planetesimal accretion rates at these radii are too small to build large embryos within 1 Myr.

The origin of fatty acids on the prebiotic Earth is important as they likely formed the encapsulating membranes of the first protocells. Carbon-rich meteorites (i.e., carbonaceous chondrites) such as Murchison and Tagish Lake are well known to contain these molecules, and their delivery to the early planet by intense early meteorite bombardments constitutes a key prebiotic source. We collect the fatty acid abundances measured in various carbonaceous chondrites from the literature and analyze them for patterns and correlations. Fatty acids in meteorites include straight-chain and branched-chain monocarboxylic and dicarboxylic acids up to 12 carbons in length---fatty acids with at least 8 carbons are required to form vesicles, and modern cell membranes employ lipids with ~12--20 carbons. To understand the origin of meteoritic fatty acids, we search the literature for abiotic fatty acid reaction pathways and create a candidate list of 11 reactions that could potentially produce these fatty acids in meteorite parent bodies. Straight-chain monocarboxylic acids (SCMA) are the dominant fatty acids in meteorites, followed by branched-chain monocarboxylic acids (BCMA). SCMA are most abundant in CM2 and Tagish Lake (ungrouped) meteorites, ranging on average from 10$^2$ ppb to 4x10$^5$ ppb, and 10$^4$ ppb to 5x10$^6$ ppb, respectively. In CM, CV, and Tagish Lake meteorites, SCMA abundances generally decrease with increasing carbon chain length. Conversely, SCMA abundances in CR meteorites peak at 5 and 6 carbons in length, and decrease on either side of this peak. This unique CR fatty acid distribution may hint at terrestrial contamination, or that certain fatty acid reactions mechanisms are active in different meteorite parent bodies (planetesimals). We identify Fischer-Tropsch-type synthesis as the most promising pathway for further analysis in the production of fatty acids in planetesimals.

Estimates of the time at which life arose on Earth make use of two types of evidence. First, astrophysical and geophysical studies provide a timescale for the formation of Earth and the Moon, for large impact events on early Earth, and for the cooling of the early magma ocean. From this evidence, we can deduce a habitability boundary, which is the earliest point at which Earth became habitable. Second, biosignatures in geological samples, including microfossils, stromatolites, and chemical isotope ratios, provide evidence for when life was actually present. From these observations we can deduce a biosignature boundary, which is the earliest point at which there is clear evidence that life existed. Studies with molecular phylogenetics and records of the changing level of oxygen in the atmosphere give additional information that helps to determine the biosignature boundary. Here, we review the data from a wide range of disciplines to summarize current information on the timings of these two boundaries. The habitability boundary could be as early as 4.5 Ga, the earliest possible estimate of the time at which Earth had a stable crust and hydrosphere, or as late as 3.9 Ga, the end of the period of heavy meteorite bombardment. The lack of consensus on whether there was a late heavy meteorite bombardment that was significant enough to prevent life is the largest uncertainty in estimating the time of the habitability boundary. The biosignature boundary is more closely constrained. Evidence from carbon isotope ratios and stromatolite fossils both point to a time close to 3.7 Ga. Life must have emerged in the interval between these two boundaries. The time taken for life to appear could, therefore, be within 200 Myr or as long as 800 Myr.

Young massive star clusters (YMCs, with M $\geq$10$^4$ M$_{\odot}$) are proposed modern-day analogues of the globular clusters (GCs) that were products of extreme star formation in the early universe. The exact conditions and mechanisms under which YMCs form remain unknown -- a fact further complicated by the extreme radiation fields produced by their numerous massive young stars. Here we show that GC-sized clusters are naturally produced in radiation-hydrodynamic simulations of isolated 10$^7$ M$_{\odot}$ Giant Molecular Clouds (GMCs) with properties typical of the local universe, even under the influence of radiative feedback. In all cases, these massive clusters grow to GC-level masses within 5 Myr via a roughly equal combination of filamentary gas accretion and mergers with several less massive clusters. Lowering the heavy-element abundance of the GMC by a factor of 10 reduces the opacity of the gas to radiation and better represents the high-redshift formation conditions of GCs. This results in higher gas accretion leading to a mass increase of the largest cluster by a factor of ~4. When combined with simulations of less massive GMCs (10$^{4-6}$ M$_{\odot}$), a clear relation emerges between the maximum YMC mass and the mass of the host GMC. Our results demonstrate that YMCs, and potentially GCs, are a simple extension of local cluster formation to more massive clouds and do not require suggested exotic formation scenarios.

The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales.

Observations show that galaxies and their interstellar media are pervaded by strong magnetic fields with energies in the diffuse component being at least comparable to the thermal and even as large or larger than the turbulent energy. Such strong magnetic fields prevent the formation of stars because patches of the interstellar medium are magnetically subcritical. Here we present the results from global numerical simulations of strongly magnetised and self-gravitating galactic discs, which show that the buoyancy of the magnetic field due to the Parker instability leads at first to the formation of giant filamentary regions. These filamentary structures become gravitationally unstable and fragment into $\sim10^5 M_{\odot}$ clouds that attract kpc long, coherent filamentary flows that build them into GMCs. Our results thus provide a solution to the long-standing problem of how the transition from sub- to supercritical regions in the interstellar medium proceeds.

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation.

We present a comprehensive body of simulations of the formation of exoplanetary populations that incorporate the role of planet traps in slowing planetary migration. The traps we include in our model are the water ice line, the disk heat transition, and the dead zone outer edge. We reduce our model parameter set to two physical parameters: the opacity of the accreting planetary atmospheres ($\kappa_{\rm{env}}$) and a measure of the efficiency of planetary accretion after gap opening ($f_{\rm{max}}$). We perform planet population synthesis calculations based on the initial observed distributions of host star and disk properties - their disk masses, lifetimes, and stellar metallicities. We find the frequency of giant planet formation scales with disk metallicity, in agreement with the observed Jovian planet frequency-metallicity relation. We consider both X-ray and cosmic ray disk ionization models, whose differing ionization rates lead to different dead zone trap locations. In both cases, Jovian planets form in our model out to 2-3 AU, with a distribution at smaller radii dependent on the disk ionization source and the setting of envelope opacity. We find that low values of $\kappa_{\rm{env}}$ (0.001-0.002 cm$^2$ g$^{-1}$) and X-ray disk ionization are necessary to obtain a separation between hot Jupiters near 0.1 AU, and warm Jupiters outside 0.6 AU, a feature present in the data. Our model also produces a large number of super Earths, but the majority are outside of 2 AU. As our model assumes a constant dust to gas ratio, we suggest that radial dust evolution must be taken into account to reproduce the observed super Earth population.

The rapid advances in observations of the different populations of exoplanets, the characterization of their host stars and the links to the properties of their planetary systems, the detailed studies of protoplanetary disks, and the experimental study of the interiors and composition of the massive planets in our solar system provide a firm basis for the next big question in planet formation theory. How do the elemental and chemical compositions of planets connect with their formation? The answer to this requires that the various pieces of planet formation theory be linked together in an end-to-end picture that is capable of addressing these large data sets. In this review, we discuss the critical elements of such a picture and how they affect the chemical and elemental make up of forming planets. Important issues here include the initial state of forming and evolving disks, chemical and dust processes within them, the migration of planets and the importance of planet traps, the nature of angular momentum transport processes involving turbulence and/or MHD disk winds, planet formation theory, and advanced treatments of disk astrochemistry. All of these issues affect, and are affected by the chemistry of disks which is driven by X-ray ionization of the host stars. We discuss how these processes lead to a coherent end-to-end model and how this may address the basic question.

The escape of ultraviolet photons from the densest regions of the interstellar medium (ISM) --- Giant Molecular Clouds (GMCs) --- is a poorly constrained parameter which is vital to understanding the ionization of the ISM and the intergalactic medium. We characterize the escape fraction, f$_{\text{esc,GMC}}$, from a suite of individual GMC simulations with masses in the range 10$^{4-6}$ M$_{\odot}$ using the adaptive-mesh refinement code FLASH. We find significantly different f$_{\text{esc,GMC}}$ depending on the GMC mass which can reach $>$90% in the evolution of 5$\times$10$^4$ and 10$^{5}$ M$_{\odot}$ clouds or remain low at $\sim$5% for most of the lifetime of more massive GMCs. All clouds show fluctuations over short, sub-Myr timescales produced by flickering HII regions. We combine our results to calculate the total escape fraction (f$_{\text{esc,tot}}$) from GMC populations in dwarf starburst and spiral galaxies by randomly drawing clouds from a GMC mass distribution (dN/dM$\propto$M$^{\alpha}$, where $\alpha$ is either -1.5 or -2.5) over fixed time intervals. We find typical f$_{\text{esc,tot}}$ values of 8% for both the dwarf and spiral models. The fluctuations of f$_{\text{esc,tot}}$, however, are much larger for the dwarf models with values as high as 90%. The photons escaping from the 5$\times$10$^4$ and 10$^{5}$ M$_{\odot}$ GMCs are the dominant contributors to f$_{\text{esc,tot}}$ in all cases. We also show that the accompanying star formation rates (SFRs) of our model ($\sim$2$\times$10$^{-2}$ and 0.73 M$_{\odot}$yr$^{-1}$) are consistent with observations of SFRs in dwarf starburst and spiral galaxies, respectively.

Prior to the origin of simple cellular life, the building blocks of RNA (nucleotides) had to form and polymerize in favourable environments on the early Earth. At this time, meteorites and interplanetary dust particles delivered organics such as nucleobases (the characteristic molecules of nucleotides) to warm little ponds whose wet-dry cycles promoted rapid polymerization. We build a comprehensive numerical model for the evolution of nucleobases in warm little ponds leading to the emergence of the first nucleotides and RNA. We couple Earth's early evolution with complex prebiotic chemistry in these environments. We find that RNA polymers must have emerged very quickly after the deposition of meteorites (< a few years). Their constituent nucleobases were primarily meteoritic in origin and not from interplanetary dust particles. Ponds appeared as continents rose out of the early global ocean but this increasing availability of "targets" for meteorites was offset by declining meteorite bombardment rates. Moreover, the rapid losses of nucleobases to pond seepage during wet periods, and to UV photodissociation during dry periods means that the synthesis of nucleotides and their polymerization into RNA occurred in just one to a few wet-dry cycles. Under these conditions, RNA polymers likely appeared prior to 4.17 billion years ago. Significance: There are two competing hypotheses for the site at which an RNA world emerged: hydrothermal vents in the deep ocean and warm little ponds. Because the former lacks wet and dry cycles, which are well known to promote polymerization (in this case, of nucleotides into RNA), we construct a comprehensive model for the origin of RNA in the latter sites. Our model advances the story and timeline of the RNA world by constraining the source of biomolecules, the environmental conditions, the timescales of reaction, and the emergence of first RNA polymers.

The process of radiative feedback in Giant Molecular Clouds (GMCs) is an important mechanism for limiting star cluster formation through the heating and ionization of the surrounding gas. We explore the degree to which radiative feedback affects early ($\lesssim$5 Myr) cluster formation in GMCs having masses that range from 10$^{4-6}$ M$_{\odot}$ using the FLASH code. The inclusion of radiative feedback lowers the efficiency of cluster formation by 20-50\% relative to hydrodynamic simulations. Two models in particular --- 5$\times$10$^4$ and 10$^5$ M$_{\odot}$ --- show the largest suppression of the cluster formation efficiency, corresponding to a factor of $\sim$2. For these clouds only, the internal energy, a measure of the energy injected by radiative feedback, exceeds the gravitational potential for a significant amount of time. We find a clear relation between the maximum cluster mass, M$_{cl,max}$, formed in a GMC of mass M$_{GMC}$; M$_{cl,max}\propto$ M$_{GMC}^{0.81}$. This scaling result suggests that young globular clusters at the necessary scale of $10^6 M_{\odot}$ form within host GMCs of masses near $\sim 5 \times 10^7 M_{\odot}$. We compare simulated cluster mass distributions to the observed embedded cluster mass function ($dlog(N)/dlog(M) \propto M^{\beta}$ where $\beta$ = -1) and find good agreement ($\beta$ = -0.99$\pm$0.14) only for simulations including radiative feedback, indicating this process is important in controlling the growth of young clusters. However, the high star formation efficiencies, which range from 16-21\%, and high star formation rates compared to locally observed regions suggest other feedback mechanisms are also important during the formation and growth of stellar clusters.

We present the next step in a series of papers devoted to connecting the composition of the atmospheres of forming planets with the chemistry of their natal evolving protoplanetary disks. The model presented here computes the coupled chemical and dust evolution of the disk and the formation of three planets per disk model. Our three canonical planet traps produce a Jupiter near 1 AU, a Hot Jupiter and a Super-Earth. We study the dependency of the final orbital radius, mass, and atmospheric chemistry of planets forming in disk models with initial disk masses that vary by 0.02 $M_\odot$ above and below our fiducial model ($M_{disk,0} = 0.1 ~M_\odot$). We compute C/O and C/N for the atmospheres formed in our 3 models and find that C/O$_{\rm planet}\sim$ C/O$_{\rm disk}$, which does not vary strongly between different planets formed in our model. The nitrogen content of atmospheres can vary in planets that grow in different disk models. These differences are related to the formation history of the planet, the time and location that the planet accretes its atmosphere, and are encoded in the bulk abundance of NH$_3$. These results suggest that future observations of atmospheric NH$_3$ and an estimation of the planetary C/O and C/N can inform the formation history of particular planetary systems.

The UV photon escape fraction from molecular clouds is a key parameter for understanding the ionization of the Interstellar Medium (ISM), and extragalactic processes, such as cosmic reionization. We present the ionizing photon flux and the corresponding photon escape fraction (f$_{esc}$) arising as a consequence of star cluster formation in a turbulent, 10$^6$ M$_{\odot}$ GMC, simulated using the code FLASH. We make use of sink particles to represent young, star-forming clusters coupled with a radiative transfer scheme to calculate the emergent UV flux. We find that the ionizing photon flux across the cloud boundary is highly variable in time and space due to the turbulent nature of the intervening gas. The escaping photon fraction remains at $\sim$5% for the first 2.5 Myr, followed by two pronounced peaks at 3.25 and 3.8 Myr with a maximum f$_{esc}$ of 30% and 37%, respectively. These peaks are due to the formation of large HII regions, that expand into regions of lower density and some of which reach the cloud surface. However, these phases are short lived and f$_{esc}$ drops sharply as the HII regions are quenched by the central cluster passing through high-density material due to the turbulent nature of the cloud. We find an average f$_{esc}$ of 15% with factor of two variations over 1 Myr timescales. Our results suggest that assuming a single value for f$_{esc}$ from a molecular cloud is in general a poor approximation, and that the dynamical evolution of the system leads to large temporal variation.

The possible meteorite parent body origin of Earth's pregenetic nucleobases is substantiated by the guanine (G), adenine (A) and uracil (U) measured in various meteorites. Cytosine (C) and thymine (T) however are absent in meteorites, making the emergence of a RNA and later RNA/DNA/protein world problematic. We investigate the meteorite parent body (planetesimal) origin of all nucleobases by computationally modeling 18 reactions that potentially contribute to nucleobase formation in such environments. Out of this list, we identify the two most important reactions for each nucleobase and find that these involve small molecules such as HCN, CO, NH3, and water that ultimately arise from the protoplanetary disks in which planetesimals are built. The primary result of this study is that cytosine is unlikely to persist within meteorite parent bodies due to aqueous deamination. Thymine has a thermodynamically favourable reaction pathway from uracil, formaldehyde and formic acid, but likely did not persist within planetesimals containing H2O2 due to an oxidation reaction with this molecule. Finally, while FT synthesis is found to be the dominant source of nucleobases within our model planetesimal, NC synthesis may still be significant under certain chemical conditions (e.g. within CR2 parent bodies). We discuss several major consequences of our results for the origin of the RNA world.

We have developed a new model for the astrochemical structure of a viscously evolving protoplanetary disk that couples an analytic description of the disk's temperature and density profile, chemical evolution, and an evolving dust distribution. We compute evolving radial distributions for a range of dust grain sizes, which depend on coagulation, fragmentation and radial drift processes. In particular we find that the water ice line plays an important role in shaping the radial distribution of the maximum grain size because ice coated grains are significantly less susceptible to fragmentation than their dry counterparts. This in turn has important effects on disk ionization and therefore on the location of dead zones. In comparison to a simple constant gas-to-dust ratio model for the dust as an example, we find that the new model predicts an outer dead zone edge that moves in by a factor of about 3 at 1 Myr (to 5 AU) and by a factor of about 14 by 3 Myr (to 0.5 AU). We show that the changing position of the dead zone and heat transition traps have important implications for the formation and trapping of planets in protoplanetary disks. Finally, we consider our results in light of recent ALMA observations of HL Tau and TW Hya.

We present an analysis of the relationship between the orientation of magnetic fields and filaments that form in 3D magnetohydrodynamic simulations of cluster-forming, turbulent molecular cloud clumps. We examine simulated cloud clumps with size scales of L ~ 2-4 pc and densities of n ~ 400-1000 cm^-3 with Alfven Mach numbers near unity. We simulated two cloud clumps of different masses, one in virial equilibrium, the other strongly gravitationally bound, but with the same initial turbulent velocity field and similar mass-to-flux ratio. We apply various techniques to analyze the filamentary and magnetic structure of the resulting cloud, including the DisPerSE filament-finding algorithm in 3D. The largest structure that forms is a 1-2 parsec-long filament, with smaller connecting sub-filaments. We find that our simulated clouds, wherein magnetic forces and turbulence are comparable, coherent orientation of the magnetic field depends on the virial parameter. Subvirial clumps undergo strong gravitational collapse and magnetic field lines are dragged with the accretion flow. We see evidence of filament-aligned flow and accretion flow onto the filament in the subvirial cloud. Magnetic fields oriented more parallel in the subvirial cloud and more perpendicular in the denser, marginally bound cloud. Radiative feedback from a 16 Msun star forming in a cluster in one of our simulations ultimately results in the destruction of the main filament, the formation of an HII region, and the sweeping up of magnetic fields within an expanding shell at the edges of the HII region.

Super Earths are the largest population of exoplanets and are seen to exhibit a rich diversity of compositions as inferred through their mean densities. Here we present a model that combines equilibrium chemistry in evolving disks with core accretion that tracks materials accreted onto planets during their formation. In doing so, we aim to explain why super Earths form so frequently and how they acquire such a diverse range of compositions. A key feature of our model is disk inhomogeneities, or planet traps, that act as barriers to rapid type-I migration. The traps we include are the dead zone, which can be caused by either cosmic ray or X-ray ionization, the ice line, and the heat transition. We find that in disks with sufficiently long lifetimes ($\gtrsim$ 4 Myr), all traps produce Jovian planets. In these disks, planet formation in the heat transition and X-ray dead zone produces hot Jupiters while the ice line and cosmic ray dead zones produce Jupiters at roughly 1 AU. Super Earth formation takes place within short-lived disks ($\lesssim$ 2 Myr), whereby the disks are photoevaporated while planets are in a slow phase of gas accretion. We find that super Earth compositions range from dry and rocky ($<$ 6 % ice by mass) to those with substantial water contents ($>$ 30 % ice by mass). The traps play a crucial role in our results, as they dictate where in the disk particular planets can accrete from, and what compositions they are able to acquire.

We present a model of the early chemical composition and elemental abundances of planetary atmospheres based on the cumulative gaseous chemical species that are accreted onto planets forming by core accretion from evolving protoplanetary disks. The astrochemistry of the host disk is computed using an ionization driven, non-equilibrium chemistry network within viscously evolving disk models. We accrete gas giant planets whose orbital evolution is controlled by planet traps using the standard core accretion model and track the chemical composition of the material that is accreted onto the protoplanet. We choose a fiducial disk model and evolve planets in 3 traps - water ice line, dead zone and heat transition. For a disk with a lifetime of 4.1 Myr we produce two Hot Jupiters (M = 1.43, 2.67 M$_{\rm Jupiter}$, r = 0.15, 0.11 AU) in the heat transition and ice line trap and one failed core (M = 0.003 M$_{\rm Jupiter}$, r =3.7 AU) in the dead zone. These planets are found with mixing ratios for CO and H$_2$O of $1.99\times 10^{-4}$, $5.0\times 10^{-4}$ respectively for both Hot Jupiters. Additionally for these planets we find CO$_2$ and CH$_4$, with mixing ratios of $1.8\times 10^{-6}\rightarrow 9.8\times 10^{-10}$ and $1.1\times 10^{-8}\rightarrow 2.3\times 10^{-10}$ respectively. These ranges correspond well with the mixing ratio ranges that have been inferred through the detection of emission spectra from Hot Jupiters by multiple authors. We compute a carbon-to-oxygen ratio of 0.227 for the ice line planet and 0.279 for the heat transition planet. These planets accreted their gas inside the ice line, hence the sub-solar C/O.

Radiative feedback is an important consequence of cluster formation in Giant Molecular Clouds (GMCs) in which newly formed clusters heat and ionize their surrounding gas. The process of cluster formation, and the role of radiative feedback, has not been fully explored in different GMC environments. We present a suite of simulations which explore how the initial gravitational boundedness, and radiative feedback, affect cluster formation. We model the early evolution ($<$ 5 Myr) of turbulent, 10$^6$ M$_{\odot}$ clouds with virial parameters ranging from 0.5 to 5. To model cluster formation, we use cluster sink particles, coupled to a raytracing scheme, and a custom subgrid model which populates a cluster via sampling an IMF with an efficiency of 20\% per freefall time. We find that radiative feedback only decreases the cluster particle formation efficiency by a few percent. The initial virial parameter plays a much stronger role in limiting cluster formation, with a spread of cluster formation efficiencies of 37\% to 71\% for the most unbound to the most bound model. The total number of clusters increases while the maximum mass cluster decreases with an increasing initial virial parameter, resulting in steeper mass distributions. The star formation rates in our cluster particles are initially consistent with observations but rise to higher values at late times. This suggests that radiative feedback alone is not responsible for dispersing a GMC over the first 5 Myr of cluster formation.

Over the past few years, astronomers have detected thousands of planets and candidate planets by observing their periodic transits in front of their host stars. Related methods might soon allow studies of the chemical imprints of life in extrasolar planetary atmospheres. Here, we address the reciprocal question, namely, from where is Earth detectable by extrasolar observers using similar methods. We explore Earth's transit zone (ETZ), the projection of a band around Earth's ecliptic onto the celestial plane, where observers can detect Earth transits across the Sun. The ETZ is between $0.520^\circ$ and $0.537^\circ$ wide due to the non-circular Earth orbit. The restricted ETZ (rETZ), where Earth transits the Sun less than 0.5 solar radii from its center, is ~$0.262^\circ$ wide. We first compile a target list of 45 K and 37 G dwarf stars inside the rETZ and within 1 kiloparsec (3260 lightyears) using the Hipparcos catalog. We then greatly enlarge the number of potential targets by constructing an analytic galactic disk model and find that ~$10^5$ K and G dwarf stars should reside within the rETZ. The ongoing GAIA space mission can potentially discover all G dwarfs among them (several $10^4$) within the next five years. Many more potentially habitable planets orbit dim, unknown M stars in the ETZ and other stars that traversed the ETZ long time ago. If any of these planets host intelligent observers, they could have identified Earth as a habitable, or even as a living, world, and we could be receiving their broadcasts today. The K2 mission, the Allen Telescope Array, the upcoming Square Kilometer Array, or the Green Bank Telescope might detect such deliberate extraterrestrial messages. Ultimately, the ETZ would be an ideal region to monitor by the Breakthrough Listen Initiatives, an upcoming survey that will constitute the most comprehensive search for extraterrestrial intelligence so far.

We present radiation hydrodynamic simulations of collapsing protostellar cores with initial masses of 30, 100, and 200 M$_{\odot}$. We follow their gravitational collapse and the formation of a massive protostar and protostellar accretion disk. We employ a new hybrid radiative feedback method blending raytracing techniques with flux-limited diffusion for a more accurate treatment of the temperature and radiative force. In each case, the disk that forms becomes Toomre-unstable and develops spiral arms. This occurs between 0.35 and 0.55 freefall times and is accompanied by an increase in the accretion rate by a factor of 2-10. Although the disk becomes unstable, no other stars are formed. In the case of our 100 and 200 M$_{\odot}$ simulation, the star becomes highly super-Eddington and begins to drive bipolar outflow cavities that expand outwards. These radiatively-driven bubbles appear stable, and appear to be channeling gas back onto the protostellar accretion disk. Accretion proceeds strongly through the disk. After 81.4 kyr of evolution, our 30 M$_{\odot}$ simulation shows a star with a mass of 5.48 M$_{\odot}$ and a disk of mass 3.3 M$_{\odot}$, while our 100 M$_{\odot}$ simulation forms a 28.8 M$_{\odot}$ mass star with a 15.8 M$_{\odot}$ disk over the course of 41.6 kyr, and our 200 M$_{\odot}$ simulation forms a 43.7 M$_{\odot}$ star with an 18 M$_{\odot}$ disk in 21.9 kyr. In the absence of magnetic fields or other forms of feedback, the masses of the stars in our simulation do not appear limited by their own luminosities.

In the Grand Tack (GT) scenario for the young solar system, Jupiter formed beyond 3.5 AU from the Sun and migrated as close as 1.5 AU until it encountered an orbital resonance with Saturn. Both planets then supposedly migrated outward for several $10^5$ yr, with Jupiter ending up at ~5 AU. The initial conditions of the GT and the timing between Jupiter's migration and the formation of the Galilean satellites remain unexplored. We study the formation of Ganymede and Callisto, both of which consist of ~50% H$_2$O and rock, in the GT scenario. We examine why they lack dense atmospheres, while Titan is surrounded by a thick N$_2$ envelope. We model an axially symmetric circumplanetary disk (CPD) in hydrostatic equilibrium around Jupiter. The CPD is warmed by viscous heating, Jupiter's luminosity, accretional heating, and the Sun. The position of the H$_2$O ice line in the CPD, which is crucial for the formation of massive moons, is computed at various solar distances. We assess the loss of Galilean atmospheres due to high-energy radiation from the young Sun. Ganymede and Callisto cannot have accreted their H$_2$O during Jupiter's supposed GT, because its CPD (if still active) was too warm to host ices and much smaller than Ganymede's contemporary orbit. From a thermal perspective, the Galilean moons might have had significant atmospheres, but these would probably have been eroded during the GT in < $10^5$ yr by solar XUV radiation. Jupiter and the Galilean moons formed beyond 4.5 (+/- 0.5) AU and prior to the proposed GT. Thereafter, Jupiter's CPD would have been dry, and delayed accretion of planetesimals should have created water-rich Io and Europa. While Galilean atmospheres would have been lost during the GT, Titan would have formed after Saturn's own tack, because Saturn still accreted substantially for ~$10^6$ yr after its closest solar approach, ending up at about 7 AU.

Carbonaceous chondrite meteorites are known for having high water and organic material contents, including amino acids. Here we address the origin of amino acids in the warm interiors of their parent bodies (planetesimals) within a few million years of their formation, and connect this with the astrochemistry of their natal protostellar disks. We compute both the total amino acid abundance pattern as well as the relative frequencies of amino acids within the CM2 (e.g. Murchison) and CR2 chondrite subclasses based on Strecker reactions within these bodies. We match the relative frequencies to well within an order of magnitude among both CM2 and CR2 meteorites for parent body temperatures $<$ 200$^{\circ}$C. These temperatures agree with 3D models of young planetesimal interiors. We find theoretical abundances of approximately 7x10$^5$ parts-per-billion (ppb), which is in agreement with the average observed abundance in CR2 meteorites of 4$\pm$7x10$^5$, but an order of magnitude higher than the average observed abundance in CM2 meteorites of 2$\pm$2x10$^4$. We find that the production of hydroxy acids could be favoured over the production of amino acids within certain meteorite parent bodies (e.g. CI1, CM2) but not others (e.g. CR2). This could be due to the relatively lower NH$_3$ abundances within CI1 and CM2 meteorite parent bodies, which leads to less amino acid synthesis. We also find that the water content in planetesimals is likely to be the main cause of variance between carbonaceous chondrites of the same subclass. We propose that amino acid abundances are primarily dependent on the ammonia and water content of planetesimals that are formed in chemically distinct regions within their natal protostellar disks.

Carbonaceous chondrites are a class of meteorite known for having a high content of water and organics. In this study, abundances of the nucleobases, i.e., the building blocks of RNA and DNA, found in carbonaceous chondrites are collated from a variety of published data and compared across various meteorite classes. An extensive review of abiotic chemical reactions producing nucleobases is then performed. These reactions are then reduced to a list of 15 individual reaction pathways that could potentially occur within meteorite parent bodies. The nucleobases guanine, adenine and uracil are found in carbonaceous chondrites in the amounts of 1$-$500 ppb. It is currently unknown which reaction is responsible for their synthesis within the meteorite parent bodies. One class of carbonaceous meteorites dominate the abundances of both amino acids and nucleobases$-$the so-called CM2 (e.g. Murchison meteorite). CR2 meteorites (e.g. Graves Nunataks) also dominate the abundances of amino acids, but are the least abundant in nucleobases. The abundances of total nucleobases in these two classes are $330 \pm250$ ppb and $16 \pm13$ ppb respectively. Guanine most often has the greatest abundances in carbonaceous chondrites with respect to the other nucleobases, but is 1$-$2 orders of magnitude less abundant in CM2 meteorites than glycine (the most abundant amino acid). Our survey of the reaction mechanisms for nucleobase formation suggests that Fischer-Tropsch synthesis (i.e. CO, H$_2$ and NH$_3$ gases reacting in the presence of a catalyst such as alumina or silica) is the most likely for conditions that characterize early states of planetesimals.

Most of the exoplanets with known masses at Earth-like distances to Sun-like stars are heavier than Jupiter, which raises the question of whether such planets are accompanied by detectable, possibly habitable moons. Here we simulate the accretion disks around super-Jovian planets and find that giant moons with masses similar to Mars can form. Our results suggest that the Galilean moons formed during the final stages of accretion onto Jupiter, when the circumjovian disk was sufficiently cool. But in contrast to other studies, with our assumptions, we show that Jupiter was still feeding from the circumsolar disk and that its principal moons cannot have formed after the complete photoevaporation of the circumsolar nebula. To counteract the steady loss of moons into the planet due to type I migration, we propose that the water ice line around Jupiter and super-Jovian exoplanets acted as a migration trap for moons. Heat transitions, however, cross the disk during the gap opening within 10^4 yr, which makes them inefficient as moon traps. This indicates a fundamental difference between planet and moon formation. We find that icy moons larger than the smallest known exoplanet can form at about 15 - 30 Jupiter radii around super-Jovian planets. Their size implies detectability by the Kepler and PLATO space telescopes as well as by the European Extremely Large Telescope.

We present an algorithm for solving the radiative transfer problem on massively parallel computers using adaptive mesh refinement and domain decomposition. The solver is based on the method of characteristics which requires an adaptive raytracer that integrates the equation of radiative transfer. The radiation field is split into local and global components which are handled separately to overcome the non-locality problem. The solver is implemented in the framework of the magneto-hydrodynamics code FLASH and is coupled by an operator splitting step. The goal is the study of radiation in the context of star formation simulations with a focus on early disc formation and evolution. This requires a proper treatment of radiation physics that covers both the optically thin as well as the optically thick regimes and the transition region in particular. We successfully show the accuracy and feasibility of our method in a series of standard radiative transfer problems and two 3D collapse simulations resembling the early stages of protostar and disc formation.

Exomoon detections might be feasible with NASA's Kepler or ESA's upcoming PLATO mission or the ground-based E-ELT. To use observational resources most efficiently we need to know where the largest, most easily detected moons can form. We explore the possibility of large exomoons by following the movement of water (H2O) ice lines in the accretion disks around young super-Jovian planets. We want to know how different heating sources in those disks affect the H2O ice lines. We simulate 2D rotationally symmetric accretion disks in hydrostatic equilibrium around super-Jovian exoplanets. The energy terms in our semi-analytical model -- (1) viscous heating, (2) planetary illumination, (3) accretional heating, and (4) stellar illumination -- are fed by precomputed planet evolution tracks. We consider planets accreting 1 to 12 Jupiter masses at distances between 1 and 20 AU to a Sun-like star. Accretion disks around Jupiter-mass planets closer than ~4.5 AU to Sun-like stars do not feature H2O ice lines, but the most massive super-Jovians can form icy satellites as close as ~3 AU to Sun-like stars. Super-Jovian planets forming beyond ~5 AU can host Mars-mass moons. We study a broad range of disk parameters for planets at 5.2 AU and find that the H2O ice lines are universally between ~15 and 30 Jupiter radii when the last generation of moons is forming. If the abundant population of super-Jovian planets at ~1 AU formed in situ, then they should lack giant icy moons because their disks did not host H2O ice in the final stages of accretion. In the more likely case that these planets migrated to their current locations from beyond a few AU, they might be orbited by large, H2O-rich moons. In this case, Mars-mass ocean moons might be common in the stellar habitable zones. Future exomoon searches can provide powerful constraints on the formation and migration history of giant exoplanets.

We use numerical simulations of turbulent cluster-forming regions to study the nature of dense filamentary structures in star formation. Using four hydrodynamic and magnetohydrodynamic simulations chosen to match observations, we identify filaments in the resulting column density maps and analyze their properties. We calculate the radial column density profiles of the filaments every 0.05 Myr and fit the profiles with the modified isothermal and pressure confined isothermal cylinder models, finding reasonable fits for either model. The filaments formed in the simulations have similar radial column density profiles to those observed. Magnetic fields provide additional pressure support to the filaments, making `puffier' filaments less prone to fragmentation than in the pure hydrodynamic case, which continue to condense at a slower rate. In the higher density simulations, the filaments grow faster through the increased importance of gravity. Not all of the filaments identified in the simulations will evolve to form stars: some expand and disperse. Given these different filament evolutionary paths, the trends in bulk filament width as a function of time, magnetic field strength, or density, are weak, and all cases are reasonably consistent with the finding of a constant filament width in different star-forming regions. In the simulations, the mean FWHM lies between 0.06 and 0.26 pc for all times and initial conditions, with most lying between 0.1 to 0.15 pc; the range in FWHMs are, however, larger than seen in typical Herschel analyses. Finally, the filaments display a wealth of substructure similar to the recent discovery of filament bundles in Taurus.

Stellar feedback, star formation and gravitational interactions are major controlling forces in the evolution of Giant Molecular Clouds (GMCs). To explore their relative roles, we examine the properties and evolution of GMCs forming in an isolated galactic disk simulation that includes both localised thermal feedback and photoelectric heating. The results are compared with the three previous simulations in this series which consists of a model with no star formation, star formation but no form of feedback and star formation with photoelectric heating in a set with steadily increasing physical effects. We find that the addition of localised thermal feedback greatly suppresses star formation but does not destroy the surrounding GMC, giving cloud properties closely resembling the run in which no stellar physics is included. The outflows from the feedback reduce the mass of the cloud but do not destroy it, allowing the cloud to survive its stellar children. This suggests that weak thermal feedback such as the lower bound expected for supernova may play a relatively minor role in the galactic structure of quiescent Milky Way-type galaxies, compared to gravitational interactions and disk shear.

We present the results of large scale, three-dimensional magneto-hydrodynamics simulations of disc-winds for different initial magnetic field configurations. The jets are followed from the source to 90 AU scale, which covers several pixels of HST images of nearby protostellar jets. Our simulations show that jets are heated along their length by many shocks. We compute the emission lines that are produced, and find excellent agreement with observations. The jet width is found to be between 20 and 30 AU while the maximum velocities perpendicular to the jet is found to be up to above 100 km/s. The initially less open magnetic field configuration simulations results in a wider, two-component jet; a cylindrically shaped outer jet surrounding a narrow and much faster, inner jet. These simulations preserve the underlying Keplerian rotation profile of the inner jet to large distances from the source. However, for the initially most open magnetic field configuration the kink mode creates a narrow corkscrew-like jet without a clear Keplerian rotation profile and even regions where we observe rotation opposite to the disc (counter-rotating). The RW Aur jet is narrow, indicating that the disc field in that case is very open meaning the jet can contain a counter-rotating component that we suggests explains why observations of rotation in this jet has given confusing results. Thus magnetized disc winds from underlying Keplerian discs can develop rotation profiles far down the jet that are not Keplerian.

We analyse simulations of turbulent, magnetised molecular cloud cores focussing on the formation of Class 0 stage protostellar discs and the physical conditions in their surroundings. We show that for a wide range of initial conditions Keplerian discs are formed in the Class 0 stage already. In particular, we show that even subsonic turbulent motions reduce the magnetic braking efficiency sufficiently in order to allow rotationally supported discs to form. We therefore suggest that already during the Class 0 stage the fraction of Keplerian discs is significantly higher than 50%, consistent with recent observational trends but significantly higher than predictions based on simulations with misaligned magnetic fields, demonstrating the importance of turbulent motions for the formation of Keplerian discs. We show that the accretion of mass and angular momentum in the surroundings of protostellar discs occurs in a highly anisotropic manner, by means of a few narrow accretion channels. The magnetic field structure in the vicinity of the discs is highly disordered, revealing field reversals up to distances of 1000 AU. These findings demonstrate that as soon as even mild turbulent motions are included, the classical disc formation scenario of a coherently rotating environment and a well-ordered magnetic field breaks down. Hence, it is highly questionable to assess the magnetic braking efficiency based on non-turbulent collapse simulation. We strongly suggest that, in addition to the global magnetic field properties, the small-scale accretion flow and detailed magnetic field structure have to be considered in order to assess the likelihood of Keplerian discs to be present.

Radiation feedback plays a crucial role in the process of star formation. In order to simulate the thermodynamic evolution of disks, filaments, and the molecular gas surrounding clusters of young stars, we require an efficient and accurate method for solving the radiation transfer problem. We describe the implementation of a hybrid radiation transport scheme in the adaptive grid-based FLASH general magnetohydrodynamics code. The hybrid scheme splits the radiative transport problem into a raytracing step and a diffusion step. The raytracer captures the first absorption event, as stars irradiate their environments, while the evolution of the diffuse component of the radiation field is handled by a flux-limited diffusion (FLD) solver. We demonstrate the accuracy of our method through a variety of benchmark tests including the irradiation of a static disk, subcritical and supercritical radiative shocks, and thermal energy equilibration. We also demonstrate the capability of our method for casting shadows and calculating gas and dust temperatures in the presence of multiple stellar sources. Our method enables radiation-hydrodynamic studies of young stellar objects, protostellar disks, and clustered star formation in magnetized, filamentary environments.

Massive exoplanets are observed preferentially around high metallicity ([Fe/H]) stars while low-mass exoplanets do not show such an effect. This so-called planet-metallicity correlation generally favors the idea that most observed gas giants at $r<10$ AU are formed via a core accretion process. We investigate the origin of this phenomenon using a semi-analystical model, wherein the standard core accretion takes place at planet traps in protostellar disks where rapid type I migrators are halted. We focus on the three major exoplanetary populations - hot-Jupiters, exo-Jupiters located at $r \simeq 1$ AU, and the low-mass planets. We show using a statistical approach that the planet-metallicity correlations are well reproduced in these models. We find that there are specific transition metallicities with values [Fe/H]$=-0.2$ to $-0.4$, below which the low-mass population dominates, and above which the Jovian populations take over. The exo-Jupiters significantly exceed the hot-Jupiter population at all observed metallicities. The low-mass planets formed via the core accretion are insensitive to metallicity, which may account for a large fraction of the observed super-Earths and hot-Neptunes. Finally, a controlling factor in building massive planets is the critical mass of planetary cores ($M_{c,crit}$) that regulates the onset of runaway gas accretion. Assuming the current data is roughly complete at [Fe/H]$>-0.6$, our models predict that the most likely value of the "mean" critical core mass of Jovian planets is $\braket{M_{c,crit}} \simeq 5 M_{\oplus}$ rather than $10 M_{\oplus}$. This implies that grain opacities in accreting envelopes should play an important role in lowering $M_{c,crit}$.

The class of meteorites called carbonaceous chondrites are examples of material from the solar system which have been relatively unchanged from the time of their initial formation. These meteorites have been classified according to the temperatures and physical conditions of their parent planetesimals. We collate available data on amino acid abundance in these meteorites and plot the concentrations of different amino acids for each meteorite within various meteorite subclasses. We plot average concentrations for various amino acids across meteorites separated by subclass and petrologic type. We see a predominance in the abundance and variety of amino acids in CM2 and CR2 meteorites. The range in temperature corresponding to these subclasses indicates high degrees of aqueous alteration, suggesting aqueous synthesis of amino acids. Within the CM2 and CR2 subclasses, we identify trends in relative frequencies of amino acids to investigate how common amino acids are as a function of their chemical complexity. These two trends (total abundance and relative frequencies) can be used to constrain formation parameters of amino acids within planetesimals. Our organization of the data supports an onion shell model for the temperature structure of planetesimals. The least altered meteorites (type 3) and their amino acids originated near cooler surface regions. The most active amino acid synthesis likely took place at intermediate depths (type 2). The most altered materials (type 1) originated furthest toward parent body cores. This region is likely too hot to either favor amino acid synthesis or for amino acids to be retained after synthesis.

Dust in protoplanetary disks is recognized as the building blocks of planets. In the core accretion scenario, the abundance of dust in disks (or metallicity) is crucial for forming cores of gas giants. We present our recent progress on the relationship between the metallicity and planet formation, wherein planet formation frequencies (PFFs) as well as the critical mass of planetary cores ($M_{c,crit}$) that can initiate gas accretion are statistically examined. We focus on three different planetary populations that are prominent for observed exoplanets in the mass-semimajor axis diagram: hot Jupiters, exo-Jupiters that are densely populated around 1 AU, and low-mass planets in tight orbits. We show that the resultant PFFs for both Jovian planets are correlated positively with the metallicity whereas low-mass planets form efficiently for a wide range of metallicities. This is consistent with the observed Planet-Metallicity correlation. Examining the statistically averaged value of $M_{c,crit}$ (defined as $<M_{c,crit}>$), we find that the correlation originates from the behavior of $<M_{c,crit}>$ that increases steadily with metallicity for two kinds of the Jovian planets while the low-mass planets obtain a rather constant value for $<M_{c,crit}>$. Such a difference in $<M_{c,crit}>$ can define transition metallicities (TMs) above which the Jovian planets gain a larger value of $<M_{c,crit}>$ than the low-mass planets, and hence gas giant formation takes place more efficiently. We find that TMs are sensitive to the important parameter that involves $M_{c,crit}$. We show, by comparing with the observations, that a most likely value of $M_{c,crit}$ is $\simeq 5M_{\oplus}$, which is smaller than the conventional value in the literature ($\simeq 10M_{\oplus}$). Our results suggest that opacities in planetary atmospheres play an important role for lowering $M_{c,crit}$.

(Abridged) In this review we focus on the observations and theory of the formation of early disks and outflows, and their connections with the first phases of planet formation. Large rotationally supported circumstellar disks, although common around more evolved young stellar objects, are rarely detected during the earliest, "Class 0" phase; however, a few excellent candidates have been discovered recently around both low and high mass protostars. In this early phase, prominent outflows are ubiquitously observed; they are expected to be associated with at least small magnetized disks. Disk formation - once thought to be a simple consequence of the conservation of angular momentum during hydrodynamic core collapse - is far more subtle in magnetized gas. In this case, the rotation can be strongly magnetically braked. Indeed, both analytic arguments and numerical simulations have shown that disk formation is suppressed in the strict ideal magnetohydrodynamic (MHD) limit for the observed level of core magnetization. We review what is known about this "magnetic braking catastrophe", possible ways to resolve it, and the current status of early disk observations. Outflows are also intimately linked to disk formation; they are a natural product of magnetic fields and rotation and are important signposts of star formation. We review new developments on early outflow generation since PPV. The properties of early disks and outflows are a key component of planet formation in its early stages and we review these major connections.

Recent studies of the nearest star-forming clouds of the Galaxy at submillimeter wavelengths with the Herschel Space Observatory have provided us with unprecedented images of the initial and boundary conditions of the star formation process. The Herschel results emphasize the role of interstellar filaments in the star formation process and connect remarkably well with nearly a decade's worth of numerical simulations and theory that have consistently shown that the ISM should be highly filamentary on all scales and star formation is intimately related to self-gravitating filaments. In this review, we trace how the apparent complexity of cloud structure and star formation is governed by relatively simple universal processes - from filamentary clumps to galactic scales. We emphasize two crucial and complementary aspects: (i) the key observational results obtained with Herschel over the past three years, along with relevant new results obtained from the ground on the kinematics of interstellar structures, and (ii) the key existing theoretical models and the many numerical simulations of interstellar cloud structure and star formation. We then synthesize a comprehensive physical picture that arises from the confrontation of these observations and simulations.

We present a model for the radiative output of star clusters in the process of star formation suitable for use in hydrodynamical simulations of radiative feedback. Gas in a clump, defined as a region whose density exceeds 10^4 cm^-3, is converted to stars via the random sampling of the Chabrier IMF. A star formation efficiency controls the rate of star formation. We have completed a suite of simulations which follow the evolution of accretion-fed clumps with initial masses ranging from 0 to 10^5 M_sol and accretion rates ranging from 10^-5 to 10^-1 M_sol yr^-1. The stellar content is tracked over time which allows the aggregate luminosity, ionizing photon rate, number of stars, and star formation rate (SFR) to be determined. For a fiducial clump of 10^4 M_sol, the luminosity is ~4x10^6 L_sol with a SFR of roughly 3x10^-3 M_sol yr^-1. We identify two regimes in our model. The accretion-dominated regime obtains the majority of its gas through accretion and is characterized by an increasing SFR while the reservoir-dominated regime has the majority of its mass present in the initial clump with a decreasing SFR. We show that our model can reproduce the expected number of O stars, which dominate the radiative output of the cluster. We find a nearly linear relationship between SFR and mass as seen in observations. We conclude that our model is an accurate and straightforward way to represent the output of clusters in hydrodynamical simulations with radiative feedback.

The rapid growth in the number of known exoplanets has revealed the existence of several distinct planetary populations in the observed mass-period diagram. Two of the most surprising are, (1) the concentration of gas giants around 1AU and (2) the accumulation of a large number of low-mass planets with tight orbits, also known as super-Earths and hot Neptunes. We have recently shown that protoplanetary disks have multiple planet traps that are characterized by orbital radii in the disks and halt rapid type I planetary migration. By coupling planet traps with the standard core accretion scenario, we showed that one can account for the positions of planets in the mass-period diagram. In this paper, we demonstrate quantitatively that most gas giants formed at planet traps tend to end up around 1 AU with most of these being contributed by dead zones and ice lines. In addition, we show that a large fraction of super-Earths and hot Neptunes are formed as "failed" cores of gas giants - this population being constituted by comparable contributions from dead zone and heat transition traps. Our results are based on the evolution of forming planets in an ensemble of disks where we vary only the lifetimes of disks as well as their mass accretion rates onto the host star. We show that a statistical treatment of the evolution of a large population of planetary cores initially caught in planet traps accounts for the existence of three distinct exoplantary populations - the hot Jupiters, the more massive planets at roughly orbital radii around 1 AU orbital, and the short period SuperEarths and hot Neptunes. There are very few evolutionary tracks that feed into the large orbital radii characteristic of the imaged Jovian planet, which agrees with recent surveys. Finally, we find that low-mass planets in tight orbits become the dominant planetary population for low mass stars ($M_* \le 0.7 M_{\odot}$).

MYStIX (Massive Young Star-Forming Complex Study in Infrared and Xray) seeks to characterize 20 OB-dominated young clusters and their environs at distances d < 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined criteria, and maps of nebular gas (particularly of hot X-ray emitting plasma) and dust. A catalog of MYStIX Probable Complex Members (MPCMs) with several hundred OB stars and > 30, 000 low mass premain sequence is assembled. This sample and related data products will be used to seek new empirical constraints on theoretical models of cluster formation and dynamics, mass segregation, OB star formation, star formation triggering on the periphery of HII regions, the survivability of protoplanetary disks in HII regions. This paper give an introduction and overview of the project, covering the data analysis methodology and application to two star forming regions, NGC 2264 and the Trifid Nebula.

Recent high-resolution observations show that protoplanetary disks have various kinds of structural properties or inhomogeneities. These are the consequence of a mixture of a number of physical and chemical processes taking place in the disks. Here, we discuss the results of our comprehensive investigations on how disk inhomogeneities affect planetary migration. We demonstrate that disk inhomogeneities give rise to planet traps - specific sites in protoplanetary disks at which rapid type I migration is halted. We show that up to three types of traps (heat transitions, ice lines and dead zones) can exist in a single disk, and that they move differently as the disk accretion rate decreases with time. We also demonstrate that the position of planet traps strongly depends on stellar masses and disk accretion rates. This indicates that host stars establish preferred (initial) scales of their planetary systems. Finally, we discuss the possible observational signatures of disk inhomogeneities.

We present the results of a study of simulated Giant Molecular Clouds (GMCs) formed in a Milky Way-type galactic disk with a flat rotation curve. This simulation, which does not include star formation or feedback, produces clouds with masses ranging between 10^4 Msun and 10^7 Msun. We compare our simulated cloud population to two observational surveys; The Boston University- Five College Radio Astronomy Observatory Galactic Ring Survey and the BIMA All-Disk Survey of M33. An analysis of the global cloud properties as well as a comparison of Larson's scaling relations is carried out. We find that simulated cloud properties agree well with the observed cloud properties, with the closest agreement occurring between the clouds at comparable resolution in M33. Our clouds are highly filamentary - a property that derives both from their formation due to gravitational instability in the sheared galactic environment, as well as to cloud- cloud gravitational encounters. We also find that the rate at which potentially star forming gas accumulates within dense regions - wherein n_{thresh} > 10^4 cm^{-3} - is 3% per 10 Myr, in clouds of roughly 10^6 Msun. This suggests that star formation rates in observed clouds are related to the rates at which gas can be accumulated into dense subregions within GMCs via filamentary flows. The most internally well-resolved clouds are chosen for listing in a catalogue of simulated GMCs; the first of its kind. The catalogued clouds are available as an extracted data set from the global simulation.

We present collapse simulations of strongly magnetised, turbulent molecular cloud cores with masses ranging from 2.6 to 1000 M_sun in order to study the influence of the initial conditions on the turbulence-induced disc formation mechanism proposed recently by Seifried et al. 2012. We find that Keplerian discs are formed in all cases independently of the core mass, the strength of turbulence, or the presence of global rotation. The discs appear within a few kyr after the formation of the protostar, are 50 - 150 AU in size, and have masses between 0.05 and a few 0.1 M_sun. During the formation of the discs the mass-to-flux ratio stays well below the critical value of 10 for Keplerian disc formation. Hence, flux-loss alone cannot explain the formation of Keplerian discs. The formation of rotationally supported discs at such early phases is rather due to the disordered magnetic field structure and due to turbulent motions in the surroundings of the discs, two effects lowering the classical magnetic braking efficiency. Binary systems occurring in the discs are mainly formed via the disc capturing mechanism rather than via disc fragmentation, which is largely suppressed by the presence of magnetic fields.

The large number of observed exoplanets ($\gtrsim $ 700) provides important constraints on their origin as deduced from the mass-period diagram of planets. The most surprising features in the diagram are 1) the (apparent) pile up of gas giants at a period of $\sim 500$ days ($\sim1$ AU) and 2) the so-called mass-period relation which indicates that planetary mass is an increasing function of orbital period. We construct the evolutionary tracks of growing planets at planet traps in evolving protoplanetary disks and show that they provide a good physical understanding of how these observational properties arise. The fundamental feature of our model is that inhomogeneities in protoplanetary disks give rise to multiple (up to 3) trapping sites for rapid (type I) planetary migration of planetary cores. The viscous evolution of disks results in the slow radial movement of the traps and their cores from large to small orbital periods. In our model, the slow inward motion of planet traps is coupled with the standard core accretion scenario for planetary growth. As planets grow, type II migration takes over. Planet growth and radial movement are ultimately stalled by the dispersal of gas disks via photoevaporation. Our model makes a number of important predictions: that distinct sub-populations of planets that reflect the properties of planet traps where they have grown result in the mass-period relation; that the presence of these sub-populations naturally explains a pile-up of planets at $\sim 1$ AU; and that evolutionary tracks from the ice line do put planets at short periods and fill an earlier claimed "planet desert" - sparse population of planets in the mass-semi-major axis diagram.

Recent observations and simulations have suggested that HII regions around massive stars may vary in their size and emitted flux on timescales short enough to be observed. This variability can have a number of causes, ranging from environmental causes to variability of the ionizing source itself. We explore the latter possibility by considering the pre-main-sequence evolution of massive protostars and conducting numerical simulations with ionizing radiation feedback using the FLASH AMR hydrodynamics code. We investigate three different models: a simple ZAMS model, a self-consistent one-zone model by Offner et al. (2009), and a model fit to the tracks computed by Hosokawa & Omukai (2009). The protostellar models show that hypercompact HII regions around massive, isolated protostars collapse or shrink from diameters of 80 or 300 AU, depending on the model choice, down to near absence during the swelling of stellar radius that accompanies the protostar's transition from a convective to a radiative internal structure. This occurs on timescales as short as ~3000 years.

We present collapse simulations of strongly magnetised, 100 M_sun, turbulent cloud cores. Around the protostars formed during the collapse Keplerian discs with typical sizes of up to 100 AU build up in contrast to previous simulations neglecting turbulence. Analysing the condensations in which the discs form, we show that the magnetic flux loss is not sufficient to explain the build-up of Keplerian discs. The average magnetic field is strongly inclined to the disc which might reduce the magnetic braking efficiency. However, the main reason for the reduced magnetic braking efficiency is the highly disordered magnetic field in the surroundings of the discs. Furthermore, due to the lack of a coherently rotating structure in the turbulent environment of the disc no toroidal magnetic field necessary for angular momentum extraction can build up. Simultaneously the angular momentum inflow remains high due to local shear flows created by the turbulent motions. We suggest that the "magnetic braking catastrophe" is an artefact of the idealised non-turbulent initial conditions and that turbulence provides a natural mechanism to circumvent this problem.

Long-lived, stable jets are observed in a wide variety of systems, from protostars, through Galactic compact objects to active galactic nuclei (AGN). Magnetic fields play a central role in launching, accelerating, and collimating the jets through various media. The termination of jets in molecular clouds or the interstellar medium deposits enormous amounts of mechanical energy and momentum, and their interactions with the external medium, as well, in many cases, as the radiation processes by which they are observed, are intimately connected with the magnetic fields they carry. This review focuses on the properties and structures of magnetic fields in long-lived jets, from their launch from rotating magnetized young stars, black holes, and their accretion discs, to termination and beyond. We compare the results of theory, numerical simulations, and observations of these diverse systems and address similarities and differences between relativistic and non-relativistic jets in protostellar versus AGN systems. On the observational side, we focus primarily on jets driven by AGN because of the strong observational constraints on their magnetic field properties, and we discuss the links between the physics of these jets on all scales.

We present collapse simulations of 100 M_{\sun}, turbulent cloud cores threaded by a strong magnetic field. During the initial collapse phase filaments are generated which fragment quickly and form several protostars. Around these protostars Keplerian discs with typical sizes of up to 100 AU build up in contrast to previous simulations neglecting turbulence. We examine three mechanisms potentially responsible for lowering the magnetic braking efficiency and therefore allowing for the formation of Keplerian discs. Analysing the condensations in which the discs form, we show that the build-up of Keplerian discs is neither caused by magnetic flux loss due to turbulent reconnection nor by the misalignment of the magnetic field and the angular momentum. It is rather a consequence of the turbulent surroundings of the disc which exhibit no coherent rotation structure while strong local shear flows carry large amounts of angular momentum. We suggest that the "magnetic braking catastrophe", i.e. the formation of sub-Keplerian discs only, is an artefact of the idealised non-turbulent initial conditions and that turbulence provides a natural mechanism to circumvent this problem.

(abridged) The nature of early Class 0/I protostellar discs is not clearly understood. Early protostellar discs are needed to drive molecular outflows and jets observed in star forming regions, but there has been some debate to how they form. From a theoretical perspective, the consequences of disc and outflow generation are crucial to understanding the very nature of how stars are assembled. We have performed 3D ideal magnetohydrodynamic (MHD) simulations of collapsing Bonnor-Ebert spheres, employing sink particles with a radius of 3.2 AU alongside an AMR grid and using a cooling function to model radiative cooling of the gas. This has allowed us to explore 2-8x10^4 yr further into the evolution of an early Class 0 disc-outflow system than previous simulations. Our outflow is precessing, kinked, turbulent, contains internal shocks and has a scale of 0.1 pc end-to-end. We form a rotationally dominated disc with a radius of 100 AU embedded inside a transient, unstable, flattened, rotating core extending out to 2000 AU. The larger flattened structure launches a low speed wind (v_r < 1.5 km/s) dominated by B_phi, while the inner disc launches a centrifugally driven jet dominated by B_p with speeds up to 20$ km/s. From the inner disk, the value of dM/dt_out/dM/dt_in ~ 0.1, wheras in the outer core dM/dt_out/dM/dt_in ~ 1.0. The inner disc becomes unstable to a warping instability due to the magnetic structure of the outflow and warps to 30 deg with respect to the z-axis by the end of the simulation. The envelope is cleared out and is less massive than the disc. We measure star formation efficiencies of eta_core=0.63 (and growing), higher than theoretical predictions. This indicates that outflows are not as efficient at expelling envelope mass as some current models estimate. We discuss the relevance of our disc misalignment concerning the formation of mis-aligned hot Jupiters.

Numerical simulations of outflows formed during the collapse of 100 M_sun cloud cores are presented. We derive a generalised criterion from MHD wind theory to analyse the launching mechanism of these outflows. The criterion is successfully applied to the whole outflow structure and cases with sub-Keplerian disc rotation. It allows us to decide whether an outflow is driven centrifugally or by the toroidal magnetic pressure. We show that quantities such as the magnetic field line inclination or the ratio of the toroidal to poloidal magnetic field alone are insufficient to determine the driving mechanism of outflows. By performing 12 runs we are able to study the influence of the initial conditions on the properties of outflows around massive protostars in detail. Our simulations reveal a strong effect of the magnetic field strength on the outflow morphology. In runs with weak fields or high rotational energies, well collimated, fast jets are observed whereas for strong fields poorly collimated, low-velocity outflows are found. We show that the occurrence of a fast jet is coupled to the existence of a Keplerian protostellar disc. Despite the very different morphologies all outflows are launched from the discs by centrifugal acceleration with the toroidal magnetic field increasingly contributing to the gas acceleration further away from the discs. The poor collimation of the outflows in runs with strong fields is a consequence of the weak hoop stresses. This in turn is caused by the slow build-up of a toroidal field due to sub-Keplerian disc rotation. The mass ejection/accretion ratios scatter around a mean of 0.3 in accordance with observational and analytical results. We suggest an evolutionary scenario for the earliest stage of massive star formation in which initially poorly collimated outflows develop which progressively get better collimated due to the generation of fast jets.

Supersonic turbulence is an essential element in understanding how structure within interstellar gas is created and shaped. In the context of star formation, many computational studies show that the mass spectrum of density and velocity fluctuations within dense clouds, as well as the distribution of their angular momenta, trace their origin to the statistical and physical properties of gas that is lashed with shock waves. In this article, we review the observations, simulations, and theories of how turbulent-like processes can account for structures we see in molecular clouds. We then compare traditional ideas of supersonic turbulence with a simpler physical model involving the effects of multiple shock waves and their interaction in the interstellar medium. Planar intersecting shock waves produce dense filaments, and generate vortex sheets that are essential to create the broad range of density and velocity structure in clouds. As an example, the lower mass behaviour of the stellar initial mass function can be traced to the tendency of a collection of shock waves to build-up a log-normal density distribution (or column density). Vorticity - which is essential to produce velocity structure over a very broad range of length scales in shocked clouds - can also be generated by the passage of curved shocks or intersecting planar shocks through such media. Two major additional physical forces affect the structure of star forming gas - gravity and feedback processes from young stars. Both of these can produce power-law tails at the high mass end of the IMF.

Radiative feedback is among the most important consequences of clustered star formation inside molecular clouds. At the onset of star formation, radiation from massive stars heats the surrounding gas, which suppresses the formation of many low-mass stars. When simulating pre-main-sequence stars, their stellar properties must be defined by a prestellar model. Different approaches to prestellar modeling may yield quantitatively different results. In this paper, we compare two existing prestellar models under identical initial conditions to gauge whether the choice of model has any significant effects on the final population of stars. The first model treats stellar radii and luminosities with a ZAMS model, while separately estimating the accretion luminosity by interpolating to published prestellar tracks. The second, more accurate prestellar model self-consistently evolves the radius and luminosity of each star under highly variable accretion conditions. Each is coupled to a raytracing-based radiative feedback code that also treats ionization. The impact of the self-consistent model is less ionizing radiation and less heating during the early stages of star formation. This may affect final mass distributions. We noted a peak stellar mass reduced by 8% from 47.3 Msun to 43.5 Msun in the evolutionary model, relative to the track-fit model. Also, the difference in mass between the two largest stars in each case is reduced from 14 Msun to 7.5 Msun. The HII regions produced by these massive stars were also seen to flicker on timescales down to the limit imposed by our timestep (< 560 years), rapidly changing in size and shape, confirming previous cluster simulations using ZAMS-based estimates for prestellar ionizing flux.

We present a model for the rotational evolution of a young, solar-mass star interacting magnetically with an accretion disk. As in a previous paper (Paper I), the model includes changes in the star's mass and radius as it descends the Hayashi track, a decreasing accretion rate, and a prescription for the angular momentum transfer between the star and disk. Paper I concluded that, for the relatively strong magnetic coupling expected in real systems, additional processes are necessary to explain the existence of slowly rotating pre-main-sequence stars. In the present paper, we extend the stellar spin model to include the effect of a spin-down torque that arises from an accretion-powered stellar wind. For a range of magnetic field strengths, accretion rates, initial spin rates, and mass outflow rates, the modeled stars exhibit rotation periods within the range of 1--10 days in the age range of 1--3 Myr. This range coincides with the bulk of the observed rotation periods, with the slow rotators corresponding to stars with the lowest accretion rates, strongest magnetic fields, and/or highest stellar wind mass outflow rates. We also make a direct, quantitative comparison between the accretion-powered stellar wind scenario and the two types of disk-locking models (namely the X-wind and Ghosh & Lamb type models) and identify some remaining theoretical issues for understanding young star spins.

Planetary migration provides a theoretical basis for the observed diversity of exoplanetary systems. We demonstrate that dust settling - an inescapable feature of disk evolution - gives even more rapid type I migration by up to a factor of about 2 than occurs in disks with fully mixed dust. On the other hand, type II migration becomes slower by a factor of 2 due to dust settling. This even more problematic type I migration can be resolved by the presence of a dead zone; the inner, high density region of a disk which features a low level of turbulence. We show that enhanced dust settling in the dead zone leaves a dusty wall at its outer edge. Back-heating of the dead zone by this wall produces a positive radial gradient for the disk temperature, which acts as a barrier for type I migration.

We present simulations of collapsing 100 M_\sun mass cores in the context of massive star formation. The effect of variable initial rotational and magnetic energies on the formation of massive stars is studied in detail. We focus on accretion rates and on the question under which conditions massive Keplerian disks can form in the very early evolutionary stage of massive protostars. For this purpose, we perform 12 simulations with different initial conditions extending over a wide range in parameter space. The equations of magnetohydrodynamics (MHD) are solved under the assumption of ideal MHD. We find that the formation of Keplerian disks in the very early stages is suppressed for a mass-to-flux ratio normalised to the critical value \mu below 10, in agreement with a series of low-mass star formation simulations. This is caused by very efficient magnetic braking resulting in a nearly instantaneous removal of angular momentum from the disk. For weak magnetic fields, corresponding to \mu > 10, large-scale, centrifugally supported disks build up with radii exceeding 100 AU. A stability analysis reveals that the disks are supported against gravitationally induced perturbations by the magnetic field and tend to form single stars rather than multiple objects. We find protostellar accretion rates of the order of a few 10^-4 M_\sun yr^-1 which, considering the large range covered by the initial conditions, vary only by a factor of ~ 3 between the different simulations. We attribute this fact to two competing effects of magnetic fields. On the one hand, magnetic braking enhances accretion by removing angular momentum from the disk thus lowering the centrifugal support against gravity. On the other hand, the combined effect of magnetic pressure and magnetic tension counteracts gravity by exerting an outward directed force on the gas in the disk thus reducing the accretion onto the protostars.

The structure of planetary systems around their host stars depends on their initial formation conditions. Massive planets will likely be formed as a consequence of rapid migration of planetesimals and low mass cores into specific trapping sites in protoplanetary discs. We present analytical modeling of inhomogeneities in protoplanetary discs around a variety of young stars, - from Herbig Ae/Be to classical T Tauri and down to M stars, - and show how they give rise to planet traps. The positions of these traps define the initial orbital distribution of multiple protoplanets. We investigate both corotation and Lindblad torques, and show that a new trap arises from the (entropy-related) corotation torque. This arises at that disc radius where disc heating changes from viscous to stellar irradiation dominated processes. We demonstrate that up to three traps (heat transitions, ice lines and dead zones) can exist in a single disc, and that they move differently as the disc accretion rate $\dot{M}$ decreases with time. The interaction between the giant planets which grow in such traps may be a crucial ingredient for establishing planetary systems. We also demonstrate that the position of planet traps strongly depends on stellar masses and disc accretion rates. This indicates that host stars establish preferred scales of planetary systems formed around them. We discuss the potential of planet traps induced by ice lines of various molecules such as water and CO, and estimate the maximum and minimum mass of planets which undergo type I migration. We finally apply our analyses to accounting for the initial conditions proposed in the Nice model for the origin of our Solar system.

We review computational approaches to understanding the origin of the Initial Mass Function (IMF) during the formation of star clusters. We examine the role of turbulence, gravity and accretion, equations of state, and magnetic fields in producing the distribution of core masses - the Core Mass Function (CMF). Observations show that the CMF is similar in form to the IMF. We focus on feedback processes such as stellar dynamics, radiation, and outflows can reduce the accreted mass to give rise to the IMF. Numerical work suggests that filamentary accretion may play a key role in the origin of the IMF.

Planetary migration is essential to explain the observed mass-period relation for exoplanets. Without some stopping mechanism, the tidal, resonant interaction between planets and their gaseous disc generally causes the planets to migrate inward so efficiently that they plunge into the host star within the gaseous disc lifetime ($\sim $ 1-3 Myrs). We investigate planetary migration by analytically calculating the migration rate and time within self-consistently computed, radiatively heated discs around M stars in which the effects of dust settling are included. We show that dust settling lowers the disc temperature and raises the gas density in the mid-plane. This inescapable evolution of disc structure speeds up type I planetary migration for lower mass bodies by up to a factor of about 2. We also examine the effects of dust settling on the gap-opening mass and type II migration, and find that the gap-opening mass is reduced by a factor of 2 and type II migration becomes slower by a factor of 2. While dust settling can somewhat alleviate the problem of planetary migration for more massive planets, the more rapid migration of low mass planets and planetary cores requires a robust slowing mechanism.

The pre-stellar cores in which low mass stars form are generally well magnetized. Our simulations show that early protostellar discs are massive and experience strong magnetic torques in the form of magnetic braking and protostellar outflows. Simulations of protostellar disk formation suggest that these torques are strong enough to suppress a rotationally supported structure from forming for near critical values of mass-to-flux. We demonstrate through the use of a 3D adaptive mesh refinement code -- including cooling, sink particles and magnetic fields -- that one produces transient 1000 AU discs while simultaneously generating large outflows which leave the core region, carrying away mass and angular momentum. Early inflow/outflow rates suggest that only a small fraction of the mass is lost in the initial magnetic tower/jet event.

We perform state-of-the-art, 3D, time-dependent simulations of magnetized disk winds, carried out to simulation scales of 60 Astronomical Units, in order to confront optical HST observations of protostellar jets. We ``observe'' the optical forbidden line emission produced by shocks within our simulated jets and compare these with actual observations. Our simulations reproduce the rich structure of time varying jets, including jet rotation far from the source, an inner (up to 400 km/s) and outer (less than 100 km/s) component of the jet, and jet widths of up to 20 Astronomical Units in agreement with observed jets. These simulations when compared with the data are able to constrain disk wind models. In particular, models featuring a disk magnetic field with a modest radial spatial variation across the disk are favored.

Planetary migration in standard models of gaseous protoplanetary disks is known to be very rapid ($\sim 10^5$ years) jeopardizing the existence of planetary systems. We present a new mechanism for significantly slowing rapid planetary migration, discovered by means of radiative transfer calculations of the thermal structure of protoplanetary disks irradiated by their central stars. Rapid dust settling in a disk's dead zone - a region with very little turbulence - leaves a dusty wall at its outer edge. We show that the back-heating of the dead zone by this irradiated wall produces a positive gradient of the disk temperature which acts as a thermal barrier to planetary migration which persists for the disk lifetime. Although we analyze in detail the migration of a Super-Earth in a low mass disk around an M star, our findings can apply to wide variety of young planetary systems. We compare our findings with other potentially important stopping mechanisms and show that there are large parameter spaces for which dead zones are likely to play the most important role for reproducing the observed mass-period relation in longer planetary periods.

In star formation, magnetic fields act as a cosmic angular momentum extractor that increases mass accretion rates onto protostars and in the process, creates spectacular outflows. However, recently it has been argued that this magnetic brake is so strong that early protostellar disks -- the cradles of planet formation -- cannot form. Our three-dimensional numerical simulations of the early stages of collapse (\lesssim 10^5 yr) of overdense star--forming clouds form early outflows and have magnetically regulated and rotationally dominated disks (inside 10 AU) with high accretion rates, despite the slip of the field through the mostly neutral gas. We find that in three dimensions, magnetic fields suppress gravitationally driven instabilities which would otherwise prevent young, well ordered disks from forming. Our simulations have surprising consequences for the early formation of disks, their density and temperature structure, the mechanism and structure of early outflows, the flash heating of dust grains through ambipolar diffusion, and the origin of planets and binary stars.

The irradiation of protoplanetary discs by central stars is the main heating mechanism for discs, resulting in their flared geometric structure. In a series of papers, we investigate the deep links between 2D self-consistent disc structure and planetary migration in irradiated discs, focusing particularly on those around M stars. In this first paper, we analyse the thermal structure of discs that are irradiated by an M star by solving the radiative transfer equation by means of a Monte Carlo code. Our simulations of irradiated hydrostatic discs are realistic and self-consistent in that they include dust settling with multiple grain sizes (N=15), the gravitational force of an embedded planet on the disc, and the presence of a dead zone (a region with very low levels of turbulence) within it. We show that dust settling drives the temperature of the mid-plane from an $r^{-3/5}$ distribution (well mixed dust models) toward an $r^{-3/4}$. The dead zone, meanwhile, leaves a dusty wall at its outer edge because dust settling in this region is enhanced compared to the active turbulent disc at larger disc radii. The disc heating produced by this irradiated wall provides a positive gradient region of the temperature in the dead zone in front of the wall. This is crucially important for slowing planetary migration because Lindblad torques are inversely proportional to the disc temperature. Furthermore, we show that low turbulence of the dead zone is self-consistently induced by dust settling, resulting in the Kelvin-Helmholtz instability (KHI). We show that the strength of turbulence arising from the KHI in the dead zone is $\alpha=10^{-5}$.

The main observational evidence for turbulence in the interstellar medium (ISM) and molecular clouds is the power-law energy spectrum for velocity fluctuations, E(k) \propto k^{\alpha}. The Kolmogorov scaling exponent, \alpha=-5/3, is typical. At the same time, the observed probability distribution function (PDF) of gas densities in both the ISM as well as in molecular clouds is a log-normal distribution, which is similar to the initial mass function (IMF) that describes the distribution of stellar masses. We examine the density and velocity structure of interstellar gas traversed by curved shock waves in the kinematic limit. We demonstrate mathematically that just a few passages of curved shock waves generically produces a log-normal density PDF. This explains the ubiquity of the log-normal PDF in many different numerical simulations. We also show that subsequent interaction with a spherical blast wave generates a power-law density distribution at high densities, qualitatively similar to the Salpeter power-law for the IMF. Finally, we show that a focused shock produces a downstream flow with energy spectrum exponent \alpha=-2. Subsequent shock passages reduce this slope, achieving \alpha \approx -5/3 after a few passages. Subsequent dissipation of energy piled up at the small scales will act to maintain the spectrum very near to the Kolomogorov value. Therefore, fully-developed turbulence may not be required to explain the observed energy spectrum and density PDF. We argue that the self-similar spherical blast wave arising from expanding HII regions or stellar winds from massive stars may ultimately be responsible for creating the high mass, power-law, tail in the IMF.

We examine the case for Quark-Novae (QNe) as possible sources for the reionization and early metal enrichment of the universe. Quark-Novae are predicted to arise from the explosive collapse (and conversion) of sufficiently massive neutron stars into quark stars. A Quark-Nova (QN) can occur over a range of time scales following the supernova event. For QNe that arise days to weeks after the supernovae, we show that dual-shock that arises as the QN ejecta encounter the supernova ejecta can produce enough photons to reionize hydrogen in most of the Inter-Galactic medium (IGM) by z ~ 6. Such events can explain the large optical depth tau_e ~ 0.1 as measured by WMAP, if the clumping factor, C, of the material being ionized is smaller than 10. We suggest a way in which a normal initial mass function (IMF) for the oldest stars can be reconciled with a large optical depth as well as the mean metallicity of the early IGM post reionization. We find that QN also make a contribution to r-process element abundances for atomic numbers A > 130. We predict that the main cosmological signatures of Quark-Novae are the gamma-ray bursts that announce their birth. These will be clustered at redshifts in the range z ~ 7-8 in our model.

Of the twenty amino acids used in proteins, ten were formed in Miller's atmospheric discharge experiments. The two other major proposed sources of prebiotic amino acid synthesis include formation in hydrothermal vents and delivery to Earth via meteorites. We combine observational and experimental data of amino acid frequencies formed by these diverse mechanisms and show that, regardless of the source, these ten early amino acids can be ranked in order of decreasing abundance in prebiotic contexts. This order can be predicted by thermodynamics. The relative abundances of the early amino acids were most likely reflected in the composition of the first proteins at the time the genetic code originated. The remaining amino acids were incorporated into proteins after pathways for their biochemical synthesis evolved. This is consistent with theories of the evolution of the genetic code by stepwise addition of new amino acids. These are hints that key aspects of early biochemistry may be universal.

Given the importance of simulating hydromagnetic processes that impact star formation, we have earlier developed a 3D adaptive mesh approach that allows us to include hydromagnetic processes during the formation and evolution of cores, discs, and stars in observed regions of star formation. In this paper, we take the next step in this program - namely - to develop a modified version of the 3D adaptive mesh refinement (AMR) code FLASH in which the ambipolar diffusion of the magnetic field in poorly ionized molecular gas is implemented. We approach the problem using a single-fluid approximation to simplify numerical calculations. In this paper, we present a series of test cases including oblique isothermal and non--isothermal C-shocks. We also present a study of the quasi-static collapse of an initial uniform, self-gravitating, magnetized sphere that is initially supported by its magnetic field against collapse (i.e. magnetically subcritical). Applications to the collapse of a pre-stellar Bonnor-Ebert sphere are presented in a companion paper.

Using numerical simulations of magnetized stellar winds, we carry out a parameter study to find the dependence of the stellar wind torque on observable parameters. We find that the power-law dependencies of the torque on parameters is significantly different than what has been used in all spin evolution models to date.

The growth of Jovian mass planets during migration in their protoplanetary disks is one of the most important problems that needs to be solved in light of observations of the exosolar planets. Studies of the migration of planets in standard gas disk models routinely show that migration is too fast to form Jovian planets, and that such migrating planetary cores generally plunge into the central stars in less than a Myr. In previous work, we have shown that a poorly ionized, less viscous region in a protoplanetary disk called a dead zone slows down the migration of fixed-mass planets. In this paper, we extend our numerical calculations to include dead zone evolution along with the disk, as well as planet formation via accretion of rocky and gaseous materials. Using our symplectic-integrator-gas dynamics code, we find that dead zones, even in evolving disks wherein migrating planets grow by accretion, still play a fundamental role in saving planetary systems. We demonstrate that Jovian planets form within 2.5 Myr for disks that are ten times more massive than a minimum mass solar nebula (MMSN) with an opacity reduction and without slowing down migration artificially. Our simulations indicate that protoplanetary disks with an initial mass comparable to the MMSN only produce Neptunian mass planets. We also find that planet migration does not help core accretion as much in the oligarchic planetesimal accretion scenario as it was expected in the runaway accretion scenario. Therefore we expect that an opacity reduction (or some other mechanisms) is needed to solve the formation timescale problem even for migrating protoplanets, as long as we consider the oligarchic growth. We also point out a possible role of a dead zone in explaining long-lived, strongly accreting gas disks.

Astrophysical jets are associated with the formation of young stars of all masses, stellar and massive black holes, and perhaps even with the formation of massive planets. Their role in the formation of planets, stars, and galaxies is increasingly appreciated and probably reflects a deep connection between the accretion flows - by which stars and black holes may be formed - and the efficiency by which magnetic torques can remove angular momentum from such flows. We compare the properties and physics of jets in both non-relativistic and relativistic systems and trace, by means of theoretical argument and numerical simulations, the physical connections between these different phenomena. We discuss the properties of jets from young stars and black holes, give some basic theoretical results that underpin the origin of jets in these systems, and then show results of recent simulations on jet production in collapsing star-forming cores as well as from jets around rotating Kerr black holes.

We present results from our numerical simulations of collapsing massive molecular cloud cores. These numerical calculations show that massive stars assemble quickly with mass accretion rates exceeding 10^-3 Msol/yr and confirm that the mass accretion during the collapsing phase is much more efficient than predicted by selfsimilar collapse solutions, dM/dt ~ c^3/G. We find that during protostellar assembly out of a non-turbulent core, the mass accretion reaches 20 - 100 c^3/G. Furthermore, we explore the self-consistent structure of bipolar outflows that are produced in our three dimensional magnetized collapse simulations. These outflows produce cavities out of which radiation pressure can be released, thereby reducing the limitations on the final mass of massive stars formed by gravitational collapse. Additional enhancement of the mass accretion rate comes from accretion along filaments that are built up by supersonic turbulent motions. Our numerical calculations of collapsing turbulent cores result in mass accretion rates as high as 10^-2 Msol/yr.

We compare the stellar wind torque calculated in a previous work (Paper II) to the spin-up and spin-down torques expected to arise from the magnetic interaction between a slowly rotating ($\sim 10$% of breakup) pre-main-sequence star and its accretion disk. This analysis demonstrates that stellar winds can carry off orders of magnitude more angular momentum than can be transferred to the disk, provided that the mass outflow rates are greater than the solar wind. Thus, the equilibrium spin state is simply characterized by a balance between the angular momentum deposited by accretion and that extracted by a stellar wind. We derive a semi-analytic formula for predicting the equilibrium spin rate as a function only of the ratio of $\dot M_{\rm w} / \dot M_{\rm a}$ and a dimensionless magnetization parameter, $\Psi \equiv B_*^2 R_*^2 (\dot M_{\rm a} v_{\rm esc})^{-1}$, where $\dot M_{\rm w}$ is the stellar wind mass outflow rate, $\dot M_{\rm a}$ the accretion rate, $B_*$ the stellar surface magnetic field strength, $R_*$ the stellar radius, and $v_{\rm esc}$ the surface escape speed. For parameters typical of accreting pre-main-sequence stars, this explains spin rates of $\sim 10$% of breakup speed for $\dot M_{\rm w} / \dot M_{\rm a} \sim 0.1$. Finally, the assumption that the stellar wind is driven by a fraction of the accretion power leads to an upper limit to the mass flow ratio of $\dot M_{\rm w} / \dot M_{\rm a} \la 0.6$.

[Abridged] In order to explain the slow rotation observed in a large fraction of accreting pre-main-sequence stars (CTTSs), we explore the role of stellar winds in torquing down the stars. For this mechanism to be effective, the stellar winds need to have relatively high outflow rates, and thus would likely be powered by the accretion process itself. Here, we use numerical magnetohydrodynamical simulations to compute detailed 2-dimensional (axisymmetric) stellar wind solutions, in order to determine the spin down torque on the star. We explore a range of parameters relevant for CTTSs, including variations in the stellar mass, radius, spin rate, surface magnetic field strength, the mass loss rate, and wind acceleration rate. We also consider both dipole and quadrupole magnetic field geometries. Our simulations indicate that the stellar wind torque is of sufficient magnitude to be important for spinning down a ``typical'' CTTS, for a mass loss rate of $\sim 10^{-9} M_\odot$ yr$^{-1}$. The winds are wide-angle, self-collimated flows, as expected of magnetic rotator winds with moderately fast rotation. The cases with quadrupolar field produce a much weaker torque than for a dipole with the same surface field strength, demonstrating that magnetic geometry plays a fundamental role in determining the torque. Cases with varying wind acceleration rate show much smaller variations in the torque suggesting that the details of the wind driving are less important. We use our computed results to fit a semi-analytic formula for the effective Alfv\'en radius in the wind, as well as the torque. This allows for considerable predictive power, and is an improvement over existing approximations.

(Abridged) We present a series of decaying turbulence simulations that represent a cluster-forming clump within a molecular cloud, investigating the role of magnetic fields on the formation of potential star-forming cores. We present an exhaustive analysis of numerical data from these simulations that includes a compilation of all of the distributions of physical properties that characterize bound cores - including their masses, radii, mean densities, angular momenta, spins, magnetizations, and mass-to-flux ratios. We also present line maps of our models that can be compared with observations. Our simulations range between 5-30 Jeans masses of gas, and are representative of molecular cloud clumps with masses between 100-1000 solar masses. The cores have mass-to-flux ratios that are generally less than that of the original cloud, and so a cloud that is initially highly supercritical can produce cores that are slightly supercritical, similar to that seen by Zeeman measurements of molecular cloud cores. Clouds that are initially only slightly supercritical will instead collapse along the field lines into sheets, and the cores that form as these sheets fragment have a different mass spectrum than what is observed. The spin rates of these cores suggests that subsequent fragmentation into multiple systems is likely. The sizes of the bound cores that are produced are typically 0.02-0.2 pc and have densities in the range 10^4-10^5 cm^{-3} in agreement with observational surveys. Finally, our numerical data allow us to test theoretical models of the mass spectrum of cores, such as the turbulent fragmentation picture of Padoan-Nordlund. We find that while this model gets the shape of the core mass spectrum reasonably well, it fails to predict the peak mass in the core mass spectrum.

Stellar winds may be important for angular momentum transport from accreting T Tauri stars, but the nature of these winds is still not well-constrained. We present some simulation results for hypothetical, hot (~1e6 K) coronal winds from T Tauri stars, and we calculate the expected emission properties. For the high mass loss rates required to solve the angular momentum problem, we find that the radiative losses will be much greater than can be powered by the accretion process. We place an upper limit to the mass loss rate from accretion-powered coronal winds of ~1e-11 solar masses per year. We conclude that accretion powered stellar winds are still a promising scenario for solving the stellar angular momentum problem, but the winds must be cool (~1e4 K) and thus are not driven by thermal pressure.

We review the theoretical efforts to understand why pre-main-sequence stars spin much more slowly than expected. The first idea put forward was that massive stellar winds may remove substantial angular momentum. Since then, it has become clear that the magnetic interaction between the stars and their accretion disks explains many of the observed emission properties. The disk locking scenario, which assumes the magnetic star-disk interaction also solves the stellar spin problem, has received the most attention in the literature to date. However, recent considerations suggest that the torques in the star-disk interaction are insufficient for disk locking to explain the slow rotators. This prompts us to revisit stellar winds, and we conclude that stellar winds, working in conjunction with magnetospheric accretion, are a promising candidate for solving the angular momentum problem. We suggest future directions for both observations and theory, to help shed light on this issue.

The tidal interaction between a disk and a planet leads to the planet's migration. A long-standing question regarding this mechanism is how to stop the migration before planets plunge into their central stars. In this paper, we propose a new, simple mechanism to significantly slow down planet migration, and test the possibility by using a hybrid numerical integrator to simulate the disk-planet interaction. The key component of the scenario is the role of low viscosity regions in protostellar disks known as dead zones, which affect planetary migration in two ways. First of all, it allows a smaller-mass planet to open a gap, and hence switch the faster type I migration to the slower type II migration. Secondly, a low viscosity slows down type II migration itself, because type II migration is directly proportional to the viscosity. We present numerical simulations of planetary migration by using a hybrid symplectic integrator-gas dynamics code. Assuming that the disk viscosity parameter inside the dead zone is (alpha=1e-4-1e-5), we find that, when a low-mass planet (e.g. 1-10 Earth masses) migrates from outside the dead zone, its migration is stopped due to the mass accumulation inside the dead zone. When a low-mass planet migrates from inside the dead zone, it opens a gap and slows down its migration. A massive planet like Jupiter, on the other hand, opens a gap and slows down inside the dead zone, independent of its initial orbital radius. The final orbital radius of a Jupiter mass planet depends on the dead zone's viscosity. For the range of alpha's noted above, this can vary anywhere from 7 AU, to an orbital radius of 0.1 AU that is characteristic of the hot Jupiters.

We present an investigation of massive star formation that results from the gravitational collapse of massive, magnetized molecular cloud cores. We investigate this by means of highly resolved, numerical simulations of initial magnetized Bonnor-Ebert-Spheres that undergo collapse and cooling. By comparing three different cases - an isothermal collapse, a collapse with radiative cooling, and a magnetized collapse - we show that massive stars assemble quickly with mass accretion rates exceeding 10^-3 Msol/yr. We confirm that the mass accretion during the collapsing phase is much more efficient than predicted by selfsimilar collapse solutions, i.e. dM/dt ~ c^3/G. We find that during protostellar assembly the mass accretion reaches 20 - 100 c^3/G. Furthermore, we determined the self-consistent structure of bipolar outflows that are produced in our three dimensional magnetized collapse simulations. These outflows produce cavities out of which radiation pressure can be released, thereby reducing the limitations on the final mass of massive stars formed by gravitational collapse. Moreover, we argue that the extraction of angular momentum by disk-threaded magnetic fields and/or by the appearance of bars with spiral arms significantly enhance the mass accretion rate, thereby helping the massive protostar to assemble more quickly.

We present a detailed computational study of the assembly of protostellar disks and massive stars in molecular clouds with supersonic turbulence. We follow the evolution of large scale filamentary structures in a cluster-forming clump down to protostellar length scales by means of very highly resolved, 3D adaptive mesh refined (AMR) simulations, and show how accretion disks and massive stars form in such environments. We find that an initially elongated cloud core which has a slight spin from oblique shocks collapses first to a filament and later develops a turbulent disk close to the center of the filament. The continued large scale flow that shocks with the filament maintains the high density and pressure within it. Material within the cooling filament undergoes gravitational collapse and an outside-in assembly of a massive protostar. Our simulations show that very high mass accretion rates of up to 10^-2 Msol/yr and high, supersonic, infall velocities result from such filamentary accretion. Accretion at these rates is higher by an order of magnitude than those found in semi-analytic studies, and can quench the radiation field of a growing massive young star.Our simulations include a comprehensive set of the important chemical and radiative processes such as cooling by molecular line emission, gas-dust interaction, and radiative diffusion in the optical thick regime, as well as H2 formation and dissociation. Therefore, we are able to probe, for the first time, the relevant physical phenomena on all scales from those characterizing the clump down to protostellar core.

We review advances in the theoretical and computational studies of disk winds, jets and outflows including: the connection between accretion and jets, the launch of jets from magnetized disks, the coupled evolution of jets and disks, the interaction of magnetized young stellar objects with their surrounding disks and the relevance to outflows, and finally, the link between jet formation and gravitational collapse. We also address the predictions that the theory makes about jet kinematics, collimation, and rotation, that have recently been confirmed by high spatial and spectral resolution observations. Disk winds have a universal character that may account for jets and outflows during the formation of massive stars as well as brown dwarfs.

Star formation is usually accompanied by outflow phenomena. There is strong evidence that these outflows and jets are launched from the protostellar disk by magneto-rotational processes. Here, we report on our three dimensional, adaptive mesh, magneto-hydrodynamic simulations of collapsing, rotating, magnetized Bonnor-Ebert-Spheres whose properties are taken directly from observations. In contrast to the pure hydro case where no outflows are seen, our present simulations show an outflow from the protodisk surface at ~ AU and a jet at ~ 0.07 AU after a strong toroidal magnetic field build up. The large scale outflow, which extends up to ~ AU at the end of our simulation, is driven by toroidal magnetic pressure (spring), whereas the jet is powered by magneto-centrifugal force (fling). At the final stage of our simulation these winds are still confined within two respective shock fronts. Furthermore, we find that the jet-wind and the disk-anchored magnetic field extracts a considerable amount of angular momentum from the protostellar disk. The initial spin of our cloud core was chosen high enough to produce a binary system. We indeed find a close binary system (separation ~ 3 R_sol) which results from the fragmentation of an earlier formed ring structure. The magnetic field strength in these protostars reaches ~ 3 kG and becomes about 3 G at 1 AU from the center in agreement with recent observational results.

Most low-mass protostellar disks evolve in clustered environments where they are affected by external radiation fields, while others evolve in more isolated star-forming regions. Assuming that the magneto-rotational instability (MRI) is the main source of viscosity, we calculate the size of a poorly ionized, MRI inactive, and hence low viscosity region - the "dead zone" - in these protostellar disks. We include disk ionization by X-rays, cosmic rays, radioactive elements and thermal collisions, recombination by molecules, metals, and grains, as well as the effect of turbulence stimulation in the dead zone by the active layers lying above it. We also calculate the gap-opening masses of planets, which are determined by a disk's viscosity and a disk aspect ratio, for disks in these environments and compare them with each other. We find that the dead zone is a robust feature of the protostellar disks that is largely independent of their environment, typically stretching out to ~ 15 AU. We analyze the possible effects of dead zones on planet formation, migration, and eccentricity evolution. We show that the gap-opening mass inside the dead zone is expected to be of the order of terrestrial and ice giant mass planets while that outside the dead zone is Jovian or super-Jovian mass planets, largely independent of the star-forming environment. We show that dead zones can significantly slow down both type I and type II planetary migration due to their lower viscosity. We also find that the growth of eccentricity of massive extrasolar planets is particularly favorable through the planet-disk interaction inside the dead zones due to the large gaps expected to be opened by planets.

We compare the angular momentum extracted by a wind from a pre-main-sequence star to the torques arising from the interaction between the star and its Keplerian accretion disk. We find that the wind alone can counteract the spin-up torque from mass accretion, solving the mystery of why accreting pre-main-sequence stars are observed to spin at less than 10% of break-up speed, provided that the mass outflow rate in the stellar winds is ~10% of the accretion rate. We suggest that such massive winds will be driven by some fraction $\epsilon$ of the accretion power. For observationally constrained typical parameters of classical T-Tauri stars, $\epsilon$ needs to be between a few and a few tens of percent. In this scenario, efficient braking of the star will terminate simultaneously with accretion, as is usually assumed to explain the rotation velocities of stars in young clusters.

(Abriged) We present a comprehensive set of axisymmetric, time-dependent simulations of jets from Keplerian disks whose mass loading as a function of disk radius is systematically changed. For a reasonable model for the density structure and injection speed of the underlying accretion disk, mass loading is determined by the radial structure of the disk's magnetic field structure. We vary this structure by using four different magnetic field configurations, ranging from the "potential" configuration (Ouyed&Pudritz 1997), to the increasingly more steeply falling Blandford&Payne (1982) and Pelletier&Pudritz (1992) models, and ending with a quite steeply raked configuration that bears similarities to the Shu X-wind model. We find that the radial distribution of the mass load has a profound effect on both the rotational profile of the underlying jet as well as the degree of collimation of its outflow velocity and magnetic field lines. We show analytically, and confirm by our simulations, that the collimation of a jet depends on its radial current distribution, which in turn is prescribed by the mass load. Models with steeply descending mass loads have strong toroidal fields, and these collimate to cylinders (this includes the Ouyed-Pudritz and Blandford-Payne outflows). On the other hand, the more gradually descending mass load profiles (the PP92 and monopolar distributions) have weaker toroidal fields, and these result in wide-angle outflows with parabolic collimation. We also present detailed structural information about jets such as their radial profiles of jet density, toroidal magnetic field, and poloidal jet speed, as well as an analysis of the bulk energetics of our different simulations.

A large body of theoretical and computational work shows that jets - modelled as magnetized disk winds - exert an external torque on their underlying disks that can efficiently remove angular momentum and act as major drivers of disk accretion. These predictions have recently been confirmed in direct HST measurements of the jet rotation and angular momentum transport in low mass protostellar systems. We review the theory of disc winds and show that their physics is universal and scales to jets from both low and high mass star forming regions. This explains the observed properties of outflows in massive star forming regions, before the central massive star generates an ultracompact HII region. We also discuss the recent numerical studies on the formation of massive accretion disks and outflows through gravitational collapse, including our own work on 3D Adaptive Mesh simulations (using the FLASH code) of the hydromagnetic collapse of an initial rotating, and cooling Bonner-Ebert sphere. Magnetized collapse gives rise to outflows. Our own simulations show that both a jet-like disk wind on sub AU scales, and a larger scale molecular outflow occur (Banerjee and Pudritz 2005).

Protoplanets accrete material from their natal protostellar disks until they are sufficiently massive to open a gap in the face of the disk's viscosity that arises from the magneto-rotational instability (MRI). By computing the ionization structure within observationally well-constrained disk models, we demonstrate that poorly ionized, low viscosity "dead zones" stretch out to 12 AU within typical disks. We find that planets of terrestrial mass robustly form within the dead zones while massive Jovian planets form beyond. Dead zones will also halt the rapid migration of planets into their central stars. Finally, we argue that the gravitational scattering of low mass planets formed in the dead zone, to larger radii by a rapidly accreting Jupiter beyond, can explain the distribution of planetary masses in our solar system.

We formulate a general, steady-state model for the torque on a magnetized star from a surrounding accretion disc. For the first time, we include the opening of dipolar magnetic field lines due to the differential rotation between the star and disc, so the magnetic topology then depends on the strength of the magnetic coupling to the disc. This coupling is determined by the effective slip rate of magnetic field lines that penetrate the diffusive disc. Stronger coupling (i.e., lower slip rate) leads to a more open topology and thus to a weaker magnetic torque on the star from the disc. In the expected strong coupling regime, we find that the spin-down torque on the star is more than an order of magnitude smaller than calculated by previous models. We also use our general approach to examine the equilibrium (`disc-locked') state, in which the net torque on the star is zero. In this state, we show that the stellar spin rate is roughly an order of magnitude faster than predicted by previous models. This challenges the idea that slowly-rotating, accreting protostars are disc locked. Furthermore, when the field is sufficiently open (e.g., for mass accretion rates > 5 x 10^{-9} M_sun / yr, for typical accreting protostars), the star will receive no magnetic spin-down torque from the disc at all. We therefore conclude that protostars must experience a spin-down torque from a source that has not yet been considered in the star-disc torque models--possibly from a stellar wind along the open field lines.

We present a detailed study of the collapse of molecular cloud cores using high resolution 3D adaptive mesh refinement (AMR) numerical simulations. In this first in a series of investigations our initial conditions consists of spherical molecular core obeying the hydrostatic Bonnor-Ebert-Profile with varying degrees of initial rotation. Our simulations cover both the formation of massive disks in which massive stars form as well as low mass disks. We use a customized version of the FLASH code whose AMR technique allows us to follow the formation of a protstellar disk and protostellar core(s) through more than ten orders in density increase while continuously resolving the local Jeans length (i.e. obeying the Truelove criterion, Truelove et al. (1997)). Our numerical simulations also incorporate the energy loss due to molecular line emission in order to obtain a more realistic picture of protostellar core and disk formation. Out initial states model system of mass 168 M_sol and 2.1 M_sol that will form high and low mass stars, respectively. We follow many features such as the development complex shock structures, and the possible fragmentation of the disk. We find that slowly rotating cores (Omega t_ff = 0.1) produce disks in which a strong bar develops but which does not fragment. Faster initial rotation rates (Omega t_ff = 0.2) result in the formation of a ring which may fragment into two star-forming cores. The size of the rings found in our simulated disks agree with the observations of similar systems.

We present the results of a series of numerical simulations of compressible, self-gravitating hydrodynamic turbulence of cluster-forming clumps in molecular clouds. We examine the role that turbulence has in the formation of gravitationally bound cores, studying the dynamical state, internal structure and bulk properties of these cores. Complex structure in turbulent clumps is formed provided that the damping time of the turbulence, $t_{damp}$ is longer than the gravitational free-fall time $t_{ff}$ in a region. We find a variety of density and infall velocity structures among the cores in the simulation, including cores that resemble the Larson-Penston collapse of an isothermal sphere ($\rho \propto r^{-2}$) as well as cores that resemble the McLaughlin-Pudritz collapse of logatropic spheres ($\rho \propto r^{-1}$). The specific angular momentum profiles range between $j \propto r^{1} - r^{2}$. The masses of the bound cores that form are well-fit by the turbulent mass spectrum of Padoan and Nordlund (2002), while the specific angular momentum distribution can be fit by a broken power law. While our hydrodynamic simulations reproduce many of the observed properties of cores, we find an upper limit for the star formation efficiency (SFE) in clusters of 40-50 per cent.

We critically examine the theory of disk locking, which assumes that the angular momentum deposited on an accreting protostar is exactly removed by torques carried along magnetic field lines connecting the star to the disk. In this letter, we consider that the differential rotation between the star and disk naturally leads to an opening (i.e., disconnecting) of the magnetic field between the two. We find that this significantly reduces the spin-down torque on the star by the disk. Thus, disk-locking cannot account for the slow rotation ($\sim$10% of breakup speed) observed in several systems and for which the model was originally developed.

In this talk, we take a new look at the theory of disk locking, which assumes that an accreting protostar rids itself of accreted angular momentum through a magnetic coupling with the accretion disk. We consider that differential rotation between the star and disk twists the field lines. For large enough twist, the magnetic field lines connecting the star and disk open and disconnect. This significantly reduces the spin-down torque on the star by the disk, and so we find that disk-locking theory predicts spin periods that are much too short to account for typical observed systems.

The final masses of Jovian planets are attained when the tidal torques that they exert on their surrounding protostellar disks are sufficient to open gaps in the face of disk viscosity, thereby shutting off any further accretion. In sufficiently well-ionized disks, the predominant form of disk viscosity originates from the Magneto-Rotational Instability (MRI) that drives hydromagnetic disk turbulence. In the region of sufficiently low ionization rate -- the so-called dead zone -- turbulence is damped and we show that lower mass planets will be formed. We considered three ionization sources (X-rays, cosmic rays, and radioactive elements) and determined the size of a dead zone for the total ionization rate by using a radiative, hydrostatic equilibrium disk model developed by Chiang et al. (2001). We studied a range of surface mass density (Sigma_{0}=10^3 - 10^5 g cm^{-2}) and X-ray energy (kT_{x}=1 - 10 keV). We also compared the ionization rate of such a disk by X-rays with cosmic rays and find that the latter dominate X-rays in ionizing protostellar disks unless the X-ray energy is very high (5 - 10 keV). Among our major conclusions are that for typical conditions, dead zones encompass a region extending out to several AU -- the region in which terrestrial planets are found in our solar system. Our results suggest that the division between low and high mass planets in exosolar planetary systems is a consequence of the presence of a dead zone in their natal protoplanetary disks. We also find that the extent of a dead zone is mainly dependent on the disk's surface mass density. Our results provide further support for the idea that Jovian planets in exosolar systems must have migrated substantially inwards from their points of origin.

We develop models for the self-similar collapse of magnetized isothermal cylinders. We find solutions for the case of a fluid with a constant toroidal flux-to-mass ratio (Gamma_phi=constant) and the case of a fluid with a constant gas to magnetic pressure ratio (beta=constant). In both cases, we find that a low magnetization results in density profiles that behave as rho ~ r^{-4} at large radii, and at high magnetization we find density profiles that behave as rho ~ r^{-2}. This density behaviour is the same as for hydrostatic filamentary structures, suggesting that density measurements alone cannot distinguish between hydrostatic and collapsing filaments--velocity measurements are required. Our solutions show that the self-similar radial velocity behaves as v_r ~ r during the collapse phase, and that unlike collapsing self-similar spheres, there is no subsequent accretion (i.e. expansion-wave) phase. We also examine the fragmentation properties of these cylinders, and find that in both cases, the presence of a toroidal field acts to strengthen the cylinder against fragmentation. Finally, the collapse time scales in our models are shorter than the fragmentation time scales. Thus, we anticipate that highly collapsed filaments can form before they are broken into pieces by gravitational fragmentation.

We have performed fully 3D simulations of the collapse of molecular cloud cores which obey the logatropic equation of state. By following the collapse of these cores from states of near hydrodynamic equilibrium, we are able to produce accretion histories which closely resemble those of observed cores. The accretion proceeds in four distinct stages: an initial period of very slow accretion; a period of vigorous accretion following the development of a central density singularity, with a mass accretion rate proportional to t^3, as predicted by self-similar models; a period of relatively stable, vigorous accretion; and finally a gradual decrease in the accretion rate once about 50% of the available mass of the molecular cloud core has been accreted. These results may explain the accretion histories of cores as they pass through the pre-protostellar, Class 0, and Class I stages.

We performed simulations of self-gravitating hydrodynamic turbulence to model the formation of filaments, clumps and cores in molecular clouds. We find that when the mass on the initial computational grid is comparable to the Jeans mass, turbulent pressure is able to prevent gravitational collapse. When the turbulence has damped away sufficiently, gravitational collapse can occur, and the resulting structure closely resembles the pre-singularity collapse of an isothermal sphere of Penston (1969). If several Jeans masses are initially placed on the grid, turbulence may not be sufficient to prevent collapse before turbulence can be significantly damped. In this case, the cores have density structures which are considerably shallower than expected for an isothermal gas, and resemble the solutions for a logatropic equation of state.

We present the extension of previous two-dimensional simulations of the time-dependent evolution of non-relativistic outflows from the surface of Keplerian accretion disks, to three dimensions. The accretion disk itself is taken to provide a set of fixed boundary conditions for the problem. The 3-D results are consistent with the theory of steady, axisymmetric, centrifugally driven disk winds up to the Alfv\'en surface of the outflow. Beyond the Alfv\'en surface however, the jet in 3-D becomes unstable to non-axisymmetric, Kelvin-Helmholtz instabilities. We show that jets maintain their long-term stability through a self-limiting process wherein the average Alfv\'enic Mach number within the jet is maintained to order unity. This is accomplished in at least two ways. First, poloidal magnetic field is concentrated along the central axis of the jet forming a ``backbone'' in which the Alfv\'en speed is sufficiently high to reduce the average jet Alfv\'enic Mach number to unity. Second, the onset of higher order Kelvin-Helmholtz ``flute'' modes (m \ge 2) reduce the efficiency with which the jet material is accelerated, and transfer kinetic energy of the outflow into the stretched, poloidal field lines of the distorted jet. This too has the effect of increasing the Alfv\'en speed, and thus reducing the Alfv\'enic Mach number. The jet is able to survive the onset of the more destructive m=1 mode in this way. Our simulations also show that jets can acquire corkscrew, or wobbling types of geometries in this relatively stable end-state, depending on the nature of the perturbations upon them. Finally, we suggest that jets go into alternating periods of low and high activity as the disappearance of unstable modes in the sub-Alfv\'enic regime enables another cycle of acceleration to super-Alfv\'enic speeds.

We present the results of fully 3-D hydrodynamic simulations of the gravitational collapse of isolated, turbulent molecular cloud cores. Starting from initial states of hydrostatic equilibrium, we follow the collapse of both singular and nonsingular logatropic cores until the central protostar has accreted > 90% of the total available mass. We find that, in the collapse of a singular core with access to a finite mass reservoir, the mass of the central protostar increases as M_acc proportional to t^4 until it has accreted about 35% of the total available mass. For nonsingular cores of fiducial masses 1, 2.5, and 5 M_solar, we find that protostellar accretion proceeds slowly prior to the formation of a singular density profile. Immediately thereafter, the accretion rate in each case increases to about 10^{-6} M_solar/yr, for cores with central temperature T_c= 10 K and truncation pressure P_s = 1.3E5 k_B K/cm^3. It remains at that level until half the available mass has been accreted. After this point, the accretion rate falls steadily as the remaining material is accreted onto the growing protostellar core. We suggest that this general behaviour of the protostellar accretion rate may be indicative of evolution from the Class 0 to the Class I protostellar phase.

In an exploration of the birthplaces of globular clusters, we present a careful examination of the formation of self-gravitating gas clouds within assembling dark matter haloes in a hierarchical cosmological model. Our high-resolution smoothed particle hydrodynamical simulations are designed to determine whether or not hypothesized supergiant molecular clouds (SGMCs) form and, if they do, to determine their physical properties and mass spectra. It was suggested in earlier work that clouds with a median mass of several 10^8 M_sun are expected to assemble during the formation of a galaxy, and that globular clusters form within these SGMCs. Our simulations show that clouds with the predicted properties are indeed produced as smaller clouds collide and agglomerate within the merging dark matter haloes of our cosmological model. We find that the mass spectrum of these clouds obeys the same power-law form observed for globular clusters, molecular clouds, and their internal clumps in galaxies, and predicted for the supergiant clouds in which globular clusters may form. We follow the evolution and physical properties of gas clouds within small dark matter haloes up to z = 1, after which prolific star formation is expected to occur. Finally, we discuss how our results may lead to more physically motivated "rules" for star formation in cosmological simulations of galaxy formation.

We model the sub-millimetre polarization patterns that are expected for filamentary clouds that are threaded by helical magnetic fields. We previously developed a three parameter model of such clouds (Fiege & Pudritz 2000a), which are described by a concentration parameter $C$, and two flux to mass ratios $\Gz$ and $\Gphi$ to specify the mass loading of the poloidal and toroidal field lines respectively. Our models provide a simple and purely geometric explanation for the well-known ``polarization hole'' effect, in which the sub-millimetre polarization percentage decreases toward the regions of peak intensity. This occurs because of a cancellation between contributions to the polarization from the ``backbone'' of poloidal flux along the filament's axis and its surrounding envelope, which is dominated by the toroidal field component. A systematic exploration of our parameter space allows us to classify the polarization patterns due to filaments aligned approximately perpendicular to the plane of the sky into three basic types. The polarization vectors are parallel to filaments when $\ratio\appleq 0.1$, where $\Bzs$ and $\Bphis$ are respectively the poloidal and toroidal magnetic field components at the outer surface of the filament. The polarization vectors are perpendicular to filaments when $\ratio\appgeq 0.33$. Intermediate cases result in polarization patterns that contain $90^\circ$ flips in the orientation of the polarization vectors. The flips are symmetric about the central axis for filaments oriented parallel to the plane of the sky, but more complicated asymmetric patters result from filaments that are inclined at some angle.

We present a new model of molecular cloud cores that originate from filamentary clouds that are threaded by helical magnetic fields. Only modest toroidal fields are required to produce elongated cores, with intrinsic axis ratios in the range 0.1 <= q <= 0.9 and mean projected axis ratios in the range 0.3 <= <p> <= 1. Thus many of our models are in good agreement with the observed shapes of cores (Myers et al 1991, Ryden 1996), which find axis ratios distributed about the mean values <p> \approx 0.5-0.6 and q \approx 0.4-0.5. We find that the Bonnor-Ebert critical mass is reduced by about 20% by the helical field threading our models. We also find that models are generally most elongated when the mass is significantly less than the Bonnor-Ebert critical mass for unmagnetized cores.

We develop a theory for filamentary molecular clouds including the effects of ordered magnetic fields, and external pressure. We first derive a form of the virial equation appropriate for filamentary clouds. By comparing with observational results collected from the literature, we find that the fields are likely helical. Secondly, we construct numerical, MHD models of filamentary clouds that agree with the observational constraints. We find that our models produce much more realistic density profiles $r\sim r^{-1.8 to -2}$ than previous models, where the density falls off as $r^{-4}$.

We review the physics of star formation, and its links with the state of the ISM in galaxies. Current obervations indicate that the preferred mode of star formation is clustered. Given that OB associations provide the dominant energy input into the ISM, deep links exist between the ISM and star formation. We present a multi-scale discussion of star formation, and attempt to create an integrated vision of these processes.

We review recent observational and theoretical results on the relationship between circumstellar accretion disks and jets in young stellar objects. We then present a theoretical framework that interprets jets as accretion-powered, centrifugally driven winds from magnetized accretion disks. Recent progress in the numerical simulation of such outflows is described. We also discuss the structure of the underlying magnetized protostellar disks, emphasizing the role that large-scale, open magnetic fields can play in angular momentum transport.

Young stars form on a wide range of scales, producing aggregates and clusters with various degrees of gravitational self-binding. The loose aggregates have a hierarchical structure in both space and time that resembles interstellar turbulence, suggesting that these stars form in only a few turbulent crossing times with positions that map out the previous gas distribution. Dense clusters, on the other hand, are often well mixed, as if self-gravitational motion has erased the initial fine structure. Nevertheless, some of the youngest dense clusters also show sub-clumping, so it may be that all stellar clustering is related to turbulence. Some of the densest clusters may also be triggered. The evidence for mass segregation of the stars inside clusters is reviewed, along with various explanations for this effect. Other aspects of the theory of cluster formation are reviewed as well, including many specific proposals for cluster formation mechanisms. The conditions for the formation of bound clusters are discussed. Critical star formation efficiencies can be as low as 10% if the gas removal process is slow and the stars are born at sub-virial speeds. Environmental conditions, particularly pressure, may affect the fraction and masses of clusters that end up bound. Globular clusters may form like normal open clusters but in conditions that prevailed during the formation of the halo and bulge, or in interacting and starburst galaxies today. Various theories for the formation of globular clusters are summarized.

In Paper I (Fiege & Pudritz, 1999), we constructed models of filamentary molecular clouds that are truncated by a realistic external pressure and contain a rather general helical magnetic field. We address the stability of our models to gravitational fragmentation and axisymmetric MHD-driven instabilities. By calculating the dominant modes of axisymmetric instability, we determine the dominant length scales and growth rates for fragmentation. We find that the role of pressure truncation is to decrease the growth rate of gravitational instabilities by decreasing the self-gravitating mass per unit length. Purely poloidal and toroidal fields also help to stabilize filamentary clouds against fragmentation. The overall effect of helical fields is to stabilize gravity-driven modes, so that the growth rates are significantly reduced below what is expected for unmagnetized clouds. However, MHD ``sausage'' instabilities are triggered in models whose toroidal flux to mass ratio exceeds the poloidal flux to mass ratio by more than a factor of $\sim 2$. We find that observed filaments appear to lie in a physical regime where the growth rates of both gravitational fragmentation and axisymmetric MHD-driven modes are at a minimum.

We study the equilibrium of pressure truncated, filamentary molecular clouds that are threaded by rather general helical magnetic fields. We first derive a new virial equation appropriate for magnetized filamentary clouds, which includes the effects of non-thermal motions and the turbulent pressure of the surrounding ISM. When compared with the data, we find that many filamentary clouds have a mass per unit length that is significantly reduced by the effects of external pressure, and that toroidal fields play a significant role in squeezing such clouds. We also develop exact numerical MHD models of filamentary molecular clouds with more general helical field configurations than have previously been considered. We also examine the effects of the equation of state by comparing ``isothermal'' filaments, with constant total (thermal plus turbulent) velocity dispersion, with equilibria constructed using a logatropic equation of state. We perform a Monte Carlo exploration of our parameter space to determine which choices of parameters result in models that agree with the available observational constraints. We find that both equations of state result in equilibria that agree with the observational results. Moreover, we find that models with helical fields have more realistic density profiles than either unmagnetized models or those with purely poloidal fields; we find that most isothermal models have density distributions that fall off as r^{-1.8} to r^{-2}, while logatropes have density profiles that range from r^{-1} to r^{-1.8}. We find that purely poloidal fields produce filaments with steep density gradients that not allowed by the observations.

New wide-field images of the Rho Ophiuchus molecular cloud at 850 and 450 microns obtained with SCUBA on the James Clerk Maxwell Telescope reveal a wide variety of large-scale features that were previously unknown. Two linear features each 4' (0.2 pc) in length extend to the north of the bright emission region containing SM1 and VLA 1623. These features may correspond to the walls of a previously unidentified outflow cavity, or the boundary of a photon dominated region powered by a nearby B star. A previously unidentified source is located in the north-east region of the image. The properties of this source (diameter ~5000 AU, mass ~0.3-1 Msun) suggest that it is a pre-protostellar core. Two arcs of emission are seen in the direction of the north-west extension of the VLA 1623 outflow. The outer arc appears relatively smooth at 850 microns and is estimated to have a mass of ~0.3 Msun, while the inner arc breaks up into a number of individual clumps, some of which are previously identified protostars.

Observational advances over the last decade reveal that star formation is associated with the simultaneous presence of gravitationally collapsing gas, bipolar outflow, and an accretion disk. Two theoretical views of star formation suppose that either stellar mass is determined from the outset by gravitational instability, or by the outflow which sweeps away the collapsing envelope of initially singular density distributions. Neither picture appears to explain all of the facts. This contribution examines some of the key issues facing star formation theory.

We present semi-analytical similarity solutions for the inside-out, expansion-wave collapse of initially virialized gas clouds with non-isothermal equations of state. Results are given for the family of negative-index polytropes, but we focus especially on the so-called logotrope, P/P_c=1+A ln(rho/rho_c). The formalism and interpretation of the present theory are extensions of those in Shu's (1977) standard model for accretion in self-gravitating isothermal spheres: a collapse front moves outwards into a cloud at rest, and the gas behind it falls back to a collapsed core, or protostar. The infalling material eventually enters free-fall, so that, at small radii, the density profiles and velocity fields have the same power-law forms in logotropic and isothermal spheres both. However, the accretion rate onto a protostar is not constant in a logotrope, but grows in proportion to t^3 during the expansion wave. Thus, low-mass stars are accreted over longer times, and high-mass stars over shorter times, than expected in isothermal clouds. This result has implications for the form and origin of the stellar IMF. We also find that infall velocities grow more slowly with time in a collapsing logotrope. These results lead to older inferred collapse ages for Class 0 protostars in general, and for the Bok globule B335 in particular.

Continuum emission at 220 and 355 GHz from the prototype Class 0 source VLA 1623 has been detected using the JCMT-CSO interferometer. Gaussian fits to the data place an upper limit of 70 AU on the half-width half-maximum radius of the emission, which implies an upper limit of ~175 AU for the cutoff radius of the circumstellar disk in the system. In the context of existing collapse models, this disk could be magnetically supported on the largest scales and have an age of ~6x10^4 yr, consistent with previous suggestions that Class 0 sources are quite young. The innermost region of the disk within ~6 AU is likely to be in centrifugal support, which is likely large enough to provide a drive for the outflow according to current theoretical models. Alternatively, if 175 AU corresponds to the centrifugal radius of the disk, the age of the system is ~2x10^5 yr, closer to age estimates for Class I sources.

We generalize the classic Bonnor-Ebert stability analysis of pressure-truncated, self-gravitating gas spheres, to include clouds with arbitrary equations of state. A virial-theorem analysis is also used to incorporate mean magnetic fields into such structures. The results are applied to giant molecular clouds (GMCs), and to individual dense cores, with an eye to accounting for recent observations of the internal velocity-dispersion profiles of the cores in particular. We argue that GMCs and massive cores are at or near their critical mass, and that in such a case the size-linewidth and mass-radius relations between them are only weakly dependent on their internal structures; any gas equation of state leads to essentially the same relations. We briefly consider the possibility that molecular clouds can be described by polytropic pressure-density relations (of either positive or negative index), but show that these are inconsistent with the apparent gravitational virial equilibrium, 2U + W = 0 of GMCs and of massive cores. This class of models would include clouds whose nonthermal support comes entirely from Alfven wave pressure. The simplest model consistent with all the salient features of GMCs and cores is a ``pure logotrope,'' in which P/P_c = 1 + A ln(rho/rho_c). Detailed comparisons with data are made to estimate the value of A, and an excellent fit to the observed dependence of velocity dispersion on radius in cores is obtained with A = 0.2.

We investigate the global stability of a differentially rotating fluid shell threaded by vertical and azimuthal magnetic fields to linear, axisymmetric perturbations. This system, which models a thick accretion disk in the vicinity of its midplane, is susceptible to the Velikhov-Chandrasekhar (VC) instability in the absence of the azimuthal field. In most cases, the azimuthal field tends to stabilize the VC instability, although strong fields (Alfv\'en speed of order the characteristic rotational speed in our incompressible model) are required for complete stabilization. Stability diagrams are constructed, indicating critical values of the two fields for instability. We find an additional strong field instability that arises when the azimuthal Alfv\'en speed exceeds the characteristic rotational speed. This instability, in the case of a freely bounded configuration, has certain similarities to the sausage instability for interpenetrating fields in plasma physics, and may be important for very massive disks or filamentary molecular clouds. An application to the L1641 region in Orion A is briefly discussed. Finally, we find that the effect of a radially varying vertical field (without an azimuthal field) is mainly stabilizing.