We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system.

AU Microscopii (AU Mic) is the second closest pre main sequence star, at a distance of 9.79 parsecs and with an age of 22 million years. AU Mic possesses a relatively rare and spatially resolved3 edge-on debris disk extending from about 35 to 210 astronomical units from the star, and with clumps exhibiting non-Keplerian motion. Detection of newly formed planets around such a star is challenged by the presence of spots, plage, flares and other manifestations of magnetic activity on the star. Here we report observations of a planet transiting AU Mic. The transiting planet, AU Mic b, has an orbital period of 8.46 days, an orbital distance of 0.07 astronomical units, a radius of 0.4 Jupiter radii, and a mass of less than 0.18 Jupiter masses at 3 sigma confidence. Our observations of a planet co-existing with a debris disk offer the opportunity to test the predictions of current models of planet formation and evolution.

Exoplanets residing close to their stars can experience evolution of both their physical structures and their orbits due to the influence of their host stars. In this work, we present a coupled analysis of dynamical tidal dissipation and atmospheric mass loss for exoplanets in XUV irradiated environments. As our primary application, we use this model to study the TRAPPIST-1 system, and place constraints on the interior structure and orbital evolution of the planets. We start by reporting on a UV continuum flux measurement (centered around $\sim1900$ Angstroms) for the star TRAPPIST-1, based on 300 ks of Neil Gehrels Swift Observatory data, and which enables an estimate of the XUV-driven thermal escape arising from XUV photo-dissociation for each planet. We find that the X-ray flaring luminosity, measured from our X-ray detections, of TRAPPIST-1 is 5.6 $\times$10$^{-4} L_{*}$, while the full flux including non-flaring periods is 6.1 $\times$10$^{-5} L_{*}$, when $L_{*}$ is TRAPPIST-1's bolometric luminosity. We then construct a model that includes both atmospheric mass-loss and tidal evolution, and requires the planets to attain their present-day orbital elements during this coupled evolution. We use this model to constrain the ratio $Q'=3Q/2k_{2}$ for each planet. Finally, we use additional numerical models implemented with the Virtual Planet Simulator \texttt{VPLanet} to study ocean retention for these planets using our derived system parameters.

The most widely-studied mechanism of mass loss from extrasolar planets is photoevaporation via XUV ionization, primarily in the context of highly irradiated planets. However, the EUV dissociation of hydrogen molecules can also theoretically drive atmospheric evaporation on low-mass planets. For temperate planets such as the early Earth, impact erosion is expected to dominate in the traditional planetesimal accretion model, but it would be greatly reduced in pebble accretion scenarios, allowing other mass loss processes to be major contributors. We apply the same prescription for photoionization to this photodissociation mechanism and compare it to an analysis of other possible sources of mass loss in pebble accretion scenarios. We find that there is not a clear path to evaporating the primordial atmosphere accreted by an early Earth analog in a pebble accretion scenario. Impact erosion could remove ~2,300 bars of hydrogen if 1% of the planet's mass is accreted as planetesimals, while the combined photoevaporation processes could evaporate ~750 bars of hydrogen. Photodissociation is likely a subdominant, but significant component of mass loss. Similar results apply to super-Earths and mini-Neptunes. This mechanism could also preferentially remove hydrogen from a planet's primordial atmosphere, thereby leaving a larger abundance of primordial water compared to standard dry formation models. We discuss the implications of these results for models of rocky planet formation including Earth's formation and the possible application of this analysis to mass loss from observed exoplanets.

In young multiple stellar systems the gravitational influence of the stars shapes the circumstellar disk, controlling accretion and the material available for planet formation. Our observations of the triple system GW Orionis provide direct evidence for disk tearing, where the gravitational torques due to misaligned disk/orbit planes warp the disk and break the inner disk into precessing rings. We image an eccentric ring that is misaligned with respect to the orbital planes and outer disk, and casts shadows on the strongly warped intermediate disk. Our ring/warp geometry constraints and the fully characterized perturber orbits make the system a potential Rosetta Stone for studying disk hydrodynamics. The ring might offer suitable conditions for planet formation, providing a mechanism for forming wide-separation planets on highly oblique orbits.

Observations indicate that members of close stellar binaries often have mass ratios close to unity, while longer-period systems exhibit a more uniform mass-ratio distribution. This paper provides a theoretical explanation for this finding by determining the tidal equilibrium states for binary star systems --- subject to the constraints of conservation of angular momentum and constant total mass. This work generalizes previous treatments by including the mass fraction as a variable in the optimization problem. The results show that the lowest energy state accessible to the system corresponds to equal mass stars on a circular orbit, where the stellar spin angular velocities are both synchronized and aligned with the orbit. These features are roughly consistent with observed properties of close binary systems. We also find the conditions required for this minimum energy state to exist: [1] The total angular momentum must exceed a critical value, [2] the orbital angular momentum must be three times greater than the total spin angular momentum, and [3] the semimajor axis is bounded from above. The last condition implies that sufficiently wide binaries are not optimized with equal mass stars, where the limiting binary separation occurs near $a_0\approx16R_\ast$.

We test whether the population of "extreme" trans-Neptunian objects (eTNOs) detected in the Y4 Dark Energy Survey (DES) data exhibit azimuthal asymmetries which might be evidence of gravitational perturbations from an unseen super-Earth in a distant orbit. By rotating the orbits of the detected eTNOs, we construct a synthetic population which, when subject to the DES selection function, reproduces the detected distribution of eTNOs in the orbital elements $a,e,$ and $i$ as well as absolute magnitude $H$, but has uniform distributions in mean anomaly $M$, longitude of ascending node $\Omega,$ and argument of perihelion $\omega.$ We then compare the detected distributions in each of $\Omega, \omega,$ and $\varpi\equiv\Omega+\omega$ to those expected from the isotropic population, using Kuiper's variant of the Kolmogorov-Smirnov test. The three angles are tested for each of 4 definitions of the eTNO population, choosing among $a>(150,250)$~AU and perihelion $q>(30,37)$~AU. These choices yield 3--7 eTNOs in the DES Y4 sample. Among the twelve total tests, two have the likelihood of drawing the observed angles from the isotropic population at $p<0.05.$ The 3 detections at $a>250, q>37$~AU, and the 4 detections at $a>250, q>30$~AU, have $\Omega$ distribution with $p=0.03$ of coming from the isotropic construction, but this is not strong evidence of anisotropy given the 12 different tests. The DES data taken on their own are thus consistent with azimuthal isotropy and do not require a "Planet 9" hypothesis. The limited sky coverage and object count mean, however, that the DES data by no means falsify this hypothesis.

Motivated by the trends found in the observed sample of extrasolar planets, this paper determines tidal equilibrium states for forming planetary systems --- subject to conservation of angular momentum, constant total mass, and fixed orbital spacing. In the low-mass limit, valid for superearth-class planets with masses of order $m_{\rm p}\sim10M_\oplus$, previous work showed that energy optimization leads to nearly equal mass planets, with circular orbits confined to a plane. The present treatment generalizes previous results by including the self-gravity of the planetary bodies. For systems with sufficiently large total mass $m_{\scriptstyle T}$ in planets, the optimized energy state switches over from the case of nearly equal mass planets to a configuration where one planet contains most of the material. This transition occurs for a critical mass threshold of approximately $m_{\scriptstyle T} \ge m_{\scriptstyle C}\sim40M_\oplus$ (where the value depends on the semimajor axes of the planetary orbits, the stellar mass, and other system properties). These considerations of energy optimization apply over a wide range of mass scales, from binary stars to planetary systems to the collection of moons orbiting the giant planets in our solar system.

The planet Mercury possesses an anomalously large iron core, and a correspondingly high bulk density. Numerous hypotheses have been proposed in order to explain such a large iron content. A long-standing idea holds that Mercury once possessed a larger silicate mantle which was removed by a giant impact early in the the Solar system's history. A central problem with this idea has been that material ejected from Mercury is typically re-accreted onto the planet after a short (~Myr) timescale. Here, we show that the primordial Solar wind would have provided sufficient drag upon ejected debris to remove them from Mercury-crossing trajectories before re-impacting the planet's surface. Specifically, the young Sun likely possessed a stronger wind, fast rotation and strong magnetic field. Depending upon the time of the giant impact, the ram pressure associated with this wind would push particles outward into the Solar system, or inward toward the Sun, on sub-Myr timescales, depending upon the size of ejected debris. Accordingly, the giant impact hypothesis remains a viable pathway toward the removal of planetary mantles, both on Mercury and extrasolar planets, particularly those close to young stars with strong winds.

Most planetary systems -- including our own -- are born within stellar clusters, where interactions with neighboring stars can help shape the system architecture. This paper develops an orbit-averaged formalism to characterize the cluster's mean-field effects as well as the physics of long-period stellar encounters. Our secular approach allows for an analytic description of the dynamical consequences of the cluster environment on its constituent planetary systems. We analyze special cases of the resulting Hamiltonian, corresponding to eccentricity evolution driven by planar encounters, as well as hyperbolic perturbations upon dissipative disks. We subsequently apply our results to the early evolution of our solar system, where the cluster's collective potential perturbs the solar system's plane, and stellar encounters act to increase the velocity dispersion of the Kuiper belt. Our results are two-fold: first, we find that cluster effects can alter the mean plane of the solar system by $\lesssim1\deg$, and are thus insufficient to explain the $\psi\approx6\deg$ obliquity of the sun. Second, we delineate the extent to which stellar flybys excite the orbital dispersion of the cold classical Kuiper belt, and show that while stellar flybys may grow the cold belt's inclination by the observed amount, the resulting distribution is incompatible with the data. Correspondingly, our calculations place an upper limit on the product of the stellar number density and residence time of the sun in its birth cluster, $\eta\,\tau\lesssim2\times10^4\,$Myr/pc$^3$.

We present a catalog of 316 trans-Neptunian bodies detected by the Dark Energy Survey (DES). These objects include 245 discoveries by DES (139 not previously published) detected in $\approx 60,000$ exposures from the first four seasons of the survey ("Y4" data). The survey covers a contiguous 5000 deg$^2$ of the southern sky in the $grizY$ optical/NIR filter set, with a typical TNO in this part of the sky being targeted by $25-30$ Y4 exposures. We describe the processes for detection of transient sources and the linkage into TNO orbits, which are made challenging by the absence of the few-hour repeat observations employed by TNO-optimized surveys. We also describe the procedures for determining detection efficiencies vs. magnitude and estimating rates of false-positive linkages. This work presents all TNOs which were detected on $\ge 6$ unique nights in the Y4 data and pass a "sub-threshold confirmation" test wherein we demand the the object be detectable in a stack of the individual images in which the orbit indicates an object should be present, but was not detected. This eliminates false positives and yields TNO detections complete to $r\lesssim 23.3$ mag with virtually no dependence on orbital properties for bound TNOs at distance $30\,{\rm AU}<d<2500\,{\rm AU}.$ The final DES TNO catalog is expected to yield $>0.3$ mag more depth, and arcs of $>4$ years for nearly all detections.

Most of the major planets in the Solar System support populations of co-orbiting bodies, known as Trojans, at their L4 and L5 Lagrange points. In contrast, Earth has only one known co-orbiting companion. This paper presents the results from a search for Earth Trojans using the DECam instrument on the Blanco Telescope at CTIO. This search found no additional Trojans in spite of greater coverage compared to previous surveys of the L5 point. Therefore, the main result of this work is to place the most stringent constraints to date on the population of Earth Trojans. These constraints depend on assumptions regarding the underlying population properties, especially the slope of the magnitude distribution (which in turn depends on the size and albedo distributions of the objects). For standard assumptions, we calculate upper limits to a 90% confidence limit on the L5 population of $N_{ET}<1$ for magnitude $H<15.5$, $N_{ET}=60-85$ for $H<19.7$, and $N_{ET}\ $= 97 for $H=20.4$. This latter magnitude limit corresponds to Trojans $\sim$300 m in size for albedo $0.15$. At H=19.7, these upper limits are consistent with previous L4 Earth Trojan constraints and significantly improve L5 constraints.

The outer Solar System contains a large number of small bodies (known as trans-Neptunian objects or TNOs) that exhibit diverse types of dynamical behavior. The classification of bodies in this distant region into dynamical classes -- sub-populations that experience similar orbital evolution -- aids in our understanding of the structure and formation of the Solar System. In this work, we propose an updated dynamical classification scheme for the outer Solar System. This approach includes the construction of a new (automated) method for identifying mean-motion resonances. We apply this algorithm to the current dataset of TNOs observed by the Dark Energy Survey (DES) and present a working classification for all of the DES TNOs detected to date. Our classification scheme yields 1 inner centaur, 19 outer centaurs, 21 scattering disk objects, 47 detached TNOs, 48 securely resonant objects, 7 resonant candidates, and 97 classical belt objects. Among the scattering and detached objects, we detect 8 TNOs with semi-major axes greater than 150 AU.

2I/Borisov is the first-ever observed interstellar comet (and the second detected interstellar object). It was discovered on 30 August 2019 and has a heliocentric orbital eccentricity of ~ 3.35, corresponding to a hyperbolic orbit that is unbound to the Sun. Given that it is an interstellar object, it is of interest to compare its properties -- such as composition and activity -- with the comets in our Solar System. This study reports low-resolution optical spectra of 2I/Borisov. The spectra were obtained by the MDM observatory Hiltner 2.4m telescope/Ohio State Multi-Object Spectrograph (on 1 and 5 November 2019). The wavelength coverage spanned from 3700A to 9200A. The dust continuum reflectance spectra of 2I/Borisov show that the spectral slope is stepper in the blue end of the spectrum (compared to the red). The spectra of 2I/Borisov clearly show CN emission at 3880A, as well as C2 emission at both 4750A and 5150A. Using a Haser model to covert the observed fluxes into estimates for the molecular production rates, we find Q(CN) = 2.4 +/- 0.2 x 10^24 s^-1, and Q(C2) = 5.5 +/- 0.4 x 10^23 s^-1 at the heliocentric distance of 2.145 au. Our Q(CN) estimate is consistent with contemporaneous observations, and the Q(C2) estimate is generally below the upper limits of previous studies. We derived the ratio Q(C2)/Q(CN) = 0.2 +/- 0.1, which indicates that 2I/Borisov is depleted in carbon chain species, but is not empty. This feature is not rare for the comets in our Solar System, especially in the class of Jupiter Family Comets.

We are undertaking a large survey of over thirty disks using the Gemini Planet Imager (GPI) to see whether the observed dust structures match spectral energy distribution (SED) predictions and have any correlation with stellar properties. GPI can observe near-infrared light scattered from dust in circumstellar environments using high-resolution Polarimetric Differential Imaging (PDI) with coronagraphy and adaptive optics. The data have been taken in J and H bands over two years, with inner working angles of 0.08'' and 0.11'' respectively. Ahead of the release of the complete survey results, here we present five objects with extended and irregular dust structures within 2'' of the central star. These objects are: FU Ori; MWC 789; HD 45677; Hen 3-365; and HD 139614. The observed structures are consistent with each object being a pre-main-sequence star with protoplanetary dust. The five objects' circumstellar environments could result from extreme youth and complex initial conditions, from asymmetric scattering patterns due to shadows cast by misaligned disks, or in some cases from interactions with companions. We see complex U_phi structures in most objects that could indicate multiple scattering or result from the illumination of companions. Specific key findings include the first high-contrast observation of MWC 789 revealing a newly-discovered companion candidate and arc, and two faint companion candidates around Hen 3-365. These two objects should be observed further to confirm whether the companion candidates are co-moving. Further observations and modeling are required to determine the causes of the structures.

The observed distribution of radii for exoplanets shows a bimodal form that can be explained by mass-loss from planetary atmospheres due to high energy radiation emitted by their host stars. The location of the minimum of this radius distribution depends on the mass-radius relation, which in turn depends on the composition of planetary cores. Current studies suggest that super-Earth and mini-Neptune planets have iron-rich and hence largely Earth-like composition cores. This paper explores how non-zero planetary magnetic fields can decrease the expected mass-loss rates from these planets. These lower mass-loss rates, in turn, affect the location of the minimum of the radius distribution and the inferred chemical composition of the planetary cores.

Current observations indicate that the planet formation process often produces multiple planet systems with nearly circular orbits, regular spacing, a narrow range of inclination angles, and similar planetary masses of order $m_{\rm p}\sim10M_\oplus$. Motivated by the observational sample, this paper determines the tidal equilibrium states for this class of extrasolar planetary systems. We start by considering two planet systems with fixed orbital spacing and variable mass ratios. The basic conjecture explored in this paper is that the planet formation process will act to distribute planetary masses in order to achieve a minimum energy state. The resulting minimum energy configuration --- subject to the constraint of constant angular momentum --- corresponds to circular orbits confined to a plane, with nearly equal planetary masses (as observed). We then generalize the treatment to include multiple planet systems, where each adjacent pair of planets attains its (local) tidal equilibrium state. The properties of observed planetary systems are close to those expected from this pairwise equilibrium configuration. In contrast, observed systems do not reside in a global minimum energy state. Both the equilibrium states of this paper and observed multi-planet systems, with planets of nearly equal mass on regularly spaced orbits, have an effective surface density of the form $\sigma\propto r^{-2}$, much steeper than most disk models.

In 2015, K2 observations of the bright (V = 8.9, K = 7.7) star HIP 41378 revealed a rich system of at least five transiting exoplanets, ranging in size from super-Earths to gas giants. The 2015 K2 observations only spanned 74.8 days, and the outer three long-period planets in the system were only detected with a single transit, so their orbital periods and transit ephemerides could not be determined at that time. Here, we report on 50.8 days of new K2 observations of HIP 41378 from summer 2018. These data reveal additional transits of the long-period planets HIP 41378 d and HIP 41378 f, yielding a set of discrete possible orbital periods for these two planets. We identify the most probable orbital periods for these two planets using our knowledge of the planets' transit durations, the host star's properties, the system's dynamics, and data from the ground-based HATNet, KELT, and WASP transit surveys. Targeted photometric follow-up during the most probable future transit times will be able to determine the planets' orbital periods, and will enable future observations with facilities like the James Webb Space Telescope. The methods developed herein to determine the most probable orbital periods will be important for long-period planets detected by the Transiting Exoplanet Survey Satellite, where similar period ambiguities will frequently arise due to the telescope's survey strategy.

Over the course of the past two decades, observational surveys have unveiled the intricate orbital structure of the Kuiper Belt, a field of icy bodies orbiting the Sun beyond Neptune. In addition to a host of readily-predictable orbital behavior, the emerging census of trans-Neptunian objects displays dynamical phenomena that cannot be accounted for by interactions with the known eight-planet solar system alone. Specifically, explanations for the observed physical clustering of orbits with semi-major axes in excess of $\sim250\,$AU, the detachment of perihelia of select Kuiper belt objects from Neptune, as well as the dynamical origin of highly inclined/retrograde long-period orbits remain elusive within the context of the classical view of the solar system. This newly outlined dynamical architecture of the distant solar system points to the existence of a new planet with mass of $m_9\sim 5-10\,M_{\oplus}$, residing on a moderately inclined orbit ($i_9\sim15-25\deg$) with semi-major axis $a_9\sim 400 - 800\,$AU and eccentricity between $e_9 \sim 0.2 - 0.5$. This paper reviews the observational motivation, dynamical constraints, and prospects for detection of this proposed object known as Planet Nine.

(abridged) Both fundamental constants that describe the laws of physics and cosmological parameters that determine the cosmic properties must fall within a range of values in order for the universe to develop astrophysical structures and ultimately support life. This paper reviews current constraints on these quantities. The standard model of particle physics contains both coupling constants and particle masses, and the allowed ranges of these parameters are discussed first. We then consider cosmological parameters, including the total energy density, the vacuum energy density, the baryon-to-photon ratio, the dark matter contribution, and the amplitude of primordial density fluctuations. These quantities are constrained by the requirements that the universe lives for a long time, emerges from the BBN epoch with an acceptable chemical composition, and can successfully produce galaxies. On smaller scales, stars and planets must be able to form and function. The stars must have sufficiently long lifetimes and hot surface temperatures. The planets must be massive enough to maintain an atmosphere, small enough to remain non-degenerate, and contain enough particles to support a complex biosphere. These requirements place constraints on the gravitational constant, the fine structure constant, and composite parameters that specify nuclear reaction rates. We consider specific instances of possible fine-tuning in stars, including the triple alpha reaction that produces carbon, as well as the effects of unstable deuterium and stable diprotons. For all of these issues, viable universes exist over a range of parameter space, which is delineated herein. Finally, for universes with significantly different parameters, new types of astrophysical processes can generate energy and support habitability.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

We report the Transiting Exoplanet Survey Satellite ($TESS$) detection of a multi-planet system orbiting the $V=10.9$ K0 dwarf TOI 125. We find evidence for up to five planets, with varying confidence. Three high signal-to-noise transit signals correspond to sub-Neptune-sized planets ($2.76$, $2.79$, and $2.94\ R_{\oplus}$), and we statistically validate the planetary nature of the two inner planets ($P_b = 4.65$ days, $P_c = 9.15$ days). With only two transits observed, we report the outer object ($P_{.03} = 19.98$ days) as a high signal-to-noise ratio planet candidate. We also detect a candidate transiting super-Earth ($1.4\ R_{\oplus}$) with an orbital period of only $12.7$ hours and a candidate Neptune-sized planet ($4.2\ R_{\oplus}$) with a period of $13.28$ days, both at low signal-to-noise. This system is amenable to mass determination via radial velocities and transit timing variations, and provides an opportunity to study planets of similar size while controlling for age and environment. The ratio of orbital periods between TOI 125 b and c ($P_c/P_b = 1.97$) is slightly smaller than an exact 2:1 commensurability and is atypical of multiple planet systems from $Kepler$, which show a preference for period ratios just $wide$ of first-order period ratios. A dynamical analysis refines the allowed parameter space through stability arguments and suggests that, despite the nearly commensurate periods, the system is unlikely to be in resonance.

We present the first images of the transition disk around the close binary system HD 34700A in polarized scattered light using the Gemini Planet Imager instrument on Gemini South. The J and H band images reveal multiple spiral-arm structures outside a large (R = 0.49" = 175 au) cavity along with a bluish spiral structure inside the cavity. The cavity wall shows a strong discontinuity and we clearly see significant non-azimuthal polarization Uphi consistent with multiple scattering within a disk at an inferred inclination ~42deg. Radiative transfer modeling along with a new Gaia distance suggest HD 37400A is a young (~5 Myr) system consisting of two intermediate-mass (~2Msun) stars surrounded by a transitional disk and not a solar-mass binary with a debris disk as previously classified. Conventional assumptions of the dust-to-gas ratio would rule out a gravitational instability origin to the spirals while hydrodynamical models using the known external companion or a hypothetical massive protoplanet in the cavity both have trouble reproducing the relatively large spiral arm pitch angles (~30deg) without fine tuning of gas temperature. We explore the possibility that material surrounding a massive protoplanet could explain the rim discontinuity after also considering effects of shadowing by an inner disk. Analysis of archival Hubble Space Telescope data suggests the disk is rotating counterclockwise as expected from the spiral arm structure and revealed a new low-mass companion at 6.45" separation. We include an appendix which sets out clear definitions of Q, U, Qphi, Uphi, correcting some confusion and errors in the literature.

In the current model of early Solar System evolution, the stable members of the Jovian and Neptunian Trojan populations were captured into resonance from the leftover reservoir of planetesimals during the outward migration of the giant planets. As a result, both Jovian and Neptunian Trojans share a common origin with the primordial disk population, whose other surviving members constitute today's trans-Neptunian object (TNO) populations. The cold classical TNOs are ultra-red, while the dynamically excited "hot" population of TNOs contains a mixture of ultra-red and blue objects. In contrast, Jovian and Neptunian Trojans are observed to be blue. While the absence of ultra-red Jovian Trojans can be readily explained by the sublimation of volatile material from their surfaces due to the high flux of solar radiation at 5AU, the lack of ultra-red Neptunian Trojans presents both a puzzle and a challenge to formation models. In this work we report the discovery by the Dark Energy Survey (DES) of two new dynamically stable L4 Neptunian Trojans,2013 VX30 and 2014 UU240, both with inclinations i >30 degrees, making them the highest-inclination known stable Neptunian Trojans. We have measured the colors of these and three other dynamically stable Neptunian Trojans previously observed by DES, and find that 2013 VX30 is ultra-red, the first such Neptunian Trojan in its class. As such, 2013 VX30 may be a "missing link" between the Trojan and TNO populations. Using a simulation of the DES TNO detection efficiency, we find that there are 162 +/- 73 Trojans with Hr < 10 at the L4 Lagrange point of Neptune. Moreover, the blue-to-red Neptunian Trojan population ratio should be higher than 17:1. Based on this result, we discuss the possible origin of the ultra-red Neptunian Trojan population and its implications for the formation history of Neptunian Trojans.

We report the discovery of a compact multi-planet system orbiting the relatively nearby (78pc) and bright ($K=8.9$) K-star, K2-266 (EPIC248435473). We identify up to six possible planets orbiting K2-266 with estimated periods of P$_b$ = 0.66, P$_{.02}$ = 6.1, P$_c$ = 7.8, P$_d$ = 14.7, P$_e$ = 19.5, and P$_{.06}$ = 56.7 days and radii of R$_P$ = 3.3 R$_{\oplus}$, 0.646 R$_{\oplus}$, 0.705 R$_{\oplus}$, 2.93 R$_{\oplus}$, 2.73 R$_{\oplus}$, and 0.90 R$_{\oplus}$, respectively. We are able to confirm the planetary nature of two of these planets (d & e) from analyzing their transit timing variations ($m_d= 8.9_{-3.8}^{+5.7} M_\oplus$ and $m_e=14.3_{-5.0}^{+6.4} M_\oplus$), confidently validate the planetary nature of two other planets (b & c), and classify the last two as planetary candidates (K2-266.02 & .06). From a simultaneous fit of all 6 possible planets, we find that K2-266 b's orbit has an inclination of 75.32$^{\circ}$ while the other five planets have inclinations of 87-90$^{\circ}$. This observed mutual misalignment may indicate that K2-266 b formed differently from the other planets in the system. The brightness of the host star and the relatively large size of the sub-Neptune sized planets d and e make them well-suited for atmospheric characterization efforts with facilities like the Hubble Space Telescope and upcoming James Webb Space Telescope. We also identify an 8.5-day transiting planet candidate orbiting EPIC248435395, a co-moving companion to K2-266.

This paper reports the discovery and orbital characterization of two extreme trans-Neptunian objects (ETNOs), 2016 QV$_{89}$ and 2016 QU$_{89}$, which have orbits that appear similar to that of a previously known object, 2013 UH$_{15}$. All three ETNOs have semi-major axes $a\approx 172$ AU and eccentricities $e\approx0.77$. The angular elements $(i,\omega,\Omega)$ vary by 6, 15, and 49 deg, respectively between the three objects. The two new objects add to the small number of TNOs currently known to have semi-major axes between 150 and 250 AU, and serve as an interesting dynamical laboratory to study the outer realm of our Solar System. Using a large ensemble of numerical integrations, we find that the orbits are expected to reside in close proximity in the $(a,e)$ phase plane for roughly 100 Myr before diffusing to more separated values. We then explore other scenarios that could influence their orbits. With aphelion distances over 300 AU, the orbits of these ETNOs extend far beyond the classical Kuiper Belt, and an order of magnitude beyond Neptune. As a result, their orbital dynamics can be affected by the proposed new Solar System member, referred to as Planet Nine in this work. With perihelion distances of 35-40 AU, these orbits are also influenced by resonant interactions with Neptune. A full assessment of any possible, new Solar System planets must thus take into account this emerging class of TNOs.

Motivated by the possible existence of other universes, this paper considers the evolution of massive stars with different values for the fundamental constants. We focus on variations in the triple alpha resonance energy and study its effects on the resulting abundances of $^{12}$C, $^{16}$O, and larger nuclei. In our universe, the $0^{+}$ energy level of carbon supports a resonant nuclear reaction that dominates carbon synthesis in stellar cores and accounts for the observed cosmic abundances. Here we define $\Delta{E}_R$ to be the change in this resonant energy level, and show how different values affect the cosmic abundances of the intermediate alpha elements. Using the state of the art computational package $MESA$, we carry out stellar evolution calculations for massive stars in the range $M_\ast$ = $15-40M_\odot$, and for a wide range of resonance energies. We also include both solar and low metallicity initial conditions. For negative $\Delta{E}_R$ , carbon yields are increased relative to standard stellar models, and such universes remain viable as long as the production of carbon nuclei remains energetically favorable, and stars remain stable, down to $\Delta{E}_R\approx-300$ keV. For positive $\Delta{E}_R$, carbon yields decrease, but significant abundances can be produced for resonance energy increments up to $\Delta{E}_R\approx+500$ keV. Oxygen yields tend to be anti-correlated with those of carbon, and the allowed range in $\Delta{E}_R$ is somewhat smaller. We also present yields for neon, magnesium, and silicon. With updated stellar evolution models and a more comprehensive survey of parameter space, these results indicate that the range of viable universes is larger than suggested by earlier studies.

Motivated by the possibility that the laws of physics could be different in other regions of space-time, we consider nuclear processes in universes where the weak interaction is either stronger or weaker than observed. We focus on the physics of both Big Bang Nucleosynthesis (BBN) and stellar evolution. For sufficiently ineffective weak interactions, neutrons do not decay during BBN, and the baryon-to-photon ratio $\eta$ must be smaller in order for protons to survive without becoming incorporated into larger nuclei. For stronger weak interactions, neutrons decay before the onset of BBN, and the early universe is left with nearly a pure hydrogen composition. We then consider stellar structure and evolution for the different nuclear compositions resulting from BBN, a wide range of weak force strengths, and the full range of stellar masses for a given universe. We delineate the range of this parameter space that supports working stars, along with a determination of the dominant nuclear reactions over the different regimes. Deuterium burning dominates the energy generation in stars when the weak force is sufficiently weak, whereas proton-proton burning into helium-3 dominates for the regime where the weak force is much stronger than in our universe. Although stars in these universes are somewhat different, they have comparable surface temperatures, luminosities, radii, and lifetimes, so that a wide range of such universes remain potentially habitable.

We report the discovery and dynamical analysis of 2015 BP$_{519}$, an extreme Trans-Neptunian Object detected detected by the Dark Energy Survey at a heliocentric distance of 55 AU and absolute magnitude Hr= 4.3. The current orbit, determined from a 1110-day observational arc, has semi-major axis $a\approx$ 450 AU, eccentricity $e\approx$ 0.92 and inclination $i\approx$ 54 degrees. With these orbital elements, 2015 BP$_{519}$ is the most extreme TNO discovered to date, as quantified by the reduced Kozai action, which is is a conserved quantity at fixed semi-major axis $a$ for axisymmetric perturbations. We discuss the orbital stability and evolution of this object in the context of the known Solar System, and find that 2015 BP$_{519}$ displays rich dynamical behavior, including rapid diffusion in semi-major axis and more constrained variations in eccentricity and inclination. We also consider the long term orbital stability and evolutionary behavior within the context of the Planet Nine Hypothesis, and find that BP$_{519}$ adds to the circumstantial evidence for the existence of this proposed new member of the Solar System, as it would represent the first member of the population of high-i, $\varpi$-shepherded TNOs.

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.

Combining visual and spectroscopic orbits of binary stars leads to a determination of the full 3D orbit, individual masses, and distance to the system. We present a full analysis of the evolved binary system $\delta$ Delphini using astrometric data from the MIRC and PAVO instruments on the CHARA long-baseline interferometer, 97 new spectra from the Fairborn Observatory, and 87 unpublished spectra from Lick Observatory. We determine the full set of orbital elements for $\delta$ Del, along with masses of $1.78 \pm 0.07$ $M_{\odot}$ and $1.62 \pm 0.07$ $M_{\odot}$ for each component, and a distance of $63.61 \pm 0.89$ pc. These results are important in two contexts: for testing stellar evolution models and defining the detection capabilities for future planet searches. We find that the evolutionary state of this system is puzzling, as our measured flux ratios, radii, and masses imply a $\sim$ 200 Myr age difference between the components using standard stellar evolution models. Possible explanations for this age discrepancy include mass transfer scenarios with a now ejected tertiary companion. For individual measurements taken over a span of 2 years we achieve $<10$ $\mu$-arcsecond precision on differential position with 10-minute observations. The high precision of our astrometric orbit suggests that exoplanet detection capabilities are within reach of MIRC at CHARA. We compute exoplanet detection limits around $\delta$ Del, and conclude that if this precision is extended to wider systems we should be able to detect most exoplanets $>2$ M$_{J}$ on orbits $>0.75$ AU around individual components of hot binary stars via differential astrometry.

We investigate a class of universes in which the weak interaction is not in operation. We consider how astrophysical processes are altered in the absence of weak forces, including Big Bang Nucleosynthesis (BBN), galaxy formation, molecular cloud assembly, star formation, and stellar evolution. Without weak interactions, neutrons no longer decay, and the universe emerges from its early epochs with a mixture of protons, neutrons, deuterium, and helium. The baryon-to-photon ratio must be smaller than the canonical value in our universe to allow free nucleons to survive the BBN epoch without being incorporated into heavier nuclei. At later times, the free neutrons readily combine with protons to make deuterium in sufficiently dense parts of the interstellar medium, and provide a power source before they are incorporated into stars. Almost all of the neutrons are incorporated into deuterium nuclei before stars are formed. As a result, stellar evolution proceeds primarily through strong interactions, with deuterium first burning into helium, and then helium fusing into carbon. Low-mass deuterium-burning stars can be long-lived, and higher mass stars can synthesize the heavier elements necessary for life. Although somewhat different from our own, such universes remain potentially habitable.

The existence of hot Jupiters has challenged theories of planetary formation since the first extrasolar planets were detected. Giant planets are generally believed to form far from their host stars, where volatile materials like water exist in their solid phase, making it easier for giant planet cores to accumulate. Several mechanisms have been proposed to explain how giant planets can migrate inward from their birth sites to short-period orbits. One such mechanism, called Kozai-Lidov migration, requires the presence of distant companions in orbits inclined by more than $\sim40$ degrees with respect to the plane of the hot Jupiter's orbit. The high occurrence rate of wide companions in hot Jupiter systems lends support to this theory for migration. However, the exact orbital inclinations of these detected planetary and stellar companions is not known, so it is not clear whether the mutual inclination of these companions is large enough for the Kozai-Lidov process to operate. This paper shows that in systems orbiting cool stars with convective outer layers, the orbits of most wide planetary companions to hot Jupiters must be well aligned with the orbits of the hot Jupiters and the spins of the host stars. For a variety of possible distributions for the inclination of the companion, the width of the distribution must be less than $\sim20$ degrees to recreate the observations with good fidelity. As a result, the companion orbits are likely well-aligned with those of the hot Jupiters, and the Kozai-Lidov mechanism does not enforce migration in these systems.

We present precise radial velocity observations of WASP-47, a star known to host a hot Jupiter, a distant Jovian companion, and, uniquely, two additional transiting planets in short-period orbits: a super-Earth in a ~19 hour orbit, and a Neptune in a ~9 day orbit. We analyze our observations from the HARPS-N spectrograph along with previously published data to measure the most precise planet masses yet for this system. When combined with new stellar parameters and reanalyzed transit photometry, our mass measurements place strong constraints on the compositions of the two small planets. We find unlike most other ultra-short-period planets, the inner planet, WASP-47 e, has a mass (6.83 +/- 0.66 Me) and radius (1.810 +/- 0.027 Re) inconsistent with an Earth-like composition. Instead, WASP-47 e likely has a volatile-rich envelope surrounding an Earth-like core and mantle. We also perform a dynamical analysis to constrain the orbital inclination of WASP-47 c, the outer Jovian planet. This planet likely orbits close to the plane of the inner three planets, suggesting a quiet dynamical history for the system. Our dynamical constraints also imply that WASP-47 c is much more likely to transit than a geometric calculation would suggest. We calculate a transit probability for WASP-47 c of about 10%, more than an order of magnitude larger than the geometric transit probability of 0.6%.

We evaluate the dynamical stability of a selection of outer solar system objects in the presence of the proposed new Solar System member Planet Nine. We use a Monte Carlo suite of numerical N-body integrations to construct a variety of orbital elements of the new planet and evaluate the dynamical stability of eight Trans-Neptunian objects (TNOs) in the presence of Planet Nine. These simulations show that some combinations of orbital elements ($a,e$) result in Planet Nine acting as a stabilizing influence on the TNOs, which can otherwise be destabilized by interactions with Neptune. These simulations also suggest that some TNOs transition between several different mean-motion resonances during their lifetimes while still retaining approximate apsidal anti-alignment with Planet Nine. This behavior suggests that remaining in one particular orbit is not a requirement for orbital stability. As one product of our simulations, we present an {\it a posteriori} probability distribution for the semi-major axis and eccentricity of the proposed Planet Nine based on TNO stability. This result thus provides additional evidence that supports the existence of this proposed planet. We also predict that TNOs can be grouped into multiple populations of objects that interact with Planet Nine in different ways: one population may contain objects like Sedna and 2012 VP$_{113}$, which do not migrate significantly in semi-major axis in the presence of Planet Nine and tend to stay in the same resonance; another population may contain objects like 2007 TG$_{422}$ and 2013 RF$_{98}$, which may both migrate and transition between different resonances.

In both stars and in the early universe, the production of deuterium is the first step on the way to producing heavier nuclei. If the strong force were slightly weaker, deuterium would not be stable, and many authors have noted that nuclesynthesis would be compromised so that helium production could not proceed through standard reaction chains. Motivated by the possibility that other regions of space-time could have different values for the fundamental constants, this paper considers stellar evolution in universes without stable deuterium and argues that such universes can remain habitable. Even in universes with no stellar nucleosynthesis, stars can form and will generate energy through gravitational contraction. We show that such stars can be sufficiently luminous and long-lived to support life. Stars with initial masses that exceed the Chandrasekhar mass cannot be supported by degeneracy pressure and explode at the end of their contraction phase. The resulting explosive nucleosynthesis can provide the universe with some heavy elements. We also explore the possibility that helium can be produced in stellar cores through a triple-nucleon reaction (roughly analogous to the triple-alpha reaction). Next we show that with even trace amounts of heavy elements --- produced through the triple-nucleon process or by explosive nucleosynthesis --- the CNO cycle can operate and allow stars to function. Finally, we consider Big Bang Nucleosynthesis without stable deuterium and find that only trace amounts of helium are produced, with even smaller abundances of other nuclei. With stars evolving through gravitational contraction, explosive nucleosynthesis, the triple-nucleon reaction, and the CNO cycle, universes with no stable deuterium are thus potentially habitable, contrary to many previous claims.

Accretion onto Classical T Tauri stars is thought to take place through the action of magnetospheric processes, with gas in the inner disk being channeled onto the star's surface by the stellar magnetic field lines. Young stars are known to accrete material in a time-variable manner and the source of this variability remains an open problem, particularly on the shortest (~ day) timescales. Using one-dimensional time-dependent numerical simulations that follow the field line geometry, we find that for plausibly realistic young stars, steady-state transonic accretion occurs naturally in the absence of any other source of variability. However, we show that if the density in the inner disk varies smoothly in time with ~ day long time-scales (e.g., due to turbulence) this complication can lead to the development of shocks in the accretion column. These shocks propagate along the accretion column and ultimately hit the star, leading to rapid, large amplitude changes in the accretion rate. We argue that when these shocks hit the star the observed time-dependence will be a rapid increase in accretion luminosity followed by a slower decline and could be an explanation for some of the short period variability observed in accreting young stars. Our one-dimensional approach bridges previous analytic work to more complicated, multi-dimensional simulations, and observations.

We report the observation and physical characterization of the possible dwarf planet \UZ\ ("DeeDee"), a dynamically detached trans-Neptunian object discovered at 92 AU. This object is currently the second-most distant known trans-Neptunian object with reported orbital elements, surpassed in distance only by the dwarf planet Eris. The object was discovered with an $r$-band magnitude of 23.0 in data collected by the Dark Energy Survey between 2014 and 2016. Its 1140-year orbit has $(a,e,i) = (109~\mathrm{AU}, 0.65, 26.8^{\circ})$. It will reach its perihelion distance of 38 AU in the year 2142. Integrations of its orbit show it to be dynamically stable on Gyr timescales, with only weak interactions with Neptune. We have performed followup observations with ALMA, using 3 hours of on-source integration time to measure the object's thermal emission in the Rayleigh-Jeans tail. The signal is detected at 7$\sigma$ significance, from which we determine a $V$-band albedo of $13.1^{+3.3}_{-2.4}\mathrm{(stat)}^{+2.0}_{-1.4}\mathrm{(sys)}$ percent and a diameter of $635^{+57}_{-61}\mathrm{(stat)}^{+32}_{-39}\mathrm{(sys)}$~km, assuming a spherical body with uniform surface properties.

Motivated by the large number of compact extrasolar planetary systems discovered by the Kepler Mission, this paper considers perturbations due to possible additional outer planets. The discovered compact systems sometimes contain multiple transiting planets, so that their orbital angular momentum vectors are tightly aligned. Since planetary orbits are susceptible to forced oscillations of their inclination angles, the highly aligned nature of these systems places constraints on possible additional (non-transiting) planets. If planets in the outer regions of these solar systems have sufficiently large mass or sufficiently small semi-major axis, they will induce the compact inner orbits to oscillate in and out of a transiting configuration. This paper considers the dynamics of the compact systems discovered to host five or more planets. In order to not perturb these systems out of a continually, mutually transiting state, additional planetary companions must generally have periastron $p>10$ AU. Specific constraints are found for each of the 18 planetary systems considered, which are obtained by marginalising over other orbital parameters using three different choices of priors for the companion properties (a uniform prior, a transit-inspired prior, and an non-transiting disk prior). A separate ensemble of numerical experiments shows that these compact systems generally cannot contain Jupiter-analogs without disrupting the observed orbits. We also consider how these constraints depend on system properties and find that the surface density of the planetary system is one of the most important variables. Finally, we provide specific results for two systems, WASP-47 and Kepler-20, for which this analysis provides interesting constraints.

In order to look for signs of on-going planet formation in young disks, we carried out the first J-band polarized emission imaging of the Herbig Ae/Be stars HD 150193, HD 163296, and HD 169142 using the Gemini Planet Imager (GPI), along with new H band observations of HD 144432. We confirm the complex "double ring" structure for the nearly face-on system HD 169142 first seen in H-band, finding the outer ring to be substantially redder than the inner one in polarized intensity. Using radiative transfer modeling, we developed a physical model that explains the full spectral energy distribution (SED) and J- and H-band surface brightness profiles, suggesting that the differential color of the two rings could come from reddened starlight traversing the inner wall and may not require differences in grain properties. In addition, we clearly detect an elongated, off-center ring in HD 163296 (MWC 275), locating the scattering surface to be 18 AU above the midplane at a radial distance of 77 AU, co-spatial with a ring seen at 1.3mm by ALMA linked to the CO snow line. Lastly, we report a weak tentative detection of scattered light for HD 150193 (MWC 863) and a non-detection for HD 144432; the stellar companion known for each of these targets has likely disrupted the material in the outer disk of the primary star. For HD 163296 and HD 169142, the prominent outer rings we detect could be evidence for giant planet formation in the outer disk or a manifestation of large-scale dust growth processes possibly related to snow-line chemistry.

Mean motion commensurabilities in multi-planet systems are an expected outcome of protoplanetary disk-driven migration, and their relative dearth in the observational data presents an important challenge to current models of planet formation and dynamical evolution. One natural mechanism that can lead to the dissolution of commensurabilities is stochastic orbital forcing, induced by turbulent density fluctuations within the nebula. While this process is qualitatively promising, the conditions under which mean motion resonances can be broken are not well understood. In this work, we derive a simple analytic criterion that elucidates the relationship among the physical parameters of the system, and find the conditions necessary to drive planets out of resonance. Subsequently, we confirm our findings with numerical integrations carried out in the perturbative regime, as well as direct N-body simulations. Our calculations suggest that turbulent resonance disruption depends most sensitively on the planet-star mass ratio. Specifically, for a disk with properties comparable to the early solar nebula with $\alpha=0.01$, only planet pairs with cumulative mass ratios smaller than $(m_1+m_2)/M\lesssim10^{-5}\sim3M_{\oplus}/M_{\odot}$ are susceptible to breaking resonance at semi-major axis of order $a\sim0.1\,$AU. Although turbulence can sometimes compromise resonant pairs, an additional mechanism (such as suppression of resonance capture probability through disk eccentricity) is required to adequately explain the largely non-resonant orbital architectures of extrasolar planetary systems.

This paper derives an upper limit on the density $\rho_{\scriptstyle\Lambda}$ of dark energy based on the requirement that cosmological structure forms before being frozen out by the eventual acceleration of the universe. By allowing for variations in both the cosmological parameters and the strength of gravity, the resulting constraint is a generalization of previous limits. The specific parameters under consideration include the amplitude $Q$ of the primordial density fluctuations, the Planck mass $M_{\rm pl}$, the baryon-to-photon ratio $\eta$, and the density ratio $\Omega_M/\Omega_b$. In addition to structure formation, we use considerations from stellar structure and Big Bang Nucleosynthesis (BBN) to constrain these quantities. The resulting upper limit on the dimensionless density of dark energy becomes $\rho_{\scriptstyle\Lambda}/M_{\rm pl}^4<10^{-90}$, which is $\sim30$ orders of magnitude larger than the value in our universe $\rho_{\scriptstyle\Lambda}/M_{\rm pl}^4\sim10^{-120}$. This new limit is much less restrictive than previous constraints because additional parameters are allowed to vary. With these generalizations, a much wider range of universes can develop cosmic structure and support observers. To constrain the constituent parameters, new BBN calculations are carried out in the regime where $\eta$ and $G=M_{\rm pl}^{-2}$ are much larger than in our universe. If the BBN epoch were to process all of the protons into heavier elements, no hydrogen would be left behind to make water, and the universe would not be viable. However, our results show that some hydrogen is always left over, even under conditions of extremely large $\eta$ and $G$, so that a wide range of alternate universes are potentially habitable.

Motivated by the possible existence of other universes, with different values for the fundamental constants, this paper considers stellar models in universes where $^8$Be is stable. Many previous authors have noted that stars in our universe would have difficulty producing carbon and other heavy elements in the absence of the well-known $^{12}$C resonance at 7.6 MeV. This resonance is necessary because $^8$Be is unstable in our universe, so that carbon must be produced via the triple alpha reaction to achieve the requisite abundance. Although a moderate change in the energy of the resonance (200 -- 300 keV) will indeed affect carbon production, an even smaller change in the binding energy of beryllium ($\sim100$ keV) would allow $^8$Be to be stable. A stable isotope with $A=8$ would obviate the need for the triple alpha process in general, and the $^{12}$C resonance in particular, for carbon production. This paper explores the possibility that $^8$Be can be stable in other universes. Simple nuclear considerations indicate that bound states can be realized, with binding energy $\sim0.1-1$ MeV, if the fundamental constants vary by a $\sim {\rm few}-10$ percent. In such cases, $^8$Be can be synthesized through helium burning, and $^{12}$C can be produced later through nuclear burning of beryllium. This paper focuses on stellar models that burn helium into beryllium; once the universe in question has a supply of stable beryllium, carbon production can take place during subsequent evolution in the same star or in later stellar generations. Using both a semi-analytic stellar structure model as well as a state-of-the-art stellar evolution code, we find that viable stellar configurations that produce beryllium exist over a wide range of parameter space. Finally, we demonstrate that carbon can be produced during later evolutionary stages.

The bulk of X-ray emission from pre-main-sequence (PMS) stars is coronal in origin. We demonstrate herein that stars on Henyey tracks in the Hertzsprung-Russell diagram have lower $\log(L_X/L_\ast)$, on average, than stars on Hayashi tracks. This effect is driven by the decay of $L_X$ once stars develop radiative cores. $L_X$ decays faster with age for intermediate mass PMS stars, the progenitors of main sequence A-type stars, compared to those of lower mass. As almost all main sequence A-type stars show no detectable X-ray emission, we may already be observing the loss of their coronae during their PMS evolution. Although there is no direct link between the size or mass of the radiative core and $L_X$, the longer stars have spent with partially convective interiors, the weaker their X-ray emission becomes. This conference paper is a synopsis of Gregory, Adams and Davies (2016).

Motivated by the current search for exomoons, this paper considers the stability of tidal equilibrium for hierarchical three-body systems containing a star, a planet, and a moon. In this treatment, the energy and angular momentum budgets include contributions from the planetary orbit, lunar orbit, stellar spin, planetary spin, and lunar spin. The goal is to determine the optimized energy state of the system subject to the constraint of constant angular momentum. Due to the lack of a closed form solution for the full three-body problem, however, we must use use an approximate description of the orbits. We first consider the Keplerian limit and find that the critical energy states are saddle points, rather than minima, so that these hierarchical systems have no stable tidal equilibrium states. We then generalize the calculation so that the lunar orbit is described by a time-averaged version of the circular restricted three-body problem. In this latter case, the critical energy state is a shallow minimum, so that a tidal equilibrium state exists. In both cases, however, the lunar orbit for the critical point lies outside the boundary (roughly half the Hill radius) where (previous) numerical simulations indicate dynamical instability. These results suggest that star-planet-moon systems have no viable long-term stable states analogous to those found for two-body systems.

The Kepler mission has revealed a great diversity of planetary systems and architectures, but most of the planets discovered by Kepler orbit faint stars. Using new data from the K2 mission, we present the discovery of a five planet system transiting a bright (V = 8.9, K = 7.7) star called HIP 41378. HIP 41378 is a slightly metal-poor late F-type star with moderate rotation (v sin(i) = 7 km/s) and lies at a distance of 116 +/- 18 from Earth. We find that HIP 41378 hosts two sub-Neptune sized planets orbiting 3.5% outside a 2:1 period commensurability in 15.6 and 31.7 day orbits. In addition, we detect three planets which each transit once during the 75 days spanned by K2 observations. One planet is Neptune sized in a likely ~160 day orbit, one is sub-Saturn sized likely in a ~130 day orbit, and one is a Jupiter sized planet in a likely ~1 year orbit. We show that these estimates for the orbital periods can be made more precise by taking into account dynamical stability considerations. We also calculate the distribution of stellar reflex velocities expected for this system, and show that it provides a good target for future radial velocity observations. If a precise orbital period can be determined for the outer Jovian planet through future observations, it will be an excellent candidate for follow-up transit observations to study its atmosphere and measure its oblateness.

Motivated by the report of a possible new planetary member of the Solar System, this work calculates cross sections for interactions between passing stars and this proposed Planet Nine. Evidence for the new planet is provided by the orbital alignment of Kuiper Belt objects, and other Solar System properties, which suggest a Neptune-mass object on an eccentric orbit with semimajor axis $a_9\approx400-1500$ AU. With such a wide orbit, Planet Nine has a large interaction cross section, and is susceptible to disruption by passing stars. Using a large ensemble of numerical simulations (several million), and Monte Carlo sampling, we calculate the cross sections for different classes of orbit-altering events: [A] scattering the planet into its proposed orbit from a smaller orbit, [B] ejecting it from the Solar System from its current orbit, [C] capturing the planet from another system, and [D] capturing a free-floating planet. Results are presented for a range of orbital elements with planetary mass $m_9=10M_{earth}$. Removing Planet Nine from the Solar System is the most likely outcome. Specifically, we obtain ejection cross sections $\sigma_{int}\sim5\times10^6$ AU$^2$ ($5\times10^4$ AU$^2$) for environments corresponding to the birth cluster (field). With these cross sections, Planet Nine is likely to be ejected if the Sun resides within its birth cluster longer than $\Delta{t} \gtrsim 100$ Myr. The probability of ejecting Planet Nine due to passing field stars is $\lesssim 3\%$ over the age of the Sun. Probabilities for producing the inferred Planet Nine orbit are low $(\lesssim 5\%)$.

This paper develops constraints on the values of the fundamental constants that allow universes to be habitable. We focus on the fine structure constant $\alpha$ and the gravitational structure constant $\alpha_G$, and find the region in the $\alpha$-$\alpha_G$ plane that supports working stars and habitable planets. This work is motivated, in part, by the possibility that different versions of the laws of physics could be realized within other universes. The following constraints are enforced: [A] long-lived stable nuclear burning stars exist, [B] planetary surface temperatures are hot enough to support chemical reactions, [C] stellar lifetimes are long enough to allow biological evolution, [D] planets are massive enough to maintain atmospheres, [E] planets are small enough in mass to remain non-degenerate, [F] planets are massive enough to support sufficiently complex biospheres, [G] planets are smaller in mass than their host stars, and [H] stars are smaller in mass than their host galaxies. This paper delineates the portion of the $\alpha$-$\alpha_G$ plane that satisfies all of these constraints. The results indicate that viable universes --- with working stars and habitable planets --- can exist within a parameter space where the structure constants $\alpha$ and $\alpha_G$ vary by several orders of magnitude. These constraints also provide upper bounds on the structure constants ($\alpha,\alpha_G$) and their ratio. We find the limit $\alpha_G/\alpha<10^{-34}$, which shows that habitable universes must have a large hierarchy between the strengths of the gravitational force and the electromagnetic force.

We present high-contrast Magellan adaptive optics (MagAO) images of HD 7449, a Sun-like star with one planet and a long-term radial velocity (RV) trend. We unambiguously detect the source of the long-term trend from 0.6-2.15 \microns ~at a separation of \about 0\fasec 54. We use the object's colors and spectral energy distribution to show that it is most likely an M4-M5 dwarf (mass \about 0.1-0.2 \msun) at the same distance as the primary and is therefore likely bound. We also present new RVs measured with the Magellan/MIKE and PFS spectrometers and compile these with archival data from CORALIE and HARPS. We use a new Markov chain Monte Carlo procedure to constrain both the mass ($> 0.17$ \msun ~at 99$\%$ confidence) and semimajor axis (\about 18 AU) of the M dwarf companion (HD 7449B). We also refine the parameters of the known massive planet (HD 7449Ab), finding that its minimum mass is $1.09^{+0.52}_{-0.19}$ \mj, its semimajor axis is $2.33^{+0.01}_{-0.02}$ AU, and its eccentricity is $0.8^{+0.08}_{-0.06}$. We use N-body simulations to constrain the eccentricity of HD 7449B to $\lesssim$ 0.5. The M dwarf may be inducing Kozai oscillations on the planet, explaining its high eccentricity. If this is the case and its orbit was initially circular, the mass of the planet would need to be $\lesssim$ 1.5 \mj. This demonstrates that strong constraints on known planets can be made using direct observations of otherwise undetectable long-period companions.

Pre-main sequence (PMS) stars of mass $\gtrsim0.35\,{\rm M}_\odot$ transition from hosting fully convective interiors to configurations with a radiative core and outer convective envelope during their gravitational contraction. This stellar structure change influences the external magnetic field topology and, as we demonstrate herein, affects the coronal X-ray emission as a stellar analog of the solar tachocline develops. We have combined archival X-ray, spectroscopic, and photometric data for $\sim$1000 PMS stars from five of the best studied star forming regions; the ONC, NGC 2264, IC 348, NGC 2362, and NGC 6530. Using a modern, PMS calibrated, spectral type-to-effective temperature and intrinsic colour scale, we deredden the photometry using colours appropriate for each spectral type, and determine the stellar mass, age, and internal structure consistently for the entire sample. We find that PMS stars on Henyey tracks have, on average, lower fractional X-ray luminosities ($L_{\rm X}/L_\ast$) than those on Hayashi tracks, where this effect is driven by changes in $L_{\rm X}$. X-ray emission decays faster with age for higher mass PMS stars. There is a strong correlation between $L_\ast$ and $L_{\rm X}$ for Hayashi track stars but no correlation for Henyey track stars. There is no correlation between $L_{\rm X}$ and radiative core mass or radius. However, the longer stars have spent with radiative cores, the less X-ray luminous they become. The decay of coronal X-ray emission from young early K to late G-type PMS stars, the progenitors of main sequence A-type stars, is consistent with the dearth of X-ray detections of the latter.

We explore the dynamics of magnetically controlled outflows from Hot Jupiters, where these flows are driven by UV heating from the central star. In these systems, some of the open field lines do not allow the flow to pass smoothly through the sonic point, so that steady-state solutions do not exist in general. This paper focuses on this type of magnetic field configuration, where the resulting flow becomes manifestly time-dependent. We consider the case of both steady heating and time-variable heating, and find the time scales for the corresponding time variations of the outflow. Because the flow cannot pass through the sonic transition, it remains subsonic and leads to so-called breeze solutions. One manifestation of the time variability is that the flow samples a collection of different breeze solutions over time, and the mass outflow rate varies in quasi-periodic fashion. Because the flow is subsonic, information can propagate inward from the outer boundary, which determines, in part, the time scale of the flow variability. This work finds the relationship between the outer boundary scale and the time scale of flow variations. In practice, the location of the outer boundary is set by the extent of the sphere of influence of the planet. The measured time variability can be used, in principle, to constrain the parameters of the system (e.g., the strengths of the surface magnetic fields).

Jupiter and Saturn play host to an impressive array of satellites, making it reasonable to suspect that similar systems of moons might exist around giant extrasolar planets. Furthermore, a significant population of such planets is known to reside at distances of several Astronomical Units (AU), leading to speculation that some moons thereof might support liquid water on their surfaces. However, giant planets are thought to undergo inward migration within their natal protoplanetary disks, suggesting that gas giants currently occupying their host star's habitable zone formed further out. Here we show that when a moon-hosting planet undergoes inward migration, dynamical interactions may naturally destroy the moon through capture into a so-called "evection resonance." Within this resonance, the lunar orbit's eccentricity grows until the moon eventually collides with the planet. Our work suggests that moons orbiting within about 10 planetary radii are susceptible to this mechanism, with the exact number dependent upon the planetary mass, oblateness and physical size. Whether moons survive or not is critically related to where the planet began its inward migration as well as the character of inter-lunar perturbations. For example, a Jupiter-like planet currently residing at 1AU could lose moons if it formed beyond 5AU. Cumulatively, we suggest that an observational census of exomoons could potentially inform us on the extent of inward planetary migration, for which no reliable observational proxy currently exists.

The Kepler Mission has detected dozens of compact planetary systems with more than four transiting planets. This sample provides a collection of close-packed planetary systems with relatively little spread in the inclination angles of the inferred orbits. A large fraction of the observational sample contains limited multiplicity, begging the question whether there is a true diversity of multi transiting systems, or if some systems merely possess high mutual inclinations, allowing them to appear as single-transiting systems in a transit-based survey. This paper begins an exploration of the effectiveness of dynamical mechanisms in exciting orbital inclination within exoplanetary systems of this class. For these tightly packed systems, we determine that the orbital inclination angles are not spread out appreciably through self-excitation. In contrast, the two Kepler multi-planet systems with additional non-transiting planets are susceptible to oscillations of their inclination angles, which means their currently observed configurations could be due to planet-planet interactions alone. We also provide constraints and predictions for the expected transit duration variations (TDVs) for each planet. In these multi-planet compact Kepler systems, oscillations of their inclination angles are remarkably hard to excite; as a result, they tend to remain continually mutually transiting (CMT-stable). We study this issue further by augmenting the planet masses and determining the enhancement factor required for oscillations to move the systems out of transit. The oscillations of inclination found here inform the recently suggested dichotomy in the sample of solar systems observed by Kepler.

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.

Using new data from the K2 mission, we show that WASP-47, a previously known hot Jupiter host, also hosts two additional transiting planets: a Neptune-sized outer planet and a super-Earth inner companion. We measure planetary properties from the K2 light curve and detect transit timing variations, confirming the planetary nature of the outer planet. We performed a large number of numerical simulations to study the dynamical stability of the system and to find the theoretically expected transit timing variations (TTVs). The theoretically predicted TTVs are in good agreement with those observed, and we use the TTVs to determine the masses of two planets, and place a limit on the third. The WASP-47 planetary system is important because companion planets can both be inferred by TTVs and are also detected directly through transit observations. The depth of the hot Jupiter's transits make ground-based TTV measurements possible, and the brightness of the host star makes it amenable for precise radial velocity measurements. The system serves as a Rosetta Stone for understanding TTVs as a planet detection technique.

Motivated by the possibility that different versions of the laws of physics could be realized within other universes, this paper delineates the galactic parameters that allow for habitable planets and revisits constraints on the amplitude $Q$ of the primordial density fluctuations. Previous work indicates that large values of $Q$ lead to galaxies so dense that planetary orbits cannot survive long enough for life to develop. Small values of $Q$ lead to delayed star formation, loosely bound galaxies, and compromised heavy element retention. This work generalizes previous treatments: [A] We consider models for the internal structure of galaxies and find the fraction of galactic real estate that allows stable, long-lived planetary orbits. [B] We perform a large ensemble of numerical simulations to estimate cross sections for the disruption of planetary orbits due to interactions with passing stars. [C] We consider disruption due to the background radiation fields produced by the galaxies. [D] One consequence of intense galactic background radiation fields is that some portion of the galaxy, denoted as the Galactic Habitable Zone, will provide the right flux levels to support habitable planets for essentially any planetary orbit. As $Q$ increases, the fraction of stars in a galaxy that allow for habitable planets decreases due to both orbital disruption and the intense background radiation. However, the outer parts of the galaxy always allow for habitable planets, so that the value of $Q$ does not have a well-defined upper limit. Moreover, some Galactic Habitable Zones are large enough to support more potentially habitable planets than the galaxies found in our universe. These results suggest that the possibilities for habitability in other universes are somewhat more favorable and far more diverse than previously imagined.

We report the discovery of the eighth and ninth known Trojans in stable orbits around Neptune's leading Lagrange point, L4. The objects 2014 QO$_{441}$ and 2014 QP$_{441}$ were detected in data obtained during the 2013-14 and 2014-15 observing seasons by the Dark Energy Survey, using the Dark Energy Camera (DECam) on the 4-meter Blanco telescope at Cerro Tololo Inter-American Observatory. Both are in high-inclination orbits (18.8$^{\circ}$ and 19.4$^{\circ}$ respectively). With an eccentricity of 0.104, 2014 QO$_{441}$ has the most eccentric orbit of the eleven known stable Neptune Trojans. Here we describe the search procedure and investigate the objects' long-term dynamical stability and physical properties.

Kepler-296 is a binary star system with two M-dwarf components separated by 0.2 arcsec. Five transiting planets have been confirmed to be associated with the Kepler-296 system; given the evidence to date, however, the planets could in principle orbit either star. This ambiguity has made it difficult to constrain both the orbital and physical properties of the planets. Using both statistical and analytical arguments, this paper shows that all five planets are highly likely to orbit the primary star in this system. We performed a Markov-Chain Monte Carlo simulation using a five transiting planet model, leaving the stellar density and dilution with uniform priors. Using importance sampling, we compared the model probabilities under the priors of the planets orbiting either the brighter or the fainter component of the binary. A model where the planets orbit the brighter component, Kepler-296A, is strongly preferred by the data. Combined with our assertion that all five planets orbit the same star, the two outer planets in the system, Kepler-296 Ae and Kepler-296 Af, have radii of 1.53 +/- 0.26 and 1.80 +/- 0.31 R_earth, respectively, and receive incident stellar fluxes of 1.40 +/- 0.23 and 0.62 +/- 0.10 times the incident flux the Earth receives from the Sun. This level of irradiation places both planets within or close to the circumstellar habitable zone of their parent star.

Winds from short-period Earth and Neptune mass exoplanets, driven by high energy radiation from a young star, may evaporate a significant fraction of a planet's mass. If the momentum flux from the evaporative wind is not aligned with the planet/star axis, then it can exert a torque on the planet's orbit. Using steady-state one-dimensional evaporative wind models we estimate this torque using a lag angle that depends on the product of the speed of the planet's upper atmosphere and a flow timescale for the wind to reach its sonic radius. We also estimate the momentum flux from time-dependent one-dimensional hydrodynamical simulations. We find that only in a very narrow regime in planet radius, mass and stellar radiation flux is a wind capable of exerting a significant torque on the planet's orbit. Similar to the Yarkovsky effect, the wind causes the planet to drift outward if atmospheric circulation is prograde (super-rotating) and in the opposite direction if the circulation is retrograde. A close-in super Earth mass planet that loses a large fraction of its mass in a wind could drift a few percent of its semi-major axis. While this change is small, it places constraints on the evolution of resonant pairs such as Kepler 36 b and c.

Most planetary systems are formed within stellar clusters, and these environments can shape their properties. This paper considers scattering encounters between solar systems and passing cluster members, and calculates the corresponding interaction cross sections. The target solar systems are generally assumed to have four giant planets, with a variety of starting states, including circular orbits with the semimajor axes of our planets, a more compact configuration, an ultra-compact state with multiple mean motion resonances, and systems with massive planets. We then consider the effects of varying the cluster velocity dispersion, the relative importance of binaries versus single stars, different stellar host masses, and finite starting eccentricities of the planetary orbits. For each state of the initial system, we perform an ensemble of numerical scattering experiments and determine the cross sections for eccentricity increase, inclination angle increase, planet ejection, and capture. This paper reports results from over 2 million individual scattering simulations. Using supporting analytic considerations, and fitting functions to the numerical results, we find a universal formula that gives the cross sections as a function of stellar host mass, cluster velocity dispersion, starting planetary orbital radius, and final eccentricity. The resulting cross sections can be used in a wide variety of applications. As one example, we revisit constraints on the birth aggregate of our Solar System due to dynamical scattering and find $N<10^4$ (consistent with previous estimates).

The recent (apparent) passage of the Voyager 1 spacecraft into interstellar space provides us with front-row seats to the complex interplay between the solar wind and the protective surrounding bubble known as heliosphere. The heliosphere extends radially out to $\sim100$ AU from the sun, and within this sphere of influence, the solar wind modulates the incoming flux of galactic cosmic rays (CRs), especially those at low energies. Newly formed stars, which support both strong magnetic fields and winds, are expected to produce analogous regions of CR exclusion, perhaps at elevated levels. Such young stars are encircled by molecular gas-rich disks, and the net removal of CRs from the circumstellar environment significantly reduces the expected CR ionization rate in the disk gas, most likely by many orders-of-magnitude. The loss of ionization reduces disk turbulence, and thereby affects both planet-formation and active chemical processes in the disk. We present models of CR exclusion and explore the implications for turbulence and for predicted chemical abundances. We also discuss means by which ALMA can be used to search for extrasolar heliosphere-analogs around young stars.

We present an observational and theoretical study of the primary ionizing agents (cosmic rays and X-rays) in the TW Hya protoplanetary disk. We use a set of resolved and unresolved observations of molecular ions and other molecular species, encompassing eleven lines total, in concert with a grid of disk chemistry models. The molecular ion constraints comprise new data from the Submillimeter Array on HCO$^+$, acquired at unprecedented spatial resolution, and data from the literature, including ALMA observations of N$_2$H$^+$. We vary the model incident CR flux and stellar X-ray spectra and find that TW Hya's HCO$^+$ and N$_2$H$^+$ emission are best fit by a moderately hard X-ray spectra, as would be expected during the "flaring" state of the star, and a low CR ionization rate, $\zeta_{\rm CR}\lesssim10^{-19}$ s$^{-1}$. This low CR rate is the first indication of the presence of CR exclusion by winds and/or magnetic fields in an actively accreting T Tauri disk system. With this new constraint, our best fit ionization structure predicts a low turbulence "dead-zone" extending from the inner edge of the disk out to $50-65$ AU. This region coincides with an observed concentration of millimeter grains, and we propose that the inner region of TW Hya is a dust (and possibly planet) growth factory as predicted by previous theoretical work.

Star formation proceeds via the collapse of a molecular cloud core over multiple dynamical timescales. Turbulence within cores results in a spatially non-uniform angular momentum of the cloud, causing a stochastic variation in orientation of the disk forming from the collapsing material. In the absence of star-disk angular momentum coupling, such disk-tilting would provide a natural mechanism for production of primordial spin-orbit misalignments in the resulting planetary systems. However, owing to high accretion rates in the embedded phase of star formation, the inner edge of the circumstellar disk extends down to the stellar surface, resulting in efficient gravitational and accretional angular momentum transfer between the star and the disk. Here, we demonstrate that the resulting gravitational coupling is sufficient to suppress any significant star-disk misalignment, with accretion playing a secondary role. The joint tilting of the star-disk system leads to a stochastic wandering of star-aligned bipolar outflows. Such wandering widens the effective opening angle of stellar outflows, allowing for more efficient clearing of the remainder of the protostar's gaseous envelope. Accordingly, the processes described in this work provide an additional mechanism responsible for sculpting the stellar Initial Mass Function (IMF).

This paper considers the stability of tidal equilibria for planetary systems in which stellar rotation provides a significant contribution to the angular momentum budget. We begin by applying classic stability considerations for two bodies to planetary systems --- where one mass is much smaller than the other. The application of these stability criteria to a subset of the Kepler sample indicates that the majority of the systems are not in a stable equilibrium state. Motivated by this finding, we generalize the stability calculation to include the quadrupole moment for the host star. In general, a stable equilibrium requires that the total system angular momentum exceeds a minimum value (denoted here as $L_X$) and that the orbital angular momentum of the planet exceeds a minimum fraction of the total. Most, but not all, of the observed planetary systems in the sample have enough total angular momentum to allow an equilibrium state. Even with the generalizations of this paper, however, most systems have too little orbital angular momentum (relative to the total) and are not in an equilibrium configuration. Finally, we consider the time evolution of these planetary systems; the results constrain the tidal quality factor of the stars and suggest that $10^6\le{Q_\ast}\le10^7$.

Hot Jupiters can experience mass loss driven by heating from UV radiation from their host stars, and this flow is often controlled by magnetic fields. More specifically, near the planetry surface, the magnetic pressure dominates the ram pressure of the outflow by several orders of magnitude. After leaving the vicinity of the planet, the flow must connect onto the background environment provided by the stellar wind and the stellar magnetic field. This contribution considers magnetically controlled planetary outflows and extends previous work by comparing two different geometries for the background magnetic field provided by the star. In the first case, stellar field is assumed to retain the form of a dipole, which is anti-aligned with the dipole field of the planet. In the second case, the stellar outflow opens up the stellar magnetic field structure so that the background field at the location of the planet is perpendicular to the planetary dipole. Using numerical simulations, we consider the launch of the planetary wind with these field configurations.

The chemical properties of protoplanetary disks are especially sensitive to their ionization environment. Sources of molecular gas ionization include cosmic rays, stellar X-rays and short-lived radionuclides, each of which varies with location in the disk. This behavior leads to a significant amount of chemical structure, especially in molecular ion abundances, which is imprinted in their submillimeter rotational line emission. Using an observationally motivated disk model, we make predictions for the dependence of chemical abundances on the assumed properties of the ionizing field. We calculate the emergent line intensity for abundant molecular ions and simulate sensitive observations with the Atacama Large Millimeter/Sub-millimeter Array (ALMA) for a disk at D=100 pc. The models readily distinguish between high ionization rates ($\zeta\gtrsim10^{-17}$ s$^{-1}$ per H$_2$) and below, but it becomes difficult to distinguish between low ionization models when $\zeta\lesssim10^{-19}$ s$^{-1}$. We find that \htdp\ emission is not detectable for sub-interstellar CR rates with ALMA (6h integration), and that \ntdp\ emission may be a more sensitive tracer of midplane ionization. HCO$^+$ traces X-rays and high CR rates ($\zeta_{\rm{CR}}\gtrsim10^{-17}$ s$^{-1}$), and provides a handle on the warm molecular ionization properties where CO is present in the gas. Furthermore, species like HCO$^+$, which emits from a wide radial region and samples a large gradient in temperature, can exhibit ring-like emission as a consequence of low-lying rotational level de-excitation near the star. This finding highlights a scenario where rings are not necessarily structural or chemical in nature, but simply a result of the underlying line excitation properties.

We consider the role magnetic fields play in guiding and controlling mass-loss via evaporative outflows from exoplanets that experience UV irradiation. First we present analytic results that account for planetary and stellar magnetic fields, along with mass-loss from both the star and planet. We then conduct series of numerical simulations for gas giant planets, and vary the planetary field strength, background stellar field strength, UV heating flux, and planet mass. These simulations show that the flow is magnetically controlled for moderate field strengths and even the highest UV fluxes, i.e., planetary surface fields $B_P\gtrsim 0.3$ gauss and fluxes $F_{UV}\sim10^{6}$ erg s$^{-1}$. We thus conclude that outflows from all hot Jupiters with moderate surface fields are magnetically controlled. The inclusion of magnetic fields highly suppresses outflow from the night-side of the planet. Only the magnetic field lines near the pole are open and allow outflow to occur. The fraction of open field lines depends sensitively on the strength (and geometry) of the background magnetic field from the star, along with the UV heating rate. The net effect of the magnetic field is to suppress the mass loss rate by (approximately) an order of magnitude. Finally, some open field lines do not allow the flow to pass smoothly through the sonic point; flow along these streamlines does not reach steady-state, resulting in time-variable mass-loss.

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 quest for Earth-like planets represents a major focus of current exoplanet research. While planets that are Earth-sized and smaller have been detected, these planets reside in orbits that are too close to their host star to allow liquid water on their surface. We present the detection of Kepler-186f, a 1.11+\-0.14 Earth radius planet that is the outermost of five planets - all roughly Earth-sized - that transit a 0.47+\-0.05 Rsun star. The intensity and spectrum of the star's radiation places Kepler-186f in the stellar habitable zone, implying that if Kepler-186f has an Earth-like atmosphere and H2O at its surface, then some of this H2O is likely to be in liquid form.

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.

This chapter concerns the long-term dynamical evolution of planetary systems from both theoretical and observational perspectives. We begin by discussing the planet-planet interactions that take place within our own Solar System. We then describe such interactions in more tightly-packed planetary systems. As planet-planet interactions build up, some systems become dynamically unstable, leading to strong encounters and ultimately either ejections or collisions of planets. After discussing the basic physical processes involved, we consider how these interactions apply to extrasolar planetary systems and explore the constraints provided by observed systems. The presence of a residual planetesimal disc can lead to planetary migration and hence cause instabilities induced by resonance crossing; however, such discs can also stabilise planetary systems. The crowded birth environment of a planetary system can have a significant impact: close encounters and binary companions can act to destabilise systems, or sculpt their properties. In the case of binaries, the Kozai mechanism can place planets on extremely eccentric orbits which may later circularise to produce hot Jupiters.

Intermediate mass stars and stellar remnants often host planets, and these dynamical systems evolve because of mass loss and tides. This paper considers the combined action of stellar mass loss and tidal dissipation on planetary orbits in order to determine the conditions required for planetary survival. Stellar mass loss is included using a so-called Jeans model, described by a dimensionless mass loss rate \gamma and an index \beta. We use an analogous prescription to model tidal effects, described here by a dimensionless dissipation rate \Gamma and two indices (q,p). The initial conditions are determined by the starting value of angular momentum parameter \eta (equivalently, the initial eccentricity) and the phase \theta of the orbit. Within the context of this model, we derive an analytic formula for the critical dissipation rate \Gamma, which marks the boundary between orbits that spiral outward due to stellar mass loss and those that spiral inward due to tidal dissipation. This analytic result \Gamma=\Gamma(\gamma,\beta,q,p,\eta,\theta) is essentially exact for initially circular orbits and holds to within an accuracy of 50% over the entire multi-dimensional parameter space, where the individual parameters vary by several orders of magnitude. For stars that experience mass loss, the stellar radius often displays quasi-periodic variations, which produce corresponding variations in tidal forcing; we generalize the calculation to include such pulsations using a semi-analytic treatment that holds to the same accuracy as the non-pulsating case. These results can be used in many applications, e.g., to predict/constrain properties of planetary systems orbiting white dwarfs.

The presence of giant gaseous planets that reside in close proximity to their host stars may be a consequence of large-scale radial migration through the proto-planetary nebulae. Within the context of this picture, significant orbital obliquities characteristic of a substantial fraction of such planets can be attributed to external torques that perturb the disks out of alignment with the spin axes of their host stars. Therefore, the acquisition of orbital obliquity exhibits sensitive dependence on the physics of disk-star interactions. Here, we analyze the primordial excitation of spin-orbit misalignment of Sun-like stars, in light of disk-star angular momentum transfer. We begin by calculating the stellar pre-main sequence rotational evolution, accounting for spin-up due to gravitational contraction and accretion as well as spin-down due to magnetic star-disk coupling. We devote particular attention to angular momentum transfer by accretion, and show that while generally subdominant to gravitational contraction, this process is largely controlled by the morphology of the stellar magnetic field (i.e. specific angular momentum accreted by stars with octupole-dominated surface fields is smaller than that accreted by dipole-dominated stars by an order of magnitude). Subsequently, we examine the secular spin-axis dynamics of disk-bearing stars, accounting for the time-evolution of stellar and disk properties and demonstrate that misalignments are preferentially excited in systems where stellar rotation is not overwhelmingly rapid. Moreover, we show that the excitation of spin-orbit misalignment occurs impulsively, through an encounter with a resonance between the stellar precession frequency and the disk-torquing frequency. Cumulatively, the model developed herein opens up a previously unexplored avenue towards understanding star-disk evolution and its consequences in a unified manner.

We present simple analytic solutions for the ionization rate $\zeta_{\rm{SLR}}$ arising from the decay of short-lived radionuclides (SLRs) within protoplanetary disks. We solve the radiative transfer problem for the decay products within the disk, and thereby allow for the loss of radiation at low disk surface densities; energy loss becomes important outside $R\gtrsim30$ for typical disk masses $M_g=0.04$ M$_\odot$. Previous studies of chemistry/physics in these disks have neglected the impact of ionization by SLRs, and often consider only cosmic rays (CRs), because of the high CR-rate present in the ISM. However, recent work suggests that the flux of CRs present in the circumstellar environment could be substantially reduced by relatively modest stellar winds, resulting in severely modulated CR ionization rates, $\zeta_{\rm{CR}}$, equal to or substantially below that of SLRs ($\zeta_{\rm{SLR}}\lesssim10^{-18}$ s$^{-1}$). We compute the net ionizing particle fluxes and corresponding ionization rates as a function of position within the disk for a variety of disk models. The resulting expressions are especially simple for the case of vertically gaussian disks (frequently assumed in the literature). Finally, we provide a power-law fit to the ionization rate in the midplane as a function of gas disk surface density and time. Depending on location in the disk, the ionization rates by SLRs are typically in the range $\zeta_{\rm{SLR}}\sim(1-10)\times10^{-19}$ s$^{-1}$.

This paper explores the effects of FUV radiation fields from external stars on circumstellar disk evolution. Disks residing in young clusters can be exposed to extreme levels of FUV flux from nearby OB stars, and observations show that disks in such environments are being actively photoevaporated. Typical FUV flux levels can be factors of \sim 10^{2} - 10^{4} higher than the interstellar value. These fields are effective in driving mass loss from circumstellar disks because they act at large radial distance from the host star, i.e., where most of the disk mass is located, and where the gravitational potential well is shallow. We combine viscous evolution (an \alpha-disk model) with an existing FUV photoevaporation model to derive constraints on disk lifetimes, and to determine disk properties as functions of time, including mass loss rates, disk masses, and radii. We also consider the effects of X-ray photoevaporation from the host star using an existing model, and show that for disks around solar-mass stars, externally-generated FUV fields are often the dominant mechanism in depleting disk material. For sufficiently large viscosities, FUV fields can efficiently photoevaporate disks over the entire range of parameter space. Disks with viscosity parameter \alpha = 10^{-3} are effectively dispersed within 1-3 Myr; for higher viscosities (\alpha = 10^{-2}) disks are dispersed within \sim 0.25 - 0.5 Myr. Furthermore, disk radii are truncated to less than \sim100 AU, which can possibly affect the formation of planets. Our model predictions are consistent with the range of observed masses and radii of proplyds in the Orion Nebula Cluster.

(Abridged) Cosmic rays (CRs) are thought to provide an important source of ionization in the outermost and densest regions of protoplanetary disks; however, it is unknown to what degree they are physically present. As is observed in the Solar System, stellar winds can inhibit the propagation of cosmic rays within the circumstellar environment and subsequently into the disk. In this work, we explore the hitherto neglected effects of cosmic ray modulation by both stellar winds and magnetic field structures and study how these processes act to reduce disk ionization rates. We construct a two-dimensional protoplanetary disk model of a T-Tauri star system, focusing on ionization from stellar and interstellar FUV, stellar X-ray photons, and cosmic rays. We show that stellar winds can power a Heliosphere-like analogue, i.e., a "T-Tauriosphere," diminishing cosmic ray ionization rates by several orders of magnitude at low to moderate CR energies (E_CR<1 GeV). We explore models of both the observed solar wind cosmic ray modulation and a highly simplified estimate for "elevated" cosmic ray modulation as would be expected from a young T-Tauri star. In the former (solar analogue) case, we estimate the ionization rate from galactic cosmic rays to be zeta_CR~(0.23-1.4) x 1e-18 s^-1. This range of values, which we consider to be the maximum CR ionization rate for the disk, is more than an order of magnitude lower than what is generally assumed in current models for disk chemistry and physics. In the latter case, i.e., for a "T-Tauriosphere," the ionization rate by cosmic rays is zeta_CR<1e-20 s^-1, which is 1000 times smaller than the interstellar value. Indeed, if winds are as efficient at cosmic ray modulation as predicted here, short-lived radionuclides (now extinct) would have provided the major source of ionization (zeta_RN~7.3e-19 s^-1) in the planet-forming zone of the young Solar Nebula.

Hill's equations arise in a wide variety of physical problems, and are specified by a natural frequency, a periodic forcing function, and a forcing strength parameter. This classic problem is generalized here in two ways: [A] to Random Hill's equations which allow the forcing strength q_k, the oscillation frequency \lambda_k, and the period \tau_k of the forcing function to vary from cycle to cycle, and [B] to Stochastic Hill's equations which contain (at least) one additional term that is a stochastic process \xi. This paper considers both random and stochastic Hill's equations with small parameter variations, so that p_k=q_k-<q_k>, \ell_k=\lambda_k-<\lambda_k>, and \xi are all O(\epsilon), where \epsilon<<1. We show that random Hill's equations and stochastic Hill's equations have the same growth rates when the parameter variations p_k and \ell_k obey certain constraints given in terms of the moments of \xi. For random Hill's equations, the growth rates for the solutions are given by the growth rates of a matrix transformation, under matrix multiplication, where the matrix elements vary from cycle to cycle. Unlike classic Hill's equations where the parameter space (the \lambda-q plane) displays bands of stable solutions interlaced with bands of unstable solutions, random Hill's equations are generically unstable. We find analytic approximations for the growth rates of the instability; for the regime where Hill's equation is classically stable, and the parameter variations are small, the growth rate \gamma = O(\epsilon^2). Using the relationship between the (\ell_k,p_k) and the \xi, this result for \gamma can be used to find growth rates for stochastic Hill's equations.

Observations indicate that intermediate mass stars, binary stars, and stellar remnants often host planets; a complete explanation of these systems requires an understanding of how planetary orbits evolve as their central stars lose mass. Motivated by these dynamical systems, this paper generalizes in two directions previous studies of orbital evolution in planetary systems with stellar mass loss: [1] Many previous treatments focus on constant mass loss rates and much of this work is carried out numerically. Here we study a class of single planet systems where the stellar mass loss rate is time dependent. The mass loss rate can be increasing or decreasing, but the stellar mass always decreases monotonically. For this class of models, we develop analytic approximations to specify the final orbital elements for planets that remain bound after the epoch of mass loss, and find the conditions required for the planets to become unbound. We also show that for some mass loss functions, planets become unbound only in the asymptotic limit where the stellar mass vanishes. [2] We consider the chaotic evolution for two planet systems with stellar mass loss. Here we focus on a model consisting of analogs of Jupiter, Saturn, and the Sun. By monitoring the divergence of initially similar trajectories through time, we calculate the Lyapunov exponents of the system. This analog solar system is chaotic in the absence of mass loss with Lyapunov time in the range 5 - 10 Myr; we find that the Lyapunov time decreases with increasing stellar mass loss rate, with a nearly linear relationship between the two time scales. Taken together, the results of this paper help provide an explanation for a wide range of dynamical evolution that occurs in solar systems with stellar mass loss.

Some theories of star formation suggest massive stars may only form in clustered environments, which would create a deficit of massive stars in low density environments. Observationally, Massey (2002) finds such a deficit in samples of the field population in the Small and Large Magellanic Clouds, with an IMF slope of {\Gamma} ~ 4. These IMF measurements represent some of the largest known deviations from the standard Salpeter IMF slope of {\Gamma}=1.35. Here, we carry out a comprehensive investigation of the mass function above 20 solar masses for the entire field population of the Small Magellanic Cloud, based on data from the Runaways and Isolated O Type Star Spectroscopic Survey of the SMC (RIOTS4). This is a spatially complete census of the entire field OB star population of the SMC obtained with the IMACS multi-object spectrograph and MIKE echelle spectrograph on the Magellan telescopes. Based on Monte-Carlo simulations of the evolved present-day mass function, we find the slope of the field IMF above 20 solar masses is {\Gamma}=2.3+/-0.4. We extend our IMF measurement to lower masses using BV photometry from the OGLE II survey. We use a statistical approach to generate a probability distribution for the mass of each star from the OGLE photometry, and we again find {\Gamma}=2.3+/-0.6 for stellar masses from 7-20 solar masses. The discovery and removal of ten runaways in our RIOTS4 sample steepens the field IMF slope to {\Gamma}=2.8+/-0.5. We discuss the possible effects of binarity and star-formation history on our results, and conclude that the steep field massive star IMF is most likely a real effect.

Motivated by the large number of extrasolar planetary systems that are near mean motion resonances, this paper explores a related type of dynamical behavior known as "nodding". Here, the resonance angle of a planetary system executes libration (oscillatory motion) for several cycles, circulates for one or more cycles, and then enters once again into libration. This type of complicated dynamics can affect our interpretation of observed planetary systems that are in or near mean motion resonance. This work shows that planetary systems in (near) mean motion resonance can exhibit nodding behavior, and outlines the portion of parameter space where it occurs. This problem is addressed using both full numerical integrations of the planetary systems and via model equations obtained through expansions of the disturbing function. In the latter approach, we identify the relevant terms that allow for nodding. The two approaches are in agreement, and show that nodding often occurs when a small body is in an external mean motion resonance with a larger planet. As a result, the nodding phenomenon can be important for interpreting observations of transit timing variations, where the existence of smaller bodies is inferred through their effects on larger, observed transiting planets. For example, in actively nodding planetary systems, both the amplitude and frequency of the transit timing variations depend on the observational time window.

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.

Observed Hot Jupiters exhibit a wide range of physical properties. For a given mass, many planets have inflated radii, while others are surprisingly compact and may harbor large central cores. Motivated by the observational sample, this paper considers possible effects from collisions of smaller rocky planets with gas giant planets. In this scenario, the Jovian planets migrate first and enter into (approximately) 4 day orbits, whereas rocky planets (mass = 0.1-20 that of Earth) migrate later and then encounter the gaseous giants. Previous work indicates that the collision rates are high for such systems. This paper calculates the trajectories of incoming rocky planets as they orbit within the gaseous planets and are subjected to gravitational, frictional, and tidal forces. These collisions always increase the metallicity of the Jovian planets. If the incoming rocky bodies survive tidal destruction and reach the central regions, they provide a means of producing large planetary cores. Both the added metallicity and larger cores act to decrease the radii of the gas giants at fixed mass. The energy released during these collisions provides the Jovian planet with an additional heat source; here we determine the radial layers where kinetic energy of the colliding body is dissipated, including the energy remaining upon impact with the existing core. This process could have long-term effects if the colliding body deposits significant energy deep in the interior, in regions of high opacity. Both Hot Jupiters and newly formed gas giants have inflated radii, large enough to allow incoming rocky planets to survive tidal disruption, enhance the central core mass, and deposit significant energy (in contrast, denser giant planets with the mass and radius of Jupiter are expected to tidally destroy incoming rocky bodies).

Young stars typically have strong magnetic fields, so that the magnetospheres of newly formed close binaries can interact, dissipate energy, and produce synchrotron radiation. The V773 Tau A binary system, a pair of T Tauri stars with a 51 day orbit, displays such a signature, with peak emission taking place near periastron. This paper proposes that the observed emission arises from the change in energy stored in the composite magnetic field of the system. We model the fields using the leading order (dipole) components and show that this picture is consistent with current observations. In this model, the observed radiation accounts for a fraction of the available energy of interaction between the magnetic fields from the two stars. Assuming antisymmetry, we compute the interaction energy $E_{\rm int}$ as a function of the stellar radii, the stellar magnetic field strengths, the binary semi-major axis, and orbital eccentricity, all of which can be measured independently of the synchrotron radiation. The variability in time and energetics of the synchrotron radiation depend on the details of the annihilation of magnetic fields through reconnection events, which generate electric fields that accelerate charged particles, and how those charged particles, especially fast electrons, are removed from the interaction region. However, the major qualitative features are well described by the background changes in the global magnetic configuration driven by the orbital motion. The theory can be tested by observing a collection of pre-main-sequence binary systems.

(abridged) Accretion from disks onto young stars is thought to follow magnetic field lines from the inner disk edge to the stellar surface. The accretion flow thus depends on the geometry of the magnetic field. This paper extends previous work by constructing a collection of orthogonal coordinate systems, including the corresponding differential operators, where one coordinate traces the magnetic field lines. This formalism allows for an (essentially) analytic description of the geometry and the conditions required for the flow to pass through sonic points. Using this approach, we revisit the problem of magnetically controlled accretion flow in a dipole geometry, and then generalize the treatment to consider magnetic fields with multiple components, including dipole, octupole, and split monopole contributions. This approach can be generalized further to consider more complex magnetic field configurations. Observations indicate that accreting young stars have substantial dipole and octupole components, and that accretion flow is transonic. If the effective equation of state for the fluid is too stiff, the flow cannot pass smoothly through the sonic points in steady state. For a multipole field of order \ell, we derive a constraint on the polytropic index, n>\ell+3/2, required for steady transonic flow to reach free-fall velocities. For octupole fields, inferred on surfaces of T Tauri stars, n>9/2, so that the flow must be close to isothermal. The inclusion of octupole field components produces higher densities at the stellar surface and smaller hot spots, which occur at higher latitudes; the magnetic truncation radius is also modified. This contribution thus increases our understanding of magnetically controlled accretion for young stellar objects and can be applied to a variety of additional astrophysical problems.

The observed population of Hot Jupiters displays a stunning variety of physical properties, including a wide range of densities and core sizes for a given planetary mass. Motivated by the observational sample, this paper studies the accretion of rocky planets by Hot Jupiters, after the Jovian planets have finished their principal migration epoch and become parked in $\sim4$-day orbits. In this scenario, rocky planets form later and then migrate inward due to torques from the remaining circumstellar disk, which also damps the orbital eccentricity. This mechanism thus represents one possible channel for increasing the core masses and metallicities of Hot Jupiters. This paper determines probabilities for the possible end states for the rocky planet: collisions with the Jovian planets, accretion onto the star, ejection from the system, and long-term survival of both planets. These probabilities depend on the mass of the Jovian planet and its starting orbital eccentricity, as well as the eccentricity damping rate for the rocky planet. Since these systems are highly chaotic, a large ensemble ($N\sim10^3$) of simulations with effectively equivalent starting conditions is required. Planetary collisions are common when the eccentricity damping rate is sufficiently low, but are rare otherwise. For systems that experience planetary collisions, this work determines the distributions of impact velocities -- both speeds and impact parameters -- for the collisions. These velocity distributions help determine the consequences of the impacts, e.g., where energy and heavy elements are deposited within the giant planets.

We present a systematic evaluation of the agreement between the observed radii of 90 well-characterized transiting extrasolar giant planets and their corresponding model radii. Our model radii are drawn from previously published calculations of core-less giant planets that have attained their asymptotic radii, and which have been tabulated for a range of planet masses and equilibrium temperatures. (We report a two-dimensional polynomial fitting function that accurately represents the models). As expected, the model radii provide a statistically significant improvement over a null hypothesis that the sizes of giant planets are completely independent of mass and effective temperature. As is well known, however, fiducial models provide an insufficient explanation; the planetary radius anomalies are strongly correlated with planetary equilibrium temperature. We find that the radius anomalies have a best-fit dependence, ${\cal R}\propto T_{\rm eff}^{\alpha}$, with $\alpha=1.4\pm0.6$. Incorporating this relation into the model radii leads to substantially less scatter in the radius correlation. The extra temperature dependence represents an important constraint on theoretical models for Hot Jupiters. Using simple scaling arguments, we find support for the hypothesis of Batygin and Stevenson (2010) that this correlation can be attributed to a planetary heating mechanism that is mediated by magnetohydrodynamic coupling between the planetary magnetic field and near-surface flow that is accompanied by ohmic dissipation at adiabatic depth. Additionally, we find that the temperature dependence is likely too strong to admit kinetic heating as the primary source of anomalous energy generation within the majority of the observed transiting planets.

Recent observations that indicate that some extrasolar planets observed in transit can experience mass loss from their surfaces. Motivated by these findings, this paper considers outflows from Hot Jupiters in the regime where the flow is controlled by magnetic fields. Given the mass loss rates estimated from current observations --- and from theoretical arguments --- magnetic fields will dominate the flow provided that field strength near the planet is greater than $\sim1$ gauss, comparable to the surface fields of the Sun and Jupiter. The problem can be separated into an inner regime, near the planet, where the outflow is launched, and an outer regime where the flow follows (primarily) stellar field lines and interacts with the stellar wind. This paper concentrates on the flow in the inner regime. For a dipole planetary field with a spatially constant background contribution, we construct a set of orthogonal coordinates that follow the field lines and determine the corresponding differential operators. Under the assumption of isothermal flow, we analytically find the conditions required for escaping material to pass smoothly through the sonic transition, and then estimate the mass outflow rates. These magnetically controlled outflows differ significantly from previous spherical models: The outflow rates are somewhat smaller, typically ${\dot M}$ $\sim 10^{9}$ g/s, and the flow is launched primarily from the polar regions of the planet. In addition, if the stellar wind is strong enough, the flow could be reversed and the planet could gain mass from the star.

A large fraction of stars forming in our galaxy are born within clusters embedded in giant molecular clouds. In these environments, the background UV radiation fields impinging upon circumstellar disks can often dominate over the radiation fields produced by each disk's central star. As a result, this background radiation can drive the evaporation of circumstellar disks and lead to the loss of planet forming potential within a cluster. This paper presents a detailed analysis of this process for clusters whose stellar membership falls within the range $100 \le N \le 1000$. For these intermediate-sized clusters, the background UV field is often dominated by the most massive stellar member. Due to the steep slope of the initial mass function, the amount of background UV light that bathes clusters of similar size displays significant variance. As a result, we perform a statistical analysis of this problem by calculating distributions of FUV flux values impinging upon star/disk systems for several cluster scenarios. We find that in the absence of dust attenuation, giant planet formation would likely be inhibited in approximately half of systems forming within intermediate-sized clusters regardless of stellar membership. In contrast, the presence of dust can significantly lower this value, with the effect considerably more pronounced in more populated clusters.

Pairs of migrating extrasolar planets often lock into mean motion resonance as they drift inward. This paper studies the convergent migration of giant planets (driven by a circumstellar disk) and determines the probability that they are captured into mean motion resonance. The probability that such planets enter resonance depends on the type of resonance, the migration rate, the eccentricity damping rate, and the amplitude of the turbulent fluctuations. This problem is studied both through direct integrations of the full 3-body problem, and via semi-analytic model equations. In general, the probability of resonance decreases with increasing migration rate, and with increasing levels of turbulence, but increases with eccentricity damping. Previous work has shown that the distributions of orbital elements (eccentricity and semimajor axis) for observed extrasolar planets can be reproduced by migration models with multiple planets. However, these results depend on resonance locking, and this study shows that entry into -- and maintenance of -- mean motion resonance depends sensitively on migration rate, eccentricity damping, and turbulence.

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.

This paper considers gravitational perturbations in geometrically thin disks with rotation curves dominated by a central object, but with substantial contributions from magnetic pressure and tension. The treatment is general, but the application is to the circumstellar disks that arise during the gravitational collapse phase of star formation. We find the dispersion relation for spiral density waves in these generalized disks and derive the stability criterion for axisymmetric $(m=0)$ disturbances (the analog of the Toomre parameter $Q_T$) for any radial distribution of the mass-to-flux ratio $\lambda$. The magnetic effects work in two opposing directions: on one hand, magnetic tension and pressure stabilize the disk against gravitational collapse and fragmentation; on the other hand, they also lower the rotation rate making the disk more unstable. For disks around young stars the first effect generally dominates, so that magnetic fields allow disks to be stable for higher surface densities and larger total masses. These results indicate that magnetic fields act to suppress the formation of giant planets through gravitational instability. Finally, even if gravitational instability can form a secondary body, it must lose an enormous amount of magnetic flux in order to become a planet; this latter requirement represents an additional constraint for planet formation via gravitational instability and places a lower limit on the electrical resistivity.

This paper derives expressions for the growth rates for the random 2 x 2 matrices that result from solutions to the random Hill's equation. The parameters that appear in Hill's equation include the forcing strength and oscillation frequency. The development of the solutions to this periodic differential equation can be described by a discrete map, where the matrix elements are given by the principal solutions for each cycle. Variations in the forcing strength and oscillation frequency lead to matrix elements that vary from cycle to cycle. This paper presents an analysis of the growth rates including cases where all of the cycles are highly unstable, where some cycles are near the stability border, and where the map would be stable in the absence of fluctuations. For all of these regimes, we provide expressions for the growth rates of the matrices that describe the solutions.

This paper reviews our current understanding of the possible birth environments of our Solar System. Since most stars form within groups and clusters, the question becomes one of determining the nature of the birth aggregate of the Sun. This discussion starts by reviewing Solar System properties that provide constraints on our environmental history. We then outline the range of star-forming environments that are available in the Galaxy, and discuss how they affect star and planet formation. The nature of the solar birth cluster is constrained by many physical considerations, including radiation fields provided by the background environment, dynamical scattering interactions, and by the necessity of producing the short-lived radioactive nuclear species inferred from meteoritic measurements. Working scenarios for the solar birth aggregate can be constructed, as discussed herein, although significant uncertainties remain.

Five new planets orbiting G and K dwarfs have emerged from the Magellan velocity survey. These companions are jovian-mass planets in eccentric (e \geq 0.24) intermediate and long-period orbits. HD 86226b orbits a solar metallicity G2 dwarf. The MP sin i mass of the planet is 1.5 MJUP, the semi-major axis is 2.6 AU, and the eccentricity 0.73. HD 129445b orbits a metal rich G6 dwarf. The minimum mass of the planet is MP sin i =1.6 MJUP, the semi-major axis is 2.9 AU, and the eccentricity 0.70. HD 164604b orbits a K2 dwarf. The MP sin i mass is 2.7 MJUP, semi-major axis is 1.3 AU, and the eccentricity is 0.24. HD 175167b orbits a metal rich G5 star. The MP sin i mass is 7.8 MJUP, the semi-major axis is 2.4 AU, and the eccentricity 0.54. HD 152079b orbits a G6 dwarf. The MP sin i mass of the planet is 3 MJUP, the semi-major axis is 3.2 AU, and the eccentricity is 0.60.

This paper investigates the dynamical evolution of embedded stellar clusters from the protocluster stage, through the embedded star-forming phase, and out to ages of 10 Myr -- after the gas has been removed from the cluster. The relevant dynamical properties of young stellar clusters are explored over a wide range of possible star formation environments using N-body simulations. Many realizations of equivalent initial conditions are used to produce robust statistical descriptions of cluster evolution including the cluster bound fraction, radial probability distributions, as well as the distributions of close encounter distances and velocities. These cluster properties are presented as a function of parameters describing the initial configuration of the cluster, including the initial cluster membership N, initial stellar velocities, cluster radii, star formation efficiency, embedding gas dispersal time, and the degree of primordial mass segregation. The results of this parameter space survey, which includes about 25,000 simulations, provide a statistical description of cluster evolution as a function of the initial conditions. We also present a compilation of the FUV radiation fields provided by these same cluster environments. The output distributions from this study can be combined with other calculations, such as disk photoevaporation models and planetary scattering cross sections, to ascertain the effects of the cluster environment on the processes involved in planet formation.

Magnetic fields are dragged in from the interstellar medium during the gravitational collapse that forms star/disk systems. Consideration of mean field magnetohydrodynamics (MHD) in these disks shows that magnetic effects produce subkeplerian rotation curves and truncate the inner disk. This letter explores the ramifications of these predicted disk properties for the migration of extrasolar planets. Subkeplerian flow in gaseous disks drives a new migration mechanism for embedded planets and modifies the gap opening processes for larger planets. This subkeplerian migration mechanism dominates over Type I migration for sufficiently small planets (m_P < 1 M_\earth) and/or close orbits (r < 1 AU). Although the inclusion of subkeplerian torques shortens the total migration time by only a moderate amount, the mass accreted by migrating planetary cores is significantly reduced. Truncation of the inner disk edge (for typical system parameters) naturally explains final planetary orbits with periods P=4 days. Planets with shorter periods P=2 days can be explained by migration during FU-Ori outbursts, when the mass accretion rate is high and the disk edge moves inward. Finally, the midplane density is greatly increased at the inner truncation point of the disk (the X-point); this enhancement, in conjunction with continuing flow of gas and solids through the region, supports the in situ formation of giant planets.

This paper presents a generalized treatment of Type I planetary migration in the presence of stochastic perturbations. In many planet-forming disks, the Type I migration mechanism, driven by asymmetric torques, acts on a short time scale and compromises planet formation. If the disk also supports MHD instabilities, however, the corresponding turbulent fluctuations produce additional stochastic torques that modify the steady inward migration scenario. This work studies the migration of planetary cores in the presence of stochastic fluctuations using complementary methods, including a Fokker-Planck approach and iterative maps. Stochastic torques have two main effects: [1] Through outward diffusion, a small fraction of the planetary cores can survive in the face of Type I inward migration. [2] For a given starting condition, the result of any particular realization of migration is uncertain, so that results must be described in terms of the distributions of outcomes. In addition to exploring different regimes of parameter space, this paper considers the effects of the outer disk boundary condition, varying initial conditions, and time-dependence of the torque parameters. For disks with finite radii, the fraction of surviving planets decreases exponentially with time. We find the survival fractions and decay rates for a range of disk models, and find the expected distribution of locations for surviving planets. For expected disk properties, the survival fraction lies in the range $0.01 < p_S < 0.1$.

This paper considers random Hill's equations in the limit where the periodic forcing function becomes a Dirac delta function. For this class of equations, the forcing strength $q_k$, the oscillation frequency $\af_k$, and the period are allowed to vary from cycle to cycle. Such equations arise in astrophysical orbital problems in extended mass distributions, in the reheating problem for inflationary cosmologies, and in periodic Schr{\"o}dinger equations. The growth rates for solutions to the periodic differential equation can be described by a matrix transformation, where the matrix elements vary from cycle to cycle. Working in the delta function limit, this paper addresses several coupled issues: We find the growth rates for the $2 \times 2$ matrices that describe the solutions. This analysis is carried out in the limiting regimes of both large $q_k \gg 1$ and small $q_k \ll 1$ forcing strength parameters. For the latter case, we present an alternate treatment of the dynamics in terms of a Fokker-Planck equation, which allows for a comparison of the two approaches. Finally, we elucidate the relationship between the fundamental parameters $(\af_k,q_k)$ appearing in the stochastic differential equation and the matrix elements that specify the corresponding discrete map. This work provides analytic -- and accurate -- expressions for the growth rates of these stochastic differential equations in both the $q_k \gg1 $ and the $q_k \ll 1$ limits.

Motivated by kinematic observations of young embedded clusters, this paper explores possible kinematic signatures produced by asphericity and departures from initial virial equilibrium in these systems. Specifically, the kinematic quantity that is measured and calculated in this study is the distribution of the line-of-sight velocities as a function of position along the cluster. Although clusters are found within a wide range of sizes, we focus on the regime with stellar membership $N \sim 10^3$. The gravitational potential of these young clusters is dominated by the gas, and the geometry of the gas distribution is generalized to include axisymmetric (and triaxial) forms. With this loss of symmetry, the kinematic results thus depend on viewing angle. This work also considers the stars to begin their trajectories with subvirial speeds, as indicated by observations of core motions in such clusters. Our results determine the conditions necessary for the kinematic signature to display interesting structure, i.e., a non-spherical potential, a viewing angle that is not along one of the principal axes, and subvirial starting conditions. We characterize the effects on this signature due to projection angle, initial stellar velocities, cluster elongation, and star formation efficiency. Finally, we compare our theoretical results to recent kinematic observations of the Orion Nebula Cluster; we find that the observations can be explained provided that the cluster is non-spherical, starts with subvirial initial velocities, and is not viewed along a principal axis.

Most stars are born in rich young stellar clusters (YSCs) embedded in giant molecular clouds. The most massive stars live out their short lives there, profoundly influencing their natal environments by ionizing HII regions, inflating wind-blown bubbles, and soon exploding as supernovae. Thousands of lower-mass pre-main sequence stars accompany the massive stars, and the expanding HII regions paradoxically trigger new star formation as they destroy their natal clouds. While this schematic picture is established, our understanding of the complex astrophysical processes involved in clustered star formation have only just begun to be elucidated. The technologies are challenging, requiring both high spatial resolution and wide fields at wavelengths that penetrate obscuring molecular material and remove contaminating Galactic field stars. We outline several important projects for the coming decade: the IMFs and structures of YSCs; triggered star formation around YSC; the fate of OB winds; the stellar populations of Infrared Dark Clouds; the most massive star clusters in the Galaxy; tracing star formation throughout the Galactic Disk; the Galactic Center region and YSCs in the Magellanic Clouds. Programmatic recommendations include: developing a 30m-class adaptive optics infrared telescope; support for high-resolution and wide field X-ray telescopes; large-aperture sub-millimeter and far-infrared telescopes; multi-object infrared spectrographs; and both numerical and analytical theory.

We report low mass companions orbiting five Solar-type stars that have emerged from the Magellan precision Doppler velocity survey, with minimum (Msini) masses ranging from 1.2 to 25 Mjup. These nearby target stars range from mildly metal-poor to metal-rich, and appear to have low chromospheric activity. The companions to the brightest two of these stars have previously been reported from the CORALIE survey. Four of these companions (HD 48265-b, HD 143361-b, HD 28185-b, HD 111232-b) are low-mass Jupiter-like planets in eccentric intermediate and long-period orbits. On the other hand, the companion to HD 43848 appears to be a long period brown dwarf in a very eccentric orbit.

This paper investigates anthropic bounds on the vacuum energy $\Lambda$. We first consider the possibility of cosmic observers existing at any random time (including the future) for constant $\Lambda$, and take into account the suppression of new structure formation as the universe approaches its eternal DeSitter (DS) geometry. Structures that collapse prior to the era of $\Lambda$-domination will lose causal contact with our Hubble volume within a finite (short) conformal time $\tau_{\ast}$. Any remnants within our Hubble volume then suffer a cosmological heat death after the universe becomes DS. The probability for finding observers by random measurements in the volume bound by the DeSitter horizon is $P \simeq 0$, and it is a simple consequence of the information loss problem for eternal DS spaces resulting from the finite and constant value of its temperature $T_{DS} \simeq \Lambda^{-1/2}$ and entropy $S = 3/(G\Lambda)$. By contrast, for geometries with $\Lambda = 0$, structures can condense and entropy production can continue without bounds at any epoch. The probability of finding observers in $\Lambda = 0$ geometries is thus overwhelming higher than in DS spaces. As a result, anthropic reasoning does not explain the small but nonzero vacuum energy observed in our universe. We also address the case where observers are considered only at a specially chosen time -- like the present epoch -- but relax the allowed values of starting density fluctuations and hence the redshift of galaxy formation. In this latter case, the bounds on a $\Lambda$ can be millions of times larger than previous estimates -- and the observed value. We thus conclude that anthropic reasoning has limited predictive power.

This paper continues previous work on the effects of turbulence on mean motion resonances in extrasolar planetary systems. Turbulence is expected to arise in the disks that form planets, and these fluctuations act to compromise resonant configurations. This paper extends previous work by considering how interactions between the planets and possible damping effects imposed by the disk affect the outcomes. These physical processes are studied using three approaches: numerical integrations of the 3-body problem with additional forcing due to turbulence, model equations that reduce the problem to stochastically driven oscillators, and Fokker-Planck equations that describe the time evolution of an ensemble of systems. With this combined approach, we elucidate the physics of how turbulence can remove extrasolar planetary systems from mean motion resonance. As expected, systems with sufficiently large damping (dissipation) can maintain resonance, in spite of turbulent forcing. In the absence of strong damping, ensembles of these systems exhibit two regimes of behavior, where the fraction of the bound states decreases as a power-law or as an exponential. Both types of behavior can be understood through the model developed herein. For systems with weak interactions between planets, the model reduces to a stochastic pendulum, and the fraction of bound states decreases as a power-law. For highly interactive systems, the dynamics are more complicated and the fraction of bound states decreases exponentially. We show how planetary interactions lead to drift terms in the Fokker-Planck equation and account for this exponential behavior. In addition to clarifying the physical processes involved, this paper strengthens the finding that turbulence implies that mean motions resonances should be rare.

We report the detection of the first two planets from the N2K Doppler planet search program at the Magellan telescopes. The first planet has a mass of M sin i = 4.96 M_Jup and is orbiting the G3 IV star HD154672 with an orbital period of 163.9 days. The second planet is orbiting the F7 V star HD205739 with an orbital period of 279.8 days and has a mass of M sin i = 1.37 M_Jup. Both planets are in eccentric orbits, with eccentricities e = 0.61 and e = 0.27, respectively. Both stars are metal rich and appear to be chromospherically inactive, based on inspection of their Ca II H and K lines. Finally, the best Keplerian model fit to HD205739b shows a trend of 0.0649 m/s/day, suggesting the presence of an additional outer body in that system.

Motivated by the possibility that the fundamental ``constants'' of nature could vary with time, this paper considers the long term evolution of white dwarf stars under the combined action of proton decay and variations in the gravitational constant. White dwarfs are thus used as a theoretical laboratory to study the effects of possible time variations, especially their implications for the future history of the universe. More specifically, we consider the gravitational constant $G$ to vary according to the parametric relation $G = G_0 (1 + t/t_\ast)^{-p}$, where the time scale $t_\ast$ is the same order as the proton lifetime. We then study the long term fate and evolution of white dwarf stars. This treatment begins when proton decay dominates the stellar luminosity, and ends when the star becomes optically thin to its internal radiation.

Motivated by the possible existence of other universes, with possible variations in the laws of physics, this paper explores the parameter space of fundamental constants that allows for the existence of stars. To make this problem tractable, we develop a semi-analytical stellar structure model that allows for physical understanding of these stars with unconventional parameters, as well as a means to survey the relevant parameter space. In this work, the most important quantities that determine stellar properties -- and are allowed to vary -- are the gravitational constant $G$, the fine structure constant $\alpha$, and a composite parameter $C$ that determines nuclear reaction rates. Working within this model, we delineate the portion of parameter space that allows for the existence of stars. Our main finding is that a sizable fraction of the parameter space (roughly one fourth) provides the values necessary for stellar objects to operate through sustained nuclear fusion. As a result, the set of parameters necessary to support stars are not particularly rare. In addition, we briefly consider the possibility that unconventional stars (e.g., black holes, dark matter stars) play the role filled by stars in our universe and constrain the allowed parameter space.

This paper considers the effects of turbulence on mean motion resonances in extrasolar planetary systems and predicts that systems rarely survive in a resonant configuration. A growing number of systems are reported to be in resonance, which is thought to arise from the planet migration process. If planets are brought together and moved inward through torques produced by circumstellar disks, then disk turbulence can act to prevent planets from staying in a resonant configuration. This paper studies this process through numerical simulations and via analytic model equations, where both approaches include stochastic forcing terms due to turbulence. We explore how the amplitude and forcing time intervals of the turbulence affect the maintenance of mean motion resonances. If turbulence is common in circumstellar disks during the epoch of planet migration, with the amplitudes indicated by current MHD simulations, then planetary systems that remain deep in mean motion resonance are predicted to be rare. More specifically, the fraction of resonant systems that survive over a typical disk lifetime of 1 Myr is of order 0.01. If mean motion resonances are found to be common, their existence would place tight constraints on the amplitude and duty cycle of turbulent fluctuations in circumstellar disks. These results can be combined by expressing the expected fraction of surviving resonant systems in the approximate form P_b = C / N_{orb}^{1/2}, where the dimensionless parameter C = 10 - 50 and where N_{orb} is the number of orbits for which turbulence is active.

We propose a model for early universe cosmology without the need for fundamental scalar fields. Cosmic acceleration and phenomenologically viable reheating of the universe results from a series of energy transitions, where during each transition vacuum energy is converted to thermal radiation. We show that this `cascading universe' can lead to successful generation of adiabatic density fluctuations and an observable gravity wave spectrum in some cases, where in the simplest case it reproduces a spectrum similar to slow-roll models of inflation. We also find the model provides a reasonable reheating temperature after inflation ends. This type of model may also be relevant for addressing the smallness of the vacuum energy today.

The rosette-shaped motion of a particle in a central force field is known to be classically solvable by quadratures. We present a new approach of describing and characterizing such motion based on the eccentricity vector of the two body problem. In general, this vector is not an integral of motion. However, the orbital motion, when viewed from the nonuniformly rotating frame defined by the orientation of the eccentricity vector, can be solved analytically and will either be a closed periodic circulation or libration. The motion with respect to inertial space is then given by integrating the argument of periapsis with respect to time. Finally we will apply the decomposition to a modern central potential, the spherical Hernquist-Newton potential, which models dark matter halos of galaxies with central black holes.

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.

Motivated by a class of orbit problems in astrophysics, this paper considers solutions to Hill's equation with forcing strength parameters that vary from cycle to cycle. The results are generalized to include period variations from cycle to cycle. The development of the solutions to the differential equation is governed by a discrete map. For the general case of Hill's equation in the unstable limit, we consider separately the case of purely positive matrix elements and those with mixed signs; we then find exact expressions, bounds, and working estimates for the growth rates. We also find exact expressions, estimates, and bounds for the infinite products of several $2 \times 2$ matrices with random variables in the matrix elements. In the limit of sharply spiked forcing terms (the delta function limit), we find analytic solutions for each cycle and for the discrete map that matches solutions from cycle to cycle; for this case we find the growth rates and the condition for instability in the limit of large forcing strength, as well as the widths of the stable/unstable zones.

This paper re-examines the problem of ambipolar diffusion as a mechanism for the production and runaway evolution of centrally condensed molecular cloud cores, a process that has been termed the gravomagneto catastrophe. Our calculation applies in the geometric limit of a highly flattened core and allows for a semi-analytic treatment of the full problem, although physical fixes are required to resolve a poor representation of the central region. A noteworthy feature of the overall formulation is that the solutions for the ambipolar diffusion portion of the evolution for negative times ($t < 0$) match smoothly onto the collapse solutions for positive times ($t > 0$). The treatment shows that the resulting cores display non-zero, but sub-magnetosonic, inward velocities at the end of the diffusion epoch, in agreement with current observations. Another important result is the derivation of an analytic relationship between the dimensionless mass to flux ratio $\lambda_0\equiv f_0^{-1}$ of the central regions produced by runaway core condensation and the dimensionless measure of the rate of ambipolar diffusion $\epsilon$. In conjunction with previous work showing that ambipolar diffusion takes place more quickly in the presence of turbulent fluctuations, i.e., that the effective value of $\epsilon$ can be enhanced by turbulence, the resultant theory provides a viable working hypothesis for the formation of isolated molecular-cloud cores and their subsequent collapse to form stars and planetary systems.

This paper constructs an analytic form for a triaxial potential that describes the dynamics of a wide variety of astrophysical systems, including the inner portions of dark matter halos, the central regions of galactic bulges, and young embedded star clusters. Specifically, this potential results from a density profile of the form $\rho (m) \propto m^{-1}$, where the radial coordinate is generalized to triaxial form so that $m^2 = x^2/a^2 + y^2/b^2 + z^2/c^2 $. Using the resulting analytic form of the potential, and the corresponding force laws, we construct orbit solutions and show that a robust orbit instability exists in these systems. For orbits initially confined to any of the three principal planes, the motion in the perpendicular direction can be unstable. We discuss the range of parameter space for which these orbits are unstable, find the growth rates and saturation levels of the instability, and develop a set of analytic model equations that elucidate the essential physics of the instability mechanism. This orbit instability has a large number of astrophysical implications and applications, including understanding the formation of dark matter halos, the structure of galactic bulges, the survival of tidal streams, and the early evolution of embedded star clusters.

The extrasolar planets discovered to date possess unexpected orbital elements. Most orbit their host stars with larger eccentricities and smaller semi-major axes than similarly sized planets in our own solar system do. It is generally agreed that the interaction between giant planets and circumstellar disks (Type II migration) drives these planets inward to small radii, but the effect of these same disks on orbital eccentricity, e, is controversial. Several recent analytic calculations suggest that disk-planet interactions can excite eccentricity, while numerical studies generally produce eccentricity damping. This paper addresses this controversy using a quasi-analytic approach, drawing on several preceding analytic studies. This work refines the current treatment of eccentricity evolution by removing several approximations from the calculation of disk torques. We encounter neither uniform damping nor uniform excitation of orbital eccentricity, but rather a function de/dt that varies in both sign and magnitude depending on eccentricity and other solar system properties. Most significantly, we find that for every combination of disk and planet properties investigated herein, corotation torques produce negative values of de/dt for some range in e within the interval [0.1, 0.5]. If corotation torques are saturated, this region of eccentricity damping disappears, and excitation occurs on a short timescale of less than 0.08 Myr. Thus, our study does not produce eccentricity excitation on a timescale of a few Myr -- we obtain either eccentricity excitation on a short time scale, or eccentricity damping on a longer time scale. Finally, we discuss the implications of this result for producing the observed range in extrasolar planet eccentricity.

(Abridged) The discovery of extrasolar planetary systems revealed an unexpected diversity of planetary systems that has revolutionized planet formation theory. A strong program of theoretical research is essential to maximize both the discovery potential and the scientific returns of future observational programs, so as to achieve a deeper understanding of the formation and evolution of planetary systems. We outline three broad categories of theoretical research: detailed studies of specific planetary systems, testing planet formation models by comparing their predictions to the observed exoplanet population, and detailed modeling of specific physical processes. We describe how such theoretical research plays an important role in analyzing observations for a wide range detection methods and contributes to understanding the Earth's place in the universe and the potential for Earth-like life beyond our solar system. In this white paper, we suggest how to maximize the scientific return of future exoplanet observations. Our recommendations include a strong theory program, support for multiple observational programs that will study a diverse set of planets and stars, significant observing time devoted to follow-up observations, and healthy collaboration between observers and theorists.

During the last three decades, evidence has mounted that star and planet formation is not an isolated process, but is influenced by current and previous generations of stars. Although cool stars form in a range of environments, from isolated globules to rich embedded clusters, the influences of other stars on cool star and planet formation may be most significant in embedded clusters, where hundreds to thousands of cool stars form in close proximity to OB stars. At the cool stars 14 meeting, a splinter session was convened to discuss the role of environment in the formation of cool stars and planetary systems; with an emphasis on the ``hot'' environment found in rich clusters. We review here the basic results, ideas and questions presented at the session. We have organized this contribution into five basic questions: what is the typical environment of cool star formation, what role do hot star play in cool star formation, what role does environment play in planet formation, what is the role of hot star winds and supernovae, and what was the formation environment of the Sun? The intention is to review progress made in addressing each question, and to underscore areas of agreement and contention.

We calculate herein the late stages of terrestrial planet accumulation around a solar type star that has a binary companion with semimajor axis larger than the terrestrial planet region. We perform more than one hundred simulations to survey binary parameter space and to account for sensitive dependence on initial conditions in these dynamical systems. As expected, sufficiently wide binaries leave the planet formation process largely unaffected. As a rough approximation, binary stars with periastron $q_B > 10$ AU have minimal effect on terrestrial planet formation within $\sim 2$ AU of the primary, whereas binary stars with $q_B \la$ 5 AU restrict terrestrial planet formation to within $\sim$ 1 AU of the primary star. Given the observed distribution of binary orbital elements for solar type primaries, we estimate that about 40 -- 50 percent of the binary population is wide enough to allow terrestrial planet formation to take place unimpeded. The large number of simulations allows for us to determine the distribution of results -- the distribution of plausible terrestrial planet systems -- for effectively equivalent starting conditions. We present (rough) distributions for the number of planets, their masses, and their orbital elements.

We explore the effects of small scale structure on the formation and equilibrium of dark matter halos in a universe dominated by vacuum energy. We present the results of a suite of four N-body simulations, two with a LCDM initial power spectrum and two with WDM-like spectra that suppress the early formation of small structures. All simulations are run into to far future when the universe is 64Gyr/h old, long enough for halos to essentially reach dynamical equilibrium. We quantify the importance of hierarchical merging on the halo mass accretion history, the substructure population, and the equilibrium density profile. We modify the mass accretion history function of Wechsler et al. (2002) by introducing a parameter, \gamma, that controls the rate of mass accretion, dln(M) / dln(a) ~ a^(-\gamma), and find that this form characterizes both hierarchical and monolithic formation. Subhalo decay rates are exponential in time with a much shorter time scale for WDM halos. At the end of the simulations, we find truncated Hernquist density profiles for halos in both the CDM and WDM cosmologies. There is a systematic shift to lower concentration for WDM halos, but both cosmologies lie on the same locus relating concentration and formation epoch. Because the form of the density profile remains unchanged, our results indicate that the equilibrium halo density profile is set independently of the halo formation process.

This paper considers general relativistic (GR) effects in currently observed extrasolar planetary systems. Although GR corrections are small, they can compete with secular interactions in these systems and thereby play an important role. Specifically, some of the observed multiple planet systems are close to secular resonance, where the dynamics is extremely sensitive to GR corrections, and these systems can be used as laboratories to test general relativity. For the three-planet solar system Upsilon Andromedae, secular interaction theory implies an 80% probability of finding the system with its observed orbital elements if GR is correct, compared with only a 2% probability in the absence of GR. In the future, tighter constraints can be obtained with increased temporal coverage.

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).

This paper constructs an analytic framework for calculating the assembly of galactic disks from the collapse of gas within dark matter halos, with the goal of determining the surface density profiles. Gas parcels (baryons) fall through the potentials of dark matter halos on nearly ballistic, zero energy orbits and collect in a rotating disk. The dark matter halos have a nearly universal form, as determined previously through numerical simulations. The calculation is first carried out for a variety of pre-collapse mass distributions and rotation profiles, including polytropic spheres in hydrostatic equilibrium with the halo potential. The resulting disk surface density profiles have nearly power-law forms, with well-defined edges. This idealized scenario is generalized to include non-spherical starting states and multiple accretion events (due to gas being added to the halo via merger events). This latter complication is explored in detail and considers a log-normal distribution for the angular momenta of the pre-collapse states. If this distribution is wide, then the composite surface density approaches a universal power-law form, independent of the shape of the constituent profiles. When the angular momentum distribution has an intermediate width, the composite surface density attains a nearly exponential form, roughly consistent with profiles of observed galaxies.

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.

This paper studies the effects of dynamical interactions among the planets in observed extrasolar planetary systems, including hypothetical additional bodies, with a focus on secular perturbations. These interactions cause the eccentricities of the planets to explore a distribution of values over time scales that are long compared to observational time baselines, but short compared to the age of the systems. The same formalism determines the eccentricity forcing of hypothetical test bodies (terrestrial planets) in these systems and we find which systems allow for potentially habitable planets. Such planets would be driven to nonzero orbital eccentricity and we derive the distribution of stellar flux experienced by the planets over the course of their orbits. The general relativistic corrections to secular interaction theory are included in the analysis and such effects are important in systems with close planets ($\sim$4 day orbits). Some extrasolar planetary systems (e.g., Upsilon Andromedae) can be used as a test of general relativity, whereas in other systems, general relativity can be used to constrain the system parameters (e.g., $\sin i \gta 0.93$ for HD160691). For the case of hot Jupiters, we discuss how the absence of observed eccentricity implies the absence of companion planets.

This paper studies the long term evolution of planetary systems containing short-period planets, including the effects of tidal circularization, secular excitation of eccentricity by companion planets, and stellar damping. For planetary systems subject to all of these effects, analytic solutions (or approximations) are presented for the time evolution of the semi-major axes and eccentricities. Secular interactions enhance the inward migration and accretion of hot Jupiters, while general relativity tends to act in opposition by reducing the effectiveness of the secular perturbations. The analytic solutions presented herein allow us to understand these effects over a wide range of parameter space and to isolate the effects of general relativity in these planetary systems.

It has been claimed recently that the black hole information-loss paradox has been resolved: the evolution of quantum states in the presence of a black hole is unitary and information preserving. We point out that, contrary to some claims in literature, information-preserving black holes still violate baryon number and any other quantum number which follows from an effective (and thus approximate) or anomalous symmetry.

We examine the properties of embedded clusters within 1 kiloparsec using new data from the Spitzer Space Telescope, as well as recent results from 2MASS and other ground-based near-infrared surveys. We use surveys of entire molecular clouds to understand the range and distribution of cluster membership, size and surface density. The Spitzer data demonstrate clearly that there is a continuum of star- forming environments, from relative isolation to dense clusters. The number of members of a cluster is correlated with the cluster radius, such that the average surface density of clusters having a few to a thousand members varies by a factor of only a few. The spatial distributions of Spitzer-identified young stellar objects frequently show elongation, low density halos, and sub-clustering. The spatial distributions of protostars resemble the distribution of dense molecular gas, suggesting that their morphologies result directly from the fragmentation of the natal gas. We also examine the effects of the cluster environments on star and planet formation. Although Far-UV and Extreme-UV radiation from massive stars can truncate disks in a few million years, fewer than half of the young stars in our sample (embedded clusters within 1 kpc) are found in regions of strong FUV and EUV fields. Typical volume densities and lifetimes of the observed clusters suggest that dynamical interactions are not an important mechanism for truncating disks on solar system size scales.

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.

This paper presents a parametric study of giant planet migration through the combined action of disk torques and planet-planet scattering. The torques exerted on planets during Type II migration in circumstellar disks readily decrease the semi-major axes, whereas scattering between planets increases the orbital eccentricities. This paper presents a parametric exploration of the possible parameter space for this migration scenario using two (initial) planetary mass distributions and a range of values for the time scale of eccentricity damping (due to the disk). For each class of systems, many realizations of the simulations are performed in order to determine the distributions of the resulting orbital elements of the surviving planets; this paper presents the results of 8500 numerical experiments. Our goal is to study the physics of this particular migration mechanism and to test it against observations of extrasolar planets. The action of disk torques and planet-planet scattering results in a distribution of final orbital elements that fills the a-e plane, in rough agreement with the orbital elements of observed extrasolar planets. In addition to specifying the orbital elements, we characterize this migration mechanism by finding the percentages of ejected and accreted planets, the number of collisions, the dependence of outcomes on planetary masses, the time spent in 2:1 and 3:1 resonances, and the effects of the planetary IMF. We also determine the distribution of inclination angles of surviving planets and the distribution of ejection speeds for exiled planets.

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$).

This paper considers the lithopanspermia hypothesis in star forming groups and clusters, where the chances of biological material spreading from one solar system to another is greatly enhanced (relative to the field) due to the close proximity of the systems and lower relative velocities. These effects more than compensate for the reduced time spent in such crowded environments. This paper uses 300,000 Monte Carlo scattering calculations to determine the cross sections for rocks to be captured by binaries and provides fitting formulae for other applications. We assess the odds of transfer as a function of the ejection speed and number of members in the birth aggregate. The odds of any given ejected meteroid being recaptured by another solar system are relatively low. Because the number of ejected rocks per system can be large, virtually all solar systems are likely to share rocky ejecta with all of the other solar systems in their birth cluster. The number of ejected rocks that carry living microorganisms is much smaller and less certain, but we estimate that several million rocks can be ejected from a biologically active solar system. For typical birth environments, the capture of life bearing rocks is expected to occur 10 -- 16,000 times per cluster (under favorable conditions), depending on the ejection speeds. Only a small fraction of the captured rocks impact the surfaces of terrestrial planets, so that only a few lithopanspermia events are expected (per cluster).

This paper explores orbits in extended mass distributions and develops an analytic approximation scheme based on epicycloids (spirograph patterns). We focus on the Hernquist potential which provides a good model for many astrophysical systems, including elliptical galaxies, dark matter halos, and young embedded star clusters. For a given potential, one can readily calculate orbital solutions as a function of energy and angular momentum using numerical methods. In contrast, this paper presents a number of analytic results for the Hernquist potential and proves a series of general constraints showing that orbits have similar properties for any extended mass distribution (including, e.g., the NFW profile). We discuss circular orbits, radial orbits, zero energy orbits, different definitions of eccentricity, analogs of Kepler's law, the definition of orbital elements, and the relation of these orbits to spirograph patterns (epicycloids). Over much of parameter space the orbits can be adequately described (with accuracy better than 10%) using the parametric equations of epicycloids, thereby providing an analytic description of the orbits. As an application of this formal development, we find a solution for the orbit of the Large Magellanic Cloud in the potential of our Galaxy.

We present an inflationary universe model which utilizes two coupled real scalar fields. The inflation field $\phi$ experiences a first order phase transition and its potential dominates the energy density of the Universe during the inflationary epoch. This field $\phi$ is initially trapped in its metastable minimum and must tunnel through a potential barrier to reach the true vacuum. The second auxiliary field $\psi$ couples to the inflaton field and serves as a catalyst to provide an abrupt end to the inflationary epoch; i.e., the $\psi$ field produces a time-dependent nucleation rate for bubbles of true $\phi$ vacuum. In this model, we find that bubbles of true vacuum can indeed percolate and we argue that thermalization of the interiors can more easily take place. The required degree of flatness (i.e., the fine tuning) in the potential of the $\psi$ field is comparable to that of other models which invoke slowly rolling fields. Pseudo Nambu-Goldstone bosons may naturally provide the flat potential for the rolling field.

In the far future of an accelerating LCDM cosmology, the cosmic web of large-scale structure consists of a set of increasingly isolated halos in dynamical equilibrium. We examine the approach of collisionless dark matter to hydrostatic equilibrium using a large N-body simulation evolved to scale factor a = 100, well beyond the vacuum--matter equality epoch, a_eq ~ 0.75, and 53/h Gyr into the future for a concordance model universe (Omega_m ~ 0.3, Omega_Lambda ~ 0.7). The radial phase-space structure of halos -- characterized at a < a_eq by a pair of zero-velocity surfaces that bracket a dynamically active accretion region -- simplifies at a > 10 a_eq when these surfaces merge to create a single zero-velocity surface, clearly defining the halo outer boundary, rhalo, and its enclosed mass, mhalo. This boundary approaches a fixed physical size encompassing a mean interior density ~ 5 times the critical density, similar to the turnaround value in a classical Einstein-deSitter model. We relate mhalo to other scales currently used to define halo mass (m200, mvir, m180b) and find that m200 is approximately half of the total asymptotic cluster mass, while m180b follows the evolution of the inner zero velocity surface for a < 2 but becomes much larger than the total bound mass for a > 3. The radial density profile of all bound halo material is well fit by a truncated Hernquist profile. An NFW profile provides a somewhat better fit interior to r200 but is much too shallow in the range r200 < r < rhalo.

The favored theoretical explanation for giant planet formation -- in both our solar system and others -- is the core accretion model (although it still has some serious difficulties). In this scenario, planetesimals accumulate to build up planetary cores, which then accrete nebular gas. With current opacity estimates for protoplanetary envelopes, this model predicts the formation of Jupiter-mass planets in 2--3 Myr at 5 AU around solar-mass stars, provided that the surface density of solids is enhanced over that of the minimum-mass solar nebula (by a factor of a few). Working within the core-accretion paradigm, this paper presents theoretical calculations which show that the formation of Jupiter-mass planets orbiting M dwarf stars is seriously inhibited at all radial locations (in sharp contrast to solar-type stars). Planet detection programs sensitive to companions of M dwarfs will test this prediction in the near future.

When stars form in small groups (N = 100 - 500 members), their circumstellar disks are exposed to little EUV radiation but a great deal of FUV radiation from massive stars in the group. This paper calculates mass loss rates for circumstellar disks exposed to external FUV radiation. Previous work treated large disks and/or intense radiation fields in which the disk radius exceeds the critical radius (supercritical disks) where the sound speed in the FUV heated layer exceeds the escape speed. This paper shows that significant mass loss still takes place for subcritical systems. Some of the gas extends beyond the disk edge (above the disk surface) to larger distances where the temperature is higher, the escape speed is lower, and an outflow develops. The evaporation rate is a sensitive function of the stellar mass and disk radius, which determine the escape speed, and the external FUV flux, which determines the temperature structure of the flow. Disks around red dwarfs are readily evaporated and shrink to disk radii of 15 AU on short time scales (10 Myr) when exposed to moderate FUV fields with $G_0$ = 3000. Although disks around solar type stars are more durable, these disks shrink to 15 AU in 10 Myr for intense FUV radiation fields with $G_0$ = 30,000; such fields exist in the central 0.7 pc of a cluster with N = 4000 stars. If our solar system formed in the presence of such strong FUV radiation fields, this mechanism could explain why Neptune and Uranus in our solar system are gas poor, whereas Jupiter and Saturn are gas rich. This mechanism for photoevaporation can also limit the production of Kuiper belt objects and can suppress giant planet formation in sufficiently large clusters, such as the Hyades, especially for disks associated with low mass stars.

Many quantum gravity theories imply that the vacuum is filled with virtual black holes. This paper explores the process in which high energy photons interact with virtual black holes and decay into gravitons and photons of lower energy. The effect requires violation (or modification) of Lorentz invariance and implies that high energy photons cannot propagate over arbitrarily large distances. For the standard Planck mass and the likely form for the interaction cross section, this quantum foam limit becomes $\dist < 450$ Mpc $(\egam/10^7 {\rm GeV})^{-5}$. For quantum gravity theories that posit a lower Planck scale, the interaction rate is larger and the limit is stronger. This paper uses extant observations of gamma rays from cosmological sources to constrain this process for varying values of the Planck mass and a range of forms for the interaction cross sections.

This paper examines how type I planet migration is affected by the presence of turbulent density fluctuations in the circumstellar disk. For type I migration, the planet does not clear a gap in the disk and its secular motion is driven by torques generated by the wakes it creates in the surrounding disk fluid. MHD turbulence creates additional density perturbations that gravitationally interact with the planet and can dominate the torques produced by the migration mechanism itself. This paper shows that conventional type I migration can be readily overwhelmed by turbulent perturbations and hence the usual description of type I migration should be modified in locations where the magnetorotational instability is active. In general, the migrating planet does not follow a smooth inward trned, but rather exhibits a random walk through phase space. Our main conclusion is that MHD turbulence will alter the time scales for type I planet migration and -- because of chaos -- requires the time scales to be described by a distribution of values.

Astronomical observations strongly suggest that our universe is now accelerating and contains a substantial admixture of dark vacuum energy. Using numerical simulations to study this newly consolidated cosmological model (with a constant density of dark energy), we show that astronomical structures freeze out in the near future and that the density profiles of dark matter halos approach the same general form. Every dark matter halo grows asymptotically isolated and thereby becomes the center of its own island universe. Each of these isolated regions of space-time approaches a universal geometry and we calculate the corresponding form of the space-time metric.

Current cosmological data indicate that our universe contains a substantial component of dark vacuum energy that is driving the cosmos to accelerate. We examine the immediate and longer term consequences of this dark energy (assumed here to have a constant density). Using analytic calculations and supporting numerical simulations, we present criteria for test bodies to remain bound to existing structures. We show that collapsed halos become spatially isolated and dynamically relax to a particular density profile with logarithmic slope steeper than -3 at radii beyond r_200. The asymptotic form of the space-time metric is then specified. We develop this scenario further by determining the effects of the accelerating expansion on the background radiation fields and individual particles. In an appendix, we generalize these results to include quintessence.

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.

Motivated by the observed correlation between black hole masses $\mbh$ and the velocity dispersion $\sigma$ of host galaxies, we develop a theoretical model of black hole formation in galactic bulges (this paper generalizes an earlier ApJ Letter). The model assumes an initial state specified by a a uniform rotation rate $\Omega$ and a density distribution of the form $\rho = \aeff^2 / 2 \pi G r^2$ (so that $\aeff$ is an effective transport speed). The black hole mass is determined when the centrifugal radius of the collapse flow exceeds the capture radius of the central black hole (for Schwarzschild geometry). This model reproduces the observed correlation between the estimated black hole masses and the velocity dispersions of galactic bulges, i.e., $\mbh \approx 10^8 M_\odot (\sigma/200 {\rm km s^{-1}})^4$, where $\sigma = \sqrt{2} \aeff$. To obtain this normalization, the rotation rate $\Omega \approx 2 \times 10^{15}$ rad/s. The model also defines a bulge mass scale $M_B$. If we identify the scale $M_B$ with the bulge mass, the model determines the ratio $\mrat$ of black hole mass to the host mass: $\mrat$ $\approx$ 0.0024 $(\sigma/200 {\rm km s^{-1}})$, again in reasonable agreement with observed values. In this scenario, supermassive black holes form quickly (in $\sim10^5$ yr) and are born rapidly rotating (with $a/M \sim 0.9$). This paper also shows how these results depend on the assumed initial conditions; the most important quantity is the initial distribution of specific angular momentum in the pre-collapse state.

This paper explores the intermediate-time dynamics of newly formed solar systems with a focus on possible mechanisms for planetary migration. We consider two limiting corners of the available parameter space -- crowded systems containing N=10 giant planets in the outer solar system, and solar systems with N=2 planets that are tidally interacting with a circumstellar disk. For a given set of initial conditions, dynamical relaxation leads to a well-defined distribution of possible solar systems. For each class of initial conditions, we perform large numbers of N-body simulations to obtain a statistical description of the possible outcomes. For N=10 planet systems, we consider several different planetary mass distributions; we also perform secondary sets of simulations to explore chaotic behavior and longer term dynamical evolution. For systems with 10 planets initially populating the radial range 5 - 30 AU, these scattering processes naturally produce planetary orbits with $a\sim1$ AU and the full range of possible eccentricity, but shorter period orbits are difficult to achieve. To account for the observed eccentric giant planets, we also explore a mechanism that combines dynamical scattering and tidal interactions with a circumstellar disk. This combined model naturally produces the observed range of semi-major axis and eccentricity. We discuss the relative merits of the different migration mechanisms for producing the observed eccentric giant planets.

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.

The Sun's gradual brightening will seriously compromise the Earth's biosphere within ~ 1E9 years. If Earth's orbit migrates outward, however, the biosphere could remain intact over the entire main-sequence lifetime of the Sun. In this paper, we explore the feasibility of engineering such a migration over a long time period. The basic mechanism uses gravitational assists to (in effect) transfer orbital energy from Jupiter to the Earth, and thereby enlarges the orbital radius of Earth. This transfer is accomplished by a suitable intermediate body, either a Kuiper Belt object or a main belt asteroid. The object first encounters Earth during an inward pass on its initial highly elliptical orbit of large (~ 300 AU) semimajor axis. The encounter transfers energy from the object to the Earth in standard gravity-assist fashion by passing close to the leading limb of the planet. The resulting outbound trajectory of the object must cross the orbit of Jupiter; with proper timing, the outbound object encounters Jupiter and picks up the energy it lost to Earth. With small corrections to the trajectory, or additional planetary encounters (e.g., with Saturn), the object can repeat this process over many encounters. To maintain its present flux of solar energy, the Earth must experience roughly one encounter every 6000 years (for an object mass of 1E22 g). We develop the details of this scheme and discuss its ramifications.

This paper argues that star forming environments should be classified into finer divisions than the traditional isolated and clustered modes. Using the observed set of galactic open clusters and theoretical considerations regarding cluster formation, we estimate the fraction of star formation that takes place within clusters. We find that less than 10% of the stellar population originates from star forming regions destined to become open clusters, confirming earlier estimates. The smallest clusters included in the observational surveys (having at least N=100 members) roughly coincide with the smallest stellar systems that are expected to evolve as clusters in a dynamical sense. We show that stellar systems with too few members N < N_\star have dynamical relaxation times that are shorter than their formation times (1-2 Myr), where the critical number of stars N_\star \approx 100. Our results suggest that star formation can be characterized by (at least) three principal modes: I. isolated singles and binaries, II. groups (N<N_\star), and III. clusters (N>N_\star). Many -- if not most -- stars form through the intermediate mode in stellar groups with 10<N<100. Such groups evolve and disperse much more rapidly than open clusters; groups also have a low probability of containing massive stars and are unaffected by supernovae and intense ultraviolet radiation fields. Because of their short lifetimes and small stellar membership, groups have relatively little effect on the star formation process (on average) compared to larger open clusters.

Using the observed properties of our solar system, in particular the isotopic compositions of meteorites and the regularity of the planetary orbits, we constrain the star formation environment of the Sun within the scenario of (external) radioactive enrichment by a massive star. This calculation yields a probability distribution for the number of stars in the solar birth aggregate. The Sun is most likely to have formed within a stellar group containing $N$ = 2000 $\pm$ 1100 members. The {\it a priori} probability of a star forming in this type of environment is ${\cal P} \approx$ 0.0085, i.e., only about 1 out of 120 solar systems are expected to form under similar conditions. We discuss additional implications of this scenario, including possible effects from the radiation fields provided by the putative cluster environment and dynamical disruption of the Kuiper Belt. The constraints of this paper place tight restrictions on the properties of the solar birth aggregate for the scenario of external enrichment by a massive star; alternately, these tight constraints slightly favor a self-enrichment scenario for the short-lived radioactive species.

We construct a model for the formation of black holes within galactic bulges. The initial state is a slowly rotating isothermal sphere, characterized by effective transport speed $\aeff$ and rotation rate $\Omega$. The black hole mass is determined when the centrifugal radius of the collapse flow exceeds the capture radius of the central black hole. This model reproduces the observed correlation between black hole masses and galactic velocity dispersions, $\mbh \approx 10^8 M_\odot (\sigma/200 \kms)^4$, where $\sigma = \sqrt{2} \aeff$. This model also predicts the ratio $\mrat$ of black hole mass to host mass: $\mrat$ $\approx$ 0.004 $(\sigma/200 \kms)$.

We construct a family of semi-analytic models for young open clusters, including a gaseous component and varying assumptions about the distribution function for the stellar component. The parameters of these models are informed by observed open clusters and general theoretical considerations regarding cluster formation. We use this framework to estimate the fraction of the stellar component that remains gravitationally bound after the gaseous component disperses. The remaining stellar fraction is a smooth function of the star formation efficiency, and depends on the distribution function of the stars. We calculate the fraction of stars remaining for this class of open cluster models and provide fitting formulae for representative cases.

We investigate the equilibrium properties of self-gravitating magnetized clouds with polytropic equations of state with negative index n. In particular, we consider scale-free isopedic configurations that have constant dimensionless spherical mass-to-flux ratio lambda_r and that may constitute ``pivotal'' states for subsequent dynamical collapse to form groups or clusters of stars. For given Gamma = 1 + 1/n, equilibria with smaller values of lambda_r are more flattened, ranging from spherical configurations with lambda_r = infinity to completely flattened states for lambda_r=1. For a given amount of support provided by the magnetic field as measured by the dimensionless parameter H_0, equilibria with smaller values of Gamma are more flattened. However, logatropic (defined by Gamma = 0) disks do not exist. The only possible scale-free isopedic equilibria with logatropic equation of state are spherical uniformly magnetized clouds.

We summarize recent observational and theoretical progress aimed at understanding the origin of the stellar initial mass function (IMF) with specific focus on galactic star-forming regions. We synthesize data from various efforts to determine the IMF in very young, partially-embedded stellar clusters and find: i) no significant variations in the low-mass IMF have been observed between different star-forming regions; and ii) the mass distributions of young stars just emerging from molecular clouds are consistent with having been drawn from the IMF derived from field stars in the solar neighborhood. These results apply only to gross characterizations of the IMF (e.g. the ratio of high to low mass stars); present observations do not rule out more subtle regional differences. Further studies are required in order to assess whether or not there is evidence for a universal turnover near the hydrogen--burning limit. We also provide a general framework for discussing theories of the IMF, and summarize recent work on several physical mechanisms which could play a role in determining the form of the stellar initial mass function.

We explore possible effects of vacuum energy on the evolution of black holes. If the universe contains a cosmological constant, and if black holes can absorb energy from the vacuum, then black hole evaporation could be greatly suppressed. For the magnitude of the cosmological constant suggested by current observations, black holes larger than $\sim 4 \times 10^{24}$ g would accrete energy rather than evaporate. In this scenario, all stellar and supermassive black holes would grow with time until they reach a maximum mass scale of $\sim 6 \times 10^{55}$ g, comparable to the mass contained within the present day cosmological horizon.

We consider the long term fate and evolution of cold degenerate stars under the action of gravity alone. Although such stars cannot emit radiation through the Hawking mechanism, the wave function of the star will contain a small admixture of black hole states. These black hole states will emit radiation and hence the star can lose its mass energy in the long term. We discuss the allowed range of possible degenerate stellar evolution within this framework.

We present a series of 2-dimensional hydrodynamic simulations of massive disks around protostars. We simulate the same physical problem using both a `Piecewise Parabolic Method' (PPM) code and a `Smoothed Particle Hydrodynamic' (SPH) code, and analyze their differences. The disks studied here range in mass from $0.05 M_*$ to $1.0 M_*$ and in initial minimum Toomre $Q$ value from 1.1 to 3.0. For this problem, the strengths of the codes overlap only in a limited fashion, but similarities exist in their predictions, including spiral arm pattern speeds and morphological features. Our results represent limiting cases (i.e. systems evolved isothermally) rather than true physical systems. Disks become active from the inner regions outward. From the earliest times, their evolution is a strongly dynamic process rather than a smooth progression toward eventual nonlinear behavior. We calculate approximate growth rates for the spiral patterns; the one-armed ($m=1$) spiral arm is not the fastest growing pattern of most disks. In our SPH simulations, disks with initial minimum $Q=1.5$ or lower break up into proto-binary or proto-planetary clumps. However, these simulations cannot follow the physics important for the flow and must be terminated before the system has completely evolved. At their termination, PPM simulations with similar initial conditions show uneven mass distributions within spiral arms, suggesting that clumping behavior might result if they were carried further. Concern that the point-like nature of SPH exaggerates clumping, that our representation of the gravitational potential in PPM is too coarse, and that our physics assumptions are too simple, suggest caution in interpretation of the clumping in both the disk and torus simulations.

Motivated by recent experimental proposals to search for extragalactic cosmic rays (including anti-matter from distant galaxies), we study particle propagation through the intergalactic medium (IGM). We first use estimates of the magnetic field strength between galaxies to constrain the mean free path for diffusion of particles through the IGM. We then develop a simple analytic model to describe the diffusion of cosmic rays. Given the current age of galaxies, our results indicate that, in reasonable models, a completely negligible number of particles can enter our Galaxy from distances greater than $\sim 100$ Mpc for relatively low energies ($E$ $< 10^6$ GeV/n). We also find that particle destruction in galaxies along the diffusion path produces an exponential suppression of the possible flux of extragalactic cosmic rays. Finally, we use gamma ray constraints to argue that the distance to any hypothetical domains of anti-matter must be roughly comparable to the horizon scale.

A number of recently discovered extrasolar planet candidates have surprisingly small orbits, which may indicate that considerable orbital migration takes place in protoplanetary systems. A natural consequence of orbital migration is for a series of planets to be accreted, destroyed, and then thoroughly mixed into the convective envelope of the central star. We study the ramifications of planet accretion for the final main sequence metallicity of the star. If maximum disk lifetimes are on the order of 10 Myr, stars with masses near 1 solar mass are predicted to have virtually no metallicity enhancement. On the other hand, early F and late A type stars with masses of 1.5--2.0 solar masses can experience significant metallicity enhancements due to their considerably smaller convection zones during the first 10 Myr of pre-main-sequence evolution. We show that the metallicities of an aggregate of unevolved F stars are consistent with an average star accreting about 2 Jupiter-mass planets from a protoplanetary disk having a 10 Myr dispersal time.

This paper outlines astrophysical issues related to the long term fate of the universe. We consider the evolution of planets, stars, stellar populations, galaxies, and the universe itself over time scales which greatly exceed the current age of the universe. (shortened version of abstract).

We study the linear evolution of small perturbations in self-gravitating fluid systems with magnetic fields. We consider wave-like perturbations to nonuniform filamentary and sheet-like hydrostatic equilibria in the presence of a uniform parallel magnetic field. Motivated by observations of molecular clouds that suggest substantial nonthermal (turbulent) pressure, we adopt equations of state that are softer than isothermal. We numerically determine the dispersion relation and the form of the perturbations in the regime of instability. The form of the dispersion relation is the same for all equations of state considered, for all magnetic field strengths, and for both geometries examined. We demonstrate the existence of a fastest growing mode for the system and study how its characteristics depend on the amount of turbulence and the strength of the magnetic field. Generally, turbulence tends to increase the rate and the length scale of fragmentation. While tending to slow the fragmentation, the magnetic field has little effect on the fragmentation length scale until reaching some threshold, above which the length scale decreases significantly. Finally, we discuss the implications of these results for star formation in molecular clouds.

Motivated by recent measurements which suggest that roughly half the mass of the galactic halo may be in the form of white dwarfs, we study the implications of such a halo. We first use current limits on the infrared background light and the galactic metallicity to constrain the allowed initial mass function (IMF) of the stellar population that produced the white dwarfs. The IMF must be sharply peaked about a characteristic mass scale $M_C \approx 2.3 M_\odot$. Since only a fraction of the initial mass of a star is incorporated into the remnant white dwarf, we argue that the mass fraction of white dwarfs in the halo is likely to be 25\% or less, and that 50\% is an extreme upper limit. We use the IMF results to place corresponding constraints on the primordial initial conditions for star formation. The initial conditions must be much more homogeneous and skewed toward higher temperatures ($T_{\rm gas} \sim$ 200 K) than the conditions which lead to the present day IMF. Next we determine the luminosity function of white dwarfs. By comparing this result with the observed luminosity function, we find that the age of the halo population must be greater than $\sim 16$ Gyr. Finally, we calculate the radiative signature of a white dwarf halo. This infrared background is very faint, but is potentially detectable with future observations.

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.

In this paper, we study infall collapse solutions for star formation in the small radius limit where the particle orbits become nearly pressure-free. We generalize previous solutions to simultaneously include the effects of both radiation pressure and angular momentum. The effects of radiation pressure can be modeled using a modified potential; for representative cases of such potentials, we obtain analytical solutions for the density and velocity fields. In general, radiation pressure limits the maximum mass of a forming star by reversing the infall when the star becomes sufficiently large. Our results imply that this maximum mass scale is given by the condition that the turnaround radius $R_R$ (the radius at which the radiation pressure force exceeds the gravitational force) exceeds the centrifugal radius $R_C$ (the angular momentum barrier). The maximum mass scale for a star forming within a rotating collapse flow with radiation pressure depends on the initial conditions, but is generally much larger than for the case of spherical infall considered previously. In particular, stars with masses $M_\ast$ $\sim 100$ $M_\odot$ can form for a fairly wide range of initial conditions.

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 discuss the physics of vortices in the circumstellar disks associated with young stellar objects. We elucidate the basic physical properties of these localized storm systems. In particular, we consider point vortices, linear vortices, the effects of self-gravity, magnetic fields, and nonlinear aspects of the problem. We find that these vortices can exist in many different forms in the disks of young stellar objects and may play a role in the formation of binary companions and/or giant planets. Vortices may enhance giant planet formation via gravitational instability by allowing dust grains (heavy elements) to settle to the center on a short timescale; the gravitational instability itself is also enhanced because the vortices also create a larger local surface density in the disk. In addition, vortices can enhance energy dissipation in disks and thereby affect disk accretion. Finally, we consider the possibility that vortices of this type exist in molecular clouds and in the disk of the galaxy itself. On all of these size scales, vortices can produce long-lived structures which may correspond to observed structures in these systems.

We present an analysis of IRAS maps of five molecular clouds: Orion, Ophiuchus, Perseus, Taurus, and Lupus. For the classification and description of these astrophysical maps, we use a newly developed technique which considers all maps of a given type to be elements of a pseudometric space. For each physical characteristic of interest, this formal system assigns a distance function (a pseudometric) to the space of all maps; this procedure allows us to measure quantitatively the difference between any two maps and to order the space of all maps. We thus obtain a quantitative classification scheme for molecular clouds. In this present study we use the IRAS continuum maps at 100$\mu$m and 60$\mu$m to produce column density (or optical depth) maps for the five molecular cloud regions given above. For this sample of clouds, we compute the ``output'' functions which measure the distribution of density, the distribution of topological components, the self-gravity, and the filamentary nature of the clouds. The results of this work provide a quantitative description of the structure in these molecular cloud regions. We then order the clouds according to the overall environmental ``complexity'' of these star forming regions. Finally, we compare our results with the observed populations of young stellar objects in these clouds and discuss the possible environmental effects on the star formation process. Our results are consistent with the recently stated conjecture that more massive stars tend to form in more ``complex'' environments.

We extend a newly developed formal system for the description of astrophysical maps. In this formalism, we consider the difference between maps to be the distance between elements of a pseudometric space (the space of all such maps). This ansatz allows us to measure quantitatively the difference between any two maps and to order the space of all maps. For each physical characteristic of interest, this technique assigns an ``output'' function to each map; the difference between the maps is then determined from the difference between their corresponding output functions. In this present study, we show that the results of this procedure are invariant under a class of transformations of the maps and the domains of the maps. In addition, we study the propagation of errors (observational uncertainties) through this formalism. We show that the uncertainties in the output functions can be controlled provided that the signal to noise ratios in the original astrophysical maps are sufficiently high. The results of this paper thus increase the effectiveness of this formal system for the description, classification, and analysis of astrophysical maps.

In this paper, we present quantitative constraints on the scalar field potential for a general class of inflationary models. (1) We first consider the reconstruction of the inflationary potential for given primordial density fluctuation spectra. Our work differs from previous work on reconstruction in that we find a semi-analytic solution for the potential for the case of density fluctuations with power-law spectra. In addition, for the case of more general spectra, we show how constraints on the density fluctuation spectra imply corresponding constraints on the potential. We present a series of figures which show how the shape of the potential depends on the shape of the perturbation spectrum and on the relative contribution of tensor modes. (2) We show that the average ratio $\rave$ of the amplitude of tensor perturbations (gravity wave perturbations) to scalar density perturbations is bounded from above: $\rave \le$ 1.6. We also show that the ratio $\rave$ is proportional to the change $\Delta \phi$ in the field: $\rave \approx 0.42 \Delta \phi/\mp$. Thus, if tensor perturbations are important for the formation of structure, then the width $\Delta \phi$ must be comparable to the Planck mass. (3) We constrain the change $\Delta V$ of the potential and the change $\Delta \phi$ of the inflation field during the portion of inflation when cosmological structure is produced. We find both upper and lower bounds for $\Delta \phi$ and for $\Delta V$. In addition, these constraints are then used to derive a bound on the scale $\Lambda$, which is the scale of the height of the potential

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.

A pseudo-Nambu-Goldstone boson, with a potential of the form $V(\phi) = \Lambda^4[1 \pm \cos(\phi/f)], naturally gives rise to inflation if $f \sim M_{Pl}$ and $\Lambda \sim M_{GUT}$. We show how this can arise in technicolor-like and superstring models, and work out an explicit string example in the context of multiple gaugino condensation models. We study the cosmology of this model in detail, and find that sufficient reheating to ensure that baryogenesis can take place requires $f > 0.3 M_{Pl}$. The primordial density fluctuation spectrum generated is a non-scale-invariant power law, $P(k) \propto k^{n_s}$, with $n_s \simeq 1 - (M^2_{Pl}/8\pi f^2)$, leading to more power on large length scales than the $n_s = 1$ Harrison-Zeldovich spectrum. The standard CDM model with $0 \la n_s \la 0.6-0.7$ could in principle explain the large-scale clustering observed in the APM and IRAS galaxy surveys as well as large-scale flows, but the COBE microwave anisotropy implies such low amplitudes (or high bias factors, $b>2$) for these CDM models that galaxy formation occurs too late to be viable; combining COBE with sufficiently early galaxy formation or the large-scale flows leads to $n_s >0.6$, or $f > 0.3 M_{Pl}$ as well. For extended and power law inflation models, this constraint is even tighter, $n_s > 0.7$; combined with other bounds on large bubbles in extended inflation, this leaves little room for most extended models.