We have used the Hubble Space Telescope to acquire new broad-band and narrow-band images of the optical line emission and red continuum associated with the HH 7-11 stellar jet in the NGC 1333 star formation region. Combining the new narrow-band images of H$\alpha$, [O~I] $\lambda$6300 and [S II] $\lambda$6716 allows us to measure electron densities and excitations at each point in the outflow with the spatial resolution of HST, while the I-band image traces out the boundary of the cavity evacuated by the outflow. Comparing these images with those taken $\sim$ 20 years ago yields high precision proper motions for all the HH objects in the outflow. HH 11 is a bullet-like clump, and emerges from the exciting source SVS 13A towards the Earth at 24 degrees to line of sight. In contrast, HH 8 and HH 10 consist of two rings of shocked gas that show no bulk proper motions even though the emitting gas is blueshifted. The HH 8 rings are expanding with time. These shocks mark places where ambient material located along the path of the jet redirects the outflow. HH 7 consists of multiple shells, and emits strongly in H$_2$ in what appears to be a terminal bow shock for the outflow, implying that the jet has yet to fully break out of its nascent cloud core. The jet largely fragments into clumps by the time it reaches HH 7. As in the case of HH 110, deflection from ambient material plays a key role in producing observable shock waves in the HH 7-11 outflow.

Young stars exhibit short-term photometric variability caused by mass accretion events from circumstellar disks, the presence of dusty warps within the inner disks, starspots that rotate across the stellar surfaces, and flares. Long-term variability also occurs owing to starspot longevity and cycles, and from changes in stellar angular momenta and activity as the stars age. We propose to observe the Carina star-forming region in different bands with a cadence of 30 minutes every night for one week per year to clarify the nature of both the short-term and long-term variability of the thousands of young stars in this region. By obtaining well-sampled multicolor lightcurves of this dense young cluster, LSST would acquire the first statistically significant data on how these objects vary on both short and long timescales. This information will allow us to relate the observed variability to stellar properties such as mass, age, binarity, and to environmental properties such as location within or exterior to the H II region, and to the presence or absence of a circumstellar disk.

Radiative cooling plays a crucial role in the dynamics of many astrophysical flows, and is particularly important in the dense shocked gas within Herbig-Haro (HH) objects and stellar jets. Simulating cooling processes accurately is necessary to compare numerical simulations with existing and planned observations of HH objects, such as those from the Hubble Space Telescope and the James Webb Space Telescope. In this paper we discuss a new, non-equilibrium cooling scheme we have implemented into the 3-D magnetohydrodynamic (MHD) code AstroBEAR. The new cooling function includes ionization, recombination, and excitation of all the important atomic species that cool below 10000 K. We tested the routine by comparing its predictions with those from the well-tested 1-D Cox-Raymond shock code (Raymond 1979). The results show thatAstroBEAR accurately tracks the ionization fraction, temperature, and other MHD variables for all low-velocity (.90 km/s) magnetized radiative shock waves. The new routine allows us to predict synthetic emission maps in all the bright forbidden and permitted lines observed in stellar jets, including H{\alpha}, [NII], [OI], and [SII]. We present an example as to how these synthetic maps facilitate a direct comparison with narrowband images of HH objects.

We examine new and pre-existing wide-field, continuum-corrected, narrowband images in H$_2$ 1-0 S(1) and Br$\gamma$ of three regions of massive star formation: IC 1396, Cygnus OB2, and Carina. These regions contain a variety of globules, pillars, and sheets, so we can quantify how the spatial profiles of emission lines behave in photodissociation regions (PDRs) that differ in their radiation fields and geometries. We have measured 450 spatial profiles of H$_2$ and Br$\gamma$ along interfaces between HII regions and PDRs. Br$\gamma$ traces photoevaporative flows from the PDRs, and this emission declines more rapidly with distance as the radius of curvature of the interface decreases, in agreement with models. As noted previously, H$_2$ emission peaks deeper into the cloud relative to Br$\gamma$, where the molecular gas absorbs far-UV radiation from nearby O-stars. Although PDRs in IC 1396, Cygnus OB2, and Carina experience orders of magnitude different levels of ionizing flux and have markedly differing geometries, all the PDRs have spatial offsets between Br$\gamma$ and H$_2$ on the order of $10^{17}$cm. There is a weak negative correlation between the offset size and the intensity of ionizing radiation and a positive correlation with the radius of curvature of the cloud. We can reproduce both the size of the offsets and the dependencies of the offsets on these other variables with simple photoevaporative flow models. Both Br$\gamma$ and H$_2$ 1-0 S(1) will undoubtedly be targeted in future JWST observations of PDRs, so this work can serve as a guide to interpreting these images.

Magnetohydrodynamic (MHD) and photoevaporative winds are thought to play an important role in the evolution and dispersal of planet-forming disks. We report the first high-resolution ($\Delta v\sim$6\kms) analysis of [S II] $\lambda$4068, [O I] $\lambda$5577, and [O I] $\lambda$6300 lines from a sample of 48 T Tauri stars. Following Simon et al. (2016), we decompose them into three kinematic components: a high-velocity component (HVC) associated with jets, and a low-velocity narrow (LVC-NC) and broad (LVC-BC) components. We confirm previous findings that many LVCs are blueshifted by more than 1.5 kms$^{-1}$ thus most likely trace a slow disk wind. We further show that the profiles of individual components are similar in the three lines. We find that most LVC-BC and NC line ratios are explained by thermally excited gas with temperatures between 5,000$-$10,000 K and electron densities $\sim10^{7}-10^{8}$ cm$^{-3}$. The HVC ratios are better reproduced by shock models with a pre-shock H number density of $\sim10^{6}-10^{7}$ cm$^{-3}$. Using these physical properties, we estimate $\dot{M}_{\rm wind}/\dot{M}_{\rm acc}$ for the LVC and $\dot{M}_{\rm jet}/\dot{M}_{\rm acc}$ for the HVC. In agreement with previous work, the mass carried out in jets is modest compared to the accretion rate. With the likely assumption that the NC wind height is larger than the BC, the LVC-BC $\dot{M}_{\rm wind}/\dot{M}_{\rm acc}$ is found to be higher than the LVC-NC. These results suggest that most of the mass loss occurs close to the central star, within a few au, through an MHD driven wind. Depending on the wind height, MHD winds might play a major role in the evolution of the disk mass.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

We present the science cases and technological discussions that came from the workshop entitled "Finding the UV-Visible Path Forward" held at NASA GSFC June 25-26, 2015. The material presented outlines the compelling science that can be enabled by a next generation space-based observatory dedicated for UV-visible science, the technologies that are available to include in that observatory design, and the range of possible alternative launch approaches that could also enable some of the science. The recommendations to the Cosmic Origins Program Analysis Group from the workshop attendees on possible future development directions are outlined.

The ~2 Myr old classical T Tauri star CI Tau shows periodic variability in its radial velocity (RV) variations measured at infrared (IR) and optical wavelengths. We find that these observations are consistent with a massive planet in a ~9-day period orbit. These results are based on 71 IR RV measurements of this system obtained over 5 years, and on 26 optical RV measurements obtained over 9 years. CI Tau was also observed photometrically in the optical on 34 nights over ~one month in 2012. The optical RV data alone are inadequate to identify an orbital period, likely the result of star spot and activity induced noise for this relatively small dataset. The infrared RV measurements reveal significant periodicity at ~9 days. In addition, the full set of optical and IR RV measurements taken together phase coherently and with equal amplitudes to the ~9 day period. Periodic radial velocity signals can in principle be produced by cool spots, hot spots, and reflection of the stellar spectrum off the inner disk, in addition to resulting from a planetary companion. We have considered each of these and find the planet hypothesis most consistent with the data. The radial velocity amplitude yields an Msin(i) of ~8.1 M_Jup; in conjunction with a 1.3 mm continuum emission measurement of the circumstellar disk inclination from the literature, we find a planet mass of ~11.3 M_Jup, assuming alignment of the planetary orbit with the disk.

Supersonic outflows from objects as varied as stellar jets, massive stars and novae often exhibit multiple shock waves that overlap one another. When the intersection angle between two shock waves exceeds a critical value, the system reconfigures its geometry to create a normal shock known as a Mach stem where the shocks meet. Mach stems are important for interpreting emission-line images of shocked gas because a normal shock produces higher postshock temperatures and therefore a higher-excitation spectrum than an oblique one does. In this paper we summarize the results of a series of numerical simulations and laboratory experiments designed to quantify how Mach stems behave in supersonic plasmas that are the norm in astrophysical flows. The experiments test analytical predictions for critical angles where Mach stems should form, and quantify how Mach stems grow and decay as intersection angles between the incident shock and a surface change. While small Mach stems are destroyed by surface irregularities and subcritical angles, larger ones persist in these situations, and can regrow if the intersection angle changes to become more favorable. The experimental and numerical results show that although Mach stems occur only over a limited range of intersection angles and size scales, within these ranges they are relatively robust, and hence are a viable explanation for variable bright knots observed in HST images at the intersections of some bow shocks in stellar jets.

High resolution observations of Young Stellar Object (YSO) jets show them to be composed of many small-scale knots or clumps. In this paper we report results of 3-D numerical simulations designed to study how such clumps interact and create morphologies and kinematic patterns seen in emission line observations. Our simulations focus on clump scale dynamics by imposing velocity differences between spherical, over-dense regions which then lead to the formation of bow shocks as faster clumps overtake slower material. We show that much of the spatial structure apparent in emission line images of jets arises from the dynamics and interactions of these bow shocks. Our simulations show a variety of time-dependent features, including bright knots associated with Mach stems where the shocks intersect, a "frothy" emission structure that arises from the presence of the Non-linear Thin Shell Instability (NTSI) along the surfaces of the bow shocks, and the merging and fragmentation of clumps. Our simulations use a new non-equilibrium cooling method to produce synthetic emission maps in H$\alpha$ and [S II]. These are directly compared with multi-epoch Hubble Space Telescope (HST) observations of Herbig-Haro (HH) jets. We find excellent agreement between features seen in the simulations and the observations in terms of both proper motion and morphologies. Thus we conclude that YSO jets may be dominated by heterogeneous structures and that interactions between these structures and the shocks they produce can account for many details of YSO jet evolution.

We show that it is possible to measure Alfv\'enic Mach numbers, defined as the shock velocity in the flow divided by the Alfv\'en velocity, for low-velocity (V$_{shock}$ $\lesssim$ 100 km$\,$s$^{-1}$) radiative shocks. The method combines observations of bright forbidden lines with a measure of the size of the cooling zone, the latter typically obtained from spatial separation between the Balmer emission lines and the forbidden lines. Because magnetic fields become compressed as gas in the postshock region cools, even relatively weak preshock magnetic fields can be detected with this method. We derive analytical formulae that explain how the spatial separations relate to emission-line ratios, and compute a large grid of radiatively-cooled shock models to develop diagnostic diagrams that can be used to derive Alfv\'enic Mach numbers in flows. Applying the method to existing data for a bright knot in the HH 111 jet, we obtain a relatively low Alfv\'enic Mach number of $\sim$ 2, indicative of a magnetized jet that has super-magnetosonic velocity perturbations within it.

We present wide-field, deep narrowband H$_2$, Br$\gamma$, H$\alpha$, [S II], [O III], and broadband I and K-band images of the Carina star formation region. The new images provide a large-scale overview of all the H$_2$ and Br$\gamma$ emission present in over a square degree centered on this signature star forming complex. By comparing these images with archival HST and Spitzer images we observe how intense UV radiation from O and B stars affects star formation in molecular clouds. We use the images to locate new candidate outflows and identify the principal shock waves and irradiated interfaces within dozens of distinct areas of star-forming activity. Shocked molecular gas in jets traces the parts of the flow that are most shielded from the intense UV radiation. Combining the H$_2$ and optical images gives a more complete view of the jets, which are sometimes only visible in H$_2$. The Carina region hosts several compact young clusters, and the gas within these clusters is affected by radiation from both the cluster stars and the massive stars nearby. The Carina Nebula is ideal for studying the physics of young H II regions and PDR's, as it contains multiple examples of walls and irradiated pillars at various stages of development. Some of the pillars have detached from their host molecular clouds to form proplyds. Fluorescent H$_2$ outlines the interfaces between the ionized and molecular gas, and after removing continuum, we detect spatial offsets between the Br$\gamma$ and H$_2$ emission along the irradiated interfaces. These spatial offsets can be used to test current models of PDRs once synthetic maps of these lines become available.

Hubble Space Telescope observations show bright knots of H$\alpha$ emission within outflowing young stellar jets. Velocity variations in the flow create secondary bow shocks that may intersect and lead to enhanced emission. When the bow shocks intersect at or above a certain critical angle, a planar shock called a Mach stem is formed. These shocks could produce brighter H$\alpha$ emission since the incoming flow to the Mach stem is parallel to the shock normal. In this paper we report first results of a study using 2-D numerical simulations designed to explore Mach stem formation at the intersection of bow shocks formed by hypersonic "bullets" or "clumps". Our 2-D simulations show how the bow shock shapes and intersection angles change as the adiabatic index $\gamma$ changes. We show that the formation or lack of a Mach stem in our simulations is consistent with the steady-state Mach stem formation theory. Our ultimate goal, which is part of an ongoing research effort, is to characterize the physical and observational consequences of bow shock intersections including the formation of Mach stems.

In this paper we explore the dynamics of radiative axisymmetric MHD jets at high resolution using AMR methods. The goal of the study is to determine both the dynamics and emission properties of such jets. To that end we have implemented microphysics enabling us to produce synthetic maps of H$\alpha$ and [S II]. The jets are pulsed either sinusoidally or randomly via a time-dependent ejection velocity which leads to a complicated structure of internal shocks and rarefactions as has been seen in previous simulations. The high resolution of our simulations allows us to explore in great detail the effect of pinch forces (due to the jet's toroidal magnetic field) within the "working surfaces" where pulses interact. We map the strong H$\alpha$ emission marking shock fronts and the strong [S II] emission inside cooling regions behind shocks as observed with high-resolution images of jets. We find that pinch forces in the stronger field cases produce additional emission regions along the axis as compared with purely hydrodynamic runs. These simulations are a first step to understanding the full 3-D emission properties of radiative MHD jets.

We consider the lifetimes of spots on the Sun and other stars from the standpoint of magnetic diffusion. While normal magnetic diffusivity predicts lifetimes of sunspots that are too large by at least two orders of magnitude, turbulent magnetic diffusivity accounts for both the functional form of the solar empirical spot-lifetime relation and for the observed sunspot lifetimes, provided that the relevant diffusion length is the supergranule size. Applying this relation to other stars, the value of turbulent diffusivity depends almost entirely on supergranule size, with very weak dependence on other variables such as magnetic field strength and density. Overall, the best observational data for other stars is consistent with the extension of the solar relation provided that stellar supergranule sizes for some stars are significantly larger than they are on the Sun.

In this review we focus on the role jets and outflows play in the star and planet formation process. Our essential question can be posed as follows: are jets/outflows merely an epiphenomenon associated with star formation or do they play an important role in mediating the physics of assembling stars both individually and globally? We address this question by reviewing the current state of observations and their key points of contact with theory. Our review of jet/outflow phenomena is organized into three length-scale domains: Source and Disk Scales ($0.1-10^2$ au) where the connection with protostellar and disk evolution theories is paramount; Envelope Scales ($10^2-10^5$ au) where the chemistry and propagation shed further light on the jet launching process, its variability and its impact on the infalling envelope; Parent Cloud Scales ($10^5-10^6$ au) where global momentum injection into cluster/cloud environments become relevant. Issues of feedback are of particular importance on the smallest scales where planet formation regions in a disk may be impacted by the presence of disk winds, irradiation by jet shocks or shielding by the winds. Feedback on envelope scales may determine the final stellar mass (core-to-star efficiency) and envelope dissipation. Feedback also plays an important role on the larger scales with outflows contributing to turbulent support within clusters including alteration of cluster star formation efficiencies (feedback on larger scales currently appears unlikely). A particularly novel dimension of our review is that we consider results on jet dynamics from the emerging field of High Energy Density Laboratory Astrophysics (HEDLA). HEDLA is now providing direct insights into the 3-D dynamics of fully magnetized, hypersonic, radiative outflows.

Regions of massive star formation offer some of the best and most easily-observed examples of radiation hydrodynamics. Boundaries where fully-ionized H II regions transition to neutral/molecular photodissociation regions (PDRs) are of particular interest because marked temperature and density contrasts across the boundaries lead to evaporative flows and fluid dynamical instabilities that can evolve into spectacular pillar-like structures. When detached from their parent clouds, pillars become ionized globules that often harbor one or more young stars. H2 molecules at the interface between a PDR and an H II region absorb ultraviolet light from massive stars, and the resulting fluoresced infrared emission lines are an ideal way to trace this boundary independent of obscuring dust. This paper presents H2 images of four regions of massive star formation that illustrate different types of PDR boundaries. The Ara OB1 star formation region contains a striking long wall that has several wavy structures which are present in H2, but the emission is not particularly bright because the ambient UV fluxes are relatively low. In contrast, the Carina star formation region shows strong H2 fluorescence both along curved walls and at the edges of spectacular pillars that in some cases have become detached from their parent clouds. The less-spectacular but more well-known Eagle Nebula has two regions that have strong fluorescence in addition to its pillars. While somewhat older than the other regions, Cyg OB2 has the highest number of massive stars of the regions surveyed and contains many isolated, fluoresced globules that have head-tail morphologies which point towards the sources of ionizing radiation. These images provide a collection of potential astrophysical analogs that may relate to ablated interfaces observed in laser experiments of radiation hydrodynamics.

We present results from an ongoing multiwavelength radial velocity (RV) survey of the Taurus-Auriga star forming region as part of our effort to identify pre--main sequence giant planet hosts. These 1-3 Myr old T Tauri stars present significant challenges to traditional RV surveys. The presence of strong magnetic fields gives rise to large, cool star spots. These spots introduce significant RV jitter which can mimic the velocity modulation from a planet-mass companion. To distinguish between spot-induced and planet-induced RV modulation, we conduct observations at ~6700 Angstroms and ~2.3 microns and measure the wavelength dependence (if any) in the RV amplitude. CSHELL observations of the known exoplanet host Gl 86 demonstrate our ability to detect not only hot Jupiters in the near infrared but also secular trends from more distant companions. Observations of nine very young stars reveal a typical reduction in RV amplitude at the longer wavelengths by a factor of ~2-3. While we can not confirm the presence of planets in this sample, three targets show different periodicities in the two wavelength regions. This suggests different physical mechanisms underlying the optical and K band variability.

Modern theoretical models of astrophysical jets combine accretion, rotation, and magnetic fields to launch and collimate supersonic flows from a central source. Near the source, magnetic field strengths must be large enough to collimate the jet requiring that the Poynting flux exceeds the kinetic-energy flux. The extent to which the Poynting flux dominates kinetic energy flux at large distances from the engine distinguishes two classes of models. In magneto-centrifugal launch (MCL) models, magnetic fields dominate only at scales $\lesssim 100$ engine radii, after which the jets become hydrodynamically dominated (HD). By contrast, in Poynting flux dominated (PFD) magnetic tower models, the field dominates even out to much larger scales. To compare the large distance propagation differences of these two paradigms, we perform 3-D ideal MHD AMR simulations of both HD and PFD stellar jets formed via the same energy flux. We also compare how thermal energy losses and rotation of the jet base affects the stability in these jets. For the conditions described, we show that PFD and HD exhibit observationally distinguishable features: PFD jets are lighter, slower, and less stable than HD jets. Unlike HD jets, PFD jets develop current-driven instabilities that are exacerbated as cooling and rotation increase, resulting in jets that are clumpier than those in the HD limit. Our PFD jet simulations also resemble the magnetic towers that have been recently created in laboratory astrophysical jet experiments.

Magnetic Towers represent one of two fundamental forms of MHD outflows. Driven by magnetic pressure gradients, these flows have been less well studied than magneto-centrifugally launched jets even though magnetic towers may well be as common. Here we present new results exploring the behavior and evolution of magnetic tower outflows and demonstrate their connection with pulsed power experimental studies and purely hydrodynamic jets which might represent the asymptotic propagation regimes of magneto-centrifugally launched jets. High-resolution AMR MHD simulations (using the AstroBEAR code) provide insights into the underlying physics of magnetic towers and help us constrain models of their propagation. Our simulations have been designed to explore the effects of thermal energy losses and rotation on both tower flows and their hydro counterparts. We find these parameters have significant effects on the stability of magnetic towers, but mild effects on the stability of hydro jets. Current-driven perturbations in the Poynting Flux Dominated (PDF) towers are shown to be amplified in both the cooling and rotating cases. Our studies of the long term evolution of the towers show that the formation of weakly magnetized central jets within the tower are broken up by these instabilities becoming a series of collimated clumps which magnetization properties vary over time. In addition to discussing these results in light of laboratory experiments, we address their relevance to astrophysical observations of young star jets and outflow from highly evolved solar type stars.

We report optical (6150 Ang) and K-band (2.3 micron) radial velocities obtained over two years for the pre-main sequence weak-lined T Tauri star Hubble I 4. We detect periodic and near-sinusoidal radial velocity variations at both wavelengths, with a semi-amplitude of 1395\pm94 m/s in the optical and 365\pm80 m/s in the infrared. The lower velocity amplitude at the longer wavelength, combined with bisector analysis and spot modeling, indicates that there are large, cool spots on the stellar surface that are causing the radial velocity modulation. The radial velocities maintain phase coherence over hundreds of days suggesting that the starspots are long-lived. This is one of the first active stars where the spot-induced velocity modulation has been resolved in the infrared.

We present new, third-epoch HST H-alpha and [S II] images of three HH jets (HH 1 and 2, HH 34, and HH 47) and compare these images with those from the previous epochs. The high-spatial resolution, coupled with a time-series whose cadence is of order both the hydrodynamical and radiative cooling timescales of the flow allows us to follow the hydrodynamical/magnetohydrodynamical evolution of an astrophysical plasma system in which ionization and radiative cooling play significant roles. Cooling zones behind the shocks are resolved, so it is possible to identify which way material flows through a given shock wave. The images show that heterogeneity is paramount in these jets, with clumps dominating the morphologies of both bow shocks and their Mach disks. This clumpiness exists on scales smaller than the jet widths and determines the behavior of many of the features in the jets. Evidence also exists for considerable shear as jets interact with their surrounding molecular clouds, and in several cases we observe shock waves as they form and fade where material emerges from the source and as it proceeds along the beam of the jet. Fine-structure within two extended bow shocks may result from Mach stems that form at the intersection points of oblique shocks within these clumpy objects. Taken together, these observations represent the most significant foray thus far into the time domain for stellar jets, and comprise one of the richest data sets in existence for comparing the behavior of a complex astrophysical plasma flows with numerical simulations and laboratory experiments.

Collimated supersonic flows in laboratory experiments behave in a similar manner to astrophysical jets provided that radiation, viscosity, and thermal conductivity are unimportant in the laboratory jets, and that the experimental and astrophysical jets share similar dimensionless parameters such as the Mach number and the ratio of the density between the jet and the ambient medium. Laboratory jets can be studied for a variety of initial conditions, arbitrary viewing angles, and different times, attributes especially helpful for interpreting astronomical images where the viewing angle and initial conditions are fixed and the time domain is limited. Experiments are also a powerful way to test numerical fluid codes in a parameter range where the codes must perform well. In this paper we combine images from a series of laboratory experiments of deflected supersonic jets with numerical simulations and new spectral observations of an astrophysical example, the young stellar jet HH 110. The experiments provide key insights into how deflected jets evolve in 3-D, particularly within working surfaces where multiple subsonic shells and filaments form, and along the interface where shocked jet material penetrates into and destroys the obstacle along its path. The experiments also underscore the importance of the viewing angle in determining what an observer will see. The simulations match the experiments so well that we can use the simulated velocity maps to compare the dynamics in the experiment with those implied by the astronomical spectra. The experiments support a model where the observed shock structures in HH 110 form as a result of a pulsed driving source rather than from weak shocks that may arise in the supersonic shear layer between the Mach disk and bow shock of the jet's working surface.

We present new adaptive optics echelle spectra of the near-infrared [Fe II] lines in the redshifted and blueshifted jets from the T Tauri star RW Aur. The spectra have an unprecedented combination of high spatial and spectral resolution that makes it possible to trace the dynamics of the flow to a projected distance of only 10 AU from the source. As noted by previous studies, the redshifted flow is much slower than its fainter blueshifted counterpart. Our observations clearly show that both the radial velocities and the emission line widths are larger closer to the source on both sides of the jet. The line widths are 20% - 30% of the jet velocity on both sides of the flow, significantly larger than would be produced by a divergent constant velocity flow. The observed line widths could arise from a layered velocity structure in the jet or from magnetic waves. A bright knot in the redshifted jet has no concomitant increase in line width, implying that it is not heated by a bow shock. Alternate heating mechanisms include planar shocks, ambipolar diffusion and magnetic reconnection.

In searches for low-mass companions to late-type stars, correlation between radial velocity variations and line bisector slope changes indicates contamination by large starspots. Two young stars demonstrate that this test is not sufficient to rule out starspots as a cause of radial velocity variations. As part of our survey for substellar companions to T Tauri stars, we identified the ~2 Myr old planet host candidates DN Tau and V836 Tau. In both cases, visible light radial velocity modulation appears periodic and is uncorrelated with line bisector span variations, suggesting close companions of several M_Jup in these systems. However, high-resolution, infrared spectroscopy shows that starspots cause the radial velocity variations. We also report unambiguous results for V827 Tau, identified as a spotted star on the basis of both visible light and infrared spectroscopy. Our results suggest that infrared follow up observations are critical for determining the source of radial velocity modulation in young, spotted stars.

We describe a new radial velocity survey of T Tauri stars and present the first results. Our search is motivated by an interest in detecting massive young planets, as well as investigating the origin of the brown dwarf desert. As part of this survey, we discovered large-amplitude, periodic, radial velocity variations in the spectrum of the weak line T Tauri star LkCa 19. Using line bisector analysis and a new simulation of the effect of star spots on the photometric and radial velocity variability of T Tauri stars, we show that our measured radial velocities for LkCa19 are fully consistent with variations caused by the presence of large star spots on this rapidly rotating young star. These results illustrate the level of activity-induced radial velocity noise associated with at least some very young stars. This activity-induced noise will set lower limits on the mass of a companion detectable around LkCa 19, and similarly active young stars.

Although several lines of evidence suggest that jets from young stars are driven magnetically from accretion disks, existing observations of field strengths in the bow shocks of these flows imply that magnetic fields play only a minor role in the dynamics at these locations. To investigate this apparent discrepancy we performed numerical simulations of expanding magnetized jets with stochastically variable input velocities with the AstroBEAR MHD code. Because the magnetic field B is proportional to the density n within compression and rarefaction regions, the magnetic signal speed drops in rarefactions and increases in the compressed areas of velocity-variable flows. In contrast, B ~ n^0.5 for a steady-state conical flow with a toroidal field, so the Alfven speed in that case is constant along the entire jet. The simulations show that the combined effects of shocks, rarefactions, and divergent flow cause magnetic fields to scale with density as an intermediate power 1 > p > 0.5. Because p > 0.5, the Alfven speed in rarefactions decreases on average as the jet propagates away from the star. This behavior is extremely important to the flow dynamics because it means that a typical Alfven velocity in the jet close to the star is significantly larger than it is in the rarefactions ahead of bow shocks at larger distances, the one place where the field is a measurable quantity. We find that the observed values of weak fields at large distances are consistent with strong fields required to drive the observed mass loss close to the star. For a typical stellar jet the crossover point inside which velocity perturbations of 30 - 40 km/s no longer produce shocks is ~ 300 AU from the source.

We present Space Telescope Imaging Spectrograph (STIS) spectral images of the HH~30 stellar jet taken through a wide slit over two epochs. The jet is unresolved spectrally, so the observations produce emission-line images for each line in the spectrum. This rich dataset shows how physical conditions in the jet vary with distance and time, produces precise proper motions of knots within the jet, resolves the jet width close to the star, and gives a spectrum of the reflected light from the disk over a large wavelength range at several positions. We introduce a new method for analyzing a set of line ratios based on minimizing a quadratic form between models and data. The method generates images of the density, temperature and ionization fraction computed using all the possible line ratios appropriately weighted. In HH 30, the density declines with distance from the source in a manner consistent with an expanding flow, and is larger by a factor of two along the axis of the jet than it is at the periphery. Ionization in the jet ranges from ~ 5% to 40%, and high ionization/excitation knots form at about 100 AU from the star and propagate outward with the flow. These high-excitation knots are not accompanied by corresponding increases in the density, so if formed by velocity variations the knots must have a strong internal magnetic pressure to smooth out density increases while lengthening recombination times.

We present moderate and high-dispersion 1-2.5 micron spectra of the 10 arcsec nebula around P Cygni, dominated by bright emission lines of [Fe II]. Observed [Fe II] line ratios disagree with theoretical transition rates in the literature, so we use the spectrum of P Cygni's nebula to constrain the atomic data for low-lying levels of [Fe II]. Of particular interest is the ratio [Fe II] 12567/16435, often used as a reddening indicator, for which we empirically derive an intrinsic value of 1.49, which is 10--40 per cent higher than previous estimates. High-dispersion spectra of [Fe II] 16435 constrain the geometry, detailed structure, and kinematics of P Cygni's nebula, which is the major product of P Cygni's outburst in 1600 A.D. We use the [N II]/[N I] line ratio to conclude that the nebula is mostly ionized with a total mass of ~0.1 Msun; more than the mass lost by the stellar wind since the eruption. For this mass, we would expect a larger infrared excess than observed. We propose that the dust which obscured the star after the outburst has since been largely destroyed, releasing Fe into the gas phase to produce the bright [Fe II] emission. The kinetic energy of this shell is ~2e46 ergs, far less than the kinetic energy released during the giant eruption of eta Car in the 1840s. The mass and kinetic energy in the nebulae of eta Car and P Cygni give insight to the range of parameters expected for extragalactic eta Car-like eruptions.

We present a proper motion study of the shock waves within the classic stellar jet HH 47 based on Hubble Space Telescope H-alpha and [S II] images of the region taken over two epochs. Individual knots within the jet and in the bow shock/Mach disk working surface of HH 47A move significantly in the five years that separate the images, and the excellent spatial resolution of HST makes it possible to measure the proper motions with enough precision to easily observe differential motions throughout the flow. The bright portion of the jet emerges at 37.5 +/- 2.5 degrees from the plane of the sky with an average velocity of 300 km/s. Dynamical ages of the shock waves in the jet range from a few decades for knots recently ejected by the source to ~ 1300 years for the faint extended bow shock HH 47D. The jet curves, but motions of knots in the jet are directed radially away from the exciting source, and velocity variability in the flow drives the shock waves that heat the jet internally. The jet orientation angle varies with time by about 15 degrees, and currently points to the northwestern portion of a cavity outlined by a reflection nebula, where a quasi-stationary shock deflects the jet. The major working surface HH 47A is more complex than a simple bow shock/Mach disk, and contains numerous clumps that move relative to one another with velocities of ~ +/- 40 km/s. Small clumps or instabilities affect the Mach disk, and dense clumps may move all the way through the working surface to cause the bumpy morphology seen at the bow shock. A localized area between the bow shock and Mach disk varies significantly between the two sets of images.

We present echelle spectroscopy of the close pre-main-sequence binary star systems DQ Tau and UZ Tau-E. Over a 16 day time interval we acquired 14 nights of spectra for DQ Tau and 12 nights of spectra for UZ Tau-E. This represents the entire phase of DQ Tau, and 63 percent of the phase of UZ Tau-E. As expected, photospheric lines such as Li I 6707 clearly split into two components as the primary and secondary orbit one another, as did the permitted line He I 5876. Unlike the photospheric features, the forbidden lines of [O I] 6300 and [O I] 5577, retain the same shape throughout the orbit. Therefore these lines must originate outside of the immediate vicinity of the two stars and any circumstellar disks that participate in the orbital motion of the stars.

We report the results of a spectroscopic survey of 20 close T Tauri binaries in the Taurus-Auriga dark cloud where the separations between primaries and their secondaries are less than the typical size of a circumstellar disk around a young star. Analysis of low-resolution and medium-resolution STIS spectra yields the stellar luminosities, reddenings, ages, masses, mass accretion rates, IR excesses, and emission line luminosities for each star in each pair. We examine the ability of IR color excesses, H-alpha equivalent widths, [O I] emission, and veiling to distinguish between weak emission and classical T Tauri stars. Four pairs have one cTTs and one wTTs; the cTTs is the primary in three of these systems. This frequency of mixed pairs among the close T Tauri binaries is similar to the frequency of mixed pairs in wider young binaries. Extinctions within pairs are usually similar; however, the secondary is more heavily reddened than the primary in some systems, where it may be viewed through the primary's disk. Mass accretion rates of primaries and secondaries are strongly correlated, and H-alpha luminosities, IR excesses, and ages also correlate within pairs. Primaries tend to have somewhat larger accretion rates than their secondaries do, and are typically slightly older than their secondaries according to three different sets of modern pre-main-sequence evolutionary tracks. Age differences for XZ Tau and FS Tau, systems embedded in reflection nebulae, are striking; the secondary in each pair is less massive but more luminous than the primary. The stellar masses of the UY Aur and GG Tau binaries measured from their rotating molecular disks are about 30% larger than the masses inferred from the spectra and evolutionary tracks.

We present an observational and theoretical study of the optical emission from nonradiative shocks in three supernova remnants: the Cygnus Loop, RCW 86 and Tycho. The spectra of these shocks are dominated by collisionally excited hydrogen Balmer lines which have both a broad component caused by proton-neutral charge exchange and a narrow component caused by excitation of cold neutrals entering the shock. In each remnant we have obtained the broad to narrow flux ratios of the Halpha and Hbeta lines and measured the Halpha broad component width. A new numerical shock code computes the broad and narrow Balmer line emission from nonradiative shocks in partially neutral gas. The Balmer line fluxes are sensitive to Lyman line trapping and the degree of electron-proton temperature equilibration. The code calculates the density, temperature and size of the postshock ionization layer and uses a Monte Carlo simulation to compute narrow Balmer line enhancement from Lyman line trapping. The initial fraction of the shock energy allocated to the electrons and protons (the equilibration) is a free parameter. Our models show that variations in electron-proton temperature equilibration and Lyman line trapping can reproduce the observed range of broad to narrow ratios. The results give 80%-100% equilibration in nonradiative portions of the NE Cygnus Loop (v_shock ~ 300 km/s), 40%-50% equilibration in nonradiative portions of RCW 86 (v_shock ~ 600 km/s) and <~ 20% equilibration in Tycho (v_shock ~ 2000 km/s). Our results suggst an inverse correlation between magnetosonic Mach number and equilibration in the observed remnants.

We have applied a new restoration technique to archival [O~I], H$\alpha$, and continuum HST images of DG~Tau. The restored [O~I] and H$\alpha$ images show that DG~Tau has a jet with a projected length of 25~AU and width $\leq$10~AU, and is already collimated at a projected distance of $\sim$~40~AU (0\farcs25) from the star. Such a narrow width and short collimation distance for a stellar jet places important constraints on theoretical models of jet formation.