A compact object moving relative to surrounding gas accretes material and perturbs the density of gas in its vicinity. In the classical picture of Bondi-Hoyle-Lyttleton accretion, the perturbation takes the form of an overdense wake behind the object, which exerts a dynamical friction drag. We use hydrodynamic simulations to investigate how the accretion rate and strength of dynamical friction are modified by the presence of outflow from the compact object. We show that the destruction of the wake by an outflow reduces dynamical friction, and reverses its sign when the outflow is strong enough, in good quantitative agreement with analytic calculations. For a strong isotropic outflow, the outcome on scales that we have simulated is a negative dynamical friction, i.e., net acceleration. For jet-like outflows driven by reprocessed accretion, both the rate of accretion and the magnitude of dynamical friction drop for more powerful jets. The accretion rate is strongly intermittent when the jet points to the same direction as the motion of the compact object. The dynamical effects of outflows may be important for the evolution of compact objects during the common envelope phase of binary systems, and for accreting compact objects and massive stars encountering AGN discs.

The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud.

Recent observations of global velocity gradients across and along molecular filaments have been interpreted as signs of gas accreting onto and along these filaments, potentially feeding star-forming cores and proto-clusters. The behavior of velocity gradients in filaments, however, has not been studied in detail, particularly on small scales (< 0.1 pc). In this paper, we present MUFASA, an efficient, robust, and automatic method to fit ammonia lines with multiple velocity components, generalizable to other molecular species. We also present CRISPy, a Python package to identify filament spines in 3D images (e.g., position-position-velocity cubes), along with a complementary technique to sort fitted velocity components into velocity-coherent filaments. In NGC 1333, we find a wealth of velocity gradient structures on a beam-resolved scale of ~0.05 pc. Interestingly, these local velocity gradients are not randomly oriented with respect to filament spines and their perpendicular, i.e., radial, component decreases in magnitude towards the spine for many filaments. Together with remarkably constant velocity gradients on larger scales along many filaments, these results suggest a scenario in which gas falling onto filaments is progressively damped and redirected to flow along these filaments.

At about 70 solar masses, the recently-discovered dark object orbited by a B-type star in the system LB-1 is difficult to understand as the end point of standard stellar evolution, except as a binary black hole (BBH). LB-1 shows a strong, broad H-alpha emission line that is best attributed to a gaseous disk surrounding the dark mass. We use the observed H-alpha line shape, particularly its wing extension, to constrain the inner radius of the disk and thereby the separation of a putative BBH. The hypothesis of a current BBH is effectively ruled out on the grounds that its merger time must be a small fraction of the current age of the B star. The hypothesis of a previous BBH that merged to create the current dark mass is also effectively ruled out by the low orbital eccentricity, due to the combination of mass loss and kick resulted from gravitational wave emission in any past merger. We conclude that the current dark mass is a single black hole produced by the highly mass-conserving, monolithic collapse of a massive star.

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.

The global climate crisis poses new risks to humanity, and with them, new challenges to the practices of professional astronomy. Avoiding the more catastrophic consequences of global warming by more than 1.5 degrees requires an immediate reduction of greenhouse gas emissions. According to the 2018 United Nations Intergovernmental Panel report, this will necessitate a 45% reduction of emissions by 2030 and net-zero emissions by 2050. Efforts are required at all levels, from the individual to the governmental, and every discipline must find ways to achieve these goals. This will be especially difficult for astronomy with its significant reliance on conference and research travel, among other impacts. However, our long-range planning exercises provide the means to coordinate our response on a variety of levels. We have the opportunity to lead by example, rising to the challenge rather than reacting to external constraints. We explore how astronomy can meet the challenge of a changing climate in clear and responsible ways, such as how we set expectations (for ourselves, our institutions, and our granting agencies) around scientific travel, the organization of conferences, and the design of our infrastructure. We also emphasize our role as reliable communicators of scientific information on a problem that is both human and planetary in scale.

We study the hydrodynamic evolution of an explosion close to the stellar surface, and give predictions for the radiation from such an event. We show that such an event will give rise to a multi-wavelength transient. We apply this model to describe a precursor burst to the peculiar supernova iPTF14hls, which occurred in 1954, sixty year before the supernova. We propose that the new generation of optical surveys might detect similar transients, and they can be used to identify supernova progenitors well before the explosion.

We present the discovery and the photometric and spectroscopic study of H-rich Type II supernova (SN) KSP-SN-2016kf (SN2017it) observed in the KMTNet Supernova Program in the outskirts of a small irregular galaxy at $z\simeq0.043$ within a day from the explosion. Our high-cadence, multi-color ($BVI$) light curves of the SN show that it has a very long rise time ($t_\text{rise}\simeq 20$ days in $V$ band), a moderately luminous peak ($M_V\simeq -$17.6 mag), a notably luminous and flat plateau ($M_V\simeq -$17.4 mag and decay slope $s\simeq0.53$ mag per 100 days), and an exceptionally bright radioactive tail. Using the color-dependent bolometric correction to the light curves, we estimate the $^{56}$Ni mass powering the observed radioactive tail to be $0.10\pm0.01$ M$_\odot$, making it a H-rich Type II SN with one of the largest $^{56}$Ni masses observed to date. The results of our hydrodynamic simulations of the light curves constrain the mass and radius of the progenitor at the explosion to be $\sim$15 M$_\odot$ (evolved from a star with an initial mass of $\sim$ 18.8 M$_\odot$) and $\sim1040$ R$_\odot$, respectively, with the SN explosion energy of $\sim 1.3\times 10^{51}$ erg s$^{-1}$. The above-average mass of the KSP-SN-2016kf progenitor, together with its low metallicity $ Z/Z_\odot \simeq0.1-0.4$ obtained from spectroscopic analysis, is indicative of a link between the explosion of high-mass red supergiants and their low-metallicity environment. The early part of the observed light curves shows the presence of excess emission above what is predicted in model calculations, suggesting there is interaction between the ejecta and circumstellar material. We further discuss the implications of the high progenitor initial mass and low-metallicity environment of KSP-SN-2016kf on our understanding of the origin of Type II SNe.

The virial ratio between kinetic and gravitational terms provides key insight into the balance of forces that confine a molecular cloud, but the clumpy and filamentary structures of resolved clouds make it difficult to evaluate this ratio in a consistent way. For clouds with resolved maps of column density as well as a line tracer, we demonstrate that the gravitational energy can be estimated directly from observations in a manner similar to the kinetic energy. This offers improved diagnostic power and consistency. Disentangling a cloud from foreground and background materials is a persistent challenge, for which we introduce a strategy based on Abel's transform. We provide proofs of principle using simulated clouds.

An over-dense wake is created by a gravitating object moving through a gaseous medium, and this wake pulls back on the object and slows it down. This is conventional dynamical friction in a gaseous medium. We argue that if the object drives a sufficiently powerful outflow, the wake is destroyed and instead an extended under-dense region is created behind the object. In this case the overall gravitational force is applied in the direction of the object's motion, producing a negative dynamical friction (NDF). Black holes in dense gas drive powerful outflows and may experience the NDF, although extensive numerical work is probably needed to demonstrate or refute this conclusively. NDF may be important for stellar-mass black holes and neutron stars inside `common envelopes' in binary systems, for stellar mass black holes inside AGN discs, or for massive black holes growing through super-Eddington accretion in early Universe.

We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transition to coherence" in velocity dispersion around its periphery. The identification of these structures provides a chance to study the coherent structures in molecular clouds statistically. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 Msun, generally smaller than previously known coherent cores identified by Goodman et al. (1998), Caselli et al. (2002), and Pineda et al. (2010). We call these structures "droplets." We find that unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor-Ebert sphere, and they have a velocity (VLSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process.

Dust barriers effectively capture the photon momentum of a central light source, but low-density channels, along with re-emission at longer wavelengths, enhance its escape. We use Monte Carlo simulations to study the effects of inhomogeneity on radiation forces imparted to a dust envelope around a central star. We survey the strength and scale of an inhomogeneous perturbation field, as well as the optical depth of its spherical reference state. We run at moderate numerical resolution, relying on our previous resolution study for calibration of the associated error. We find that inhomogeneities matter most when their scale exceeds the characteristic mean free path. As expected, they tend to reduce the net radiation force and extend its range; however, there is significant variance among realizations. Within our models, force integrals correlate with the emergent spectral energy distribution, given a specified set of dust properties. A critical issue is the choice of integral measures of the radiation force: for strong deviations from spherical symmetry the relevant measures assess radial forces relative to the cloud centre of mass. Of these, we find the virial term due to radiation to be the least stochastic of several integral measures in the presence of inhomogeneities.

We address in detail the radiation forces on spherical dust envelopes around luminous stars, and numerical solutions for these forces, as a first step toward more general dust geometries. Two physical quantities, a normalized force and a force-averaged radius, suffice to capture the overall effects of radiation forces; these combine to yield the radiation term in the virial theorem. In addition to the optically thin and thick regimes, the wavelength dependence of dust opacity allows for an intermediate case in which starlight is easily trapped but infrared radiation readily escapes. Scattering adds a non-negligible force in this intermediate regime. We address all three regimes analytically and provide approximate formulae for the force parameters, for arbitrary optical depth and inner dust temperature. Turning to numerical codes, we examine the convergence properties of the Monte Carlo code Hyperion run in Cartesian geometry. We calibrate both Hyperion and our analytical estimate using the DUSTY code, run in spherical geometry. We find that Monte Carlo codes tend to underestimate the radiation force when the mean free path of starlight is not well resolved, as this causes the inner dust temperature, and therefore the inner Rosseland opacity, to be too low. We briefly discuss implications for more complicated radiation transfer problems.

Early light from core-collapse supernovae, now detectable in high-cadence surveys, holds clues to a star and its environment just before it explodes. However, effects that alter the early light have not been fully explored. We highlight the possibility of non-radial flows at the time of shock breakout. These develop in sufficiently non-spherical explosions if the progenitor is not too diffuse. When they do develop, non-radial flows limit ejecta speeds and cause ejecta-ejecta collisions. We explore these phenomena and their observational implications, using global, axisymmetric, non-relativistic FLASH simulations of simplified polytropic progenitors, which we scale to representative stars. We develop a method to track photon production within the ejecta, enabling us to estimate band-dependent light curves from adiabatic simulations. Immediate breakout emission becomes hidden as an oblique flow develops. Non-spherical effects lead the shock-heated ejecta to release a more constant luminosity at a higher, evolving color temperature at early times, effectively mixing breakout light with the early light curve. Collisions between non-radial ejecta thermalize a small fraction of the explosion energy; we address emission from these collisions in a subsequent paper.

Star cluster formation is unlikely to be a sudden event: instead, matter will flow to a cluster's formation site over an extended period, even as stars form and inject energy to the region. A cluster's gaseous precursor must persist under the competing influences of accretion and feedback for several crossing times, insofar as star formation is a slow process. The new-born stellar cluster should therefore preserve a memory of this competition. Using analytical approximations we assess the dynamical state of the gas, mapping regimes in which various types of feedback are weak or strong. Protostellar outflows, radiation pressure, and ionized gas pressure are accounted for. Comparison to observations shows that feedback is often incapable of expelling gas in the more massive, rapidly-accreting clusters, but feedback may nevertheless starve accretion by acting on larger scales.

Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe Nebula, to achieve a comprehensive view of the kinematic properties of the source. These information were derived by simultaneously fitting ammonia inversion transition lines (1,1) and (2,2). Our analysis unveils the imprint of protostellar feedback, such as increasing line widths, temperature and turbulent motions in our molecular data. Combined with complementary observations of dust thermal emission, we estimate that the core is gravitationally bound following a virial analysis. If the core is not contracting, another source of internal pressure, most likely the magnetic field, is supporting it against gravitational collapse and limits its star formation efficiency.

We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores.

We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure confined.

We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH$_3$ emission to dust continuum emission from Herschel, and find that the two tracers correspond closely. NH$_3$ is present in over 60\% of lines-of-sight with $A_V \gtrsim 7$ mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH$_3$ is detected toward ~ 40\% of lines-of-sight with $A_V \gtrsim 7$ mag. Moreover, we find that the NH$_3$ emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH$_3$ column densities through forward modeling of the hyperfine structure of the NH$_3$ (1,1) and (2,2) lines. We show that the NH$_3$ velocity dispersion, ${\sigma}_v$, and gas kinetic temperature, $T_K$, vary systematically between the regions included in this release, with an increase in both the mean value and spread of ${\sigma}_v$ and $T_K$ with increasing star formation activity. The data presented in this paper are publicly available.

Wave-driven outflows and non-disruptive explosions have been implicated in pre-supernova outbursts, supernova impostors, LBV eruptions, and some narrow-line and superluminous supernovae. To model these events, we investigate the dynamics of stars set in motion by strong acoustic pulses and wave trains, focusing here on nonlinear wave propagation, shock formation, and an early phase of the development of a weak shock. We identify the shock formation radius, showing that a heuristic estimate based on crossing characteristics matches an exact expansion around the wave front and verifying both with numerical experiments. Our general analytical condition for shock formation applies to one-dimensional motions within any static environment, including both eruptions and implosions, and can easily be extended to non-stationary flows. We also consider the early phase of shock energy dissipation. We find that waves of super-Eddington acoustic luminosity always create shocks, rather than damping by radiative diffusion. Therefore, shock formation is integral to super-Eddington outbursts.

Hydrostatic models of Wolf-Rayet stars typically contain low-density outer envelopes that inflate the stellar radii by a factor of several and are capped by a denser shell of gas. Inflated envelopes and density inversions are hallmarks of envelopes that become super-Eddington as they cross the iron-group opacity peak, but these features disappear when mass loss is sufficiently rapid. We re-examine the structures of steady, spherically symmetric wind solutions that cross a sonic point at high optical depth, identifying the physical mechanism by which outflow affects the stellar structure, and provide an improved analytical estimate for the critical mass loss rate above which extended structures are erased. Weak-flow solutions below this limit resemble hydrostatic stars even in supersonic zones; however, we infer that these fail to successfully launch optically thick winds. Wolf-Rayet envelopes will therefore likely correspond to the strong, compact solutions. We also find that wind solutions with negligible gas pressure are stably stratified at and below the sonic point. This implies that convection is not the source of variability in Wolf-Rayet stars, as has been suggested; but, acoustic instabilities provide an alternative explanation. Our solutions are limited to high optical depths by our neglect of Doppler enhancements to the opacity, and do not account for acoustic instabilities at high Eddington factors; yet they provide useful insights into Wolf-Rayet stellar structures.

During star cluster formation, ongoing mass accretion is resisted by stellar feedback in the form of protostellar outflows from the low-mass stars and photo-ionization and radiation pressure feedback from the massive stars. We model the evolution of cluster-forming regions during a phase in which both accretion and feedback are present, and use these models to investigate how star cluster formation might terminate. Protostellar outflows are the strongest form of feedback in low-mass regions, but these cannot stop cluster formation if matter continues to flow in. In more massive clusters, radiation pressure and photo-ionization rapidly clear the cluster-forming gas when its column density is too small. We assess the rates of dynamical mass ejection and of evaporation, while accounting for the important effect of dust opacity on photo-ionization. Our models are consistent with the census of protostellar outflows in NGC 1333 and Serpens South, and with the dust temperatures observed in regions of massive star formation. Comparing observations of massive cluster-forming regions against our model parameter space, and against our expectations for accretion-driven evolution, we infer that massive-star feedback is a likely cause of gas disruption in regions with velocity dispersions less than a few kilometers per second, but that more massive and more turbulent regions are too strongly bound for stellar feedback to be disruptive.

We present the first fully calibrated H$_2$, 1-0 S(1) image of the entire 30 Doradus nebula. The observations were conducted using the NOAO Extremely Wide-Field Infrared Imager on the CTIO 4-meter Blanco Telescope. Together with a NEWFIRM Br$\gamma$ image of 30 Doradus, our data reveal the morphologies of the warm molecular gas and ionized gas in 30 Doradus. The brightest H$_2$-emitting area, which extends from the northeast to the southwest of R136, is a photodissociation region viewed face-on, while many clumps and pillar features located at the outer shells of 30 Doradus are photodissociation regions viewed edge-on. Based on the morphologies of H$_2$, Br$\gamma$, $^{12}$CO, and 8$\mu$m emission, the H$_2$ to Br$\gamma$ line ratio and Cloudy models, we find that the H$_2$ emission is formed inside the photodissociation regions of 30 Doradus, 2 - 3 pc to the ionization front of the HII region, in a relatively low-density environment $<$ 10$^4$ cm$^{-3}$. Comparisons with Br$\gamma$, 8$\mu$m, and CO emission indicate that H$_2$ emission is due to fluorescence, and provide no evidence for shock excited emission of this line.

Non-spherical explosions develop non-radial flows as the pattern of shock emergence progresses across the stellar surface. In supernovae these flows can limit ejecta speeds, stifle shock breakout emission, and cause collisions outside the star. Similar phenomena occur in stellar and planetary collisions, tidal disruption events, accretion-induced collapses, and propagating detonations. We present two-dimensional, nested-grid Athena simulations of non-radial shock emergence in a frame comoving with the breakout pattern, focusing on the adiabatic, non-relativistic limit in a plane stratified envelope. We set boundary conditions using a known self-similar solution and explore the role of box size and resolution on the result. The shock front curves toward the stellar surface, and exhibits a kink from which weak discontinuities originate. Flow around the point of shock emergence is neither perfectly steady nor self-similar. Waves and vortices, which are not predominantly due to grid effects, emanate from this region. The post-shock flow is deflected along the stellar surface, and its pressure disturbs the stellar atmosphere upstream of the emerging shock. We use the numerical results and their analytical limits to predict the effects of radiation transfer and gravity, which are not included in our simulations.

During a stellar tidal disruption event (TDE), an accretion disk forms as stellar debris returns to the disruption site and circularizes. Rather than being confined within the circularizing radius, the disk can spread to larger radii to conserve angular momentum. A spreading disk is a source of matter for re-accretion at rates which can exceed the later stellar fallback rate, although a disk wind can suppress its contribution to the central black hole accretion rate. A spreading disk is detectible through a break in the central accretion rate history, or, at longer wavelengths, by its own emission. We model the evolution of TDE disk size and accretion rate, by accounting for the time-dependent fallback rate, for the influence of wind losses in the early, advective stage, and for the possibility of thermal instability for accretion rates intermediate between the advection-dominated and gas-pressure dominated states. The model provides a dynamic basis for modeling TDE light curves. All or part of a young TDE disk will precess as a solid body due to Lense-Thirring effect, and precession may manifest itself as quasi-periodic modulation of light curve. The precession period increases with time. Applying our results to the jetted TDE candidate Swift J1644+57, whose X-ray light curve shows numerous quasi-periodic dips, we argue that the data best fit a scenario in which a main-sequence star was fully disrupted by an intermediate mass black hole on an orbit significantly inclined from the black hole equator, with the apparent jet shutoff at t= 500 d corresponding to a disk transition from the advective state to the gas-pressure dominated state.

In a non-spherical stellar explosion, non-radial motions become important near the stellar surface. For realistic deviations from spherical symmetry, non-radial flow dramatically alters the dynamics and emission of shock emergence on a significant fraction of the surface. The breakout flash is stifled, ejecta speeds are limited, and matter is cast sideways. Non-radial ejection allows for collisions outside the star, which may engender a new type of transient. Strongly oblique breakouts are most easily produced in compact stellar progenitors, such as white dwarfs and stripped-envelope core collapse supernovae. We study the shock structure and post-shock acceleration using conservation laws, a similarity analysis, and an approximate theory for oblique shocks. The shock is likely to extend vertically from the stellar surface, then kink before joining a deep asymptotic solution. Outflow from the region crossed by an oblique shock is probably unsteady and may affect the surface ahead of the main shock. We comment on the implications for several notable explosions in which the non-spherical dynamics described in this paper are likely to play an important role. We also briefly consider relativistic and superluminal pattern speeds.

We examine the dynamics of accelerating normal shocks in stratified planar atmospheres, providing accurate fitting formulae for the scaling index relating shock velocity to the initial density and for the post-shock acceleration factor as functions of the polytropic and adiabatic indices which parameterize the problem. In the limit of a uniform initial atmosphere there are analytical formulae for these quantities. In the opposite limit of a very steep density gradient the solutions match the outcome of shock acceleration in exponential atmospheres.

Optical and infrared emission lines from HII regions are an important diagnostic used to study galaxies, but interpretation of these lines requires significant modeling of both the internal structure and dynamical evolution of the emitting regions. Most of the models in common use today assume that HII region dynamics are dominated by the expansion of stellar wind bubbles, and have neglected the contribution of radiation pressure to the dynamics, and in some cases also to the internal structure. However, recent observations of nearby galaxies suggest that neither assumption is justified, motivating us to revisit the question of how HII region line emission depends on the physics of winds and radiation pressure. In a companion paper we construct models of single HII regions including and excluding radiation pressure and winds, and in this paper we describe a population synthesis code that uses these models to simulate galactic collections of HII regions with varying physical parameters. We show that the choice of physical parameters has significant effects on galactic emission line ratios, and that in some cases the line ratios can exceed previously claimed theoretical limits. Our results suggest that the recently-reported offset in line ratio values between high-redshift star-forming galaxies and those in the local universe may be partially explained by the presence of large numbers of radiation pressured-dominated HII regions within them.

The emission line ratios [OIII]5007/H-beta and [NII]6584/H-alpha have been adopted as an empirical way to distinguish between the fundamentally different mechanisms of ionization in emission-line galaxies. However, detailed interpretation of these diagnostics requires calculations of the internal structure of the emitting HII regions, and these calculations depend on the assumptions one makes about the relative importance of radiation pressure and stellar winds. In this paper we construct a grid of quasi-static HII region models to explore how choices about these parameters alter HII regions' emission line ratios. We find that, when radiation pressure is included in our models, HII regions reach a saturation point beyond which further increases in the luminosity of the driving stars does not produce any further increase in effective ionization parameter, and thus does not yield any further alteration in an HII region's line ratio. We also show that, if stellar winds are assumed to be strong, the maximum possible ionization parameter is quite low. As a result of this effect, it is inconsistent to simultaneously assume that HII regions are wind-blown bubbles and that they have high ionization parameters; some popular HII region models suffer from this inconsistency. Our work in this paper provides a foundation for a companion paper in which we embed the model grids we compute here within a population synthesis code that enables us to compute the integrated line emission from galactic populations of HII regions.

The Milky Way's most luminous, young and massive (M > 10^4 Msun) star clusters and OB associations have largely evaded detection despite knowledge of their surrounding H II regions. We search for these clusters and associations within the 40 star forming complexes from Rahman & Murray in the 13 most luminous WMAP free-free emission sources of the Galaxy. Selecting for objects with the dust-reddened colors of OB stars, we identify new candidate associations using the 2MASS point source catalog. In 40 star forming complexes searched, 22 contain cluster/association candidates with sizes and masses in the range of 3' - 26' and 10^{2.3} - 10^{5} Msun. Of the 22 candidates, at least 7 have estimated masses > 10^4 Msun, doubling the number of such massive clusters known in the Galaxy. Applying our method to a statistically similar set of test locations, we estimate that 3.0 +/- 0.6 of our 22 candidate associations are unrelated to the star forming complexes. In addition, the apparent extinctions of our candidate associations correlate well with the predictions from a Galactic model. These facts, along with the clear detection of a known OB association and the previous spectral verification of one cluster found by this method, validate our method. Only one of the searched WMAP sources remains without a candidate. In 8 of the most luminous WMAP sources, the candidate associations can account for the observed free-free flux. With our new compilation, the Galactic census of young, massive stellar associations may now be about half complete.

The ionization parameter U is potentially useful for measuring radiation pressure feedback from massive star clusters, as it reflects the radiation-to-gas-pressure ratio and is readily derived from mid-infrared line ratios. We consider several effects which determine the apparent value of U in HII regions and galaxies. An upper limit is set by the compression of gas by radiation pressure. The pressure from stellar winds and the presence of neutral clumps both reduce U for a given radiation intensity. The most intensely irradiated regions are selectively dimmed by internal dust absorption of ionizing photons, inducing observational bias on galactic scales. We explore these effects analytically and numerically, and use them to interpret previous observational results. We find that radiation confinement sets the upper limit log_10 U = -1 seen in individual regions. Unresolved starbursts display a maximum value of ~ -2.3. While lower, this is also consistent with a large portion of their HII regions being radiation dominated, given the different technique used to interpret unresolved regions, and given the bias caused by dust absorption. We infer that many individual, strongly illuminated regions cannot be dominated by stellar winds, and that even when averaged on galactic scales, shocked wind pressures cannot be large compared to radiation pressure. Therefore, most HII regions cannot be adiabatic wind bubbles. Our models imply a metallicity dependence in the physical structure and dust attenuation of radiation-dominated regions, both of which should vary strongly across a critical metallicity of about one-twentieth solar.

Young OB associations with masses greater than 10$^4$ M$_{\sun}$ have been inferred to exist in the Galaxy but have largely evaded detection. Recently, a candidate OB association has been identified within the most luminous star forming complex in the Galaxy, the Dragonfish Nebula. We identify 18 young, massive stars with near-infrared spectroscopy from a sample of 50 members within the candidate OB association, including 15 O-type, and three Luminous Blue Variables or Wolf-Rayet stars. This number matches the expected yield of massive stars from the candidate association, confirming its existence and ability to power the parent star forming complex. These results demonstrate the existence of a 10$^5$ M$_{\sun}$ OB association, more powerful than any previously known in the Galaxy, comparable in mass only to Westerlund 1. Further, the results also validate the color selection method used to identify the association, adding credence to others discovered in the same way.

The plateaus observed in about one half of the early X-ray afterglows are the most puzzling feature in gamma-ray bursts (GRBs) detected by Swift. By analyzing the temporal and spectral indices of a large X-ray plateau sample, we find that 55% can be explained by external, forward shock synchrotron emission produced by a relativistic ejecta coasting in a \rho ~ r^{-2}, wind-like medium; no energy injection into the shock is needed. After the ejecta collects enough medium and transitions to the adiabatic, decelerating blastwave phase, it produces the post-plateau decay. For those bursts consistent with this model, we find an upper limit for the initial Lorentz factor of the ejecta, \Gamma_0 \leq 46 (\epsilon_e/0.1)^{-0.24} (\epsilon_B/0.01)^{0.17}; the isotropic equivalent total ejecta energy is E_{iso} ~ 10^{53} (\epsilon_e/0.1)^{-1.3} (\epsilon_B/0.01)^{-0.09} (t_b/10^4 s) erg, where \epsilon_e and \epsilon_B are the fractions of the total energy at the shock downstream that are carried by electrons and the magnetic field, respectively, and t_b is the end of the plateau. Our finding supports Wolf-Rayet stars as the progenitor stars of some GRBs. It raises intriguing questions about the origin of an intermediate-\Gamma_0 ejecta, which we speculate is connected to the GRB jet emergence from its host star. For the remaining 45% of the sample, the post-plateau decline is too rapid to be explained in the coasting-in-wind model, and energy injection appears to be required.

We present virial models for the global evolution of giant molecular clouds. Focusing on the presence of an accretion flow, and accounting for the amount of mass, momentum, and energy supplied by accretion and star formation feedback, we are able to follow the growth, evolution, and dispersal of individual giant molecular clouds. Our model clouds reproduce the scaling relations observed in both galactic and extragalactic clouds. We find that accretion and star formation contribute contribute roughly equal amounts of turbulent kinetic energy over the lifetime of the cloud. Clouds attain virial equilibrium and grow in such a way as to maintain roughly constant surface densities, with typical surface densities of order 50 - 200 Msun pc^-2, in good agreement with observations of giant molecular clouds in the Milky Way and nearby external galaxies. We find that as clouds grow, their velocity dispersion and radius must also increase, implying that the linewidth-size relation constitutes an age sequence. Lastly, we compare our models to observations of giant molecular clouds and associated young star clusters in the LMC and find good agreement between our model clouds and the observed relationship between H ii regions, young star clusters, and giant molecular clouds.

The Milky Way harbors giant H II regions which may be powered by star complexes more luminous than any Galactic OB association known. Being across the disk of the Galaxy, however, these brightest associations are severely extinguished and confused. We present a search for one such association toward the most luminous H II region in the recent catalog by Murray and Rahman, which, at ~9.7 kpc, has recombination rate of ~7x10^{51} sec$^{-1}$. Prior searches have identified only small scale clustering around the rim of this shell-like region, but the primary association has not previously been identified. We apply a near-infrared color selection and find an overdensity of point sources toward its southern central part. The colors and magnitudes of these excess sources are consistent with O- and early B-type stars at extinctions $0.96 < A_{K} < 1.2$, and they are sufficiently numerous (406 +/- 102 after subtraction of field sources) to ionize the surrounding H II region, making this a candidate for the most luminous OB association in the Galaxy. We reject an alternate theory, in which the apparent excess is caused by localized extinction, as inconsistent with source demographics.

We conduct a survey of numerical simulations to probe the structure and appearance of non-radiative black hole accretion flows like the supermassive black hole at the Galactic centre. We find a generic set of solutions, and make specific predictions for currently feasible rotation measure (RM) observations, which are accessible to current instruments including the EVLA, GMRT and ALMA. The slow time variability of the RM is a key quantitative signature of this accretion flow. The time variability of RM can be used to quantitatively measure the nature of the accretion flow, and to differentiate models. Sensitive measurements of RM can be achieved using RM synthesis or using pulsars. Our energy conserving ideal magneto-hydrodynamical simulations, which achieve high dynamical range by means of a deformed-mesh algorithm, stretch from several Bondi radii to about one thousandth of that radius, and continue for tens of Bondi times. Magnetized flows which lack outward convection possess density slopes around -1, almost independent of physical parameters, and are more consistent with observational constraints than are strongly convective flows We observe no tendency for the flows to become rotationally supported in their centres, or to develop steady outflow. We support these conclusions with formulae which encapsulate our findings in terms of physical and numerical parameters. We discuss the relation of these solutions to other approaches. The main potential uncertainties are the validity of ideal MHD and the absence of a fully relativistic inner boundary condition. The RM variability predictions are testable with current and future telescopes.

We characterize the infall rate onto protostellar systems forming in self-gravitating radiation-hydrodynamic simulations. Using two dimensionless parameters to determine disks' susceptability to gravitational fragmentation, we infer limits on protostellar system multiplicity and the mechanism of binary formation. We show that these parameters give robust predictions even in the case of marginally resolved protostellar disks. We find that protostellar systems with radiation feedback predominately form binaries via turbulent fragmentation, not disk instability, and we predict turbulent fragmentation is the dominant channel for binary formation for low-mass stars. We clearly demonstrate that systems forming in simulations including radiative feedback have fundamentally different parameters than those in purely hydrodynamic simulations.

The Gould Belt Legacy Survey (GBS) on the JCMT has observed a region of 260 square arcminutes in 12CO J=3--2 emission, and a 190 square arcminute subset of this in 13CO and C18O towards the Serpens molecular cloud. We examine the global velocity structure of the non-outflowing gas, and calculate excitation temperatures and opacities. The large scale mass and energetics of the region are evaluated, with special consideration for high velocity gas. We find the cloud to have a mass of 203 solar masses, and to be gravitationally bound, and that the kinetic energy of the outflowing gas is approximately seventy percent of the turbulent kinetic energy of the cloud. We identify compact outflows towards some of the submillimetre Class 0/I sources in the region

We investigate how the removal of interstellar material by stellar feedback limits the efficiency of star formation in molecular clouds and how this determines the shape of the mass function of young star clusters. In particular, we derive relations between the power-law exponents of the mass functions of the clouds and clusters in the limiting regimes in which the feedback is energy-driven and momentum-driven, corresponding to minimum and maximum radiative losses and likely to bracket all realistic cases. We find good agreement between the predicted and observed exponents, especially for momentum-driven feedback, provided the protoclusters have roughly constant mean surface density, as indicated by observations of the star-forming clumps within molecular clouds. We also consider a variety of specific feedback mechanisms, concluding that HII regions inflated by radiation pressure predominate in massive protoclusters, a momentum-limited process when photons can escape after only a few interactions with dust grains. We then present a first estimate of the star formation efficiency in this case, finding that it depends on the masses and sizes of the protoclusters only through their mean surface density, thus ensuring consistency between the observed power-law exponents of the mass functions of the clouds and clusters. The numerical value of this efficiency is also consistent with observations.

We evaluate the role of radiation pressure in the dynamics of HII regions. We first determine under what conditions radiation pressure is significant in comparison to gas pressure and show that, while radiation pressure is generally unimportant for HII regions driven by a handful of massive stars, it is dominant for the larger HII regions produced by the massive star clusters found near the Galactic center and in starburst environments. We then provide a solution for the problem of how HII regions expand when radiation pressure influences their behavior. Finally, we compare radiation-dominated HII regions to other sources of stellar feedback, and argue that HII regions are probably the primary mechanism for regulating the formation of massive star clusters.

We study rapidly accreting, gravitationally unstable disks with a series of global, three dimensional, numerical experiments using the code ORION. In this paper we conduct a numerical parameter study focused on protostellar disks, and show that one can predict disk behavior and the multiplicity of the accreting star system as a function of two dimensionless parameters which compare the disk's accretion rate to its sound speed and orbital period. Although gravitational instabilities become strong, we find that fragmentation into binary or multiple systems occurs only when material falls in several times more rapidly than the canonical isothermal limit. The disk-to-star accretion rate is proportional to the infall rate, and governed by gravitational torques generated by low-m spiral modes. We also confirm the existence of a maximum stable disk mass: disks that exceed ~50% of the total system mass are subject to fragmentation and the subsequent formation of binary companions.

Most analytic work to date on protostellar disks has focused on those in isolation from their environments. However, observations are now beginning to probe the earliest, most embedded phases of star formation, during which disks are rapidly accreting from their parent cores and cannot be modeled in isolation. We present a simple, one-zone model of protostellar accretion disks with high mass infall rates. Our model combines a self-consistent calculation of disk temperatures with an approximate treatment of angular momentum transport via two mechanisms. We use this model to survey the properties of protostellar disks across a wide range of stellar masses and evolutionary times, and make predictions for disks' masses, sizes, spiral structure, and fragmentation that will be directly testable by future large-scale surveys of deeply embedded disks. We define a dimensionless accretion-rotation parameter which, in conjunction with the disk's temperature, controls the disk evolution. We track the dominant mode of angular momentum transport, and demonstrate that for stars with final masses greater than roughly one solar mass, gravitational instabilities are the most important mechanism as most of the mass accumulates. We predict that binary formation through disk fission, fragmentation of the disk into small objects, and spiral arm strength all increase in importance to higher stellar masses.

Recent advances in our understanding of massive star formation have made clear the important role of protostellar disks in mediating accretion. Here we describe a simple, semi-analytic model for young, deeply embedded, massive accretion disks. Our approach enables us to sample a wide parameter space of stellar mass and environmental variables, providing a means to make predictions for a variety of sources that next generation telescopes like ALMA and the EVLA will observe. Moreover we include, at least approximately, multiple mechanisms for angular momentum transport, a comprehensive model for disk heating and cooling, and a realistic estimate for the angular momentum in the gas reservoir. We make predictions for the typical sizes, masses, and temperatures of the disks, and describe the role of gravitational instabilities in determining the binarity fraction and upper mass cut-off.

This paper describes a James Clerk Maxwell Telescope (JCMT) legacy survey that has been awarded roughly 500 hrs of observing time to be carried out from 2007 to 2009. In this survey we will map with SCUBA-2 (Submillimetre Common User Bolometer Array 2) almost all of the well-known low-mass and intermediate-mass star-forming regions within 0.5 kpc that are accessible from the JCMT. Most of these locations are associated with the Gould Belt. From these observations we will produce a flux-limited snapshot of star formation near the Sun, providing a legacy of images, as well as point-source and extended-source catalogues, over almost 700 square degrees of sky. The resulting images will yield the first catalogue of prestellar and protostellar sources selected by submillimetre continuum emission, and should increase the number of known sources by more than an order of magnitude. We will also obtain CO maps with the array receiver HARP (Heterodyne Array Receiver Programme), in three CO isotopologues, of a large typical sample of prestellar and protostellar sources. We will then map the brightest hundred sources with the SCUBA-2 polarimeter (POL-2), producing the first statistically significant set of polarization maps in the submillimetre. The images and source catalogues will be a powerful reference set for astronomers, providing a detailed legacy archive for future telescopes, including ALMA, Herschel and JWST.

Protostellar outflows crisscross the regions of star cluster formation, stirring turbulence and altering the evolution of the forming cluster. We model the stirring of turbulent motions by protostellar outflows, building on an observation that the scaling law of supersonic turbulence implies a momentum cascade analogous to the energy cascade in Kolmogorov turbulence. We then generalize this model to account for a diversity of outflow strengths, and for outflow collimation, both of which enhance turbulence. For a single value of its coupling coefficient the model is consistent with turbulence simulations by Li & Nakamura and, plausibly, with observations of the NGC 1333 cluster-forming region. Outflow-driven turbulence is strong enough to stall collapse in cluster-forming regions for several crossing times, relieving the mismatch between star formation and turbulent decay rates. The predicted line-width-size scaling implies radial density indices between -1 and -2 for regions supported by outflow-driven turbulence, with a tendency for steeper profiles in regions that are more massive or have higher column densities.

We examine whether massive-star accretion disks are likely to fragment due to self-gravity. Rapid accretion and high angular momentum push these disks toward fragmentation, whereas viscous heating and the high protostellar luminosity stabilize them. We find that for a broad range of protostar masses and for reasonable accretion times, massive disks larger than ~150 AU are prone to fragmentation. We develop an analytical estimate for the angular momentum of accreted material, extending the analysis of Matzner and Levin (2005) to account for strongly turbulent initial conditions. In a core-collapse model, we predict that disks are marginally prone to fragmentation around stars of about four to 15 solar masses -- even if we adopt conservative estimates of the disks' radii and tendency to fragment. More massive stars are progressively more likely to fragment, and there is a sharp drop in the stability of disk accretion at the very high accretion rates expected above 110 solar masses. Fragmentation may starve accretion in massive stars, especially above this limit, and is likely to create swarms of small, coplanar companions.

We present semi-analytic dynamical models for giant molecular clouds evolving under the influence of HII regions launched by newborn star clusters. In contrast to previous work, we neither assume that clouds are in virial or energetic equilibrium, nor do we ignore the effects of star formation feedback. The clouds, which we treat as spherical, can expand and contract homologously. Photoionization drives mass ejection; the recoil of cloud material both stirs turbulent motions and leads to an effective confining pressure. The balance between these effects and the decay of turbulent motions through isothermal shocks determines clouds' dynamical and energetic evolution. We find that for realistic values of the rates of turbulent dissipation, photoevaporation, and energy injection by HII regions, the massive clouds where most molecular gas in the Galaxy resides live for a few crossing times, in good agreement with recent observational estimates that large clouds in local group galaxies survive roughly 20-30 Myr. During this time clouds remain close to equilibrium, with virial parameters of 1-3 and column densities near 10^22 H atoms cm^-2, also in agreement with observed cloud properties. Over their lives they convert 5-10% of their mass into stars, after which point most clouds are destroyed when a large HII region unbinds them. In contrast, small clouds like those found in the solar neighborhood only survive ~1 crossing time before being destroyed.

We argue that gravitational instability of typical protostellar disks is not a viable mechanism for the fragmentation into multiple systems -- binary stars, brown dwarf companions, or gas giant planets -- except at periods above roughly 20,000 years. Our conclusion is based on a comparison between prior numerical work on disk self-gravity by Gammie (2001) with our own analytical models for the dynamical and thermal state of protostellar disks. For this purpose we first develop a simple theory for the initial conditions of low-mass star formation, accounting for the effect of turbulence on the characteristic mass, accretion rate, and angular momentum of collaping cores. Second, we examine the criterion for fragmentation to occur during star formation, concentrating on the self-gravitational instabilities of protostellar accretion disks in their main accretion phase. Self-gravitational instabilities are strongly dependent on the thermal state of the disk, and we find that the combination of viscous heating and stellar irradiation quenches fragmentation due to Toomre's local instability. Global instability of the disk may be required to process mass supply, but is unlikely to produce fragments. Our findings help to explain the dearth of substellar companions to stellar type stars: the brown dwarf desert.

We investigate the potential of the upcoming LOBSTER space observatory (due circa 2009) to detect soft X-ray flashes from shock breakout in supernovae, primarily from Type II events. LOBSTER should discover many SN breakout flashes, although the number is sensitive to the uncertain distribution of extragalactic gas columns. X-ray data will constrain the radii of their progenitor stars far more tightly than can be accomplished with optical observations of the SN light curve. We anticipate the appearance of blue supergiant explosions (SN 1987A analogs), which will uncover a population of these underluminous events. We consider also how the mass, explosion energy, and absorbing column can be constrained from X-ray observables alone and with the assistance of optically-determined distances. These conclusions are drawn using known scaling relations to extrapolate, from previous numerical calculations, the LOBSTER response to explosions with a broad range of parameters. We comment on a small population of flashes with 0.2 < z < 0.8 that should exist as transient background events in XMM, Chandra, and ROSAT integrations.

We present a scenario for non-radiative accretion onto the supermassive black hole at the galactic center. Conducting MHD simulations with 1400^3 grid zones that break the axial and reflection symmetries of earlier investigations and extend inward from the Bondi radius, we find a quasi-hydrostatic radial density profile rho propto r^{-0.72} with superadiabatic gradient corresponding to an n~0.72 polytrope. Buoyancy generated by magnetic dissipation is resisted by the same fields so effectively that energy is advected inward: a state of magnetically-frustrated convection. This scenario is consistent with observational constraints on energetics andouter boundary conditions.

I constrain a possible supernova origin for gamma-ray bursts by modeling the dynamical interaction between a relativistic jet and a stellar envelope surrounding it. The delay in observer's time introduced by the jet traversing the envelope should not be long compared to the duration of gamma-ray emission; also, the jet should not be swallowed by a spherical explosion it powers. The only stellar progenitors that comfortably satisfy these constraints, if one assumes that jets move ballistically within their host stars, are compact carbon-oxygen or helium post-Wolf-Rayet stars (type Ic or Ib supernovae); type II supernovae are ruled out. Notably, very massive stars do not appear capable of producing the observed bursts at any redshift unless the stellar envelope is stripped prior to collapse. The presence of a dense stellar wind places an upper limit on the Lorentz factor of the jet in the internal shock model; however, this constraint may be evaded if the wind is swept forward by a photon precursor. Shock breakout and cocoon blowout are considered individually; neither presents a likely source of precursors for cosmological GRBs. These envelope constraints could conceivably be circumvented if jets are laterally pressure-confined while traversing the outer stellar envelope. If so, jets responsible for observed GRBs must either have been launched from a region several hundred kilometers wide, or have mixed with envelope material as they travel. A phase of pressure confinement and mixing would imprint correlations among jets that may explain observed GRB variability-luminosity and lag-luminosity correlations.

We argue that massive stars are the dominant sources of energy for the turbulent motions within giant molecular clouds, and that the primary agent of feedback is the expansion of H II regions within the cloud volume. This conclusion is suggested by the low efficiency of star formation and corroborated by dynamical models of H II regions. We evaluate the turbulent energy input rate in clouds more massive than one third of a million solar masses, for which gravity does not significantly affect the expansion of H II regions. Such clouds achieve a balance between the decay of turbulent energy and its regeneration in H II regions; summed over clouds, the implied ionizing luminosity and star formation rate are roughly consistent with the Galactic total. H II regions also photoevaporate their clouds: we derive cloud destruction times somewhat shorter than those estimated by Williams and McKee. The upper mass limit for molecular clouds in the Milky Way may derive from the fact that larger clouds would destroy themselves in less than one crossing time. The conditions within starburst galaxies do not permit giant molecular clouds to be supported or destroyed by H II regions. This should lead to rapid cloud collapse and the efficient formation of massive star clusters, explaining some aspects of the starburst phenomenon.

Supernova (SN) 1998bw and gamma-ray burst (GRB) 980425 offer the first direct evidence that supernovae are the progenitors of some GRBs. However, this burst was unusually dim, smooth and soft compared to other bursts with known afterglows. Whether it should be considered a prototype for cosmological GRBs depends largely on whether the supernova explosion and burst were asymmetrical or can be modeled as spherical. We address this question by treating the acceleration of the supernova shock in the outermost layers of the stellar envelope, the transition to relativistic flow, and the subsequent expansion (and further acceleration) of the ejecta into the surrounding medium. We find that GRB 980425 could plausibly have been produced by a collision between the relativistic ejecta from SN 1998bw and the star's pre-supernova wind; the model requires no significant asymmetry. This event therefore belongs to a dim subclass of GRBs and is not a prototype for jet-like cosmological GRBs.

The spatiotemporal coincidence of supernova (SN) 1998bw and gamma-ray burst (GRB) 980425 and this supernova's unusual optical and radio properties have prompted many theoretical models that produce GRBs from supernovae. We review the salient features of our simple, spherical model in which an energetic supernova explosion shock accelerates a small fraction of the progenitor's stellar envelope to mildly relativistic velocities. This material carries sufficient energy to produce a weak GRB and a bright radio supernova through an external shock against a dense stellar wind.

We investigate the acceleration of shock waves to relativistic velocities in the outer layers of exploding stars. By concentrating the explosion energy in the outermost ejecta, such trans-relativistic blast waves can serve as the progenitors of gamma-ray bursts (GRBs); in particular, the ``baryon-loading'' problem that plagues many models of GRBs is circumvented. We present physically motivated and numerically validated analytic expressions to describe trans-relativistic blast waves in supernovae. We find that relativistic ejecta are enhanced in more centrally condensed envelopes, e.g., for radiative envelopes, when the luminosity approaches the Eddington limit. We present convenient formulae for estimating the relativistic ejecta from a given progenitor. We apply our analytic and numerical methods to a model of SN 1998bw, finding significantly enhanced relativistic ejecta compared to previous studies. We propose that GRB 980425 is associated with SN 1998bw and may have resulted from an approximately spherical explosion producing ~10^-6 M_sun of mildly relativistic ejecta with mean Lorentz factor ~2, which then interacted with a dense circumstellar wind with mass loss rate ~few x 10^-4 M_sun/yr. A highly asymmetric explosion is not required. An extreme model of ``hypernova'' explosions in massive stars is able to account for the energetics and relativistic ejecta velocities required by many of the observed cosmological GRBs. We present simplified models of explosions resulting from accretion-induced collapse of white dwarfs and phase transitions of neutron stars. While we find increased energies in relativistic ejecta compared to previous studies, these explosions are unlikely to be observed at cosmological distances with current detectors. (abridged)

Using a quantitative model for bipolar outflows driven by hydromagnetic protostellar winds, we calculate the efficiency of star formation assuming that available gas is either converted into stars or ejected in outflows. We estimate the efficiency of a single star formation event in a protostellar core, finding 25%-70% for cores with various possible degrees of flattening. The core mass function and the stellar initial mass function have similar slopes, because the efficiency is not sensitive to its parameters. We then consider the disruption of gas from a dense molecular clump in which a cluster of young stars is being born. In both cases, we present analytical formulae for the efficiencies that compare favorably against observations and, for clusters, against numerical simulations. We predict efficiencies in the range 30%-50% for the regions that form clusters of low-mass stars. In our model, star formation and gas dispersal happen concurrently. We neglect the destructive effects of massive stars: our results are therefore upper limits to the efficiency in regions more massive than about 3000 Msun.

Magnetic stresses collimate protostellar winds into a common distribution of force with angle. Sweeping into the ambient medium, such winds drive bipolar molecular outflows whose properties are insensitive to the distribution of ambient gas and to the details of how the wind is launched, and how its intensity varies over time. Moreover, these properties are in accord with the commonly observed features of outflows. This model is simple enough to permit a quantitative study of the feedback effects from low-mass star formation. It predicts the rate at which star-forming gas is ejected by winds, and hence the efficiency with which stars form. Applied to individual star formation, it relates the stellar initial mass function to the distribution of pre-stellar cores. Applied to cluster formation, it indicates whether the resulting stellar system will remain gravitationally bound. Using the energy injection and mass ejection implied by this model, we investigate the dynamical evolution of a molecular clump as a stellar cluster forms within it. This depends critically on the rate at which turbulence decays: it may involve equilibrium star formation (slow decay), overstable oscillations, or collapse (fast decay).

Rossby waves (r-modes) have been suggested as a means to regulate the spin periods of young or accreting neutron stars, and also to produce observable gravitational wave radiation. R-modes involve primarily transverse, incompressive motions of the star's fluid core. However, neutron stars gain crusts early in their lives: therefore, r-modes also imply shear in the fluid beneath the crust. We examine the criterion for this shear layer to become turbulent, and derive the rate of dissipation in the turbulent regime. Unlike dissipation from a viscous boundary layer, turbulent energy loss is nonlinear in mode energy and can therefore cause the mode to saturate at amplitudes typically much less than unity. This energy loss also reappears as heat below the crust. We study the possibility of crust melting as well as its implications for the spin evolution of low-mass X-ray binaries. Lastly, we identify some universal features of the spin evolution that may have observational consequences.

We demonstrate that magnetically-collimated protostellar winds will sweep ambient material into thin, radiative, momentum-conserving shells whose features reproduce those commonly observed in bipolar molecular outflows. We find the typical position-velocity and mass-velocity relations to occur in outflows in a wide variety of ambient density distributions, regardless of the time histories of their driving winds.

We examine the relation between presupernova stellar structure and the distribution of ejecta in core-collapse supernovae, assuming adiabatic, spherically symmetric flow. We develop a simple yet accurate formula for the blastwave shock velocity, and demonstrate that the entire final density distribution can be approximated with simple models for the final pressure distribution, along with the approximate shock-deposited entropy, in a way that matches the results of simulations. We find that the distribution of density in a star's ejecta depends on whether its outer envelope is radiative or convective, and if convective, on the composition structure of the star; simple approximate forms are presented for red and blue supergiant ejecta. Our models are most accurate for the high-velocity ejecta from the periphery of a star, where the shock dynamics are predictable. We present formulae for the final density distribution of this material, for both radiative and efficiently convective envelopes. These formulae limit to the well-known planar, self-similar solutions for mass shells approaching the stellar surface. But, the assumption of adiabatic flow fails at low optical depth, so this planar limit need not be attained. Formulae are given for the observable properties of the X-ray burst accompanying shock emergence, and their dependence on the parameters of the explosion. Motivated by the relativistic expansion recently inferred by Kulkarni et al. (1998) for the synchrotron shell around SN1998bw, we estimate the criterion for relativistic mass ejection and the rest mass of relativistic ejecta.

We present an analytic calculation of the thermonuclear depletion of the light elements lithium, beryllium, and boron in fully convective, low-mass stars. Under the presumption that the pre--main-sequence star is always fully mixed during contraction, we find that the burning of these rare light elements can be computed analytically, even when the star is degenerate. Using the effective temperature as a free parameter, we constrain the properties of low-mass stars from observational data, independently of the uncertainties associated with modeling their atmospheres and convection. Our analytic solution explains the dependence of the age at a given level of elemental depletion on the stellar effective temperature, nuclear cross sections, and chemical composition. Most importantly, our results allow observers to translate lithium non-detections in young cluster members into a model-independent minimum age for that cluster. Using this procedure, we have found lower limits to the ages of the Pleiades (100 Myr) and Alpha Persei (60 Myr) clusters. Recent experimental work on the low energy resonance in the ^10B(p,\alpha)^7Be reaction has greatly enhanced estimates of the destruction rate of ^10B, making it possible for stars with M>0.1 M_sun to deplete both ^10B and ^11B before reaching the main sequence. Moreover, there is an interesting range of masses, 0.085 M_sun < M < 0.13 M_sun, where boron depletion occurs on the main sequence in less than a Hubble time, providing a potential ``clock'' for dating low-mass stars.

We present an analytic calculation of the thermonuclear depletion of lithium in contracting, fully convective, pre-main sequence stars of mass M < 0.5 M_sun. Previous numerical work relies on still-uncertain physics (atmospheric opacities and convection, in particular) to calculate the effective temperature as a unique function of stellar mass. We assume that the star's effective temperature, T_eff, is fixed during Hayashi contraction and allow its actual value to be a free parameter constrained by observation. Using this approximation, we compute lithium burning analytically and explore the dependence of lithium depletion on T_eff, M, and composition. Our calculations yield the radius, age, and luminosity of a pre-main sequence star as a function of lithium depletion. This allows for more direct comparisons to observations of lithium depleted stars. Our results agree with those numerical calculations that explicitly determine stellar structure during Hayashi contraction. In agreement with Basri, Marcy, and Graham (1996), we show that the absence of lithium in the Pleiades star HHJ 3 implies that it is older than 100 Myr. We also suggest a generalized method for dating galactic clusters younger than 100 Myr (i.e., those with contracting stars of M > 0.08 M_sun) and for constraining the masses of lithium depleted stars.