Stellar bow shocks result from relative motions between stars and their environment. The interaction of the stellar wind and radiation with gas and dust in the interstellar medium produces curved arcs of emission at optical, infrared, and radio wavelengths. We recently proposed a new two-dimensional classification scheme for the shape of such bow shocks, which we here apply to three very different observational datasets: mid-infrared arcs around hot OB stars; far-infrared arcs around luminous cool stars; and H alpha emission-line arcs around proplyds and other young stars in the Orion Nebula. For OB stars, the average shape is consistent with simple thin-shell models for the interaction of a spherical wind with a parallel stream, but the diversity of observed shapes is many times larger than such models predict. We propose that this may be caused by time-dependent oscillations in the bow shocks, due to either instabilities or wind variability. Cool star bow shocks have markedly more closed wings than hot star bow shocks, which may be due to the dust emission arising in the shocked stellar wind instead of the shocked interstellar medium. The Orion Nebula arcs, on the other hand, have both significantly more open wings and significantly flatter apexes than the hot star bow shocks. We test several possible explanations for this difference (divergent ambient stream, low Mach number, observational biases, and influence of collimated jets), but the evidence for each is inconclusive.

Dust waves and bow waves result from the action of a star's radiation pressure on a stream of dusty plasma that flows past it. They are an alternative mechanism to hydrodynamic bow shocks for explaining the curved arcs of infrared emission seen around some stars. When gas and grains are perfectly coupled, for a broad class of stellar parameters, wind-supported bow shocks predominate when the ambient density is below 100 per cubic cm. At higher densities radiation-supported bows can form, tending to be optically thin bow waves around B stars, or optically thick bow shocks around early O stars. For OB stars with particularly weak stellar winds, radiation-supported bows become more prevalent.

Astrophysical bow shocks are a common result of the interaction between two supersonic plasma flows, such as winds or jets from stars or active galaxies, or streams due to the relative motion between a star and the interstellar medium. For cylindrically symmetric bow shocks, we develop a general theory for the effects of inclination angle on the apparent shape. We propose a new two-dimensional classification scheme for bow shapes, which is based on dimensionless geometric ratios that can be estimated from observational images. The two ratios are related to the flatness of the bow's apex, which we term planitude and the openness of its wings, which we term alatude. We calculate the expected distribution in the planitude-alatude plane for a variety of simple geometrical and physical models: quadrics of revolution, wilkinoids, cantoids, and ancantoids. We further test our methods against numerical magnetohydrodynamical simulations of stellar bow shocks and find that the apparent planitude and alatude measured from infrared dust continuum maps serve as accurate diagnostics of the shape of the contact discontinuity, which can be used to discriminate between different physical models. We present an algorithm that can determine the planitude and alatude from observed bow shock emission maps with a precision of 10 to 20%.

Stellar bow shocks, bow waves, and dust waves all result from the action of a star's wind and radiation pressure on a stream of dusty plasma that flows past it. The dust in these bows emits prominently at mid-infrared wavelengths in the range 8 to 60 micron. We propose a novel diagnostic method, the tau-eta diagram, for analyzing these bows, which is based on comparing the fractions of stellar radiative energy and stellar radiative momentum that is trapped by the bow shell. This diagram allows the discrimination of wind-supported bow shocks, radiation-supported bow waves, and dust waves in which grains decouple from the gas. For the wind-supported bow shocks, it allows the stellar wind mass-loss rate to be determined. We critically compare our method with a previous method that has been proposed for determining wind mass-loss rates from bow shock observations. This comparison points to ways in which both methods can be improved and suggests a downward revision by a factor of two with respect to previously reported mass-loss rates. From a sample of 23 mid-infrared bow-shaped sources, we identify at least 4 strong candidates for radiation-supported bow waves, which need to be confirmed by more detailed studies, but no strong candidates for dust waves.

Dust waves are a result of gas-grain decoupling in a stream of dusty plasma that flows past a luminous star. The radiation field is sufficiently strong to overcome the collisional coupling between grains and gas at a "rip-point", where the ratio of radiation pressure to gas pressure exceeds a critical value of roughly 1000. When the rip point occurs outside the hydrodynamic bow shock, a separate dust wave may form, decoupled from the gas shell, which can either be drag-confined or inertia-confined, depending on the stream density and relative velocity. In the drag-confined case, there is a minimum stream velocity of roughly 60 km/s that allows a steady-state stagnant drift solution for the dust wave apex. For lower relative velocities, the dust dynamics close to the axis exhibit a limit cycle behavior (rip and snap back) between two different radii. Strong coupling of charged grains to the plasma's magnetic field can modify these effects, but for a quasi-parallel field orientation the results are qualitatively similar to the non-magnetic case. For a quasi-perpendicular field, on the other hand, the formation of a decoupled dust wave is strongly suppressed.

We present the results of long-slit echelle spectroscopy and deep narrow-band imaging of the Owl Nebula (NGC 3587), obtained at the \textit{Observatorio Astron\'omico Nacional, San Pedro M\'artir}. These data allow us to construct an iso-velocity data cube and develop a 3-D morpho-kinematic model. We find that, instead of the previously assumed bipolar dumbbell shape, the inner cavity consists of multi-polar fingers within an overall tripolar structure. We identify three additional planetary nebulae that show very similar morphologies and kinematics to the Owl, and propose that these constitute a new class of \textit{strigiform} (owl-like) nebulae. Common characteristics of the strigiform nebulae include a double-shell (thin outside thick) structure, low-luminosity and high-gravity central stars, the absence of a present-day stellar wind, and asymmetric inner cavities, visible in both optical and mid-infrared emission lines, that show no evidence for surrounding bright rims. The origin of the cavities is unclear, but they may constitute relics of an earlier stage of evolution when the stellar wind was active.

In order to study the nature, origin, and impact of turbulent velocity fluctuations in the ionized gas of the Orion Nebula, we apply a variety of statistical techniques to observed velocity cubes. The cubes are derived from high resolving power ($R \approx 40,000$) longslit spectroscopy of optical emission lines that span a range of ionizations. From Velocity Channel Analysis (VCA), we find that the slope of the velocity power spectrum is consistent with predictions of Kolmogorov theory between scales of 8 and 22 arcsec (0.02 to 0.05 pc). The outer scale, which is the dominant scale of density fluctuations in the nebula, approximately coincides with the autocorrelation length of the velocity fluctuations that we determine from the second order velocity structure function. We propose that this is the principal driving scale of the turbulence, which originates in the autocorrelation length of dense cores in the Orion molecular filament. By combining analysis of the non-thermal line widths with the systematic trends of velocity centroid versus ionization, we find that the global champagne flow and smaller scale turbulence each contribute in equal measure to the total velocity dispersion, with respective root-mean-square widths of 4-5 km/s. The turbulence is subsonic and can account for only one half of the derived variance in ionized density, with the remaining variance provided by density gradients in photoevaporation flows from globules and filaments. Intercomparison with results from simulations implies that the ionized gas is confined to a thick shell and does not fill the interior of the nebula.

We report the first direct imaging of protoplanetary disks in the star-forming region of Carina, the most distant, massive cluster in which disks have been imaged. Using the Atacama Large Millimeter/sub-millimeter Array (ALMA), disks are observed around two young stellar objects (YSOs) that are embedded inside evaporating gaseous globules and exhibit jet activity. The disks have an average size of 120 AU and total masses of 30 and 50 M_Jup. Given the measured masses, the minimum timescale required for planet formation (~1-2 Myr) and the average age of the Carina population (~1-4 Myr), it is plausible that young planets are present or their formation is currently ongoing in these disks. The non-detection of millimeter emission above the 4sigma threshold (~7 M_Jup) in the core of the massive cluster Trumpler~14, an area containing previously identified proplyd candidates, suggest evidence for rapid photo-evaporative disk destruction in the cluster's harsh radiation field. This would prevent the formation of giant gas planets in disks located in the cores of Carina's dense sub-clusters, whereas the majority of YSO disks in the wider Carina region remain unaffected by external photo-evaporation.

Photoionization produces supra-thermal electrons, electrons with much more energy than is found in a thermalized gas at electron temperatures characteristic of nebulae. The presence of these high energy electrons may solve the long-standing t^2/ADF puzzle, the observations that abundances obtained from recombination and collisionally excited lines do not agree, and that different temperature indicators give different results, if they survive long enough to affect diagnostic emission lines. The presence of these non-Maxwellian distribution electrons is usually designated by the term kappa. Here we use well-established methods to show that the distance over which heating rates change are much longer than the distance supra thermal electrons can travel, and that the timescale to thermalize these electrons are much shorter than the heating or cooling timescales. These estimates establish that supra thermal electrons will have disappeared into the Maxwellian velocity distribution long before they affect the collisionally excited forbidden and recombination lines that are used for deriving abundances relative to hydrogen. The electron velocity distribution in nebulae should be closely thermal.

Recent Hubble Space Telescope images have allowed the determination with unprecedented accuracy of motions and changes of shocks within the inner Orion Nebula. These originate from collimated outflows from very young stars, some within the ionized portion of the nebula and others within the host molecular cloud. We have doubled the number of Herbig-Haro objects known within the inner Orion Nebula. We find that the best-known Herbig-Haro shocks originate from a relatively few stars, with the optically visible X-ray source COUP 666 driving many of them. While some isolated shocks are driven by single collimated outflows, many groups of shocks are the result of a single stellar source having jets oriented in multiple directions at similar times. This explains the feature that shocks aligned in opposite directions in the plane of the sky are usually blue shifted because the redshifted outflows pass into the optically thick Photon Dominated Region behind the nebula. There are two regions from which optical outflows originate for which there are no candidate sources in the SIMBAD data base.

We investigate the scale dependence of fluctuations inside a realistic model of an evolving turbulent HII region and to what extent these may be studied observationally. We find that the multiple scales of energy injection from champagne flows and the photoionization of clumps and filaments leads to a flatter spectrum of fluctuations than would be expected from top-down turbulence driven at the largest scales. The traditional structure function approach to the observational study of velocity fluctuations is shown to be incapable of reliably determining the velocity power spectrum of our simulation. We find that a more promising approach is the Velocity Channel Analysis technique of Lazarian & Pogosyan (2000), which, despite being intrinsically limited by thermal broadening, can successfully recover the logarithmic slope of the velocity power spectrum to a precision of +-0.1 from high resolution optical emission line spectroscopy.

We have been able to compare with astrometric precision AstroDrizzle processed images of NGC 6720 (the Ring Nebula) made using two cameras on the Hubble Space Telescope. The time difference of the observations was 12.925 yrs. This large time-base allowed determination of tangential velocities of features within this classic planetary nebula. Individual features were measured in [N II] images as were the dark knots seen in silhouette against background nebular [O III] emission. An image magnification and matching technique was also used to test the accuracy of the usual assumption of homologous expansion. We found that homologous expansion does apply, but the rate of expansion is greater along the major axis of the nebula, which is intrinsically larger than the minor axis. We find that the dark knots expand more slowly that the nebular gas, that the distance to the nebula is 720 pc +/-30%, and the dynamic age of the Ring Nebula is about 4000 yrs. The dynamic age is in agreement with the position of the central star on theoretical curves for stars collapsing from the peak of the Asymptotic Giant Branch to being white dwarfs.

This is a summary of the 2013 release of the plasma simulation code Cloudy. Cloudy models the ionization, chemical, and thermal state of material that may be exposed to an external radiation field or other source of heating, and predicts observables such as emission and absorption spectra. It works in terms of elementary processes, so is not limited to any particular temperature or density regime. This paper summarizes advances made since the last major review in 1998. Much of the recent development has emphasized dusty molecular environments, improvements to the ionization / chemistry solvers, and how atomic and molecular data are used. We present two types of simulations to demonstrate the capability of the code. We consider a molecular cloud irradiated by an X-ray source such as an Active Nucleus and show how treating EUV recombination lines and the full SED affects the observed spectrum. A second example illustrates the very wide range of particle and radiation density that can be considered.

We have performed a detailed analysis of the Ring Nebula (NGC 6720) using HST WFC3 images and derived a new 3-D model. Existing high spectral resolution spectra played an important supplementary role in our modeling. It is shown that the Main Ring of the nebula is an ionization-bounded irregular non-symmetric disk with a central cavity and perpendicular extended lobes pointed almost towards the observer. The faint outer halos are determined to be fossil radiation, i.e. radiation from gas ionized in an earlier stage of the nebula when it was not ionization bounded. The narrow-band WFC3 filters that isolate some of the emission-lines are affected by broadening on their short wavelength side and all the filters were calibrated using ground-based spectra. The filter calibration results are presented in an appendix.

We have performed a detailed analysis of the electron temperature and density in the the Ring Nebula using the calibrated HST WFC3 images described in the preceding paper. The electron temperature (Te) determined from [N II] and [O III] rises slightly and monotonically towards the central star. The observed equivalent width (EW) in the central region indicates that Te rises as high as 13000 K. In contrast, the low EW's in the outer regions are largely due to scattered diffuse Galactic radiation by dust. The images allowed determination of unprecedented small scale variations in Te. These variations indicate that the mean square area temperature fluctuations are significantly higher than expected from simple photoionization. The power producing these fluctuations occurs at scales of less than 3.5E15 cm. This scale length provides a strong restriction on the mechanism causing the large t^2 values observed.

Photoevaporating protoplanetary discs (proplyds) in the vicinity of hot massive stars, such as those found in Orion, are important objects of study for the fields of star formation, early disc evolution, planetary formation, and H II region astrophysics. Their element abundances are largely unknown, unlike those of the main-sequence stars or the host Orion nebula. We present a spectroscopic analysis of the Orion proplyd HST 10, based on integral field observations with the Very Large Telescope/FLAMES fibre array with 0.31" x 0.31" spatial pixels. The proplyd and its vicinity are imaged in a variety of emission lines across a 6.8" x 4.3" area. The reddening, electron density and temperature are mapped out from various line diagnostics. The abundances of helium, and eight heavy elements are measured relative to hydrogen using the direct method based on the [O III] electron temperature. The abundance ratios of O/H and S/H are derived without resort to ionization correction factors. We construct dynamic photoevaporation models of HST 10 with the Cloudy microphysics code that validate the oxygen and sulfur abundances. With the exception of [O I] 6300-A and [S II] 4069-A, the model fit is satisfactory for all spectral lines arising from the proplyd. The models show that the classic ionization correction factor for neon significantly underestimates (0.4 dex) this element's abundance in the low ionization conditions of HST 10. Apart from iron, whose gas-phase abundance is ~0.3 dex lower than in the local Orion nebula, most other elements in the proplyd do not show substantially different gas-phase abundances from the nebula. The abundances of carbon, oxygen and neon in HST 10 are practically the same as those in B-type stars in Orion.

The optical [N I] doublet near 5200 {\AA} is anomalously strong in a variety of emission-line objects. We compute a detailed photoionization model and use it to show that pumping by far-ultraviolet (FUV) stellar radiation previously posited as a general explanation applies to the Orion Nebula (M42) and its companion M43; but, it is unlikely to explain planetary nebulae and supernova remnants. Our models establish that the observed nearly constant equivalent width of [N I] with respect to the dust-scattered stellar continuum depends primarily on three factors: the FUV to visual-band flux ratio of the stellar population; the optical properties of the dust; and the line broadening where the pumping occurs. In contrast, the intensity ratio [N I]/H{\beta} depends primarily on the FUV to extreme-ultraviolet ratio, which varies strongly with the spectral type of the exciting star. This is consistent with the observed difference of a factor of five between M42 and M43, which are excited by an O7 and B0.5 star respectively. We derive a non-thermal broadening of order 5 km/s for the [N I] pumping zone and show that the broadening mechanism must be different from the large-scale turbulent motions that have been suggested to explain the line-widths in this H II region. A mechanism is required that operates at scales of a few astronomical units, which may be driven by thermal instabilities of neutral gas in the range 1000 to 3000 K. In an appendix, we describe how collisional and radiative processes are treated in the detailed model N I atom now included in the Cloudy plasma code.

We present results from integral field spectroscopy with PMAS. The observed field contains: five protoplanetary discs (also known as proplyds), the high-velocity jet HH 514 and a bowshock. Spatial distribution maps are obtained for different emission line fluxes, the c(H{\beta}) coefficient, electron densities and temperatures, ionic abundances of different ions from collisionally excited lines (CELs), C2+ and O2+ abundances from recombination lines (RLs) and the abundance discrepancy factor of O2+, ADF(O2+). We find that collisional de-excitation has a major influence on the line fluxes in the proplyds. If this is not properly accounted for then physical conditions deduced from commonly used line ratios will be in error, leading to unreliable chemical abundances for these objects. We obtain the intrinsic emission of the proplyds 177-341, 170-337 and 170-334 by a direct subtraction of the background emission, though the last two present some background contamination due to their small sizes. A detailed analysis of 177-341 spectra reveals the presence of high-density gas (3.8\times10^5 cm^-3) in contrast to the typical values observed in the background gas of the nebula (3800 cm^-3). We also explore how the background subtraction could be affected by the possible opacity of the proplyd. We construct a physical model for the proplyd 177-341 finding a good agreement between the predicted and observed line ratios. Finally, we find that the use of reliable physical conditions returns an ADF(O2+) about zero for the intrinsic spectra of 177-341, while the background emission presents the typical ADF(O2+) observed in the Orion Nebula. We conclude that the presence of high-density ionized gas is severely affecting the abundances determined from CELs and, therefore, those from RLs should be considered as a better approximation to the true abundances.

We present the results of radiation-magnetohydrodynamic simulations of the expansion of HII regions and surrounding photodissociation regions in turbulent, magnetised, molecular clouds on scales of up to 4 parsecs, including the effects of ionising and non-ionising ultraviolet radiation and x rays from young star clusters. We find that HII region expansion reduces the disordered component of the B field, imposing a large-scale order on the field around its border. The field in the neutral gas lies along the ionisation front, while the field in the ionised gas tends to be perpendicular to this. The highest pressure compressed neutral/molecular gas is driven towards approximate equipartition between thermal/magnetic/turbulent energy densities, whereas lower pressure neutral/molecular gas divides into quiescent, magnetically dominated regions, and, on the other hand, turbulent, demagnetised regions. The ionised gas shows approximate thermal/turbulent equipartition, but with magnetic energy densities 1 to 3 orders of magnitude lower. A high velocity dispersion (approx 8 km/s) is maintained in the ionised gas throughout our simulations, despite the mean expansion velocity being significantly lower. The B field does not significantly brake the HII region expansion on the length and timescales accessible to our simulations, but it does tend to suppress the small-scale fragmentation and radiation-driven implosion of neutral/molecular gas that forms globules and pillars at the edge of the HII region. However, the relative luminosity of ionising and non-ionising radiation has a much larger influence than the presence or absence of the B field. When the radiation field is relatively soft (as in the case of a lower mass cluster, with earliest spectral type of B0.5), then fragmentation is less vigorous and a thick, relatively smooth PDR forms. Movies available at http://youtube.com/user/divBequals0

In planetary nebulae, abundances of oxygen and other heavy elements derived from optical recombination lines are systematically higher than those derived from collisionally excited lines. We investigate the hypothesis that the destruction of solid bodies may produce pockets of cool, high-metallicity gas that could explain these abundance discrepancies. Under the assumption of maximally efficient radiative ablation, we derive two fundamental constraints that the solid bodies must satisfy in order that their evaporation during the planetary nebula phase should generate a high enough gas phase metallicity. A local constraint implies that the bodies must be larger than tens of meters, while a global constraint implies that the total mass of the solid body reservoir must exceed a few hundredths of a solar mass. This mass greatly exceeds the mass of any population of comets or large debris particles expected to be found orbiting evolved low- to intermediate-mass stars. We therefore conclude that contemporaneous solid body destruction cannot explain the observed abundance discrepancies in planetary nebulae. However, similar arguments applied to the sublimation of solid bodies during the preceding asymptotic giant branch (AGB) phase do not lead to such a clear-cut conclusion. In this case, the required reservoir of volatile solids is only one ten-thousandth of a solar mass, which is comparable to the most massive debris disks observed around solar-type stars, implying that this mechanism may contribute to abundance discrepancies in at least some planetary nebulae, so long as mixing of the high metallicity gas is inefficient.

We compare the accuracy of various methods for determining the transfer of the diffuse Lyman continuum in HII regions, by comparing them with a high-resolution discrete-ordinate integration. We use these results to suggest how, in multidimensional dynamical simulations, the diffuse field may be treated with acceptable accuracy without requiring detailed transport solutions. The angular distribution of the diffuse field derived from the numerical integration provides insight into the likely effects of the diffuse field for various material distributions.

The Orion Bar is one of the nearest and best-studied photodissociation or photon-dominated regions (PDRs). Observations reveal the presence of H2 lines from vibrationally or rotationally excited upper levels that suggest warm gas temperatures (400 to 700 K). However, standard models of PDRs are unable to reproduce such warm rotational temperatures. In this paper we attempt to explain these observations with new comprehensive models which extend from the H+ region through the Bar and include the magnetic field in the equation of state. We adopt the model parameters from our previous paper which successfully reproduced a wide variety of spectral observations across the Bar. In this model the local cosmic-ray density is enhanced above the galactic background, as is the magnetic field, and which increases the cosmic-ray heating elevating the temperature in the molecular region. The pressure is further enhanced above the gas pressure in the H+ region by the momentum transferred from the absorbed starlight. Here we investigate whether the observed H2 lines can be reproduced with standard assumptions concerning the grain photoelectric emission. We also explore the effects due to the inclusion of recently computed H2 + H2, H2 + H and H2 + He collisional rate coefficients.

We present the first three-dimensional radiation-magnetohydrodynamic simulations of the photoionisation of a dense, magnetised molecular globule by an external source of ultraviolet radiation. We find that, for the case of a strong ionising field, significant deviations from the non-magnetic evolution are seen when the initial magnetic field threading the globule has an associated magnetic pressure that is greater than one hundred times the gas pressure. In such a strong-field case, the photoevaporating globule will adopt a flattened or "curled up" shape, depending on the initial field orientation, and magnetic confinement of the ionised photoevaporation flow can lead to recombination and subsequent fragmentation during advanced stages of the globule evolution. We find suggestive evidence that such magnetic effects may be important in the formation of bright, bar-like emission features in H II regions. We include simple but realistic fits to heating and cooling rates in the neutral and molecular gas in the vicinity of a high-mass star cluster and show that the frequently used isothermal approximation can lead to an overestimate of the importance of gravitational instability in the radiatively imploded globule. For globules within 2 parsecs of a high-mass star cluster, we find that heating by stellar x rays prevents the molecular gas from cooling below 50 K.

We have combined recent Hubble Space Telescope WFPC2 images in the [O III] 5007 and [N II] 6583 lines with similar images made 9.557 years earlier to determine the motion of the Ring Nebula within the plane of the sky. Scaled ratio images argue for homologous expansion, that is, larger velocities scale with increasing distance from the central star. The rather noisy pattern of motion of individual features argues for the same conclusion and that the silhouetted knots move at the same rate as the surrounding gas. These tangential velocities are combined with information from a recent high resolution radial velocity study to determine a dynamic distance, which is in basic agreement with the distance determined from the parallax of the central star. We have also obtained very high signal to noise ratio moderate resolution spectra (9.4 Angstrom) along the major and minor axes of the nebula and from this determined the electron temperatures and density in the multiple ionization zones present. These results confirm the status of the Ring Nebula as one of the older planetary nebulae, with a central star transitioning to the white dwarf cooling curve. (Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA Contract No. NAS 5-26555 and the San Pedro Martir Observatory operated by the Universidad Nacional Autonoma de Mexico.)

The three dimensional structure of the brightest part of the Orion Nebula is assessed in the light of published and new data. We find that the widely accepted model of a concave blister of ionized material needs to be altered in the southwest direction from the Trapezium, where we find that the Orion-S feature is a separate cloud of very optically thick molecules within the body of ionized gas, which is probably the location of the multiple embedded sources that produce the outflows that define the Orion-S star formation region. Evidence for this cloud comes from the presence of H2CO lines in absorption in the radio continuum and discrepancies in the extinction derived from radio-optical and optical only emission. We present an equilibrium Cloudy model of the Orion-S cloud, which successfully reproduces many observed properties of this feature. We also report the discovery of an open-sided shell of [O III] surrounding the Trapezium stars, revealed through emission line ratio images and the onset of radiation shadows beyond some proplyds. We show that the observed properties of the shell are consistent with it being a stationary structure, produced by shock interactions between the ambient nebular gas and the high-velocity wind from theta^1 Ori C. We examine the implications of the recently published evidence for a large blueshifted velocity of theta^1 Ori C with respect to the Orion Molecular Cloud, which could mean that this star has only recently begun to photoionize the Orion Nebula. We show that current observations of the Nebula do not rule out such a possibility, so long as the ionization front has propagated into a pre-existing low-density region. In addition, a young age for the Nebula would help explain the presence of nearby proplyds with a short mass-loss timescale to photoablation.

We have used widely spaced in time Hubble Space Telescope images to determine tangential velocities of features associated with outflows from young stars. These observations were supplemented by groundbased telescope spectroscopy and from the resultant radial velocities, space velocities were determined for many outflows. Numerous new moving features were found and grouped into known and newly assigned Herbig Haro objects. It was found that stellar outflow is highly discontinuous, as frequently is the case, with long-term gaps of a few hundred years and that these outflow periods are marked by staccato bursts over periods of about ten years. Although this has been observed in other regions, the Orion Nebula Cluster presents the richest display of this property. Most of the large scale Herbig Haro objects in the brightest part of the Orion Nebula appear to originate from a small region northeast of the strong Orion-S radio and infrared sources. With the possible exception of HH 203, we are not able to identify specific stellar sources, but do identify candidate sources for several other bright Herbig Haro objects. We find that there are optical features in the BN-KL region that can be related to the known large scale outflow that originates there. We find additional evidence for this outflow originating 500 to 1000 years ago.

We present an atlas of three-dimensional (position-position-velocity) spectra of the Orion Nebula in optical emission lines from a variety of different ionization stages: [O I] 6300, [S II] 6716,6731, [N II] 6584, [S III] 6312, H alpha 6563, and [O III] 5007. These transitions provide point to point information about the physical structure and kinematics of the nebula at an effective resolution of 3'' x 2'' x 10 km/s, clearly showing the large scale behavior of the ionized gas and the presence of localized phenomena such as Herbig-Haro outflows. As an example application of the atlas, we present a statistical analysis of the widths of the H alpha, [O III], and [N II] lines that permits a determination of the mean electron temperature in the nebula of (9200 +/- 400) K. We also find, in contradiction to previous claims, that the non-thermal line broadening is not significantly different between recombination lines and collisional lines.

The hydrogen ionization and dissociation front around an ultraviolet radiation source should merge when the ratio of ionizing photon flux to gas density is sufficiently low and the spectrum is sufficiently hard. This regime is particularly relevant to the molecular knots that are commonly found in evolved planetary nebulae, such as the Helix Nebula, where traditional models of photodissociation regions have proved unable to explain the high observed luminosity in H_2 lines. In this paper we present results for the structure and steady-state dynamics of such advection-dominated merged fronts, calculated using the Cloudy plasma/molecular physics code. We find that the principal destruction processes for H_2 are photoionization by extreme ultraviolet radiation and charge exchange reactions with protons, both of which form H_2^+, which rapidly combines with free electrons to undergo dissociative recombination. Advection moves the dissociation front to lower column densities than in the static case, which vastly increases the heating in the partially molecular gas due to photoionization of He^0, H_2, and H^0. This causes a significant fraction of the incident bolometric flux to be re-radiated as thermally excited infrared H_2 lines, with the lower excitation pure rotational lines arising in 1000 K gas and higher excitation H_2 lines arising in 2000 K gas, as is required to explain the H_2 spectrum of the Helix cometary knots.

We critically examine a scenario for the enrichment of the interstellar medium (ISM) in which supernova ejecta follow a long (10^8 yr) journey before falling back onto the galactic disk in the form of metal-rich ``droplets'', These droplets do not become fully mixed with the interstellar medium until they become photoionized in HII regions. We investigate the hypothesis that the photoionization of these highly metallic droplets can explain the observed ``abundance discrepancy factors'' (ADFs), which are found when comparing abundances derived from recombination lines and from collisionally excited lines, both in Galactic and extragalactic HII regions. We derive bounds of 10^{13}--10^{15} cm on the droplet sizes inside HII regions in order that (1) they should not have already been detected by direct imaging of nearby nebulae, and (2) they should not be too swiftly destroyed by diffusion in the ionized gas. From photoionization modelling we find that, if this inhomogeneous enrichment scenario holds, then the recombination lines strongly overestimate the metallicities of the fully mixed HII regions. The abundances derived from collisionally excited lines also suffer some bias, although to a much lesser extent. In the absence of any recipe for correcting these biases, we recommend the discarding of all objects showing large ADFs from studies of galactic chemical evolution. These biases must also be kept in mind when comparing the galactic abundance gradients for elements derived from recombination lines with those derived from collisionally excited lines. Finally, we propose a set of observations that could be undertaken to test our scenario and improve our understanding of element mixing in the ISM.

We argue that the discrepancies observed in HII regions between abundances derived from optical recombination lines (ORLs) and collisionally excited lines (CELs) might well be the signature of a scenario of the enrichment of the interstellar medium (ISM) proposed by Tenorio-Tagle (1996). In this scenario, the fresh oxygen released during massive supernova explosions is confined within the hot superbubbles as long as supernovae continue to explode. Only after the last massive supernova explosion, the metal-rich gas starts cooling down and falls on the galaxy within metal-rich droplets. Full mixing of these metal-rich droplets and the ISM occurs during photoionization by the next generations of massive stars. During this process, the metal-rich droplets give rise to strong recombination lines of the metals, leading to the observed ORL-CEL discrepancy. (The full version of this work is submitted to Astronomy and Astrophysics.)

[Abridged] We use new HST and archived images to clarify the nature of the knots in the Helix Nebula. We employ published far infrared spectrophotometry and existing 2.12 micron images to establish that the population distribution of the lowest ro-vibrational states of H2 is close to the distribution of a gas in LTE at 988 +- 119 K. We derive a total flux from the nebula in H2 lines and compare this with the power available from the central star for producing this radiation. We establish that neither soft X-rays nor FUV radiation has enough energy to power the H2 radiation, only the stellar EUV radiation shortward of 912 Angstrom does. Advection of material from the cold regions of the knots produces an extensive zone where both atomic and molecular hydrogen are found, allowing the H2 to directly be heated by Lyman continuum radiation, thus providing a mechanism that can explain the excitation temperature and surface brightness of the cusps and tails. New images of the knot 378-801 reveal that the 2.12 micron cusp and tail lie immediately inside the ionized atomic gas zone. This firmly establishes that the "tail" structure is an ionization bounded radiation shadow behind the optically thick core of the knot. A unique new image in the HeII 4686 Angstrom line fails to show any emission from knots that might have been found in the He++ core of the nebula. We also re-examined high signal-to-noise ratio ground-based telescope images of this same inner region and found no evidence of structures that could be related to knots.

Multiple optical outflows are known to exist in the vicinity of the active star formation region called Orion-South (Orion-S). We have mapped the velocity of low ionization features in the brightest part of the Orion Nebula, including Orion-S, and imaged the entire nebula with the Hubble Space Telescope. These new data, combined with recent high resolution radio maps of outflows from the Orion-S region, allow us to trace the origin of the optical outflows. It is confirmed that HH 625 arises from the blueshifted lobe of the CO outflow from 136-359 in Orion-S while it is likely that HH 507 arises from the blueshifted lobe of the SiO outflow from the nearby source 135-356. It is likely that redshifted lobes are deflected within the photon dominated region behind the optical nebula. This leads to a possible identification of a new large shock to the southwest from Orion-S as being driven by the redshifted CO outflow arising from 137-408. The distant object HH 400 is seen to have two even further components and these all are probably linked to either HH 203, HH 204, or HH 528. Distant shocks on the west side of the nebula may be related to HH 269. The sources of multiple bright blueshifted Herbig-Haro objects (HH 202, HH 203, HH 204, HH 269, HH 528) remain unidentified, in spite of earlier claimed identifications. Some of this lack of identification may arise from the fact that deflection in radial velocity can also produce a change in direction in the plane of the sky. The best way to resolve this open question is through improved tangential velocities of low ionization features arising where the outflows first break out into the ionized nebula.

High-dispersion echelle spectroscopy in optical forbidden lines of O^0, S^+, and S^2+ is used to construct velocity-resolved images and electron density maps of the inner region of the Orion nebula with a resolution of 10 km s-1 x 3" x 2". Among the objects and regions newly discovered in this study are (1) the Diffuse Blue Layer: an extended layer of moderately blue-shifted, low-density, low-ionization emission in the southeast region of the nebula; (2) the Red Bay: a region to the east of the Trapezium where the usual correlation between velocity and ionization potential is very weak, and where the emitting layer is very thick; and (3) HH 873: a new redshifted jet to the southwest of the Trapezium.

We present numerical radiation-hydrodynamic simulations of the evolution of HII regions formed in an inhomogeneous medium resulting from turbulence simulations. We find that the filamentary structure of the underlying density distribution produces a highly irregular shape for the ionized region, in which the ionization front escapes to large distances in some directions within 80,000 years. In other directions, on the other hand, neutral gas in the form of dense globules persists within 1 parsec of the central star for the full duration of our simulation (400,000 years). Divergent photoablation flows from these globules maintain a root-mean-squared velocity in the ionized gas that is close to the ionized sound speed. Simulated images in optical emission lines show morphologies that are in strikingly detailed agreement with those observed in real HII regions.

This review covers the dynamic processes that are important in the evolution and structure of galactic HII regions, concentrating on an elementary presentation of the physical concepts and recent numerical simulations of HII region evolution in a non-uniform medium. The contents are as follows: (1) The equations (Euler equations; Radiative transfer; Rate equations; How to avoid the dynamics; How to avoid the atomic physics). (2) Physical concepts (Static photoionization equilibrium; Ionization front propagation; Structure of a D-type front; Photoablation flows; Other ingredients - Stellar winds, Radiation pressure, Magnetic fields, Instabilities). (3) HII region evolution (Early phases: hypercompact and ultracompact regions; Later phases: compact and extended regions; Clumps and turbulence).

We present hydrodynamical simulations of the photoevaporation of a cloud with large-scale density gradients, giving rise to an ionized, photoevaporation flow. The flow is found to be approximately steady during the large part of its evolution, during which it can resemble a "champagne flow" or a "globule flow" depending on the curvature of the ionization front. The distance from source to ionization front and the front curvature uniquely determine the structure of the flow, with the curvature depending on the steepness of the lateral density gradient in the neutral cloud. We compare these simulations with both new and existing observations of the Orion nebula and find that a model with a mildly convex ionization front can reproduce the profiles of emission measure, electron density, and mean line velocity for a variety of emitting ions on scales of 10^{17} to 10^{18} cm. The principal failure of our model is that we cannot explain the large observed widths of the [O I] 6300 Angstrom line that forms at the ionization front.

We have conducted a combined observational and theoretical investigation of the ubiquitous knots in the Helix Nebula (NGC 7293). We have constructed a combined hydrodynamic+radiation model for the ionized portion of these knots and have accurately calculated a static model for their molecular regions. Imaging observations in optical emission lines were made with the Hubble Space Telescope's STIS spectrograph, operating in a "slitless" mode, complemented by WFPC2 images in several of the same lines. The NICMOS camera was used to image the knots in molecular hydrogen. These observations, when combined with other studies of molecular hydrogen and CO provide a complete characterization of the knots. They possess dense molecular cores of densities about 1,000,000 per cubic centimeter surrounded on the central star side by a zone of hot molecular hydrogen. The temperature of the molecular hydrogen emitting layer defies explanation either through detailed calculations for radiative equilibrium or for simplistic calculations for shock excitation. Further away from the core is the ionized zone, whose peculiar distribution of emission lines is explained by the expansion effects of material flowing through this region. The shadowed region behind the core is the source of most of the CO emission from the knot and is of the low temperature expected for a radiatively heated molecular region.

We present a method for including steady-state gas flows in the plasma physics code Cloudy, which was previously restricted to modeling static configurations. The numerical algorithms are described in detail, together with an example application to plane-parallel ionization-bounded HII regions. As well as providing the foundation for future applications to more complex flows, we find the following specific results regarding the effect of advection upon ionization fronts in HII regions: 1. Significant direct effects of advection on the global emission properties occur only when the ionization parameter is lower than is typical for HII regions. 2. The overheating of partially ionized gas in the front is not large, even for supersonic (R-type) fronts. 3. The most significant morphological signature of advective fronts is an electron density spike that occurs at the ionization front. Observational evidence for such a spike is found in images of the Orion bar. 4. Plane-parallel, weak-D fronts are found to show at best a shallow correlation between mean velocity and ionization potential for optical emission lines even when the flow velocity closely approaches the ionized sound speed.

Presented here are high angular resolution MERLIN 5 GHz (6 cm) continuum observations of the binary proplyd system, LV 1 in the Orion nebula, which consists of proplyd 168--326SE and its binary proplyd companion 168--326NW (separation 0.4 arcsec). Accurate astrometric alignment allows a detailed comparison between these data and published HST PC Halpha and [Oiii] images. Thermal radio sources coincide with the two proplyds and originate in the ionized photoevaporating flows seen in the optical emission lines. Flow velocities of approx 50 km/s from the ionized proplyd surfaces and \geq 100 km/s from a possible micro-jet have been detected using the Manchester Echelle spectrometer. A third radio source is found to coincide with a region of extended, high excitation, optical line emission that lies between the binary proplyds 168--326SE/326NW . This is modelled as a bowshock due to the collision of the photoevaporating flows from the two proplyds. Both a thermal and a non-thermal origin for the radio emission in this collision zone are considered.

We describe an efficient Monte Carlo algorithm for a restricted class of scattering problems in radiation transfer. This class includes many astrophysically interesting problems, including the scattering of ultraviolet and visible light by grains. The algorithm correctly accounts for multiply-scattered light. We describe the algorithm, present a number of important optimizations, and explicity show how the algorithm can be used to estimate quantities such as the emergent and mean intensity. We present two test cases, examine the importance of the optimizations, and show that this algorithm can be usefully applied to optically-thin problems, a regime sometimes considered limited to explicit single-scattering plus attenuation approximations.

We explain the Ha emission of the cometary knots in the Helix Nebula (NGC 7293) with an analytical model that describes the emission of the head of the globules as a photoevaporated flow produced by the incident ionizing radiation of the central star.We compare these models with the Ha emission obtained from the HST archival images of the Helix Nebula. From a comparison of the Ha emission with the predictions of the analytical model we obtain a rate of ionizing photons from the central star of about 5e45 s^-1, which is consistent with estimates based on the total Hb flux of the nebula. We also model the tails of the cometary knots as a photoevaporated wind from a neutral shadow region produced by the diffuse ionizing photon field of the nebula. A comparison with the HST images allows us to obtain a direct determination of the value of the diffuse ionizing flux. We compare the ratio of diffuse to direct stellar flux as a function of radius inside an HII region with those obtained from the observational data through the analytical tail and head wind model. The agreement of this model with the values determined from the observations of the knots is excellent.

We present the results of spectroscopy of four bright proplyds in the Orion Nebula obtained at a velocity resolution of 6 km/s. After careful isolation of the proplyd spectra from the confusing nebular radiation, the emission line profiles are compared with those predicted by realistic dynamic/photoionization models of the objects. The spectral line widths show a clear correlation with ionization potential, which is consistent with the free expansion of a transonic, ionization-stratified, photoevaporating flow. Fitting models of such a flow simultaneously to our spectra and HST emission line imaging provides direct measurements of the proplyd size, ionized density and outflow velocity. These measurements confirm that the ionization front in the proplyds is approximately D-critical and provide the most accurate and robust estimate to date of the proplyd mass loss rate. Values of 0.7E-6 to 1.5E-6 Msun/year are found for our spectroscopic sample, although extrapolating our results to a larger sample of proplyds implies that 0.4E-6 Msun/year is more typical of the proplyds as a whole. In view of the reported limits on the masses of the circumstellar disks within the proplyds, the length of time that they can have been exposed to ionizing radiation should not greatly exceed 10,000 years - a factor of 30 less than the mean age of the proplyd stars. We review the various mechanisms that have been proposed to explain this situation, and conclude that none can plausibly work unless the disk masses are revised upwards by a substantial amount.

Many low-mass stars in the Orion nebula are associated with very compact (~= 1 arcsec) emission knots, known variously as proplyds, PIGs or LV knots. Some of these knots are teardrop-shaped, with ``tails'' pointing away from the massive star Theta 1 Ori C, which is the principal exciting star of the nebula. We discuss models of such knots, which invoke the interaction of the fast stellar wind from Theta 1 Ori C with a transonic photoevaporated flow from the surface of an accretion disk around a young low-mass star. We review previous analytic work and compare the results of the model with the observed brightnesses, morphologies and emission line profiles of the knots, as well as presenting new results from numerical hydrodynamical simulations.

The Doppler shifts of optical emission lines which have been scattered by surrounding dust and electrons can provide useful information about the kinematics, geometry and physical conditions of astrophysical flows. In principle, the scatterers can provide views of the line-emitting gas from different directions, allowing the 3-d velocity field of the emitting gas to be determined and revealing sources which are hidden from direct view. Unfortunately, as in a Hall of Mirrors, the resultant multiple images can be confusing and hard to interpret. In general, the scattering geometry will be unknown, which makes it difficult to disentangle the effect of the motion of the scatterers from that of the motion of the emitting source. In this situation, spectropolarimetric observations can be a great help, since the scattered light will be partially polarised to a degree dependent on the angle of scattering and on the details of the scattering process. Analysis of the polarisation profiles of the scattered emission lines can then determine the relative orientation and velocity of the source and the scatterers. Such techniques are applied to the upstream scattered emission lines in the HH 1 jet, where it is shown that spectropolarimetric measurements would allow an independent determination of the bowshock velocity with respect to the upstream matter, which would test ``multiple outburst'' models of the jet.

Detailed comparisons are made between observations of scattered light upstream of the head of the HH~1 jet and predictions of simple scattering models. It is shown that, in order to unambiguously determine the velocity of the head of the jet (bow shock) with respect to the upstream dust, existing spectroscopic observations are insufficient and that spectropolarimetric observations of the scattered light are necessary. Such an independent measure of the bow shock velocity is important in order to test ``multiple outflow'' theories of Herbig-Haro jets. It is also shown 2that the scattering dust must have a very forward-throwing scattering phase function ($\langle\cos\theta\rangle\msim 0.7$) and slight evidence is found for a dust-gas ratio that is higher than average.

A general theoretical framework is developed for interpreting spectropolarimetric observations of optically thin emission line scattering from small dust particles. Spatially integrated and spatially resolved line profiles of both scattered intensity and polarization are calculated analytically from a variety of simple kinematic models. These calculations will provide a foundation for further studies of emission line scattering from dust and electrons in such diverse astrophysical environments as Herbig-Haro objects, symbiotic stars, starburst galaxies and active galactic nuclei.