External cosmic rays impinging upon circumstellar disks associated with young stellar objects provide an important source of ionization, and as such, play an important role in disk evolution and planet formation. However, these incoming cosmic rays are affected by a variety of physical processes internal to stellar/disk systems, including modulation by turbulent magnetic fields. Globally, these fields naturally provide both a funneling effect, where cosmic rays from larger volumes are focused into the disk region, and a magnetic mirroring effect, where cosmic rays are repelled due to the increasing field strength. This paper considers cosmic ray propagation in the presence of a turbulent spiral magnetic field, analogous to that produced by the Solar wind. The interaction of this wind with the interstellar medium defines a transition radius, analogous to the Heliopause, which provides the outer boundary to this problem. We construct a new coordinate system where one coordinate follows the spiral magnetic field lines and consider magnetic perturbations to the field in the perpendicular directions. The presence of magnetic turbulence replaces the mirroring points with a distribution of values and moves the mean location outward. Our results thus help quantify the degree to which cosmic ray fluxes are reduced in circumstellar disks by the presence of magnetic field structures that are shaped by stellar winds. The new coordinate system constructed herein should also be useful in other astronomical applications.

Quasars at high redshift provide direct information on the mass growth of supermassive black holes and, in turn, yield important clues about how the Universe evolved since the first (Pop III) stars started forming. Yet even basic questions regarding the seeds of these objects and their growth mechanism remain unanswered. The anticipated launch of eROSITA and ATHENA is expected to facilitate observations of high-redshift quasars needed to resolve these issues. In this paper, we compare accretion-based supermassive black hole growth in the concordance LCDM model with that in the alternative Friedmann-Robertson Walker cosmology known as the R_h=ct universe. Previous work has shown that the timeline predicted by the latter can account for the origin and growth of the > 10^9 M_sol highest redshift quasars better than that of the standard model. Here, we significantly advance this comparison by determining the soft X-ray flux that would be observed for Eddington-limited accretion growth as a function of redshift in both cosmologies. Our results indicate that a clear difference emerges between the two in terms of the number of detectable quasars at redshift z > 6, raising the expectation that the next decade will provide the observational data needed to discriminate between these two models based on the number of detected high-redshift quasar progenitors. For example, while the upcoming ATHENA mission is expected to detect ~0.16 (i.e., essentially zero) quasars at z ~ 7 in R_h=ct, it should detect ~160 in LCDM---a quantitatively compelling difference.

The measured properties of the epoch of reionization (EoR) show that reionization probably began around z ~ 12-15 and ended by z=6. In addition, a careful analysis of the fluctuations in the cosmic microwave background indicate a scattering optical depth tau ~ 0.066+/-0.012 through the EoR. In the context of LCDM, galaxies at intermediate redshifts and dwarf galaxies at higher redshifts now appear to be the principal sources of UV ionizing radiation, but only for an inferred (ionizing) escape fraction f_ion ~ 0.2, which is in tension with other observations that suggest a value as small as ~ 0.05. In this paper, we examine how reionization might have progressed in the alternative Friedmann-Robertson Walker cosmology known as the R_h=ct Universe, and determine the value of f_ion required with this different rate of expansion. We find that R_h=ct accounts quite well for the currently known properties of the EoR, as long as its fractional baryon density falls within the reasonable range 0.026 < Omega_b < 0.037. This model can also fit the EoR data with f_ion ~ 0.05, but only if the Lyman continuum photon production is highly efficient and Omega_b ~ 0.037. These results are still preliminary, however, given their reliance on a particular form of the star-formation rate density, which is still uncertain at very high redshifts. It will also be helpful to reconsider the EoR in R_h=ct when complete structure formation models become available.

Radioactive nuclei play an important role in planetary evolution by providing an internal heat source, which affects planetary structure and helps facilitate plate tectonics. A minimum level of nuclear activity is thought to be necessary --- but not sufficient --- for planets to be habitable. Extending previous work that focused on short-lived nuclei, this paper considers the delivery of long-lived radioactive nuclei to circumstellar disks in star forming regions. Although the long-lived nuclear species are always present, their abundances can be enhanced through multiple mechanisms. Most stars form in embedded cluster environments, so that disks can be enriched directly by intercepting ejecta from supernovae within the birth clusters. In addition, molecular clouds often provide multiple episodes of star formation, so that nuclear abundances can accumulate within the cloud; subsequent generations of stars can thus receive elevated levels of radioactive nuclei through this distributed enrichment scenario. This paper calculates the distribution of additional enrichment for $^{40}$K, the most abundant of the long-lived radioactive nuclei. We find that distributed enrichment is more effective than direct enrichment. For the latter mechanism, ideal conditions lead to about 1 in 200 solar systems being directly enriched in $^{40}$K at the level inferred for the early solar nebula (thereby doubling the abundance). For distributed enrichment from adjacent clusters, about 1 in 80 solar systems are enriched at the same level. Distributed enrichment over the entire molecular cloud is more uncertain, but can be even more effective.

Most star formation in the Galaxy takes place in clusters, where the most massive members can affect the properties of other constituent solar systems. This paper considers how clusters influence star formation and forming planetary systems through nuclear enrichment from supernova explosions, where massive stars deliver short-lived radioactive nuclei (SLRs) to their local environment. The decay of these nuclei leads to both heating and ionization, and thereby affects disk evolution, disk chemistry, and the accompanying process of planet formation. Nuclear enrichment can take place on two spatial scales: [1] Within the cluster itself ($\ell\sim1$pc), the SLRs are delivered to the circumstellar disks associated with other cluster members. [2] On the next larger scale ($\ell\sim2-10$pc), SLRs are injected into the background molecular cloud; these nuclei provide heating and ionization to nearby star-forming regions, and to the next generation of disks. For the first scenario, we construct the expected distributions of radioactive enrichment levels provided by embedded clusters. Clusters can account for the SLR mass fractions inferred for the early Solar Nebula, but typical SLR abundances are lower by a factor of $\sim10$. For the second scenario, we find that distributed enrichment of SLRs in molecular clouds leads to comparable abundances. For both the direct and distributed enrichment processes, the masses of $^{26}$Al and $^{60}$Fe delivered to individual circumstellar disks typically fall in the range $10-100pM_\odot$ (where $1pM_\odot=10^{-12}M_\odot$). The corresponding ionization rate due to SLRs typically falls in the range $\zeta_{SLR}\sim1-5\times10^{-19}$ sec$^{-1}$. This ionization rate is smaller than that due to cosmic rays, $\zeta_{CR}\sim10^{-17}$ sec$^{-1}$, but will be important in regions where cosmic rays are attenuated (e.g., disk mid-planes).

The magnetic fields associated with young stellar objects are expected to have an hour-glass geometry, i.e., the magnetic field lines are pinched as they thread the equatorial plane surrounding the forming star but merge smoothly onto a background field at large distances. With this field configuration, incoming cosmic rays experience both a funneling effect that acts to enhance the flux impinging on the circumstellar disk and a magnetic mirroring effect that acts to reduce that flux. To leading order, these effects nearly cancel out for simple underlying magnetic field structures. However, the environments surrounding young stellar objects are expected to be highly turbulent. This paper shows how the presence of magnetic field fluctuations affects the process of magnetic mirroring, and thereby changes the flux of cosmic rays striking circumstellar disks. Turbulence has two principle effects: 1) The (single) location of the magnetic mirror point found in the absence of turbulence is replaced with a wide distribution of values. 2) The median of the mirror point distribution moves outward for sufficiently large fluctuation amplitudes (roughly when $\delta B/B_0>0.2$ at the location of the turbulence-free mirror point); the distribution becomes significantly non-gaussian in this regime as well. These results may have significant consequences for the ionization fraction of the disk, which in turn dictates the efficiency with which disk material can accrete onto the central object. A similar reduction in cosmic ray flux can occur during the earlier protostellar stages; the decrease in ionization can help alleviate the magnetic braking problem that inhibits disk formation.

How and where cosmic rays are produced, and how they diffuse through various turbulent media, represent fundamental problems in astrophysics with far reaching implications, both in terms of our theoretical understanding of high-energy processes in the Milky Way and beyond, and the successful interpretation of space-based and ground based GeV and TeV observations. For example, recent and ongoing detections, e.g., by Fermi (in space) and HESS (in Namibia), of $\gamma$-rays produced in regions of dense molecular gas hold important clues for both processes. In this paper, we carry out a comprehensive numerical investigation of relativistic particle acceleration and transport through turbulent magnetized environments in order to derive broadly useful scaling laws for the energy diffusion coefficients.

Together, the Fermi-LAT and HESS have revealed the presence of an unusual GeV-TeV source coincident with Sgr A* at the Galactic center. Its high-energy emission appears to be bimodal, hinting at an energizing process more sophisticated than mere shock acceleration. It has been suggested that this may be evidence of strong, rapid variability, as required if Sgr A*'s emission were responsible for the fluorescent X-ray echos detected in nearby molecular clouds. In this {\it Letter}, however, we show that stochastic acceleration in a more realistic two-phase environment surrounding the central black hole can accommodate Sgr A*'s high-energy spectrum quite well. We therefore suggest that the Fermi-HESS data alone do not necessarily provide evidence for strong variability in Sgr A*.

Most star formation in our galaxy occurs within embedded clusters, and these background environments can affect the star and planet formation processes occurring within them. In turn, young stellar members can shape the background environment and thereby provide a feedback mechanism. This work explores one aspect of stellar feedback by quantifying the background X-ray radiation fields produced by young stellar objects. Specifically, the distributions of X-ray luminosities and X-ray fluxes produced by cluster environments are constructed as a function of cluster membership size $N$. Composite flux distributions, for given distributions of cluster sizes $N$, are also constructed. The resulting distributions are wide and the X-ray radiation fields are moderately intense, with the expected flux levels exceeding the cosmic and galactic X-ray backgrounds by factors of $\sim10-1000$ (for energies 0.2 -- 15 keV). For circumstellar disks that are geometrically thin and optically thick, the X-ray flux from the background cluster dominates that provided by a typical central star in the outer disk where $r \ga 9 - 14$ AU. In addition, the expectation value of the ionization rate provided by the cluster X-ray background is $\zeta_X\sim8\times10^{-17}$ s$^{-1}$, about 4 -- 8 times larger than the canonical value of the ionization rate from cosmic rays. These elevated flux levels in clusters indicate that X-rays can affect ionization, chemistry, and heating in circumstellar disks and in the material between young stellar objects.

The presence of relativistic particles at the center of our galaxy is evidenced by the diffuse TeV emission detected from the inner $\sim$$2^\circ$ of the Galaxy. Although it is not yet entirely clear whether the origin of the TeV photons is due to hadronic or leptonic interactions, the tight correlation of the intensity distribution with the distribution of molecular gas along the Galactic ridge strongly points to a pionic-decay process involving relativistic protons. In earlier work, we concluded that point-source candidates, such as the supermassive black hole Sagittarius A* (identified with the HESS source J1745-290), or the pulsar wind nebulae dispersed along the Galactic plane, could not account for the observed diffuse TeV emission from this region. Motivated by this result, we consider here the feasibility that the cosmic rays populating the Galactic Center (GC) region are accelerated in situ by magnetic turbulence. Our results indicate that even in a highly conductive environment, this mechanism is efficient enough to energize protons within the intercloud medium to the $\ga $ TeV energies required to produce the HESS emission.

The diffuse TeV emission detected from the inner $\sim2^\circ$ of the Galaxy appears to be strongly correlated with the distribution of molecular gas along the Galactic ridge. Although it is not yet entirely clear whether the origin of the TeV photons is due to hadronic or leptonic interactions, the tight correlation of the intensity distribution with the molecular gas strongly points to a pionic-decay process involving relativistic protons. But the spectrum of the TeV radiation---a power law with index $\alpha\approx -2.3$---cannot be accommodated easily with the much steeper distribution of cosmic rays seen at Earth. In earlier work, we examined the possible sources of these relativistic protons and concluded that neither the supermassive black hole Sagittarius A* (identified with the HESS source J1745-290), nor several pulsar wind nebulae dispersed along the Galactic plane, could produce a TeV emission profile morphologically similar to that seen by HESS. We concluded from this earlier study that only relativistic protons accelerated throughout the inter-cloud medium could account for the observed diffuse TeV emission from this region. In this paper, we develop a model for diffusive proton acceleration driven by a turbulent Alfv\'enic magnetic field present throughout the gaseous medium. Though circumstantial, this appears to be the first evidence that at least some cosmic rays are accelerated diffusively within the inner $\sim300$ pc of the Galaxy.

The propagation of high-energy cosmic rays through giant molecular clouds constitutes a fundamental process in astronomy and astrophysics. The diffusion of cosmic-rays through these magnetically turbulent environments is often studied through the use of energy-dependent diffusion coefficients, although these are not always well motivated theoretically. Now, however, it is feasible to perform detailed numerical simulations of the diffusion process computationally. While the general problem depends upon both the field structure and particle energy, the analysis may be greatly simplified by dimensionless analysis. That is, for a specified purely turbulent field, the analysis depends almost exclusively on a single parameter -- the ratio of the maximum wavelength of the turbulent field cells to the particle gyration radius. For turbulent magnetic fluctuations superimposed over an underlying uniform magnetic field, particle diffusion depends on a second dimensionless parameter that characterizes the ratio of the turbulent to uniform magnetic field energy densities. We consider both of these possibilities and parametrize our results to provide simple quantitative expressions that suitably characterize the diffusion process within molecular cloud environments. Doing so, we find that the simple scaling laws often invoked by the high-energy astrophysics community to model cosmic-ray diffusion through such regions appear to be fairly robust for the case of a uniform magnetic field with a strong turbulent component, but are only valid up to $\sim 50$ TeV particle energies for a purely turbulent field. These results have important consequences for the analysis of cosmic-ray processes based on TeV emission spectra associated with dense molecular clouds.

Molecular clouds at the Galactic center (GC) have environments considerably different from their disk counterparts. The GC may therefore provide important clues about how the environment affects star formation. Interestingly, while the inner 50 parsecs of our Galaxy include a remarkable population of high-mass stars, the initial mass function (IMF) appears to be consistent with a Salpeter slope down to ~ 1 solar mass. We show here that the loss of turbulent pressure due to ambipolar diffusion and the damping of Alfven and fast MHD waves can lead to the formation of dense condensations exceeding their Jeans limit. The fragmentation and subsequent collapse of these condensations is similar to the diffusion-driven protostellar collapse mechanism expected to occur within nearby "regular" molecular clouds. As such, a Salpeter IMF at the GC is not surprising, though the short dynamical timescales associated with the GC molecular clouds may help explain the lower star formation efficiency observed from this region.

A self-similar formalism for the study of the gravitational collapse of molecular gas provides an important theoretical framework from which to explore the dynamics of star formation. Motivated by the presence of elongated and filamentary structures observed in giant molecular clouds, we build upon the existing body of work on cylindrical self-similar collapse flows by including dynamic equations of state that are different from the effective equation of state that produces the initial density distribution. We focus primarily on the collapse of initial states for which the gas is at rest and everywhere overdense from its corresponding hydrostatic equilibrium profile by a factor $\Lambda$, and apply our results toward the analysis of star formation within dense, elongated molecular cores. An important aspect of this work is the determination of the mass infall rates over a range of the parameters which define the overall state of the gas -- the overdensity parameter $\Lambda$, the index $\Gamma$ of the static equation of state, and the index $\gamma$ of the dynamic equation of state. While most of the parameter space explored in this work leads to solutions for which the underlying equations do not become singular, we do include a discussion on how to treat cases for which solutions pass smoothly through the singular surface. In addition, we also present a different class of collapse solutions for the special case $\gamma = 1$.

The High-Energy Stereoscopic System (HESS) has detected intense diffuse TeV emission correlated with the distribution of molecular gas along the galactic ridge at the centre of our Galaxy. Earlier HESS observations of this region had already revealed the presence of several point sources at these energies, one of them (HESS J1745-290) coincident with the supermassive black hole Sagittarius A*. It is still not entirely clear what the origin of the TeV emission is, nor even whether it is due to hadronic or leptonic interactions. It is reasonable to suppose, however, that at least for the diffuse emission, the tight correlation of the intensity distribution with the molecular gas indicates a pionic-decay process involving relativistic protons. In this paper, we explore the possible source(s) of energetic hadrons at the galactic centre, and their propagation through a turbulent medium. We conclude that though Sagittarius A* itself may be the source of cosmic rays producing the emission in HESS J1745-290, it cannot be responsible for the diffuse emission farther out. A distribution of point sources, such as pulsar wind nebulae dispersed along the galactic plane, similarly do not produce a TeV emission profile consistent with the HESS map. We conclude that only a relativistic proton distribution accelerated throughout the inter-cloud medium can account for the TeV emission profile measured with HESS.

A large fraction of stars form within young embedded clusters, and these environments produce a substantial ultraviolet (UV) background radiation field, which can provide feedback on the star formation process. To assess the possible effects of young stellar clusters on the formation of their constituent stars and planets, this paper constructs the expected radiation fields produced by these clusters. We include both the observed distribution of cluster sizes $N$ in the solar neighborhood and an extended distribution that includes clusters with larger $N$. The paper presents distributions of the FUV and EUV luminosities for clusters with given stellar membership $N$, distributions of FUV and EUV luminosity convolved over the expected distribution of cluster sizes $N$, and the corresponding distributions of FUV and EUV fluxes. These flux distributions are calculated both with and without the effects of extinction. Finally, we consider the effects of variations in the stellar initial mass function on these radiation fields. Taken together, these results specify the distributions of radiation environments that forming solar systems are expected to experience.

Molecular clouds interacting with supernova remnants may be subject to a greatly enhanced irradiation by cosmic rays produced at the shocked interface between the ejecta and the molecular gas. Over the past decade, broad-band observations have provided important clues about these relativistic particles and indicate that they may dominate over the locally observed cosmic-ray population by a significant amount. In this paper, we estimate the enhancement and find that the cosmic ray energy density can be up to $\sim$1000 times larger in the molecular cloud than in the field. This enhancement can last for a few Myr and leads to a corresponding increase in the ionization fraction, which has important consequences for star formation. Ionization fractions in] molecular cloud cores determine, in part, the rate of ambipolar diffusion, an important process in core formation and pre-collapse evolution. Ionization fractions in newly formed circumstellar disks affect the magneto-rotational instability mechanism, which in turn affects the rate of disk accretion. As estimated here, the increased ionization acts to increase the ambipolar diffusion time by a factor of $\sim30$ and thereby suppresses star formation. In contrast, the increased ionization fraction reduces the sizes of dead zones in accretion disks (by up to an order of magnitude) and thus increases disk accretion rates (by a comparable factor).

Stochastic acceleration of electrons interacting resonantly with a turbulent magnetic field in a small accretion torus appears to be the likely mechanism responsible for much of Sagittarius A*'s millimeter and shorter wavelength spectrum. The longer wavelength radiation is produced at larger radii by electrons either diffusing from smaller scales or accelerated in situ. An important prediction of this model is the ejection of a significant flux of relativistic protons from a magnetic-field-dominated acceleration site into the wind-shocked medium surrounding the black hole. Recently, several air Cerenkov telescopes, notably HESS, have detected TeV emission from the Galactic center, with characteristics hinting at a p-p-induced pion decay process for the \gamma-ray emission. Given (1) the size of this acceleration region measured in the radio band and (2) the wind-injected ISM mapped with Chandra using the diffuse X-rays, it is feasible to test the idea that protons accelerated within \~20 Schwarzschild radii of the black hole produce the TeV emission farther out. We show a fraction of TeV protons scattering about once within ~3 pc of Sagittarius A* and the proton power (~10^37 erg s^-1) produced in concert with the 7 mm radio emission matches the TeV luminosity well. This model explains why the TeV source does not vary on a timescale of a year or less. The particle cascade generated by the p-p scatterings also produces bremsstrahlung, inverse Compton, and synchrotron emission at longer wavelengths from secondary particles. We compare these with current measurements and demonstrate that GLAST will detect this source during its one-year all-sky survey.

This paper explores the stability of an Earth-like planet orbiting a solar-mass star in the presence of a stellar companion using ~ 400,000 numerical integrations. Given the chaotic nature of the systems being considered, we perform a statistical analysis of the ensuing dynamics for ~500 orbital configurations defined by the following set of orbital parameters: the companion mass; the companion eccentricity; the companion periastron; and the planet's inclination angle relative to the stellar binary plane. Specifically, we generate a large sample of survival times for each orbital configuration through the numerical integration of N >> 1 equivalent experiments (e.g., with the same orbital parameters but randomly selected initial orbital phases). We then construct distributions of survival time using the variable mu_s = log tau_s (where tau_s is in years) for each orbital configuration. The primary objective of this work is twofold. First, we use the mean of the distributions to gain a better understanding of what orbital configurations, while unstable, have sufficiently long survival times to make them interesting to the study of planet habitability. Second, we calculate the width, skew, and kurtosis of each mu_s distribution and look for general features that may aid further understanding and numerical exploration of these chaotic systems.

In recent years, apparent anisotropies in the ~EeV cosmic ray (CR) flux arriving at Earth from the general direction of the galactic center have been reported from the analysis of AGASA and SUGAR data. The more recently commissioned Auger Observatory has not confirmed these results. HESS has now detected an unusually soft TeV source roughly coincident with the location of the previously claimed CR anisotropy. In this paper, we develop a model for the TeV emission from this object, consistent with observations at other wavelengths, and examine the circumstances under which it might have contributed to the $\sim$ EeV cosmic ray spectrum. We find that the supernova remnant G8.7-0.1 can plausibly account for all the known radiative characteristics of HESS J1804-216, but that it can accelerate cosmic rays only up to an energy $\sim 10^5$ GeV. On the other hand, the pulsar (PSR J1803-2137) embedded within this remnant can in principle inject EeV protons into the surrounding medium, but it cannot account for the broadband spectrum of HESS J1804-216. We therefore conclude that although G8.7-0.1 is probably the source of TeV photons originating from this direction, there is no compelling theoretical motivation for expecting a cosmic ray anisotropy at this location. However, if G8.7-0.1 is indeed correctly identified with HESS J1804-216, it should also produce a $\sim$ GeV flux detectable in a one-year all sky survey by GLAST.

This paper studies the dynamical evolution of young stellar clusters with $N$ = 100 - 1000 members. We use N-body simulations to explore how evolution depends on system size $N$ and the initial conditions. Motivated by recent observations of extremely young systems, this study compares subvirial and virial starting states. Multiple realizations of equivalent cases (100 simulations per case) are used to build up a robust statistical description of these systems, e.g., distributions of closest approaches, mass profiles, and distributions of radial locations. These results provide a framework from which to assess the effects of these clusters on star and planet formation. The distributions of radial positions are used in conjunction with distributions of FUV luminosities (also calculated here) to determine the radiation exposure of circumstellar disks. The distributions of closest approaches are used in conjunction with scattering cross sections (calculated here from $10^5$ scattering experiments) to determine the probability of solar system disruption. We also use the nearby cluster NGC 1333 as a test case. Our main conclusion is that clusters in this size range have only a modest effect on forming planetary systems: Interaction rates are low so that the typical solar system experiences a single encounter within 1000 AU. Radiation exposure is low, with median FUV flux $G_0$ = 900, so that photoevaporation of disks is only important beyond 30 AU. Given the low interaction rates and modest radiation levels, we suggest that solar system disruption is a rare event in these clusters.

Motivated by recent observations that show starless molecular cloud cores exhibit subsonic inward velocities, we revisit the collapse problem for polytropic gaseous spheres. In particular, we provide a generalized treatment of protostellar collapse in the spherical limit and find semi-analytic (self-similar) solutions, corresponding numerical solutions, and purely analytic calculations of the mass infall rates (the three approaches are in good agreement). This study focuses on collapse solutions that exhibit nonzero inward velocities at large radii, as observed in molecular cloud cores, and extends previous work in four ways: (1) The initial conditions allow nonzero initial inward velocity. (2) The starting states can exceed the density of hydrostatic equilibrium so that the collapse itself can provide the observed inward motions. (3) We consider different equations of state, especially those that are softer than isothermal. (4) We consider dynamic equations of state that are different from the effective equation of state that produces the initial density distribution. This work determines the infall rates over a wide range of parameter space, as characterized by four variables: the initial inward velocity $\vin$, the overdensity $\overdense$ of the initial state, the index $\Gamma$ of the static equation of state, and the index $\gamma$ of the dynamic equation of state. For the range of parameter space applicable to observed cores, the resulting infall rate is about a factor of two larger than found in previous theoretical studies (those with hydrostatic initial conditions and $\vin = 0$).

The EGRET supernova remnants (SNR's) are all expanding into nearby dense molecular clouds, powering a shock at the interface where protons and electrons accelerate to relativistic energies. A viable mechanism for the emission of gamma$-rays in these sources is the decay of neutral pions created in collisions between the relativistic hadrons and protons in the ambient medium. But neutral pion decay alone cannot reproduce the whole high-energy spectrum, particularly below 100 MeV. A pion-decay scenario thus requires a lepton population to "fill in" the lower part of the spectrum via bremsstrahlung emission. This population, however, is constrained by the SNR radio spectrum. Taking our cue from the behavior of Sgr A East, an energetic EGRET SNR at the Galactic center, we here examine the role played in these sources by secondary leptons--electrons and positrons produced in proton-proton scattering events and the ensuing particle cascades. We show that while secondary leptons cannot account for the gamma-rays below 100 MeV, they can account for the hard radio spectra observed from the EGRET SNR's. Thus, it appears that both primary and secondary leptons may be important contributors to the overall broadband emission from these sources, but if so, must radiate most of their energy in different parts of the SNR-cloud environment. We show that shock acceleration in dense cores being overtaken by the expanding SNR shell can naturally lead to such a scenario.

Recent analysis of data sets from two extensive air shower cosmic ray detectors shows tantalizing evidence of an anisotropic overabundance of cosmic rays towards the Galactic Center (GC) that ``turns on'' around $10^{18}$ eV. We demonstrate that the anisotropy could be due to neutrons created at the Galactic Center through charge-exchange in proton-proton collisions, where the incident, high energy protons obey an $\sim E^{-2}$ power law associated with acceleration at a strong shock. We show that the normalization supplied by the gamma-ray signal from EGRET GC source 3EG J1746-2851 -- ascribed to pp induced neutral pion decay at GeV energies -- together with a very reasonable spectral index of 2.2, predicts a neutron flux at $\sim 10^{18}$ eV fully consistent with the extremely high energy cosmic ray data. Likewise, the normalization supplied by the very recent GC data from the HESS air-Cerenkov telescope at \~TeV energies is almost equally-well compatible with the $\sim 10^{18}$ eV cosmic ray data. Interestingly, however, the EGRET and HESS data appear to be themselves incompatible. We consider the implications of this discrepancy. We discuss why the Galactic Center environment can allow diffusive shock acceleration at strong shocks up to energies approaching the ankle in the cosmic ray spectrum. Finally, we argue that the shock acceleration may be occuring in the shell of Sagittarius A East, an unusual supernova remnant located very close to the Galactic Center. If this connection between the anisotropy and Sagittarius A East could be firmly established it would be the first direct evidence for a particular Galactic source of cosmic rays up to energies near the ankle.

This paper explores the stability of an Earth-like planet orbiting a solar mass star in the presence of an outer-lying intermediate mass companion. The overall goal is to estimate the fraction of binary systems that allow Earth-like planets to remain stable over long time scales. We numerically determine the planet's ejection time $\tauej$ over a range of companion masses ($M_C$ = 0.001 -- 0.5 $M_\odot$), orbital eccentricities $\epsilon$, and semi-major axes $a$. This suite of $\sim40,000$ numerical experiments suggests that the most important variables are the companion's mass $M_C$ and periastron distance $\rmin$ = $a(1-\epsilon)$ to the primary star. At fixed $M_C$, the ejection time is a steeply increasing function of $\rmin$ over the range of parameter space considered here (although the ejection time has a distribution of values for a given $\rmin$). Most of the integration times are limited to 10 Myr, but a small set of integrations extend to 500 Myr. For each companion mass, we find fitting formulae that approximate the mean ejection time as a function of $\rmin$. These functions can then be extrapolated to longer time scales. By combining the numerically determined ejection times with the observed distributions of orbital parameters for binary systems, we estimate that (at least) 50 percent of binaries allow an Earth-like planet to remain stable over the 4.6 Gyr age of our solar system.

Sgr A East appears to be a single, mixed morphology 10,000 year-old supernova remnant (SNR) at the Galactic Center. It also appears to belong to a class of remnants that have been observed and detected at 1720 MHz, the transition frequency of OH maser emission. However, if the EGRET source 3EG J1746-2852 conincident with the Galactic center is itself associated with this object, it would endow it with a gamma-ray luminosity almost two orders of magnitude greater than that of the other EGRET-detected SNR's. We here reconsider the viability of a pion-production mechanism as the source of the broadband emission observed from Sgr A East, and show that what connects these objects -- and ultimately also accounts for their different gamma-ray emissivity -- is the very important interaction between the expanding SNR shell and the surrounding molecular cloud environment. The singularly high gamma-ray luminosity of Sgr A East, as well as its unusually steep radio spectral index, can thereby be attributed to the high-density, strong magnetized environment in which it is located.

PAPER WITHDRAWN. The recent detection of Sgr A* in the X-ray band, together with the radio polarization measurements conducted over the past few years, offer the best constraints yet for understanding the nature of the emitting gas within several Schwarzschild radii ($r_S$) of this supermassive black hole candidate at the Galactic Center. Earlier, we showed that the sub-mm radiation from this source may be associated with thermal synchrotron emission from an inner Keplerian region within the circularization radius of the accreting plasma. In this paper, we extend this analysis in a very important way, by calculating the implied high-energy emission of Sgr A* associated with the orbiting, hot, magnetized gas. We find that for the accretion rate inferred from the fits to the sub-mm data, the dominant contribution to Sgr A*'s X-ray flux is due to self-Comptonization of the radio photons, rather than from bremsstrahlung. The latter is a two-body process, which would produce significant X-ray emission only at much higher accretion rates. This picture leads to the testable prediction that the physical conditions within the inner $\sim5r_S$ are variable on a time scale of order a year. In particular, the accretion rate $\dot M$ appears to have changed by about 15% between the sub-mm measurements in 1996 and 1999. Given that the radio and self-Comptonized fluxes are strongly correlated in this picture, the upcoming second generation Chandra observations of Sgr A* may provide the direct evidence required to test this model.

Previous treatments of ambipolar diffusion in star-forming molecular clouds do not consider the effects of fluctuations in the fluid fields about their mean values. This paper generalizes the ambipolar diffusion problem in molecular cloud layers to include such fluctuations. Because magnetic diffusion is a nonlinear process, fluctuations can lead to an enhancement of the ambipolar diffusion rate. In addition, the stochastic nature of the process makes the ambipolar diffusion time take on a distribution of different values. In this paper, we focus on the case of long wavelength fluctuations and find that the rate of ambipolar diffusion increases by a significant factor $\Lambda \sim 1 - 10$. The corresponding decrease in the magnetic diffusion time helps make ambipolar diffusion more consistent with observations.

In this paper, we explore the possibility that the radiative properties of the most compact region in NGC 6251* may be understood in the same sense as Sgr A*, though with some telling differences that may hint at the nature of jet formation. We show that observations of this object with ASCA, ROSAT, HST and VLBI together may be hinting at a picture in which Bondi-Hoyle accretion from an ambient ionized medium feeds a standard disk accreting at ~ 4.0*10^{22} g s^{-1}. Somewhere near the event horizon, this plasma is heated to >10^{11} K, where it radiates via thermal synchrotron (producing a radio component) and self-Comptonization (accounting for a nonthermal X-ray flux). This temperature is much greater than its virial value and the hot cloud expands at roughly the sound speed (~0.1c), after which it begins to accelerate on a parsec scale to relativistic velocities. In earlier work, the emission from the extended jet has been modeled successfully using nonthermal synchrotron self-Compton processes, with a self-absorbed inner core. In the picture we are developing here, the initial ejection of matter is associated with a self-absorbed thermal radio component that dominates the core emission on the smallest scales. The nonthermal particle distributions responsible for the emission in the extended jet are then presumably energized, e.g., via shock acceleration, within the expanding, hot gas. The power associated with this plasma represents an accretion efficiency of about 0.54, requiring dissipation in a prograde disk around a rapidly spinning black hole (with spin parameter a~1).

Maps of the Galactic electron-positron annihilation radiation show evidence for three distinct and significant features: (1) a central bulge source, (2) emission in the Galactic plane, and (3) an enhancement of emission at positive latitudes above the Galactic Center. In this paper, we explore the possibility that Sgr A East, a very prominent radio structure surrounding the Galactic nucleus, may be a significant contributer to the central bulge feature. The motivation for doing so stems from a recently proposed link between this radio object and the EGRET gamma-ray source 2EG J1746-2852. If this association is correct, then Sgr A East is also expected to be a source of copious positron production. The results presented here show that indeed Sgr A East must have produced a numerically significant population of positrons, but also that most of them have not yet had sufficient time to thermalize and annihilate. As such, Sgr A East by itself does not appear to be the dominant current source of annihilation radiation, but it will be when the positrons have cooled sufficiently and they have become thermalized. This raises the interesting possibility that the bulge component may be due to the relics of earlier explosive events like the one that produced Sgr A East.

Sgr A East is a very prominent elongated shell structure surrounding (though off-centered from) the Galactic nucleus. Its energetics (~4x10^52 ergs), based on the power required to carve out the radio synchrotron remnant within the surrounding dense molecular cloud, appear to be extreme compared to the total energy (~10^51 ergs) released in a typical supernova (SN) explosion. Yet it shares several characteristics in common with SN remnants (SNRs), the most significant of which is the ~0.1-10 GeV gamma-ray spectrum measured by EGRET, if we associate the Galactic center source 2EGJ1746-2852 with this nonthermal shell. We here show that the highest-energy component in Sgr A East's spectrum, like that of SNRs, can be fitted with the gamma-rays produced in neutral pion decays. Further, we demonstrate in a self-consistent manner that the leptons released in the associated charged pion decays produce an e+e- distribution that can mimic a power-law with index ~3, like that inferred from the VLA data for this source. These relativistic electrons and positrons also radiate by bremsstrahlung, and inverse Compton scattering with the intense IR and UV fields from the nucleus. We show that the overall emissivity calculated in this way may account for Sgr A East's broadband spectrum ranging from GHz frequencies all the way to TeV energies, where Whipple has thus far set an upper limit to the flux corresponding to a 2.5-sigma signal above the noise.

EGRET (20 MeV to 30 GeV) on board the Compton GRO has observed high-energy emission from about 40-50 Active Galactic Nuclei. Theoretical models of this emission based on the upscattering of thermal disk photons by cooling, relativistic electrons can successfully account for the EGRET observations, but they predict a considerably greater X-ray flux than that actually observed in a majority of these sources. This inconsistency may be an indication that the particles are energized during the Compton scattering process, since the X-ray emission is produced by the lowest energy electrons, whose density may be relatively small due to the acceleration. Such a situation may arise as a result of resistive field generation in electromagnetic acceleration schemes, which we here explore. A key feature of this model is the assumed existence of a current associated with the azimuthal component B_phi of the underlying magnetic field by a slight imbalance in the energy distributions of outwardly moving, relativistic electrons and protons produced at the disk surface via shock acceleration. The generation of an electric field (via magnetic field line reconnection) is thus required to maintain the current in the presence of a resistivity induced by the radiative drag on the relativistic electrons. We show that the resulting spectrum can exhibit a significant deficit of X-rays compared with gamma-rays. In addition, due to the uni-directional flow of the current associated with B_phi, this model would predict that the electrons are energized relative to the protons in the outflow only on one side of the disk. They should be decelerated on the reverse side. As such, we would anticipate that any given blazar should have a ~50% probability of being gamma-bright, which appears to be consistent with the observed ratio.

Continuum low-frequency radio observations of the Galactic Center reveal the presence of two prominent radio sources, Sgr A East and its surrounding Halo, containing non-thermal particle distributions with power-law indices around 2.5-3.3 and 2.4, respectively. The central 1-2 pc region is also a source of intense (stellar) UV and (dust-reprocessed) far-IR radiation that bathes these extended synchrotron-emitting structures. A recent detection of gamma-rays (2EGJ1746-2852) from within around 1 degree of the Galactic Center by EGRET onboard the Compton GRO shows that the emission from this environment extends to very high energies. We suggest that inverse Compton scatterings between the power-law electrons inferred from the radio properties of Sgr A East and its Halo, and the UV and IR photons from the nucleus, may account for the possibly diffuse gamma-ray source as well. We show that both particle distributions may be contributing to the gamma-ray emission, though their relevant strength depends on the actual physical properties (such as the magnetic field intensity) in each source. If this picture is correct, the high-energy source at the Galactic Center is extended over several arcminutes, which can be tested with thenext generation of gamma-ray and hard X-ray missions.

We present models for the initial mass function (IMF) for stars forming within molecular clouds. These models use the idea that stars determine their own masses through the action of powerful stellar outflows. This concept allows us to calculate a semi-empirical mass formula (SEMF), which provides the transformation between initial conditions in molecular clouds and the final masses of forming stars. For a particular SEMF, a given distribution of initial conditions predicts a corresponding IMF. We consider several different descriptions for the distribution of initial conditions in star forming molecular clouds. We first consider the limiting case in which only one physical variable -- the effective sound speed -- determines the initial conditions. In this limit, we use observed scaling laws to determine the distribution of sound speed and the SEMF to convert this distribution into an IMF. We next consider the opposite limit in which many different independent physical variables play a role in determining stellar masses. In this limit, the central limit theorem shows that the IMF approaches a log-normal form. Realistic star forming regions contain an intermediate number of relevant variables; we thus consider intermediate cases between the two limits. Our results show that this picture of star formation and the IMF naturally produces stellar mass distributions that are roughly consistent with observations. This paper thus provides a calculational framework to construct theoretical models of the IMF.

We study the linear evolution of small perturbations in self-gravitating fluid systems in two spatial dimensions; we consider both cylindrical and cartesian (i.e., slab) geometries. The treatment is general, but the application is to molecular clouds. We consider a class of equations of state which heuristically take into account the presence of turbulence; in particular, we consider equations of state which are {\it softer} than isothermal. We take the unperturbed cloud configuration to be in hydrostatic equilibrium. We find a class of wave solutions which propagate along a pressure supported cylinder (or slab) and have finite (trapped) spatial distributions in the direction perpendicular to the direction of propagation. Our results indicate that the dispersion relations for these two dimensional waves have similar forms for the two geometries considered here. Both cases possess a regime of instability and a fastest growing mode. We also find the (perpendicular) form of the perturbations for a wide range of wavelengths. Finally, we discuss the implications of our results for star formation and molecular clouds. The mass scales set by instabilities in both molecular cloud filaments and sheets are generally much larger than the masses of stars. However, these instabilities can determine the length scales for the initial conditions for protostellar collapse.

We study nonlinear wave phenomena in self-gravitating fluid systems, with a particular emphasis on applications to molecular clouds. This paper presents analytical results for one spatial dimension. We show that a large class of physical systems can be described by theories with a ``charge density'' $q(\rho)$; this quantity replaces the density on the right hand side of the Poisson equation for the gravitational potential. We use this formulation to prove general results about nonlinear wave motions in self-gravitating systems. We show that in order for stationary waves to exist, the total charge (the integral of the charge density over the wave profile) must vanish. This ``no-charge'' property for solitary waves is related to the capability of a system to be stable to gravitational perturbations for arbitrarily long wavelengths. We find necessary and sufficient conditions on the charge density for the existence of solitary waves and stationary waves. We then apply these results to study nonlinear waves in Jeans-type theories, Yukawa theories, and a pseudo-two-dimensional treatment. We study the allowed types of wave behavior for all of these models. Finally, we discuss the implications of this work for molecular cloud structure.

Highly dynamic magnetospheric perturbations in neutron star environments can naturally account for the features observed in Gamma-ray Burst spectra. The source distribution, however, appears to be extragalactic. Although noncatastrophic isotropic emission mechanisms may be ruled out on energetic and timing arguments, MHD processes can produce strongly anisotropic gamma-rays with an observable flux out to distances of around 1-2 Gpc. Here we show that sheared Alfven waves propagating along open magnetospheric field lines at the poles of magnetized neutron stars transfer their energy dissipationally to the current sustaining the field misalignment and thereby focus their power into a spatial region around 1000 times smaller than that of the crustal disturbance. This produces a strong (observable) flux enhancement along certain directions. We apply this model to a source population of ``turned-off'' pulsars that have nonetheless retained their strong magnetic fields and have achieved alignment at a period of greater than 5 seconds.