Atmospheric Dynamics of Short-period Extra Solar Gas Giant Planets I: Dependence of Night-Side Temperature on Opacity

More than two dozen short-period Jupiter-mass gas giant planets have been discovered around nearby solar-type stars in recent years, several of which undergo transits, making them ideal for the detection and characterization of their atmospheres. Here we adopt a three-dimensional radiative hydrodynamical numerical scheme to simulate atmospheric circulation on close-in gas giant planets. In contrast to the conventional GCM and shallow water algorithms, this method does not assume quasi hydrostatic equilibrium and it approximates radiation transfer from optically thin to thick regions with flux-limited diffusion. In the first paper of this series, we consider synchronously-spinning gas giants. We show that a full three-dimensional treatment, coupled with rotationally modified flows and an accurate treatment of radiation, yields a clear temperature transition at the terminator. Based on a series of numerical simulations with varying opacities, we show that the night-side temperature is a strong indicator of the opacity of the planetary atmosphere. Planetary atmospheres that maintain large, interstellar opacities will exhibit large day-night temperature differences, while planets with reduced atmospheric opacities due to extensive grain growth and sedimentation will exhibit much more uniform temperatures throughout their photosphere's. In addition to numerical results, we present a four-zone analytic approximation to explain this dependence.Comment: 35 Pages, 13 Figure

Tidal Barrier and the Asymptotic Mass of Proto Gas-Giant Planets

Extrasolar planets found with radial velocity surveys have masses ranging from several Earth to several Jupiter masses. While mass accretion onto protoplanetary cores in weak-line T-Tauri disks may eventually be quenched by a global depletion of gas, such a mechanism is unlikely to have stalled the growth of some known planetary systems which contain relatively low-mass and close-in planets along with more massive and longer period companions. Here, we suggest a potential solution for this conundrum. In general, supersonic infall of surrounding gas onto a protoplanet is only possible interior to both of its Bondi and Roche radii. At a critical mass, a protoplanet's Bondi and Roche radii are equal to the disk thickness. Above this mass, the protoplanets' tidal perturbation induces the formation of a gap. Although the disk gas may continue to diffuse into the gap, the azimuthal flux across the protoplanets' Roche lobe is quenched. Using two different schemes, we present the results of numerical simulations and analysis to show that the accretion rate increases rapidly with the ratio of the protoplanet's Roche to Bondi radii or equivalently to the disk thickness. In regions with low geometric aspect ratios, gas accretion is quenched with relatively low protoplanetary masses. This effect is important for determining the gas-giant planets' mass function, the distribution of their masses within multiple planet systems around solar type stars, and for suppressing the emergence of gas-giants around low mass stars

On disc driven inward migration of resonantly coupled planets with application to the system around GJ876

We consider two protoplanets gravitationally interacting with each other and a protoplanetary disc. The two planets orbit interior to a tidally maintained disc cavity while the disc interaction indices inward migration. When the migration is slow enough, the more rapidly migrating outer protoplanet approaches and becomes locked in a 2:1 commensurability with the inner one. This is maintained in subsequent evolution. We study this evolution using a simple anaytic model, full hydrodynamic 2D simulations of the disc planet system and longer time N body integrations incorporating simple prescriptions for the effect of the disc on the planet orbits. The eccentricity of the protoplanets are found to be determined by the migration rate induced in the outer planet orbit by the external disc. We apply our results to the recently discovered resonant planets around GJ876. Simulation shows that a disc with parameters expected for protoplanetary discs causes trapping in the 2:1 commensurability when the planets orbit in an inner cavity and that eccentricities in the observed range may be obtained.Comment: 8 pages, 5 figures, submitted to A&A on 30/03/200

3D-MHD simulations of an accretion disk with star-disk boundary layer

We present global 3D MHD simulations of geometrically thin but unstratified accretion disks in which a near Keplerian disk rotates between two bounding regions with initial rotation profiles that are stable to the MRI. The inner region models the boundary layer between the disk and an assumed more slowly rotating central, non magnetic star. We investigate the dynamical evolution of this system in response to initial vertical and toroidal fields imposed in a variety of domains contained within the near Keplerian disk. Cases with both non zero and zero net magnetic flux are considered and sustained dynamo activity found in runs for up to fifty orbital periods at the outer boundary of the near Keplerian disk. Simulations starting from fields with small radial scale and with zero net flux lead to the lowest levels of turbulence and smoothest variation of disk mean state variables. For our computational set up, average values of the Shakura & Sunyaev (1973) $\alpha$ parameter in the Keplerian disk are typically $0.004\pm 0.002.$ Magnetic field eventually always diffuses into the boundary layer resulting in the build up of toroidal field inward angular momentum transport and the accretion of disk material. The mean radial velocity, while exhibiting large temporal fluctuations is always subsonic. Simulations starting with net toroidal flux may yield an average $\alpha \sim 0.04.$ While being characterized by one order of magnitude larger average $\alpha$, simulations starting from vertical fields with large radial scale and net flux may lead to the formation of persistent non-homogeneous, non-axisymmetric magnetically dominated regions of very low density.Comment: Accepted for publication in Ap

On the observability of bow shocks of Galactic runaway OB stars

Massive stars that have been ejected from their parent cluster and supersonically sailing away through the interstellar medium (ISM) are classified as exiled. They generate circumstellar bow shock nebulae that can be observed. We present two-dimensional, axisymmetric hydrodynamical simulations of a representative sample of stellar wind bow shocks from Galactic OB stars in an ambient medium of densities ranging from n_ISM=0.01 up to 10.0/cm3. Independently of their location in the Galaxy, we confirm that the infrared is the most appropriated waveband to search for bow shocks from massive stars. Their spectral energy distribution is the convenient tool to analyze them since their emission does not depend on the temporary effects which could affect unstable, thin-shelled bow shocks. Our numerical models of Galactic bow shocks generated by high-mass (~40 Mo) runaway stars yield H$\alpha$ fluxes which could be observed by facilities such as the SuperCOSMOS H-Alpha Survey. The brightest bow shock nebulae are produced in the denser regions of the ISM. We predict that bow shocks in the field observed at Ha by means of Rayleigh-sensitive facilities are formed around stars of initial mass larger than about 20 Mo. Our models of bow shocks from OB stars have the emission maximum in the wavelength range 3 <= lambda <= 50 micrometer which can be up to several orders of magnitude brighter than the runaway stars themselves, particularly for stars of initial mass larger than 20 Mo.Comment: 13 pages, 12 figures. Accepted to MNRAS (2016

The Vertical Structure of Planet-induced Gaps in Proto-Planetary Discs

Giant planets embedded in circumstellar discs are expected to open gaps in these discs. We examine the vertical structure of the gap edges. We find that the planet excites spiral arms with significant (Mach number of a half) vertical motion of the gas, and discuss the implications of these motions. In particular, the spiral arms will induce strong vertical stirring of the dust, making the edge appeared `puffed up' relative to the bulk of the disc. Infra-red observations (sensitive to dust) would be dominated by the light from the thick inner edge of the disc. Sub-millimetre observations (sensitive to gas velocities) would appear to be hot in `turbulent' motions (actually the ordered motion caused by the passage of the spiral arms), but cold in chemistry. Resolved sub-millimetre maps of circumstellar discs might even be able to detect the spiral arms directly.Comment: Revision adds new data, and corrects physical intepretatio

Dynamical evolution and leading order gravitational wave emission of Riemann-S binaries

An approximate strategy for studying the evolution of binary systems of extended objects is introduced. The stars are assumed to be polytropic ellipsoids. The surfaces of constant density maintain their ellipsoidal shape during the time evolution. The equations of hydrodynamics then reduce to a system of ordinary differential equations for the internal velocities, the principal axes of the stars and the orbital parameters. The equations of motion are given within Lagrangian and Hamiltonian formalism. The special case when both stars are axially symmetric fluid configurations is considered. Leading order gravitational radiation reaction is incorporated, where the quasi-static approximation is applied to the internal degrees of freedom of the stars. The influence of the stellar parameters, in particular the influence of the polytropic index $n$, on the leading order gravitational waveforms is studied.Comment: 31 pages, 7 figures, typos correcte

Orbital migration and the frequency of giant planet formation

We present a statistical study of the post-formation migration of giant planets in a range of initial disk conditions. For given initial conditions we model the evolution of giant planet orbits under the influence of disk, stellar, and mass loss torques. We determine the mass and semi-major axis distribution of surviving planets after disk dissipation, for various disk masses, lifetimes, viscosities, and initial planet masses. The majority of planets migrate too fast and are destroyed via mass transfer onto the central star. Most surviving planets have relatively large orbital semi-major axes of several AU or larger. We conclude that the extrasolar planets observed to date, particularly those with small semi-major axes, represent only a small fraction (~25% to 33%) of a larger cohort of giant planets around solar-type stars, and many undetected giant planets must exist at large (>1-2 AU) distances from their parent stars. As sensitivity and completion of the observed sample increases with time, this distant majority population of giant planets should be revealed. We find that the current distribution of extrasolar giant planet masses implies that high mass (more than 1-2 Jupiter masses) giant planet formation must be relatively rare. Finally, our simulations imply that the efficiency of giant planet formation must be high: at least 10% and perhaps as many as 80% of solar-type stars possess giant planets during their pre-main sequence phase. These predictions, including those for pre-main sequence stars, are testable with the next generation of ground- and space-based planet detection techniquesComment: 25 pages, 5 figures. Double-space, single-column format to show long equations. Accepted for publication in A&

Black hole formation via hypercritical accretion during common envelope evolution

Neutron stars inspiralling into a stellar envelope can accrete at rates vastly exceeding the Eddington limit if the flow develops pressures high enough to allow neutrinos to radiate the released gravitational energy. It has been suggested that this hypercritical mode of accretion leads inevitably to the formation of stellar mass black holes during common envelope evolution. We study the hydrodynamics of this flow at large radii (R >> R_ns), and show that for low Mach number flows, in two dimensions, modest density gradients in the stellar envelope suffice to produce a hot, advection dominated accretion disk around the accreting object. The formation of outflows from such a disk is highly probable, and we discuss the impact of the resultant mass loss and feedback of energy into the envelope for the survival of the neutron star. Unless outflows are weaker than those inferred for well observed accreting systems, we argue that in most cases insufficient accretion occurs to force collapse to a black hole before the envelope has been ejected. This conclusions is of interest for black hole formation in general, for some models of gamma ray bursts, and for predictions of the event rate in future LIGO observations.Comment: ApJ, submitte