ABSTRACT The role of radiation pressure in shaping exoplanet photoevaporation remains a topic of contention. Radiation pressure from the exoplanet’s host star has been proposed as a mechanism to drive the escaping atmosphere into a ‘cometary’ tail and explain the high velocities observed in systems where mass-loss is occurring. In this paper, we present results from high-resolution 3D hydrodynamic simulations of a planet similar to HD 209458b. We self-consistently launch a wind flowing outwards from the planet by calculating the ionization and heating resulting from incident high-energy radiation, and account for radiation pressure. We first present a simplified calculation, setting a limit on the Lyman-α flux required to drive the photoevaporated planetary material to larger radii and line-of-sight velocities. We then present the results of our simulations, which confirm the limits determined by our analytic calculation. We thus demonstrate that, within the limits of our hydrodynamic simulation and for the Lyman-α fluxes expected for HD 209458, radiation pressure is unlikely to significantly affect photoevaporative winds or to explain the high velocities at which wind material is observed, though further possibilities remain to be investigated.
more »
« less
Effects of charge exchange on the evaporative wind of HD 209458b
The role of charge exchange in shaping exoplanet photoevaporation remains a topic of contention. Exchange of electrons between stellar wind protons from the exoplanet’s host star and neutral hydrogen from the planet’s wind has been proposed as a mechanism to create ‘energetic neutral atoms’ (ENAs), which could explain the high absorption line velocities observed in systems where mass-loss is occurring. In this paper, we present results from three-dimensional hydrodynamic simulations of the mass-loss of a planet similar to HD 209458b. We self-consistently launch a planetary wind by calculating the ionization and heating resulting from incident high-energy radiation, inject a stellar wind into the simulation, and allow electron exchange between the stellar and planetary winds. We predict the potential production of ENAs by the wind–wind interaction analytically, and then present the results of our simulations, which confirm the analytic limits. Within the limits of our hydrodynamic simulation, we find that charge exchange with the stellar wind properties examined here is unable to explain the absorption observed at high Doppler velocities.
more » « less- Award ID(s):
- 1813298
- NSF-PAR ID:
- 10373640
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 517
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 1724-1736
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Interstellar neutral atoms propagating into the heliosphere experience charge exchange with the supersonic solar wind (SW) plasma, generating ions that are picked up by the SW. These pickup ions (PUIs) constitute ∼25% of the proton number density by the time they reach the heliospheric termination shock (HTS). Preferential acceleration of PUIs at the HTS leads to a suprathermal, kappa-like PUI distribution in the heliosheath, which may be further heated in the heliosheath by traveling shocks or pressure waves. In this study, we utilize a dynamic, 3D magnetohydrodynamic model of the heliosphere to show that dynamic heating of PUIs at the HTS and in the inner heliosheath (IHS), as well as a background source of energetic neutral atoms (ENAs) from outside the heliopause, can explain the heliospheric ENA signal observed by the Interstellar Boundary Explorer (IBEX) in the Voyager 2 direction. We show that the PUI heating process at the HTS is characterized by a polytropic index larger than 5/3, likely ranging between
γ ∼ 2.3 and 2.7, depending on the time in solar cycle 24 and SW conditions. The ENA fluxes at energies >1.5 keV show large-scale behavior in time with the solar cycle and SW dynamic pressure, whereas ENAs < 1.5 keV primarily exhibit random-like fluctuations associated with SW transients affecting the IHS. We find that ≲20% of the ENAs observed at ∼0.5–6 keV come from other sources, likely from outside the heliopause as secondary ENAs. This study offers the first model replication of the intensity and evolution of IBEX-Hi ENA observations from the outer heliosphere. -
Interacting binaries are of general interest as laboratories for investigating the physics of accretion, which gives rise to the bulk of high-energy radiation in the Galaxy. They allow us to probe stellar evolution processes that cannot be studied in single stars. Understanding the orbital evolution of binaries is essential in order to model the formation of compact binaries. Here we focus our attention on studying orbital evolution driven by angular momentum loss through stellar winds in massive binaries. We run a suite of hydrodynamical simulations of binary stars hosting one mass losing star with varying wind velocity, mass ratio, wind velocity profile and adiabatic index, and compare our results to analytic estimates for drag and angular momentum loss. We find that, at leading order, orbital evolution is determined by the wind velocity and the binary mass ratio. Small ratios of wind to orbital velocities and large accreting companion masses result in high angular momentum loss and a shrinking of the orbit. For wider binaries and binaries hosting lighter mass-capturing companions, the wind mass-loss becomes more symmetric, which results in a widening of the orbit. We present a simple analytic formula that can accurately account for angular momentum losses and changes in the orbit, which depends on the wind velocity and mass ratio. As an example of our formalism, we compare the effects of tides and winds in driving the orbital evolution of high mass X-ray binaries, focusing on Vela X-1 and Cygnus X-1 as examples.more » « less
-
more » « less
ABSTRACT Although galactic winds play a critical role in regulating galaxy formation, hydrodynamic cosmological simulations do not resolve the scales that govern the interaction between winds and the ambient circumgalactic medium (CGM). We implement the Physically Evolved Wind (PhEW) model of Huang et al. in the gizmo hydrodynamics code and perform test cosmological simulations with different choices of model parameters and numerical resolution. PhEW adopts an explicit subgrid model that treats each wind particle as a collection of clouds that exchange mass and metals with their surroundings and evaporate by conduction and hydrodynamic instabilities as calibrated on much higher resolution cloud scale simulations. In contrast to a conventional wind algorithm, we find that PhEW results are robust to numerical resolution and implementation details because the small scale interactions are defined by the model itself. Compared to our previous wind simulations with the same resolution, our PhEW simulations are in better agreement with low-redshift galactic stellar mass functions at M* < 1011M⊙ because PhEW particles shed mass to the CGM before escaping low mass haloes. PhEW radically alters the CGM metal distribution because PhEW particles disperse metals to the ambient medium as their clouds dissipate, producing a CGM metallicity distribution that is skewed but unimodal and is similar between cold and hot gas. While the temperature distributions and radial profiles of gaseous haloes are similar in simulations with PhEW and conventional winds, these changes in metal distribution will affect their predicted UV/X-ray properties in absorption and emission.
-
more » « less
Abstract Atmospheric escape from close-in exoplanets is thought to be crucial in shaping observed planetary populations. Recently, significant progress has been made in observing this process in action through excess absorption in-transit spectra and narrowband light curves. We model the escape of initially homogeneous planetary winds interacting with a stellar wind. The ram pressure balance of the two winds governs this interaction. When the impingement of the stellar wind on the planetary outflow is mild or moderate, the planetary outflow expands nearly spherically through its sonic surface before forming a shocked boundary layer. When the confinement is strong, the planetary outflow is redirected into a cometary tail before it expands to its sonic radius. The resultant transmission spectra at the He 1083 nm line are accurately represented by a 1D spherical wind solution in cases of mild to moderate stellar wind interaction. In cases of strong stellar wind interaction, the degree of absorption is enhanced and the cometary tail leads to an extended egress from transit. The crucial features of the wind–wind interaction are, therefore, encapsulated in the light curve of He 1083 nm equivalent width as a function of time. The possibility of extended He 1083 nm absorption well beyond the optical transit carries important implications for planning out-of-transit observations that serve as a baseline for in-transit data.
