Source author record

Peter Bodenheimer

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

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

27 published item(s)

preprint2022arXiv

Mixing of Condensable Constituents with H-He During the Formation & Evolution of Jupiter

Simulations of Jupiter's formation are presented that incorporate mixing of H-He with denser material entering the planet as solids. Heavy compounds and gas mix substantially when the planet becomes roughly as massive as Earth, because incoming planetesimals can fully vaporize. Supersaturation of vaporized silicates causes the excess to sink as droplets, but water remains at higher altitudes. Because the mean molecular weight decreases rapidly outward, some of the compositional inhomogeneities produced during formation can survive for billions of years. After 4.57 Gyr, our Jupiter model retains compositional gradients; proceeding outwards one finds: i) an inner heavy-element core, the outer part derived from hot supersaturated rain-out; ii) a composition-gradient region, containing most of the heavy elements, where H-He abundance increases outward, reaching about 0.9 mass fraction at 0.3 of the radius, with silicates enhanced relative to water in the lower parts and depleted in the upper parts; iii) a uniform composition region (neglecting He immiscibility) that is enriched over protosolar and contains most of the planet's mass; and iv) an outer region where cloud formation (condensation) of heavy constituents occurs. This radial compositional profile has heavy elements more broadly distributed than predicted by classical formation models, but less diluted than suggested by Juno-constrained gravity models. The compositional gradients in the region containing the bulk of the heavy elements prevent convection, in both our models and those fitting current gravity, resulting in a hot interior where much of the accretion energy remains trapped.

preprint2020arXiv

Detailed Calculations of the Efficiency of Planetesimal Accretion in the Core-Accretion Model

We present results of a detailed study of the rate of the accretion of planetesimals by a growing proto-Jupiter in the core-accretion model. Using a newly developed code, we accurately combine a detailed three-body trajectory calculation with gas drag experienced during the passage of planetesimals in the protoplanet's envelope. We find that the motion of planetesimals is excited to the extent that encounters with the proto-planetary envelope become so fast that ram pressure breaks up the planetesimals in most encounters. As a result, the accretion rate is largely independent of the planetesimal size and composition. For the case we explored of a planet forming at 5.2 AU from the Sun in a disk with a solid surface density of 6 g/cm^2 (Lozovsky et al. 2017) the accretion rate we compute differs in several respects from that assumed by those authors. We find that only 4-5 M_Earth is accreted in the first 1.5x10^6 years before the onset of rapid gas accretion. Most of the mass, some 10 M_Earth, is accreted simultaneously with this rapid gas accretion. In addition, we find that the mass accretion rate remains small, but non-zero for at least a million years after this point, and an additional 0.3-0.4 M_Earth is accreted during that time. This late accretion, together with a rapid infall of gas could lead to the accreted material being mixed throughout the outer regions, and may account for the enhancement of high-Z material in Jupiter's envelope.

preprint2016arXiv

In Situ and Ex Situ Formation Models of Kepler 11 Planets

We present formation simulations of the six Kepler 11 planets. Models assume either in situ or ex situ assembly, the latter with migration, and are evolved to the estimated age of the system, 8 Gyr. Models combine detailed calculations of both the gaseous envelope and the condensed core structures, including accretion of gas and solids, of the disk's viscous and thermal evolution, including photo-evaporation and disk-planet interactions, and of the planets' evaporative mass loss after disk dispersal. Planet-planet interactions are neglected. Both sets of simulations successfully reproduce measured radii, masses, and orbital distances of the planets, except for the radius of Kepler 11b, which loses its entire gaseous envelope shortly after formation. Gaseous (H+He) envelopes account for < 18% of the planet masses, and between 35 and 60% of the planet radii. In situ models predict a very massive inner disk, whose solids' surface density (sigma_Z) varies from over 1e4 to 1e3 g/cm2 at stellocentric distances 0.1 < r < 0.5 AU. Initial gas densities would be in excess of 1e5 g/cm2 if solids formed locally. Given the high disk temperatures (> 1000 K), planetary interiors can only be composed of metals and highly refractory materials. Sequestration of hydrogen by the core and subsequent outgassing is required to account for the observed radius of Kepler 11b. Ex situ models predict a relatively low-mass disk, whose initial sigma_Z varies from 10 to 5 g/cm2 at 0.5 < r < 7 AU and whose initial gas density ranges from 1e3 to 100 g/cm2. All planetary interiors are expected to be rich in H2O, as core assembly mostly occurs exterior to the ice condensation front. Kepler 11b is expected to have a steam atmosphere, and H2O is likely mixed with H+He in the envelopes of the other planets. Results indicate that Kepler 11g may not be more massive than Kepler 11e.

preprint2014arXiv

Accretion and Evolution of ~2.5 Earth-mass Planets with Voluminous H/He Envelopes

Formation of planets in the Neptune size range with low-mass, but voluminous, H_2/He gaseous envelopes is modeled by detailed numerical simulations according to the core-nucleated accretion scenario. Formation locations ranging from 0.5 to 4 AU from a star of 1 solar mass are considered. The final planets have heavy-element cores of 2.2--2.5 Earth masses and envelopes in the range 0.037--0.16 Earth masses. After the formation process, which lasts 2 Myr or less, the planets evolve at constant mass up to an age of several Gyr. For assumed equilibrium temperatures of 250, 500, and 1000 K, their calculated final radii are compared with those observed by the Kepler spacecraft. For the particular case of Kepler-11 f, we address the question whether it could have formed in situ or whether migration from a formation location farther out in the disk is required.

preprint2014arXiv

Growth of Jupiter: Enhancement of Core Accretion by a Voluminous Low-Mass Envelope

We present calculations of the early stages of the formation of Jupiter via core nucleated accretion and gas capture. The core begins as a seed body of about 350 kilometers in radius and orbits in a swarm of planetesimals whose initial radii range from 15 meters to 50 kilometers. The evolution of the swarm accounts for growth and fragmentation, viscous and gravitational stirring, and for drag-assisted migration and velocity damping. During this evolution, less than 9% of the mass is in planetesimals smaller than 1 kilometer in radius; < ~25% is in planetesimals with radii between 1 and 10 kilometers; and < ~7% is in bodies with radii larger than 100 kilometers. Gas capture by the core substantially enhances the size-dependent cross-section of the planet for accretion of planetesimals. The calculation of dust opacity in the planet's envelope accounts for coagulation and sedimentation of dust particles released as planetesimals are ablated. The calculation is carried out at an orbital semi-major axis of 5.2 AU and the initial solids' surface density is 10 g/cm^2 at that distance. The results give a core mass of nearly 7.3 Earth masses (Mearth) and an envelope mass of approximately 0.15 Mearth after about 4e5 years, at which point the envelope growth rate surpasses that of the core. The same calculation without the envelope yields a core of only about 4.4 Mearth.

preprint2014arXiv

The Formation of Uranus & Neptune: Challenges and Implications For Intermediate-Mass Exoplanets

In this paper we investigate the formation of Uranus and Neptune, according to the core-nucleated accretion model, considering formation locations ranging from 12 to 30 AU from the Sun, and with various disk solid-surface densities and core accretion rates. It is shown that in order to form Uranus-like and Neptune-like planets in terms of final mass and solid-to-gas ratio, very specific conditions are required. We also show that when recently proposed high solid accretion rates are assumed, along with solid surface densities about 10 times those in the minimum-mass solar nebula, the challenge in forming Uranus and Neptune at large radial distances is no longer the formation timescale, but is rather finding agreement with the final mass and composition of these planets. In fact, these conditions are more likely to lead to gas-giant planets. Scattering of planetesimals by the forming planetary core is found to be an important effect at the larger distances. Our study emphasizes how (even slightly) different conditions in the protoplanetary disk and the birth environment of the planetary embryos can lead to the formation of very different planets in terms of final masses and compositions (solid-to-gas ratios), which naturally explains the large diversity of intermediate-mass exoplanets.

preprint2013arXiv

Deuterium Burning in Massive Giant Planets and Low-Mass Brown Dwarfs formed by Core-Nucleated Accretion

Formation of bodies near the deuterium-burning limit is considered by detailed numerical simulations according to the core-nucleated giant planet accretion scenario. The objects, with heavy-element cores in the range 5-30 Mearth, are assumed to accrete gas up to final masses of 10-15 Jupiter masses (Mjup). After the formation process, which lasts 1-5 Myr and which ends with a 'cold-start', low-entropy configuration, the bodies evolve at constant mass up to an age of several Gyr. Deuterium burning via proton capture is included in the calculation, and we determined the mass, M50, above which more than 50% of the initial deuterium is burned. This often-quoted borderline between giant planets and brown dwarfs is found to depend only slightly on parameters, such as core mass, stellar mass, formation location, solid surface density in the protoplanetary disk, disk viscosity, and dust opacity. The values for M50 fall in the range 11.6-13.6 Mjup, in agreement with previous determinations that do not take the formation process into account. For a given opacity law during the formation process, objects with higher core masses form more quickly. The result is higher entropy in the envelope at the completion of accretion, yielding lower values of M50. For masses above M50, during the deuterium-burning phase, objects expand and increase in luminosity by 1 to 3 orders of magnitude. Evolutionary tracks in the luminosity-versus-time diagram are compared with the observed position of the companion to Beta Pictoris.

preprint2013arXiv

Numerical Evidence for Dark Star Formation: A Comment on "Weakly Interacting Massive Particle Dark Matter and First Stars: Suppression of Fragmentation in Primordial Star Formation" by Smith et al. 2012, ApJ 761, 154

(abridged) This comment is intended to show that simulations by Smith et al. (S12) support the Dark Star (DS) scenario and even remove some potential obstacles. Our previous work illustrated that the initial hydrogen densities of the first equilibrium DSs are high, ~10^{17}/cm^3 for the case of 100 GeV WIMPs, with a stellar radius of ~2-3 AU. Subsequent authors have somehow missed the fact that equilibrium DSs have the high densities they do. S12 have numerically simulated the effect of dark matter annihilation on the contraction of a protostellar gas cloud en route to forming the first stars. They show results at a density ~5 10^{14}/cm^3, slightly higher than the value at which annihilation heating prevails over cooling. However, they are apparently unable to reach the ~10^{17}/cm^3 density of our hydrostatic DS solutions. We are in complete agreement with their physical result that the gas keeps collapsing to densities > 5 10^{14}/cm^3, as it must before equilibrium DSs can form. However we are in disagreement with some of the words in their paper which imply that DSs never come to exist. It seems to us that S12 supports the DS scenario. They use the sink particle approach to treat the gas that collapses to scales smaller than their resolution limit. We argue that their sink is effectively a DS, or contains one. An accretion disk forms as more mass falls onto the sink, and the DS grows. S12 not only confirm our predictions about DS in the range where the simulations apply, but also solve a potential obstruction to DS formation by showing that dark matter annihilation prevents the fragmentation of the collapsing gas. Whereas fragmentation might perturb the dark matter away from the DS and remove its power source, instead S12 show that further sinks, if any, form only far enough away as to leave the DS undisturbed in the comfort of its dark matter surroundings.

preprint2013arXiv

Three-Dimensional Radiation-Hydrodynamics Calculations of the Envelopes of Young Planets Embedded in Protoplanetary Disks

We perform global three-dimensional (3D) radiation-hydrodynamics calculations of the envelopes surrounding young planetary cores of 5, 10, and 15 Earth masses, located in a protoplanetary disk at 5 and 10 AU from a solar-mass star. We apply a nested-grid technique to resolve the thermodynamics of the disk at the orbital-radius length scale and that of the envelope at the core-radius length scale. The gas is modeled as a solar mixture of molecular and atomic hydrogen, helium, and their ions. The equation of state accounts for both gas and radiation, and gas energy includes contributions from rotational and vibrational states of molecular hydrogen and from ionization of atomic species. Dust opacities are computed from first principles, applying the full Mie theory. One-dimensional (1D) calculations of planet formation are used to supplement the 3D calculations by providing energy deposition rates in the envelope due to solids accretion. We compare 1D and 3D envelopes and find that masses and gas accretion rates agree within factors of 2, and so do envelope temperatures. The trajectories of passive tracers are used to define the size of 3D envelopes, resulting in radii much smaller than the Hill radius and smaller than the Bondi radius. The moments of inertia and angular momentum of the envelopes are determined and the rotation rates are derived from the rigid-body approximation, resulting in slow bulk rotation. We find that the polar flattening is < ~0.05. The dynamics of the accretion flow is examined by tracking the motion of tracers that move into the envelope. The anisotropy of this flow is characterized in terms of both its origin and impact site at the envelope surface. Gas merges with the envelope preferentially at mid- to high latitudes.

preprint2011arXiv

Formation and Structure of Low Density Exo-Neptunes

Kepler has found hundreds of Neptune-size (2-6 R_Earth) planet candidates within 0.5 AU of their stars. The nature of the vast majority of these planets is not known because their masses have not been measured. Using theoretical models of planet formation, evolution and structure, we explore the range of minimum plausible masses for low-density exo-Neptunes. We focus on highly irradiated planets with T_eq>=500K. We consider two separate formation pathways for low-mass planets with voluminous atmospheres of light gases: core nucleated accretion and outgassing of hydrogen from dissociated ices. We show that Neptune-size planets at T_eq=500K with masses as small as a few times that of Earth can plausibly be formed core nucleated accretion coupled with subsequent inward migration. We also derive a limiting low-density mass-radius relation for rocky planets with outgassed hydrogen envelopes but no surface water. Rocky planets with outgassed hydrogen envelopes typically have computed radii well below 3 R_Earth. For both planets with H/He envelopes from core nucleated accretion and planets with outgassed hydrogen envelopes, we employ planet interior models to map the range of planet mass--envelope mass--equilibrium temperature parameter space that is consistent with Neptune-size planet radii. Atmospheric mass loss mediates which corners of this parameter space are populated by actual planets and ultimately governs the minimum plausible mass at a specified transit radius. We find that Kepler's 2-6 R_Earth planet candidates at T_eq=500--1000K could potentially have masses less than ~4 M_Earth. Although our quantitative results depend on several assumptions, our qualitative finding that warm Neptune-size planets can have masses substantially smaller than those given by interpolating the masses and radii of planets within our Solar System is robust.

preprint2010arXiv

Composition of Massive Giant Planets

The two current models for giant planet formation are core accretion and disk instability. We discuss the core masses and overall planetary enrichment in heavy elements predicted by the two formation models, and show that both models could lead to a large range of final compositions. For example, both can form giant planets with nearly stellar compositions. However, low-mass giant planets, enriched in heavy elements compared to their host stars, are more easily explained by the core accretion model. The final structure of the planets, i.e., the distribution of heavy elements, is not firmly constrained in either formation model.

preprint2010arXiv

Formation of Jupiter using opacities based on detailed grain physics

Numerical simulations, based on the core-nucleated accretion model, are presented for the formation of Jupiter at 5.2 AU in 3 primordial disks with three different assumed values of the surface density of solid particles. The grain opacities in the envelope of the protoplanet are computed using a detailed model that includes settling and coagulation of grains and that incorporates a recalculation of the grain size distribution at each point in time and space. We generally find lower opacities than the 2% of interstellar values used in previous calculations [Hubickyj, O., Bodenheimer, P., Lissauer, J. J., 2005. Icarus 179, 415--431; Lissauer, J. J., Hubickyj, O., D'Angelo, G., Bodenheimer, P., 2009. Icarus 199, 338-350]. These lower opacities result in more rapid heat loss from and more rapid contraction of the protoplanetary envelope. For a given surface density of solids, the new calculations result in a substantial speedup in formation time as compared with those previous calculations. Formation times are calculated to be 1.0, 1.9, and 4.0 Myr, and solid core masses are found to be 16.8, 8.9, and 4.7 Earth masses, for solid surface densities, sigma, of 10, 6, and 4 grams per squared centimeter, respectively. For sigma=10 and sigma=6 g/cm^2, respectively, these formation times are reduced by more than 50% and more than 80% compared with those in a previously published calculation with the old approximation to the opacity.

preprint2010arXiv

Gyroscopic pumping of large-scale flows in stellar interiors, and application to Lithium Dip stars

The maintenance of large-scale differential rotation in stellar convective regions by rotationally influenced convective stresses also drives large-scale meridional flows by angular--momentum conservation. This process is an example of ``gyroscopic pumping'', and has recently been studied in detail in the solar context. An important question concerns the extent to which these gyroscopically pumped meridional flows penetrate into nearby stably stratified (radiative) regions, since they could potentially be an important source of non-local mixing. Here we present an extensive study of the gyroscopic pumping mechanism, using a combination of analytical calculations and numerical simulations both in Cartesian geometry and in spherical geometry. The various methods, when compared with one another, provide physical insight into the process itself, as well as increasingly sophisticated means of estimating the gyroscopic pumping rate. As an example of application, we investigate the effects of this large-scale mixing process on the surface abundances of the light elements Li and Be for stars in the mass range 1.3-1.5 solar masses (so-called ``Li-dip stars''). We find that gyroscopic pumping is a very efficient mechanism for circulating material between the surface and the deep interior, so much in fact that it over-estimates Li and Be depletion by orders of magnitude for stars on the hot side of the dip.However, when the diffusion of chemical species back into the surface convection zone is taken into account, a good fit with observed surface abundances of Li and Be as a function of stellar mass in the Hyades cluster can be found for reasonable choices of model parameters.

preprint2010arXiv

Supermassive Dark Stars: Detectable in JWST

The first phase of stellar evolution in the history of the Universe may be Dark Stars, powered by dark matter heating rather than by nuclear fusion. Weakly Interacting Massive Particles, which may be their own antipartners, collect inside the first stars and annihilate to produce a heat source that can power the stars for millions to billions of years. In this paper we show that these objects can grow to be supermassive dark stars (SMDS) with masses $\gtrsim (10^5-10^7) \msun$. The growth continues as long as dark matter heating persists, since dark stars are large and cool (surface temperature $\lesssim 5\times 10^4$K) and do not emit enough ionizing photons to prevent further accretion of baryons onto the star. The dark matter may be provided by two mechanisms: (1) gravitational attraction of dark matter particles on a variety of orbits not previously considered, and (2) capture of WIMPs due to elastic scattering. Once the dark matter fuel is exhausted, the SMDS becomes a heavy main sequence star; these stars eventually collapse to form massive black holes that may provide seeds for supermassive black holes in the Universe. SMDS are very bright, with luminosities exceeding $(10^9-10^{11}) L_\odot$. We demonstrate that for several reasonable parameters, these objects will be detectable with JWST. Such an observational discovery would confirm the existence of a new phase of stellar evolution powered by dark matter.

preprint2010arXiv

The Effects of Metallicity, and Grain Growth and Settling on the Early Evolution of Gaseous Protoplanets

Giant protoplanets formed by gravitational instability in the outer regions of circumstellar disks go through an early phase of quasi-static contraction during which radii are large and internal temperatures are low. The main source of opacity in these objects is dust grains. We investigate two problems involving the effect of opacity on the evolution of planets of 3, 5, and 7 M_J. First, we pick three different overall metallicities for the planet and simply scale the opacity accordingly. We show that higher metallicity results in slower contraction as a result of higher opacity. It is found that the pre-collapse time scale is proportional to the metallicity. In this scenario, survival of giant planets formed by gravitational instability is predicted to be more likely around low-metallicity stars, since they evolve to the point of collapse to small size on shorter time scales. But metal-rich planets, as a result of longer contraction times, have the best opportunity to capture planetesimals and form heavy-element cores. Second, we investigate the effects of opacity reduction as a result of grain growth and settling, for the same three planetary masses and for three different values of overall metallicity. When these processes are included, the pre-collapse time scale is found to be of order 1000 years for the three masses, significantly shorter than the time scale calculated without these effects. In this case the time scale is found to be relatively insensitive to planetary mass and composition. However, the effects of planetary rotation and accretion of gas and dust, which could increase the timescale, are not included in the calculation. The short time scale we find would preclude metal enrichment by planetesimal capture, as well as heavy-element core formation, over a large range of planetary masses and metallicities.

preprint2009arXiv

Dark Stars: a new look at the First Stars in the Universe

We have proposed that the first phase of stellar evolution in the history of the Universe may be Dark Stars (DS), powered by dark matter heating rather than by nuclear fusion, and in this paper we examine the history of these DS. The power source is annihilation of Weakly Interacting Massive Particles (WIMPs) which are their own antiparticles. These WIMPs are the best motivated dark matter (DM) candidates and may be discovered by ongoing direct or indirect detection searches (e.g. FERMI /GLAST) or at the Large Hadron Collider at CERN. A new stellar phase results, powered by DM annihilation as long as there is DM fuel, from millions to billions of years. We build up the dark stars from the time DM heating becomes the dominant power source, accreting more and more matter onto them. We have included many new effects in the current study, including a variety of particle masses and accretion rates, nuclear burning, feedback mechanisms, and possible repopulation of DM density due to capture. Remarkably, we find that in all these cases, we obtain the same result: the first stars are very large, 500-1000 times as massive as the Sun; as well as puffy (radii 1-10 A.U.), bright ($10^6-10^7 L_\odot$), and cool ($T_{surf} < $10,000 K) during the accretion. These results differ markedly from the standard picture in the absence of DM heating. Hence DS should be observationally distinct from standard Pop III stars. In addition, DS avoid the (unobserved) element enrichment produced by the standard first stars. Once the dark matter fuel is exhausted, the DS becomes a heavy main sequence star; these stars eventually collapse to form massive black holes that may provide seeds for the supermassive black holes and intermediate black holes, and explain ARCADE data.

preprint2009arXiv

Dark Stars: A New Study of the FIrst Stars in the Universe

We have proposed that the first phase of stellar evolution in the history of the Universe may be Dark Stars (DS), powered by dark matter heating rather than by nuclear fusion. Weakly Interacting Massive Particles, which may be their own antipartners, collect inside the first stars and annihilate to produce a heat source that can power the stars. A new stellar phase results, a Dark Star, powered by dark matter annihilation as long as there is dark matter fuel, with lifetimes from millions to billions of years. We find that the first stars are very bright ($\sim 10^6 L_\odot$) and cool ($T_{surf} < 10,000$K) during the DS phase, and grow to be very massive (500-1000 times as massive as the Sun). These results differ markedly from the standard picture in the absence of DM heating, in which the maximum mass is about 140$M_\odot$ and the temperatures are much hotter ($T_{surf} > 50,000$K); hence DS should be observationally distinct from standard Pop III stars. Once the dark matter fuel is exhausted, the DS becomes a heavy main sequence star; these stars eventually collapse to form massive black holes that may provide seeds for supermassive black holes observed at early times as well as explanations for recent ARCADE data and for intermediate black holes.

preprint2009arXiv

Dark Stars: Begynnelsen

The first phase of stellar evolution in the history of the universe may be Dark Stars, powered by dark matter heating rather than by fusion. Weakly interacting massive particles, which are their own antiparticles, can annihilate and provide an important heat source for the first stars in the universe. This and the following contribution present the story of Dark Stars. In this first part, we describe the conditions under which dark stars form in the early universe: 1) high dark matter densities, 2) the annihilation products get stuck inside the star, and 3) dark matter heating wins over all other cooling or heating mechanisms.

preprint2009arXiv

Dark Stars: Död och Återuppståndelse

The first phase of stellar evolution in the history of the universe may be Dark Stars, powered by dark matter heating rather than by fusion. Weakly interacting massive particles, which are their own antiparticles, can annihilate and provide an important heat source for the first stars in the universe. This and the previous contribution present the story of Dark Stars. In this second part, we describe the structure of Dark Stars and predict that they are very massive ($\sim 800 M_\odot$), cool (6000 K), bright ($\sim 10^6 L_\odot$), long-lived ($\sim 10^6$ years), and probable precursors to (otherwise unexplained) supermassive black holes. Later, once the initial dark matter fuel runs out and fusion sets in, dark matter annihilation can predominate again if the scattering cross section is strong enough, so that a Dark Star is born again.

preprint2009arXiv

Determination of the Interior Structure of Transiting Planets in Multiple-Planet Systems

Tidal dissipation within a short-period transiting extrasolar planet perturbed by a companion object can drive orbital evolution of the system to a so-called tidal fixed point, in which the apsidal lines of the transiting planet and its perturber are aligned, and for which variations in the orbital eccentricities of both planet and perturber are damped out. Significant contributions to the apsidal precession rate are made by the secular planet-planet interaction, by general relativity, and by the gravitational quadropole fields created by the transiting planet's tidal and rotational distortions. The fixed-point orbital eccentricity of the inner planet is therefore a strong function of the planet's interior structure. We illustrate these ideas in the specific context of the recently discovered HAT-P-13 exo-planetary system, and show that one can already glean important insights into the physical properties of the inner transiting planet. We present structural models of the planet, which indicate that its observed radius can be maintained for a one-parameter sequence of models that properly vary core mass and tidal energy dissipation in the interior. We use an octopole-order secular theory of the orbital dynamics to derive the dependence of the inner planet's eccentricity, on its tidal Love number. We find that the currently measured eccentricity, implies 0.116 < k2_{b} < 0.425, 0 M_{Earth}<M_{core}<120 M_{Earth}$, and Q_{b} < 300,000. Improved measurement of the eccentricity will soon allow for far tighter limits to be placed on all three of these quantities, and will provide an unprecedented probe into the interior structure of an extrasolar planet.

preprint2009arXiv

The Formation of Uranus and Neptune in Solid-Rich Feeding Zones: Connecting Chemistry and Dynamics

The core accretion theory of planet formation has at least two fundamental problems explaining the origins of Uranus and Neptune: (1) dynamical times in the trans-Saturnian solar nebula are so long that core growth can take > 15 Myr, and (2) the onset of runaway gas accretion that begins when cores reach 10 Earth masses necessitates a sudden gas accretion cutoff just as the ice giant cores reach critical mass. Both problems may be resolved by allowing the ice giants to migrate outward after their formation in solid-rich feeding zones with planetesimal surface densities well above the minimum-mass solar nebula. We present new simulations of the formation of Uranus and Neptune in the solid-rich disk of Dodson-Robinson et al. (2009) using the initial semimajor axis distribution of the Nice model (Gomes et al. 2005; Morbidelli et al. 2005; Tsiganis et al. 2005), with one ice giant forming at 12 AU and the other at 15 AU. The innermost ice giant reaches its present mass after 3.8-4.0 Myr and the outermost after 5.3-6 Myr, a considerable time decrease from previous one-dimensional simulations (e.g. Pollack et al. 1996). The core masses stay subcritical, eliminating the need for a sudden gas accretion cutoff. Our calculated carbon mass fractions of 22% are in excellent agreement with the ice giant interior models of Podolak et al. (1995) and Marley et al. (1995). Based on the requirement that the ice giant-forming planetesimals contain >10% mass fractions of methane ice, we can reject any solar system formation model that initially places Uranus and Neptune inside the orbit of Saturn. We also demonstrate that a large population of planetesimals must be present in both ice giant feeding zones throughout the lifetime of the gaseous nebula.

preprint2009arXiv

Tidal and Magnetic Interactions between a Hot Jupiter and its Host Star in the Magnetospheric Cavity of a Protoplanetary Disk

We present a simplified model to study the orbital evolution of a young hot Jupiter inside the magnetospheric cavity of a proto-planetary disk. The model takes into account the disk locking of stellar spin as well as the tidal and magnetic interactions between the star and the planet. We focus on the orbital evolution starting from the orbit in the 2:1 resonance with the inner edge of the disk, followed by the inward and then outward orbital migration driven by the tidal and magnetic torques as well as the Roche-lobe overflow of the tidally inflated planet. The goal in this paper is to study how the orbital evolution inside the magnetospheric cavity depends on the cavity size, planet mass, and orbital eccentricity. In the present work, we only target the mass range from 0.7 to 2 Jupiter masses. In the case of the large cavity corresponding to the rotational period ~ 7 days, the planet of mass >1 Jupiter mass with moderate initial eccentricities (>~ 0.3) can move to the region < 0.03 AU from its central star in 10^7 years, while the planet of mass <1 Jupiter mass cannot. We estimate the critical eccentricity beyond which the planet of a given mass will overflow its Roche radius and finally lose all of its gas onto the star due to runaway mass loss. In the case of the small cavity corresponding to the rotational period ~ 3 days, all of the simulated planets lose all of their gas even in circular orbits. Our results for the orbital evolution of young hot Jupiters may have the potential to explain the absence of low-mass giant planets inside ~ 0.03 AU from their dwarf stars revealed by transit surveys.

preprint2008arXiv

Dark Stars: the First Stars in the Universe may be powered by Dark Matter Heating

A new line of research on Dark Stars is reviewed, which suggests that the first stars to exist in the universe were powered by dark matter heating rather than by fusion. Weakly Interacting Massive Particles, which may be there own antipartmers, collect inside the first stars and annihilate to produce a heat source that can power the stars. A new stellar phase results, a Dark Star, powered by dark matter annihilation as long as there is dark matter fuel.

preprint2008arXiv

Models of Jupiter's Growth Incorporating Thermal and Hydrodynamic Constraints

[Abridged] We model the growth of Jupiter via core nucleated accretion, applying constraints from hydrodynamical processes that result from the disk-planet interaction. We compute the planet's internal structure using a Henyey-type stellar evolution code. The planet's interactions with the protoplanetary disk are calculated using 3-D hydrodynamic simulations. Previous models of Jupiter's growth have taken the radius of the planet to be approximately one Hill sphere radius, Rhill. However, 3-D hydrodynamic simulations show that only gas within 0.25Rhill remains bound to the planet, with the more distant gas eventually participating in the shear flow of the protoplanetary disk. Therefore in our new simulations, the planet's outer boundary is placed at the location where gas has the thermal energy to reach the portion of the flow not bound to the planet. We find that the smaller radius increases the time required for planetary growth by ~5%. Thermal pressure limits the rate at which a planet less than a few dozen times as massive as Earth can accumulate gas from the protoplanetary disk, whereas hydrodynamics regulates the growth rate for more massive planets. Within a moderately viscous disk, the accretion rate peaks when the planet's mass is about equal to the mass of Saturn. In a less viscous disk hydrodynamical limits to accretion are smaller, and the accretion rate peaks at lower mass. To account for disk dissipation, we perform some of our simulations of Jupiter's growth within a disk whose surface gas density decreases on a timescale of 3Myr. According to our simulations, proto-Jupiter's distended and thermally-supported envelope was too small to capture the planet's current retinue of irregular satellites as advocated by Pollack et al. (1979).

preprint2008arXiv

Stellar Structure of Dark Stars: a first phase of Stellar Evolution due to Dark Matter Annihilation

Dark Stars are the very first phase of stellar evolution in the history of the universe: the first stars to form (typically at redshifts $z \sim 10-50$) are powered by heating from dark matter (DM) annihilation instead of fusion (if the DM is made of particles which are their own antiparticles). We find equilibrium polytropic configurations for these stars; we start from the time DM heating becomes important ($M \sim 1-10 M_\odot$) and build up the star via accretion up to 1000 M$_\odot$. The dark stars, with an assumed particle mass of 100 GeV, are found to have luminosities of a few times $10^6$ L$_\odot$, surface temperatures of 4000--10,000 K, radii $\sim 10^{14}$ cm, lifetimes of at least $ 0.5$ Myr, and are predicted to show lines of atomic and molecular hydrogen. Dark stars look quite different from standard metal-free stars without DM heating: they are far more massive (e.g. $\sim 800 M_\odot$ for 100 GeV WIMPs), cooler, and larger, and can be distinguished in future observations, possibly even by JWST or TMT.

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.

preprint2003arXiv

On the Surface Heating of Synchronously-Spinning Short-Period Jovian Planets

We consider the atmospheric flow on short-period extra-solar planets through two-dimensional numerical simulations of hydrodynamics with radiation transfer. One side is always exposed to the irradiation from the host star. The other is always in shadow. The temperature of the day side is determined by the equilibrium which the planetary atmosphere establishes with stellar radiation. Part of the thermal energy deposited on the day side is advected to the night side by a current. The radiation transfer, the night-side temperature distribution and by this the spectroscopic signature of the planet are sensitive functions of the atmospheric opacity. If the atmosphere contains grains with an abundanceand size distribution comparable to that of the interstellar medium, shallow heating occurs on the day side and the night side cools well below the day side. The temperature difference decreases as the abundance of grains is reduced. A simple analytic model of the dissipation of the circulation flow and associated kinetic heating is considered. This heating effect occurs mostly near the photosphere, not deep enough to significantly affect the size of planets. We show that the surface irradiation suppresses convection near the photospheric region on the day side. In some cases convection zones appear near the surface on the night side. This structural modification may influence the response and dissipation of tidal disturbances and alter the circularization and synchronization time scales.