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

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

Why Is Supercritical Disk Accretion Feasible?

Although the occurrence of steady supercritical disk accretion onto a black hole has been speculated about since the 1970s, it has not been accurately verified so far. For the first time, we previously demonstrated it through two-dimensional, long-term radiation-hydrodynamic simulations. To clarify why this accretion is possible, we quantitatively investigate the dynamics of a simulated supercritical accretion flow with a mass accretion rate of ~10^2 L_E/c^2 (with L_E and c being, respectively, the Eddington luminosity and the speed of light). We confirm two important mechanisms underlying supercritical disk accretion flow, as previously claimed, one of which is the radiation anisotropy arising from the anisotropic density distribution of very optically thick material. We qualitatively show that despite a very large radiation energy density, E_0>10^2L_E/(4 pi r^2 c) (with r being the distance from the black hole), the radiative flux F_0 cE_0/tau could be small due to a large optical depth, typically tau 10^3, in the disk. Another mechanism is photon trapping, quantified by vE_0, where v is the flow velocity. With a large |v| and E_0, this term significantly reduces the radiative flux and even makes it negative (inward) at r<70r_S, where r_S is the Schwarzschild radius. Due to the combination of these effects, the radiative force in the direction along the disk plane is largely attenuated so that the gravitational force barely exceeds the sum of the radiative force and the centrifugal force. As a result, matter can slowly fall onto the central black hole mainly along the disk plane with velocity much less than the free-fall velocity, even though the disk luminosity exceeds the Eddington luminosity. Along the disk rotation axis, in contrast, the strong radiative force drives strong gas outflows.Comment: 8 pages, 7 figures, accepted for publication in Ap

Potential Vorticity Evolution of a Protoplanetary Disk with An Embedded Protoplanet

We present two-dimensional inviscid hydrodynamic simulations of a protoplanetary disk with an embedded planet, emphasizing the evolution of potential vorticity (the ratio of vorticity to density) and its dependence on numerical resolutions. By analyzing the structure of spiral shocks made by the planet, we show that progressive changes of the potential vorticity caused by spiral shocks ultimately lead to the excitation of a secondary instability. We also demonstrate that very high numerical resolution is required to both follow the potential vorticity changes and identify the location where the secondary instability is first excited. Low-resolution results are shown to give the wrong location. We establish the robustness of a secondary instability and its impact on the torque onto the planet. After the saturation of the instability, the disk shows large-scale non-axisymmetry, causing the torque on the planet to oscillate with large amplitude. The impact of the oscillating torque on the protoplanet's migration remains to be investigated.Comment: 17 pages total with 9 figures (Fig.4,5,9 are in .jpg), accepted to Ap

Motion and gravitational radiation of a binary system consisting of an oscillating and rotating coplanar dusty disk and a point-like object

A binary system composed of an oscillating and rotating coplanar dusty disk and a point mass is considered. The conservative dynamics is treated on the Newtonian level. The effects of gravitational radiation reaction and wave emission are studied to leading quadrupole order. The related waveforms are given. The dynamical evolution of the system is determined semi-analytically exploiting the Hamiltonian equations of motion which comprise the effects both of the Newtonian tidal interaction and the radiation reaction on the motion of the binary system in elliptic orbits. Tidal resonance effects between orbital and oscillatory motions are considered in the presence of radiation damping.Comment: 26 pages, 8 figure

Dynamics of Planetesimals due to Gas Drag from an Eccentric Precessing Disk

We analyze the dynamics of individual kilometer-size planetesimals in circumstellar orbits of a tight binary system. We include both the gravitational perturbations of the secondary star and a non-linear gas drag stemming from an eccentric gas disk with a finite precession rate. We consider several precession rates and eccentricities for the gas, and compare the results with a static disk in circular orbit. The disk precession introduces three main differences with respect to the classical static case: (i) The equilibrium secular solutions generated by the gas drag are no longer fixed points in the averaged system, but limit cycles with frequency equal to the precession rate of the gas. The amplitude of the cycle is inversely dependent on the body size, reaching negligible values for $\sim 50$ km size planetesimals. (ii) The maximum final eccentricity attainable by small bodies is restricted to the interval between the gas eccentricity and the forced eccentricity, and apsidal alignment is no longer guaranteed for planetesimals strongly coupled with the gas. (iii) The characteristic timescales of orbital decay and secular evolution decrease significantly with increasing precession rates, with values up to two orders of magnitude smaller than for static disks. Finally, we apply this analysis to the $\gamma$-Cephei system and estimate impact velocities for different size bodies and values of the gas eccentricity. For high disk eccentricities, we find that the disk precession decreases the velocity dispersion between different size planetesimals, thus contributing to accretional collisions in the outer parts of the disk. The opposite occurs for almost circular gas disks, where precession generates an increase in the relative velocities.Comment: 11 pages, 9 figures. Accepted in MNRA

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

Model for nucleation in GaAs homoepitaxy derived from first principles

The initial steps of MBE growth of GaAs on beta 2-reconstructed GaAs(001) are investigated by performing total energy and electronic structure calculations using density functional theory and a repeated slab model of the surface. We study the interaction and clustering of adsorbed Ga atoms and the adsorption of As_2 molecules onto Ga atom clusters adsorbed on the surface. The stable nuclei consist of bound pairs of Ga adatoms, which originate either from dimerization or from an indirect interaction mediated through the substrate reconstruction. As_2 adsorption is found to be strongly exothermic on sites with a square array of four Ga dangling bonds. Comparing two scenarios where the first As_2 gets incorporated in the incomplete surface layer, or alternatively in a new added layer, we find the first scenario to be preferable. In summary, the calculations suggest that nucleation of a new atomic layer is most likely on top of those surface regions where a partial filling of trenches in the surface has occurred before.Comment: 8 pages, 14 figures, Submitted to Phys. Rev. B (December 15, 1998). Other related publications can be found at http://www.fhi-berlin.mpg.de/th/paper.htm

Super-critical Accretion Flows around Black Holes: Two-dimensional, Radiation-pressure-dominated Disks with Photon-trapping

The quasi-steady structure of super-critical accretion flows around a black hole is studied based on the two-dimensional radiation-hydrodynamical (2D-RHD) simulations. The super-critical flow is composed of two parts: the disk region and the outflow regions above and below the disk. Within the disk region the circular motion as well as the patchy density structure are observed, which is caused by Kelvin-Helmholtz instability and probably by convection. The mass-accretion rate decreases inward, roughly in proportion to the radius, and the remaining part of the disk material leaves the disk to form outflow because of strong radiation pressure force. We confirm that photon trapping plays an important role within the disk. Thus, matter can fall onto the black hole at a rate exceeding the Eddington rate. The emission is highly anisotropic and moderately collimated so that the apparent luminosity can exceed the Eddington luminosity by a factor of a few in the face-on view. The mass-accretion rate onto the black hole increases with increase of the absorption opacity (metalicity) of the accreting matter. This implies that the black hole tends to grow up faster in the metal rich regions as in starburst galaxies or star-forming regions.Comment: 16 pages, 12 figures, accepted for publication in ApJ (Volume 628, July 20, 2005 issue