Engineering topics

Arras, Phil

Publications and source records attributed to Arras, Phil.

8 records

DOE OSTI2017

Tidal Dissipation in WASP-12

WASP-12 is a hot Jupiter system with an orbital period of P = 1.1 days, making it one of the shortest-period giant planets known. Recent transit timing observations by Maciejewski et al. and Patra et al. found a decreasing period with P/| P-dot | = 3.2 Myr. This has been interpreted as evidence of either orbital decay due to tidal dissipation or a long-term oscillation of the apparent period due to apsidal precession. Here, we consider the possibility that it is orbital decay. We show that the parameters of the host star are consistent with either a M {sub *} ≃ 1.3 M {sub ⊙} main sequence star or a M {sub *} ≃ 1.2 M {sub ⊙} subgiant. We find that if the star is on the main sequence, the tidal dissipation is too inefficient to explain the observed P-dot . However, if it is a subgiant, the tidal dissipation is significantly enhanced due to nonlinear wave-breaking of the dynamical tide near the star’s center. The subgiant models have a tidal quality factor Q′{sub ∗}≃2×10{sup 5} and an orbital decay rate that agrees well with the observed P-dot . It would also explain why the planet survived for ≃3 Gyr while the star was on the main sequence and yet is now inspiraling on a 3 Myr timescale. Although this suggests that we are witnessing the last ∼0.1% of the planet’s life, the probability of such a detection is a few percent given the observed sample of ≃30 hot Jupiters in P < 3-day orbits around M {sub *} > 1.2 M {sub ⊙} hosts.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Using Long-term Millisecond Pulsar Timing to Obtain Physical Characteristics of the Bulge Globular Cluster Terzan 5

Over the past decade, the discovery of three unique stellar populations and a large number of confirmed pulsars within the globular cluster Terzan 5 has raised questions over its classification. Using the long-term radio pulsar timing of 36 ms pulsars in the cluster core, we provide new measurements of key physical properties of the system. As Terzan 5 is located within the galactic bulge, stellar crowding and reddening make optical and near-infrared observations difficult. Pulsar accelerations, however, allow us to study the intrinsic characteristics of the cluster independent of reddening and stellar crowding and probe the mass density profile without needing to quantify the mass-to-light ratio. Relating the spin and orbital periods of each pulsar to the acceleration predicted by a King model, we find a core density of 1.58{sub −0.13}{sup +0.13} × 10{sup 6} M{sub ⊙} pc{sup −3}, a core radius of 0.16{sub −0.01}{sup +0.01} pc, a pulsar density profile of n∝r{sup −3.14{sub −}{sub 0}{sub .}{sub 5}{sub 3}{sup +{sup 0{sup .{sup 5{sup 2}}}}}}, and a total mass of M{sub T}(R{sub ⊥}<1.0 pc) ≃ 3.0 × 10{sup 5} M{sub ⊙}, assuming a cluster distance of 5.9 kpc. Using this information, we argue against Terzan 5 being a disrupted dwarf galaxy and discuss the possibility of it being a fragment of the Milky Way’s proto-bulge. We also discuss whether low-mass pulsars were formed via electron-capture supernovae or exist in a core full of heavy white dwarfs and hard binaries. Finally, we provide an upper limit for the mass of a possible black hole at the core of the cluster of M{sub BH}≃3×10{sup 4} M{sub ⊙}.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

AXISYMMETRIC SIMULATIONS OF HOT JUPITER–STELLAR WIND HYDRODYNAMIC INTERACTION

Gas giant exoplanets orbiting at close distances to the parent star are subjected to large radiation and stellar wind fluxes. In this paper, hydrodynamic simulations of the planetary upper atmosphere and its interaction with the stellar wind are carried out to understand the possible flow regimes and how they affect the Lyα transmission spectrum. Following Tremblin and Chiang, charge exchange reactions are included to explore the role of energetic atoms as compared to thermal particles. In order to understand the role of the tail as compared to the leading edge of the planetary gas, the simulations were carried out under axisymmetry, and photoionization and stellar wind electron impact ionization reactions were included to limit the extent of the neutrals away from the planet. By varying the planetary gas temperature, two regimes are found. At high temperature, a supersonic planetary wind is found, which is turned around by the stellar wind and forms a tail behind the planet. At lower temperatures, the planetary wind is shut off when the stellar wind penetrates inside where the sonic point would have been. In this regime mass is lost by viscous interaction at the boundary between planetary and stellar wind gases. Absorption by cold hydrogen atoms is large near the planetary surface, and decreases away from the planet as expected. The hot hydrogen absorption is in an annulus and typically dominated by the tail, at large impact parameter, rather than by the thin leading edge of the mixing layer near the substellar point.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Magnetohydrodynamic simulations of hot jupiter upper atmospheres

Two-dimensional simulations of hot Jupiter upper atmospheres including the planet's magnetic field are presented. The goal is to explore magnetic effects on the layer of the atmosphere that is ionized and heated by stellar EUV radiation, and the imprint of these effects on the Lyα transmission spectrum. The simulations are axisymmetric, isothermal, and include both rotation and azimuth-averaged stellar tides. Mass density is converted to atomic hydrogen density through the assumption of ionization equilibrium. The three-zone structure—polar dead zone (DZ), mid-latitude wind zone (WZ), and equatorial DZ—found in previous analytic calculations is confirmed. For a magnetic field comparable to that of Jupiter, the equatorial DZ, which is confined by the magnetic field and corotates with the planet, contributes at least half of the transit signal. For even stronger fields, the gas escaping in the mid-latitude WZ is found to have a smaller contribution to the transit depth than the equatorial DZ. Transmission spectra computed from the simulations are compared to Hubble Space Telescope Space Telescope Imaging Spectrograph and Advanced Camera for Surveys data for HD 209458b and HD 189733b, and the range of model parameters consistent with the data is found. The central result of this paper is that the transit depth increases strongly with magnetic field strength when the hydrogen ionization layer is magnetically dominated, for dipole magnetic field B {sub 0} ≳ 10 G. Hence transit depth is sensitive to magnetic field strength, in addition to standard quantities such as the ratio of thermal to gravitational binding energies. Another effect of the magnetic field is that the planet loses angular momentum orders of magnitude faster than in the non-magnetic case, because the magnetic field greatly increases the lever arm for wind braking of the planet's rotation. Spin-down timescales for magnetized models of HD 209458b that agree with the observed transit depth can be as short as ≅ 30 Myr, which is much shorter than the age of the system.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

MODULES FOR EXPERIMENTS IN STELLAR ASTROPHYSICS (MESA): PLANETS, OSCILLATIONS, ROTATION, AND MASSIVE STARS

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{sub Sun} 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

H{alpha} ABSORPTION IN TRANSITING EXOPLANET ATMOSPHERES

Absorption of stellar H{alpha} by the upper atmosphere of the planet HD 189733b has recently been detected by Jensen et al. Motivated by this observation, we have developed a model for atomic hydrogen in the n = 2 state and compared the resulting H{alpha} line profile to the observations. The model atmosphere is in hydrostatic balance, as well as thermal and photoionization equilibrium. Collisional and radiative transitions are included in the determination of the n = 2 state level population. We find that H{alpha} absorption is dominated by an optical depth {tau} {approx} 1 shell, composed of hydrogen in the metastable 2s state that is located below the hydrogen ionization layer. The number density of the 2s state within the shell is found to vary slowly with radius, while that of the 1s state falls rapidly. Thus while the Ly{alpha} absorption, for a certain wavelength, occurs inside a relatively well defined impact parameter, the contribution to H{alpha} absorption is roughly uniform over the entire atomic hydrogen layer. The model can approximately reproduce the observed Ly{alpha} and H{alpha} integrated transit depths for HD 189733b by using an ionization rate enhanced over that expected for the star by an order of magnitude. For HD 209458b, we are unable to explain the asymmetric H{alpha} line profile observed by Jensen et al., as the model produces a symmetric line profile with transit depth comparable to that of HD 189733b. In an appendix, we study the effect of the stellar Ly{alpha} absorption on the net cooling rate.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

An instability due to the nonlinear coupling of p-modes to g-modes: Implications for coalescing neutron star binaries

A weakly nonlinear fluid wave propagating within a star can be unstable to three-wave interactions. The resonant parametric instability is a well-known form of three-wave interaction in which a primary wave of frequency ω {sub a} excites a pair of secondary waves of frequency ω {sub b} + ω {sub c} ≅ ω {sub a}. Here we consider a nonresonant form of three-wave interaction in which a low-frequency primary wave excites a high-frequency p-mode and a low-frequency g-mode such that ω {sub b} + ω {sub c} >> ω {sub a}. We show that a p-mode can couple so strongly to a g-mode of similar radial wavelength that this type of nonresonant interaction is unstable even if the primary wave amplitude is small. As an application, we analyze the stability of the tide in coalescing neutron star binaries to p-g mode coupling. We find that the equilibrium tide and dynamical tide are both p-g unstable at gravitational wave frequencies f {sub gw} ≳ 20 Hz and drive short wavelength p-g mode pairs to significant energies on very short timescales (much less than the orbital decay time due to gravitational radiation). Resonant parametric coupling to the tide is, by contrast, either stable or drives modes at a much smaller rate. We do not solve for the saturation of the p-g instability and therefore we cannot say precisely how it influences the evolution of neutron star binaries. However, we show that if even a single daughter mode saturates near its wave breaking amplitude, the p-g instability of the equilibrium tide will (1) induce significant orbital phase errors (Δφ ≳ 1 radian) that accumulate primarily at low frequencies (f {sub gw} ≲ 50 Hz) and (2) heat the neutron star core to a temperature of T ∼ 10{sup 10} K. Since there are at least ∼100 unstable p-g daughter pairs, Δφ and T are potentially much larger than these values. Tides might therefore significantly influence the gravitational wave signal and electromagnetic emission from coalescing neutron star binaries at much larger orbital separations than previously thought.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

NONLINEAR TIDES IN CLOSE BINARY SYSTEMS

We study the excitation and damping of tides in close binary systems, accounting for the leading-order nonlinear corrections to linear tidal theory. These nonlinear corrections include two distinct physical effects: three-mode nonlinear interactions, i.e., the redistribution of energy among stellar modes of oscillation, and nonlinear excitation of stellar normal modes by the time-varying gravitational potential of the companion. This paper, the first in a series, presents the formalism for studying nonlinear tides and studies the nonlinear stability of the linear tidal flow. Although the formalism we present is applicable to binaries containing stars, planets, and/or compact objects, we focus on non-rotating solar-type stars with stellar or planetary companions. Our primary results include the following: (1) The linear tidal solution almost universally used in studies of binary evolution is unstable over much of the parameter space in which it is employed. More specifically, resonantly excited internal gravity waves in solar-type stars are nonlinearly unstable to parametric resonance for companion masses M' {approx}> 10-100 M{sub Circled-Plus} at orbital periods P Almost-Equal-To 1-10 days. The nearly static 'equilibrium' tidal distortion is, however, stable to parametric resonance except for solar binaries with P {approx}< 2-5 days. (2) For companion masses larger than a few Jupiter masses, the dynamical tide causes short length scale waves to grow so rapidly that they must be treated as traveling waves, rather than standing waves. (3) We show that the global three-wave treatment of parametric instability typically used in the astrophysics literature does not yield the fastest-growing daughter modes or instability threshold in many cases. We find a form of parametric instability in which a single parent wave excites a very large number of daughter waves (N Almost-Equal-To 10{sup 3}[P/10 days] for a solar-type star) and drives them as a single coherent unit with growth rates that are a factor of Almost-Equal-To N faster than the standard three-wave parametric instability. These are local instabilities viewed through the lens of global analysis; the coherent global growth rate follows local rates in the regions where the shear is strongest. In solar-type stars, the dynamical tide is unstable to this collective version of the parametric instability for even sub-Jupiter companion masses with P {approx}< a month. (4) Independent of the parametric instability, the dynamical and equilibrium tides excite a wide range of stellar p-modes and g-modes by nonlinear inhomogeneous forcing; this coupling appears particularly efficient at draining energy out of the dynamical tide and may be more important than either wave breaking or parametric resonance at determining the nonlinear dissipation of the dynamical tide.

79 ASTRONOMY AND ASTROPHYSICS↗