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.

Formation of the planets in the Kepler-36 system is modeled by detailed numerical simulations according to the core-nucleated accretion scenario. The standard model is updated to include the dissolution of accreting rocky planetesimals in the gaseous envelope of the planet, leading to substantial enrichment of the envelope mass in heavy elements and a non-uniform composition with depth. For Kepler-36 c, models involving in situ formation and models involving orbital migration are considered. The results are compared with standard formation models. The calculations include the formation (accretion) phase, as well as the subsequent cooling phase, up to the age of Kepler-36 (7 Gyr). During the latter phase, mass loss induced by stellar XUV radiation is included. In all cases, the results fit the measured mass, 7.84 M$_\oplus$, and radius, 3.68 R$_\oplus$, of Kepler-36 c. Two parameters are varied to obtain these fits: the disk solid surface density at the formation location, and the "efficiency" factor in the XUV mass loss rate. The updated models are hotter and therefore less dense in the silicate portion of the planet and in the overlying layers of H/He, as compared with standard models. The lower densities mean that only about half as much H/He is needed to be accreted to fit the present-day mass and radius constraints. For Kepler-36 b, an updated in situ calculation shows that the entire H/He envelope is lost, early in the cooling phase, in agreement with observation.

The composition of Jupiter and the primordial distribution of the heavy elements are determined by its formation history. As a result, in order to constrain the primordial internal structure of Jupiter the growth of the core and the deposition and settling of accreted planetesimals must be followed in detail. In this paper we determine the distribution of the heavy elements in proto-Jupiter and determine the mass and composition of the core. We find that while the outer envelope of proto-Jupiter is typically convective and has an homogeneous composition, the innermost regions have compositional gradients. In addition, the existence of heavy elements in the envelope leads to much higher internal temperatures (several times 10$^4$ K) than in the case of a hydrogen-helium envelope. The derived core mass depends on the actual definition of the core: if the core is defined as the region in which the heavy-element mass fraction is above some limit (say 0.5), then it can be much more massive ($\sim$ 15 M$_{\oplus}$) and more extended (10\% of the planet's radius) than in the case where the core is just the region with 100\% heavy elements. In the former case Jupiter's core also consists of hydrogen and helium. Our results should be taken into account when constructing internal structure models of Jupiter and when interpreting the upcoming data from the {\it Juno} (NASA) mission.

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.

Hot Jupiters, giant extrasolar planets with orbital periods shorter than ~10 days, have long been thought to form at large radial distances, only to subsequently experience long-range inward migration. Here, we propose that in contrast with this picture, a substantial fraction of the hot Jupiter population formed in situ via the core accretion process. We show that under conditions appropriate to the inner regions of protoplanetary disks, rapid gas accretion can be initiated by Super-Earth type planets, comprising 10-20 Earth masses of refractory composition material. An in situ formation scenario leads to testable consequences, including the expectation that hot Jupiters should frequently be accompanied by additional low-mass planets with periods shorter than ~100 days. Our calculations further demonstrate that dynamical interactions during the early stages of planetary systems' lifetimes should increase the inclinations of such companions, rendering transits rare. High-precision radial velocity monitoring provides the best prospect for their detection.

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.

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.

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.

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.

The large number of detected giant exoplanets offers the opportunity to improve our understanding of the formation mechanism, evolution, and interior structure of gas giant planets. The two main models for giant planet formation are core accretion and disk instability. There are substantial differences between these formation models, including formation timescale, favorable formation location, ideal disk properties for planetary formation, early evolution, planetary composition, etc. First, we summarize the two models including their substantial differences, advantages, and disadvantages, and suggest how theoretical models should be connected to available (and future) data. We next summarize current knowledge of the internal structures of solar- and extrasolar- giant planets. Finally, we suggest the next steps to be taken in giant planet exploration.

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.

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

We present calculations of thermal evolution of Hot Jupiters with various masses and effective temperatures under Ohmic dissipation. The resulting evolutionary sequences show a clear tendency towards inflated radii for effective temperatures that give rise to significant ionization of alkali metals in the atmosphere, compatible with the trend of the data. The degree of inflation shows that Ohmic dissipation, along with the likely variability in heavy element content can account for all of the currently detected radius anomalies. Furthermore, we find that in absence of a massive core, low-mass hot Jupiters can over-flow their Roche-lobes and evaporate on Gyr time-scales, possibly leaving behind small rocky cores.

We revisit the calculation of the Ohmic dissipation in a hot Jupiter presented in Laine et al. (2008) by considering more realistic interior structures, stellar obliquity, and the resulting orbital evolution. In this simplified approach, the young hot Jupiter of one Jupiter mass is modelled as a diamagnetic sphere with a finite resistivity, orbiting across tilted stellar magnetic dipole fields in vacuum. Since the induced Ohmic dissipation occurs mostly near the planet's surface, we find that the dissipation is unable to significantly expand the young hot Jupiter. Nevertheless, the planet inside a small co-rotation orbital radius can undergo orbital decay by the dissipation torque and finally overfill its Roche lobe during the T Tauri star phase. The stellar obliquity can evolve significantly if the magnetic dipole is parallel/anti-parallel to the stellar spin. Our results are validated by the general torque-dissipation relation in the presence of the stellar obliquity. We also run the fiducial model in Laine et al. (2008) and find that the planet's radius is sustained at a nearly constant value by the Ohmic heating, rather than being thermally expanded to the Roche radius as suggested by the authors.

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.

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.

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.

The first stars to form in the history of the universe may have been powered by dark matter annihilation rather than by fusion. This new phase of stellar evolution may have lasted millions to billions of years. These dark stars can grow to be very large, > 10^5 solar masses, and are relatively cool (~10^4 K). They are also very bright, being potentially detectable in the upcoming James Webb Space Telescope or even the Hubble Space Telescope. Once the dark matter runs out, the dark stars have a short fusion phase, before collapsing into black holes (BH). The resulting BH could serve as seeds for the (unexplained) supermassive black holes at high redshift and at the centers of galaxies.

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.

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.

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.

The final composition of giant planets formed as a result of gravitational instability in the disk gas depends on their ability to capture solid material (planetesimals) during their 'pre-collapse' stage, when they are extended and cold, and contracting quasi-statically. The duration of the pre-collapse stage is inversely proportional roughly to the square of the planetary mass, so massive protoplanets have shorter pre-collapse timescales and therefore limited opportunity for planetesimal capture. The available accretion time for protoplanets with masses of 3, 5, 7, and 10 Jupiter masses is found to be 7.82E4, 2.62E4, 1.17E4 and 5.67E3 years, respectively. The total mass that can be captured by the protoplanets depends on the planetary mass, planetesimal size, the radial distance of the protoplanet from the parent star, and the local solid surface density. We consider three radial distances, 24, 38, and 68 AU, similar to the radial distances of the planets in the system HR 8799, and estimate the mass of heavy elements that can be accreted. We find that for the planetary masses usually adopted for the HR 8799 system, the amount of heavy elements accreted by the planets is small, leaving them with nearly stellar compositions.

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.

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.

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.

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.

The transiting "hot Saturn" HD 149026b, which has the highest mean density of any confirmed planet in the Neptune-Jupiter mass range, has challenged theories of planet formation since its discovery in 2005. Previous investigations could not explain the origin of the planet's 45-110 Earth-mass solid core without invoking catastrophes such as gas giant collisions or heavy planetesimal bombardment launched by neighboring planets. Here we show that HD 149026b's large core can be successfully explained by the standard core accretion theory of planet formation. The keys to our reconstruction of HD 149026b are (1) applying a model of the solar nebula to describe the protoplanet nursery; (2) placing the planet initially on a long-period orbit at Saturn's heliocentric distance of 9.5 AU; and (3) adjusting the solid mass in the HD 149026 disk to twice that of the solar nebula in accordance with the star's heavy element enrichment. We show that the planet's migration into its current orbit at 0.042 AU is consistent with our formation model. Our study of HD 149026b demonstrates that it is possible to discover the growth history of any planet with a well-defined core mass that orbits a solar-type star.

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.

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.

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.

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.

To date, there is no core accretion simulation that can successfully account for the formation of Uranus or Neptune within the observed 2-3 Myr lifetimes of protoplanetary disks. Since solid accretion rate is directly proportional to the available planetesimal surface density, one way to speed up planet formation is to take a full accounting of all the planetesimal-forming solids present in the solar nebula. By combining a viscously evolving protostellar disk with a kinetic model of ice formation, we calculate the solid surface density in the solar nebula as a function of heliocentric distance and time. We find three effects that strongly favor giant planet formation: (1) a decretion flow that brings mass from the inner solar nebula to the giant planet-forming region, (2) recent lab results (Collings et al. 2004) showing that the ammonia and water ice lines should coincide, and (3) the presence of a substantial amount of methane ice in the trans-Saturnian region. Our results show higher solid surface densities than assumed in the core accretion models of Pollack et al. (1996) by a factor of 3 to 4 throughout the trans-Saturnian region. We also discuss the location of ice lines and their movement through the solar nebula, and provide new constraints on the possible initial disk configurations from gravitational stability arguments.

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

We present two new in situ core accretion simulations of Saturn with planet formation timescales of 3.37 Myr (model S0) and 3.48 Myr (model S1), consistent with observed protostellar disk lifetimes. In model S0, we assume rapid grain settling reduces opacity due to grains from full interstellar values (Podolak 2003). In model S1, we do not invoke grain settling, instead assigning full interstellar opacities to grains in the envelope. Surprisingly, the two models produce nearly identical formation timescales and core/atmosphere mass ratios. We therefore observe a new manifestation of core accretion theory: at large heliocentric distances, the solid core growth rate (limited by Keplerian orbital velocity) controls the planet formation timescale. We argue that this paradigm should apply to Uranus and Neptune as well.

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.

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 the universe. This talk presents the story of these Dark Stars. We make predictions that the first stars 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 DM fuel runs out and fusion sets in, DM annihilation can predominate again if the scattering cross section is strong enough, so that a Dark Star is born again.

Traditional thermal evolution models of giant planets employ arbitrary initial conditions selected more for computational expediency than physical accuracy. Since the initial conditions are eventually forgotten by the evolving planet, this approach is valid for mature planets, if not young ones. To explore the evolution at young ages of jovian mass planets we have employed model planets created by one implementation of the core accretion mechanism as initial conditions for evolutionary calculations. We find that young jovian planets are smaller, cooler, and several to 100 times less luminous than predicted by earlier models. Furthermore the time interval during which the young jupiters are fainter than expected depends on the mass of planet. Jupiter mass planets (1 M_J) align with the conventional model luminosity in as little at 20 million years, but 10 M_J planets can take up to 1 billion years to match commonly cited luminosities, given our implementation of the core accretion mechanism. If our assumptions, especially including our treatment of the accretion shock, are correct, then young jovian planets are substantially fainter at young ages than currently believed. These results have important consequences both for detection strategies and for assigning masses to young jovian planets based on observed luminosities.

We present evidence that stars with planets exhibit statistically significant silicon and nickel enrichment over the general metal-rich population. We also present simulations which predict silicon enhancement of planet hosts within the context of the core-accretion hypothesis for giant planet formation. Because silicon and oxygen are both alpha elements, [Si/Fe] traces [O/Fe], so the silicon enhancement in planet hosts predicts that these stars are oxygen-rich as well. We present new numerical simulations of planet formation by core accretion that establish the timescale on which a Jovian planet reaches rapid gas accretion, t_rga, as a function of solid surface density sigma_solid: (t_rga / 1 Myr) = (sigma_solid / 25.0 g cm^{-2})^{-1.44}. This relation enables us to construct Monte Carlo simulations that predict the fraction of star-disk systems that form planets as a function of [Fe/H], [Si/Fe], disk mass, outer disk radius and disk lifetime. Our simulations reproduce both the known planet-metallicity correlation and the planet-silicon correlation reported in this paper. The simulations predict that 16% of Solar-type stars form Jupiter-mass planets, in agreement with 12% predicted from extrapolation of the observed planet frequency-semimajor axis distribution. Although a simple interpretation of core accretion predicts that the planet-silicon correlation should be much stronger than the planet-nickel correlation, we observe the same degree of silicon and nickel enhancement in planet hosts. If this result persists once more planets have been discovered, it might indicate a complexity in the chemistry of planet formation beyond the simple accumulation of solids in the core accretion theory.

We study the formation of a giant gas planet by the core--accretion gas--capture process, with numerical simulations, under the assumption that the planetary core forms in the center of an anti-cyclonic vortex. The presence of the vortex concentrates particles of centimeter to meter size from the surrounding disk, and speeds up the core formation process. Assuming that a planet of Jupiter mass is forming at 5 AU from the star, the vortex enhancement results in considerably shorter formation times than are found in standard core--accretion gas--capture simulations. Also, formation of a gas giant is possible in a disk with mass comparable to that of the minimum mass solar nebula.

Here we show preliminary calculations of the cooling and contraction of a 2 MJ planet. These calculations, which are being extended to 1-10 MJ, differ from other published "cooling tracks" in that they include a core accretion-gas capture formation scenario, the leading theory for the formation of gas giant planets. We find that the initial post-accretionary intrinsic luminosity of the planet is ~3 times less than previously published models which use arbitrary initial conditions. These differences last a few tens of millions of years. Young giant planets are intrinsically fainter than has been previously appreciated. We also discuss how uncertainties in atmospheric chemistry and the duration of the formation time of giant planets lead to challenges in deriving planetary physical properties from comparison with tabulated model values.

Doppler measurements from Subaru and Keck have revealed radial velocity variations in the V=8.15, G0IV star HD 149026 consistent with a Saturn-Mass planet in a 2.8766 day orbit. Photometric observations at Fairborn Observatory have detected three complete transit events with depths of 0.003 mag at the predicted times of conjunction. HD 149026 is now the second brightest star with a transiting extrasolar planet. The mass of the star, based on interpolation of stellar evolutionary models, is 1.3 +/- 0.1 solar masses; together with the Doppler amplitude, K=43.3 m s^-1, we derive a planet mass Msin(i)=0.36 Mjup, and orbital radius of 0.042 AU. HD 149026 is chromospherically inactive and metal-rich with spectroscopically derived [Fe/H]=+0.36, Teff=6147 K, log g=4.26 and vsin(i)=6.0 km s^-1. Based on Teff and the stellar luminosity of 2.72 Lsun, we derive a stellar radius of 1.45 Rsun. Modeling of the three photometric transits provides an orbital inclination of 85.3 +/- 1.0 degrees and (including the uncertainty in the stellar radius) a planet radius of 0.725 +/- 0.05 Rjup. Models for this planet mass and radius suggest the presence of a ~67 Mearth core composed of elements heavier than hydrogen and helium. This substantial planet core would be difficult to construct by gravitational instability.

We present the first results from simulations of processes leading to planet formation in protoplanetary disks with different metallicities. For a given metallicity, we construct a two-dimensional grid of disk models with different initial masses and radii ($M_0$, $R_0$). For each disk, we follow the evolution of gas and solids from an early evolutionary stage, when all solids are in the form of small dust grains, to the stage when most solids have condensed into planetesimals. Then, based on the core accretion - gas capture scenario, we estimate the planet-bearing capability of the environment defined by the final planetesimal swarm and the still evolving gaseous component of the disk. We define the probability of planet-formation, $P_p$, as the normalized fractional area in the ($M_0$, $\log R_0$) plane populated by disks that have formed planets inside 5 AU. With such a definition, and under the assumption that the population of planets discovered at $R$ $<$ 5 AU is not significantly contaminated by planets that have migrated from $R$ $>$ 5 AU, our results agree fairly well with the observed dependence between the probability that a star harbors a planet and the star's metal content. The agreement holds for the disk viscosity parameter $\alpha$ ranging from $10^{-3}$ to $10^{-2}$, and it becomes much poorer when the redistribution of solids relative to the gas is not allowed for during the evolution of model disks.

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.

Several short-period Jupiter-mass planets have been discovered around nearby solar-type stars. During the circularization of their orbits, the dissipation of tidal disturbance by their host stars heats the interior and inflates the sizes of these planets. Based on a series of internal structure calculations for giant planets, we examine the physical processes which determine their luminosity-radius relation. In order for young or intensely heated gas giant planets to attain quasi-hydrostatic equilibria, with sizes comparable to or larger than two Jupiter radii, their interiors must have sufficiently high temperature and low density such that degeneracy effects are relatively weak compared to those in a mature and compact Jupiter. Consequently, the polytropic index monotonically increases whereas the central temperature increases and then decreases with the planets' size. These effects, along with a temperature-sensitive opacity for the radiative surface layers of giant planets, cause the power index of the luminosity's dependence on radius to decrease with increasing radius. For planets larger than twice Jupiter's radius, this index is sufficiently small that they become unstable to tidal inflation. We make comparisons between cases of uniform heating and cases in which the heating is concentrated in various locations within the giant planet. Based on these results we suggest that accurate measurement of the sizes of close-in young Jupiters can be used to probe their internal structure under the influence of tidal heating.

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.

We have computed evolutionary models for extrasolar planets which range in mass from 0.1 to 3.0 Jovian Masses, and which range in equilibrium temperature from 113 K to 2000 K. We present four sequences of models, designed to show the structural effects of a solid core and of internal heating due to the conversion of kinetic to thermal energy at pressures of tens of bars. The model planetary radii are intended for comparisons with radii derived from observations of transiting extrasolar planets. To provide such comparisons, we expect that of order 10 transiting planets with orbital periods less than 200 days can be detected around bright (V<10) main-sequence stars for which accurate, well-sampled radial velocity measurements can be readily accumulated. Through these observations, structural properties of the planets will be derivable, particularly for low-mass, high-temperature planets. Implications regarding the transiting companion to OGLE-TR-56 recently announced by Konacki et al. are discussed. With regard to the confirmed transiting planet, HD 209458b, we find, in accordance with other recent calculations, that models without internal heating predict a radius that is 35 percent smaller than the observed radius. We explore the possibility that HD 209458b owes its large size to dissipation of energy arising from ongoing tidal circularization of the orbit. We show that residual scatter in the current radial velocity data set for HD 209458b is consistent with the presence of an as-of-yet undetected second companion, and that further radial velocity monitoring of the star is indicated.

We investigate the possibility of substantial inflation of short-period Jupiter-mass planets, as a result of their internal tidal dissipation associated with the synchronization and circularization of their orbits. We employ the simplest prescription based on an equilibrium model with a constant lag angle for all components of the tide. We show that for young Jupiter-mass planets, with a period less than 3 days, an initial radius about 2 Jupiter radii, and an orbital eccentricity greater than 0.2, the energy dissipated during the circularization of their orbits is sufficiently intense and protracted to inflate their sizes up to their Roche radii.

In this paper we present the global baroclinic instability as a source for vigorous turbulence leading to angular momentum transport in Keplerian accretion disks. We show by analytical considerations and three-dimensional radiation hydro simulations that, in particular, protoplanetary disks have a negative radial entropy gradient, which makes them baroclinic. Two-dimensional numerical simulations show that a baroclinic flow is unstable and produces turbulence. These findings are tested for numerical effects by performing a simulation with a barotropic initial condition which shows that imposed turbulence rapidly decays. The turbulence in baroclinic disks transports angular momentum outward and creates a radially inward bound accretion of matter. Potential energy is released and excess kinetic energy is dissipated. Finally the reheating of the gas supports the radial entropy gradient, forming a self consistent process. We measure accretion rates in our 2D and 3D simulations of dotM= - 1E-9 - 1E-7 Msun/yr and viscosity parameters of alpha = 1E-4 - 1E-2, which fit perfectly together and agree reasonably with observations. The turbulence creates pressure waves, Rossby waves, and vortices in the R-phi plane of the disk. We demonstrate in a global simulation that these vortices tend to form out of little background noise and to be long-lasting features, which have already been suggested to lead to the formation of planets.

Formation of planets in the 47 UMa system is followed in an evolving protoplanetary disk composed of gas and solids. The evolution of the disk is calculated from an early stage, when all solids, assumed to be high-temperature silicates, are in the dust form, to the stage when most solids are locked in planetesimals. The simulation of planetary evolution starts with a solid embryo of ~1 Earth mass, and proceeds according to the core accretion -- gas capture model. Orbital parameters are kept constant, and it is assumed that the environment of each planet is not perturbed by the second planet. It is found that conditions suitable for both planets to form within several Myr are easily created, and maintained throughout the formation time, in disks with $\alpha \approx 0.01$. In such disks, a planet of 2.6 Jupiter masses (the minimum for the inner planet of the 47 UMa system) may be formed at 2.1 AU from the star in \~3 Myr, while a planet of 0.89 Jupiter masses (the minimum for the outer planet) may be formed at 3.95 AU from the star in about the same time. The formation of planets is possible as a result of a significant enhancement of the surface density of solids between 1.0 and 4.0 AU, which results from the evolution of a disk with an initially uniform gas-to-dust ratio of 167 and an initial radius of 40 AU.

The rotational properties of centrally condensed, turbulent molecular cloud cores with velocity fields that are characterized by Gaussian random fields are investigated. It is shown that the observed line width-size relationship can be reproduced if the velocity power spectrum is a power-law with P(k)=k**n and n=-3 to -4. The line-of-sight velocity maps of these cores show velocity gradients that can be interpreted as rotation. For n=-4, both, the deduced values of the angular velocity Omega=1.6 km/s/pc * (R/0.1 pc)**0.5 and the scaling relations between Omega and the core radius R are in very good agreement with the observations. As a result of the dominance of long wavelength modes, the cores also have a net specific angular momentum with an average value of j=7*(10**20)*(R/0.1 pc)**(1.5) cm**2/s with a large spread. Their internal dimensionless rotational parameter is beta=0.03, independent of the scale radius R. In general, the line-of-sight velocity gradient of an individual turbulent core does not provide a good estimate of its internal specific angular momentum. We find however that the distribution of the specific angular momenta of a large sample of cores which are described by the same power spectrum can be determined very accurately from the distribution of their line-of-sight velocity gradients Omega using the simple formula j=p*Omega*R*R where p depends on the density distribution of the core and has to be determined from a Monte-Carlo study. Our results show that for centrally condensed cores the intrinsic angular momentum is overestimated by a factor of 2-3 if p=0.4 is used.

The evolution of binary systems consisting of an asymptotic giant branch star of mass equal to 3 M_sun or 5 M_sun, and a main sequence star of mass equal to 0.4 M_sun or 0.6 M_sun with orbital periods > 200 days has been followed from the onset through the late stages of the common envelope phase. Using a nested grid technique, the three-dimensional hydrodynamical simulations of an asymptotic giant branch star with radii approximately 1 A.U. indicate that a significant fraction of the envelope gas is unbound (about 31% and 23% for binaries of 3 M_sun and 0.4 M_sun, and 5 M_sun and 0.6 M_sun respectively) by the ends of the simulations, and that the efficiency of the mass ejection process is about 40%. While the original volume of the giant is virtually evacuated in the late stages, most of the envelope gas remains marginally bound on the grid. At the ends of our simulations, when the orbital decay timescale exceeds about 5 years, the giant core and companion orbit one another with a period of about 1 day (2.4 days for a binary involving a more evolved giant), although this is an upper limit to the final orbital period. For a binary of 5 M_sun and 0.4 M_sun, the common envelope may not be completely ejected.

Hydrodynamical calculations in three space dimensions of the collapse of an isothermal, rotating 1 M\sol protostellar cloud are presented. The initial density stratification is a power law with density $\rho \propto r^{-p}$, with $p=1$. The case of the singular isothermal sphere ($p=2$) is not considered; however $p=1$ has been shown observationally to be a good representation of the density distribution in molecular cloud cores just before the beginning of collapse. The collapse is studied with two independent numerical methods, an SPH code with 200,000 particles, and a finite-difference code with nested grids which give high spatial resolution in the inner regions. Although previous numerical studies have indicated that such a power-law distribution would not result in fragmentation into a binary system, both codes show, in contrast, that multiple fragmentation does occur in the central regions of the protostar. Thus the process of binary formation by fragmentation is shown to be consistent with the fact that a large fraction of young stars are observed to be in binary or multiple systems.

Hydrodynamical calculations in three space dimensions of the collapse of an isothermal, centrally condensed, rotating 1 M\sol protostellar cloud are presented. A numerical algorithm involving nested subgrids is used to resolve the region where fragmentation occurs in the central part of the protostar. A previous calculation by Boss, which produced a hierarchical multiple system, is evolved further, at comparable numerical resolution, and the end result is a binary, with more than half of the mass of the original cloud, whose orbital separation increases with time as a result of accretion of high-angular momentum material and as a result of merging with fragments that have formed farther out. Repeating the calculation with significantly higher resolution, we find that a sequence of binaries can be induced by fragmentation of circumbinary disks. The stability of the resulting multiple system is investigated Using $n$-body calculations, which indicate that such a system would transform on a short time scale into a more stable hierarchical structure. The outermost and most massive binary which forms in the high-resolution run has properties similar to that of the binary found in the low-resolution calculation. Thus the basic outcome is shown to be independent of the numerical spatial resolution. The high-resolution run, in addition, leads to the formation of a system of smaller fragments, which might be important for the understanding of the origin of close binaries with low-mass components and of low mass single stars.