The habitability of the surface of any planet is determined by a complex evolution of its interior, surface, and atmosphere. The electromagnetic and particle radiation of stars drive thermal, chemical and physical alteration of planetary atmospheres, including escape. Many known extrasolar planets experience vastly different stellar environments than those in our Solar system: it is crucial to understand the broad range of processes that lead to atmospheric escape and evolution under a wide range of conditions if we are to assess the habitability of worlds around other stars. One problem encountered between the planetary and the astrophysics communities is a lack of common language for describing escape processes. Each community has customary approximations that may be questioned by the other, such as the hypothesis of H-dominated thermosphere for astrophysicists, or the Sun-like nature of the stars for planetary scientists. Since exoplanets are becoming one of the main targets for the detection of life, a common set of definitions and hypotheses are required. We review the different escape mechanisms proposed for the evolution of planetary and exoplanetary atmospheres. We propose a common definition for the different escape mechanisms, and we show the important parameters to take into account when evaluating the escape at a planet in time. We show that the paradigm of the magnetic field as an atmospheric shield should be changed and that recent work on the history of Xenon in Earth's atmosphere gives an elegant explanation to its enrichment in heavier isotopes: the so-called Xenon paradox.

Observing the dynamics of compact astrophysical objects provides insights into their inner workings and allows to probe physics under extreme conditions. The immediate vicinity of an active supermassive black hole with its event horizon, photon ring, accretion disk, and relativistic jets is a perfect place to study general relativity, magneto-hydrodynamics, and high energy plasma physics. The recent observations of the black hole shadow of M87* with Very Long Baseline Interferometry (VLBI) by the Event Horizon Telescope (EHT) open the possibility to investigate dynamical processes there on timescales of days. In this regime, radio astronomical imaging algorithms are brought to their limits. Compared to regular radio interferometers, VLBI networks have fewer antennas. The resulting sparser sampling of the Fourier sky can only be partly compensated by co-adding observations from different days, as the source changes. Here, we present an imaging algorithm that copes with the data scarcity and the source's temporal evolution, while simultaneously providing uncertainty quantification on all results. Our algorithm views the imaging task as a Bayesian inference problem of a time-varying flux density, exploits the correlation structure between time frames, and reconstructs a whole, $(2+1+1)$-dimensional time-variable and spectral-resolved image at once. We apply the method to the EHT observation of M87* and validate our approach on synthetic data. The obtained first time-resolved reconstruction of M87* indicates varying structures on and outside the emission ring on a time scale of days.

The space-borne gravitational wave interferometer, LISA, is expected to detect signals from numerous binary white dwarfs. At small orbital separation, rapid rotation and large tidal bulges may allow for the stellar internal structure to be probed through such observations. Finite-size effects are encoded in quantities like the moment of inertia ($I$), tidal Love number (Love), and quadrupole moment ($Q$). The universal relations among them (I-Love-Q relations) can be used to reduce the number of parameters in the gravitational-wave templates. We here study I-Love-Q relations for more realistic white dwarf models than used in previous studies. In particular, we extend previous works by including (i) differential rotation and (ii) internal temperature profiles taken from detailed stellar evolution calculations. We use the publicly available stellar evolution code MESA to generate cooling models of both low- and high-mass white dwarfs. We show that differential rotation causes the I-Q relation (and similarly the Love-Q relation) to deviate from that of constant rotation. We also find that the introduction of finite temperatures causes the white dwarf to move along the zero-temperature mass sequence of I-Q values, moving towards values that suggest a lower mass. We further find that after only a few Myrs, high-mass white dwarfs are well-described by the zero-temperature model, suggesting that the relations with zero-temperature may be good enough in most practical cases. Low-mass, He-core white dwarfs with thick hydrogen envelopes may undergo long periods of H burning which sustain the stellar temperature and allow deviations from the I-Love-Q relations for longer times.

We analyze the full $\mathit{Kepler}$ phase curve of KOI-964, a binary system consisting of a hot white dwarf on an eclipsing orbit around an A-type host star. Using all 18 quarters of long-cadence photometry, we carry out a joint light curve fit and obtain improved phase curve amplitudes, occultation depths, orbital parameters, and transit ephemeris over the previous results of Carter et al. 2011. A periodogram of the residuals from the phase curve fit reveals an additional stellar variability signal from the host star with a characteristic period of $0.620276\pm0.000011$ days and a full amplitude of $24\pm2$ ppm. We also present new Keck/HIRES radial velocity observations which we use to measure the orbit and obtain a mass ratio of $q=0.106\pm0.012$. Combining this measurement with the results of a stellar isochrone analysis, we find that the masses of the host star and white dwarf companion are $2.23\pm0.12\,M_{\odot}$ and $0.236^{+0.028}_{-0.027}\,M_{\odot}$, respectively. The effective temperatures of the two components are $9940^{+260}_{-230}$ K and $15080\pm400$ K, respectively, and we determine the age of the system to be $0.21^{+0.11}_{-0.08}$ Gyr. We use the measured system properties to compute predicted phase curve amplitudes and find that while the measured Doppler boosting and mutual illumination components agree well with theory, the ellipsoidal distortion amplitude is significantly underestimated. We detail possible explanations for this discrepancy, including interactions between the dynamical tide of the host star and the tidal bulge and possible non-synchronous rotation of the host star.

The data reduction procedure for radio interferometers can be viewed as a combined calibration and imaging problem. We present an algorithm that unifies cross-calibration, self-calibration, and imaging. Being a Bayesian method, that algorithm does not only calculate an estimate of the sky brightness distribution, but also provides an estimate of the joint uncertainty which entails both the uncertainty of the calibration and the one of the actual observation. The algorithm is formulated in the language of information field theory and uses Metric Gaussian Variational Inference (MGVI) as the underlying statistical method. So far only direction-independent antenna-based calibration is considered. This restriction may be released in future work. An implementation of the algorithm is contributed as well.

Turbulent motions in the convective envelope of red giants excite a rich spectrum of solar-like oscillation modes. Observations by CoRoT and Kepler have shown that the mode amplitudes increase dramatically as the stars ascend the red giant branch, i.e., as the frequency of maximum power, $\nu_\mathrm{max}$, decreases. Most studies nonetheless assume that the modes are well described by the linearized fluid equations. We investigate to what extent the linear approximation is justified as a function of stellar mass $M$ and $\nu_\mathrm{max}$, focusing on dipole mixed modes with frequency near $\nu_\mathrm{max}$. A useful measure of a mode's nonlinearity is the product of its radial wavenumber and its radial displacement, $k_r \xi_r$ (i.e., its shear). We show that $k_r \xi_r \propto \nu_\mathrm{max}^{-9/2}$, implying that the nonlinearity of mixed modes increases significantly as a star evolves. The modes are weakly nonlinear ($k_r \xi_r > 10^{-3}$) for $\nu_\mathrm{max} \lesssim 150 \, \mu\mathrm{Hz}$ and strongly nonlinear ($k_r \xi_r > 1$) for $\nu_\mathrm{max} \lesssim 30 \, \mu\mathrm{Hz}$, with only a mild dependence on $M$ over the range we consider ($1.0 - 2.0 M_\odot$). A weakly nonlinear mixed mode can excite secondary waves in the stellar core through the parametric instability, resulting in enhanced, but partial, damping of the mode. By contrast, a strongly nonlinear mode breaks as it propagates through the core and is fully damped there. Evaluating the impact of nonlinear effects on observables such as mode amplitudes and linewidths requires large mode network simulations. We plan to carry out such calculations in the future and investigate whether nonlinear damping can explain why some red giants exhibit dipole modes with unusually small amplitudes, known as depressed modes.

The orbital decay of binaries containing a primary sub-giant or red giant star and a stellar or substellar companion is investigated. The tide raised in the primary by the companion leads to an exchange of angular momentum between the orbit and the stellar spin, causing the orbit to contract and the primary to spin up. The rate of orbital decay is computed including both the equilibrium tide, damped by turbulent viscosity in the convective envelope, and the dynamical tide, assumed to be a traveling internal-gravity wave in the radiative core. For close binaries, the tidal forcing period is expected to be much shorter than the eddy turnover timescale in the convective envelope, and the prescription for "reduced" viscosity is an important consideration. The dynamical tide tends to dominate for the closest orbits, while the equilibrium tide dominates for more distant orbits, with the crossover point depending on the stellar mass. The spin up of the primary to synchronous rotation occurs for sufficiently massive secondaries, and this greatly slows the orbital decay until the Darwin instability occurs. A parameter survey is presented for orbital decay as a function of primary and secondary mass, as well as turbulent viscosity prescription. These results are summarized with analytic formulae and numerical results for the age-dependent critical separation, $a_{\rm crit}$, inside of which orbital decay is rapid, and few systems are expected to be observed. The calculations of $a_{\rm crit}$ are compared with APOGEE binaries, as well as solar mass exoplanet host stars.

We study the impact of nonlinear bulk viscosity due to Urca reactions driven by tidally-induced fluid motion during binary neutron star inspiral. Fluid compression is computed for low radial order oscillation modes through an adiabatic, time-dependent solution of the mode amplitudes. Optically thin neutrino emission and heating rates are then computed from this adiabatic fluid motion. Calculations use direct and modified Urca reactions operating in a $M=1.4\, M_\odot$ neutron star, which is constructed using the Skyrme Rs equation of state. We find that the energy pumped into low order oscillation modes is not efficiently thermalized even by direct Urca reactions, with core temperatures reaching only $T \simeq 10^8\, {\rm K}$ during the inspiral. Although this is an order of magnitude larger than the heating due to shear viscosity considered by previous studies, it reinforces the result that the stars are quite cold at merger. Upon excitation of the lowest order g-mode, the chemical potential imbalance reaches $\beta > 1\, {\rm MeV}$ at orbital frequencies $\nu_{\rm orb} > 200\, {\rm Hz}$, implying significant charged-current optical depths and Fermi blocking. To asses the importance of neutrino degeneracy effects, the neutrino transfer equation is solved in the static approximation for the three-dimensional density distribution, and the reaction rates are then computed including Fermi-blocking. We find that the heating rate is suppressed by factors of a few for $\nu_{\rm orb} > 200\, {\rm Hz}$. The spectrum of emitted $\nu_e$ and $\bar{\nu}_e$, including radiation transfer effects, is presented for a range of orbital separations.

Motivated by the discovery of a handful of pulsating, extremely low mass white dwarfs (ELM WDs, mass $M \lesssim 0.17\, M_\odot$) which likely have WD companions, this paper discusses binary formation models for these systems. Formation of an ELM WD by unstable mass transfer (MT) or a common envelope (CE) event is unlikely. Stable Roche-lobe overflow with conservative MT produces only $M \gtrsim 0.2\, M_\odot$. This paper discusses the formation of ELM WD using angular momentum losses by magnetic braking. The initially more massive star becomes the companion WD through a CE event. The initially less massive star becomes the ELM WD as an evolved donor in a Cataclysmic Variable binary. Evolutionary models are constructed using the Modules for Experiments in Stellar Astrophysics (MESA), with ELM WD progenitors in the range $M_d=1.0-1.5\, M_\odot$ and WD companions in the range $M_{\rm a} \simeq 0.4-0.9\, M_\odot$. A prescription to reduce magnetic braking for thin surface convection zones is included. Upon the thinning of the evolved donor's envelope, the donor star shrinks out of contact and MT ceases, revealing the ELM WD. Systems with small masses have previously been suggested as possible AM CVN's. Systems with large masses, up to the limit $M \simeq 0.18\, M_\odot$ at which shell flashes occur on the WD cooling track, will tend to expand out to orbital periods $P_{\rm orb} \gtrsim 15\, {\rm hr}$. In between this range, ELM WD may become pulsators both as pre-WD and on the WD cooling track. Brickhill's criterion for convective mode driving is used to estimate the location of the blue edge of the g-mode instability strip.

Data from radio interferometers provide a substantial challenge for statisticians. It is incomplete, noise-dominated and originates from a non-trivial measurement process. The signal is not only corrupted by imperfect measurement devices but also from effects like fluctuations in the ionosphere that act as a distortion screen. In this paper we focus on the imaging part of data reduction in radio astronomy and present RESOLVE, a Bayesian imaging algorithm for radio interferometry in its new incarnation. It is formulated in the language of information field theory. Solely by algorithmic advances the inference could be sped up significantly and behaves noticeably more stable now. This is one more step towards a fully user-friendly version of RESOLVE which can be applied routinely by astronomers.

This paper presents a detailed hydrostatic model of the upper atmosphere of HD 189733b, with the goal of constraining its temperature, particle densities, and radiation field over the pressure range $10^{-4}-10\, \mu \rm bar$, where the observed H$\alpha$ transmission spectrum is produced. The atomic hydrogen level population is computed including both collisional and radiative transition rates. The Ly$\alpha$ resonant scattering is computed using a Monte-Carlo simulation. The model transmission spectra are in broad agreement with the data. Excitation of the H(2$\ell$) population is mainly by Ly$\alpha$ radiative excitation due to the large Ly$\alpha$ intensity. The density of H(2$\ell$) is nearly flat over two decades in pressure, and is optically thick to H$\alpha$. Additional models computed for a range of the stellar Lyman continuum (LyC) flux suggest that the variability in H$\alpha$ transit depth may be due to the variability in the stellar LyC. Since metal lines provide the dominant cooling of this part of the atmosphere, the atmosphere structure is sensitive to the density of species such as Mg and Na which may themselves be constrained by observations. Since the H$\alpha$ and Na D lines have comparable absorption depths, we argue that the center of the Na D lines are also formed in the atomic layer where the H$\alpha$ line is formed.

WASP-12 is a hot Jupiter system with an orbital period of $P= 1.1\textrm{ day}$, making it one of the shortest-period giant planets known. Recent transit timing observations by Maciejewski et al. (2016) and Patra et al. (2017) find a decreasing period with $P/|\dot{P}| = 3.2\textrm{ 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_\ast \simeq 1.3 M_\odot$ main sequence star or a $M_\ast \simeq 1.2 M_\odot$ subgiant. We find that if the star is on the main sequence, the tidal dissipation is too inefficient to explain the observed $\dot{P}$. 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_\ast'\simeq 2\times10^5$ and an orbital decay rate that agrees well with the observed $\dot{P}$. It would also explain why the planet survived for $\simeq 3\textrm{ 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 $\sim 0.1\%$ of the planet's life, the probability of such a detection is a few percent given the observed sample of $\simeq 30$ hot Jupiters in $P<3\textrm{ day}$ orbits around $M_\ast>1.2 M_\odot$ hosts.

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 millisecond 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\times$10$^6$ M$_\odot$ pc$^{-3}$, a core radius of 0.16 pc, a pulsar density profile $n\propto r^{-3.14}$, and a total mass of M$_{\rm T}$($R_\perp<$1.0 pc)$\simeq3.0\times$10$^5$ M$_\odot$ 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 Terzan 5 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 3.0$\times$10$^4$ M$_\odot$.

Some close-in gaseous exoplanets are nearly in Roche-lobe contact, and previous studies show tidal decay can drive hot Jupiters into contact during the main sequence of their host stars. Improving upon a previous model, we present a revised model for mass transfer in a semi-detached binary system that incorporates an extended atmosphere around the donor and allows for an arbitrary mass ratio. We apply this new formalism to hypothetical, confirmed, and candidate planetary systems to estimate mass loss rates and compare with models of evaporative mass loss. Overflow may be significant for hot Neptunes out to periods of $\sim$ 2 days, while for hot Jupiters, it may only be important inward of 0.5 days. We find that CoRoT-24 b may be losing mass at a rate of more than an Earth mass in a Gyr. The hot Jupiter WASP-12 b may lose an Earth mass in a Myr, while the putative planet orbiting a T-Tauri star PTFO8-8695 might shed its atmosphere in a few Myrs. We point out that the orbital expansion that can accompany mass transfer may be less effective than previously considered because the gas accreted by the host star removes some of the system's angular momentum from the orbit, but simple scaling arguments suggest that the Roche-lobe overflow might remain stable. Consequently, the recently discovered small planets in ultra-short-periods ($<$ 1 day) may not be the remnants of hot Jupiters/Neptunes. The new model presented here has been incorporated into Modules for Experiments in Stellar Astrophysics (MESA).

Many gaseous exoplanets in short-period orbits are on the verge or are in the process of Roche-lobe overflow (RLO). Moreover, orbital stability analysis shows tides can drive many hot Jupiters to spiral inevitably toward their host stars. Thus, the coupled processes of orbital evolution and RLO likely shape the observed distribution of close-in exoplanets and may even be responsible for producing some of the short-period rocky planets. However, the exact outcome for an overflowing planet depends on its internal response to mass loss, and the accompanying orbital evolution can act to enhance or inhibit RLO. In this study, we apply the fully-featured and robust Modules for Experiments in Stellar Astrophysics (MESA) suite to model RLO of short-period gaseous planets. We show that, although the detailed evolution may depend on several properties of the planetary system, it is largely determined by the core mass of the overflowing gas giant. In particular, we find that the orbital expansion that accompanies RLO often stops and reverses at a specific maximum period that depends on the core mass. We suggest that RLO may often strand the remnant of a gas giant near this orbital period, which provides an observational prediction that can corroborate the hypothesis that short-period gas giants undergo RLO. We conduct a preliminary comparison of this prediction to the observed population of small, short-period planets and find some planets in orbits that may be consistent with this picture. To the extent that we can establish some short-period planets are indeed the remnants of gas giants, that population can elucidate the properties of gas giant cores, the properties of which remain largely unconstrained.

It has been suggested that hot stellar wind gas in a bow shock around an exoplanet is sufficiently opaque to absorb stellar photons and give rise to an observable transit depth at optical and UV wavelengths. In the first part of this paper, we use the CLOUDY plasma simulation code to model the absorption from X-ray to radio wavelengths by 1-D slabs of gas in coronal equilibrium with varying densities ($10^{4}-10^{8} \, {\rm cm^{-3}}$) and temperatures ($2000-10^{6} \ {\rm K}$) illuminated by a solar spectrum. For slabs at coronal temperatures ($10^{6} \ {\rm K}$) and densities even orders of magnitude larger than expected for the compressed stellar wind ($10^{4}-10^{5} \, {\rm cm^{-3}}$), we find optical depths orders of magnitude too small ($> 3\times10^{-7}$) to explain the $\sim3\%$ UV transit depths seen with Hubble. Using this result and our modeling of slabs with lower temperatures ($2000-10^4 {\rm K}$), the conclusion is that the UV transits of WASP-12b and HD 189733b are likely due to atoms originating in the planet, as the stellar wind is too highly ionized. A corollary of this result is that transport of neutral atoms from the denser planetary atmosphere outward must be a primary consideration when constructing physical models. In the second part of this paper, additional calculations using CLOUDY are carried out to model a slab of planetary gas in radiative and thermal equilibrium with the stellar radiation field. Promising sources of opacity from the X-ray to radio wavelengths are discussed, some of which are not yet observed.

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 Lyman-alpha 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.

A number of White Dwarfs (WDs) in cataclysmic binaries have shown brightness variations consistent with non-radial oscillations as observed in isolated WDs. A few objects have been well-characterized with photometric campaigns in the hopes of gleaning information about the mass, spin, and possibly internal structural characteristics. The novel aspect of this work is the possiblity to measure or constrain the interior structure and spin rate of WDs which have spent gigayears accreting material from their companion, undergoing thousands of nova outbursts in the process. In addition, variations in the surface temperature affect the site of mode driving, and provide unique and challenging tests for mode driving theories previously applied to isolated WD's. Having undergone long-term accretion, these WDs are expected to have been spun up. Spin periods in the range 60-100 seconds have been measured by other means for two objects, GW Lib and V455 And. Compared to typical mode frequencies, the spin frequency may be similar or higher, and the Coriolis force can no longer be treated as a small perturbation on the fluid motions. We present the results of a non-perturbative calculation of the normal modes of these WDs, using interior thermal structures appropriate to accreting systems. This includes a discussion of the surface brightness distributions, which are strongly modified from the non-rotating case. Using the measured spin period of approximately 100 seconds, we show that the observed pulsations from GW Lib are consistent with the three lowest azimuthal order rotationally modified modes that have the highest frequency in the stellar frame. The high frequencies are needed for the convective driving, but are then apparently shifted to lower frequencies by a combination of their pattern motion and the WD rotation.

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 $\omega_a$ excites a pair of secondary waves of frequency $\omega_b+\omega_c \simeq \omega_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 $\omega_b+\omega_c >>\omega_a$. We show that a p-mode can couple so strongly to a g-mode of similar radial wavelength that this type of 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_{gw} > 20 Hz and drive p-g pairs to significant energies on very short timescales (much less than the orbital decay time). 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 instability and therefore cannot say precisely how it influences 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: (i) induces significant orbital phase errors ($\Delta \phi$ > 1 radian) that accumulate primarily at low frequencies (f_{gw} < 50 Hz) and (ii) heats the neutron star core to T~10^{10} K. Since there are >100 unstable daughters, $\Delta \phi$ and T are potentially much larger than these values. Tides might therefore significantly influence the gravitational wave signal and electromagnetic emission from neutron star binaries at much larger orbital separations than previously thought.

We examine the dynamics of resonance locking in detached, tidally interacting binary systems. In a resonance lock, a given stellar or planetary mode is trapped in a highly resonant state for an extended period of time, during which the spin and orbital frequencies vary in concert to maintain the resonance. This phenomenon is qualitatively similar to resonance capture in planetary dynamics. We show that resonance locks can accelerate the course of tidal evolution in eccentric systems and also efficiently couple spin and orbital evolution in circular binaries. Previous analyses of resonance locking have not treated the mode amplitude as a fully dynamical variable, but rather assumed the adiabatic (i.e. Lorentzian) approximation valid only in the limit of relatively strong mode damping. We relax this approximation, analytically derive conditions under which the fixed point associated with resonance locking is stable, and further check these analytic results using numerical integrations of the coupled mode, spin, and orbital evolution equations. These show that resonance locking can sometimes take the form of complex limit cycles or even chaotic trajectories. We provide simple analytic formulae that define the binary and mode parameter regimes in which resonance locks of some kind occur (stable, limit cycle, or chaotic). We briefly discuss the astrophysical implications of our results for white dwarf and neutron star binaries as well as eccentric stellar binaries.

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$\alpha$ 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, mid-latitude wind zone, and equatorial dead zone -- found in previous analytic calculations is confirmed. For a magnetic field comparable to that of Jupiter, the equatorial dead zone, 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 wind zone is found to have a smaller contribution to the transit depth than the equatorial dead zone. Transmission spectra computed from the simulations are compared to HST STIS and ACS 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_0 > 10\ {\rm 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.

We calculate the tidal response of helium and carbon/oxygen (C/O) white dwarf (WD) binaries inspiraling due to gravitational wave emission. We show that resonance locks, previously considered in binaries with an early-type star, occur universally in WD binaries. In a resonance lock, the orbital and spin frequencies evolve in lockstep, so that the tidal forcing frequency is approximately constant and a particular normal mode remains resonant, producing efficient tidal dissipation and nearly synchronous rotation. We show that analogous locks between the spin and orbital frequencies can occur not only with global standing modes, but even when damping is so efficient that the resonant tidal response becomes a traveling wave. We derive simple analytic formulas for the tidal quality factor Q and tidal heating rate during a g-mode resonance lock, and verify our results numerically. We find that Q ~ 10^7 for orbital periods ~ 1 - 2 hr in C/O WDs, and Q ~ 10^9 for P_orb ~ 3 - 10 hr in helium WDs. Typically tidal heating occurs sufficiently close to the surface that the energy should be observable as surface emission. Moreover, near an orbital period of ~ 10 min, the tidal heating rate reaches ~ 10^{-2} L_\sun, rivaling the luminosities of our fiducial WD models. Recent observations of the 13-minute double-WD binary J0651 are roughly consistent with our theoretical predictions. Tides naturally tend to generate differential rotation; however, we show that the fossil magnetic field strength of a typical WD can maintain solid-body rotation down to at least P_orb ~ 10 min even in the presence of a tidal torque concentrated near the WD surface.

Absorption of stellar H\alpha\ by the upper atmosphere of the planet HD189733b 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\ ~ 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 HD189733b 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.

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 Msun 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 MESA's performance scaling 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 (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. [Abridged]

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 effects: three-mode nonlinear interactions and nonlinear excitation of modes by the time-varying gravitational potential of the companion. This paper presents the formalism for studying nonlinear tides and studies the nonlinear stability of the linear tidal flow. Although the formalism is applicable to binaries containing stars, planets, or compact objects, we focus on solar type stars with stellar or planetary companions. Our primary results include: (1) The linear tidal solution often used in studies of binary evolution is unstable over much of the parameter space in which it is employed. More specifically, resonantly excited gravity waves are unstable to parametric resonance for companion masses M' > 10-100 M_Earth at orbital periods P = 1-10 days. The nearly static equilibrium tide is, however, parametrically stable except for solar binaries with P < 2-5 days. (2) For companion masses larger than a few Jupiter masses, the dynamical tide causes waves to grow so rapidly that they must be treated as traveling waves rather than standing waves. (3) We find a novel form of parametric instability in which a single parent wave excites a very large number of daughter waves (N = 10^3[P / 10 days]) and drives them as a single coherent unit with growth rates that are ~N times faster than the standard three wave parametric instability. (4) Independent of the parametric instability, 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.

We investigate the pulsational stability of massive (M >~ 120 Msun) main sequence stars of a range of metallicities, including primordial, Population III stars. We include a formulation of convective damping motivated by numerical simulations of the interaction between convection and periodic shear flows. We find that convective viscosity is likely strong enough to stabilize radial pulsations whenever nuclear-burning (the epsilon-mechanism) is the dominant source of driving. This suggests that massive main sequence stars with Z <~ 2 x 10^-3 are pulsationally stable and are unlikely to experience pulsation-driven mass loss on the main sequence. These conclusions are, however, sensitive to the form of the convective viscosity and highlight the need for further high-resolution simulations of the convection-oscillation interaction. For more metal-rich stars (Z >~ 2 x 10^-3), the dominant pulsational driving arises due to the kappa-mechanism arising from the iron-bump in opacity and is strong enough to overcome convective damping. Our results highlight that even for oscillations with periods a few orders of magnitude shorter than the outer convective turnover time, the "frozen-in" approximation for the convection-oscillation interaction is inappropriate, and convective damping should be taken into account when assessing mode stability.

We develop a general framework for interpreting and analyzing high-precision lightcurves from eccentric stellar binaries. Although our methods are general, we focus on the recently discovered Kepler system KOI-54, a face-on binary of two A stars with $e=0.83$ and an orbital period of 42 days. KOI-54 exhibits strong ellipsoidal variability during its periastron passage; its lightcurve also contains ~20 pulsations at perfect harmonics of the orbital frequency, and another ~10 nonharmonic pulsations. Analysis of such data is a new form of asteroseismology in which oscillation amplitudes and phases rather than frequencies contain information that can be mined to constrain stellar properties. We qualitatively explain the physics of mode excitation and the range of harmonics expected to be observed. To quantitatively model observed pulsation spectra, we develop and apply a linear, tidally forced, nonadiabatic stellar oscillation formalism including the Coriolis force. We produce temporal power spectra for KOI-54 that are semi-quantitatively consistent with the observations. Both stars in the KOI-54 system are expected to be rotating pseudosynchronously, with resonant nonaxisymmetric modes providing a key contribution to the total torque; such resonances provide a possible explanation for the two largest-amplitude harmonic pulsations observed in KOI-54, although we find quantitative problems with this interpretation. We show in detail that the nonharmonic pulsations observed in KOI-54 can be produced by nonlinear three-mode coupling. The methods developed in this paper can be generalized in the future to determine the best-fit stellar parameters given pulsation data. We also derive an analytic model of KOI-54's ellipsoidal variability, including both tidal distortion and stellar irradiation, which can be used to model other similar systems.

Close-in, massive exoplanets raise significant tides in their stellar hosts. We compute the radial velocity (RV) signal due to this fluid motion in the equilibrium tide approximation. The predicted radial velocities in the observed sample of exoplanets exceed 1 m/s for 17 systems, with the largest predicted signal being ~30 m/s for WASP-18 b. Tidally-induced RV's are thus detectable with present methods. Both tidal fluid flow and the epicyclic motion of a slightly eccentric orbit produce an RV signal at twice the orbital frequency. If care is not taken, the tidally induced RV may, in some cases, be confused with a finite orbital eccentricity. Indeed, WASP-18 b is reported to have an eccentric orbit with small e=0.009 and pericenter longitude \omega=-\pi/2. Whereas such a close alignment of the orbit and line of sight to the observer requires fine tuning, this phase in the RV signal is naturally explained by the tidal velocity signature of an e=0 orbit. Additionally, the equilibrium tide estimate for the amplitude is in rough agreement with the data. Thus the reported eccentricity for WASP-18 b is instead likely a signature of the tidally-induced RV in the stellar host. Measurement of both the orbital and tidal velocity for non-transiting planets may allow planet mass and inclination to be separately determined solely from radial velocity data. We suggest that high precision fitting of RV data should include the tidal velocity signal in those cases where it may affect the determination of orbital parameters.

We have computed the fate of exoplanet companions around main sequence stars to explore the frequency of planet ingestion by their host stars during the red giant branch evolution. Using published properties of exoplanetary systems combined with stellar evolution models and Zahn's theory of tidal friction, we modeled the tidal decay of the planets' orbits as their host stars evolve. Most planets currently orbiting within 2 AU of their star are expected to be ingested by the end of their stars' red giant branch ascent. Our models confirm that many transiting planets are sufficiently close to their parent star that they will be accreted during the main sequence lifetime of the star. We also find that planet accretion may play an important role in explaining the mysterious red giant rapid rotators, although appropriate planetary systems do not seem to be plentiful enough to account for all such rapid rotators. We compare our modeled rapid rotators and surviving planetary systems to their real-life counterparts and discuss the implications of this work to the broader field of exoplanets.

(Abridged) The upper atmospheres of close-in gas giant exoplanets are subjected to intense heating/tidal forces from their parent stars. Atomic/ionized hydrogen (H) layers are sufficiently rarefied that magnetic pressure may dominate gas pressure for expected planetary magnetic field strength. We examine the magnetospheric structure using a 3D isothermal magnetohydrodynamic model that includes: a static "dead zone" near the magnetic equator containing magnetically confined gas; a "wind zone" outside the magnetic equator in which thermal pressure gradients and the magneto-centrifugal-tidal effect give rise to transonic outflow; and a region near the poles where sufficiently strong tidal forces may suppress transonic outflow. Using dipole field geometry, we estimate the size of the dead zone to be ~1-10 planetary radii for a range of parameters. To understand appropriate base conditions for the 3D isothermal model, we compute a 1D thermal model in which photoelectric heating from the stellar Lyman continuum is balanced by collisionally-excited Lyman {\alpha} cooling. This 1D model exhibits a H layer with temperatures T=5000-10000K down to pressures of 10-100 nbar. Using the 3D isothermal model, we compute H column densities and Lyman {\alpha} transmission spectra for parameters appropriate to HD 209458b. Line-integrated transit depths of 5-10% can be achieved for the above base conditions. Strong magnetic fields increase the transit signal while decreasing the mass loss, due to higher covering fraction and density of the dead zone. In our model, most of the transit signal arises from magnetically confined gas, some of which may be outside the L1 equipotential. Hence the presence of gas outside the L1 equipotential does not directly imply mass loss. Lastly, we discuss the domain of applicability for the magnetic wind model described in this paper and in the Roche-lobe overflow model.

Diffusive nuclear burning of H by an underlying material capable of capturing protons can readily consume H from the surface of neutron stars (NSs) during their early cooling history. In the absence of subsequent accretion, it will be depleted from the photosphere. We now extend diffusive nuclear burning to He, motivated by the recent observation by Ho \& Heinke of a carbon atmosphere on the NS in the Cassiopeia A supernova remnant. We calculate the equilibrium structure of He on an underlying $\alpha$ capturing material, accounting for thermal, mass defect, and Coulomb corrections on the stratification of material with the same zeroth order $\mu_e = A/Z$. We show that Coulomb corrections dominate over thermal and mass defect corrections in the highly degenerate part of the envelope. We also show that the bulk of the He sits deep in the envelope rather than near the surface. Thus, even if the photospheric He abundance is low, the total He column could be substantially larger than the photospheric column, which may have implications for rapid surface evolution ($\approx 1$ yr timescales) of neutron stars. When nuclear reactions are taken into account, we find that for base temperatures $\gtrsim 1.6 \times 10^8$ K, He is readily captured onto C. As these high temperatures are present during the early stages of NS evolution, we expect that the primordial He is completely depleted from the NS surface like the case for primordial H. We also find that magnetic fields $\lesssim 10^{12}$ G do not affect our conclusions. Armed with the results of this work and our prior efforts, we expect that primordial H and He are depleted, and so any observed H or He on the surfaces of these NS must be due to subsequent accretion (with or without spallation). If this subsequent accretion can be prevented, the underlying mid-Z material would be exposed.

Helium core white dwarfs (WDs) with mass M <~ 0.20 Msun undergo several Gyrs of stable hydrogen burning as they evolve. We show that in a certain range of WD and hydrogen envelope masses, these WDs may exhibit g-mode pulsations similar to their passively cooling, more massive carbon/oxygen core counterparts, the ZZ Cetis. Our models with stably burning hydrogen envelopes on helium cores yield g-mode periods and period spacings longer than the canonical ZZ Cetis by nearly a factor of 2. We show that core composition and structure can be probed using seismology since the g-mode eigenfunctions predominantly reside in the helium core. Though we have not carried out a fully nonadiabatic stability analysis, the scaling of the thermal time in the convective zone with surface gravity highlights several low-mass helium WDs that should be observed in search of pulsations: NLTT 11748, SDSS J0822+2753, and the companion to PSR J1012+5307. Seismological studies of these He core WDs may prove especially fruitful, as their luminosity is related (via stable hydrogen burning) to the hydrogen envelope mass, which eliminates one model parameter.

Asynchronous rotation and orbital eccentricity lead to time-dependent irradiation of the close-in gas giant exoplanets -- the hot Jupiters. This time-dependent surface heating gives rise to fluid motions which propagate throughout the planet. We investigate the ability of this "thermal tide" to produce a quadrupole moment which can couple to the stellar gravitational tidal force. While previous investigations discussed planets with solid surfaces, here we focus on entirely fluid planets in order to understand gas giants with small cores. The Coriolis force, thermal diffusion and self-gravity of the perturbations are ignored for simplicity. First, we examine the response to thermal forcing through analytic solutions of the fluid equations which treat the forcing frequency as a small parameter. In the "equilibrium tide" limit of zero frequency, fluid motion is present but does not induce a quadrupole moment. In the next approximation, finite frequency corrections to the equilibrium tide do lead to a nonzero quadrupole moment, the sign of which torques the planet {\it away} from synchronous spin. We then numerically solve the boundary value problem for the thermally forced, linear response of a planet with neutrally stratified interior and stably stratified envelope. The numerical results find quadrupole moments in agreement with the analytic non-resonant result at sufficiently long forcing period. Surprisingly, in the range of forcing periods of 1-30 days, the induced quadrupole moments can be far larger than the analytic result due to response of internal gravity waves which propagate in the radiative envelope. We discuss the relevance of our results for the spin, eccentricity and thermal evolution of hot Jupiters.

Motivated by the comments of Goodman (2009) on our paper concerning thermal tides (Arras and Socrates 2009a), we have studied an idealized problem to understand the global response of a completely fluid gas giant planet to thermal forcing at the surface (Arras and Socrates 2009b). Our findings disagree with the main claims in Goodman (2009). We find that significant quadrupole moments can indeed be induced as a result of thermal forcing. Furthermore, we find that it is possible for the orientation of the quadrupoles to be such that the planet is torqued away from synchronous rotation. Given these results, we believe our proposed thermal tide mechanism (Arras and Socrates 2009a) provides a viable scenario for generating steady-state asychronous rotation, inflated radii and possibly eccentric orbits of the hot Jupiters.

We study the orbits, tidal heating and mass loss from satellites around close-in gas giant exoplanets. The focus is on large satellites which are potentially observable by their transit signature. We argue that even Earth-size satellites around hot Jupiters may be immune to destruction by orbital decay; detection of such a massive satellite would strongly constrain theories of tidal dissipation in gas giants, in a manner complementary to orbital circularization. The star's gravity induces significant periodic eccentricity in the satellite's orbit. The resulting tidal heating rates, per unit mass, are far in excess of Io's and dominate radioactive heating out to planet orbital periods of months for reasonable satellite tidal $Q$. Inside planet orbital periods of about a week, tidal heating can completely melt the satellite. Lastly, we compute an upper limit to the satellite mass loss rate due to thermal evaporation from the surface, valid if the satellite's atmosphere is thin and vapor pressure is negligible. Using this upper limit, we find that although rocky satellites around hot Jupiters with orbital periods less than a few days can be significantly evaporated in their lifetimes, detectable satellites suffer negligible mass loss at longer orbital periods.

Rapid rotation in field red giant stars is a relatively rare but well-studied phenomenon; here we investigate the potential role of planet accretion in spinning up these stars. Using Zahn's theory of tidal friction and stellar evolution models, we compute the decay of a planet's orbit into its evolving host star and the resulting transfer of angular momentum into the stellar convective envelope. This experiment assesses the frequency of planet ingestion and rapid rotation on the red giant branch (RGB) for a sample of 99 known exoplanet host stars. We find that the known exoplanets are indeed capable of creating rapid rotators; however, the expected fraction due to planet ingestion is only ~10% of the total seen in surveys of present-day red giants. Of the planets ingested, we find that those with smaller initial semimajor axes are more likely to create rapid rotators because these planets are accreted when the stellar moment of inertia is smallest. We also find that many planets may be ingested prior to the RGB phase, contrary to the expectation that accretion would generally occur when the stellar radii expand significantly as giants. Finally, our models suggest that the rapid rotation signal from ingested planets is most likely to be seen on the lower RGB, which is also where alternative mechanisms for spin-up, e.g., angular momentum dredged up from the stellar core, do not operate. Thus, rapid rotators on the lower RGB are the best candidates to search for definitive evidence of systems that have experienced planet accretion.

Time-dependent insolation in a planetary atmosphere induces a mass quadrupole upon which the stellar tidal acceleration can exert a force. This "thermal tide" force can give rise to secular torques on the planet and orbit as well as radial forces causing eccentricity evolution. We apply this idea to the close-in gas giant exoplanets ("hot Jupiters"). The response of radiative atmospheres is computed in a hydrostatic model which treats the insolation as a time-dependent heat source, and solves for thermal radiation using flux-limited diffusion. Fully nonlinear numerical simulations are compared to solutions of the linearized equations, as well as analytic approximations, all of which are in good agreement. We find generically that thermal tide density perturbations {\it lead} the semi-diurnal forcing. As a result thermal tides can generate asynchronous spin and eccentricity. Our results are as follows: (1) Departure from synchronous spin is significant for hot Jupiters, and increases with orbital period. (2) Ongoing gravitational tidal dissipation in spin equilibrium leads to steady-state internal heating rates up to $\sim 10^{28} {\rm erg\ s^{-1}}$. If deposited sufficiently deep, these heating rates may explain the anomalously large radii of many hot Jupiters in terms of a "tidal main sequence" where cooling balances tidal heating. At fixed stellar type, planet mass and tidal $Q$, planetary radius increases strongly toward the star inside orbital periods $\la 2$ weeks. (3) There exists a narrow window in orbital period where small eccentricities, $e$, grow exponentially with a large rate. This window may explain the $\sim 1/4$ of hot Jupiters which should have been circularized by the gravitational tide long ago, but are observed to have significant nonzero $e$.(Abridged)

We study a non-dissipative hydrodynamical mechanism that can stabilize the spin of the accretor in an ultra-compact double white dwarf binary. This novel synchronization mechanism relies on a nonlinear wave interaction spinning down the background star. The essential physics of the synchronization mechanism is summarized as follows. As the compact binary coalesces due to gravitational wave emission, the largest star eventually fills its Roche lobe and accretion starts. The accretor then spins up due to infalling material and eventually reaches a spin frequency where a normal mode of the star is resonantly driven by the gravitational tidal field of the companion. If the resonating mode satisfies a set of specific criteria, which we elucidate in this paper, it exchanges angular momentum with the background star at a rate such that the spin of the accretor locks at this resonant frequency, even though accretion is ongoing. Some of the accreted angular momentum that would otherwise spin up the accretor is fed back into the orbit through this resonant tidal interaction. Two modes capable of stabilizing the accretor's spin are the l=4,m=2 and l=5,m=3 CFS unstable hybrid r-modes, which stabilize the spin of the accretor at frequency 2.6 and 1.5 times the binary's orbital frequency respectively. Since the stabilization mechanism relies on continuously driving a mode at resonance, its lifetime is limited since eventually the mode amplitude saturates due to non-linear mode-mode coupling. Rough estimates of the lifetime of the effect lie from a few orbits to millions of years.

Hundreds of substellar companions to solar-type stars will be discovered with the Kepler satellite. Kepler's extreme photometric precision gives access to low-amplitude stellar variability contributed by a variety of physical processes. We discuss in detail the periodic flux modulations arising from the tidal force on the star due to a substellar companion. An analytic expression for the variability is derived in the equilibrium-tide approximation. We demonstrate analytically and through numerical solutions of the linear, nonadiabatic stellar oscillation equations that the equilibrium-tide formula works extremely well for stars of mass <1.4 Msun with thick surface convection zones. More massive stars with largely radiative envelopes do not conform to the equilibrium-tide approximation and can exhibit flux variations $\ga$10 times larger than naive estimates. Over the full range of stellar masses considered, we treat the oscillatory response of the convection zone by adapting a prescription that A. J. Brickhill developed for pulsating white dwarfs. Compared to other sources of periodic variability, the ellipsoidal lightcurve has a distinct dependence on time and system parameters. We suggest that ellipsoidal oscillations induced by giant planets may be detectable from as many as ~100 of the 10^5 Kepler target stars. (Abridged)

(Abridged) The Cataclysmic Variable (CV) population harbors a diverse range of donor stars and accreting white dwarfs (WDs). A range of WD masses is expected, from low mass Helium core WDs, to massive WDs which have previously accreted at rates high enough for Hydrogen to burn steadily. Furthermore, a wide range of Helium enrichment is expected in the accreted material depending on the degree to which the donor star is evolved. We investigate the impact of this diversity on the range of effective temperatures ($T_{\rm eff}$) for which g-modes are unstable. The critical $T_{\rm eff}$ below which modes are unstable ("blue edge") depends on both surface gravity, $g$, and He abundance, $Y$. The Hydrogen/first Helium ionization instability strip is more sensitive to $g$ than $Y$. We find that (for solar composition envelopes), relative to a fiducial WD mass $0.6 M_\odot$, the blue edge for a $0.4 M_\odot$ He core WD shifts downward by $\approx 1000 {\rm K}$, while that for a massive $\approx 1.2 M_\odot$ WD shifts upward by $\approx 2000 {\rm K}$. The second Helium ionization instability strip exhibits strong dependences on both $g$ and $Y$. Surprisingly, increasing $Y$ by only 10% relative to solar creates an instability strip near $15,000 {\rm K}$. Hence CV's below the period gap with evolved donor stars of $Y\ga 0.4$ may have an "intermediate" instability strip well outside of the isolated DA and DB variables. This "intermediate" instability strip also occurs for low mass He WD with solar composition envelopes. The lack of pulsations in CV's with $T_{\rm eff}$ in the pure Hydrogen ZZ Ceti instability strip is also easily explained.

Quasi-periodic oscillations (QPOs) in the X-ray lightcurves of accreting neutron star and black hole binaries have been widely interpreted as being due to standing wave modes in accretion disks. These disks are thought to be highly turbulent due to the magnetorotational instability (MRI). We study wave excitation by MRI turbulence in the shearing box geometry. We demonstrate that axisymmetric sound waves and radial epicyclic motions driven by MRI turbulence give rise to narrow, distinct peaks in the temporal power spectrum. Inertial waves, on the other hand, do not give rise to distinct peaks which rise significantly above the continuum noise spectrum set by MRI turbulence, even when the fluid motions are projected onto the eigenfunctions of the modes. This is a serious problem for QPO models based on inertial waves.

Accretion flows with pressure gradients permit the existence of standing waves which may be responsible for observed quasi-periodic oscillations (QPO's) in X-ray binaries. We present a comprehensive treatment of the linear modes of a hydrodynamic, non-self-gravitating, polytropic slender torus, with arbitrary specific angular momentum distribution, orbiting in an arbitrary axisymmetric spacetime with reflection symmetry. We discuss the physical nature of the modes, present general analytic expressions and illustrations for those which are low order, and show that they can be excited in numerical simulations of relativistic tori. The mode oscillation spectrum simplifies dramatically for near Keplerian angular momentum distributions, which appear to be generic in global simulations of the magnetorotational instability. We discuss our results in light of observations of high frequency QPO's, and point out the existence of a new pair of modes which can be in an approximate 3:2 ratio for arbitrary black hole spins and angular momentum distributions, provided the torus is radiation pressure dominated. This mode pair consists of the axisymmetric vertical epicyclic mode and the lowest order axisymmetric breathing mode.

We consider the thermal structure and radii of strongly irradiated gas giant planets over a range in mass and irradiating flux. The cooling rate of the planet is sensitive to the surface boundary condition, which depends on the detailed manner in which starlight is absorbed and energy redistributed by fluid motion. We parametrize these effects by imposing an isothermal boundary condition $T \equiv T_{\rm deep}$ below the photosphere, and then constrain $T_{\rm deep}$ from the observed masses and radii. We compute the dependence of luminosity and core temperature on mass, $T_{\rm deep}$ and core entropy, finding that simple scalings apply over most of the relevant parameter space. These scalings yield analytic cooling models which exhibit power-law behavior in the observable age range $0.1-10 {\rm Gyr}$, and are confirmed by time-dependent cooling calculations. We compare our model to the radii of observed transiting planets, and derive constraints on $T_{\rm deep}$. Only HD 209458 has a sufficiently accurate radius measurement that $T_{\rm deep}$ is tightly constrained; the lower error bar on the radii for other planets is consistent with no irradiation. More accurate radius and age measurements will allow for a determination of the correlation of $T_{\rm deep}$ with the equilibrium temperature, informing us about both the greenhouse effect and day-night asymmetries.

Recent observations show that the thermal X-ray spectra of many isolated neutron stars are featureless and in some cases (e.g., RX J1856.5-3754) well fit by a blackbody. Such a perfect blackbody spectrum is puzzling since radiative transport through typical neutron star atmospheres causes noticeable deviation from blackbody. Previous studies have shown that in a strong magnetic field, the outermost layer of the neutron star may be in a condensed solid or liquid form because of the greatly enhanced cohesive energy of the condensed matter. The critical temperature of condensation increases with the magnetic field strength, and can be as high as 10^6 K (for Fe surface at B \sim 10^{13} G or H surface at B \sim a few times 10^{14} G). Thus the thermal radiation can directly emerge from the degenerate metallic condensed surface, without going through a gaseous atmosphere. Here we calculate the emission properties (spectrum and polarization) of the condensed Fe and H surfaces of magnetic neutron stars in the regimes where such condensation may be possible. For a smooth condensed surface, the overall emission is reduced from the blackbody by less than a factor of 2. The spectrum exhibits modest deviation from blackbody across a wide energy range, and shows mild absorption features associated with the ion cyclotron frequency and the electron plasma frequency in the condensed matter. The roughness of the solid condensate (in the Fe case) tends to decrease the reflectivity of the surface, and make the emission spectrum even closer to blackbody. We discuss the implications of our results for observations of dim, isolated neutron stars and magnetars.

We follow the time dependent thermal evolution of a white dwarf (WD) undergoing sudden accretion in a dwarf nova outburst, using both simulations and analytic estimates. The post-outburst lightcurve clearly separates into early times when the WD flux is high, and late times when the flux is near the quiescent level. The break between these two regimes, occurring at a time of order the outburst duration, corresponds to a thermal diffusion wave reaching the base of the freshly accreted layer. Our principal result is that long after the outburst, the fractional flux perturbation about the quiescent flux decays as a power law with time (and {\it not} as an exponential). We use this result to construct a simple fitting formula that yields estimates for both the quiescent flux and the accreted column, i.e. the total accreted mass divided by WD surface area. The WD mass is not well constrained by the late time lightcurve alone, but it can be inferred if the accreted mass is known from observations. We compare our work with the well-studied outburst of WZ Sge, finding that the cooling is well described by our model, giving an effective temperature $T_{\rm eff}=14,500 {\rm K}$ and accreted column $\Delta y\approx10^6 {\rm g cm^{-2}}$, in agreement with the modeling of Godon et al. To reconcile this accreted column with the accreted mass inferred from the bolometric accretion luminosity, a large WD mass $\gtrsim1.1M_\odot$ is needed. Our power law result is a valuable tool for making quick estimates of the outburst properties. We show that fitting the late time lightcurve with this formula yields a predicted column within 20% of that estimated from our full numerical calculations.

We compute the rate of diffusive nuclear burning for hydrogen on the surface of a "magnetar" (Soft Gamma-Ray Repeater or Anomalous X-Ray Pulsar). We find that hydrogen at the photosphere will be burned on an extremely rapid timescale of hours to years, depending on composition of the underlying material. Improving on our previous studies, we explore the effect of a maximally thick "inert" helium layer, previously thought to slow down the burning rate. Since hydrogen diffuses faster in helium than through heavier elements, we find this helium buffer actually increases the burning rate for magnetars. We compute simple analytic scalings of the burning rate with temperature and magnetic field for a range of core temperature. We conclude that magnetar photospheres are very unlikely to contain hydrogen. This motivates theoretical work on heavy element atmospheres that are needed to measure effective temperature from the observed thermal emission and constrains models of AXPs that rely on magnetar cooling through thick light element envelopes.

Pulsation modes have recently been observed in a handful of white dwarf (WD) primaries of cataclysmic variables, allowing an interesting new probe into the structure of accreting WD's. We briefly discuss the seismology of these objects, how stellar properties may be inferred from the observed mode frequencies, and mode driving mechanisms. For one pulsator, GW Lib, we have shown that a WD mass $M=1.05 M_\odot$ and accreted envelope mass $M_{\rm env}=0.4\times 10^{-4} M_\odot$ give the best match to the observed pulsation periods. A first exploration of mode driving favors $T_{\rm eff} = 14 000 K$ and a massive WD, but more work is necessary.

We present a first analysis of the g-mode oscillation spectrum for the white dwarf (WD) primary of GW Lib, a faint cataclysmic variable (CV). Stable periodicities have been observed from this WD for a number of years, but their interpretation as stellar pulsations has been hampered by a lack of theoretical models appropriate to an accreting WD. Using the results of Townsley and Bildsten, we construct accreting models for the observed effective temperature and approximate mass of the WD in GW Lib. We compute g-mode frequencies for a range of accreted layer masses, Macc, and long term accretion rates, <Mdot>. If we assume that the observed oscillations are from l=1 g-modes, then the observed periods are matched when M ~= 1.02 Msun, Macc ~= 0.31 x 10^-4 Msun and <Mdot> ~= 7.3 x 10^-11 Msun/yr. Much more sensitive observations are needed to discover more modes, after which we will be able to more accurately measure these parameters and constrain or measure the WD's rotation rate.

We present calculations of magnetic field evolution by the Hall effect and Ohmic decay in neutron star crusts. We first discuss when and for which neutron stars the Hall effect is important relative to Ohmic decay. We compute the evolution of an initial field distribution by Ohmic decay, and give approximate analytic formulas for both the surface and interior fields as a function of time. We then discuss the subsequent evolution by the Hall effect, giving several examples which illustrate how the field configuration evolves. Hall wave eigenfunctions are computed, including the effect of the large density change across the crust. We estimate the response of the crust to the magnetic stresses induced by Hall waves, and give a detailed discussion of the boundary conditions at the solid-liquid interface. Finally, we discuss the implications for the Hall cascade proposed by Goldreich & Reisenegger.

We calculate the coupled thermal evolution and magnetic field decay in relativistic model neutron stars threaded by superstrong magnetic fields (B > 10^{15} G). Our main goal is to evaluate how such ``magnetars'' evolve with time and how field decay modifies the transitions to core superfluidity and cooling dominated by surface X-ray emission. Observations of a thermal X-ray spectral component and fast timing noise place strong constraints on the presence of a superfluid core. We find that the transition to core superfluidity can be significantly delayed by field decay in the age range ~ 10^3-10^5 yrs. The mechanism of Hall drift is related to the stability of the core magnetic field, and to currents flowing outward through the crust. The heating effect is enhanced if it is continuous rather than spasmodic. Condensation of a heavy element layer at the surface is shown to cause only modest changes in the outward conduction of heat.

We describe a numerical method to solve the magnetohydrodynamic (MHD) equations. The fluid variables are updated along each direction using the flux conservative, 2nd order, total variation diminishing (TVD), upwind scheme of Jin and Xin. The magnetic field is updated separately in two-dimensional advection-constraint steps. The electromotive force (EMF) is computed in the advection step using the TVD scheme, and this same EMF is used immediately in the constraint step in order to preserve \grad.B=0 without the need to store intermediate fluxes. Operator splitting is used to extend the code to three dimensions, and Runge-Kutta is used to get second order accuracy in time. The advantages of this code are high resolution per grid cell, second order accuracy in space and time, enforcement of the \grad.B=0 constraint to machine precision, no memory overhead, speed, and simplicity. A 3-D Fortran implementation less than 400 lines long is made freely available. We also implemented a fully scalable message-passing parallel MPI version. We present tests of the code on MHD waves and shocks.

Rossby waves (r-modes) in rapidly rotating neutron stars are unstable because of the emission of gravitational radiation. We study saturation of this instability by nonlinear transfer of energy to stellar "inertial" oscillation modes. We present detailed calculations of stellar inertial modes in the WKB limit, their linear damping by bulk and shear viscosity, and the nonlinear coupling forces among these modes. The saturation amplitude is derived in the extreme limits of strong or weak driving by radiation reaction, as compared to the damping rate of low order inertial modes. We find the saturation energy is {\it extremely small}, at least four orders of magnitude smaller than that found by previous investigators. We discuss the consequences of this result for spin evolution of young neutron stars, and neutron stars being spun up by accretion in Low Mass X-ray Binaries.We also discuss the detection of these gravitational waves by LIGO.

We develop the formalism required to study the nonlinear interaction of modes in rotating Newtonian stars in the weakly nonlinear regime. The formalism simplifies and extends previous treatments. At linear order, we elucidate and extend slightly a formalism due to Schutz, show how to decompose a general motion of a rotating star into a sum over modes, and obtain uncoupled equations of motion for the mode amplitudes under the influence of an external force. Nonlinear effects are added perturbatively via three-mode couplings. We describe a new, efficient way to compute the coupling coefficients, to zeroth order in the stellar rotation rate, using spin-weighted spherical harmonics. We apply this formalism to derive some properties of the coupling coefficients relevant to the nonlinear interactions of unstable r-modes in neutron stars, postponing numerical integrations of the coupled equations of motion to a later paper. From an astrophysical viewpoint, the most interesting result of this paper is that many couplings of r-modes to other rotational modes (modes with zero frequencies in the non-rotating limit) are small: either they vanish altogether because of various selection rules, or they vanish to lowest order in the angular velocity. In zero-buoyancy stars, the coupling of three r-modes is forbidden entirely and the coupling of two r-modes to one hybrid rotational mode vanishes to zeroth order in rotation frequency. In incompressible stars, the coupling of any three rotational modes vanishes to zeroth order in rotation frequency.

We study spectroscopically determined iron abundances of 642 solar-type stars to search for the signature of accreted iron-rich material. We find that the metallicity [Fe/H] of a subset of 466 main sequence stars, when plotted as a function of stellar mass, mimics the pattern seen in lithium abundances in open clusters. Using Monte Carlo models we find that, on average, these stars have accreted about 0.4 Earth masses of iron while on the main sequence. A much smaller sample of 19 stars in the Hertzsprung gap, which are slightly evolved and whose convection zones are significantly more massive, have lower average [Fe/H], and their metallicity shows no clear variation with stellar mass. These findings suggest that terrestrial-type material is common around solar type stars.

Rossby waves (r-modes) have been suggested as a means to regulate the spin periods of young or accreting neutron stars, and also to produce observable gravitational wave radiation. R-modes involve primarily transverse, incompressive motions of the star's fluid core. However, neutron stars gain crusts early in their lives: therefore, r-modes also imply shear in the fluid beneath the crust. We examine the criterion for this shear layer to become turbulent, and derive the rate of dissipation in the turbulent regime. Unlike dissipation from a viscous boundary layer, turbulent energy loss is nonlinear in mode energy and can therefore cause the mode to saturate at amplitudes typically much less than unity. This energy loss also reappears as heat below the crust. We study the possibility of crust melting as well as its implications for the spin evolution of low-mass X-ray binaries. Lastly, we identify some universal features of the spin evolution that may have observational consequences.

We establish constraints on the mass and abundance of black holes in the Galactic halo by determining their impact on globular clusters which are conventionally considered to be little evolved. Using detailed Monte Carlo simulations and simple analytic estimates, we conclude that, at Galactocentric radius R~8 kpc, black holes with masses M_bh >~(1-3) x 10^6 M_sun can comprise no more than a fraction f_bh ~ 0.025-0.05 of the total halo density. This constraint significantly improves those based on disk heating and dynamical friction arguments as well as current lensing results. At smaller radius, the constraint on f_bh strengthens, while, at larger radius, an increased fraction of black holes is allowed.

We study neutrino-nucleon scattering and absorption in a dense, magnetized nuclear medium. These are the most important sources of neutrino opacity governing the cooling of a proto-neutron star in the first tens of seconds after its formation. Because the weak interaction is parity violating, the absorption and scattering cross-sections depend asymmetrically on the directions of the neutrino momenta with respect to the magnetic field. We develop the moment formalism of neutrino transport in the presence of such asymmetric opacities and derive explicit expressions for the neutrino flux and other angular moments of the Boltzmann transport equation. For a given neutrino species, there is a drift flux of neutrinos along the magnetic field in addition to the usual diffusive flux. This drift flux depends on the deviation of the neutrino distribution function from thermal equilibrium. Hence, despite the fact that the neutrino cross-sections are asymmetric throughout the star, asymmetric neutrino flux can be generated only in the outer region of the proto-neutron star where the neutrino distribution deviates significantly from thermal equilibrium. In addition to the asymmetric absorption opacity arising from nucleon polarization, we find the contribution of the electron (or positron) ground state Landau level. For neutrinos of energy less than a few times the temperature, this is the dominant source of asymmetric opacity. Lastly, we discuss the implication of our result to the origin of pulsar kicks: in order to generate kick velocity of a few hundred km/s from asymmetric neutrino emission using the parity violation effect, the proto-neutron star must have a dipole magnetic field of at least $10^{15}-10^{16}$ G.

We explore the possibility that massive black holes comprise a significant fraction of the dark matter of our galaxy by studying the dissolution of galactic globular clusters bombarded by them. In our simulations, we evolve the clusters along a sequence of King models determined by changes of state resulting from collisions with the black holes. The results divide naturally into regimes of `small' and `large' black hole mass. `Small' black holes do not destroy clusters in single collisions; their effect is primarily cumulative, leading to a relation between $\mbh$ and $\fhalo$, the fraction of the halo in black holes of mass $\mbh$, which is $\fhalo\mbh < $ constant (up to logarithmic corrections). For $\fhalo=1$, we find $\mbh \simless 10^{3} \msun$ by requiring survival of the same clusters studied by Moore (1993), who neglected cluster evolution, mass loss, and stochasticity of energy inputs in his estimates, but reached a similar conclusion. `Large' black holes may not penetrate a cluster without disrupting it; their effect is mainly catastrophic (close collisions), but also partly cumulative (distant collisions). In the large $\mbh$ limit, $\fhalo$ (but not $\mbh$) can be constrained by computing the probability that a cluster survives a combination of close, destructive encounters and distant, nondestructive encounters. We find that it is unlikely that $\fhalo \simgreat 0.3$ by requiring 50 per cent survival probability for Moore's clusters over $10^{10}$ years.

In magnetized proto-neutron stars, neutrino cross sections depend asymmetrically on the neutrino momenta due to parity violation. However, these asymmetric opacities do not induce any asymmetric flux in the bulk interior of the star where neutrinos are nearly in thermal equilibrium. Consequently, parity violation in neutrino absorption and scattering can only give rise to asymmetric neutrino flux above the neutrino-matter decoupling layer. The kick velocity is substantially reduced from previous estimates, requiring a dipole field $B \sim 10^{16}$~G to get $v_{kick}$ of order a few hundred km~s$^{-1}$.