We extend our earlier work on turbulence-induced relative velocity between equal-size particles (Pan and Padoan, Paper I) to particles of arbitrarily different sizes. The Pan and Padoan (PP10) model shows that the relative velocity between different particles has two contributions, named the generalized shear and acceleration terms, respectively. The generalized shear term represents the particles' memory of the spatial flow velocity difference across the particle distance in the past, while the acceleration term is associated with the temporal flow velocity difference on individual particle trajectories. Using the simulation of Paper I, we compute the root-mean-square relative velocity, <w^2>^1/2, as a function of the friction times, tau_p1 and tau_p2, of the two particles, and show that the PP10 prediction is in satisfactory agreement with the data, confirming its physical picture. For a given tau_p1 below the Lagrangian correlation time of the flow, T_L, <w^2>^1/2 as a function of tau_p2 shows a dip at tau_p2~tau_p1, indicating tighter velocity correlation between similar particles. Defining a ratio f=tau_pl/tau_ph, with tau_pl and tau_ph the friction times of the smaller and larger particles, we find that <w^2>^1/2 increases with decreasing f due to the generalized acceleration contribution, which dominates at f<1/4. At a fixed f, our model predicts that <w^2>^1/2 scales as tau_ph^1/2 for tau_p,h in the inertial range of the flow, stays roughly constant for T_L <tau_ph < T_L/f, and finally decreases as tau_ph^-1/2 for tau_ph>>T_L/f. The acceleration term is independent of the particle distance, r, and thus reduces the r-dependence of <w^2>^1/2 in the bidisperse case.

Motivated by its important role in the collisional growth of dust particles in protoplanetary disks, we investigate the probability distribution function (PDF) of the relative velocity of inertial particles suspended in turbulent flows. Using the simulation from our previous work, we compute the relative velocity PDF as a function of the friction timescales, tau_p1 and tau_p2, of two particles of arbitrary sizes. The friction time of particles included in the simulation ranges from 0.1 tau_eta to 54T_L, with tau_eta and T_L the Kolmogorov time and the Lagrangian correlation time of the flow, respectively. The relative velocity PDF is generically non-Gaussian, exhibiting fat tails. For a fixed value of tau_p1, the PDF is the fattest for equal-size particles (tau_p2~tau_p1), and becomes thinner at both tau_p2<tau_ p1 and tau_p2>tau_p1. Defining f as the friction time ratio of the smaller particle to the larger one, we find that, at a given f in 1/2<f<1, the PDF fatness first increases with the friction time, tau_ph, of the larger particle, peaks at tau_ph~tau_eta, and then decreases as tau_ph increases further. For 0<f<1/4, the PDF shape becomes continuously thinner with increasing tau_ph. The PDF is nearly Gaussian only if tau_ph is sufficiently large (>>T_L). These features are successfully explained by the Pan & Padoan model. Using our simulation data and some simplifying assumptions, we estimated the fractions of collisions resulting in sticking, bouncing, and fragmentation as a function of the dust size in protoplanerary disks, and argued that accounting for non-Gaussianity of the collision velocity may help further alleviate the bouncing barrier problem.

We conduct a comprehensive theoretical and numerical investigation of the pollution of pristine gas in turbulent flows, designed to provide new tools for modeling the evolution of the first generation of stars. The properties of such Population III (Pop III) stars are thought to be very different than later generations, because cooling is dramatically different in gas with a metallicity below a critical value Z_c, which lies between ~10^-6 and 10^-3 solar value. Z_c is much smaller than the typical average metallicity, <Z>, and thus the mixing efficiency of the pristine gas in the interstellar medium plays a crucial role in the transition from Pop III to normal star formation. The small critical value, Z_c, corresponds to the far left tail of the probability distribution function (PDF) of the metallicity. Based on closure models for the PDF formulation of turbulent mixing, we derive equations for the fraction of gas, P, lying below Z_c, in compressible turbulence. Our simulation data shows that the evolution of the fraction P can be well approximated by a generalized self-convolution model, which predicts dP/dt = -n/tau_con P (1-P^(1/n)), where n is a measure of the locality of the PDF convolution and the timescale tau_con is determined by the rate at which turbulence stretches the pollutants. Using a suite of simulations with Mach numbers ranging from M = 0.9 to 6.2, we provide accurate fits to n and tau_con as a function of M, Z_c/<Z>, and the scale, L_p, at which pollutants are added to the flow. For P>0.9, mixing occurs only in the regions surrounding the pollutants, such that n=1. For smaller P, n is larger as mixing becomes more global. We show how the results can be used to construct one-zone models for the evolution of Pop III stars in a single high-redshift galaxy, as well as subgrid models for tracking the evolution of the first stars in large cosmological simulations.

The first generation of stars had very different properties than later stellar generations, as they formed from a "pristine" gas that was free of heavy elements. Normal star formation took place only after the first stars polluted the surrounding turbulent interstellar gas, increasing its local heavy element concentration, Z, beyond a critical value, Z_c (10^-8 < Z_c <10^-5). Motivated by this astrophysical problem, we investigate the fundamental physics of the pollution of pristine fluid elements in isotropic compressible turbulence. Turbulence stretches the pollutants, produces concentration structures at small scales, and brings the pollutants and the unpolluted flow in closer contact. Our theoretical approach employs the probability distribution function (PDF) method for turbulent mixing. We adopt three PDF closure models and derive evolution equations for the pristine fraction from the models. To test and constrain the theoretical models, we conduct numerical simulations for decaying passive scalars in isothermal turbulent flows with Mach numbers of 0.9 and 6.2, and compute the mass fraction, P(Z_c, t), of the flow with Z < Z_c. In the Mach 0.9 flow, the evolution of P(Z_c, t)$ is well described by a continuous convolution model and dP(Z_c, t)/dt = P(Z_c, t) ln[P(Z_c, t)]/tau_con, if the mass fraction of the polluted flow is larger than ~ 0.1. If the initial pollutant fraction is smaller than ~ 0.1, an early phase exists during which the pristine fraction follows an equation from a nonlinear integral model: dP(Z_c, t)/dt = P(Z_c, t) [P(Z_c, t)-1]/tau_int. The timescales tau_con and tau_int are measured from our simulations. When normalized to the flow dynamical time, the decay of P(Z_ c, t) in the Mach 6.2 flow is slower than at Mach 0.9, and we show that P(Z_c, t) in the Mach 6.2 flow can be well fit using a formula from a generalized version of the self-convolution model.

We study clustering of inertial particles in turbulent flows and discuss its applications to dust particles in protoplanetary disks. Using numerical simulations, we compute the radial distribution function (RDF), which measures the probability of finding particle pairs at given distances, and the probability density function of the particle concentration. The clustering statistics depend on the Stokes number, $St$, defined as the ratio of the particle friction timescale, $\tau_{\rm p} $, to the Kolmogorov timescale in the flow. In the dissipation range, the clustering intensity strongly peaks at $St \simeq 1$, and the RDF for $St \sim 1$ shows a fast power-law increase toward small scales, suggesting that turbulent clustering may considerably enhance the particle collision rate. Clustering at inertial-range scales is of particular interest to the problem of planetesimal formation. At these scales, the strongest clustering is from particles with $\tau_{\rm p}$ in the inertial range. Clustering of these particles occurs primarily around a scale where the eddy turnover time is $\sim\tau_{\rm p}$. Particles of different sizes tend to cluster at different locations, leading to flat RDFs between different particles at small scales. In the presence of multiple particle sizes, the overall clustering strength decreases as the particle size distribution broadens. We discuss particle clustering in recent models for planetesimal formation. We point out that, in the model based on turbulent clustering of chondrule-size particles, the probability of finding strong clusters that can seed planetesimals may have been significantly overestimated. We discuss various clustering mechanisms in simulations of planetesimal formation by gravitational collapse of dense clumps of meter-size particles, in particular the contribution from turbulent clustering due to the limited numerical resolution.

Stellar astrospheres--the plasma cocoons carved out of the interstellar medium by stellar winds--are continually influenced by their passage through the fluctuating interstellar medium (ISM). Inside dense interstellar regions, an astrosphere may be compressed to a size smaller than the liquid-water habitable zone distance. Habitable planets then enjoy no astrospheric buffering from the full flux of Galactic cosmic rays and interstellar dust and gas, a situation we call ``descreening.'' Recent papers (Yeghikyan and Fahr, Pavlov et al.) have suggested such global consequences as severe ozone depletion and glaciation. Using a ram-pressure balance model that includes gravitational focusing of the interstellar flow, we compute the size of the astrosphere in the apex direction as a function of parent star mass. We derive a dependence on the parent-star mass M due to gravitational focusing for densities larger than about 100 (M/M_\odot)^{-2} cm^{-3}. We calculate the interstellar densities required to descreen planets in the habitable zone of solar- and subsolar-mass stars and find a critical descreening density of roughly 600 (M/M_\odot)^{-2} cm^{-3} for the Sun's velocity relative to the local ISM. Finally, we estimate from ISM observations the frequency of descreening encounters as 1--10 Gyr^{-1} for solar-type stars and 10^2--10^9 times smaller for M stars. Given this disparity, we conclude that M star habitable-zone planets are virtually never exposed to the severe effects discussed by Yeghikyan and Fahr and Pavlov et al.

We examine the effects of turbulent intermittency on the deflagration to detonation transition (DDT) in Type Ia supernovae. The Zel'dovich mechanism for DDT requires the formation of a nearly isothermal region of mixed ash and fuel that is larger than a critical size. We primarily consider the hypothesis by Khokhlov et al. and Niemeyer and Woosley that the nearly isothermal, mixed region is produced when the flame makes the transition to the distributed regime. We use two models for the distribution of the turbulent velocity fluctuations to estimate the probability as a function of the density in the exploding white dwarf that a given region of critical size is in the distributed regime due to strong local turbulent stretching of the flame structure. We also estimate lower limits on the number of such regions as a function of density. We find that the distributed regime, and hence perhaps DDT, occurs in a local region of critical size at a density at least a factor of 2-3 larger than predicted for mean conditions that neglect intermittency. This factor brings the transition density to be much larger than the empirical value from observations in most situations. We also consider the intermittency effect on the more stringent conditions for DDT by Lisewski et al. and Woosley. We find that a turbulent velocity of $10^8$ cm/s in a region of size $10^6$ cm, required by Lisewski et al., is rare. We expect that intermittency gives a weaker effect on the Woosley model with stronger criterion. The predicted transition density from this criterion remains below $10^7$ g/cm$^3$ after accounting for intermittency using our intermittency models.

Smaller terrestrial planets (< 0.3 Earth masses) are less likely to retain the substantial atmospheres and ongoing tectonic activity probably required to support life. A key element in determining if sufficiently massive "sustainably habitable" planets can form is the availability of solid planet-forming material. We use dynamical simulations of terrestrial planet formation from planetary embryos and simple scaling arguments to explore the implications of correlations between terrestrial planet mass, disk mass, and the mass of the parent star. We assume that the protoplanetary disk mass scales with stellar mass as Mdisk ~ f Mstar^h, where f measures the relative disk mass, and 1/2 < h < 2, so that disk mass decreases with decreasing stellar mass. We consider systems without Jovian planets, based on current models and observations for M stars. We assume the mass of a planet formed in some annulus of a disk with given parameters is proportional to the disk mass in that annulus, and show with a suite of simulations of late-stage accretion that the adopted prescription is surprisingly accurate. Our results suggest that the fraction of systems with sufficient disk mass to form > 0.3 Earth mass habitable planets decreases for low-mass stars for every realistic combination of parameters. This "habitable fraction" is small for stellar masses below a mass in the interval 0.5 to 0.8 Solar masses, depending on disk parameters, an interval that excludes most M stars. Radial mixing and therefore water delivery are inefficient in lower-mass disks commonly found around low-mass stars, such that terrestrial planets in the habitable zones of most low-mass stars are likely to be small and dry.

Solar hard X-ray flares can expose astronauts on lunar and deep space extravehicular activities (EVAs) to dangerous acute biological doses. We combine calculations of radiative transfer through shielding materials with subsequent transfer through tissue to show that hazardous doses, taken as >= 0.1 Gy, should occur with a probability of about 10% per 100 hours of accumulated EVA inside current spacesuits. The rapid onset and short duration of X-ray flares and the lack of viable precursor events require strategies for quick retreat, in contrast to solar proton events, which usually take hours to deliver significant fluence and can often be anticipated by flares or other light-speed precursors. Our results contrast with the view that only particle radiation poses dangers for human space exploration. Heavy-element shields provide the most efficient protection from X-ray flares, since X-rays produce no significant secondary radiation. We calculate doses due to X-ray flares behind aluminum shields and estimate the required shield masses to accompany EVA rovers.

Stable, hydrogen-burning, M dwarf stars comprise about 75% of all stars in the Galaxy. They are extremely long-lived and because they are much smaller in mass than the Sun (between 0.5 and 0.08 MSun), their temperature and stellar luminosity are low and peaked in the red. We have re-examined what is known at present about the potential for a terrestrial planet forming within, or migrating into, the classic liquid-surface-water habitable zone close to an M dwarf star. Observations of protoplanetary disks suggest that planet-building materials are common around M dwarfs, but N-body simulations differ in their estimations of the likelihood of potentially-habitable, wet planets residing within their habitable zones, which are only ~ 1/5 to 1/50 of the width of that for a G star. Particularly in light of the claimed detection of the planets with masses as small as 5.5 and 7.5 MEarth orbiting M stars, there seems no reason to exclude the possibility of terrestrial planets. Tidally locked synchronous rotation within the narrow habitable zone doesn't necessarily lead to atmospheric collapse, and active stellar flaring may not be as much of an evolutionarily disadvantageous factor as has previously been supposed. We conclude that M dwarf stars may indeed be viable hosts for planets on which the origin and evolution of life can occur. A number of planetary processes such as cessation of geothermal activity, or thermal and non-thermal atmospheric loss processes may limit the duration of planetary habitability to periods far shorter than the extreme lifetime of the M dwarf star. Nevertheless, it makes sense to include M dwarf stars in programs that seek to find habitable worlds and evidence of life.

Motivated by an interpretation of $z \sim 3$ objects by Jimenez and Haiman (2006), we examine processes that control the fraction of primordial ($Z = 0$) gas, and so primordial stars, in high-SFR Lyman break galaxies. A primordial fraction different from 1 or 0 requires microscopic diffusion catalyzed by a velocity field with timescale comparable to the duration of star formation. The only process we found that satisfies this requirement for LBGs without fine-tuning is turbulence-enhanced mixing induced by exponential stretching and compressing of metal-rich ejecta. The time-dependence of the primordial fraction for this model is calculated. We show that conclusions for all the models discussed here are virtually independent of the IMF, including extremely top-heavy IMFs.

Putative organisms on the Martian surface would be exposed to potentially high doses of ionizing radiation during strong solar X-ray flares. We extrapolate the observed flare frequency-energy release scaling relation to releases much larger than seen so far for the sun, an assumption supported by observations of flares on other solar- and subsolar-mass main sequence stars. We calculate the surficial reprocessed X-ray spectra using a Monte Carlo code we have developed. Biological doses from indirect genome damage are calculated for each parameterized flare spectrum by integration over the X-ray opacity of water. We estimate the mean waiting time for solar flares producing a given biological dose of ionizing radiation on Mars and compare with lethal dose data for a wide range of terrestrial organisms. These timescales range from decades for significant human health risk to 0.5 Myr for D. radiodurans lethality. Such doses require total flare energies of 10^33--10^38 erg, the lower range of which has been observed for other stars. Flares are intermittent bursts, so acute lethality will only occur on the sunward hemisphere during a sufficiently energetic flare, unlike low-dose-rate, extended damage by cosmic rays. We estimate the soil and CO_2 ice columns required to provide 1/e shielding as 4--9 g cm^-2, depending on flare mean energy and atmospheric column density. Topographic altitude variations give a factor of two variation in dose for a given flare. Life in ice layers that may exist ~ 100 g cm^-2 below the surface would be well protected.

Coupled one-dimensional photochemical-climate calculations have been performed for hypothetical Earth-like planets around M dwarfs. Visible, near-infrared and thermal-infrared synthetic spectra of these planets were generated to determine which biosignature gases might be observed by a future, space-based telescope. Our star sample included two observed active M dwarfs, AD Leo and GJ 643, and three quiescent model stars. The spectral distribution of these stars in the ultraviolet generates a different photochemistry on these planets. As a result, the biogenic gases CH4, N2O, and CH3Cl have substantially longer lifetimes and higher mixing ratios than on Earth, making them potentially observable by space-based telescopes. On the active M-star planets, an ozone layer similar to Earth's was developed that resulted in a spectroscopic signature comparable to the terrestrial one. The simultaneous detection of O2 (or O3) and a reduced gas in a planet's atmosphere has been suggested as strong evidence for life. Planets circling M stars may be good locations to search for such evidence.

Peaks and lulls in the star formation rate (SFR) over the history of the Galaxy produce plateaux and declines in the present day mass function (PDMF) where the main-sequence lifetime overlaps the age and duration of the SFR variation. These PDMF features can be misinterpreted as the form of the intrinsic stellar initial mass function (IMF) if the star formation rate is assumed to be constant or slowly varying with time. This effect applies to all regions that have formed stars for longer than the age of the most massive stars, including OB associations, star complexes, and especially galactic field stars. Related problems may apply to embedded clusters. Evidence is summarized for temporal SFR variations from parsec scales to entire galaxies, all of which should contribute to inferred IMF distortions. We give examples of various star formation histories to demonstrate the types of false IMF structures that might be seen. These include short-duration bursts, stochastic histories with log-normal amplitude distributions, and oscillating histories with various periods and phases. The inferred IMF should appear steeper than the intrinsic IMF over mass ranges where the stellar lifetimes correspond to times of decreasing SFRs; shallow portions of the inferred IMF correspond to times of increasing SFRs. If field regions are populated by dispersed clusters and defined by their low current SFRs, then they should have steeper inferred IMFs than the clusters. The SFRs required to give the steep field IMFs in the LMC and SMC are determined. Structure observed in several determinations of the Milky Way field star IMF can be accounted for by a stochastic and bursty star formation history.

Low-energy cosmic rays (CRs) are confined by self-generated MHD waves in the mostly neutral ISM. We show that the CR transport equation can be expressed as a continuity equation for the CR number density involving an effective convection velocity. Assuming balance between wave growth and ion-neutral damping, this equation gives a steady-state condition n_cr ~ n_i^1/2 up to a critical density for free streaming. This relation naturally accounts for the heretofore unexplained difference in CR ionization rates derived for dense diffuse clouds (McCall et al. 2003) and dark clouds, and predicts large spatial variations in the CR heating rate and pressure.

I track the history of star count estimates of the Milky Way field star and open cluster IMFs, concentrating on the neglected mass range from 1 to 15 M${_\odot}$. The prevalent belief in a universal IMF appears to be without basis for this mass range. Two recent estimates of the field star IMF using different methods and samples give values of the average logarithmic slope $\Gamma$ between -1.7 and -2.1 in the mass range 1.1 to 4 M${_\odot}$. Two older estimates between 2 and 15 M${_\odot}$ disagree severely; the field IMF in this range is essentially unknown from star counts. Variations in $\Gamma$ among open cluster IMFs in this mass range have not decreased despite numerous detailed studies, even for studies using homogeneous data and reduction procedures and including only clusters with a significant mass range. These cluster variations \textit{might} be due to the combined effects of sampling, systematic errors, stellar evolution uncertainties, dynamical evolution, and unresolved binaries. If so, then the cluster data are consistent with a universal IMF, but are also consistent with sizeable variations. The cluster data do not allow an estimate of an average IMF or $\Gamma$ because the average depends on the choice of weighting procedure and other effects. If the spread in cluster IMFs is in excess of the effects listed above, real IMF variations must occur that do not depend much on physical conditions explored so far. The complexity of the star formation process seen in observations and simulations suggests that large realization-to-realization differences might be expected, in which case an individual cluster IMF would be in part the product of evolutionary contingency in star formation, and the function of interest is the probability distribution of IMF parameters.

(Abridged) The propagation of ionizing radiation through model atmospheres of terrestrial-like exoplanets is studied for a large range of column densities and incident photon energies using a Monte Carlo code we have developed to treat Compton scattering and photoabsorption. Incident spectra from parent star flares, supernovae, and gamma-ray bursts are modeled and compared to energetic particles in importance. We find that terrestrial-like exoplanets with atmospheres thinner than about 100 g cm^-2 transmit and reprocess a significant fraction of incident gamma-rays, producing a characteristic, flat surficial spectrum. Thick atmospheres (>~ 100 g cm^-2) efficiently block even gamma-rays, but nearly all incident energy is redistributed into diffuse UV and visible aurora-like emission, increasing the effective atmospheric transmission by many orders of magnitude. Depending on the presence of molecular UV absorbers and atmospheric thickness, up to 10% of the incident energy can reach the surface as UV reemission. For the Earth, between 2 x 10^-3 and 4 x 10^-2 of the incident flux reaches the ground in the biologically effective 200--320 nm range, depending on O_2/O_3 shielding. Finally, we suggest that transient atmospheric ionization layers can be frequently created at low altitudes. We conclude that these events can produce frequent fluctuations in atmospheric ionization levels and surficial UV fluxes on terrestrial-like planets.

In this paper, we investigate some chemokinematical properties of the Milky Way disk, by using a sample composed by 424 late-type dwarfs. We show that the velocity dispersion of a stellar group correlates with the age of this group, according to a law proportional to t^0.26, where t is the age of the stellar group. The temporal evolution of the vertex deviation is considered in detail. It is shown that the vertex deviation does not seem to depend strongly on the age of the stellar group. Previous studies in the literature seem to not have found it due to the use of statistical ages for stellar groups, rather than individual ages. The possibility to use the orbital parameters of a star to derive information about its birthplace is investigated, and we show that the mean galactocentric radius is likely to be the most reliable stellar birthplace indicator. However, this information cannot be presently used to derive radial evolutionary constraints, due to an intrinsic bias present in all samples constructed from nearby stars. An extensive discussion of the secular and stochastic heating mechanisms commonly invoked to explain the age-velocity dispersion relation is presented. We suggest that the age-velocity dispersion relation could reflect the gradual decrease in the turbulent velocity dispersion from which disk stars form, a suggestion originally made by Tinsley and Larson (1978) and supported by several more recent disk evolution calculations. A test to distinguish between the two types of models using high-redshift galaxies is proposed.

Interstellar turbulence has implications for the dispersal and mixing of the elements, cloud chemistry, cosmic ray scattering, and radio wave propagation through the ionized medium. This review discusses the observations and theory of these effects. Metallicity fluctuations are summarized, and the theory of turbulent transport of passive tracers is reviewed. Modeling methods, turbulent concentration of dust grains, and the turbulent washout of radial abundance gradients are discussed. Interstellar chemistry is affected by turbulent transport of various species between environments with different physical properties and by turbulent heating in shocks, vortical dissipation regions, and local regions of enhanced ambipolar diffusion. Cosmic rays are scattered and accelerated in turbulent magnetic waves and shocks, and they generate turbulence on the scale of their gyroradii. Radio wave scintillation is an important diagnostic for small scale turbulence in the ionized medium, giving information about the power spectrum and amplitude of fluctuations. The theory of diffraction and refraction is reviewed, as are the main observations and scintillation regions.

This two-part review summarizes interstellar turbulence and its implications. The first part begins with diagnostics and energy sources. Turbulence theory is considered in detail, including the basic fluid equations, solenoidal and compressible modes, global inviscid quadratic invariants, scaling arguments, energy transfer and cascade, velocity probability distributions, and turbulent pressure. Theories of magnetohydrodynamic turbulence, including collisionless MHD turbulence, are discussed, emphasizing various proposals for energy spectra. A large variety of numerical simulations of interstellar turbulence are reviewed. Models have reproduced many of the basic features of the observed scaling relations, predicted fast decay rates for supersonic MHD turbulence, and derived probability distribution functions for density. Thermal instabilities and thermal phases have a new interpretation in a turbulent medium. Large-scale models with various combinations of self-gravity, magnetic fields, supernovae, and star formation are beginning to resemble the observed interstellar medium in morphology and statistical properties. The role of self-gravity in turbulent gas is clarified, leading to new paradigms for the formation of star clusters, the stellar mass function, the origin of stellar rotation and binary stars, and the effects of magnetic fields.

(Abridged) We show that sizeable fractions of incident ionizing radiation from stochastic astrophysical sources can be redistributed to biologically and chemically important UV wavelengths, a significant fraction of which can reach the surface. This redistribution is mediated by secondary electrons, resulting from Compton scattering and X-ray photoabsorption, with energies low enough to excite atmospheric molecules and atoms, resulting in a rich aurora-like spectrum. We calculate the fraction of energy redistributed into biologically and chemically important wavelength regions for spectra characteristic of stellar flares and supernovae using a Monte-Carlo transport code written for this problem and then estimate the fraction of this energy that is transmitted from the atmospheric altitudes of redistribution to the surface for a few illustrative cases. Redistributed fractions are found to be of order 1%, even in the presence of an ozone shield. This result implies that planetary organisms will be subject to mutationally significant, if intermittent, fluences of UV-B and harder radiation even in the presence of a narrow-band UV shield like ozone. We also calculate the surficial transmitted fraction of ionizing radiation and redistributed ultraviolet radiation for two illustrative evolving Mars atmospheres whose initial surface pressures were 1 bar. Our results suggest that coding organisms on planets orbiting low-mass stars (and on the early Earth) may evolve very differently than on contemporary Earth, with diversity and evolutionary rate controlled by a stochastically varying mutation rate and frequent hypermutation episodes.

Combining results for the local cosmic rate and mean peak luminosity of GRBs with the cosmic history of the star formation rate, we provide estimates for the local GRB rate per unit blue luminosity in galaxies. We find a typical GRB rate per unit B luminosity of 2.4x10^-17 h_{70}^2/Lsun/yr. The corresponding mean rate in the Milky Way is 5.5x10^-7 h_{70}^2/yr. We conclude: 1) the ratio of supernova rate to isotropic equivalent GRB rate is large: more than 6000 SNIbc per GRB or 30,000 SNII per GRB. GRBs could arise in a large fraction of black hole-forming events only with collimation in the range 0.01 - 0.001 and a steep enough slope of the IMF; 2) GRBs cannot account for the majority of large HI holes observed in galaxies; 3) the probability that the solar system was exposed to a fluence large enough to melt the chondrules during the first 10^7 yr of solar system history is negligibly small; 4) Even for very opaque atmospheres, a significant fraction of the GRB energy is transmitted as UV lines due to excitation by secondary electrons. For eukaryotic-like organisms in thin atmospheres (e.g. contemporary Mars), or for UV line exposure in thick atmospheres (e.g. Earth), biologically significant events occur at a rate of about 100--500 /Gyr. The direct contribution of these "jolts" to mutational evolution may, however, be negligible because of the short duration of the GRBs. Evolutionary effects due to partial sterilizations and to longer-lived disruptions of atmospheric chemistry should be more important. (Abridged)

The star formation rate in galaxies should be related to the fraction of gas that can attain densities large enough for gravitational collapse. In galaxies with a turbulent interstellar medium, this fraction is controlled by the effective barotropic index $gamma = dlog P/dlog (rho)$ which measures the turbulent compressibility. When the cooling timescale is smaller than the dynamical timescale, gamma can be evaluated from the derivatives of cooling and heating functions, using the condition of thermal equilibrium. We present calculations of gamma for protogalaxies in which the metal abundance is so small that H_2 and HD cooling dominates. For a heating rate independent of temperature and proportional to the first power of density, the turbulent gas is relatively "hard", with $gamma >= 1$, at large densities, but moderately "soft", $gamma <= 0.8$, at densities below around $10^4 cm^(-3)$. At low temperatures the density probability distribution should fall ra pidly for densities larger than this value, which corresponds physically to the critical density at which collisional and radiative deexcitation rate s of HD are equal. The densities attained in turbulent protogalaxies thus depend on the relatively large deuterium abundance in our universe. We expect the same physical effect to occur in higher metallicity gas with different coolants. The case in which adiabatic (compressional) heating due to cloud collapse dominates is also discussed, and suggests a criterion for the maximum mass of Population III stars.

We discuss the temperature distribution in a two-dimensional, thermally unstable numerical simulation of the warm and cold gas in the Galactic disk, including the magnetic field, self-gravity, the Coriolis force, stellar energy injection and a realistic cooling function. We find that ~50% of the turbulent gas mass has temperatures in what would be the thermally unstable range if thermal instability were to be considered alone. This appears to be a consequence of there being many other forces at play than just thermal pressure. We also point out that a bimodal temperature pdf is a simple consequence of the form of the interstellar cooling function and is not necessarily a signature of discontinuous phase transitions.

According to recent models, gamma-ray bursts apparently explode in a wide variety of ambient densities ranging from ~ 10^{-3} to 30 cm^{-3}. The lowest density environments seem, at first sight, to be incompatible with bursts in or near molecular clouds or with dense stellar winds and hence with the association of gamma-ray bursts with massive stars. We argue that low ambient density regions naturally exist in areas of active star formation as the interiors of superbubbles. The evolution of the interior bubble density as a function of time for different assumptions about the evaporative or hydrodynamical mass loading of the bubble interior is discussed. We present a number of reasons why there should exist a large range of inferred afterglow ambient densities whether gamma-ray bursts arise in massive stars or some version of compact star coalescence. We predict that many gamma-ray bursts will be identified with X-ray bright regions of galaxies, corresponding to superbubbles, rather than with blue localized regions of star formation. Massive star progenitors are expected to have their own circumstellar winds. The lack of evidence for individual stellar winds associated with the progenitor stars for the cases with afterglows in especially low density environments may imply low wind densities and hence low mass loss rates combined with high velocities. If gamma-ray bursts are associated with massive stars, this combination might be expected for compact progenitors with atmospheres dominated by carbon, oxygen or heavier elements, that is, progenitors resembling Type Ic supernovae.

We estimate the frequency of intermittent hypermutation events and disruptions of planetary/satellite photochemistry due to ultraviolet radiation from core collapse supernova explosions. Calculations are presented for planetary systems in the local Milky Way, including the important moderating effects of vertical Galactic structure and UV absorption by interstellar dust. The events are particularly frequent for satellites of giant gas planets at \gtrsim 5-10 AU distance from solar-type parent stars, or in the conventional habitable zones for planets orbiting spectral type K and M parent stars, with rates of significant jolts about 10^3 - 10^4 per Gyr. The steep source spectra and existing data on UVA and longer-wavelength radiation damage in terrestrial organisms suggest that the mutational effects may operate even on planets with ozone shields. We argue that the mutation doubling dose for UV radiation should be much smaller than the mean lethal dose, using terrestrial prokaryotic organisms as our model, and that jolts may lead to important real-time evolutionary episodes if the jolt durations are longer than about a week, corresponding to several hundred generation times, or much less if the equivalent of mutator genes exist in extraterrestrial organisms. Longer-term phylogenetic effects are likely if atmospheric photochemical disturbances lead to niche creation or destruction in relevant habitats.

The effects of wind-driven star formation feedback on the spatio-temporal organization of stars and gas in galaxies is studied using two-dimensional intermediate-representational quasi-hydrodynamical simulations. The model retains only a reduced subset of the physics, including mass and momentum conservation, fully nonlinear fluid advection, inelastic macroscopic interactions, threshold star formation, and momentum forcing by winds from young star clusters on the surrounding gas. Expanding shells of swept-up gas evolve through the action of fluid advection to form a ``turbulent'' network of interacting shell fragments whose overall appearance is a web of filaments (in two dimensions). A new star cluster is formed whenever the column density through a filament exceeds a critical threshold based on the gravitational instability criterion for an expanding shell, which then generates a new expanding shell after some time delay. A filament- finding algorithm is developed to locate the potential sites of new star formation. The major result is the dominance of multiple interactions between advectively-distorted shells in controlling the gas and star morphology, gas velocity distribution and mass spectrum of high mass density peaks, and the global star formation history. The gas morphology observations of gas in the LMC and in local molecular clouds. The frequency distribution of present-to-past average global star formation rate, the distribution of gas velocities in filaments (found to be exponential), and the cloud mass spectra (estimated using a structure tree method), are discussed in detail.

Multifractal scaling (MFS) refers to structures that can be described as a collection of interwoven fractal subsets which exhibit power-law spatial scaling behavior with a range of scaling exponents (concentration, or singularity, strengths) and dimensions. The existence of MFS implies an underlying multiplicative (or hierarchical, or cascade) process. Panoramic column density images of several nearby star- forming cloud complexes, constructed from IRAS data and justified in an appendix, are shown to exhibit such multifractal scaling, which we interpret as indirect but quantitative evidence for nested hierarchical structure. The relation between the dimensions of the subsets and their concentration strengths (the "multifractal spectrum'') appears to satisfactorily order the observed regions in terms of the mixture of geometries present: strong point-like concentrations, line- like filaments or fronts, and space-filling diffuse structures. This multifractal spectrum is a global property of the regions studied, and does not rely on any operational definition of "clouds.'' The range of forms of the multifractal spectrum among the regions studied implies that the column density structures do not form a universality class, in contrast to indications for velocity and passive scalar fields in incompressible turbulence, providing another indication that the physics of highly compressible interstellar gas dynamics differs fundamentally from incompressible turbulence. (Abstract truncated)

A chromospheric age distribution of 552 late-type dwarfs is transformed into a star formation history by the application of scale height corrections, stellar evolutionary corrections and volume corrections. We show that the disk of our Galaxy has experienced enhanced episodes of star formation at 0-1 Gyr, 2-5 Gyr and 7-9 Gyr ago, although the reality of the latter burst is still uncertain. The star sample birthsites are distributed over a very large range of distances because of orbital diffusion, and so give an estimate of the global star formation rate. These results are compared with the metal-enrichment rate, given by the age-metallicity relation, with the expected epochs of close encounters between our Galaxy and the Magellanic Clouds, and with previous determinations of the star formation history. Simulations are used to examine the age-dependent smearing of the star formation history due to age uncertainties, and the broadening of the recovered features, as well as to measure the probability level that the history derived to be produced by statistical fluctuations of a constant star formation history. We show with more than 98% level of confidence that the the star formation of our Galaxy was not constant.

The age-metallicity relation of the solar neighbourhood is studied using a sample of 552 late-type dwarfs. This sample was built from the intersection of photometric catalogues with chromospheric activity surveys of the Mount Wilson group. For these stars, metallicities were estimated from uvby data, and ages were calculated from their chromospheric emission levels using a new metallicity-dependent chromospheric activity-age relation developed by Rocha-Pinto & Maciel (1998). A careful estimate of the errors in the chromospheric age is made. The errors in the chromospheric indices are shown to include partially the effects of the stellar magnetic cycles, although a detailed treatment of this error is still beyond our knowledge. It is shown that the results are not affected by the presence of unresolved binaries in the sample. We derive an age-metallicity relation which confirms the mean trend found by previous workers. The mean metallicity shows a slow, steady increase with time, amounting at least 0.56 dex in 15 Gyr. The initial metallicity of the disk is around -0.70 dex, in agreement with the G dwarf metallicity distribution. According to our data, the intrinsic cosmic dispersion in metal abundances is around 0.13 dex, a factor of two smaller than that found by Edvardsson et al. (1993). We show that chromospheric ages are compatible with isochrone ages, within the expected errors, so that the difference in the scatter cannot be caused by the accuracy of our ages and metallicities. This reinforces some suggestions that the Edvarsson et al.'s sample is not suitable to the determination of the age-metallicity relation.

We propose that galactic shocks propagating through interstellar density fluctuations provide a mechanism for the intermittent replenishment, or "pumping," of the supersonic motions and internal density enhancements observed pervasively within cool atomic and molecular interstellar structures, without necessarily requiring the presence of self-gravity, magnetic fields, or young stars. The shocks are assumed to be due to a variety of galactic sources on a range of scales. An analytic result for the kinematic vorticity generated by a shock passing through a radially-stratified two-dimensional isobaric model cloud is derived, assuming that the Mach number is not so large that the cloud is disrupted, and neglecting the shock curvature and cloud distortion. Two-dimensional lattic gas hydrodynamic simulations at modest Mach numbers were used to verify the analytic result. The induced internal velocities are initially a significant fraction of the shock speed divided by the square root of the density contrast, accounting for both the observed linewidth amplitudes and the apparent cloud-to-cloud linewidth-density scaling. The linewidth-size relation could then be interpreted in terms of the well-known power spectrum of a system of shocks. The induced vortical energy should quickly be converted to compressible and MHD modes, and so would be difficult to observe directly, even though it would still be the power source for the other modes. The shockpump thus produces density structure without the necessity of any sort of instability. We argue that the shockpump should lead to nested shock-induced structures, providing a cascade mechanism for supersonic "turbulence" and a physical explanation for the fractal-like structure of the cool interstellar medium.

The star formation rate history of the Milky Way is derived using the chromospheric age distribution for 552 stars in the solar neighborhood. The stars sample birthsites distributed over a very large range of distances because of orbital diffusion, and so give an estimate of the global star formation rate history. The derivation incorporates the metallicity dependence of chromospheric emission at a given age, and corrections to account for incompleteness, scale height-age correlations, and stellar evolutionary effects. The results show that the Milky Way has undergone fluctuations in the global star formation rate with amplitudes greater than a factor of 2-3 on timescales less than 0.2-1 Gyr. The actual history is likely to be more bursty than found here because of the smearing effect of age uncertainties. There is some evidence for a slow secular increase in the star formation rate, perhaps a record of the accumulation history of our galaxy. A smooth nearly-constant star formation rate history is strongly ruled out, confirming the result first discovered by Barry (1988) using a smaller sample and a different age calibration. This result suggests that galaxies can fluctuate coherently on large scales.

We investigate numerically the role of thermal instability (TI) as a generator of density structures in the interstellar medium (ISM), both by itself and in the context of a globally turbulent medium. Simulations of the instability alone show that the condenstion process which forms a dense phase (``clouds'') is highly dynamical, and that the boundaries of the clouds are accretion shocks, rather than static density discontinuities. The density histograms (PDFs) of these runs exhibit either bimodal shapes or a single peak at low densities plus a slope change at high densities. Final static situations may be established, but the equilibrium is very fragile: small density fluctuations in the warm phase require large variations in the density of the cold phase, probably inducing shocks into the clouds. This result suggests that such configurations are highly unlikely. Simulations including turbulent forcing show that large- scale forcing is incapable of erasing the signature of the TI in the density PDFs, but small-scale, stellar-like forcing causes erasure of the signature of the instability. However, these simulations do not reach stationary regimes, TI driving an ever-increasing star formation rate. Simulations including magnetic fields, self-gravity and the Coriolis force show no significant difference between the PDFs of stable and unstable cases, and reach stationary regimes, suggesting that the combination of the stellar forcing and the extra effective pressure provided by the magnetic field and the Coriolis force overwhelm TI as a density-structure generator in the ISM. We emphasize that a multi-modal temperature PDF is not necessarily an indication of a multi-phase medium, which must contain clearly distinct thermal equilibrium phases.

The series of papers by Chandrasekhar and Munch in the 1950s were concerned with the use of statistical models to infer the properties of interstellar clouds based on the observed spatial brightness fluctuations of the Milky Way. The present paper summarizes the subsequent influence of this work, concentrating on the departure from their earlier discrete cloud model to a continuous stochastic model in Paper V of the series. The contrast between the two models anticipated and parallels current tensions in the interpretation of interstellar structure, as well as intergalactic Lyman alpha clouds. The case of interstellar structure is discussed in some detail. Implications concerning the reification of models and the ability of scientific abstraction to model complex phenomena are also briefly discussed.

The relation between the star formation rate and the kinetic energy increase in a region containing a large number of stellar sources is investigated as a possible prescription for star formation feedback in larger scale galaxy evolution simulations, and in connection with observed scaling relations for molecular clouds, extragalactic giant H II regions, and starburst galaxies. The kinetic energy increase is not simply proportional to the source input rate, but depends on the competition between stellar power input and dissipation due to interactions between structures formed and driven by the star formation. A simple one-zone model is used to show that, in a steady-state, the energy increase should be proportional to the two-thirds power of the stellar energy injection rate, with additional factors depending on the mean density of the region and mean column density of fragments. The scaling relation is tested using two-dimensional pressureless hydrodynamic simulations of wind-driven star formation, in which star formation occurs according to a threshold condition on the column density through a shell, and a large number of shells are present at any one time. The morphology of the simulations resembles an irregular network or web of dynamically-interacting filaments. A set of 16 simulations in which different parameters were varied agree remarkably with the simple analytical prescription for the scaling relation. Converting from wind power of massive stars to Lyman continuum luminosity shows that the cluster wind model for giant H II regions may still be viable.

The probability density functions (pdfs) of molecular line centroid velocity fluctuations and fluctuation differences at different spatial lags are estimated for several nearby molecular clouds with active internal star formation. The data consist of over 75,000 $^{13}$CO line profiles divided among twelve spatially and/or kinematically distinct regions. Although three regions (all in Mon R2) appear nearly Gaussian, the others show strong evidence for non-Gaussian, often nearly exponential, centroid velocity pdfs, possibly with power law contributions in the far tails. Evidence for nearly exponential centroid pdfs in the neutral HI component of the ISM is also presented, based on older optical and radio observations. These results are in contrast to pdfs found in isotropic incompressible turbulence experiments and simulations. Furthermore, no evidence is found for the scaling of difference pdf kurtosis with Reynolds number which is seen in incompressible turbulence, and the spatial distribution of high-amplitude velocity differences shows little indication of the filamentary appearance predicted by decay simulations dominated by vortical interactions. The variation with lag of the difference pdf moments is presented as a constraint on future simulations.

Observations have not yielded convincing results concerning the form of the stellar initial mass function (IMF) or its variations in space and time, so it is proposed that theoretical models may provide useful guidance. Several classes of theoretical models for the physics controlling the IMF are reviewed, emphasizing the dependence of the results on the probability distribution of velocity and density fluctuations, the spatial density distribution, the softness of the equation of state as controlled by the dominant heating and cooling processes, the ambient radiation field, and the star formation rate itself. One particularly relevant result is the ability of turbulent interactions to form high-density condensations, which should vary with decreasing redshift as the dominant cooling mechanism changes from H$_2$, to fine structure lines, and finally to CO and gas-grain collisions, due to corresponding changes in the effective polytropic exponent. A second theme is the ability of some of these models to produce a characteristic or turnover mass in the IMF which has nothing to do with the Jeans mass. Further development of such models should provide IMFs whose functional dependence can be used as input for simulations of galaxy formation, in which the energy input from massive stars plays a crucial role, and for calculations of the consequences of protogalactic star formation on the subequent evolution of the universe. The probable irrelevance of the concepts of pressure-confined clouds and the formation of condensations by thermal instability in a turbulent interstellar medium is also discussed.

We examine the idea that diffuse and giant molecular clouds and their substructure form as density fluctuations induced by large scale interstellar turbulence. We do this by investigating the topology of various fields in realistic simulations of the ISM. We find that a) the velocity field is continuous across threshold-defined cloud boundaries; b) such cloud boundaries are rather arbitrary, with no correspondence to any actual physical boundary, such as a density discontinuity; c) abrupt velocity jumps are coincident with the density maxima; d) the volume and surface kinetic terms in the Eulerian Virial Theorem for a cloud ensemble are comparable in general; e) the magnetic field exhibits bends and reversals highly correlated with similar density features. These results suggest that clouds are formed by colliding gas streams. Within this framework, we argue that thermal pressure equilibrium is irrelevant for cloud confinement in a turbulent medium, since inertial motions can still distort or disrupt a cloud. Turbulent pressure confinement appears self-defeating, because turbulence contains large-scale motions which necessarily distort cloud boundaries. Density-weighted velocity histograms show similar FWHMs and similar multi-component structure to those of observational line profiles, though the histogram features do not correspond to isolated "clumps", but rather to extended regions throughout a cloud. We argue that the results presented here may be also applicable to small scales with larger densities (molecular clouds and cores) and suggest that quasi- hydrostatic configurations cannot be produced from turbulent fluctuations unless the thermodynamic behavior of the flow becomes nearly adiabatic at late stages of collapse. We expect this to occur only at protostellar densities.

Several recent observational studies have shown that the clustering of young stars in local star-forming regions, and of Cepheids in the LMC, can be described by a power law two-point correlation function. We show by numerical simulations that the observed range in power law slopes can be accounted for by a model in which stellar winds drive expanding shells that are subjected to nonlinear fluid advection and interactions with other shells, and in which star formation occurs when a threshold shell column density is exceeded. The models predict how the power law slope should depend on the maximum age of the stellar sample and the average star formation rate, although a number of effects preclude a comparison with currently-available data. We also show how stellar migration flattens the power law slope below a scale that depends on the velocity dispersion and age of the sample, an effect which may explain the secondary breaks in the observed correlation functions of some regions at large separations. Problems with using the correlation function as a descriptor of clustering structure for statistically inhomogeneous data sets are discussed.

A survey of results concerning the IMF derived from star counts is presented, including work up to, but not including, that presented in these proceedings. The situation regarding low-mass stars in the field and in clusters, high-mass stars and intermediate-mass stars in clusters and associations of the Milky Way and LMC, pre-main sequence objects in visible and embedded clusters, and the IMF in galaxies more distant than the Magellanic Clouds is discussed, with an emphasis on the sources of uncertainty. Most of these uncertainties, especially radial mass segregation and unresolved binaries, would steepen the true IMF relative to the apparent IMF. Several cases of apparently large variations in cluster IMFs are pointed out, and a graphical comparison of results for about 60 clusters shows a spread of at least unity in the logarithmic IMF index for all mass ranges above about 1 $M_{\sun}$. I conclude that either: 1. The uncertainties are so large that very little can be said about an average IMF or IMF variations; or 2. If the observations are taken at face value, there are strong indications of IMF variations, which do not seem to correlate with obvious environmental conditions like metallicity or stellar density. If there is an average IMF, I suggest that it is steepest at intermediate masses. If the variations are real, they offer a useful test of theoretical models.

We investigate the form of the one-point probability distribution function (pdf) for the density field of the interstellar medium using numerical simulations that successively reduce the number of physical processes included. Two-dimensional simulations of self-gravitating supersonic MHD and hydrodynamic turbulence, and of decaying Burgers turbulence, produce in all cases filamentary density structures and a power-law density pdf with logarithmic slope around -1.7. This suggests that the functional form of the pdf and the general filamentary morphology are the signature of the nonlinear advection operator. These results do not support previous claims that the pdf is lognormal. A series of 1D simulations of forced supersonic polytropic turbulence is used to resolve the discrepancy. They suggest that the pdf is lognormal only for effective polytropic indices $\gamma=1$ (or nearly lognormal for $\gamma\not=1$ if the Mach number is sufficiently small), while power laws develop at high densities if $\gamma<1$. We evaluate the polytropic index for conditions relevant to the cool interstellar medium using published cooling functions and different heating sources, finding that a lognormal pdf may occur at densities between 10$^3$ and at least 10$^4$ cm$^{-3}$. Several applications are examined. First, we question a recent derivation of the IMF from the density pdf by Padoan, Nordlund & Jones because a) the pdf does not contain spatial information, and b) their derivation produces the most massive stars in the voids of the density distribution. Second, we illustrate how a distribution of ambient densities can alter the predicted form of the size distribution of expanding shells. Finally, a brief comparison is made with the density pdfs found in cosmological simulations.

It is often assumed that galaxies cannot generate large-scale coherent star- forming activity without some organizing agent, such as spiral density waves, bars, large-scale instabilities, or external perturbations due to encounters with other galaxies. We present simulations of a simple model of star formation in which local spatial couplings lead to large-scale coherent, and even synchronized, patterns of star formation without any explicit propagation or any separate organizing agent. At a given location, star formation is assumed to occur when the gas velocity dispersion falls below a critical value dependent on the density. Young stars inject energy into the gas in their neighborhood, increasing the velocity dispersion and inhibiting the instability. A dissipation function continually "cools" the gas. The stability of this local inhibitory feedback model is examined both analytically and numerically. A large number of two-dimensional simulations are used to examine the effect of spatial couplings due to energy injection into neighboring regions. We find that several distinct types of behavior can be demarcated in a phase diagram whose parameter axes are the density (assumed constant in most models) and spatial coupling strength. These "phases" include, with decreasing density, a spatially homogeneous steady state, oscillatory "islands," traveling waves of star formation or global synchronization, and scattered "patches" of star formation activity. The coherence effects are explained in terms of the ability of the energy injected near a star formation site to introduce phase correlations in the subsequent cooling curves of neighboring regions.

The frequency distributions of sizes of ``clouds" and ``clumps" within clouds are significantly flatter for extinction surveys than for CO spectral line surveys, even for comparable size ranges. A possible explanation is the blocking of extinction clouds by larger foreground clouds (occlusion), which should not affect spectral line surveys much because clouds are resolved in velocity space along a given line of sight. We present a simple derivation of the relation between the true and occluded size distributions, assuming clouds are uniformly distributed in space or the distance to a cloud comples is much greater than the size of the complex. Because the occlusion is dominated by the largest clouds, we find that occlusion does not affect the measured size distribution except for sizes comparable to the largest size, implying that occlusion is not responsible for the discrepancy if the range in sizes of the samples is large. However, we find that the range in sizes for many of the published observed samples is actually quite small, which suggests that occlusion does affect the extinction sample and/or that the discrepancy could arise from the different operational definitions and selection effects involved in the two samples. Size and mass spectra from an IRAS survey (Wood \etal\ 1994) suggest that selection effects play a major role in all the surveys. We conclude that a reliable determination of the ``true" size and mass spectra of clouds will require spectral line surveys with very high signal-to-noise and sufficient resolution and sampling to cover a larger range of linear sizes, as well as careful attention to selection effects.

Probability density functions (pdfs) of $^{13}CO$ emission line centroid (line-of-sight, intensity-weighted average) velocities are presented for several densely sampled molecular clouds as quantitative descriptors of their underlying dynamics. Although some are approximately Gaussian in form, most of the pdfs exhibit relatively broader, often nearly exponential, tails, similar to the pdfs of velocity {\em differences} and {\em derivatives} (but not the velocity field itself) found in experiments and numerical simulations of incompressible turbulence. The broad pdf tails found in the present work are also similar to those found in decades-old measurements of interstellar velocity pdfs using atomic line centroids, and to the excess wing emission recently found in individual molecular line profiles. Some possible interpretations of the observed deviations are briefly discussed, although none of these account for the nearly exponential tails.