No Arabic abstract
To shed more light on the nature of the observed Ly{alpha} absorption during transits of HD 209458b and to quantify the major mechanisms responsible for the production of fast hydrogen atoms (the so called energetic neutral atoms, ENAs) around the planet, 2D hydrodynamic multifluid modeling of the expanding planetary upper atmosphere, which is driven by stellar XUV, and its interaction with the stellar wind has been performed. The model selfconsistently describes the escaping planetary wind, taking into account the generation of ENAs due to particle acceleration by the radiation pressure and by the charge exchange between the stellar wind protons and planetary atoms. The calculations in a wide range of stellar wind parameters and XUV flux values showed that under typical Sun-like star conditions, the amount of generated ENAs is too small, and the observed absorption at the level of 6-8 percent can be attributed only to the non-resonant natural line broadening. For lower XUV fluxes, e.g., during the activity minima, the number of planetary atoms that survive photoionization and give rise to ENAs increases, resulting in up to 10-15 percent absorption at the blue wing of the Lya line, caused by resonant thermal line broadening. A similar asymmetric absorption can be seen under the conditions realized during coronal mass ejections, when sufficiently high stellar wind pressure confines the escaping planetary material within a kind of bowshock around the planet. It was found that the radiation pressure in all considered cases has a negligible contribution to the production of ENAs and the corresponding absorption.
The absorption of stellar radiation observed by the HD209458b in resonant lines of OI and CII has not yet been satisfactorily modeled. In our previous 2D simulations we have shown that the hydrogen-dominated upper atmosphere of HD209458b, heated by XUV radiation, expands supersonically beyond the Roche lobe and drags the heavier species along with it. Assuming solar abundances, OI and CII particles accelerated by tidal forces to velocities up to 50 km/s should produce the absorption due to Doppler resonance mechanism at the level of 6-10%, consistent with the observations. Since the 2D geometry does not take into account the Coriolis force in the planet reference frame, the question remained to which extent the spiraling of the escaping planetary material and its actually achieved velocity may influence the conclusions made on the basis of 2D modeling. In the present paper we apply for the first time in the study of HD209458b a global 3D hydrodynamic multi-fluid model that self-consistently describes the formation and expansion of the escaping planetary wind, affected by the tidal and Coriolis forces, as well as by the surrounding stellar wind. The modeling results confirm our previous findings that the velocity and density of the planetary flow are sufficiently high to produce the absorption in HI, OI, and CII resonant lines at the level close to the in-transit observed values. The novel finding is that the matching of the absorption measured in MgII and SiIII lines requires at least 10 times lower abundances of these elements than the Solar system values.
Warm Neptune GJ3470b has been recently observed in 23S-23P transition of metastable helium, yielding absorption of about 1% in Doppler velocity range of [-40; 10] km/s. Along with previous detection of absorption in Ly{alpha} with depth of 20-40% in the blue and red wings of the line, it offers a complex target for simulation and testing of the current models. Obtained results suggest that absorption in both these lines comes from interaction of expanding upper planetary atmosphere with stellar plasma wind, allowing to constrain the stellar plasma parameters and the helium abundance in planet atmosphere.
We present the first three-dimensional magnetohydrodynamic (MHD) simulations of the atmosphere of HD 209458b which self-consistently include reduction of winds due to the Lorentz force and Ohmic heating. We find overall wind structures similar to that seen in previous models of hot Jupiter atmospheres, with strong equatorial jets and meridional flows poleward near the day side and equatorward near the night side. Inclusion of magnetic fields slows those winds and leads to Ohmic dissipation. We find wind slowing ranging from 10%-40% for reasonable field strengths. We find Ohmic dissipation rates ~10^17 W at 100 bar, orders of magnitude too small to explain the inflated radius of this planet. Faster wind speeds, not achievable in these anelastic calculations, may be able to increase this value somewhat, but likely will not be able to close the gap necessary to explain the inflated radius. We demonstrate that the discrepancy between the simulations presented here and previous models is due to inadequate treatment of magnetic field geometry and evolution. Induced poloidal fields become much larger than those imposed, highlighting the need for a self-consistent MHD treatment of these hot atmospheres.
We perform a systematic combined photometric and kinematic analysis of a sample of globular clusters under different relaxation conditions, based on their core relaxation time (as listed in available catalogs), by means of two well-known families of spherical stellar dynamical models. Systems characterized by shorter relaxation time scales are expected to be better described by isotropic King models, while less relaxed systems might be interpreted by means of non-truncated, radially-biased anisotropic f^( u) models, originally designed to represent stellar systems produced by a violent relaxation formation process and applied here for the first time to the study of globular clusters. The comparison between dynamical models and observations is performed by fitting simultaneously surface brightness and velocity dispersion profiles. For each globular cluster, the best-fit model in each family is identified, along with a full error analysis on the relevant parameters. Detailed structural properties and mass-to-light ratios are also explicitly derived. We find that King models usually offer a good representation of the observed photometric profiles, but often lead to less satisfactory fits to the kinematic profiles, independently of the relaxation condition of the systems. For some less relaxed clusters, f^( u) models provide a good description of both observed profiles. Some derived structural characteristics, such as the total mass or the half-mass radius, turn out to be significantly model-dependent. The analysis confirms that, to answer some important dynamical questions that bear on the formation and evolution of globular clusters, it would be highly desirable to acquire larger numbers of accurate kinematic data-points, well distributed over the cluster field.
There is evidence that the transiting planet HD 209458b has a large exosphere of neutral hydrogen, based on a 15% decrement in Lyman-alpha flux that was observed by Vidal-Madjar et al. during transits. Here we report upper limits on H-alpha absorption by the exosphere. The results are based on optical spectra of the parent star obtained with the Subaru High Dispersion Spectrograph. Comparison of the spectra taken inside and outside of transit reveals no exospheric H-alpha signal greater than 0.1% within a 5.1A band (chosen to have the same Delta_lambda/lambda as the 15% Ly-alpha absorption). The corresponding limit on the column density of n=2 neutral hydrogen is N_2 <~ 10^9 cm^{-2}. This limit constrains proposed models involving a hot (~10^4 K) and hydrodynamically escaping exosphere.