We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA star. Improvements in MESA star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 M ? stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA star solves the fully coupled stellar structure and composition equations, and we show how this has improved the scaling of MESA's calculational speed on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit that packages all the required components needed to form a unified, maintained, and well-validated build environment for MESA. We also highlight a few tools developed by the community for rapid visualization of MESA star results.

/pdf/modules-for-experiments-in-stellar-astrophysics-mesa-planets-29j324a9bq.pdf

Modules for Experiments in Stellar Astrophysics (MESA): Planets, Oscillations, Rotation, and Massive Stars

A collection of three papers in this issue, tackling seemingly unrelated planetary phenomena, marks a notable unification of Solar System dynamics. The three problems covered are the hard-to-explain orbits of giant planets, the evolution of the orbits of Jupiter's Trojan asteroids, and the cause of the ‘Late Heavy Bombardment’ that peppered the Moon with meteors, comets and asteroids some 700 million years after the planets were formed. Key to all these events, on this new model, was a rapid migration of the giant planets (Saturn, Jupiter, Neptune and Uranus) after a long period of stability within the Solar System. Planetary formation theories1,2 suggest that the giant planets formed on circular and coplanar orbits. The eccentricities of Jupiter, Saturn and Uranus, however, reach values of 6 per cent, 9 per cent and 8 per cent, respectively. In addition, the inclinations of the orbital planes of Saturn, Uranus and Neptune take maximum values of ∼2 degrees with respect to the mean orbital plane of Jupiter. Existing models for the excitation of the eccentricity of extrasolar giant planets3,4,5 have not been successfully applied to the Solar System. Here we show that a planetary system with initial quasi-circular, coplanar orbits would have evolved to the current orbital configuration, provided that Jupiter and Saturn crossed their 1:2 orbital resonance. We show that this resonance crossing could have occurred as the giant planets migrated owing to their interaction with a disk of planetesimals6,7. Our model reproduces all the important characteristics of the giant planets' orbits, namely their final semimajor axes, eccentricities and mutual inclinations.

http://www.ifa.hawaii.edu/~meech/a740/2006/spring/papers/nature03539.pdf

Origin of the orbital architecture of the giant planets of the Solar System.

Stellar multiplicity is an ubiquitous outcome of the star formation process Characterizing the frequency and main characteristics of multiple systems and their dependencies on primary mass and environment is therefore a powerful tool to probe this process While early attempts were fraught with selection biases and limited completeness, instrumentation breakthroughs in the last two decades now enable robust analyses In this review, we summarize our current empirical knowledge of stellar multiplicity for Main Sequence stars and brown dwarfs, as well as among populations of Pre-Main Sequence stars and embedded protostars Clear trends as a function of both primary mass and stellar evolutionary stage are identified that will serve as a comparison basis for numerical and analytical models of star formation

Stellar Multiplicity

Gravitational interaction between a planet and a three-dimensional isothermal gaseous disk is studied. In the present paper we mainly examine the torque on a planet and the resultant radial migration of the planet. A planet excites density waves at Lindblad and corotation resonances and experiences a negative torque by the density waves, which causes a rapid inward migration of the planet during its formation in a protoplanetary disk. We formulate the linear wave excitation in three-dimensional isothermal disks and calculate the torques of Lindblad resonances and corotation resonances. For corotation resonances, a general torque formula is newly derived, which is also applicable to two-dimensional disks. The new formula succeeds in reproducing numerical results on the corotation torques, which do not agree with the previously well-known formula. The net torque of the inner and the outer Lindblad resonances (i.e., the differential Lindblad torque) is caused by asymmetry such as the radial pressure gradient and the scale height variation. In three-dimensional disks, the differential Lindblad torques are generally smaller than those in two-dimensional disks. Especially, the effect of a pressure gradient becomes weak. The scale height variation, which is a purely three-dimensional effect, makes the differential Lindblad torque decrease. As a result, the migration time of a planet is obtained as of the order of 106 yr for an Earth-size planet at 5 AU for a typical disk model, which is longer than the result of two-dimensional calculation by the factor of 2 or 3. The reflected waves from disk edges, which are neglected in the torque calculation, can further weaken the disk-planet interaction.

Three-Dimensional Interaction between a Planet and an Isothermal Gaseous Disk. I. Corotation and Lindblad Torques and Planet Migration

We present stellar parameters and metallicities, obtained from a detailed spectroscopic analysis, for a large sample of 98 stars known to be orbited by planetary mass companions (almost all known targets), as well as for a volume-limited sample of 41 stars not known to host any planet. For most of the stars the stellar parameters are revised versions of the ones presented in our previous work. However, we also present parameters for 18 stars with planets not previously published, and a compilation of stellar parameters for the remaining 4 planet-hosts for which we could not obtain a spectrum. A comparison of our stellar parameters with values of T eff , log g, and [Fe/H] available in the literature shows a remarkable agreement. In particular, our spectroscopic log g values are now very close to trigonometric log g estimates based on Hipparcos parallaxes. The derived [Fe/H] values are then used to confirm the previously known result that planets are more prevalent around metal-rich stars. Furthermore, we confirm that the frequency of planets is a strongly rising function of the stellar metallicity, at least for stars with [Fe/H] > 0. While only about 3% of the solar metallicity stars in the CORALIE planet search sample were found to be orbited by a planet, this number increases to more than 25% for stars with [Fe/H] above +0.3. Curiously, our results also suggest that these percentages might remain relatively constant for values of [Fe/H] lower than about solar, increasing then linearly with the mass fraction of heavy elements. These results are discussed in the context of the theories of planetary formation.

/pdf/spectroscopic-fe-h-for-98-extra-solar-planet-host-stars-wcdhm4n4nu.pdf

M. Podolak

Papers

Formation of the Giant Planets by Concurrant Accretion of Solids and Gas