BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Fred C. Adams

Fred C. Adams appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-phastro-ph.COastro-ph.GAmath-phmath.MPgr-qcastro-ph.HEastro-ph.IMMachine Learningnlin.AO

Catalog footprint

What is connected

80works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

80 published item(s)

preprint2022arXiv

A Collision Mechanism for the Removal of Earth's Trojan Asteroids

Due to their strong resonances with their host planet, Trojan asteroids can remain in stable orbits for billions of years. As a result, they are powerful probes for constraining the dynamical and chemical history of the solar system. Although we have detected thousands of Jupiter Trojans and dozens of Neptune Trojans, there are currently no known long-term stable Earth Trojans. Dynamical simulations show that the parameter space for stable Earth Trojans in substantial, so their apparent absence poses a mystery. This work uses a large ensemble of N-body simulations to explore how the Trojan population dynamically responds if Earth suffers large collisions, such as those thought to have occurred to form the Moon and/or to have given Earth its Late Veneer. We show that such collisions can be highly disruptive to the primordial Trojan population, and could have eliminated it altogether. More specifically, if Earth acquired the final 1\% of its mass through ${\cal O}(10)$ collisions, then only $\sim1\%$ of the previously bound Trojan population would remain.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gemini-LIGHTS: Herbig Ae/Be and massive T-Tauri protoplanetary disks imaged with Gemini Planet Imager

We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of transitional, pre-transitional, and full disks. Our approach explicitly did not favor well-known, "famous" disks or those observed by ALMA, resulting in a less-biased sample suitable to probe the major stages of disk evolution during planet formation. Our optimized data reduction allowed polarized flux as low as 0.002% of the stellar light to be detected, and we report polarized scattered light around 80% of our targets. We detected point-like companions for 47% of the targets, including 3 brown dwarfs (2 confirmed, 1 new), and a new super-Jupiter mass candidate around V1295 Aql. We searched for correlations between the polarized flux and system parameters, finding a few clear trends: presence of a companion drastically reduces the polarized flux levels, far-IR excess correlates with polarized flux for non-binary systems, and systems hosting disks with ring structures have stellar masses $<$ 3 Msun. Our sample also included four hot, dusty "FS CMa" systems and we detected large-scale ($>100$ au) scattered light around each, signs of extreme youth for these enigmatic systems. Science-ready images are publicly available through multiple distribution channels using a new FITS file standard jointly developed with members of the VLT/SPHERE team.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Inferring Late Stage Enrichment of Exoplanet Atmospheres from Observed Interstellar Comets

The discovery of the first two interstellar objects implies that, on average, every star contributes a substantial amount of material to the galactic population by ejecting such bodies from the host system. Since scattering is a chaotic process, a comparable amount of material should be injected into the inner regions of each system that ejects comets. For comets that are transported inwards and interact with planets, this Letter estimates the fraction of material that is accreted or outward-scattered as a function of planetary masses and orbital parameters. These calculations indicate that planets with escape velocities smaller than their current day orbital velocities will efficiently accrete comets. We estimate the accretion efficiency for members of the current census of extrasolar planets, and find that planetary populations including but not limited to hot and warm Jupiters, sub-Neptunes and super-Earths can efficiently capture incoming comets. This cometary enrichment may have important ramifications for post-formation atmospheric composition and chemistry. As a result, future detections and compositional measurements of interstellar comets will provide direct measurements of material that potentially enriched a sub-population of the extrasolar planets. Finally, we estimate the efficiency of this enrichment mechanism for extrasolar planets that will be observed with the $\textit{James Webb Space Telescope}$ (JWST). With JWST currently operational and these observations imminently forthcoming, it is of critical importance to investigate how enrichment from interstellar comet analogues may affect the interpretations of exoplanet atmospheric compositions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Testing Self-Organized Criticality Across the Main Sequence using Stellar Flares from TESS

Self-organized criticality describes a class of dynamical systems that maintain themselves in an attractor state with no intrinsic length or time scale. Fundamentally, this theoretical construct requires a mechanism for instability that may trigger additional instabilities locally via dissipative processes. This concept has been invoked to explain nonlinear dynamical phenomena such as featureless energy spectra that have been observed empirically for earthquakes, avalanches, and solar flares. If this interpretation proves correct, it implies that the solar coronal magnetic field maintains itself in a critical state via a delicate balance between the dynamo-driven injection of magnetic energy and the release of that energy via flaring events. All-sky high-cadence surveys like the Transiting Exoplanet Survey Satellite (TESS) provide the necessary data to compare the energy distribution of flaring events in stars of different spectral types to that observed in the Sun. We identified $\sim 10^6$ flaring events on $\sim 10^5$ stars observed by TESS at 2-minute cadence. By fitting the flare frequency distribution for different mass bins, we find that all main sequence stars exhibit distributions of flaring events similar to that observed in the Sun, independent of their mass or age. This may suggest that stars universally maintain a critical state in their coronal topologies via magnetic reconnection events. If this interpretation proves correct, we may be able to infer properties of magnetic fields, interior structure, and dynamo mechanisms for stars that are otherwise unresolved point sources.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Theoretical and Observational Evidence for Coriolis Effects in Coronal Magnetic Fields Via Direct Current Driven Flaring Events

All stars produce explosive surface events such as flares and coronal mass ejections. These events are driven by the release of energy stored in coronal magnetic fields, generated by the stellar dynamo. However, it remains unclear if the energy deposition in the magnetic fields is driven by direct or alternating currents. Recently, we presented observational measurements of the flare intensity distributions for a sample of $\sim10^5$ stars across the main sequence observed by $\textit{TESS}$, all of which exhibited power-law distributions similar to those observed in the Sun, albeit with varying slopes. Here we investigate the mechanisms required to produce such a distribution of flaring events via direct current energy deposition, in which coronal magnetic fields braid, reconnect, and produce flares. We adopt a topological model for this process which produces a power-law distribution of energetic flaring events. We expand this model to include the Coriolis effect, which we demonstrate produces a shallower distribution of flare energies in stars that rotate more rapidly (corresponding to a weaker decline in occurrence rates toward increasing flare energies). We present tentative evidence for the predicted rotation-power-law index correlation in the observations. We advocate for future observations of stellar flares that would improve our measurements of the power-law exponents, and yield key insights into the underlying dynamo mechanisms that underpin the self-similar flare intensity distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Coupled Analysis of Atmospheric Mass Loss and Tidal Evolution in XUV Irradiated Exoplanets: the TRAPPIST-1 Case Study

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Giant Planet Candidate Transiting a White Dwarf

Astronomers have discovered thousands of planets outside the solar system, most of which orbit stars that will eventually evolve into red giants and then into white dwarfs. During the red giant phase, any close-orbiting planets will be engulfed by the star, but more distant planets can survive this phase and remain in orbit around the white dwarf. Some white dwarfs show evidence for rocky material floating in their atmospheres, in warm debris disks, or orbiting very closely, which has been interpreted as the debris of rocky planets that were scattered inward and tidally disrupted. Recently, the discovery of a gaseous debris disk with a composition similar to ice giant planets demonstrated that massive planets might also find their way into tight orbits around white dwarfs, but it is unclear whether the planets can survive the journey. So far, the detection of intact planets in close orbits around white dwarfs has remained elusive. Here, we report the discovery of a giant planet candidate transiting the white dwarf WD 1856+534 (TIC 267574918) every 1.4 days. The planet candidate is roughly the same size as Jupiter and is no more than 14 times as massive (with 95% confidence). Other cases of white dwarfs with close brown dwarf or stellar companions are explained as the consequence of common-envelope evolution, wherein the original orbit is enveloped during the red-giant phase and shrinks due to friction. In this case, though, the low mass and relatively long orbital period of the planet candidate make common-envelope evolution less likely. Instead, the WD 1856+534 system seems to demonstrate that giant planets can be scattered into tight orbits without being tidally disrupted, and motivates searches for smaller transiting planets around white dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A triple star system with a misaligned and warped circumstellar disk shaped by disk tearing

Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple star system GWOrionis, finding evidence for disk tearing. Our images show an eccentric ring that is misaligned with the orbital planes and the outer disk. The ring casts shadows on a strongly warped intermediate region of the disk. If planets can form within the warped disk, disk tearing could provide a mechanism for forming wide-separation planets on oblique orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

ARMADA I: Triple Companions Detected in B-Type Binaries alpha Del and nu Gem

Ground-based optical long-baseline interferometry has the power to measure the orbits of close binary systems at ~10 micro-arcsecond precision. This precision makes it possible to detect "wobbles" in the binary motion due to the gravitational pull from additional short period companions. We started the ARrangement for Micro-Arcsecond Differential Astrometry (ARMADA) survey with the MIRC-X instrument at the CHARA array for the purpose of detecting giant planets and stellar companions orbiting individual stars in binary systems. We describe our observations for the survey, and introduce the wavelength calibration scheme that delivers precision at the tens of micro-arcseconds level for <0.2 arcsecond binaries. We test our instrument performance on a known triple system kappa Peg, and show that our survey is delivering a factor of 10 better precision than previous similar surveys. We present astrometric detections of tertiary components to two B-type binaries: a 30-day companion to alpha Del, and a 50-day companion to nu Gem. We also collected radial velocity data for alpha Del with the Tennessee State University Automated Spectroscopic Telescope at Fairborn Observatory. We are able to measure the orbits and masses of all three components in these systems. We find that the previously published RV orbit for the inner pair of nu Gem is not consistent with our visual orbit. The precision achieved for these orbits suggests that our ARMADA survey will be successful at discovering new compact triple systems to A/B-type binary systems, leading to better statistics of hierarchical system architectures and formation history.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Detection of Diatomic Carbon in 2I/Borisov

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamics of Planetary Systems Within Star Clusters: Aspects of the Solar System's Early Evolution

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°$, and are thus insufficient to explain the $ψ\approx6°$ 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, $η\,τ\lesssim2\times10^4\,$Myr/pc$^3$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Energy Optimization in Binary Star Systems: Explanation for Equal Mass Members in Close 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$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Energy Optimization in Extrasolar Planetary Systems: The Transition from Peas-in-a-Pod to Runaway Growth

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Search for L5 Earth Trojans with DECam

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Testing the isotropy of the Dark Energy Survey's extreme trans-Neptunian objects

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 $Ω,$ and argument of perihelion $ω.$ We then compare the detected distributions in each of $Ω, ω,$ and $\varpi\equivΩ+ω$ 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 $Ω$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The First Habitable Zone Earth-sized Planet from TESS. I: Validation of the TOI-700 System

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Solar wind prevents re-accretion of debris after Mercury's giant impact

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Trans-Neptunian objects found in the first four years of the Dark Energy Survey

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Irregular dust features around intermediate-mass young stars with GPI: signs of youth or misaligned disks?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Evidence for Color Dichotomy in the Primordial Neptunian Trojan Population

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Constraints on Alternate Universes: Stars and habitable planets with different fundamental constants

This paper develops constraints on the values of the fundamental constants that allow universes to be habitable. We focus on the fine structure constant $α$ and the gravitational structure constant $α_G$, and find the region in the $α$-$α_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 $α$-$α_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 $α$ and $α_G$ vary by several orders of magnitude. These constraints also provide upper bounds on the structure constants ($α,α_G$) and their ratio. We find the limit $α_G/α<10^{-34}$, which shows that habitable universes must have a large hierarchy between the strengths of the gravitational force and the electromagnetic force.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Five Planets Transiting a Ninth Magnitude Star

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Interaction Cross Sections and Survival Rates for Proposed Solar System Member Planet Nine

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 $σ_{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 $Δ{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\%)$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

MagAO Imaging of Long-period Objects (MILO). I. A Benchmark M Dwarf Companion Exciting a Massive Planet around the Sun-like Star HD 7449

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The coronal evolution of pre-main-sequence stars

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The influence of radiative core growth on coronal X-ray emission from pre-main sequence stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Stability of Tidal Equilibrium for Hierarchical Star-Planet-Moon Systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cross Sections for Planetary Systems Interacting with Passing Stars and Binaries

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hot Jupiter Breezes: Time-dependent Outflows from Extrasolar Planets

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Oscillations of Relative Inclination Angles in Compact Extrasolar Planetary Systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Planets in Other Universes: Habitability constraints on density fluctuations and galactic structure

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resonant Removal of Exomoons During Planetary Migration

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Five Planets in the Kepler-296 Binary System All Orbit the Primary: A Statistical and Analytical Analysis

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Torque on an exoplanet from an anisotropic evaporative wind

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Alignment of Protostars and Circumstellar Disks During the Embedded Phase

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An Earth-sized Planet in the Habitable Zone of a Cool Star

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constraining the X-ray and Cosmic Ray Ionization Chemistry of the TW Hya Protoplanetary Disk: Evidence for a Sub-interstellar Cosmic Ray Rate

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, $ζ_{\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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Distributions of Short-Lived Radioactive Nuclei Produced by Young Embedded Stellar Clusters

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 $ζ_{SLR}\sim1-5\times10^{-19}$ sec$^{-1}$. This ionization rate is smaller than that due to cosmic rays, $ζ_{CR}\sim10^{-17}$ sec$^{-1}$, but will be important in regions where cosmic rays are attenuated (e.g., disk mid-planes).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Effects of Turbulence on Cosmic Ray Propagation in Protostars and Young Star/Disk Systems

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exclusion of Cosmic Rays in Protoplanetary Disks. II. Chemical Gradients and Observational Signatures

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 ($ζ\gtrsim10^{-17}$ s$^{-1}$ per H$_2$) and below, but it becomes difficult to distinguish between low ionization models when $ζ\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 ($ζ_{\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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetically controlled mass loss from extrasolar planets in close orbits

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetically Controlled Outflows from Planets

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Stability of Extrasolar Planetary Systems and other Closely Orbiting Pairs

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Wind-driven Exclusion of Cosmic Rays in the Protoplanetary Disk Environment

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evolution of Planetary Orbits with Stellar Mass Loss and Tidal Dissipation

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 γand an index β. We use an analogous prescription to model tidal effects, described here by a dimensionless dissipation rate Γand two indices (q,p). The initial conditions are determined by the starting value of angular momentum parameter η(equivalently, the initial eccentricity) and the phase θof the orbit. Within the context of this model, we derive an analytic formula for the critical dissipation rate Γ, 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 Γ=Γ(γ,β,q,p,η,θ) 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evolution of Planetary Systems with Time Dependent Stellar Mass Loss

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exclusion of Cosmic Rays in Protoplanetary Disks: Stellar and Magnetic Effects

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hill's Equation with Small Fluctuations: Cycle to Cycle Variations and Stochastic Processes

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 λ_k, and the period τ_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 ξ. This paper considers both random and stochastic Hill's equations with small parameter variations, so that p_k=q_k-<q_k>, \ell_k=λ_k-<λ_k>, and ξare all O(ε), where ε<<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 ξ. 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 λ-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 γ= O(ε^2). Using the relationship between the (\ell_k,p_k) and the ξ, this result for γcan be used to find growth rates for stochastic Hill's equations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Magnetic and Gravitational Disk-Star Interactions: An Interdependence of PMS Stellar Rotation Rates and Spin-Orbit Misalignments

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Radionuclide Ionization in Protoplanetary Disks: Calculations of Decay Product Radiative Transfer

We present simple analytic solutions for the ionization rate $ζ_{\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, $ζ_{\rm{CR}}$, equal to or substantially below that of SLRs ($ζ_{\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 $ζ_{\rm{SLR}}\sim(1-10)\times10^{-19}$ s$^{-1}$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Long-Term Dynamical Evolution of Planetary Systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Viscous Evolution and Photoevaporation of Circumstellar Disks due to External FUV Radiation Fields

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 α-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 α= 10^{-3} are effectively dispersed within 1-3 Myr; for higher viscosities (α= 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effects of Collisions with Rocky Planets on the Properties of Hot Jupiters

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mean Motion Resonances in Exoplanet Systems: An Investigation Into Nodding Behavior

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Accretion of Rocky Planets by Hot Jupiters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Magnetic Interactions in Pre-Main-Sequence Binaries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Magnetically Controlled Accretion Flows onto Young Stellar Objects

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Magnetically Controlled Outflows from Hot Jupiters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the Anomalous Radii of the Transiting Extrasolar 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}^α$, with $α=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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An Investigation of the Loss of Planet-Forming Potential in Intermediate Sized Young Embedded Star Clusters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effects of Turbulence, Eccentricity Damping, and Migration Rate on the Capture of Planets into Mean Motion Resonance

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Five Long-period Extrasolar Planets in Eccentric orbits from the Magellan Planet Search Program

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hill's Equation with Random Forcing Parameters: Determination of Growth Rates through Random Matrices

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stability of Magnetized Disks and Implications for Planet Formation

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Birth Environment of the Solar System

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dynamical Evolution of Young Embedded Clusters: A Parameter Space Survey

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

General Analysis of Type I Planetary Migration with Stochastic Perturbations

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hill's Equation with Random Forcing Parameters: The Limit of Delta Function Barriers

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{ö}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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Kinematic Signatures of Subvirial Initial Conditions in Young Clusters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Migration of Extrasolar Planets: Effects from X-Wind Accretion Disks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Low Mass Companions for Five Solar-Type Stars from the Magellan Planet Search Program

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Mean Motion Resonances in Extrasolar Planetary Systems with Turbulence, Interactions, and Damping

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Terrestrial Planet Formation Around Individual Stars Within Binary Star Systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The Asymptotic Form of Cosmic Structure: Small Scale Power and Accretion History

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, γ, that controls the rate of mass accretion, dln(M) / dln(a) ~ a^(-γ), 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.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Core-Accretion Model Predicts Few Jovian-Mass Planets Orbiting Red Dwarfs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The Ultimate Halo Mass in a LCDM Universe

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.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Future Evolution of Structure in an Accelerating Universe

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.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Constraints on the Intergalactic Transport of Cosmic Rays

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.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Linear Graviational Instability of Filamentary and Sheet-like Molecular Clouds with Magnetic Fields

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.

0 citations0 reviewsNo code linkedOpen access