Planetary systems with sufficiently small orbital spacings can experience planetary mergers and ejections. The branching ratio of mergers vs ejections depends sensitively on the treatment of planetary close encounters. Previous works have adopted a simple "sticky-sphere" prescription, whose validity is questionable. We apply both smoothed particle hydrodynamics and $N$-body integrations to investigate the fluid effects in close encounters between gas giants and the long-term evolution of closely-packed planetary systems. Focusing on parabolic encounters between Jupiter-like planets with $M_J$ and $2M_J$, we find that quick mergers occur when the impact parameter $r_p$ (the pericenter separation between the planets) is less than $2R_J$, and the merger conserved 97% of the initial mass. Strong tidal effects can affect the "binary-planet" orbit when $r_p$ is between $2R_J$ and $4R_J$. We quantify these effects using a set of fitting formulae that can be implemented in $N$-body codes. We run a suite of $N$-body simulations with and without the formulae for systems of two giant planets initially in unstable, nearly circular coplanar orbits. The fluid (tidal) effects significantly increase the branching ratio of planetary mergers relative to ejections by doubling the effective collision radius. While the fluid effects do not change the distributions of semi-major axis and eccentricity of each type of remnant planets (mergers vs surviving planets in ejections), the overall orbital properties of planet scattering remnants are strongly affected due to the change in the branching ratio. We also find that the merger products have broad distributions of spin magnitudes and obliquities.

Several white dwarfs with atmospheric metal pollution have been found to host small planetary bodies (planetesimals) orbiting near the tidal disruption radius. We study the physical properties and dynamical origin of these bodies under the hypothesis that they underwent high-eccentricity migration from initial distances of several astronomical units. We examine two plausible mechanisms for orbital migration and circularization: tidal friction and ram-pressure drag in a compact disc. For each mechanism, we derive general analytic expressions for the evolution of the orbit that can be rescaled for various situations. We identify the physical parameters that determine whether a planetesimal's orbit can circularize within the appropriate time-scale and constrain these parameters based on the properties of the observed systems. For tidal migration to work, an internal viscosity similar to that of molten rock is required, and this may be naturally produced by tidal heating. For disc migration to operate, a minimal column density of the disc is implied; the inferred total disc mass is consistent with estimates of the total mass of metals accreted by polluted WDs.

Large planetary spin-orbit misalignments (obliquities) may strongly influence atmospheric circulation and tidal heating in the planet. A promising avenue to generate obliquities is via spin-orbit resonances, where the spin and orbital precession frequencies of the planet cross each other as the system evolves in time. One such mechanism involves a dissipating (mass-losing) protoplanetary disk that drives orbital precession of an interior planet. We study this scenario analytically in this paper, and obtain the mapping between the general initial spin orientation and the final obliquity. We show that (i) under adiabatic evolution (i.e. the disk dissipates at a sufficiently slow rate), the final planetary obliquity as a function of the initial spin orientation bifurcates into distinct tracks governed by interactions with the resonance; (ii) under nonadiabatic evolution, a broad range of obliquities can be excited. We obtain analytical expressions for the final obliquities for various regimes of parameter space. The dynamical system studied in this paper is an example of "Colombo's top", and our analysis and results can be adapted to other applications.

We study the dynamics of merging compact binaries near a rotating supermassive black hole (SMBH) in a hierarchical triple configuration. We include various general relativistic effects that couple the inner orbit, the outer orbit and the spin of the SMBH. During the orbital decay due to gravitational radiation, the inner binary can encounter an "apsidal precession resonance" and experience eccentricity excitation. This resonance occurs when the apsidal precession rate of the inner binary matches that of the outer binary, with the precessions driven by both Newtonian interactions and various post-Newtonian effects. The eccentricity excitation requires the outer orbit to have a finite eccentricity, and is most effective for triples with small mutual inclinations, in contrast to the well-studied Lidov-Kozai effect. The resonance and the associated eccentricity growth may occur while the binary emits gravitational waves in the low-frequency band, and may be detectable by future space-based gravitational wave detectors.

In binaries composed of either early-type stars or white dwarfs, the dominant tidal process involves the excitation of internal gravity waves (IGWs), which propagate towards the stellar surface, and their dissipation via nonlinear wave breaking. We perform 2D hydrodynamical simulations of this wave breaking process in a stratified, isothermal atmosphere. We find that, after an initial transient phase, the dissipation of the IGWs naturally generates a sharp critical layer, separating the lower stationary region (with no mean flow) and the upper "synchronized" region (with the mean flow velocity equal to the horizontal wave phase speed). While the critical layer is steepened by absorption of these waves, it is simultaneously broadened by Kelvin-Helmholtz instabilities such that, in steady state, the critical layer width is determined by the Richardson criterion. We study the absorption and reflection of incident waves off the critical layer and provide analytical formulae describing its long-term evolution. The result of this study is important for characterizing the evolution of tidally heated white dwarfs and other binary stars.

In planetary systems with sufficiently small inter-planet spacing, close encounters can lead to planetary collisions/mergers or ejections. We study the spin property of the merger products of two giant planets in a statistical manner using numerical simulations and analytical modeling. Planetary collisions lead to rapidly rotating objects and a broad range of obliquities. We find that, under typical conditions for two-planet scatterings, the distributions of spin magnitude $S$ and obliquity $\Theta_{\rm SL}$ of the merger products have simple analytical forms: $f_{S} \propto S$ and $f_{\cos\Theta_{\rm SL}} \propto (1-\cos^2\Theta_{\rm SL})^{-1/2}$. Though parameter studies, we determine the regime of validity for the analytical distributions of spin and obliquity. Since planetary mergers is a major outcome of planet-planet scatterings, observational search for the spin/obliquity signatures of exoplanets would provide important constraints on the dynamical history of planetary systems.

The CHIME/FRB collaboration recently reported the detection of a 16 day periodicity in the arrival times of radio bursts from FRB 180916.J0158+65. We study the possibility that the observed periodicity arises from free precession of a magnetized neutron star, and put constraints on different components of the star's magnetic fields. Using a simple geometric model, where radio bursts are emitted from a rotating neutron star magnetosphere, we show that the emission pattern as a function of time can match that observed from FRB 180916.J0158+65.

While the radiation mechanism of fast radio bursts (FRBs) is unknown, coherent curvature radiation and synchrotron maser are promising candidates. We find that both radiation mechanisms work for a neutron star (NS) central engine with $B\gtrsim 10^{12}$ G, while for the synchrotron maser, the central engine can also be an accreting black hole (BH) with $B\gtrsim 10^{12}$ G and a white dwarf (WD) with $B\sim 10^8-10^9$ G. We study the electromagnetic counterparts associated with such central engines, i.e., nebulae for repeating FRBs and afterglows for non-repeating FRBs. In general, the energy spectrum and flux density of the counterpart depend strongly on its size and total injected energy. We apply the calculation to the nebula of FRB 121102 and find that the persistent radio counterpart requires the average energy injection rate into the nebula to be between $2.7\times10^{39}~{\rm erg/s}$ and $1.5\times10^{44}~{\rm erg/s}$, and the minimum injected energy be $6.0\times10^{47}~{\rm erg}$ in around $7$ yr. Consequently, we find that for FRB 121102 and its nebula: (1) WD and accretion BH central engines are disfavored; (2) a rotation-powered NS central engine works when $1.2\times10^{12}~{\rm G}\lesssim B\lesssim 7.8\times10^{14}~{\rm G}$ with initial period $P<180$ ms, but the radio emission must be more efficient than that in typical giant pulses of radio pulsars; and (3) a magnetic-powered NS central engine works when its internal magnetic field $B\gtrsim 10^{16}$ G. We also find that the radio-emitting electrons in the nebula could produce a significant rotation measure (RM), but cannot account for the entire observed RM of FRB 121102.

We carry out 2D viscous hydrodynamics simulations of circumbinary disk (CBD) accretion using {\footnotesize AREPO}. We resolve the accretion flow from a large-scale CBD down to the streamers and disks around individual binary components. Extending our recent studies \citep{mun19}, we consider circular binaries with various mass ratios ($0.1\leq q_{\rm{b}}\leq1$) and study accretion from ``infinite'', steady-supply disks and from finite-sized, viscously spreading tori. For ``infinite'' disks, a global steady state can be reached, and the accretion variability has a dominant frequency ${\sim}0.2\Omega_{\rm{b}}$ for $q_{\rm{b}}>0.5$ and $\Omega_{\rm{b}}$ for $q_{\rm{b}}<0.5$, ($\Omega_{\rm{b}}$ is the binary angular frequency). We find that the accretion ``eigenvalue'' $l_0$ -- the net angular momentum transfer from the disk to the binary per unit accreted mass -- is always positive and falls in the range ($0.65$-$0.85)a_{\rm b}^2\Omega_{\rm{b}}$ (with $a_{\rm{b}}$ the binary separation), depending weakly on the mass ratio and viscosity. This leads to binary expansion when $q_{\rm{b}}\gtrsim0.3$. Accretion from a finite torus can be separated into two phases: an initial transient phase, corresponding to the filling of the binary cavity, followed by a viscous pseudo-stationary phase, during which the torus viscously spreads and accretes onto the binary. In the viscous phase, the net torque on the binary per unit accreted mass is close to $l_0$, the value derived for ``infinite'' disks. We conclude that similar-mass binaries accreting from CBDs gain angular momentum and expand over long time scales. This result significantly impacts the coalescence of supermassive binary black holes and newly formed binary stars. We offer a word of caution against conclusions drawn from simulations of transient accretion onto empty circumbinary cavities.

Tidal dissipation due to convective turbulent viscosity shapes the evolution of a variety of astrophysical binaries. For example, this type of dissipation determines the rate of orbital circularization in a binary with a post-main sequence star that is evolving toward a common envelope phase. Viscous dissipation can also influence binaries with solar-type stars, or stars with a close-in giant planet. In general, the effective viscosity in a convective stellar envelope depends on the tidal forcing frequency $\omega_{\rm tide}$; when $\omega_{\rm tide}$ is larger than the turnover frequency of convective eddies, the viscosity is reduced. Previous work has focused on binaries in nearly circular orbits. However, for eccentric orbits, the tidal potential has many forcing frequencies. In this paper, we develop a formalism for computing tidal dissipation that captures the effect of frequency-dependent turbulent viscosity and is valid for arbitrary binary eccentricities. We also present an alternative simpler formulation that is suitable for very high eccentricities. We apply our formalisms to a red giant branch (RGB) star model and a solar-type star model. We find that a range of pseudosynchronous rotation rates are possible for both stellar models, and the pseudosynchronous rate can differ from the prediction of the commonly-used weak tidal friction theory by up to a factor of a few. We also find that tidal decay and circularization due to turbulent viscosity can be a few orders of magnitude faster than predicted by weak tidal friction in RGB stars on eccentric, small pericentre orbits, but is suppressed by a few orders of magnitude in solar-type stars due to viscosity reduction.

We present 2D hydrodynamical simulations of hot Jupiters orbiting near the inner edge of protoplanetary discs. We systemically explore how the accretion rate at the inner disc edge is regulated by a giant planet of different mass, orbital separation and eccentricity. We find that a massive (with planet-to-star mass ratio $\gtrsim 0.003$) eccentric ($e_p\gtrsim 0.1$) planet drives a pulsed accretion at the inner edge of the disc, modulated at one or two times the planet's orbital frequency. The amplitude of accretion variability generally increases with the planet mass and eccentricity, although some non-monotonic dependences are also possible. Applying our simulation results to the T Tauri system CI Tau, where a young hot Jupiter candidate has been detected, we show that the observed luminosity variability in this system can be explained by pulsed accretion driven by an eccentric giant planet.

Many warm Jupiters (WJs) have substantial eccentricities, which are linked to their formation and migration histories. This paper explores eccentricity excitation of WJs due to planet-planet scattering, beginning with 3-4 planets in unstable orbits, with the innermost planet placed in the range (0.1 - 1)AU. Such a setup is consistent with either in-situ formation or arrival at sub-AU orbits due to disk migration. Most previous N-body experiments have focused on "cold" Jupiters at several AU, where scattering results in planet ejections, efficiently exciting the eccentricities of surviving planets. In contrast, scattering at sub-AU distances results in a mixture of collisions and ejections, and the final eccentricities of surviving planets are unclear. We conduct scattering experiments for a range of planet masses and initial spacings, including the effect of general relativistic apsidal precession, and systematically catalogue the scattering outcomes and properties of surviving planets. A comparable number of one-planet and two-planet systems are produced. Two-planet systems arise exclusively through planet-planet collisions, and tend to have low eccentricities/mutual inclinations and compact configurations. One-planet systems arise through a combination of ejections and collisions, resulting in higher eccentricities. The observed eccentricity distribution of solitary WJs (lacking detection of a giant planet companion) is consistent with roughly 60% of the systems having undergone in-situ scattering, and the remaining experiencing a quiescent history.

We study the general relativitic (GR) effects induced by a spinning supermassive black hole on the orbital and spin evolution of a merging black hole binary (BHB) in a hierarchical triple system. A sufficiently inclined outer orbit can excite Lidov-Kozai eccentricity oscillations in the BHB and induce its merger. These GR effects generate extra precessions on the BHB orbits and spins, significantly increasing the inclination window for mergers and producing a wide range of spin orientations when the BHB enters LIGO band. This "GR-enhanced" channel may play an important role in BHB mergers.

We develop a simplified model for studying the long-term evolution of giant planets in protoplanetary discs. The model accounts for the eccentricity evolution of the planets and the dynamics of eccentric discs under the influences of secular planet-disc interactions and internal disc pressure, self-gravity and viscosity. Adopting the ansatz that the disc precesses coherently with aligned apsides, the eccentricity evolution equations of the planet-disc system reduce to a set of linearized ODEs, which allows for fast computation of the evolution of planet-disc eccentricities over long timescales. Applying our model to "giant planet + external disc" systems, we are able to reproduce and explain the secular behaviours found in previously published hydrodynamical simulations. We re-examine the possibility of eccentricity excitation (due to secular resonance) of multiple planets embedded in a dispersing disc, and find that taking into account the dynamics of eccentric discs can significantly affect the evolution of the planets' eccentricities.

We study the dynamical signatures of black hole (BH) and neutron star (NS) binary mergers via Lidov-Kozai oscillations induced by tertiary companions in hierarchical triple systems. For each type of binaries (BH-BH and BH-NS), we explore a wide range of binary/triple parameters that lead to binary mergers, and determine the distributions of eccentricity ($e_\mathrm{m}$) and spin-orbit misalignment angle ($\theta_\rm {sl}^\rm f$) when the binary enters the LIGO/VIRGO band. We use the double-averaged (over both orbits) and single-averaged (over the inner orbit) secular equations, as well as N-body integration, to evolve systems with different hierarchy levels, including the leading-order post-Newtonian effect, spin-orbit coupling and gravitational radiation. We find that for merging BH-BH binaries with comparable masses, about $7\%$ have $e_\mathrm{m}>0.1$ and $0.7\%$ have $e_{\rm m}>0.9$. The majority of the mergers have significant eccentricities in the LISA band. The BH spin evolution and $\theta_\rm {sl}^\rm f$ are correlated with the orbital evolution and $e_{\rm m}$. Mergers with $e_{\rm m} \lesssim 10^{-3}$ have a distribution of $\theta_\rm {sl}^\rm f$ that peaks around $90^\circ$ (and thus favoring a projected binary spin parameter $\chi_{\rm eff}\sim 0$), while mergers with larger $e_{\rm m}$ have a more isotropic spin-orbit misalignments. For typical BH-NS binaries, strong octuple effects lead to more mergers with non-negligible $e_\mathrm{m}$ (with $\sim 18\%$ having $e_\mathrm{m}>0.1$ and $2.5\%$ having $e_\mathrm{m}>0.9$), and the final BH spin axis tends to be randomly orientated. Measurements or constraints on eccentric mergers and $\theta_\rm {sl}^\rm f$ from LIGO/VIRGO and LISA would provide useful diagnostics on the dynamical formation of merging BH or NS binaries in triples.

Dynamically formed compact object binaries may still be eccentric while in the LIGO/Virgo band. For a neutron star (NS) in an eccentric binary, the fundamental modes (f-modes) are excited at pericenter, transferring energy from the orbit to oscillations in the NS. We model this system by coupling the evolution of the NS f-modes to the orbital evolution of the binary as it circularizes and moves toward coalescence. NS f-mode excitation generally speeds up the orbital decay and advances the phase of the gravitational wave signal from the system. We calculate how this effect changes the timing of pericenter passages and examine how the cumulative phase shift before merger depends on the initial eccentricity of the system. This phase shift can be much larger for highly eccentric mergers than for circular mergers, and can be used to probe the NS equation of state.

After the Tidal Disruption Event (TDE) of a star around a SuperMassive Black Hole (SMBH), the bound stellar debris rapidly forms an accretion disk. If the accretion disk is not aligned with the spinning SMBH's equatorial plane, the disk will be driven into Lense-Thirring precession around the SMBH's spin axis, possibly affecting the TDE's light curve. We carry out an eigenmode analysis of such a disk to understand how the disk's warp structure, precession, and inclination evolution are influenced by the disk's and SMBH's properties. We find an oscillatory warp may develop as a result of strong non-Keplarian motion near the SMBH. The global disk precession frequency matches the Lense-Thirring precession frequency of a rigid disk around a spinning black hole within a factor of a few when the disk's accretion rate is high, but deviates significantly at low accretion rates. Viscosity aligns the disk with the SMBH's equatorial plane over timescales of days to years, depending on the disk's accretion rate, viscosity, and SMBH's mass. We also examine the effect of fall-back material on the warp evolution of TDE disks, and find that the fall-back torque aligns the TDE disk with the SMBH's equatorial plane in a few to tens of days for the parameter space investigated. Our results place constraints on models of TDE emission which rely on the changing disk orientation with respect to the line of sight to explain observations.

Recent studies suggest that ultra-short period planets (USPs), Earth-sized planets with sub-day periods, constitute a statistically distinct sub-sample of {\it Kepler} planets: USPs have smaller radii ($1-1.4R_\oplus$) and larger mutual inclinations with neighboring planets than nominal {\it Kepler} planets, and their period distribution is steeper than longer-period planets. We study a "low-eccentricity" migration scenario for the formation of USPs, in which a low-mass planet with initial period of a few days maintains a small but finite eccentricity due to secular forcings from exterior companion planets, and experiences orbital decay due to tidal dissipation. USP formation in this scenario requires that the initial multi-planet system have modest eccentricities ($\gtrsim 0.1$) or angular momentum deficit. During the orbital decay of the inner-most planet, the system can encounter several apsidal and nodal precession resonances that significantly enhance eccentricity excitation and increase the mutual inclination between the inner planets. We develop an approximate method based on eccentricity and inclination eigenmodes to efficiently evolve a large number of multi-planet systems over Gyr timescales in the presence of rapid (as short as $\sim 100$~years) secular planet-planet interactions and other short-range forces. Through a population synthesis calculation, we demonstrate that the "low-$e$ migration" mechanism can naturally produce USPs from the large population of {\it Kepler} multis under a variety of conditions, with little fine tuning of parameters. This mechanism favors smaller inner planets with more massive and eccentric companion planets, and the resulting USPs have properties that are consistent with observations.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

Relativistic jets are observed only in the low/hard and intermediate states of X-ray binaries (XRBs), and are switched off in the thermal state, but they appear to be present in both low-luminosity and luminous active galactic nuclei (AGNs). It is widely believed that strong large-scale magnetic fields is a crucial ingredient in jet production; such fields can be attained only through efficient advection from the outer disc. We suggest that geometrically thin accretion discs with magnetic outflows are present in luminous radio-loud AGNs; this is likely because the interstellar medium provides both mass and sufficient magnetic flux to the outer disc. Most angular momentum of such disc is removed by the outflows, and the radial velocity of the disc is significantly increased compared to viscous drift velocity. This facilitates efficient magnetic field advection through the disc to produce a strong field near the black hole in luminous AGNs, which helps launch relativistic jets. In XRBs, the magnetic fields of the gas from companion stars are too weak to drive outflows from outer discs. Jets are therefore switched off in the thermal state due to inefficient magnetic field advection in the disc.

The population of giant planets on short-period orbits can potentially be explained by some flavours of high-eccentricity migration. In this paper we investigate one such mechanism involving "secular chaos", in which secular interactions between at least three giant planets push the inner planet to a highly eccentric orbit, followed by tidal circularization and orbital decay. In addition to the equilibrium tidal friction, we incorporate dissipation due to dynamical tides that are excited inside the giant planet. Using the method of Gaussian rings to account for planet-planet interactions, we explore the conditions for extreme eccentricity excitation via secular chaos and the properties of hot Jupiters formed in this migration channel. Our calculations show that once the inner planet reaches a sufficiently large eccentricity, dynamical tides quickly dissipate the orbital energy, producing an eccentric warm Jupiter, which then decays in semi-major axis through equilibrium tides to become a hot Jupiter. Dynamical tides help the planet avoid tidal disruption, increasing the chance of forming a hot Jupiter, although not all planets survive the process. We find that the final orbital periods generally lie in the range of 2-3 days, somewhat shorter than those of the observed hot Jupiter population. We couple the planet migration to the stellar spin evolution to predict the final spin-orbit misalignments. The distribution of the misalignment angles we obtain shows a lack of retrograde orbits compared to observations. Our results suggest that high-eccentricity migration via secular chaos can only account for a fraction of the observed hot Jupiter population.

We carry out 2D viscous hydrodynamical simulations of circumbinary accretion using the AREPO code. We self-consistently compute the accretion flow over a wide range of spatial scales, from the circumbinary disk (CBD) far from the central binary, through accretion streamers, to the disks around individual binary components, resolving the flow down to 2% of the binary separation. We focus on equal-mass binaries with arbitrary eccentricities. We evolve the flow over long (viscous) timescales until a quasi-steady is reached, in which the mass supply rate at large distances $\dot{M}_0$ (assumed constant) equals the time-averaged mass transfer rate across the disk and the total mass accretion rate onto the binary components. This quasi-steady state allows us to compute the secular angular momentum transfer rate onto the binary, $\langle\dot{J}_b\rangle$, and the resulting orbital evolution. Through direct computation of the gravitational and accretion torques on the binary, we find that $\langle\dot{J}_b\rangle$ is consistently positive (i.e., the binary gains angular momentum), with $l_0\equiv\langle\dot{J}_b\rangle/\dot M_0$ in the range of $(0.4-0.8)a_b^2\Omega_b$, depending on the binary eccentricity (where $a_b,~\Omega_b$ are the binary semi-major axis and angular frequency); we also find that this $\langle\dot{J}_b\rangle$ is equal to the net angular momentum current across the CBD, indicating that global angular momentum balance is achieved in our simulations. We compute the time-averaged rate of change of the binary orbital energy for eccentric binaries, and thus obtain the secular rates $\langle\dot a_b\rangle$ and $\langle \dot{e}_b\rangle$. In all cases, $\langle\dot{a}_b\rangle$ is positive, i.e., the binary expands while accreting. We discuss the implications of our results for the merger of supermassive binary black holes and for the formation of close stellar binaries.

We study the orbital evolution of black hole (BH) binaries in quadruple systems, where the tertiary binary excites large eccentricity in the BH binary through Lidov-Kozai (LK) oscillations, causing the binary BHs to merge via gravitational radiation. For typical BH binaries with masses $m_{1,2}\simeq 20M_\odot-30M_\odot$ and initial semimajor axis $a_0\sim100$ AU (such that the binaries have no chance of merging by themselves within $\sim10^{10}$ yrs), we show that binary-binary interactions can significantly increase the LK window for mergers (the range of companion inclinations that allows the BH binary to merge within 10~Gyrs). This increase arises from a secular resonance between the LK oscillation of the BH binary and the nodal precession of the outer (binary-binary) orbit driven by the tertiary binary. Therefore, in the presence of tertiary binary, the BH merger fraction is increased to $10-30\%$, an order of magnitude larger than the merger fraction found in similar triple systems. Our analysis (with appropriate scalings) can be easily adapted to other configurations of systems, such as relatively compact BH binaries and moderately hierarchical triples, which may generate even higher merger fractions. Since the occurrence rate of stellar quadruples in the galactic fields is not much smaller than that of stellar triples, our result suggests that dynamically induced BH mergers in quadruple systems may be an important channel of producing BH mergers observed by LIGO/VIRGO.

The prevalence of binary stars at close separations implies that many of these systems will interact or merge during the binary's lifetime. This paper presents hydrodynamic simulations of the scenario of binary coalescence through unstable mass transfer, which drives the pair to closer separations. When the donor star does not rotate synchronously with respect to the orbit, dynamical tidal waves are excited in its envelope. We show that resonance crossings with high azimuthal-order $(m\sim3$ to $6$) fundamental modes induce a visible "polygram" distortion to the star. As the binary orbit tightens, the system sweeps through resonance with modes of decreasing azimuthal order, which are selectively excited. We compare our hydrodynamic simulations to predictions from linear theory of resonant mode excitation. The linear theory provides an estimate of mode amplitudes to within a factor of two, even as the oscillations become quite non-linear as the stars coalesce. We estimate that a wave with 10% radial amplitude generates approximately 1% photometric variability; this may be detectible if such a binary coalescence is caught in action by future photometric all-sky surveys.

High-eccentricity migration is an important channel for the formation of hot Jupiters (HJs). In particular, Lidov-Kozai (LK) oscillations of orbital eccentricity/inclination induced by a distant planetary or stellar companion, combined with tidal friction, have been shown to produce HJs on Gyr timescales, provided that efficient tidal dissipation operates in the planet. We re-examine this scenario with the inclusion of dynamical tides. When the planet's orbit is in a high-eccentricity phase, the tidal force from the star excites oscillatory f-modes and r-modes in the planet. For sufficiently large eccentricity and small pericentre distance, the mode can grow chaotically over multiple pericentre passages and eventually dissipate non-linearly, drawing energy from the orbit and rapidly shrinking the semi-major axis. We study the effect of such chaotic tides on the planet's orbital evolution. We find that this pathway produces very eccentric ($e\gtrsim 0.9$) warm Jupiters (WJs) on short timescales (a few to 100 Myrs). These WJs efficiently circularize to become HJs due to their persistently small pericentre distances. Chaotic tides can also save some planets from tidal disruption by truncating the LK eccentricity oscillations, significantly increasing the HJ formation fraction for a range of planet masses and radii. Using a population synthesis calculation, we determine the characteristics of WJs and HJs produced in this scenario, including the final period distribution, orbital inclinations and stellar obliquities. Chaotic tides endow LK migration with several favorable features to explain observations of HJs. We expect that chaotic tides are also important in other flavours of high-$e$ migration.

Black hole (BH) mergers driven by gravitational perturbations of external companions constitute an important class of formation channels for merging BH binaries detected by LIGO. We have studied the orbital and spin evolution of binary BHs in triple systems, where the tertiary companion excites large eccentricity in the inner binary through Lidov-Kozai oscillations, causing the binary to merge via gravitational radiation. Using the single-orbital-averaged and double-orbital-averaged secular dynamics equations of triples, we determine the merger window (the range of companion inclinations that allows the inner binary to merge) and the merger fraction as a function of various system parameters. For typical BH binaries, the merger fraction increases rapidly with the perturber eccentricity $e_{\rm out}$ because of the octupole perturbation, reaching $10-20\%$ at $e_{\rm out}=0.9$. We derive the analytical expressions for the merger window and merger fraction for systems with negligible octupole effect, and apply them to neutron star binary mergers in triples. We also follow the spin evolution of the BHs during the companion-induced orbital decay. Starting from aligned spin axes, a wide range of final spin-orbit misalignment angle $\theta_{\rm sl}^{\rm f}$ can be generated when the binary enters the LIGO band. For systems where the octupole effect is small (such as equal mass binaries or $e_{\rm out}\sim 0$), the distribution of $\theta_{\rm sl}^{\rm f}$ peaks around $90^\circ$. As the octuple effect increases, a more isotropic distribution of final spin axis is produced. Overall, merging BH binaries produced by Lidov-Kozai oscillations in triples exhibit a unique distribution of the effective spin parameter $\chi_{\rm eff}$; this may be used to distinguish this formation channel from others.

Stellar spin-orbit misalignments (obliquities) in hot Jupiter systems have been extensively probed. Such obliquities may reveal clues about hot Jupiter dynamical histories. Common explanations for generating obliquities include high-eccentricity migration and primordial disk misalignment. This paper investigates another mechanism for producing stellar spin-orbit misalignments in systems hosting a close-in planet with an external, modestly inclined companion. Spin-orbit misalignment may be excited due to a secular resonance, occurring when the precession rate of the stellar spin axis (driven by the inner planet) becomes comparable to the nodal precession rate of the inner planet (driven by the companion). Due to the spin-down of the host star via magnetic braking, this resonance may be achieved during the star's main-sequence lifetime for a wide range of planet masses and orbital architectures. Obliquity excitation is accompanied by a decrease in mutual inclination between the inner planet and perturber, and can thus erase high inclinations. For hot Jupiters, the stellar spin axis is strongly coupled to the orbital axis, and obliquity excitation by a giant planet companion requires a strong perturber, usually located within 1-2 AU. For warm Jupiters, the spin and orbital axes are more weakly coupled, and the resonance may be achieved for distant giant planet perturbers (at several to tens of AU). Since warm Jupiters have a high occurrence rate of distant planetary companions with appropriate properties for resonant obliquity excitation, stellar obliquities in warm Jupiter systems may be common, particularly for warm Jupiters orbiting cool stars that have undergone significant spin-down.

The mass-period or radius-period distribution of close-in exoplanets shows a paucity of intermediate mass/size (sub-Jovian) planets with periods ~< 3 days. We show that this sub-Jovian desert can be explained by the photoevaporation of highly irradiated sub-Neptunes and the tidal disruption barrier for gas giants undergoing high-eccentricity migration. The distinctive triangular shape of the sub-Jovain desert result from the fact that photoevaporation is more effective closer to the host star, and that in order for a gas giant to tidally circularise closer to the star without tidal disruption it needs to be more massive. Our work indicates that super-Earths/mini-Neptunes and hot-Jupiters had distinctly separate formation channels and arrived at their present locations at different times.

A number of multiplanet systems are observed to contain planets very close to mean motion resonances, although there is no significant pileup of precise resonance pairs. We present theoretical and numerical studies on the outcome of capture into first-order mean motion resonances (MMRs) using a parametrized planet migration model that takes into account nonlinear eccentricity damping due to planet-disk interaction. This parametrization is based on numerical hydrodynamical simulations and is more realistic than the simple linear parametrization widely used in previous analytic studies. We find that nonlinear eccentricity damping can significantly influence the stability and outcome of resonance capture. In particular, the equilibrium eccentricity of the planet captured into MMRs become larger, and the captured MMR state tends to be more stable compared to the prediction based on the simple migration model. In addition, when the migration is sufficiently fast or/and the planet mass ratio is sufficiently small, we observe a novel phenomenon of eccentricity overshoot, where the planet's eccentricity becomes very large before settling down to the lower equilibrium value. This can lead to the ejection of the smaller planet if its eccentricity approaches unity during the overshoot. This may help explain the lack of low-mass planet companion of hot Jupiters when compared to warm Jupiters.

Several recent studies have suggested that circumstellar disks in young stellar binaries may be driven into misalignement with their host stars due to secular gravitational interactions between the star, disk and the binary companion. The disk in such systems is twisted/warped due to the gravitational torques from the oblate central star and the external companion. We calculate the disk warp profile, taking into account of bending wave propagation and viscosity in the disk. We show that for typical protostellar disk parameters, the disk warp is small, thereby justifying the "flat-disk" approximation adopted in previous theoretical studies. However, the viscous dissipation associated with the small disk warp/twist tends to drive the disk toward alignment with the binary or the central star. We calculate the relevant timescales for the alignment. We find the alignment is effective for sufficiently cold disks with strong external torques, especially for systems with rapidly rotating stars, but is ineffective for the majority of star-disk-binary systems. Viscous warp driven alignment may be necessary to account for the observed spin-orbit alignment in multi-planet systems if these systems are accompanied by an inclined binary companion.

Many hot Jupiter (HJ) systems have been observed to have their stellar spin axis misaligned with the planet's orbital angular momentum axis. The origin of this spin-orbit misalignment and the formation mechanism of HJs remain poorly understood. A number of recent works have suggested that gravitational interactions between host stars, protoplanetary disks, and inclined binary companions may tilt the stellar spin axis with respect to the disk's angular angular momentum axis, producing planetary systems with misaligned orbits. These previous works considered idealized disk evolution models and neglected the gravitational influence of newly formed planets. In this paper, we explore how disk photoevaporation and planet formation and migration affect the inclination evolution of planet-star-disk-binary systems. We take into account planet-disk interactions and the gravitational spin-orbit coupling between the host star and the planet. We find that the rapid depletion of the inner disk via photoevaporation reduces the excitation of stellar obliquities. Depending on the formation and migration history of HJs, the spin-orbit coupling between the star and the planet may reduces and even completely suppress the excitation of stellar obliquities. Our work constrains the formation/migration history of HJs. On the other hand, planetary systems with "cold" Jupiters or close-in super-earths may experience excitation of stellar obliquities in the presence of distant inclined companions.

Compact multi-planet systems containing super-Earths or sub-Neptunes, commonly found around solar-type stars, may be surrounded by external giant planet or stellar companions, which can shape the architechture and observability of the inner systems. We present a comprehensive study on the evolution of the inner planetary system subject to the gravitational influence of an eccentric, misaligned outer perturber. Analytic results are derived for the inner planet eccentricities ($e_i$) and mutual inclination ($\theta_{12}$) of the "2-planet + perturber" system, calibrated with numerical secular and N-body integrations, as a function of the perturber mass $m_p$, semi-major axis $a_p$ and inclination angle $\theta_p$. We find that the dynamics of the inner system is determined by the dimensionless parameter $\epsilon_{12}$, given by the ratio between the differential precession rate driven by the perturber and the mutual precession rate of the inner planets. Loosely packed systems (corresponding to $\epsilon_{12} \gg 1$) are more susceptible to eccentricity/inclination excitations by the perturber than tightly packed inner systems (with $\epsilon_{12} \ll 1$) (or singletons), although resonance may occur around $\epsilon_{12}\sim 1$, leading to large $e_i$ and $\theta_{12}$. Dynamical instability may set in for inner planet systems with large excited eccentricities and mutual inclinations. We present a formalism to extend our analytical results to general inner systems with $N>2$ planets and apply our results to constrain possible external companions to the Kepler-11 system. Eccentricity and inclination excitation by external companions may help explain the observational trend that systems with fewer transiting planets are dynamically hotter than those with more transiting planets.

Highly eccentric binary systems appear in many astrophysical contexts, ranging from tidal capture in dense star clusters, precursors of stellar disruption by massive black holes, to high-eccentricity migration of giant planets. In a highly eccentric binary, the tidal potential of one body can excite oscillatory modes in the other during a pericenter passage, resulting in energy exchange between the modes and the binary orbit. These modes exhibit one of three behaviors over multiple passages: low-amplitude oscillations, large amplitude oscillations corresponding to a resonance between the orbital frequency and the mode frequency, and chaotic growth. We study these phenomena with an iterative map, fully exploring how the mode evolution depends on the pericenter distance and other parameters. In addition, we show that the dissipation of mode energy results in a quasi-steady state, with gradual orbital decay punctuated by resonances, even in systems where the mode amplitude would initially grow stochastically. A newly captured star around a black hole can experience significant orbital decay and heating due to the chaotic growth of the mode amplitude and dissipation. A giant planet pushed into a high-eccentricity orbit may experience a similar effect and become a hot or warm Jupiter.

Consider a planet with its orbital angular momentum axis aligned with the spin axis of its host star. To what extent does an inclined distant companion (giant planet or binary star) affect this alignment? We provide an analytic, quantitative answer and apply it to hot Jupiter systems, for which misalignments between the orbital axis and the stellar spin axis have been detected. We also show how similar consideration can be applied to multi-planet systems with distant companions (such as Kepler-56). The result of this paper provides a simple method to assess the dynamical role played by external companions on spin-orbit misalignments in exoplanetary systems.

We investigate the migration of a low-mass ($\lesssim 10 M_\oplus$) planet near the inner edge of a protoplanetary disc using two-dimensional viscous hydrodynamics simulations. We employ an inner boundary condition representing the truncation of the disc at the stellar corotation radius. As described by Tsang (2011), wave reflection at the inner disc boundary modifies the Type I migration torque on the planet, allowing migration to be halted before the planet reaches the inner edge of the disc. For low-viscosity discs ($\alpha \lesssim 10^{-3}$), planets may be trapped with semi-major axes as large as $3-5$ times the inner disc radius. In general, planets are trapped closer to the inner edge as either the planet mass or the disc viscosity parameter $\alpha$ increases, and farther from the inner edge as the disc thickness is increased. This planet trapping mechanism may impact the formation and migration history of close-in compact multiplanet systems.

We study the chaotic and secular evolution of hierarchical quadruple systems in the $3+1$ configuration, focusing on the evolution of mutual inclination of the inner binaries as the system undergoes coupled Lidov-Kozai (LK) oscillations. We include short-range forces (SRF; such as those due to tidal and rotational distortions) that control the eccentricity excitation of the inner binary. The evolution of mutual inclination is described, a priori, by two dimensionless parameters, $\pazocal{R}_0$, the ratio between the inner and outer LK time-scales and $\epsilon_{SRF}$, the ratio between the SRF precession and the inner LK precession rates. We find that the chaotic zones for the mutual inclination depend mainly on $\pazocal{R}_0$, while $\epsilon_{SRF}$ controls mainly the range of eccentricity excitation. The mutual inclination evolves chaotically for $1\lesssim \pazocal{R}_0\lesssim 10$, leading to large misalignments. For $0.4 \lesssim \pazocal{R}_0 \lesssim 0.8$, the system could be weakly excited and produce bimodal distribution of mutual inclination angles. Our results can be applied to exomoons-planets in stellar binaries and Warm/Hot Jupiters in stellar triples. Such systems could develop large mutual inclination angles if the inner binary is tight enough, and also high eccentricities, depending of the strength of the short-range forces. Future detections of tilted Warm/Hot Jupiters and exomoons could put our mechanism under observational tests.

It is usually thought that viscous torque works to align a circumbinary disk with the binary's orbital plane. However, recent numerical simulations suggest that the disk may evolve to a configuration perpendicular to the binary orbit ("polar alignment") if the binary is eccentric and the initial disk-binary inclination is sufficiently large. We carry out a theoretical study on the long-term evolution of inclined disks around eccentric binaries, calculating the disk warp profile and dissipative torque acting on the disk. For disks with aspect ratio $H/r$ larger than the viscosity parameter $\alpha$, bending wave propagation effectively makes the disk precess as a quasi-rigid body, while viscosity acts on the disk warp and twist to drive secular evolution of the disk-binary inclination. We derive a simple analytic criterion (in terms of the binary eccentricity and initial disk orientation) for the disk to evolve toward polar alignment with the eccentric binary. When the disk has a non-negligible angular momentum compared to the binary, the final "polar alignment" inclination angle is reduced from $90^\circ$. For typical protoplanetary disk parameters, the timescale of the inclination evolution is shorter than the disk lifetime, suggesting that highly-inclined disks and planets may exist orbiting eccentric binaries.

We study the effect of external companion on the orbital and spin evolution of merging black-hole (BH) binaries. An sufficiently close by and inclined companion can excite Lidov-Kozai (LK) eccentricity oscillations in the binary, thereby shortening its merger time. During such LK-enhanced orbital decay, the spin axis of the BH generally exhibits chaotic evolution, leading to a wide range ($0^\circ$-$180^\circ$) of final spin-orbit misalignment angle from an initially aligned configuration. For systems that do not experience eccentricity excitation, only modest ($\lesssim 20^\circ$) spin-orbit misalignment can be produced, and we derive an analytic expression for the final misalignment using the principle of adiabatic invariance. The spin-orbit misalignment directly impacts the gravitational waveform, and can be used to constrain the formation scenarios of BH binaries and dynamical influences of external companions.

Recent studies have proposed that most warm Jupiters (WJs, giant planets with semi-major axes in the range of 0.1-1 AU) probably form in-situ, or arrive in their observed orbits through disk migration. However, both in-situ formation and disk migration, in their simplest flavors, predict WJs to be in low-eccentricity orbits, in contradiction with many observed WJs that are moderately eccentric (e=0.2-0.7). This paper examines the possibility that the WJ eccentricities are raised by secular interactions with exterior giant planet companions, following in-situ formation or migration on a circular orbit. Eccentricity growth may arise from an inclined companion (through Lidov-Kozai cycles), or from an eccentric, nearly coplanar companion (through apsidal precession resonances). We quantify the necessary conditions (in terms of the eccentricity, semi-major axis and inclination) for external perturbers of various masses to raise the WJ eccentricity. We also consider the sample of eccentric WJs with detected outer companions, and for each system, identify the range of mutual inclinations needed to generate the observed eccentricity. For most systems, we find that relatively high inclinations (at least $\sim 40^\circ$) are needed so that Lidov-Kozai cycles are induced; the observed outer companions are typically not sufficiently eccentric to generate the observed WJ eccentricity in a low-inclination configuration. The results of this paper place constraints on possibly unseen external companions to eccentric WJs. Observations that probe mutual inclinations of giant planet systems will help clarify the origin of eccentric WJs and the role of external companions.

In coalescing neutron star (NS) binaries, tidal force can resonantly excite low-frequency (< 500 Hz) oscillation modes in the NS, transferring energy between the orbit and the NS. This resonant tide can induce phase shift in the gravitational waveforms, and potentially provide a new window of studying NS interior using gravitational waves. Previous works have considered tidal excitations of pure g-modes (due to stable stratification of the star) and pure inertial modes (due to Coriolis force), with the rotational effect treated in an approximate manner. However, for realistic NSs, the buoyancy and rotational effects can be comparable, giving rise to mixed inertial-gravity modes. We develop a non-perturbative numerical spectral code to compute the frequencies and tidal coupling coefficients of these modes. We then calculate the phase shift in the gravitational waveform due to each resonance during binary inspiral. We adopt polytropic NS models with a parameterized stratification. We derive relevant scaling relations and survey how the phase shift depends on various properties of the NS. We find that for canonical NSs (with mass M = 1.4M_sun and radius R = 10 km) and modest rotation rates (< 300 Hz), the gravitational wave phase shift due to a resonance is generally less than 0.01 radian. But the phase shift is a strong function of R and M, and can reach a radian or more for low-mass NSs with larger radii (R > 15 km). Significant phase shift can also be produced when the combination of stratification and rotation gives rise to a very low frequency (< 20 Hz in the inertial frame) modified g-mode. We also find that some inertial modes can be strongly affected by stratification, and that the m = 1 r-mode, previously identified to have a small but finite inertial-frame frequency based on the Cowling approximation, in fact has essentially zero frequency, and therefore cannot be excited.

We study the secular dynamics of compact binaries (consisting of white dwarfs, neutron stars or black holes) with tertiary companions in hierarchical triple systems. As the inner binary (with initially negligible eccentricity) undergoes orbital decay due to gravitational radiation, its eccentricity can be excited by gravitational forcing from the tertiary. This excitation occurs when the triple system passes through an "apsidal precession resonance", when the precession rate of the inner binary, driven by the gravitational perturbation of the external companion and general relativity, matches the precession rate of the outer binary. The eccentricity excitation requires the outer companion to be on an eccentric orbit, with the mutual inclination between the inner and outer orbits less than $\sim 40^\circ$. Gravitational wave (GW) signals from the inner binary can be significantly modified as the system evolves through the apsidal precession resonance. For some system parameters (e.g., a white dwarf binary with a brown dwarf tertiary), the resonance can happen when the binary emits GWs in the $10^{-4}-10^{-1}$~Hz range (the sensitivity band of LISA).

Hierarchical quadruple systems arise naturally in stellar binaries and triples that harbour planets. Examples are hot Jupiters (HJs) in stellar triple systems, and planetary companions to HJs in stellar binaries. The secular dynamical evolution of these systems is generally complex, with secular chaotic motion possible in certain parameter regimes. The latter can lead to extremely high eccentricities and, therefore, strong interactions such as efficient tidal evolution. These interactions are believed to play an important role in the formation of HJs through high-eccentricity migration. Nevertheless, a deeper understanding of the secular dynamics of these systems is still lacking. Here we study in detail the secular dynamics of a special case of hierarchical quadruple systems in either the `2+2' or `3+1' configurations. We show how the equations of motion can be cast in a form representing a perturbed hierarchical three-body system, in which the outer orbital angular momentum vector is precessing steadily around a fixed axis. In this case, we show that eccentricity excitation can be significantly enhanced when the precession period is comparable to the Lidov-Kozai oscillation time-scale of the inner orbit. This arises from an induced large mutual inclination between the inner and outer orbits driven by the precession of the outer orbit, even if the initial mutual inclination is small. We present a simplified semi-analytic model that describes the latter phenomenon.

Recent SPH simulations by Martin et al. (2014) suggest a circumstellar gaseous disk may exhibit coherent eccentricity-inclination oscillations due to the tidal forcing of an inclined binary companion, in a manner that resembles Lidov-Kozai oscillations in hierarchical triple systems. We carry out linear stability analysis for the eccentricity growth of circumstellar disks in binaries, including the effects of gas pressure and viscosity and secular (orbital-averaged) tidal force from the inclined companion. We find that the growth of disk eccentricity depends on the dimensionless ratio ($S$) between $c_{\rm s}^2$ (the disk sound speed squared) and the tidal torque acting on the disk (per unit mass) from the companion. For $S\ll 1$, the standard Lidov-Kozai result is recovered for a thin disk annulus: eccentricity excitation occurs when the mutual inclination $I$ between the disk and binary lies between $39^\circ$ and $141^\circ$. As $S$ increases, the inclination window for eccentricity growth generally becomes narrower. For $S\gtrsim$ a few, eccentricity growth is suppressed for all inclination angles. Surprisingly, we find that for $S\sim 1$ and certain disk density/pressure profiles, eccentricity excitation can occur even when $I$ is much less than $39^\circ$.

Rocky planets orbiting M-dwarf stars in the habitable zone tend to be driven to synchronous rotation by tidal dissipation, potentially causing difficulties for maintaining a habitable climate on the planet. However, the planet may be captured into asynchronous spin-orbit resonances, and this capture may be more likely if the planet has a sufficiently large intrinsic triaxial deformation. We derive the analytic expression for the maximum triaxiality of a rocky planet, with and without a liquid envelope, as a function of the planet's radius, density, rigidity and critical strain of fracture. The derived maximum triaxiality is consistent with the observed triaxialities for terrestrial planets in the solar system, and indicates that rocky planets in the habitable zone of M-dwarfs can in principle be in a state of asynchronous spin-orbit resonances.

It is well known that asteroids and comets fall into the Sun. Metal pollution of white dwarfs and transient spectroscopic signatures of young stars like $\beta$-Pic provide growing evidence that extra solar planetesimals can attain extreme orbital eccentricities and fall onto their parent stars. We aim to develop a general, practically implementable, semi-analytical theory of secular eccentricity excitation of small bodies in mean motion resonances with an eccentric planet valid for arbitrary values of the eccentricities and including the short-range force due to General Relativity. Our semi-analytic model for the restricted planar three-body problem does not make use of any series expansion and therefore is valid for any values of eccentricities and semi-major axes ratios. The model is based on the application of the adiabatic principle, which is valid when the precession period of the longitude of pericenter of the planetesimal is much longer than the libration period in the mean motion resonance. This holds down to vanishingly small eccentricities in resonances of order larger than 1. We provide a Mathematica notebook with the implementation of the model allowing direct use to the interested reader. We confirm that the 4:1 mean motion resonance with a moderately eccentric planet is the most powerful one to lift the eccentricity of planetesimals from nearly circular orbits to star-grazing ones. However, if the planet is too eccentric, we find that this resonances becomes unable to pump the planetesimal's eccentricity to very high value. The inclusion of the General Relativity effect imposes a condition on the mass of the planet to drive the planetsimals into star-grazing orbits. For a planetesimal at $\sim$1 AU around a solar-mass star (or white dwarf), we find a threshold planetary mass of about 17 Earth masses. We finally derive an analytical formula for this critical mass.

Many protostellar gapped and binary discs show misalignments between their inner and outer discs; in some cases, $\sim70$ degree misalignments have been observed. Here we show that these misalignments can be generated through a "secular precession resonance" between the nodal precession of the inner disc and the precession of the gap-opening (stellar or massive planetary) companion. An evolving protostellar system may naturally cross this resonance during its lifetime due to disc dissipation and/or companion migration. If resonance crossing occurs on the right timescale, of order a few Myrs, characteristic for young protostellar systems, the inner and outer discs can become highly misaligned ($\gtrsim 60$ degrees). When the primary star has a mass of order a solar mass, generating a significant misalignment typically requires the companion to have a mass of $\sim 0.01-0.1$ M$_\odot$ and an orbital separation of tens of AU. The recently observed companion in the cavity of the gapped, highly misaligned system HD 142527 satisfies these requirements, indicating that a previous resonance crossing event misaligned the inner and outer discs. Our scenario for HD 142527's misaligned discs predicts that the companion's orbital plane is aligned with the outer disc's; this prediction should be testable with future observations as the companion's orbit is mapped out. Misalignments observed in several other gapped disc systems could be generated by the same secular resonance mechanism.

Recent observations of Kepler multi-planet systems have revealed a number of systems with planets very close to second-order mean motion resonances (MMRs, with period ratio $1:3$, $3:5$, etc.) We present an analytic study of resonance capture and its stability for planets migrating in gaseous disks. Resonance capture requires slow convergent migration of the planets, with sufficiently large eccentricity damping timescale $T_e$ and small pre-resonance eccentricities. We quantify these requirements and find that they can be satisfied for super-Earths under protoplanetary disk conditions. For planets captured into resonance, an equilibrium state can be reached, in which eccentricity excitation due to resonant planet-planet interaction balances eccentricity damping due to planet-disk interaction. We show that this "captured" equilibrium can be overstable, leading to partial or permanent escape of the planets from the resonance. In general, the stability of the captured state depends on the inner to outer planet mass ratio $q=m_1/m_2$ and the ratio of the eccentricity damping times. The overstability growth time is of order $T_e$, but can be much larger for systems close to the stability threshold. For low-mass planets undergoing type I (non-gap opening) migration, convergent migration requires $q \lesssim 1$, while the stability of the capture requires $q\gtrsim 1$. These results suggest that planet pairs stably captured into second-order MMRs have comparable masses. This is in contrast to first-order MMRs, where a larger parameter space exists for stable resonance capture. We confirm and extend our analytical results with $N$-body simulations, and show that for overstable capture, the escape time from the MMR can be comparable to the time the planets spend migrating between resonances.

Eclipsing binaries are observed to have a range of eccentricities and spin-orbit misalignments (stellar obliquities). Whether such properties are primordial, or arise from post-formation dynamical interactions remains uncertain. This paper considers the scenario in which the binary is the inner component of a hierarchical triple stellar system, and derives the requirements that the tertiary companion must satisfy in order to raise the eccentricity and obliquity of the inner binary. Through numerical integrations of the secular octupole-order equations of motion of stellar triples, coupled with the spin precession of the oblate primary star due to the torque from the secondary, we obtain a simple, robust condition for producing spin-orbit misalignment in the inner binary: In order to excite appreciable obliquity, the precession rate of the stellar spin axis must be smaller than the orbital precession rate due to the tertiary companion. This yields quantitative requirements on the mass and orbit of the tertiary. We also present new analytic expressions for the maximum eccentricity and range of inclinations allowing eccentricity excitation (Lidov-Kozai window), for stellar triples with arbitrary masses and including the non-Keplerian potentials introduced by general relativity, stellar tides and rotational bulges. The results of this paper can be used to place constraints on unobserved tertiary companions in binaries that exhibit high eccentricity and/or spin-orbit misalignment, and will be helpful in guiding efforts to detect external companions around stellar binaries. As an application, we consider the eclipsing binary DI Herculis, and identify the requirements that a tertiary companion must satisfy to produce the observed spin-orbit misalignment.

A white dwarf (WD) captured into a high-eccentricity orbit around a massive black hole (MBH) may undergo many pericenter passages before tidal disruption. During these passages, the tidal potential of the MBH excites internal oscillations or waves in the WD, and the dissipation of these oscillations can significantly influence the physical properties of the WD prior to its disruption. We calculate the amplitude of the tidally excited gravity (buoyancy) waves in the WD as a function of the pericenter distance and eccentricity for realistic WD models, under the assumption that these outgoing gravity waves are efficiently dissipated in outer layers of the WD by non-linear effects or radiative damping. We obtain fitting formulae for the tidal energy and angular momentum transfer rates as well as the tidal heating rate. We find that these dynamical tides are much weaker than gravitational radiation in driving the orbital decay of the WD-MBH binary, and they are also inefficient in changing the WD spin during the orbital evolution. Incorporating our computed tidal dissipation rate into a MESA-based WD evolution code, we find that tidal heating can lead to appreciable brightening of the WD and may induce runaway fusion in the hydrogen envelope well before the WD undergoes tidal disruption.

We carry out numerical simulations of circumbinary discs, solving the viscous hydrodynamics equations on a polar grid covering an extended disc outside the binary co-orbital region. We use carefully controlled outer boundary conditions and long-term integrations to ensure that the disc reaches a quasi-steady state, in which the time-averaged mass accretion rate onto the binary, $\langle\dot{M}\rangle$, matches the mass supply rate at the outer disc. We focus on binaries with comparable masses and a wide range of eccentricities ($e_\mathrm{B}$). For $e_\mathrm{B} \lesssim 0.05$, the mass accretion rate of the binary is modulated at about $5$ times the binary period; otherwise it is modulated at the binary period. The inner part of the circumbinary disc ($r \lesssim 6 a_\mathrm{B}$) generally becomes coherently eccentric. For low and high $e_\mathrm{B}$, the disc line of apsides precesses around the binary, but for intermediate $e_\mathrm{B}$ ($0.2 - 0.4$), it instead becomes locked with that of the binary. By considering the balance of angular momentum transport through the disc by advection, viscous stress, and gravitational torque, we determine the time-averaged net angular momentum transfer rate to the binary, $\langle\dot{J}\rangle$. The specific angular momentum, $l_0 = \langle\dot{J}\rangle/\langle\dot{M}\rangle$, depends non-monotonically on $e_\mathrm{B}$. Contrary to previous claims, we find that $l_0$ is positive for most $e_\mathrm{B}$, implying that the binary receives net angular momentum, which may cause its separation to grow with time. The minimum $l_0$ occurs at intermediate $e_\mathrm{B}$ ($0.2 - 0.4$), corresponding to the regime where the inner eccentric disc is apsidally aligned with the binary.

The {\it Kepler} mission has detected thousands of planetary systems with 1-7 transiting planets packed within 0.7~au from their host stars. There is an apparent excess of single-transit planet systems that cannot be explained by transit geometries alone, when a single planetary mutual inclination dispersion is assumed. This suggests that the observed compact planetary systems have at least two different architectures. We present a scenario where the "Kepler dichotomy" may be explained by the action of an external giant planet or stellar companion misaligned with the inner multi-planet system. The external companion excites mutual inclinations of the inner planets, causing such systems to appear as "Kepler singles" in transit surveys. We derive approximate analytic expressions (in various limiting regimes), calibrated with numerical calculations, for the mutual inclination excitations for various planetary systems and perturber properties (mass $m_p$, semi-major axis $a_p$ and inclination $\theta_p$). In general, the excited mutual inclination increases with $m_p/a_p^3$ and $\theta_p$, although secular resonances may lead to large mutual inclinations even for small $\theta_p$. We discuss the implications of our results for understanding the dynamical history of transiting planet systems with known external perturbers.

Newly formed planets (or brown dwarfs) may possess disks or rings that occupy an appreciable fraction of the planet's Hill sphere and extend beyond the Laplace radius, where the tidal torque from the host star dominates over the torque from the oblate planet. Such a disk/ring can exhibit unique, detectable transit signatures, provided that the disk/ring is significantly misaligned with the orbital plane of the planet. There exists tentative evidence for an extended ring system around the young K5 star 1 SWASP J140747-354542. We present a general theoretical study of the inclination (warp) profile of circumplanetary disks under the combined influences of the tidal torque from the central star, the torque from the oblate planet and the self-gravity of the disk. We calculate the steady-state warp profile ("generalized Laplace Surface") and investigate the condition for coherent precession of the disk. We find that to maintain non-negligible misalignment between the extended outer disk and the planet's orbital plane, and to ensure coherent disk precession, the disk surface density must be sufficiently large so that the self-gravity torque overcomes the tidal torque from the central star. Our analysis and quantitative results can be used to constrain the parameters of transiting circumplanetary disks that may be detected in the future.

Bailey et al.~(2016) and Gomes et al.~(2016) recently suggested that the 6 degree misalignment between the Sun's rotational equator and the orbital plane of the major planets may be produced by the forcing from the hypothetical Planet Nine on an inclined orbit. Here we present a simple but accurate calculation of the effect, which provides a clear description of how the Sun's spin orientation depends on the property of Planet Nine in this scenario.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

Many exoplanetary systems containing hot Jupiters (HJs) exhibit significant misalignment between the spin axes of the host stars and the orbital angular momentum axes of the planets ("spin-orbit misalignment"). High-eccentricity migration involving Lidov-Kozai oscillations of the planet's orbit induced by a distant perturber is a possible channel for producing such misaligned HJ systems. Previous works have shown that the dynamical evolution of the stellar spin axis during the high-$e$ migration plays a dominant role in generating the observed spin-orbit misalignment. Numerical studies have also revealed various patterns of the evolution of the stellar spin axis leading to the final misalignment. Here we develop an analytic theory to elucidate the evolution of spin-orbit misalignment during the Lidov-Kozai migration of planets in stellar binaries. Secular spin-orbit resonances play a key role in the misalignment evolution. We include the effects of short-range forces and tidal dissipation, and categorize the different possible paths to spin-orbit misalignment as a function of various physical parameters (e.g. planet mass and stellar rotation period). We identify five distinct spin-orbit evolution paths and outcomes, only two of which are capable of producing retrograde orbits. We show that these paths to misalignment and the outcomes depend only on two dimensionless parameters, which compare the stellar spin precession frequency with the rate of change of the planet's orbital axis, and the Lidov-Kozai oscillation frequency. Our analysis reveals a number of novel phenomena for the stellar spin evolution, ranging from bifurcation, adiabatic advection, to fully chaotic evolution of spin-orbit angles.

Oscillations between photons and axion-like particles (ALP) travelling in intergalactic magnetic fields have been invoked to explain a number of astrophysical phenomena, or used to constrain ALP properties using observations. One example is the anomalous transparency of the universe to TeV gamma-rays. The intergalactic magnetic field is usually modeled as patches of coherent domains, each with a uniform magnetic field, but the field orientation changes randomly from one domain to the next ("discrete-$\varphi$ model"). We show in this paper that in more realistic situations, when the magnetic field direction varies continuously along the propagation path, the photon-to-ALP conversion probability $P$ can be significantly different from the discrete-$\varphi$ model. In particular, $P$ has a distinct dependence on the photon energy and ALP mass, and can be as large as 100 percent. This result may affect previous constraints on ALP properties based on ALP-photon propagation in intergalactic magnetic fields.

We present numerical simulations of circumbinary accretion onto eccentric and circular binaries using the moving-mesh code AREPO. This is the first set of simulations to tackle the problem of binary accretion using a finite-volume scheme on a freely moving mesh, which allows for accurate measurements of accretion onto individual stars for arbitrary binary eccentricity. While accretion onto a circular binary shows bursts with period of ~5 times the binary period P_b,accretion onto an eccentric binary is predominantly modulated at the period ~1P_b. For an equal-mass circular binary, the accretion rates onto individual stars are quite similar to each other, following the same variable pattern in time. By contrast, for eccentric binaries, one of the stars can accrete at a rate 10-20 times larger than its companion. This "symmetry breaking" between the stars, however, alternates over timescales of order 200 P_b, and can be attributed to a slowly precessing, eccentric circumbinary disk. Over longer timescales, the net accretion rates onto individual stars are the same, reaching a quasi-steady state with the circumbinary disk. These results have important implications for the accretion behavior of binary T-Tauri stars and supermassive binary black holes.

Planets around binary stars and those in multiplanet systems may experience resonant eccentricity excitation and disruption due to perturbations from a distant stellar companion. This "evection resonance" occurs when the apsidal precession frequency of the planet, driven by the quadrupole associated with the inner binary or the other planets, matches the orbital frequency of the external companion. We develop an analytic theory to study the effects of evection resonance on circumbinary planets and multiplanet systems. We derive the general conditions for effective eccentricity excitation or resonance capture of the planet as the system undergoes long-term evolution. Applying to circumbinary planets, we show that inward planet migration may lead to eccentricity growth due to evection resonance with an external perturber, and planets around shrinking binaries may not survive the resonant eccentricity growth. On the other hand, significant eccentricity excitation in multiplanet systems occurs in limited parameter space of planet and binary semimajor axes, and requires the planetary migration to be sufficiently slow.

Lidov-Kozai oscillations of planets in stellar binaries, combined with tidal dissipation, can lead to the formation of hot Jupiters (HJs) or tidal disruption of planets. Recent population synthesis studies have found that the fraction of systems resulting in HJs (F_HJ) depends strongly on the planet mass, host stellar type and tidal dissipation strength, while the total migration fraction F_mig = F_ HJ + F_dis (including both HJ formation and tidal disruption) exhibits much weaker dependence. We present an analytical method for calculating F_HJ and F_mig in the Lidov-Kozai migration scenario. The key ingredient of our method is to determine the critical initial planet-binary inclination angle that drives the planet to reach sufficiently large eccentricity for efficient tidal dissipation or disruption. This calculation includes the effects of the octupole potential and short-range forces on the planet. Our analytical method reproduces the resulting planet migration/disruption fractions from population synthesis, and can be easily implemented for various planet, stellar/companion types, and for different distributions of initial planetary semi-major axes, binary separations and eccentricities. We extend our calculations to planets in the super-Earth mass range and discuss the conditions for such planets to survive Lidov-Kozai migration and form close-in rocky planets.

Observed hot Jupiter (HJ) systems exhibit a wide range of stellar spin-orbit misalignment angles. The origin of these HJs remains unclear. This paper investigates the inward migration of giant planets due to Lidov-Kozai (LK) oscillations induced by a distant (100-1000 AU) stellar companion. We conduct a large population synthesis study, including the octupole gravitational potential from the stellar companion, mutual precession of the host stellar spin axis and planet orbital axis, tidal dissipation in the planet, and stellar spin-down in the host star due to magnetic braking. We consider a range of planet masses ($0.3-5\,M_J$) and initial semi-major axes ($1-5$AU), different properties for the host star, and varying tidal dissipation strengths. The fraction of systems that result in HJs depends on planet mass and stellar type, with $f_{\rm HJ} = 1-4\%$ (depending on tidal dissipation strength) for $M_p=1\,M_J$, and larger (up to $8\%$) for more massive planets. The production efficiency of "hot Saturns" ($M_p=0.3M_J$) is much lower, because most migrating planets are tidally disrupted. We find that the fraction of systems that result in either HJ formation or tidal disruption, $f_{\rm mig} \simeq 11-14\%$ is roughly constant, having little variation with planet mass, stellar type and tidal dissipation strength. This "universal" migration fraction can be understood qualitatively from analytical migration criteria based on the properties of octupole LK oscillations. The distribution of final HJ stellar obliquities exhibits a complex dependence on the planet mass and stellar type. For $M_p = (1-3)M_J$, the distribution is always bimodal, with peaks around $30^{\circ}$ and $130^{\circ}$. The distribution for $5M_J$ planets depends on host stellar type, with a preference for low obliquities for solar-type stars, and higher obliquities for more massive ($1.4M_{\odot}$) stars.

The physical mechanism of angular momentum transport in poorly ionized regions of protoplanetary discs, the dead zones (DZs), is not understood. The presence of a DZ naturally leads to conditions susceptible to the Rossby wave instability (RWI), which produces vortices and spiral density waves that may revive the DZ and be responsible for observed large-scale disc structures. We present a series of two-dimensional hydrodynamic simulations to investigate the role of the RWI in DZs, including its impact on the long-term evolution of the disc and its morphology. The nonlinear RWI can generate Reynolds stresses (effective $\alpha$ parameter) as large as $0.01 - 0.05$ in the DZ, helping to sustain quasi-steady accretion throughout the disc. It also produces novel disc morphologies, including azimuthal asymmetries with $m = 1, 2$, and atypical vortex shapes. The angular momentum transport strength and morphology are most sensitive to two parameters: the radial extent of the DZ and the disc viscosity. The largest Reynolds stresses are produced when the radial extent of the DZ is less than its distance to the central star. Such narrow DZs lead to a single vortex or two coherent antipodal vortices in the quasi-steady state. The edges of wider DZs evolve separately, resulting in two independent vortices and reduced angular momentum transport efficiency. In either case, we find that, because of the Reynolds stresses generated by the nonlinear RWI, gravitational instability is unlikely to play a role in angular momentum transport across the DZ, unless the accretion rate is sufficiently high.

PSR J1357$-$6429 is a young and energetic radio pulsar detected in X-rays and $\gamma$-rays. It powers a compact pulsar wind nebula with a jet visible in X-rays and a large scale plerion detected in X-ray and TeV ranges. Previous multiwavelength studies suggested that the pulsar has a significant proper motion of about 180 mas yr$^{-1}$ implying an extremely high transverse velocity of about 2000 km s$^{-1}$. In order to verify that, we performed radio-interferometric observations of PSR J1357$-$6429 with the the Australia Telescope Compact Array (ATCA) in the 2.1 GHz band. We detected the pulsar with a mean flux density of $212\pm5$ $\mu$Jy and obtained the most accurate pulsar position, RA = 13:57:02.525(14) and Dec = $-$64:29:29.89(15). Using the new and archival ATCA data, we did not find any proper motion and estimated its 90 per cent upper limit $\mu < 106$ mas yr$^{-1}$. The pulsar shows a highly polarised single pulse, as it was earlier observed at 1.4 GHz. Spectral analysis revealed a shallow spectral index $\alpha_{\nu}$ = $0.5 \pm 0.1$. Based on our new radio position of the pulsar, we disclaim its optical counterpart candidate reported before.

Recent observations have shown that circumstellar and circumbinary discs in young stellar binaries are often misaligned with respect to the binary orbital plane. We analyze the tidal truncation of such misaligned discs due to torques applied to the disc at the Lindblad resonances from the tidal forcings of the binary. We consider eccentric binaries with arbitrary binary-disc inclination angles. We determine the dependence of the tidal forcing strengths on the binary parameters and show that they are complicated non-monotonic functions of eccentricity and inclination. We adopt a truncation criterion determined by the balance between resonant torque and viscous torque, and use it to calculate the outer radii of circumstellar discs and the inner radii of circumbinary discs. Misaligned circumstellar discs have systematically larger outer radii than aligned discs, and are likely to fill their Roche lobes if inclined by more than $45^\circ - 90^\circ$, depending on the binary mass ratio and disc viscosity parameter. Misaligned circumbinary discs generally have smaller inner radii than aligned discs, but the details depend sensitively on the binary and disc parameters.

The discovery of transiting circumbinary planets by the Kepler mission suggests that planets can form efficiently around binary stars. None of the stellar binaries currently known to host planets has a period shorter than 7 days, despite the large number of eclipsing binaries found in the Kepler target list with periods shorter than a few days. These compact binaries are believed to have evolved from wider orbits into their current configurations via the so-called Lidov-Kozai migration mechanism, in which gravitational perturbations from a distant tertiary companion induce large-amplitude eccentricity oscillations in the binary, followed by orbital decay and circularization due to tidal dissipation in the stars. Here we explore the orbital evolution of planets around binaries undergoing orbital decay by this mechanism. We show that planets may survive and become misaligned from their host binary, or may develop erratic behavior in eccentricity, resulting in their consumption by the stars or ejection from the system as the binary decays. Our results suggest that circumbinary planets around compact binaries could still exist, and we offer predictions as to what their orbital configurations should be like.

We present analytic calculations of the electromagnetic torques acting on a magnetic neutron star rotating in vacuum, including near-zone torques associated with the inertia of dipole and quadrupole magnetic fields. We incorporate these torques into the rotational dynamics of a rigid-body neutron star, and show that the effects of the inertial torque can be understood as a modification of the moment of inertia tensor of the star. We apply our rotational dynamics equation to the Crab pulsar, including intrinsic distortions of the star and various electromagnetic torques, to investigate the possibility that the counter-alignment of the magnetic inclination angle, as suggested by recent observations, could be explained by pulsar precession. We find that if the effective principal axis of the pulsar is nearly aligned with either the magnetic dipole axis or the rotation axis, then precession may account for the observed counter-alignment over decade timescales. Over the spindown timescale of the pulsar, the magnetic inclination angle always decreases.

We examine the role played by viscosity in the excitation of global oscillation modes (both axisymmetric and non-axisymmetric) in accretion discs around black holes using two-dimensional hydrodynamic simulations. The turbulent viscosity is modeled by the $\alpha$-ansatz, with different equations of state. We consider both discs with transonic radial inflows across the innermost stable circular orbit, and stationary discs truncated by a reflecting wall at their inner edge, representing a magnetosphere. In transonic discs, viscosity can excite several types of global oscillation modes. These modes are either axisymmetric with frequencies close to multiples of the maximum radial epicyclic frequency $\kappa_\mathrm{max}$, non-axisymmetric with frequencies close to multiples of of the innermost stable orbit frequency $\Omega_\mathrm{ISCO}$, or hybrid modes whose frequencies are linear combinations of these two frequencies. Small values of the viscosity parameter $\alpha$ primarily produce non-axisymmetric modes, while axisymmetric modes become dominant for large $\alpha$. The excitation of these modes may be related to an instability of the sonic point, at which the radial infall speed is equal to the sound speed of the gas. In discs with a reflective inner boundary, we explore the effect of viscosity on trapped $p$-modes which are intrinsically overstable due to the corotation resonance effect. The effect of viscosity is either to reduce the growth rates of these modes, or to completely suppress them and excite a new class of higher frequency modes. The latter requires that the dynamic viscosity scales positively with the disc surface density, indicating that it is a result of the classic viscous overstability effect.

We study the effect of large-scale magnetic fields on the non-axisymmetric inertial-acoustic modes (also called p-modes) trapped in the innermost regions of accretion discs around black holes (BHs). These global modes could provide an explanation for the high-frequency quasi-periodic oscillations (HFQPOs) observed in BH X-ray binaries. There may be observational evidence for the presence of such large-scale magnetic fields in the disks since episodic jets are observed in the same spectral state when HFQPOs are detected. We find that a large-scale poloidal magnetic field can enhance the corotational instability and increase the growth rate of the purely hydrodynamic overstable p-modes. In addition, we show that the frequencies of these overstable p-modes could be further reduced by such magnetic fields, making them agree better with observations.

We present analytical expressions for the tidal Love numbers of a giant planet with a solid core and a fluid envelope. We model the core as a uniform, incompressible, elastic solid, and the envelope as a non-viscous fluid satisfying the $n=1$ polytropic equation of state. We discuss how the Love numbers depend on the size, density, and shear modulus of the core. We then model the core as a viscoelastic Maxwell solid and compute the tidal dissipation rate in the planet as characterized by the imaginary part of the Love number $k_2$. Our results improve upon existing calculations based on planetary models with a solid core and a uniform ($n=0$) envelope. Our analytical expressions for the Love numbers can be applied to study tidal distortion and viscoelastic dissipation of giant planets with solid cores of various rheological properties, and our general method can be extended to study tidal distortion/dissipation of super-earths.

The Lidov-Kozai (LK) mechanism plays an important role in the secular evolution of many hierarchical triple systems. The standard LK mechanism consists of large-amplitude oscillations in eccentricity and inclination of a binary subject to the quadrupole potential from an outer perturber. Recent work has shown that when the octupole terms are included in the potential, the inner binary can reach more extreme eccentricities as well as undergo orientation flips. It is known that pericenter precessions due to short-range effects, such as General Relativity and tidal and rotational distortions, can limit the growth of eccentricity and even suppress standard (quadrupolar) LK oscillations, but their effect on the octupole-level LK mechanism has not been fully explored. In this paper, we systematically study how these short-range forces affect the extreme orbital behaviour found in octupole LK cycles. In general, the influence of the octupole potential is confined to a range of initial mutual inclinations Itot centered around 90deg (when the inner binary mass ratio is <<1), with this range expanding with increasing octupole strength. We find that, while the short-range forces do not change the width and location of this "window of influence", they impose a strict upper limit on the maximum achievable eccentricity. This limiting eccentricity can be calculated analytically, and its value holds even for strong octupole potential and for the general case of three comparable masses. Short-range forces also affect orbital flips, progressively reducing the range of Itot within which flips are possible as the intensity of these forces increases.

This paper provides an introduction to a number of astrophysics problems related to strong magnetic fields. The first part deals with issues related to atoms, condensed matter and high-energy processes in very strong magnetic fields, and how these issues influence various aspects of neutron star astrophysics. The second part deals with classical astrophysical effects of magnetic fields: Even relatively "weak" fields can play a strong role in various astrophysical problems, ranging from stars, accretion disks and outflows, to the formation and merger of compact objects.

Many exoplanetary systems containing hot Jupiters are found to possess significant misalignment between the spin axis of the host star and the planet's orbital angular momentum axis. A possible channel for producing such misaligned hot Jupiters involves Lidov-Kozai oscillations of the planet's orbital eccentricity and inclination driven by a distant binary companion. In a recent work (Storch, Anderson & Lai 2014), we have shown that a proto-hot Jupiter undergoing Lidov-Kozai oscillations can induce complex, and often chaotic, evolution of the spin axis of its host star. Here we explore the origin of the chaotic spin behavior and its various features in an idealized non-dissipative system where the secular oscillations of the planet's orbit are strictly periodic. Using Hamiltonian perturbation theory, we identify a set of secular spin-orbit resonances in the system, and show that overlaps of these resonances are responsible for the onset of wide-spread chaos in the evolution of stellar spin axis. The degree of chaos in the system depends on the adiabaticity parameter $\epsilon$, proportional to the ratio of the Lidov-Kozai nodal precession rate and the stellar spin precession rate, and thus depends on the planet mass, semi-major axis and the stellar rotation rate. For systems with zero initial spin-orbit misalignment, our theory explains the occurrence (as a function of $\epsilon$) of "periodic islands" in the middle of a "chaotic ocean" of spin evolution, and the occurrence of restricted chaos in middle of regular/periodic spin evolution. Finally, we discuss a novel "adiabatic resonance advection" phenomenon, in which the spin-orbit misalignment, trapped in a resonance, gradually evolves as the adiabaticity parameter slowly changes. This phenomenon can occur for certain parameter regimes when tidal decay of the planetary orbit is included.

Many exoplanetary systems containing hot Jupiters are observed to have highly misaligned orbital axes relative to the stellar spin axes. Kozai-Lidov oscillations of orbital eccentricity/inclination induced by a binary companion, in conjunction with tidal dissipation, is a major channel for the production of hot Jupiters. We demonstrate that gravitational interaction between the planet and its oblate host star can lead to chaotic evolution of the stellar spin axis during Kozai cycles. As parameters such as the planet mass and stellar rotation period vary, periodic islands can appear in an ocean of chaos, in a manner reminiscent of other dynamical systems. In the presence of tidal dissipation, the complex spin evolution can leave an imprint on the final spin-orbit misalignment angles.

Warped accretion discs are expected in many protostellar binary systems. In this paper, we study the long-term evolution of disc warp and precession for discs with dimensionless thickness $H/r$ larger than their viscosity parameter $\alpha$, such that bending waves can propagate and dominate the warp evolution. For small warps, these discs undergo approximately rigid-body precession. We derive analytical expressions for the warp/twist profiles of the disc and the alignment timescale for a variety of models. Applying our results to circumbinary discs, we find that these discs align with the orbital plane of the binary on a timescale comparable to the global precession time of the disc, and typically much smaller than its viscous timescale. We discuss the implications of our finding for the observations of misaligned circumbinary discs (such as KH 15D) and circumbinary planetary systems (such as Kepler-413); these observed misalignments provide useful constraints on the uncertain aspects of the disc warp theory. On the other hand, we find that circumstellar discs can maintain large misalignments with respect to the plane of the binary companion over their entire lifetime. We estimate that inclination angles larger than $\sim 20^\circ$ can be maintained for typical disc parameters. Overall, our results suggest that while highly misaligned circumstellar discs in binaries are expected to be common, such misalignments should be rare for circumbinary discs. These expectations are consistent with current observations of protoplanetary discs and exoplanets in binaries, and can be tested with future observations.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

Tidal interactions play an important role in the evolution and ultimate fate of compact white dwarf (WD) binaries. Not only do tides affect the pre-merger state (such as temperature and rotation rate) of the WDs, but they may also determine which systems merge and which undergo stable mass transfer. In this paper, we attempt to quantify the effects of rotation on tidal angular momentum transport in binary stars, with specific calculations applied to WD stellar models. We incorporate the effect of rotation using the traditional approximation, in which the dynamically excited gravity waves within the WDs are transformed into gravito-inertial Hough waves. The Coriolis force has only a minor effect on prograde gravity waves, and previous results predicting the tidal spin-up and heating of inspiraling WDs are not significantly modified. However, rotation strongly alters retrograde gravity waves and inertial waves, with important consequences for the tidal spin-down of accreting WDs. We identify new dynamical tidal forcing terms that arise from a proper separation of the equilibrium and dynamical tide components; these new forcing terms are very important for systems near synchronous rotation. Additionally, we discuss the impact of Stokes drift currents on the wave angular momentum flux. Finally, we speculate on how tidal interactions will affect super-synchronously rotating WDs in accreting systems.

A number of radio pulsars exhibit intriguing mode-switching behavior. Recent observations of PSR B0943+10 revealed correlated radio and X-ray mode switches, providing a new avenue for understanding this class of objects. The large X-ray pulse fraction observed during the radio quiet phase (Q mode) was previously interpreted as a result of changing obscuration of X-rays by dense magnetosphere plasma. We show that the large X-ray pulse fraction can be explained by including the beaming effect of a magnetic atmosphere, while remaining consistent with the dipole field geometry constrained by radio observations. We also explore a more extreme magnetic field configuration, where a magnetic dipole displaced from the center of the star produces two magnetic polar caps of different sizes and magnetic field strengths. These models are currently consistent with data in radio and X-rays and can be tested or constrained by future X-ray observations.

Disk accretion onto magnetic stars occurs in a variety of systems, including accreting neutron stars (with both high and low magnetic fields), white dwarfs, and protostars. We review some of the key physical processes in magnetosphere-disk interaction, highlighting the theoretical uncertainties. We also discuss some applications to the observations of accreting neutron stars and protostellar systems, as well as possible connections to protoplanetary disks and exoplanets.

We study the interactions between a protostar and its circumstellar disc under the influence of an external binary companion to determine the evolution of the mutual stellar spin - disc misalignment angle. The gravitational torque on the disc from an inclined binary makes the disc precess around the binary axis, while the star-disc interaction torque due to the rotation-induced quadrupole makes the stellar spin and the disc angular momentum axes precess around each other. A significant star-disc misalignment angle can be generated from a small initial value as the star-disc system evolves in time such that the two precession frequencies cross each other. This ``secular resonance'' behavior can be understood in a geometric way from the precession dynamics of spin and disc angular momenta. We derive the conditions for such resonance to occur, and find that they can be satisfied for reasonable protostar-disc-binary parameters. The evolution of star-disc inclination is also affected by mass accretion and by magnetic star-disc interaction torques, which can either promote or reduce star-disc misalignment. In general, as long as the initial binary-disc inclination is greater than a few degrees, a variety of star-disc misalignment angles can be generated within the disc lifetimes. We discuss the implications of our results for the stellar spin orientations in binaries, for the alignments/misalignments of protostellar discs and debris discs, and for the stellar obliquities in exoplanetary systems. In particular, if hot Jupiters are produced by the Kozai effect induced by an external stellar companion, then it is likely that ``primordial'' star-disc misalignments are already generated by the star-disc-binary interactions. Even for systems where the Kozai effect is suppressed, misaligned planets may still be produced during the protoplanetary disc phase.

We study the possibility of tidal dissipation in the solid cores of giant planets and its implication for the formation of hot Jupiters through high-eccentricity migration. We present a general framework by which the tidal evolution of planetary systems can be computed for any form of tidal dissipation, characterized by the imaginary part of the complex tidal Love number, ${\rm Im}[{\tilde k}_2(\omega)]$, as a function of the forcing frequency $\omega$. Using the simplest viscoelastic dissipation model (the Maxwell model) for the rocky core and including the effect of a nondissipative fluid envelope, we show that with reasonable (but uncertain) physical parameters for the core (size, viscosity and shear modulus), tidal dissipation in the core can accommodate the tidal-Q constraint of the Solar system gas giants and at the same time allows exoplanetary hot Jupiters to form via tidal circularization in the high-e migration scenario. By contrast, the often-used weak friction theory of equilibrium tide would lead to a discrepancy between the Solar system constraint and the amount of dissipation necessary for high-e migration. We also show that tidal heating in the rocky core can lead to modest radius inflation of the planets, particularly when the planets are in the high-eccentricity phase ($e\sim 0.6$) during their high-e migration. Finally, as an interesting by-product of our study, we note that for a generic tidal response function ${\rm Im}[{\tilde k}_2(\omega)]$, it is possible that spin equilibrium (zero torque) can be achieved for multiple spin frequencies (at a given $e$), and the actual pseudo-synchronized spin rate depends on the evolutionary history of the system.

X-ray polarimetry, sometimes alone, and sometimes coupled to spectral and temporal variability measurements and to imaging, allows a wealth of physical phenomena in astrophysics to be studied. X-ray polarimetry investigates the acceleration process, for example, including those typical of magnetic reconnection in solar flares, but also emission in the strong magnetic fields of neutron stars and white dwarfs. It detects scattering in asymmetric structures such as accretion disks and columns, and in the so-called molecular torus and ionization cones. In addition, it allows fundamental physics in regimes of gravity and of magnetic field intensity not accessible to experiments on the Earth to be probed. Finally, models that describe fundamental interactions (e.g. quantum gravity and the extension of the Standard Model) can be tested. We describe in this paper the X-ray Imaging Polarimetry Explorer (XIPE), proposed in June 2012 to the first ESA call for a small mission with a launch in 2017 but not selected. XIPE is composed of two out of the three existing JET-X telescopes with two Gas Pixel Detectors (GPD) filled with a He-DME mixture at their focus and two additional GPDs filled with pressurized Ar-DME facing the sun. The Minimum Detectable Polarization is 14 % at 1 mCrab in 10E5 s (2-10 keV) and 0.6 % for an X10 class flare. The Half Energy Width, measured at PANTER X-ray test facility (MPE, Germany) with JET-X optics is 24 arcsec. XIPE takes advantage of a low-earth equatorial orbit with Malindi as down-link station and of a Mission Operation Center (MOC) at INPE (Brazil).

Recent observations have revealed evidence for the global oscillations of Jupiter and Saturn, which can potentially provide a new window into the interior structure of giant planets. Motivated by these observations, we study the non-radial oscillation modes of giant planets containing a solid core. Our calculations include the elastic response of the core and consider a wide range of possible values of the core shear modulus. While the elasticity of the core only slightly changes the frequencies of acoustic modes, which reside mostly in the fluid envelope, it adds two new classes of shear modes that are largely confined to the core. We also calculate the effects of the Coriolis force on the planetary oscillation modes. In addition to changing the mode frequencies, the Coriolis force can cause the shear modes to mix with the f-modes. Such mixing occurs when the frequencies of the shear mode and the f-mode are close to each other, and results in "mixed modes" with similar properties that are slightly split in frequency. We discuss our results in light of the recent work by Hedman & Nicholson (2013), which revealed the presence of density waves in Saturn's C-ring that appear to be excited by the gravitational perturbations associated with the f-mode oscillations within Saturn. We find that the fine splitting in wave frequencies observed in the rings can in principle be explained by the rotation-induced mixing between core shear modes and f-modes, possibly indicating the presence of a solid core within Saturn. However, in our current calculations, which assume rigid-body rotation and include only first-order rotational effects, significant fine-tuning in the planetary model parameters is needed in order to achieve these mode mixings and to explain the observed fine frequency splitting. We briefly discuss other effects that may modify the f-modes and facilitate mode mixing.

We study the dynamics of spiral waves and oscillation modes in relativistic rotating discs around black holes. Generalizing the Newtonian theory, we show that wave absorption can take place at the corotation resonance, where the pattern frequency of the wave matches the background disc rotation rate. We derive the general relativistic expression for the disc vortensity (vorticity divided by surface density), which governs the behaviour of density perturbation near corotation. Depending on the gradient of the generalized disc vortensity, corotational wave absorption can lead to the amplification or damping of the spiral wave. We apply our general theory of relativistic wave dynamics to calculate the non-axisymmetric inertial-acoustic modes (also called p-modes) trapped in the inner-most region of a black hole accretion disc. Because general relativity changes the profiles of the radial epicyclic frequency and disc vortensity near the inner disc edge close to the black hole, these p-modes can become overstable under appropriate conditions. We present the numerical results of the frequencies and growth rates of p-modes for various black hole spin and model disc parameters (the surface density profile and sound speed), and discuss their implications for understanding the enigmatic high-frequency quasi-periodic oscillations observed in black hole X-ray binaries.

The Kepler satellite has discovered a number of transiting planets around close binary stars. These circumbinary systems have highly aligned planetary and binary orbits. In this paper, we explore how the mutual inclination between the planetary and binary orbits may reflect the physical conditions of the assembly of protoplanetary disks and the interaction between protostellar binaries and circumbinary disks. Given the turbulent nature of star-forming molecular clouds, it is possible that the gas falling onto the outer region of a circumbinary disk and the central protostellar binary have different axes of rotation. Thus, the newly assembled circumbinary disk can be misaligned with respect to the binary. However, the gravitational torque from the binary produces a warp and twist in the disk, and the back-reaction torque tends to align the disk and the binary orbital plane. We present a new, analytic calculation of this alignment torque, and show that the binary-disk inclination angle can be reduced appreciably after the binary accretes a few percent of its mass from the disk. Our calculation suggests that in the absence of other disturbances, circumbinary disks and planets around close (sub-AU) stellar binaries, for which mass accretion onto the proto-binary is very likely to have occurred, are expected to be highly aligned with the binary orbits, while disks and planets around wide binaries can be misaligned. Measurements of the mutual inclinations of circumbinary planetary systems can provide a clue to the birth environments of such systems.

In this document, we describe the scientific potential of blazar observations with a X-ray polarimetry mission like GEMS (Gravity and Extreme Magnetism SMEX). We describe five blazar science investigations that such a mission would enable: (i) the structure and the role of magnetic fields in AGN jets, (ii) analysis of the polarization of the synchrotron X-ray emission from AGN jets, (iii) discrimination between synchrotron self-Compton and external Compton models for blazars with inverse Compton emission in the X-ray band, (iv) a precision study of the polarization properties of the X-ray emission from Cen-A, (v) tests of Lorentz Invariance based on X-ray polarimetric observations of blazars. We conclude with a discussion of a straw man observation program and recommended accompanying multiwavelength observations.

We examine the expected X-ray polarization properties of neutron-star X-ray sources of various types, e.g., accretion and rotation powered pulsars, magnetars, and low-mass X-ray binaries. We summarize the model calculations leading to these expected properties. We describe how a comparison of these with their observed properties, as inferred from GEMS data, will probe the essential dynamical, electromagnetic, plasma, and emission processes in neutron-star binaries, discriminate between models of these processes, and constrain model parameters. An exciting goal is the first observational demonstration in this context of the existence of vacuum resonance, a fundamental quantum electrodynamical phenomenon first described in the 1930s.

We present here a summary of the scientific goals behind the Gravity and Extreme Magnetism SMEX (GEMS) X-ray polarimetry mission's black hole (BH) observing program. The primary targets can be divided into two classes: stellar-mass galactic BHs in accreting binaries, and super-massive BHs in the centers of active galactic nuclei (AGN). The stellar-mass BHs can in turn be divided into various X-ray spectral states: thermal-dominant (disk), hard (radio jet), and steep power-law (hot corona). These different spectral states are thought to be generated by different accretion geometries and emission mechanisms. X-ray polarization is an ideal tool for probing the geometry around these BHs and revealing the specific properties of the accreting gas.

The Rossby wave instability, associated with density bumps in differentially rotating discs, may arise in several different astrophysical contexts, such as galactic or protoplanetary discs. While the linear phase of the instability has been well studied, the nonlinear evolution and especially the saturation phase remain poorly understood. In this paper, we test the non-linear saturation mechanism analogous to that derived for wave-particle interaction in plasma physics. To this end we perform global numerical simulations of the evolution of the instability in a two-dimensional disc. We confirm the physical mechanism for the instability saturation and show that the maximum amplitude of vorticity can be estimated as twice the linear growth rate of the instability. We provide an empirical fitting formula for this growth rate for various parameters of the density bump. We also investigate the effects of the azimuthal mode number of the instability and the energy leakage in the spiral density waves. Finally, we show that our results can be extrapolated to 3D discs.

Merging neutron star (NS) binaries may be detected by ground-based gravitational wave (GW) interferometers (e.g. LIGO/VIRGO) within this decade and may also generate electromagnetic radiation detectable by wide-field, fast imaging telescopes that are coming online. The GWs can provide new constraint on the NS equation of state (including mass-radius relation and the related nuclear symmetry energy through resonant g-modes). This paper reviews various hydrodynamical (including equilibrium and dynamical/resonant tides) and electrodynamical processes in coalescing NS binaries, with focus on the pre-merger phase.

The physical origin of high-frequency QPOs (HFQPOs) in black-hole X-ray binaries remains an enigma despite many years of detailed observational studies. Although there exists a number of models for HFQPOs, many of these are simply "notions" or "concepts" without actual calculation derived from fluid or disk physics. Future progress requires a combination of numerical simulations and semi-analytic studies to extract physical insights. We review recent works on global oscillation modes in black-hole accretion disks, and explain how, with the help of general relativistic effects, the energy stored in the disk differential rotation can be pumped into global spiral density modes in the disk, making these modes grow to large amplitudes under certain conditions ("corotational instability"). These modes are robust in the presence of disk magnetic fields and turbulence. The computed oscillation mode frequencies are largely consistent with the observed values for HFQPOs in BH X-ray binaries. The approximate 2:3 frequency ratio is also expected from this model. The connection of HFQPOs with other disk properties (such as production of episodic jets) is also discussed.

Tidal dissipation in compact white dwarf (WD) binary systems significantly influences the physical conditions (such as surface temperature and rotation rate) of the WDs prior to mass transfer or merger. In these systems, the dominant tidal effects involve the excitation of gravity waves and their dissipation in the outer envelope of the star. We calculate the amplitude of tidally excited gravity waves in low-mass (0.3M_\odot) helium-core (He) WDs as a function of the tidal forcing frequency \omega. Like carbon-oxygen (CO) WDs studied in our previous paper, we find that the dimensionless tidal torque F(\omega) (inversely proportional to the effective tidal quality factor) has an erratic dependence on \omega. On average, F(\omega) scales approximately as \omega^6, and is several orders of magnitude smaller for He WDs than for CO WDs. We find that tidal torques can begin to synchronize the WD rotation when the orbital period is less than about a hour, although a nearly constant asynchronization is maintained even at small periods. We examine where the tidally excited gravity waves experience non-linear breaking or resonant absorption at a critical layer, allowing us to estimate the location and magnitude of tidal heating in the WD envelope. We then incorporate tidal heating in the MESA stellar evolution code, calculating the physical conditions of the WD as a function of orbital period for different WD models. We find that tidal heating makes a significant contribution to the WD luminosity for short-period (~10 min) systems such as SDSS J0651+2844. We also find that for WDs containing a hydrogen envelope, tidal heating can trigger runaway hydrogen shell burning, leading to a nova-like event before the onset of mass transfer.

We present two-dimensional inviscid hydrodynamic simulations of overstable inertial-acoustic oscillation modes (p-modes) in black-hole accretion discs. These global spiral waves are trapped in the inner-most region of the disc, and are driven overstable by wave absorption at the corotation resonance ($r_c$) when the gradient of the background disc vortensity (vorticity divided by surface density) at $r_c$ is positive and the disc inner boundary is sufficiently reflective. Previous linear calculations have shown that the growth rates of these modes can be as high as 10% of the rotation frequency at the disc inner edge. We confirm these linear growth rates and the primary disc oscillation frequencies in our simulations when the mode amplitude undergoes exponential growth. We show that the mode growth saturates when the radial velocity perturbation becomes comparable to the disc sound speed. During the saturation stage, the primary disc oscillation frequency differs only slightly (by less than a few percent) from the linear mode frequency. Sharp features in the fluid velocity profiles at this stage suggest that the saturation results from nonlinear wave steepening and mode-mode interactions.

We study the effect of large-scale magnetic fields on the non-axisymmetric Rossby wave instability (RWI) in accretion discs. The instability develops around a density bump, which is likely present in the transition region between the active zone and dead zone of protoplanetary discs. Previous works suggest that the vortices resulting from the RWI may facilitate planetesimal formation and angular momentum transport. We consider discs threaded by a large-scale poloidal magnetic field, with a radial field component at the disc surface. Such field configurations may lead to the production of magnetic winds or jets. In general, the magnetic field can affect the RWI even when it is sub-thermal (plasma $\beta\sim 10$). For infinitely thin discs, the instability can be enhanced by about 10 percent. For discs with finite thickness, with a radial gradient of the magnetic field strength, the RWI growth rate can increase significantly (by a factor of $\sim 2$) as the field approaches equipartition ($\beta \sim 1$). Our result suggests that the RWI can continue to operate in discs that produce magnetic winds.

We study the tidal excitation of gravity modes (g-modes) in compact white dwarf binary systems with periods ranging from minutes to hours. As the orbit of the system decays via gravitational radiation, the orbital frequency increases and sweeps through a series of resonances with the g-modes of the white dwarf. At each resonance, the tidal force excites the g-mode to a relatively large amplitude, transferring the orbital energy to the stellar oscillation. We calculate the eigenfrequencies of g-modes and their coupling coefficients with the tidal field for realistic non-rotating white dwarf models. Using these mode properties, we numerically compute the excited mode amplitude in the linear approximation as the orbit passes though the resonance, including the backreaction of the mode on the orbit. We also derive analytical estimates for the mode amplitude and the duration of the resonance, which accurately reproduce our numerical results for most binary parameters. We find that the g-modes can be excited to a dimensionless (mass-weighted) amplitude up to 0.1, with the mode energy approaching $10^{-3}$ of the gravitational binding energy of the star. This suggests that thousands of years prior to the binary merger, the white dwarf may be heated up significantly by tidal interactions. However, more study is needed since the physical amplitudes of the excited oscillation modes become highly nonlinear in the outer layer of the star, which can reduce the mode amplitude attained by tidal excitation.

Compact binary white dwarfs (WDs) undergoing orbital decay due to gravitational radiation can experience significant tidal heating prior to merger. In these WDs, the dominant tidal effect involves the excitation of outgoing gravity waves in the inner stellar envelope and the dissipation of these waves in the outer envelope. As the binary orbit decays, the WDs are synchronized from outside in (with the envelope synchronized first, followed by the core). We examine the deposition of tidal heat in the envelope of a Carbon-Oxygen WD and study how such tidal heating affects the structure and evolution of the WD. We show that significant tidal heating can occur in the star's degenerate hydrogen layer. This layer heats up faster than it cools, triggering runaway nuclear fusion. Such "tidal novae" may occur in all WD binaries containing a CO WD, at orbital periods between 5 min and 20 min, and precede the final merger by 10^5-10^6 years.

In compact white dwarf (WD) binary systems (with periods ranging from minutes to hours), dynamical tides involving the excitation and dissipation of gravity waves play a dominant role in determining the physical conditions of the WDs prior to mass transfer or binary merger. We calculate the amplitude of the tidally excited gravity waves as a function of the tidal forcing frequency \omega=2(\Omega-\Omega_s) (where \Omega is the orbital frequency and \Omega_s is the spin frequency) for several realistic carbon-oxygen WD models, assuming that the waves are efficiently dissipated in the outer layer of the star by nonlinear effects or radiative damping. The mechanism of wave excitation in WDs is complex due to the sharp features associated with composition changes inside the WD, and in our WD models gravity waves are launched just below the helium-carbon boundary. We find that the tidal torque on the WD and the related tidal energy transfer rate, \dot E_{\rm tide}, depend on \omega in an erratic way. On average, \dot E_{\rm tide} scales approximately as \Omega^5\omega^5 for a large range of tidal frequencies. We also study the effects of dynamical tides on the long-term evolution of WD binaries. Above a critical orbital frequency \Omega_c, corresponding to an orbital period of order one hour (depending on WD models), dynamical tides efficiently drive \Omega_s toward \Omega, although a small, almost constant degree of asynchronization (\Omega-\Omega_s\sim {\rm constant}) is maintained even at the smallest binary periods. While the orbital decay is always dominated by gravitational radiation, the tidal energy transfer can induce significant phase error in the low-frequency gravitational waveforms, detectable by the planned LISA project. Tidal dissipation may also lead to significant heating of the WD envelope and brightening of the system long before binary merger.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

Observations of hot Jupiters around solar-type stars with very short orbital periods (~day) suggest that tidal dissipation in such stars is not too efficient so that these planets can survive against rapid orbital decay. This is consistent with recent theoretical works, which indicate that the tidal Q of planet-hosting stars can indeed be much larger than the values inferred from stellar binaries. On the other hand, recent measurements of Rossiter-McLaughlin effects in transiting hot Jupiter systems not only reveal that many such systems have misaligned stellar spin with respect to the orbital axis, but also show that systems with cooler host stars tend to have aligned spin and orbital axes. Winn et al. suggested that this obliquity - temperature correlation may be explained by efficient damping of stellar obliquity due to tidal dissipation in the star. This explanation, however, is in apparent contradiction with the survival of these short-period hot Jupiters. We show that in the solar-type parent stars of close-in exoplanetary systems, the effective tidal Q governing the damping of stellar obliquity can be much smaller than that governing orbital decay. This is because for misaligned systems, the tidal potential contains a Fourier component with frequency equal to the stellar spin frequency (in the rotating frame of the star). This component can excite inertial waves in the convective envelope of the star, and the dissipation of inertial waves then leads to a spin-orbit alignment torque, but not orbital decay. By contrast, for aligned systems, such inertial wave excitation is forbidden since the tidal forcing frequency is much larger than the stellar spin frequency. We derive a general effective tidal evolution theory for misaligned binaries, taking account of different tidal responses and dissipation rates for different tidal forcing components.

The unipolar induction DC circuit model, originally developed by Goldreich & Lynden-Bell for the Jupiter-Io system, has been applied to different types of binary systems in recent years. We show that there exists an upper limit to the magnetic interaction torque and energy dissipation rate in such model. This arises because when the resistance of the circuit is too small, the large current flow severely twists the magnetic flux tube connecting the two binary components, leading to breakdown of the circuit. Applying this limit, we find that in coalescing neutron star binaries, magnetic interactions produce negligible correction to the phase evolution of the gravitational waveform, even for magnetar-like field strengths. However, energy dissipation in the binary magnetosphere may still give rise to electromagnetic radiation prior to the final merger. For ultra-compact white dwarf binaries, we find that DC circuit does not provide adequate energy dissipation to explain the observed X-ray luminosities of several sources. For exoplanetary systems containing close-in Jupiters or super-Earths, magnetic torque and dissipation are negligible, except possibly during the early T Tauri phase, when the stellar magnetic field is stronger than 10^3G.

Recent observations of the central compact object in the Kesteven 79 supernova remnant show that this neutron star (NS) has a weak dipole magnetic field (a few x 10^{10} G) but an anomalously large (~ 64%) pulse fraction in its surface X-ray emission. We explore the idea that a substantial sub-surface magnetic field exists in the NS crust, which produces diffuse hot spots on the stellar surface due to anisotropic heat conduction, and gives rise to the observed X-ray pulsation. We develop a general-purpose method, termed "Temperature Template with Full Transport" (TTFT), that computes the synthetic pulse profile of surface X-ray emission from NSs with arbitrary magnetic field and surface temperature distributions, taking into account magnetic atmosphere opacities, beam pattern, vacuum polarization, and gravitational light bending. We show that a crustal toroidal magnetic field of order a few x 10^{14} G or higher, varying smoothly across the crust, can produce sufficiently distinct surface hot spots to generate the observed pulse fraction in the Kes 79 NS. This result suggests that substantial sub-surface magnetic fields, much stronger than the "visible" dipole fields, may be buried in the crusts of some young NSs, and such hidden magnetic fields can play an important role in their observational manifestations. The general TTFT tool we have developed can also be used for studying radiation from other magnetic NSs.

We study global non-axisymmetric oscillation modes and instabilities in magnetosphere- disc systems, as expected in neutron star X-ray binaries and possibly also in accreting black hole systems. Our two-dimensional magnetosphere-disc model consists of a Keplerian disc in contact with an uniformly rotating magnetosphere with low plasma density. Two types of global overstable modes exist in such systems, the interface modes and the disc inertial-acoustic modes. We examine various physical effects and parameters that influence the properties of these oscillation modes, particularly their growth rates, including the magnetosphere field configuration, the velocity and density contrasts across the magnetosphere-disc interface, the rotation profile (with Newtonian or General Relativistic potential), the sound speed and magnetic field of the disc. The interface modes are driven unstable by Rayleigh-Taylor and Kelvin-Helmholtz in- stabilities, but can be stabilized by the toroidal field (through magnetic tension) and disc differential rotation (through finite vorticity). General relativity increases their growth rates by modifying the disc vorticity outside the magnetosphere boundary. The interface modes may also be affected by wave absorption associated with corotation resonance in the disc. In the presence of a magnetosphere, the inertial-acoustic modes are effectively trapped at the innermost region of the relativistic disc just outside the interface. They are driven unstable by wave absorption at the corotation resonance, but can be stabilized by modest disc magnetic fields. The overstable oscillation modes studied in this paper have characteristic properties that make them possible candidates for the quasi-periodic oscillations observed in X-ray binaries.

The Rossby wave instability (RWI) is a promising mechanism for producing large-scale vortices in protoplanetary discs. The instability operates around a density bump in the disc, and the resulting vortices may facilitate planetesimal formation and angular momentum transfer in the disc dead zone. Most previous works on the RWI deal with two-dimensional (height-integrated) discs. However, vortices may have different dynamical behaviours in 3D than in 2D. Recent numerical simulations of the RWI in 3D global discs by Meheut et al. have revealed intriguing vertical structure of the vortices, including appreciable vertical velocities. In this paper we present a linear analysis of the RWI in 3D global models of isothermal discs. We calculate the growth rates of the Rossby modes (of various azimuthal wave numbers m = 2 - 6) trapped around the fiducial density bump and carry out 3D numerical simulations to compare with our linear results. We show that the 3D RWI growth rates are only slightly smaller than the 2D growth rates, and the velocity structures seen in the numerical simulations during the linear phase are in agreement with the velocity eigenfunctions obtained in our linear calculations. This numerical benchmark shows that numerical simulations can accurately describe the instability. The angular momentum transfer rate associated with Rossby vortices is also studied.

Recent observation of the tidally-excited stellar oscillations in the main-sequence binary KOI-54 by the KEPLER satellite provides a unique opportunity for studying dynamical tides in eccentric binary systems. We develop a general theory of tidal excitation of oscillation modes of rotating binary stars, and apply our theory to tidally excited gravity modes (g-modes) in KOI-54. The strongest observed oscillations, which occur at 90 and 91 times the orbital frequency, are likely due to prograde m=2 modes (relative to the stellar spin axis) locked in resonance with the orbit. The remaining flux oscillations with frequencies that are integer multiples of the orbital frequency are likely due to nearly resonant m=0 g-modes; such axisymmetric modes generate larger flux variations compared to the m=2 modes, assuming that the spin inclination angle of the star is comparable to the orbital inclination angle. We examine the process of resonance mode locking under the combined effects of dynamical tides on the stellar spin and orbit and the intrinsic stellar spindown. We show that KOI-54 can naturally evolve into a state in which at least one m=2 mode is locked in resonance with the orbital frequency. Our analysis provides an explanation for the fact that only oscillations with frequencies less than 90-100 times the orbital frequency are observed. We have also found evidence from the published KEPLER result that three-mode nonlinear coupling occurs in the KOI-54 system. We suggest that such nonlinear mode coupling may explain the observed oscillations that are not harmonics of the orbital frequency.

We quantify the systematic impact of the non-spherical shape of transiting planets and brown dwarfs, due to tidal forces and rotation, on the observed transit depth. Such a departure from sphericity leads to a bias in the derivation of the transit radius from the light curve and affects the comparison with planet structure and evolution models which assume spherical symmetry. As the tidally deformed planet projects its smallest cross section area during the transit, the measured effective radius is smaller than the one of the unperturbed spherical planet. This effect can be corrected by calculating the theoretical shape of the observed planet. We derive simple analytical expressions for the ellipsoidal shape of a fluid object (star or planet) accounting for both tidal and rotational deformations and calibratre it with fully numerical evolution models in the 0.3Mjup-75Mjup mass range. Our calculations yield a 20% effect on the transit depth, i.e. a 10% decrease of the measured radius, for the extreme case of a 1Mjup planet orbiting a Sun-like star at 0.01AU. For the closest planets detected so far (< 0.05 AU), the effect on the radius is of the order of 1 to 10%, by no means a negligible effect, enhancing the puzzling problem of the anomalously large bloated planets. These corrections must thus be taken into account for a correct determination of the radius from the transit light curve. Our analytical expressions can be easily used to calculate these corrections, due to the non-spherical shape of the planet, on the observed transit depth and thus to derive the planet's real equilibrium radius. They can also be used to model ellipsoidal variations of the stellar flux now detected in the CoRoT and Kepler light curves. We also derive directly usable analytical expressions for the moment of inertia, oblateness and Love number (k_2) of a fluid planet as a function of its mass.

We present the results of a low-resolution spectral abundance study of 25 stars in the Bootes I dwarf spheroidal (dSph) galaxy. The data were obtained with the LRIS instrument at Keck Observatory, and allow us to measure [Fe/H], [C/Fe], and [alpha/Fe] for each star. We find both a large spread in metallicity (2.1 dex in [Fe/H]) as well as the low average metallicity in this system, <[Fe/H]>=-2.59, matching previous estimates. This sample includes a newly discovered extremely metal-poor star, with [Fe/H]=-3.8, that is one of the most metal-poor stars yet found in a dSph. We compare the metallicity distribution function of Bootes I to analytic chemical evolution models. While the metallicity distribution function of Bootes I is best fit by an Extra Gas chemical evolution model, leaky-box models also provide reasonable fits. We also find that the [alpha/Fe] distribution and the carbon-enhanced metal-poor fraction of our sample (12%) are reasonable matches to Galactic halo star samples in the same metallicity range, indicating that at these low metallicities, systems like the Bootes I ultra-faint dSph could have been contributors to the Galactic halo.

The magnetosphere of a pulsar is composed of relativistic plasmas streaming along the magnetic field lines and corotating with the pulsar. We study the intrinsic Faraday rotation in the pulsar magnetosphere by critically examining the wave modes and the variations of polarization properties for the circularly polarized natural modes under various assumptions about the magnetosphere plasma properties. Since it is difficult to describe analytically the Faraday rotation effect in such a plasma, we use numerical integrations to study the wave propagation effects in the corotating magnetosphere. Faraday rotation effect is identified among other propagation effects, such as wave mode coupling and the cyclotron absorption. In a highly symmetrical electron-positron pair plasma, the Faraday rotation effect is found to be negligible. Only for asymmetrical plasmas, such as the electron-ion streaming plasma, can the Faraday rotation effect become significant, and the Faraday rotation angle is found to be approximately proportional to $\lambda^{0.5}$ instead of the usual $\lambda^2$-law. For such electrons-ion plasma of pulsar magnetosphere, the induced rotation measure becomes larger at higher frequencies, and should have opposite signs for the emissions from opposite magnetic poles.

We present analysis of high-resolution spectra of a sample of stars in the globular cluster M5 (NGC 5904). The sample includes stars from the red giant branch (seven stars), the red horizontal branch (two stars), and the asymptotic giant branch (eight stars), with effective temperatures ranging from 4000 K to 6100 K. Spectra were obtained with the HIRES spectrometer on the Keck I telescope, with a wavelength coverage from 3700 to 7950 angstroms for the HB and AGB sample, and 5300 to 7600 angstroms for the majority of the RGB sample. We find offsets of some abundance ratios between the AGB and the RGB branches. However, these discrepancies appear to be due to analysis effects, and indicate that caution must be exerted when directly comparing abundance ratios between different evolutionary branches. We find the expected signatures of pollution from material enriched in the products of the hot hydrogen burning cycles such as the CNO, Ne-Na, and Mg-Al cycles, but no significant differences within these signatures among the three stellar evolutionary branches especially when considering the analysis offsets. We are also able to measure an assortment of neutron-capture element abundances, from Sr to Th, in the cluster. We find that the neutron-capture signature for all stars is the same, and shows a predominately r-process origin. However, we also see evidence of a small but consistent extra s-process signature that is not tied to the light-element variations, pointing to a pre-enrichment of this material in the protocluster gas.

Magnetic interactions between a protostar and its accretion disc tend to induce warping in the disc and produce secular changes in the stellar spin direction, so that the spin axis may not always be perpendicular to the disc. This may help explain the recently observed spin-orbit misalignment in a number of exoplanetary systems. We study the dynamics of warped protoplanetary discs under the combined effects of magnetic warping/precession torques and internal stresses in the disc, including viscous damping of warps and propagation of bending waves. We show that when the outer disc axis is misaligned with the stellar spin axis, the disc evolves towards a warped steady-state on a timescale that depends on the disc viscosity or the bending wave propagation speed, but in all cases is much shorter than the timescale for the spin evolution (of order of a million years). Moreover, for the most likely physical parameters characterizing magnetic protostars, circumstellar discs and their interactions, the steady-state disc has a rather small warp, such that the whole disc lies approximately in a single plane determined by the outer disc boundary conditions, although more extreme parameters may give rise to larger disc warps. In agreement with our recent analysis (Lai et al. 2010) based on flat discs, we find that the back-reaction magnetic torques of the slightly warped disc on the star can either align the stellar spin axis with the disc axis or push it towards misalignment, depending on the parameters of the star-disc system. This implies that newly formed planetary systems may have a range of inclination angles between the stellar spin axis and the symmetry axis of the planetary orbits.

Recent observations have shown that in many exoplanetary systems the spin axis of the parent star is misaligned with the planet's orbital axis. These have been used to argue against the scenario that short-period planets migrated to their present-day locations due to tidal interactions with their natal discs. However, this interpretation is based on the assumption that the spins of young stars are parallel to the rotation axes of protostellar discs around them. We show that the interaction between a magnetic star and its circumstellar disc can (but not always) have the effect of pushing the stellar spin axis away from the disc angular momentum axis toward the perpendicular state and even the retrograde state. Planets formed in the disc may therefore have their orbital axes misaligned with the stellar spin axis, even before any additional planet-planet scatterings or Kozai interactions take place. In general, magnetosphere--disc interactions lead to a broad distribution of the spin--orbit angles, with some systems aligned and other systems misaligned.

Recent hydrodynamical simulations have shown that differentially rotating neutron stars formed in core-collapse supernovae may develop global non-axisymmetric instabilities even when $T/|W|$ (the ratio of the rotational kinetic energy $T$ to the gravitational potential energy $|W|$) is relatively small (less than 0.1). Such low-$T/|W|$ instability can give rise to efficient gravitational wave emission from the proto-neutron star. We investigate how this instability is affected by magnetic fields using a cylindrical stellar model. Wave absorption at the corotation resonance plays an important role in facilitating the hydrodynamic low-$T/|W|$ instability. In the presence of a toroidal magnetic field, the corotation resonance is split into two magnetic resonances where wave absorptions take place. We show that the toroidal magnetic field suppresses the low-$T/|W|$ instability when the total magnetic energy $W_{\rm B}$ is of order $0.2\,T$ or larger, corresponding to toroidal fields of a few $\times 10^{16}$ G or stronger. Although poloidal magnetic fields do not influence the instability directly, they can affect the instability by generating toroidal fields through linear winding of the initial poloidal field and magneto-rotational instability. We show that an initial poloidal field with strength as small as $10^{14}$ G may suppress the low-$T/|W|$ instability.

We describe a three hour inquiry activity involving converging lenses and telescopes as part of a semester-long astronomy lab course for non-science majors at Hartnell Community College in Salinas, CA. Students were shown several short demonstrations and given the chance to experiment with the materials, after which there was a class discussion about the phenomena they observed. Students worked in groups of 2-4 to design their own experiments to address a particular question of interest to them and then presented their findings to the class. An instructor-led presentation highlighted the students' discoveries and the lab's content goals, followed by a short worksheet-based activity that guided them in applying their new knowledge to build a simple telescope using two converging lenses. The activity was successful in emphasizing communication skills and giving students opportunities to engage in the process of science in different ways. One of the biggest challenges in designing this activity was covering all of the content given the short amount of time available. Future implementations may have more success by splitting the lab into two sessions, one focusing on converging lenses and the other focusing on telescopes.

We present a method for the determination of [alpha/Fe] ratios from low-resolution (R = 2000) SDSS/SEGUE stellar spectra. By means of a star-by-star comparison with degraded spectra from the ELODIE spectral library and with a set of moderately high-resolution (R = 15,000) and medium-resolution (R = 6000) spectra of SDSS/SEGUE stars, we demonstrate that we are able to measure [alpha/Fe] from SDSS/SEGUE spectra (with S/N > 20/1) to a precision of better than 0.1 dex, for stars with atmospheric parameters in the range Teff = [4500, 7000] K, log g = [1.5, 5.0], and [Fe/H] = [-1.4, +0.3], over the range [alpha/Fe] = [-0.1, +0.6]. For stars with [Fe/H] < -1.4, our method requires spectra with slightly higher signal-to-noise to achieve this precision (S/N > 25/1). Over the full temperature range considered, the lowest metallicity star for which a confident estimate of [alpha/Fe] can be obtained from our approach is [Fe/H] ~ -2.5; preliminary tests indicate that a metallicity limit as low as [Fe/H] ~ -3.0 may apply to cooler stars. As a further validation of this approach, weighted averages of [alpha/Fe] obtained for SEGUE spectra of likely member stars of Galactic globular clusters (M15, M13, and M71) and open clusters (NGC 2420, M67, and NGC 6791) exhibit good agreement with the values of [alpha/Fe] from previous studies. The results of the comparison with NGC 6791 imply that the metallicity range for the method may extend to ~ +0.5.

Hot accretion tori around a compact object are known to be susceptible to a global hydrodynamical instability, the so-called Papaloizou-Pringle (PP) instability, arising from the interaction of non-axisymmetric waves across the corotation radius, where the wave pattern speed matches the fluid rotation rate. However, accretion tori produced in various astrophysical situations (e.g., collapsars and neutron star binary mergers) are likely to be highly magnetized. We study the effect of magnetic fields on the PP instability in incompressible tori with various magnetic strengths and structures. In general, toroidal magnetic fields have significant effects on the PP instability: For thin tori (with the fractional width relative to the outer torus radius much less than unity), the instability is suppressed at large field strengths with the corresponding toroidal Alfven speed $v_{A\phi}\go 0.2r\Omega$ (where $\Omega$ is the flow rotation rate). For thicker tori (with the fractional width of order 0.4 or larger), which are hydrodynamically stable, the instability sets in for sufficiently strong magnetic fields (with $v_{A\phi}\go 0.2 r\Omega$). Our results suggest that highly magnetized accretion tori may be subjected to global instability even when it is stable against the usual magneto-rotational instability.

Low-order, non-axisymmetric p-modes (also referred as inertial-acoustic modes) trapped in the inner-most region of hydrodynamic accretion discs around black holes, are plausible candidates for high-frequency quasi-periodic oscillations (QPOs) observed in a number of accreting black-hole systems. These modes are subject to global instabilities due to wave absorption at the corotation resonance (where the wave pattern frequency $\omega/m$ equals the disc rotation rate $\Omega$), when the fluid vortensity, $\zeta=\kappa^2/(2\Omega\Sigma)$ (where $\kappa$ and $\Sigma$ are the radial epicyclic frequency and disc surface density, respectively), has a positive gradient. We investigate the effects of disc magnetic fields on the wave absorption at corotation and the related wave super-reflection of the corotation barrier, and on the overstability of disc p-modes. For discs with a pure toroidal field, the corotation resonance is split into two magnetic resonances, where the wave frequency in the corotating frame of the fluid, $\tomega=\omega-m\Omega$, matches the slow magnetosonic wave frequency. Significant wave energy/angular momentum absorption occurs at both magnetic resonances, but with opposite signs. The combined effect of the two magnetic resonances is to reduce the super-reflection and the growth rate of the overstable p-modes. We show that even a subthermal toroidal field may suppress the overstability of hydrodynamic (B=0) p-modes. For accretion discs with mixed (toroidal and vertical) magnetic fields, two additional Alfven resonances appear, where $\tomega$ matches the local Alfven wave frequency. They further reduce the growth rate of p-modes. Our results suggest that in order for the non-axisymmetric p-modes to be a viable candidate for the observed high-frequency QPOs, the disc magnetic field must be appreciably subthermal, or other mode excitation mechanisms are at work.

We study mass transfer by Roche lobe overflow in close-in exoplanetary systems. The planet's atmospheric gas passes through the inner Lagrangian point and flows along a narrow stream, accelerating to 100-$200\kms$ velocity before forming an accretion disk. We show that the cylinder-shaped accretion stream can have an area (projected in the plane of the sky) comparable to that of the planet and a significant optical depth to spectral line absorption. Such a "transiting cylinder" may produce an earlier ingress of the planet transit, as suggested by recent HST observations of the WASP-12 system. The asymmetric disk produced by the accretion stream may also lead to time-dependent obscuration of the star light and apparent earlier ingress. We also consider the interaction of the stellar wind with the planetary magnetosphere. Since the wind speed is subsonic/sub-Alfvenic and comparable to the orbital velocity of the planet, the head of the magnetopause lies eastward relative to the substellar line (the line joining the planet and the star). The gas around the magnetopause may, if sufficiently compressed, give rise to asymmetric ingress/egress during the planet transit, although more works are needed to evaluate this possibility.

We study the effect of corotation resonance on the inertial-acoustic oscillations (p-modes) of black-hole accretion discs. Previous works have shown that for barotropic flows (where the pressure depends only on the density), wave absorption at the corotation resonance can lead to mode growth when the disc vortensity, $\zeta=\kappa^2/(2\Omega\Sigma)$ (where $\Omega, \kappa, \Sigma$ are the rotation rate, radial epicyclic frequency and surface density of the disc, respectively), has a positive gradient at the corotation radius. Here we generalize the analysis of the corotation resonance effect to non-barotropic fluids. We show that the mode instability criterion is modified by the finite radial Brunt-V\"as\"al\"a frequency of the disc. We derive an analytic expression for the reflectivity when a density wave impinges upon the corotation barrier, and calculate the frequencies and growth rates of global p-modes for disc models with various $\alpha$-viscosity parameterizations. We find that for disc fluids with constant adiabatic index $\Gamma$, super-reflection and mode growth depend on the gradient of the effective vortensity, $\zeta_{\rm eff} = \zeta/S^{2/\Gamma}$ (where $S \equiv P/\Sigma^{\Gamma}$ measures the entropy): when $d\zeta_{\rm eff}/dr > 0$ at the corotation radius, wave absorption leads to amplification of the p-mode. Our calculations show that the lowest-order p-modes with azimuthal wave number $m=2, 3, 4,...$ have the largest growth rates, with the frequencies approximately in (but distinct from) the $2:3:4...$ commensurate ratios. We discuss the implications of our results for the high-frequency quasi-periodic oscillations observed in accreting black-hole systems.

We present numerical simulations of electron-positron pair cascades in the magnetospheres of magnetic neutron stars for a wide range of surface fields (B_p = 10^{12}--10^{15} G), rotation periods (0.1--10 s), and field geometries. This has been motivated by the discovery in recent years of a number of radio pulsars with inferred magnetic fields comparable to those of magnetars. Evolving the cascade generated by a primary electron or positron after it has been accelerated in the inner gap of the magnetosphere, we follow the spatial development of the cascade until the secondary photons and pairs leave the magnetosphere, and we obtain the pair multiplicity and the energy spectra of the cascade pairs and photons under various conditions. Going beyond previous works, which were restricted to weaker fields (B < a few x 10^{12} G), we have incorporated in our simulations detailed treatments of physical processes that are potentially important (especially in the high field regime) but were either neglected or crudely treated before, including photon splitting with the correct selection rules for photon polarization modes, one-photon pair production into low Landau levels for the e^+e^-, and resonant inverse Compton scattering from polar cap hot spots. We discuss the implications of our results for the radio pulsar death line and for the hard X-ray emission from magnetized neutron stars.

We study the propagation effects of radio waves in a pulsar magnetosphere, composed of relativistic electron-positron pair plasmas streaming along the magnetic field lines and corotating with the pulsar. We critically examine the various physical effects that can potentially influence the observed wave intensity and polarization, including resonant cyclotron absorption, wave mode coupling due to pulsar rotation, wave propagation through quasi-tangential regions (where the photon ray is nearly parallel to the magnetic field) and mode circularization due to the difference in the electron/positron density/velocity distributions. We numerically integrate the transfer equations for wave polarization in the rotating magnetosphere, taking account of all the propagation effects in a self-consistent manner. For typical magnetospheric plasma parameters produced by pair cascade, we find that the observed radio intensity and polarization profiles can be strongly modified by the propagation effects. For relatively large impact parameter (the minimum angle between the magnetic dipole axis and the line of sight), the polarization angle profile is similar to the prediction from the Rotating Vector Model, except for a phase shift and an appreciable circular polarization. For smaller impact parameter, the linear polarization position angle may exhibit a sudden $90^o$ jump due to the quasi-tangential propagation effect, accompanied by complex circular polarization profile. Some applications of our results are discussed, including the origin of non-gaussion pulse profiles, the relationship between the position angle profile and circular polarization in conal-double pulsars, and the orthogonal polarization modes.

We review the polarization properties of X-ray emission from highly magnetized neutron stars, focusing on emission from the stellar surfaces. We discuss how x-ray polarization can be used to constrain neutron star magnetic field and emission geometry, and to probe strong-field quantum electrodynamics and possibly constrain the properties of axions.

As part of a program to measure abundance ratios in stars beyond 15 kpc from the Galactic center, we have discovered a metal-poor star in the outer halo with a unique chemical signature. We originally identified it in the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey as a distant metal-poor star. We obtained a follow-up spectrum using the Echellette Spectrograph and Imager (ESI) at the Keck 2 telescope, and measure [Fe/H] =-3.17, [Mg/Fe] =-0.10 and [Ca/Fe] =+1.11. This is one of the largest over-abundances of Ca measured in any star to date; the extremely low value of [Mg/Ca]=-1.21 is entirely unique. To have found such an unusual star in our small sample of 27 targets suggests that there may be previously unobserved classes of stars yet to be found in situ in the Galactic halo.

X-ray photons emitted from the surface or atmosphere of a magnetized neutron star is highly polarized. However, the observed polarization may be modified due to photon propagation through the star's magnetosphere. For photon frequencies much larger than the typical radio frequency, vacuum birefringence due to strong-field quantum electrodynamics dominates over the plasma effect. We study the evolution of photon polarization in the magnetized QED vacuum of a neutron star magnetosphere, paying particular attention to the propagation effect across the quasi-tangential (QT) point, where the photon momentum is nearly aligned with the magnetic field. In agreement with previous studies, we find that in most regions of the magnetosphere, the photon polarization modes are decoupled due to vacuum birefringence, and therefore a large net linear polarization can be expected when the radiation escapes the magnetosphere. However, we show that X-ray polarization may change significantly when the photon passes through the QT region. When averaging over a finite emission area, the net effect of QT propagation is to reduce the degree of linear polarization; the reduction factor depends on the photon energy, magnetic field strength, geometry, rotation phase and the emission area, and can be more than a factor of two. We derive the general conditions under which the QT propagation effect is important, and provide an easy-to-use prescription to account for the QT effect for most practical calculations of X-ray polarization signals from magnetic neutron stars. For a neutron star with a dipole magnetic field, the QT effect can be important for emission from the polar cap for certain magnetic field and energy ranges, and is negligible for emission from the entire stellar surface.

The NSF's Astronomy and Astrophysics Postdoctoral Fellowship (AAPF) is exceptional among the available postdoctoral awards in Astronomy and Astrophysics. The fellowship is one of the few that allows postdoctoral researchers to pursue an original research program, of their own design, at the U.S. institution of their choice. However, what makes this fellowship truly unique is the ability of Fellows to lead an equally challenging, original educational program simultaneously. The legacy of this singular fellowship has been to encourage and advance leaders in the field who are equally as passionate about their own research as they are about sharing that research and their passion for astronomy with students and the public. In this positional paper we address the importance of fellowships like the AAPF to the astronomical profession by identifying the science and educational contributions that Fellows have made to the community. Further, we recommend that fellowships that encourage leading postdoctoral researchers to also become leaders in Astronomy education be continued and expanded.

We highlight recent theoretical and observational progress in several areas of neutron star astrophysics, and discuss the prospect for advances in the next decade.

The characteristic physical timescales near stellar-mass compact objects are measured in milliseconds. These timescales -- the free-fall time, the fastest stable orbital period, and stellar spin periods -- encode the fundamental physical properties of compact objects: mass, radius, and angular momentum. The characteristic temperature of matter in the vicinity of neutron stars is such that the principal electromagnetic window into their realms is the X-ray band. Because of these connections to the fundamental properties of neutron stars, X-ray timing studies remain today the most direct means of probing their structure and dynamics. While current X-ray observatories have revealed many relevant and fascinating phenomena, they lack the sensitivity to fully exploit them to uncover the fundamental properties of compact objects and their extreme physics. With this white paper, we summarize and highlight the science opportunities that will accompany an order-of-magnitude improvement in X-ray timing sensitivity, a goal attainable in the coming decade.

We study global non-axisymmetric oscillation modes trapped near the inner boundary of an accretion disc. Observations indicate that some of the quasi-periodic oscillations (QPOs) observed in the luminosities of accreting compact objects (neutron stars, black holes and white dwarfs) are produced in the inner-most regions of accretion discs or boundary layers. Two simple models are considered in this paper: The magnetosphere-disc model consists of a thin Keplerian disc in contact with a uniformly rotating magnetosphere with and low plasma density, while the star-disc model involves a Keplerian disc terminated at the stellar atomosphere with high density and small density scale height. We find that the interface modes at the magnetosphere-disc boundary are generally unstable due to Rayleigh-Taylor and/or Kelvin-Helmholtz instabilities. However, differential rotation of the disc tends to suppress Rayleigh-Taylor instability and a sufficiently high disc sound speed (or temperature) is needed to overcome this suppression and to attain net mode growth. On the other hand, Kelvin-Helmholtz instability may be active at low disc sound speeds. We also find that the interface modes trapped at the boundary between a thin disc and an unmagnetized star do not suffer Rayleigh-Taylor or Kelvin-Helmholtz instability, but can become unstable due to wave leakage to large disc radii and, for sufficiently steep disc density distributions, due to wave absorption at the corotation resonance in the disc. The non-axisymmetric interface modes studied in this paper may be relevant to the high-frequency QPOs observed in some X-ray binaries and in cataclysmic variables.

Diskoseismic c-modes in accretion discs have been invoked to explain low-frequency variabilities observed in black-hole X-ray binaries. These modes are trapped in the inner-most region of the disc and have frequencies much lower than the rotation frequency at the disc inner radius. We show that because the trapped waves can tunnel through the evanescent barrier to the corotational wave zone, the c-modes are damped due to wave absorption at the corotation resonance. We calculate the corotational damping rates of various c-modes using the WKB approximation. The damping rate varies widely depending on the mode frequency, the black hole spin parameter and the disc sound speed, and is generally much less than 10% of the mode frequency. A sufficiently strong excitation mechanism is needed to overcome this corotational damping and make the mode observable.

We study the global stability of non-axisymmetric p-modes (also called inertial-acoustic modes) trapped in the inner-most regions of accretion discs around black holes. We show that the lowest-order (highest-frequency) p-modes, with frequencies $\omega=(0.5-0.7) m\Omega_{\rm ISCO}$, can be overstable due to general relativistic effects, according to which the radial epicyclic frequency is a non-monotonic function of radius near the black hole. The mode is trapped inside the corotation resonance radius and carries a negative energy. The mode growth arises primarily from wave absorption at the corotation resonance, and the sign of the wave absorption depends on the gradient of the disc vortensity. When the mode frequency is sufficiently high, such that the slope of the vortensity is positive at corotation positive wave energy is absorbed at the resonance, leading to the growth of mode amplitude. We also study how the rapid radial inflow at the inner edge of the disc affects the mode trapping and growth. Our analysis of the behavior of the fluid perturbations in the transonic flow near the ISCO indicates that, while the inflow tends to damp the mode, the damping effect is sufficiently small under some conditions so that net mode growth can still be achieved. We further clarify the role of the Rossby wave instability and show that it does not operate for black hole accretion discs with smooth-varying vortensity profiles. Overstable non-axisymmetric p-modes driven by the corotational instability provide a plausible explanation for the high-frequency (> 100 Hz) quasi-periodic oscillations (HFQPOs) observed from a number of black-hole X-ray binaries in the very high state. The absence of HFQPOs in the soft (thermal) state may result from mode damping due to the radial infall at the ISCO.

The origin of the rapid quasi-periodic variabilities observed in a number of accreting black hole X-ray binaries is not understood. It has been suggested that these variabilities are associated with diskoseismic oscillation modes of the black hole accretion disk. In particular, in a disk with no magnetic field, the so-called g-modes (inertial oscillations) can be self-trapped at the inner region of the disk due to general relativistic effects. Real accretion disks, however, are expected to be turbulent and contain appreciable magnetic fields. We show in this paper that even a weak magnetic field (with the magnetic energy much less than the thermal energy) can modify or "destroy" the self-trapping zone of disk g-modes, rendering their existence questionable in realistic black hole accretion disks. The so-called corrugation modes (c-modes) are also strongly affected when the poloidal field approaches equal-partition. On the other hand, acoustic oscillations (p-modes), which do not have vertical structure, are not affected qualitatively by the magnetic field, and therefore may survive in a turbulent, magnetic disk.

We report the discovery of a strong magnetic field in the unique pulsating carbon-atmosphere white dwarf SDSS J142625.71+575218.3. From spectra gathered at the MMT and Keck telescopes, we infer a surface field of B_s ~1.2 MG, based on obvious Zeeman components seen in several carbon lines. We also detect the presence of a Zeeman-splitted He I 4471 line, which is an indicator of the presence of a non-negligible amount of helium in the atmosphere of this Hot DQ star. This is important for understanding its pulsations, as nonadabatic theory reveals that some helium must be present in the envelope mixture for pulsation modes to be excited in the range of effective temperature where the target star is found. Out of nearly 200 pulsating white dwarfs known today, this is the first example of a star with a large detectable magnetic field. We suggest that SDSS J142625.71+575218.3 is the white dwarf equivalent of a roAp star.

We present the results of an abundance analysis for a sample of stars with $-4<$[Fe/H]$<-2$. The data were obtained with the HIRES spectrograph at Keck Observatory. The set includes 28 stars, with effective temperature ranging from 4800 to 6600 K. For 13 stars with [Fe/H]$<-2.6$, including nine with [Fe/H]$<-3.0$ and one with [Fe/H]$=-4.0$, these are the first reported detailed abundances. For the most metal-poor star in our sample, CS 30336-049, we measure an abundance pattern that is very similar to stars in the range [Fe/H]$\sim-3.5$, including a normal C+N abundance. We also find that it has very low but measurable Sr and Ba, indicating some neutron-capture activity even at this low of a metallicity. We explore this issue further by examining other very neutron-capture-deficient stars, and find that at the lowest levels, [Ba/Sr] exhibits the ratio of the main r-process. We also report on a new r-process-enhanced star, CS 31078-018. This star has [Fe/H]$=-2.85$, [Eu/Fe]$=1.23$, and [Ba/Eu]$=-0.51$. CS 31078-018 exhibits an ``actinide boost'', i.e. much higher [Th/Eu] than expected and at a similar level to CS 31082-001. Our spectra allow us to further constrain the abundance scatter at low metallicities, which we then use to fit to the zero-metallicity Type II supernova yields of Heger & Woosley (2008). We find that supernovae with progenitor masses between 10 and 20 M$_{\odot}$ provide the best matches to our abundances.

In differentially rotating discs with no self-gravity, density waves cannot propagate around the corotation, where the wave pattern rotation speed equals the fluid rotation rate. Waves incident upon the corotation barrier may be super-reflected (commonly referred to as corotation amplifier), but the reflection can be strongly affected by wave absorptions at the corotation resonance/singularity. The sign of the absorption is related to the Rossby wave zone very near the corotation radius. We derive the explicit expressions for the complex reflection and transmission coefficients, taking into account wave absorption at the corotation resonance. We show that for generic discs, this absorption plays a much more important role than wave transmission across the corotation barrier. Depending on the sign of the gradient of the specific vorticity of the disc the corotation resonance can either enhance or diminish the super-reflectivity, and this can be understood in terms of the location of the Rossby wave zone relative to the corotation radius. Our results provide the explicit conditions (in terms of disc thickness, rotation profile and specific vorticity gradient) for which super-reflection can be achieved. Global overstable disc modes may be possible for discs with super-reflection at the corotation barrier.

The accretion disk around a rotating magnetic star (neutron star, white dwarf or T Tauri star) is subjected to periodic vertical magnetic forces from the star, with the forcing frequency equal to the stellar spin frequency or twice the spin frequency. This gives rise bending waves in the disk that may influence the variabilities of the system. We study the excitation, propagation and dissipation of these waves using a hydrodynamical model coupled with a generic model description of the magnetic forces. The $m=1$ bending waves are excited at the Lindblad/vertical resonance, and propagate either to larger radii or inward toward the corotation resonance where dissipation takes place. While the resonant torque is negligible compared to the accretion torque, the wave nevertheless may reach appreciable amplitude and can cause or modulate flux variabilities from the system. We discuss applications of our result to the observed quasi-periodic oscillations from various systems, in particular neutron star low-mass X-ray binaries.

For sufficiently strong magnetic fields and/or low temperatures, the neutron star surface may be in a condensed state with little gas or plasma above it. Such surface condensation can significantly affect the thermal emission from isolated neutron stars, and may lead to the formation of a charge-depleted acceleration zone ("vacuum gap") in the magnetosphere above the stellar polar cap. Using the latest results on the cohesive property of magnetic condensed matter, we quantitatively determine the conditions for surface condensation and vacuum gap formation in magnetic neutron stars. We find that condensation can occur if the thermal energy kT of the neutron star surface is less than about 8% of its cohesive energy Q_s, and that a vacuum gap can form if the neutron star's rotation axis and magnetic moment point in opposite directions and kT is less than about 4% of Q_s. Thus, vacuum gap accelerators may exist for some neutron stars. Motivated by this result, we also study the physics of pair cascades in the vacuum gap model for photon emission by accelerating electrons and positrons due to both curvature radiation and resonant/nonresonant inverse Compton scattering. Our calculations of the condition of cascade-induced vacuum breakdown and the related pulsar death line/boundary generalize previous works to the superstrong field regime. We find that inverse Compton scatterings do not produce a sufficient number of high energy photons in the gap and thus do not lead to pair cascades for most neutron star parameters. We discuss the implications of our results for the recent observations of neutron star thermal radiation as well as for the detection/non-detection of radio emission from high-B pulsars and magnetars.

Recent observations of thermally emitting isolated neutron stars revealed spectral features that could be interpreted as radiative transitions of He in a magnetized neutron star atmosphere. We present Hartree-Fock calculations of the polarization-dependent photoionization cross sections of the He atom in strong magnetic fields ranging from 10^12 G to 10^14 G. Convenient fitting formulae for the cross sections are given as well as related oscillator strengths for various bound-bound transitions. The effects of finite nucleus mass on the radiative absorption cross sections are examined using perturbation theory.

We report high-resolution spectroscopy of 125 field stars previously observed as part of the Sloan Digital Sky Survey and its program for Galactic studies, the Sloan Extension for Galactic Understanding and Exploration (SEGUE). These spectra are used to measure radial velocities and to derive atmospheric parameters, which we compare with those reported by the SEGUE Stellar Parameter Pipeline (SSPP). The SSPP obtains estimates of these quantities based on SDSS ugriz photometry and low-resolution (R = 2000) spectroscopy. For F- and G-type stars observed with high signal-to-noise ratios (S/N), we empirically determine the typical random uncertainties in the radial velocities, effective temperatures, surface gravities, and metallicities delivered by the SSPP to be 2.4 km/s, 130 K (2.2%), 0.21 dex, and 0.11 dex, respectively, with systematic uncertainties of a similar magnitude in the effective temperatures and metallicities. We estimate random errors for lower S/N spectra based on numerical simulations.

We present initial results measuring the companion fraction of metal-poor stars ([Fe/H]$<-$2.0). We are employing the Lick Observatory planet-finding system to make high-precision Doppler observations of these objects. The binary fraction of metal-poor stars provides important constraints on star formation in the early Galaxy (Carney et al. 2003). Although it has been shown that a majority of solar metallicity stars are in binaries, it is not clear if this is the case for metal-poor stars. Is there a metallicity floor below which binary systems do not form or become rare? To test this we are determining binary fractions at metallicities below [Fe/H]$=-2.0$. Our measurments are not as precise as the planet finders', but we are still finding errors of only 50 to 300 m/s, depending on the signal-to-noise of a spectrum and stellar atmosphere of the star. At this precision we can be much more complete than previous studies in our search for stellar companions.

Hall drift, i. e., transport of magnetic flux by the moving electrons giving rise to the electrical current, may be the dominant effect causing the evolution of the magnetic field in the solid crust of neutron stars. It is a nonlinear process that, despite a number of efforts, is still not fully understood. We use the Hall induction equation in axial symmetry to obtain some general properties of nonevolving fields, as well as analyzing the evolution of purely toroidal fields, their poloidal perturbations, and current-free, purely poloidal fields. We also analyze energy conservation in Hall instabilities and write down a variational principle for Hall equilibria. We show that the evolution of any toroidal magnetic field can be described by Burgers' equation, as previously found in plane-parallel geometry. It leads to sharp current sheets that dissipate on the Hall time scale, yielding a stationary field configuration that depends on a single, suitably defined coordinate. This field, however, is unstable to poloidal perturbations, which grow as their field lines are stretched by the background electron flow, as in instabilities earlier found numerically. On the other hand, current-free poloidal configurations are stable and could represent a long-lived crustal field supported by currents in the fluid stellar core.

We present carbon and strontium abundances for 100 metal-poor stars measured from R$\sim $7000 spectra obtained with the Echellette Spectrograph and Imager at the Keck Observatory. Using spectral synthesis of the G-band region, we have derived carbon abundances for stars ranging from [Fe/H]$=-1.3$ to [Fe/H]$=-3.8$. The formal errors are $\sim 0.2$ dex in [C/Fe]. The strontium abundance in these stars was measured using spectral synthesis of the resonance line at 4215 {\AA}. Using these two abundance measurments along with the barium abundances from our previous study of these stars, we show it is possible to identify neutron-capture-rich stars with our spectra. We find, as in other studies, a large scatter in [C/Fe] below [Fe/H]$ = -2$. Of the stars with [Fe/H]$<-2$, 9$\pm$4% can be classified as carbon-rich metal-poor stars. The Sr and Ba abundances show that three of the carbon-rich stars are neutron-capture-rich, while two have normal Ba and Sr. This fraction of carbon enhanced stars is consistent with other studies that include this metallicity range.

We study the wave propagation modes in the relativistic streaming pair plasma of the magnetospheres of pulsars and magnetars, focusing on the effect of vacuum polarization. We show that the combined plasma and vacuum polarization effects give rise to a vacuum resonance, where ``avoided mode crossing'' occurs between the extraordinary mode and the (superluminous) ordinary mode. When a photon propagates from the vacuum-polarization-dominated region at small radii to the plasma-dominated region at large radii, its polarization state may undergo significant change across the vacuum resonance. We map out the parameter regimes (e.g., field strength, plasma density and Lorentz factor) under which the vacuum resonance occurs and examine how wave propagation is affected by the resonance. Some possible applications of our results are discussed, including high-frequency radio emission from pulsars and possibly magnetars, and optical/IR emission from neutron star surfaces and inner magnetospheres.

The cohesive energy of condensed matter in strong magnetic fields is a fundamental quantity characterizing magnetized neutron star surfaces. The cohesive energy refers to the energy required to pull an atom out of the bulk condensed matter at zero pressure. Theoretical models of pulsar and magnetar magnetospheres depend on the cohesive properties of the surface matter in strong magnetic fields. For example, depending on the cohesive energy of the surface matter, an acceleration zone ("polar gap") above the polar cap of a pulsar may or may not form. Also, condensation of the neutron star surface, if it occurs, can significantly affect thermal emission from isolated neutron stars. We describe our calculations of the cohesive property of matter in strong magnetic fields, and discuss the implications of our results to the recent observations of neutron star surface emission as well as to the detection/non-detection of radio emission from magnetars.

We calculate the transition rates between proton Landau levels due to non-radiative and radiative Coulomb collisions in an electron-proton plasma with strong magnetic field B. Both electron-proton collisions and proton-proton collisions are considered. The roles of the first-order cyclotron absorption and second-order free-free absorption and scattering in determining the line strength and shape as well as the continuum are analysed in detail. We solve the statistical balance equation for the populations of proton Landau levels. For temperatures \sim 10^6-10^7 K, the deviations of the proton populations from LTE are appreciable at density \rho < 0.1 B_{14}^{3.5} g cm^{-3}, where B_{14}=B/(10^{14} G). We present general formulae for the plasma emissivity and absorption coefficents under a wide range of physical conditions. Our results are useful for studying the possibility and the conditions of proton/ion cyclotron line formation in the near vicinity of highly magnetized neutron stars.

We present new ab initio calculations of the electronic structure of various atoms and molecules in strong magnetic fields ranging from B=10^12 G to 2x10^15 G, appropriate for radio pulsars and magnetars. For these field strengths, the magnetic forces on the electrons dominate over the Coulomb forces, and to a good approximation the electrons are confined to the ground Landau level. Our calculations are based on the density functional theory, and use a local magnetic exchange-correlation function which is tested to be reliable in the strong field regime. Numerical results of the ground-state energies are given for H_N (up to N=10), He_N (up to N=8), C_N (up to N=5) and Fe_N (up to N=3), as well as for various ionized atoms. Fitting formulae for the B-dependence of the energies are also given. In general, as N increases, the binding energy per atom in a molecule, |E_N|/N, increases and approaches a constant value. For all the field strengths considered in this paper, hydrogen, helium, and carbon molecules are found to be bound relative to individual atoms (although for B less than a few x 10^12 G, the relative binding between C and C_2 is small). Iron molecules are not bound at B<10^13 G, but become energetically more favorable than individual atoms at larger field strengths.

We present new, ab initio calculations of the electronic structure of one-dimensional infinite chains and three-dimensional condensed matter in strong magnetic fields ranging from B=10^12 G to 2x10^15 G, appropriate for observed magnetic neutron stars. At these field strengths, the magnetic forces on the electrons dominate over the Coulomb forces, and to a good approximation the electrons are confined to the ground Landau level. Our calculations are based on the density functional theory, and use a local magnetic exchange-correlation function appropriate in the strong field regime. The band structures of electrons in different Landau orbitals are computed self-consistently. Numerical results of the ground-state energies and electron work functions are given for one-dimensional chains of H, He, C, and Fe. Fitting formulae for the B-dependence of the energies are also provided. For all the field strengths considered in this paper, hydrogen, helium, and carbon chains are found to be bound relative to individual atoms (although for B less than a few x 10^12 G, the relative binding between C and C_infinity is small). Iron chains are significantly bound for B>10^14 G and are weakly bound if at all at B<10^13 G. We also study the cohesive property of three-dimensional condensed matter of H, He, C, and Fe at zero pressure, constructed from interacting chains in a body-centered tetragonal lattice. Such three-dimensional condensed matter is found to be bound relative to individual atoms, with the cohesive energy increasing rapidly with increasing B.

Recent experiments suggest that polarized photons may couple significantly to pseudoscalar particles such as axions. We study the possible observational signatures of axion-photon coupling for radiation from magnetic stars, with particular focus on neutron stars. We present general methods for calculating the axion-photon conversion probability during propagation through a varying magnetized vacuum as well as across an inhomogeneous atmosphere. Partial axion-photon conversion may take place in the vacuum region outside the neutron star. Strong axion-photon mixing occurs due to a resonance in the atmosphere, and depending on the axion coupling strength and other parameters, significant axion-photon conversion can take place at the resonance. Such conversions may produce observable effects on the radiation spectra and polarization signals from the star. We also apply our results to axion-photon propagation in the Sun and in magnetic white dwarfs. We find that there is no appreciable conversion of solar axions to photons during the propagation.

Recent observations of pulsar wind nebulae and radio polarization profiles revealed a tendency of the alignment between the spin and velocity directions in neutron stars. We study the condition for spin-kick alignment using a toy model, in which the kick consists of many off-centered, randomly-oriented thrusts. Both analytical considerations and numerical simulations indicate that spin-kick alignment cannot be easily achieved if the proto-neutron star does not possess some initial angular momentum, contrary to some previous claims. To obtain the observed spin-kick misalignment angle distribution, the initial spin period of the neutron star must be smaller than the kick timescale. Typically, an initial period of a hundred milliseconds or less is required.

Theoretical modeling of surface emission from magnetized neutron stars (NSs) requires proper treatment of QED effects, in particular the effect of photon mode conversion due to the ``vacuum resonance'' between plasma and vacuum polarization. Previous NS atmosphere models incorporated this effect approximately, using transfer equations for the photon modes which are inadequate near the vacuum resonance, particularly for field strengths around $B_l\simeq 7\times 10^{13}$ G, where the vacuum resonance occurs near the photosphere. In this paper, we provide an accurate treatment of the QED-induced mode conversion effect in NS atmosphere models, employing both the modal radiative transfer equations, coupled with an accurate mode conversion probability at the resonance, and the full evolution equations for the photon Stokes parameters. In doing so, we are able to quantitatively calculate the effects of vacuum polarization on the atmosphere emission spectra, beam patterns, and polarizations for the entire range of field strengths, $B=10^{12}-10^{15}$ G. We find that for NSs with $B\ga 2 B_l$, vacuum polarization reduces the widths of spectral features, and softens the hard spectral tail typical of magnetized atmosphere models. For $B\la B_l/2$, vacuum polarization does not change the emission spectra, but can significantly affect the polarization signals. We show that vacuum polarization induces a unique energy-dependent linear polarization signature, and that circular polarization can be generated in the magnetospheres of rapidly rotating NSs. We discuss the implications of our results for recent observations of thermally emitting isolated NSs and magnetars, as well as the prospects for future spectral and polarization observations.

There has recently been growing evidence for the existence of neutron stars possessing magnetic fields with strengths that exceed the quantum critical field strength of $4.4 \times 10^{13}$ G, at which the cyclotron energy equals the electron rest mass. Such evidence has been provided by new discoveries of radio pulsars having very high spin-down rates and by observations of bursting gamma-ray sources termed magnetars. This article will discuss the exotic physics of this high-field regime, where a new array of processes becomes possible and even dominant, and where familiar processes acquire unusual properties. We review the physical processes that are important in neutron star interiors and magnetospheres, including the behavior of free particles, atoms, molecules, plasma and condensed matter in strong magnetic fields, photon propagation in magnetized plasmas, free-particle radiative processes, the physics of neutron star interiors, and field evolution and decay mechanisms. Application of such processes in astrophysical source models, including rotation-powered pulsars, soft gamma-ray repeaters, anomalous X-ray pulsars and accreting X-ray pulsars will also be discussed. Throughout this review, we will highlight the observational signatures of high magnetic field processes, as well as the theoretical issues that remain to be understood.

We study the effect of resonant tidal excitation of inertial modes in neutron stars during binary inspiral. For spin frequencies less than 100 Hz, the phase shift in the gravitational waveform associated with the resonance is small and does not affect the matched filtering scheme for gravitational wave detection. For higher spin frequencies, the phase shift can become significant. Most of the resonances take place at orbital frequencies comparable to the spin frequency, and thus significant phase shift may occur only in the high-frequency band (hundreds of Hertz) of gravitational wave. The exception is a single odd-paity $m=1$ mode, which can be resonantly excited for misaligned spin-orbit inclinations, and may occur in the low-frequency band (tens of Hertz) of gravitational wave and induce significant (>> 1 radian) phase shift.

We study the effect of pulsar rotation on timing of binary pulsars, with particular emphasis on the double pulsar system J0737-3039. Special relativistic aberration due to the orbital motion of pulsar changes both the longitude and colatitude of the emission direction with respect to the pulsar spin axis. The former gives rise to a shift of the arrival time of the pulse centroid (this is the conventional ``longitudinal'' aberration delay), the latter results in a distortion (contraction or dilation) of the pulse profile on the orbital timescale. The amplitude of pulse distortion depends inversely on the variation of polarization position angle across the pulse. For small angle between the pulsar magnetic and spin axes, as inferred for PSR J0737-3039A from polarimetric observations, the pulse distortion is significant (\sim 1%) and the associated ``latitudinal'' aberration delay is much larger than the longitudinal one. By monitoring the arrival time of separate pulse components as a function of pulsar orbital phase, the latitudinal aberration delay may be easily measured with the current timing precision, thus constraining the spin geometry of the system. The latitudinal delay can also be detected by monitoring system's orbital parameters on the geodetic precession timescale. Because of the near edge-on orbital orientation of the PSR J0737-3039 system, general relativistic bending of pulsar A's radio beam near its superior conjunction also introduces spin-dependent time delays of similar order of magnitude as the aberration delays. In addition, light bending splits the pulse profile into two variable components, corresponding to two gravitationally lensed images of the source. Detection of lensing effects is challenging, but may be possible with existing technology.

The measurement of the Shapiro time delay in binary pulsar systems with highly inclined orbit can be affected both by the motion of the pulsar's companion because of the finite time it takes a photon to cross the binary, and by the gravitational light bending if the orbit is sufficiently edge-on relative to the line of sight. Here we calculate the effect of retardation due to the companion's motion on various time delays in pulsar binaries, including the Shaipro delay, the geometric lensing delay, and the lens-induced delays associated with the pulsar rotation. Our results can be applied to systems so highly inclined that near conjunction gravitational lensing of the pulsar radiation by the companion becomes important (the recently discovered double pulsar system J0737-3039 may exemplify such a system). To the leading order, the effect of retardation is to shift all the delay curves backward in time around the orbit conjunction, without affecting the shape and amplitude of the curves. The time shift is of order the photon orbit crossing time, and ranges from a second to a few minutes for the observed binary pulsar systems. In the double pulsar system J0737-3039, the motion of the companion may also affect the interpretation of the recent correlated interstellar scintillation measurements. Finally, we show that lensing sets an upper limit on the magnitude of the frame-dragging time delay caused by the companion's spin, and makes this delay unobservable in stellar-mass binary pulsar systems.

We study the excitation of density and bending waves and the associated angular momentum transfer in gaseous disks with finite thickness by a rotating external potential. The disk is assumed to be isothermal in the vertical direction and has no self-gravity. The disk perturbations are decomposed into different modes, each characterized by the azimuthal index m and the vertical index n, which specifies the nodal number of the density perturbation along the disk normal direction. The n=0 modes correspond to the two-dimensional density waves previously studied by Goldreich & Tremaine and others. In a three-dimensional disk, waves can be excited at both Lindblad resonances (for modes with n=0,1,2...) and vertical resonances (for the $n\ge 1$ modes only). The torque on the disk is positive for waves excited at outer Linblad/vertical resonances and negative at inner Lindblad/vertical resonances. While the n=0 modes are evanescent around corotation, the $n\ge 1$ modes can propagate into the corotation region where they are damped and deposit their angular momenta. We have derived analytical expressions for the amplitudes of different wave modes excited at Lindblad and/or vertical resonances and the resulting torques on the disk. It is found that for $n\ge 1$, angular momentum transfer through vertical resonances is much more efficient than Lindblad resonances. This implies that in some situations (e.g., a circumstellar disk perturbed by a planet in an inclined orbit), vertical resonances may be an important channel of angular momentum transfer between the disk and the external potential. We have also derived new formulae for the angular momentum deposition at corotation and studied wave excitations at disk boundaries.

We study the equation of state, polarization and radiation properties for nonideal, strongly magnetized plasmas which compose outer envelopes of magnetic neutron stars. Detailed calculations are performed for partially ionized hydrogen atmospheres and for condensed hydrogen or iron surfaces of these stars. This is a companion paper to astro-ph/0511803

We study observational constraints on neutron star (NS) kicks for isolated pulsars and for neutron stars in binary systems. We are particularly interested in the evidence of kick-spin alignment/misalignment and its dependence on the neutron star initial spin period. For several young pulsars, X-ray observations of compact nebulae showed that pulsar proper motion is aligned with the spin direction as defined by the symmetry axis of the nebula. We also critically examine the measurements of the proper motion and the projected spin axis from a large sample of pulsars with well-calibrated polarization data. We find that among the two dozen pulsars for which reliable measurements are available, there is a significant correlation between the spin axis and the proper motion. For various NS binaries, including double NS systems, binaries with massive main-sequence star companion and binaries with massive white-dwarf companion, we obtain constraints on the kick magnitudes and directions from the observed orbital characteristics of the system. We find that the kick velocity is misaligned with the the NS spin axis in a number of systems, and the NS spin period (when available) in these systems is generally longer than several hundreds milliseconds. These constraints, together with spin-kick alignment observed in many isolated pulsars, suggest that the kick timescale is hundreds of milliseconds to 1 s, so that spin-kick alignment or misalignment can be obtained depending on the initial spin period of the NS. We discuss the implication of our result for various NS kick mechanisms.

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 present our view of the main physical ingredients determining the evolution of neutron star magnetic fields. This includes the basic properties of neutron star matter, possible scenarios for the origin of the magnetic field, constraints and mechanisms for its evolution, and a discussion of our recent work on the Hall drift.

We investigate the effects of gravitational lensing in the binary pulsar system J0737-3039. Current measurement of the orbital inclination allows the millisecond pulsar (A) to pass very close (at R_{min}=4000 km) in projection to the companion pulsar (B), with R_{min} comparable to the Einstein radius (2600 km). For this separation at the conjunction, lensing causes small (about 10%) magnification of the pulsar A signal on a timescale of several seconds, and displaces the pulsar image on the sky plane by about 1200 km. More importantly, lensing introduces a correction (of several microsec) to the conventional Shapiro delay formula used in pulsar timing analysis, and gives rise to a geometric time delay together with the delays associated with the pulsar spin period. These lensing effects can influence the determination of the system parameters by both timing and scintillation studies. Given the current uncertainty in the orbital inclination, more extreme manifestations of lensing (e.g. magnification by a factor of several) are possible. We compare our predictions with the existing observations and discuss the possibility of detecting gravitational lensing signatures in the system. The anomalously high point in A's lightcurve close to superior conjunction might be caused by gravitational lensing.

We construct hydrogen atmosphere models for strongly magnetized neutron stars in thermodynamic equilibrium, taking into account partial ionization. The presence of bound states affects the equation of state, absorption coefficients, and polarizability tensor of a strongly magnetized plasma. Therefore the partial ionization influences the polarization vectors and opacities of normal electromagnetic waves, and thus the spectra of outgoing radiation. Here we review a model suitable for the most typical neutron-star atmospheres and focus on the problems that remain to be solved for its extension to other atmospheric parameters.

We present chemical abundances for 110 stars identified in objective-prism surveys as candidates to be very metal-poor. The abundances are derived from high S/N, intermediate-resolution spectra obtained with the Keck Observatory Echelle Spectrometer and Imager. An additional 25 stars with well-determined abundances ranging from [Fe/H]$=-1.5$ and -3.2 were observed and the results used to help calibrate our analysis and determine the accuracy of our abundance determinations. Abundances for the program stars were measured for Fe, Mg, Ca, Ti, Cr and Ba with an accuracy of approximately 0.3 dex. 53 of the stars in our sample have [Fe/H]<= -2, 22 have [Fe/H]<= -2.5 and 13 stars have [Fe/H]<= -2.9. Surprisingly, approximately one third of the sample is relatively metal rich with [Fe/H]>-1.5. In addition to identifying a number of extremely metal-poor stars, this study also shows that moderate-resolution spectra obtained with the Keck Echelle Spectrometer and Imager yield relatively accurate abundances for stars as faint as V=14 in modest exposure time (~20 minutes). This capability will prove useful if the so-far elusive stars at [Fe/H]<-4 turn out to be mostly fainter than V=15.

Polarizability tensor of a strongly magnetized plasma and the polarization vectors and opacities of normal electromagnetic waves are studied for the conditions typical of neutron star atmospheres, taking account of partial ionization effects. Vacuum polarization is also included using a new set of fitting formulae that are accurate for wide range of field strengths. The full account of the coupling of the quantum mechanical structure of the atoms to their center-of-mass motion across the magnetic field is shown to be crucial for the correct evaluation of the polarization properties and opacities of the plasma. The self-consistent treatment of the polarizability and absorption coefficients proves to be necessary if the ionization degree of the plasma is low, which can occur in the atmospheres of cool or ultramagnetized neutron stars. Atmosphere models and spectra are presented to illustrate the importance of such self-consistent treatment.

We study several effects that influence the strength of the proton cyclotron and atomic features in the thermal spectra of magnetic neutron stars. We show that it is possible for vacuum polarization to strongly suppress the spectral lines when the magnetic field B\ga 10^14 G. For weaker fields (B\la 7\times 10^13 G), the surface spectrum is unaffected by vacuum polarization; thus the proton cyclotron absorption line can have a large equivalent width. Using an approximate calculation of the synthetic spectra, we show that variation of magnetic fields over the neutron star surface leads to broadening of the line. The recently detected absorption features from the isolated neutron stars RX J1308.6+2127, RX J1605.3+3249, and RX J0720.4-3125 can plausibly be explained by the proton cyclotron resonance, with possible blending due to atomic transitions, in the atmosphere of the star.

The inner region of the accretion disk around a magnetized star (T Tauri star, white dwarf or neutron star) is subjected to magnetic torques that induce warping and precession of the disk. These torques arise from the interaction between the stellar field and the induced electric currents in the disk. We carry out numerical simulations of the nonlinear evolution of warped, viscous accretion disks driven by the magnetic torques. We show that the disk can develop into a highly warped steady state in which the disk attains a fixed (warped) shape and precesses rigidly. The warp is most pronounced at the disk inner radius (near the magnetosphere boundary). As the system parameters (such as accretion rate) change, the disk can switch between a completely flat state (warping stable) and a highly warped state. The precession of warped disks may be responsible for a variety of quasi-periodic oscillations or radiation flux variabilities observed in many different systems, including young stellar objects and X-ray binaries.

Observations over the last decade have shown that neutron stars receive a large kick velocity (of order a hundred to a thousand km/s) at birth. The physical origin of the kicks and the related supernova asymmetry is one of the central unsolved mysteries of supernova research. We review the physics of different kick mechanisms, including hydrodynamically driven, neutrino -- magnetic field driven, and electromagnetically driven kicks. The viabilities of the different kick mechanisms are directly related to the other key parameters characterizing nascent neutron stars, such as the initial magnetic field and the initial spin. Recent observational constraints on kick mechanisms are also discussed.

We construct partially ionized hydrogen atmosphere models for magnetized neutron stars in radiative equilibrium with surface fields B=10^12-5 \times 10^14 G and effective temperatures T_eff \sim a few \times 10^5-10^6 K. These models are based on the latest equation of state and opacity results for magnetized, partially ionized hydrogen plasmas that take into account various magnetic and dense medium effects. The atmospheres directly determine the characteristics of thermal emission from isolated neutron stars. For the models with B=10^12-10^13 G, the spectral features due to neutral atoms lie at extreme UV and very soft X-ray energy bands and therefore are difficult to observe. However, the continuum flux is also different from the fully ionized case, especially at lower energies. For the superstrong field models (B\ga 10^14 G), we show that the vacuum polarization effect not only suppresses the proton cyclotron line as shown previously, but also suppresses spectral features due to bound species; therefore spectral lines or features in thermal radiation are more difficult to observe when the neutron star magnetic field is \ga 10^14 G.

In the atmospheric plasma of a strongly magnetized neutron star, vacuum polarization can induce a Mikheyev-Smirnov-Wolfenstein like resonance across which a X-ray photon may (depending on its energy) convert from one mode into the other, with significant changes in opacities and polarizations. We show that this vacuum resonance effect gives rise to an unique energy-dependent polarization signature in the surface emission from neutron stars. The detection of polarized X-rays from neutron stars can provide a direct probe of strong-field quantum electrodynamics and constrain the neutron star magnetic field and geometry.

Current theoretical models for the outflows/jets from AGN, X-ray binaries and young stellar objects involve large-scale magnetic fields threading an underlying accretion disk. We suggest that such a disk is subjected to warping instability and retrograde precession driven by magnetic torques associated with the outflow. The growth timescale for the disk warp and the precession period are of order the radial infall time of the disk. These effects may be relevant to jet precession and other variabilities observed in many systems.

Thermal surface emissions have now been detected from more than a dozen isolated neutron stars, including radio pulsars, radio-quiet neutron stars and magnetars. These detections can potentially provide important information on the interior physics, magnetic fields, and surface composition of neutron stars. Understanding the properties of matter and radiative transfer in strong magnetic fields is essential for the proper interpretation of the observations. We review current theoretical works on modeling magnetized neutron star atmospheres/surface layers, discussing some of the novel properties of matter and radiative transfer in strong magnetic fields. Of particular interest is the effect of the strong-field vacuum polarization, which dramatically changes the radiative transfer and the emergent X-ray spectra from magnetars.

We present a theoretical study of radiative transfer in strongly magnetized electron-ion plasmas, focusing on the effect of vacuum polarization due to quantum electrodynamics. This study is directly relevant to thermal radiation from the surfaces of highly magnetized neutron stars, which have been detected in recent years. Strong-field vacuum polarization modifies the photon propagation modes in the plasma, and induces a ``vacuum resonance'' at which a polarized X-ray photon propagating outward in the neutron star atmosphere can convert from a low-opacity mode to a high-opacity mode and vice versa. The effectiveness of this mode conversion depends on the photon energy and the atmosphere density gradient. For a wide range of field strengths, $7\times 10^{13}< B < 10^{16}$ G, the vacuum resonance lies between the photospheres of the two photon modes, and the emergent radiation spectrum from the neutron star is significantly modified by the vacuum resonance. (For lower field strengths, only the polarization spectrum is affected.) Under certain conditions, which depend on the field strength, photon energy and propagation direction, the vacuum resonance is accompanied by the phenomenon of mode collapse (at which the two photon modes become degenerate) and the breakdown of Faraday depolarization. Thus, the widely used description of radiative transfer based on photon modes is not adequate to treat the vacuum polarization effect rigorously. We study the evolution of polarized X-rays across the vacuum resonance and derive the transfer equation for the photon intensity matrix (Stokes parameters), taking into account the effect of birefringence of the plasma-vacuum medium, free-free absorption, and scatterings by electrons and ions.

We study the effect of vacuum polarization on the atmosphere structure and radiation spectra of neutron stars with surface magnetic fields B=10^14-10^15 G, as appropriate for magnetars. Vacuum polarization modifies the dielectric property of the medium and gives rise to a resonance feature in the opacity; this feature is narrow and occurs at a photon energy that depends on the plasma density. Vacuum polarization can also induce resonant conversion of photon modes via a mechanism analogous to the MSW mechanism for neutrino oscillation. We construct atmosphere models in radiative equilibrium with an effective temperature of a few \times 10^6 K by solving the full radiative transfer equations for both polarization modes in a fully ionized hydrogen plasma. We discuss the subtleties in treating the vacuum polarization effects in the atmosphere models and present approximate solutions to the radiative transfer problem which bracket the true answer. We show from both analytic considerations and numerical calculations that vacuum polarization produces a broad depression in the X-ray flux at high energies (a few keV \la E \la a few tens of keV) as compared to models without vacuum polarization; this arises from the density dependence of the vacuum resonance feature and the large density gradient present in the atmosphere. Thus the vacuum polarization effect softens the high energy tail of the thermal spectrum, although the atmospheric emission is still harder than the blackbody spectrum because of the non-grey opacities. We also show that the depression of continuum flux strongly suppresses the equivalent width of the ion cyclotron line and therefore makes the line more difficult to observe.

Recent observations have revealed many new puzzles related to core-collapse supernovae, including the formation of magnetars and black holes and their possible GRB connections. We review our current understanding of the origin of pulsar kicks and supernova asymmetry. It is argued that neutron star kicks are intimately connected to the other fundamental parameters of young neutron stars, such as the initial spin and magnetic field strength.

We construct hydrogen atmosphere models for magnetized neutron stars in radiative equilibrium with surface fields B=10^12-5x10^14 G and effective temperatures T_eff a few x 10^5-10^6 K by solving the full radiative transfer equations for both polarization modes in the magnetized hydrogen plasma. The atmospheres directly determine the characteristics of thermal emission from isolated neutron stars. We study the effects of vacuum polarization and bound atoms on the atmosphere structure and spectra. For the lower magnetic field models (B 10^12 G), the spectral features due to neutral atoms lie at extreme UV and very soft X-ray energies and therefore are not likely to be observed. However, the continuum flux is also different from the fully ionized case, especially at lower energies. For the higher magnetic field models, we find that vacuum polarization softens the high energy tail of the thermal spectrum. We show that this depression of continuum flux strongly suppresses not only the proton cyclotron line but also spectral features due to bound species; therefore spectral lines or features in thermal radiation are more difficult to observe when the neutron star magnetic field is > 10^14 G.

The inner region of the accretion disk around a magnetized star is subjected to magnetic torques that induce warping and precession of the disk. These torques arise from interactions between the stellar field and the induced electric currents in the disk. These novel magnetic effects give rise to some ``exotic'' stellar variabilities, and may play an important role in explaining a number of puzzling behaviors related to disk accretion onto magnetic stars, such as mHz QPOs in X-ray pulsars, long-term periodicities of X-ray binaries (including precession of jets), low-Frequency (10-50 Hz) QPO's in low-mass X-ray binaries, and photometric variabilities of T Tauri stars.

Recent studies of the 55 Cancri system suggest the existence of three planets with periods of ~15 days, ~45 days, and ~5500 days (Marcy et al. 2002). The inner two planets are near the 3:1 mean motion commensurability and it is likely that these two planets became trapped in the resonance while farther from the star and migrated together. As the innermost planet begins to dissipate energy through tides the planets break out of the resonance. The final state of the system gives important information about its past history, such as the migration timescale that led to capture.

We have recently shown that in a highly magnetized neutron star atmospheric plasma, vacuum polarization can induce resonant conversion of photon polarization modes via a mechanism analogous to MSW neutrino oscillation. In a recent paper Ozel has dismissed this mode conversion effect as ``mistakes''. Here we explain why our arguments/calculations of this effect are correct.

We construct atmosphere models for strongly magnetized neutron stars with surface fields $B\sim 10^{12}-10^{15}$ G and effective temperatures $T_{\rm eff}\sim 10^6-10^7$ K. The atmospheres directly determine the characteristics of thermal emission from isolated neutron stars, including radio pulsars, soft gamma-ray repeaters, and anomalous X-ray pulsars. In our models, the atmosphere is composed of pure hydrogen or helium and is assumed to be fully ionized. The radiative opacities include free-free absorption and scattering by both electrons and ions computed for the two photon polarization modes in the magnetized electron-ion plasma. Since the radiation emerges from deep layers in the atmosphere with $\rho\ga 10^2$ g/cm$^3$, plasma effects can significantly modify the photon opacities by changing the properties of the polarization modes. In the case where the magnetic field and the surface normal are parallel, we solve the full, angle-dependent, coupled radiative transfer equations for both polarization modes. We also construct atmosphere models for general field orientations based on the diffusion approximation of the transport equations and compare the results with models based on full radiative transport. In general, the emergent thermal radiation exhibits significant deviation from blackbody, with harder spectra at high energies. The spectra also show a broad feature ($\Delta E/\Ebi\sim 1$) around the ion cyclotron resonance $\Ebi=0.63 (Z/A)(B/10^{14}{G})$ keV, where $Z$ and $A$ are the atomic charge and atomic mass of the ion, respectively; this feature is particularly pronounced when $\Ebi\ga 3k\Teff$. Detection of the resonance feature would provide a direct measurement of the surface magnetic fields on magnetars.

It is known that vacuum polarization can modify the photon propagation modes in the atmospheric plasma of a strongly magnetized neutron star. A resonance occurs when the effect of vacuum polarization on the photon modes balances that of the plasma. We show that a photon (with energy $E\go$ a few keV) propagating outward in the atmosphere can convert from one polarization mode into another as it traverses the resonant density, $\rho_{res}\simeq Y_e^{-1}\eta^{-2}(B/10^{14} G)^2(E/1 keV)^2$ g cm$^{-3}$, where $Y_e$ is the electron fraction, and $\eta\sim 1$ is a slowly varying function of the magnetic field $B$. The physics of this mode conversion is analogous to the Mikheyev-Smirnov-Wolfenstein mechanism for neutrino oscillation. Because the two photon modes have vastly different opacities in the atmosphere, this vacuum-induced mode conversion can significantly affect radiative transport and surface emission from strongly magnetized neutron stars.

The inner region of the accretion disk around a magnetized neutron star is subjected to magnetic torques that induce warping and precession of the disk. These torques arise from interactions between the stellar field and the induced electric currents in the disk. We carry out a global analysis of warping/precession modes in a viscous accretion disk, and show that under a wide range of conditions typical of accreting X-ray pulsars, the magnetic warping torque can overcome viscous damping and make the mode grow. The warping/precession modes are concentrated near the inner edge of the disk (at the magnetosphere-disk boundary), and can give rise to variabilities or quasi-periodic oscillations (QPOs) in the X-ray/UV/optical fluxes from X-ray pulsars. We examine the observed properties of mHz QPOs in several systems (such as 4U 1626-67), and suggest that some hitherto unexplained QPOs are the results of magnetically driven disk warping/precession.

An accretion disk around a rotating magnetized star is subjected to magnetic torques which induce disk warping and precession. These torques arise from interactions between the stellar field and the induced surface currents on the disk. Applying these new effects to weakly magnetized ($B\sim 10^7-10^9$G) neutron stars in low-mass X-ray binaries, we study the global hydrodynamical warping/precession modes of the disk under the combined influences of relativistic frame dragging, classical precession due to the oblateness of the neutron star, and the magnetic torques. Under quite general conditions, the magnetic warping torque can overcome the ``Bardeen-Petterson'' viscous damping and makes the modes grow. The modes are confined to the inner region of the disk, and have frequencies close to the sum of the Lense-Thirring frequency, the classical precession frequency, and the magnetically driven precession frequency evaluated at the inner disk radius $r_{in}$. As $\dot M$ increases, the mode frequency is reduced relative to the total precession frequency at $r_{in}$ since the mode is less concentrated around $r_{in}$ due to the increasing viscous stress associated with the large $\dot M$. Because of this, and because the magnetically driven precession is retrograde and depends strongly on $\dot M$, the mode frequency can have a non-monotonic dependence on $\dot M$. This may account for several observed features of low-frequency (10--60 Hz) quasi-periodic oscillations (LFQPOs) in low-mass X-ray binaries which are otherwise difficult to explain, such as the flattening/turnover in the LFQPO frequency--$\dot M$ correlation or in the LFQPO frequency--kHz QPO frequency correlation (e.g., as seen clearly in GX 17+2).

The properties of matter are significantly modified by strong magnetic fields, $B>>2.35\times 10^9$ Gauss ($1 G =10^{-4} Tesla$), as are typically found on the surfaces of neutron stars. In such strong magnetic fields, the Coulomb force on an electron acts as a small perturbation compared to the magnetic force. The strong field condition can also be mimicked in laboratory semiconductors. Because of the strong magnetic confinement of electrons perpendicular to the field, atoms attain a much greater binding energy compared to the zero-field case, and various other bound states become possible, including molecular chains and three-dimensional condensed matter. This article reviews the electronic structure of atoms, molecules and bulk matter, as well as the thermodynamic properties of dense plasma, in strong magnetic fields, $10^9G << B < 10^{16}G$. The focus is on the basic physical pictures and approximate scaling relations, although various theoretical approaches and numerical results are also discussed. For the neutron star surface composed of light elements such as hydrogen or helium, the outermost layer constitutes a nondegenerate, partially ionized Coulomb plasma if $B<<10^{14}G$, and may be in the form of a condensed liquid if the magnetic field is stronger (and temperature $<10^6$ K). For the iron surface, the outermost layer of the neutron star can be in a gaseous or a condensed phase depending on the cohesive property of the iron condensate.

The secular instability and nonlinear evolution of the m=2 f-mode (bar-mode) driven by gravitational radiation reaction in a rapidly rotating, newly formed neutron star are reviewed. There are two types of rotating bars which generate quite different gravitational waveforms: those with large internal rotation relative to the bar figure (Dedekind-like bars) have GW frequency sweeping downward during the evolution, and those with small internal rotation (Jacobi-like) have GW frequency sweeping upward. Various sources of viscosity (which affects the instability) in hot nuclear matter are reexamined, and the possible effect of star-disk coupling on the bar-mode instability is also discussed.

Observational advances over the last decade have left little doubt that neutron stars received a large kick velocity (of order a few hundred to a thousand km/s) at birth. The physical origin of the kicks and the related supernova asymmetry is one of the central unsolved mysteries of supernova research. We review the physics of different kick mechanisms, including hydrodynamically driven, neutrino -- magnetic field driven, and electromagnetically driven kicks. The viabilities of the different kick mechanisms are directly related to the other key parameters characterizing nascent neutron stars, such as the initial magnetic field and the initial spin. Recent observational constraints on kick mechanisms are also discussed.

We study implications for the apparent alignment of the spin axes, proper-motions, and polarization vectors of the Crab and Vela pulsars. The spin axes are deduced from recent Chandra X-ray Observatory images that reveal jets and nebular structure having definite symmetry axes. The alignments indicate these pulsars were born either in isolation or with negligible velocity contributions from binary motions. We examine the effects of rotation and the conditions under which spin-kick alignment is produced for various models of neutron star kicks. If the kick is generated when the neutron star first forms by asymmetric mass ejection or/and neutrino emission, then the alignment requires that the protoneutron star possesses an original spin with period $P_s$ much less than the kick timescale, thus spin-averaging the kick forces. The kick timescale ranges from 100 ms to 10 s depending on whether the kick is hydrodynamically driven or neutrino-magnetic field driven. For hydrodynamical models, spin-kick alignment further requires the rotation period of an asymmetry pattern at the radius near shock breakout (>100 km) to be much less than ~100 ms; this is difficult to satisfy unless rotation plays a dynamically important role in the core collapse and explosion ($P_s\lo 1$ ms). Aligned kick and spin vectors are inherent to the slow process of asymmetric electromagnetic radiation from an off-centered magnetic dipole. We reassess the viability of this effect, correcting a factor of 4 error in Harrison and Tademaru's calculation that increases the size of the effect. To produce a kick velocity of order a few hundred km/s requires that the neutron star be born with an initial spin close to 1 ms and that spindown due to r-mode driven gravitational radiation be inefficient compared to standard magnetic braking.

Recent work has shown that a young, rapidly rotating neutron star loses angular momentum to gravitational waves generated by unstable r-mode oscillations. We study the spin evolution of a young, magnetic neutron star including both the effects of gravitational radiation and magnetic braking (modeled as magnetic dipole radiation). Our phenomenological description of nonlinear r-modes is similar to, but distinct from, that of Owen et al. (1998) in that our treatment is consistent with the principle of adiabatic invariance in the limit when direct driving and damping of the mode are absent. We show that, while magnetic braking tends to increase the r-mode amplitude by spinning down the neutron star, it nevertheless reduces the efficiency of gravitational wave emission from the star. For B >= 10^14 (\nus/300 Hz)^2 G, where \nus is the spin frequency, the spindown rate and the gravitational waveforms are significantly modified by the effect of magnetic braking. We also estimate the growth rate of the r-mode due to electromagnetic (fast magnetosonic) wave emission and due to Alfven wave emission in the neutron star magnetosphere. The Alfven wave driving of the r-mode becomes more important than the gravitational radiation driving when B >= 10^13 (\nus/150 Hz)^3 G; the electromagnetic wave driving of the r-mode is much weaker. Finally, we study the properties of local Rossby-Alfven waves inside the neutron star and show that the fractional change of the r-mode frequency due to the magnetic field is of order 0.5 (B/10^16 G)^2 (\nus/100 Hz)^-2.

Self-similar solutions provide good descriptions for the gravitational collapse of spherical clouds or stars when the gas obeys a polytropic equation of state, $p=K\rho^\gamma$ (with $\gamma\le 4/3$). We study the behaviors of nonradial perturbations in the similarity solutions of Larson, Penston and Yahil, which describe the evolution of the collapsing cloud prior to core formation. Our global stability analysis reveals the existence of unstable bar-modes ($l=2$) when $\gamma\le 1.09$. In particular, for the collapse of isothermal spheres, which applies to the early stages of star formation, the $l=2$ density perturbation relative to the background, $\delta\rho({\bf r},t)/\rho(r,t)$, increases as $(t_0-t)^{-0.352}\propto \rho_c(t)^{0.176}$, where $t_0$ denotes the epoch of core formation, and $\rho_c(t)$ is the cloud central density. Thus, the isothermal cloud tends to evolve into an ellipsoidal shape (prolate bar or oblate disk, depending on initial conditions) as the collapse proceeds. In the context of Type II supernovae, core collapse is described by the $\gamma\simeq 1.3$ equation of state, and our analysis indicates that there is no growing mode (with density perturbation) in the collapsing core before the proto-neutron star forms, although nonradial perturbations can grow during the subsequent accretion of the outer core and envelope onto the neutron star. We also carry out a global stability analysis for the self-similar expansion-wave solution found by Shu, which describes the post-collapse accretion (``inside-out'' collapse) of isothermal gas onto a protostar. We show that this solution is unstable to perturbations of all $l$'s, although the growth rates are unknown.

It is no longer necessary to `sell' the idea of pulsar kicks, the notion that neutron stars receive a large velocity (a few hundred to a thousand km s$^{-1}$) at birth. However, the origin of the kicks remains mysterious. We review the physics of different kick mechanisms, including hydrodynamically driven, neutrino and magnetically driven kicks.

The inner region of the accretion disk onto a rotating magnetized central star (neutron star, white dwarf or T Tauri star) is subjected to magnetic torques which induce warping and precession of the disk. The origin of these torques lies in the interaction between the (induced) surface current on the disk and the horizontal magnetic field (parallel to the disk) produced by the inclined magnetic dipole. Under quite general conditions, there exists a magnetic warping instability in which the magnetic torque drives the disk plane away from the equatorial plane of the star toward a state where the disk normal vector is perpendicular to the spin axis. Viscous stress tends to suppress the warping instability at large radii, but the magnetic torque always dominates as the disk approaches the magnetosphere boundary. The magnetic torque also drives the tilted inner disk into retrograde precession around the stellar spin axis. Moreover, resonant magnetic forcing on the disk can occur which may affect the dynamics of the disk. The magnetically driven warping instability and precession may be related to a number observational puzzles, including: (1) Spin evolution (torque reversal) of accreting X-ray pulsars; (2) Quasi-periodic oscillations in low-mass X-ray binaries; (3) Super-orbital periods in X-ray binaries; (4) Photometric period variations of T Tauri stars.

When a self-gravitating spherical gas cloud collapses or accretes onto a central mass, the inner region of the cloud develops a density profile $\rho\propto r^{-3/2}$ and the velocity approaches free-fall. We show that in this region, nonspherical perturbations grow with decreasing radius. In the linear regime, the tangential velocity perturbation increases as $r^{-1}$, while the Lagrangian density perturbation, $\Delta\rho/\rho$, grows as $r^{-1/2}$. Faster growth occurs if the central collapsed object maintains a finite multiple moment, in which case $\Delta\rho/\rho$ increases as $r^{-l}$, where $l$ specifies the angular degree of the perturbation. These scaling relations are different from those obtained for the collapse of a homogeneous cloud. Our numerical calculations indicate that nonspherical perturbations are damped in the subsonic region, and that they grow and approach the asymptotic scalings in the supersonic region. The implications of our results to asymmetric supernova collapse and to black hole accretion are briefly discussed.

We present a semi-analytic study of the equilibrium models of close binary systems containing a fluid star (mass $m$ and radius $R_0$) and a Kerr black hole (mass $M$) in circular orbit. We consider the limit $M\gg m$ where spacetime is described by the Kerr metric. The tidally deformed star is approximated by an ellipsoid, and satisfies the polytropic equation of state. The models also include fluid motion in the stellar interior, allowing binary models with nonsynchronized stellar spin (as expected for coalescing neutron star-black hole binaries) to be constructed. Tidal disruption occurs at orbital radius $r_{\rm tide}\sim R_0(M/m)^{1/3}$, but the dimensionless ratio $\hat r_{\rm tide}=r_{\rm tide}/[R_0(M/m)^{1/3}]$ depends on the spin parameter of the black hole as well as on the equation of state and the internal rotation of the star. We find that the general relativistic tidal field disrupts the star at a larger $\hat r_{\rm tide}$ than the Newtonian tide; the difference is particularly prominent if the disruption occurs in the vicinity of the black hole's horizon. In general, $\hat r_{\rm tide}$ is smaller for a (prograde rotating) Kerr black hole than for a Schwarzschild black hole. We apply our results to coalescing black hole-neutron star and black hole-white dwarf binaries. The tidal disruption limit is important for characterizing the expected gravitational wave signals and is relevant for determining the energetics of gamma ray bursts which may result from such disruption.

Gravitational radiation tends to drive gravity modes in rotating neutron stars unstable. For an inviscid star, the instability sets in when the rotation frequency is about 0.7 times the corresponding mode frequency of the nonrotating star. Neutron stars with spin frequencies $\go 100$ Hz are susceptible to this instability, with growth time of order years. However, it is likely that viscous dissipation suppresses the instability except for a narrow range of temperatures around $10^9$ K. We also show that the viscosity driven instability of g-modes is absent.

It has recently been suggested that the maximum observed quasi-periodic oscillation (QPO) frequencies, $\nu_{max}$, for several low-mass X-ray binaries, particularly 4U 1820-30, correspond to the orbital frequency at the inner-most stable orbit of the accretion disk. This would imply that the neutron stars in these systems have masses $\go 2~M_\odot$, considerably larger than any well-measured neutron star mass. We suggest that the levelling off of $\nu_{QPO}$ may be also understood in terms of a steepening magnetic field which, although possibly dipolar at the stellar surface, is altered substantially by disk accretion, and presents a ``wall'' to the accretion flow that may be outside the innermost stable orbit. General relativistic effects add to the flattening of the $\nu_{QPO}-\mdot$ relation at frequencies below the Kepler frequency at the innermost stable orbit. We offer two other possible ways to reconcile the low value of $\nu_{max}$ ($\approx 1060$ Hz for 4U 1820-30) with a moderate neutron star mass, $\approx 1.4\msun$: at sufficiently large $\mdot$, either (i) the disk terminates in a very thin boundary layer near the neutron star surface, or (ii) $\nu_{QPO}$ is not the orbital frequency right at the inner edge of the disk, but rather at a somewhat larger radius, where the emissivity of the disk peaks.

In a coalescing neutron star-neutron star (NS-NS) or neutron star-black hole (NS-BH) binary, oscillation modes of the NS can be resonantly excited by the companion during the final minutes of the inspiral. The resonant energy transfer between the orbit and NS speeds up or slows down the inspiral and induces a phase change in the emitted gravitational waves from the binary. A tidal resonance, (jk,m), occurs when the mode frequency equals m times the orbital frequency. For f-mode resonance to occur before coalescence, the NS must have rapid rotation, with spin frequency nu_s>710Hz for (22,2)-resonance and nu_s>570Hz for (33,3)-resonance (M=1.4Mo and R=10km; however, for R=15km, these spin frequencies become 330Hz and 250Hz). Because of strong tidal coupling, f-mode resonances induce a large change in the number of orbital cycles, N_orb, with maximum N_orb~10-1000 for (22,2)-resonance and N_orb~1 for (33,3)- resonance. Such resonant effects, if present, must be included in constructing waveform templates used in searching for gravitational wave signals. Higher order f-mode resonances can occur at slower rotation rates, but N_orb<0.1. For the dominant g-mode (22,2)-resonance, even modest rotation (nu_s<100Hz) can enhance the resonant effect on the orbit by shifting resonance to a smaller orbital frequency. However, because of the weak coupling, N_orb lies in the range 10^-3-10^-2 (depending on the NS EOS) and is probably negligible for the purpose of detecting gravitational waves. R-mode resonances require misaligned spin-orbit inclinations, and the dominant resonances correspond to (22,3) and (22,1). Since tidal coupling depends strongly on rotation rate, N_orb<10^-2(R/10km)^(10)(M/1.4Mo)^(-20/3) is negligible for canonical NS parameters but can be appreciable if the NS radius is larger.

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.

The original paper describes neutrino cooling calculations in the presence of asymmetric neutrino opacities due to parity violation in strong magnetic fields. In early April, through the work of P. Arras and D. Lai at Cornell (see astro-ph/9806285) we became aware that the paper contains a serious flaw: although the neutrino scattering opacities depend asymmetrically on the initial and final neutrino momenta, detailed balance requires that there be no cumulative effect from multiple scatterings in the bulk interior of the proto-neutron star where local thermodynamic equilibrium applies to a good approximation. This paper was therefore withdrawn from ApJ in April. Also in April, we submitted an erratum to ApJ Lett. pointing out this error. (The erratum will appear in the July 1, 1998 issue of ApJ Lett.)

In proto-neutron stars with strong magnetic fields, the cross section for $\nu_e$ ($\bar\nu_e$) absorption on neutrons (protons) depends on the local magnetic field strength due to the quantization of energy levels for the $e^-$ ($e^+$) produced in the final state. If the neutron star possesses an asymmetric magnetic field topology in the sense that the magnitude of magnetic field in the north pole is different from that in the south pole, then asymmetric neutrino emission may be generated. We calculate the absorption cross sections of $\nue$ and $\bnue$ in strong magnetic fields as a function of the neutrino energy. These cross sections exhibit oscillatory behaviors which occur because new Landau levels for the $e^-$ ($e^+$) become accessible as the neutrino energy increases. By evaluating the appropriately averaged neutrino opacities, we demonstrate that the change in the local neutrino flux due to the modified opacities is rather small. To generate appreciable kick velocity ($\sim 300$ km~s$^{-1}$) to the newly-formed neutron star, the difference in the field strengths at the two opposite poles of the star must be at least $10^{16}$~G. We also consider the magnetic field effect on the spectral neutrino energy fluxes. The oscillatory features in the absorption opacities give rise to modulations in the emergent spectra of $\nu_e$ and $\bar\nu_e$.

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}$.

The inner regions of accretion disks of weakly magnetized neutron stars are affected by general relativity and stellar magnetic fields. Even for field strengths sufficiently small so that there is no well-defined magnetosphere surrounding the neutron star, there is still a region in the disk where magnetic field stress plays an important dynamical role. We construct magnetic slim disk models appropriate for neutron stars in low-mass X-ray binaries (LMXBs) which incorporate both effects (GR and magnetic fields). The B-field--disk interaction is treated in a phenomenological manner, allowing for both closed and open field configurations. We show that even for surface magnetic fields as weak as $10^7-10^8$ G, the sonic point of the accretion flow can be significantly modified from the pure GR value (near $6M$). We derive an approximate analytical expression for the sonic radius and show that it mainly depends on the surface field strength $B_0$ and mass accretion rate $\dot M$ through the ratio $b^2\propto B_0^2/\dot M$. The sonic radius thus obtained approaches the usual Alfven radius for high $b^2$, and asymptotes to $6M$ as $b^2\to 0$. We therefore suggest that for neutron stars in LMXBs, the distinction between the disk sonic radius and the magnetospheric radius may not exist. We apply our theoretical results to the kHz QPOs observed in the X-ray fluxes of LMXBs. If these QPOs are associated with the orbital frequency at the inner radius of the disk, then the QPO frequencies and their correlation with mass accretion rate can provide useful constraints on the nature of the magnetic field -- disk interactions as well as on the structure of magnetic fields in LMXBs. Current observational data may suggest that the magnetic fields in LMXBs have complex topology.

In proto-neutron stars with strong magnetic fields, the neutrino-nucleon scattering/absorption cross sections depend on the direction of neutrino momentum with respect to the magnetic field axis, a manifestation of parity violation in weak interactions. We study the deleptonization and thermal cooling (via neutrino emission) of proto-neutron stars in the presence of such asymmetric neutrino opacities. Significant asymmetry in neutrino emission is obtained due to multiple neutrino-nucleon scatterings. For an ordered magnetic field threading the neutron star interior, the fractional asymmetry in neutrino emission is about $0.006 (B/10^{14}G)$, corresponding to a pulsar kick velocity of about $200 (B/10^{14}G)$ km/s for a total radiated neutrino energy of $3\times 10^{53}$ erg.

Recent timing observations have revealed the presence of orbital precession due to spin-orbit coupling and rapid orbital decay due to dynamical tidal interaction in the PSR J0045-7319/B-star binary system. They can be used to put concrete constraints on the age, initial spin and velocity of the neutron star.

We study the effect of rotation on the excitation of internal oscillation modes of a star by the external gravitational potential of its companion. Unlike the nonrotating case, there are difficulties with the usual mode decomposition for rotating stars because of the asymmetry between modes propagating in the direction of rotation and those propagating opposite to it. For an eccentric binary system, we derive general expressions for the energy transfer $\Delta E_s$ and the corresponding angular momentum transfer $\Delta J_s$ in a periastron passage when there is no initial oscillation present in the star. Except when a nearly precise orbital resonance occurs (i.e., the mode frequency equals multiple of the orbital frequency), $\Delta E_s$ is very close to the steady-state mode energy in the tide in the presence of dissipation. It is shown that stellar rotation can change the strength of dynamical tide significantly. In particular, retrograde rotation with respect to the orbit increases the energy transfer by bringing lower-order g-modes (or f-mode for convective stars), which couple more strongly to the tidal potential, into closer resonances with the orbital motion of the companion. We apply our general formalism to the problems of tidal capture binary formation and the orbital evolution of PSR J0045-7319/B-star binary. Stellar rotation changes the critical impact parameter for binary capture. Although the enhancement (by retrograde rotation) in the capture cross section is at most $\sim 20%$, the probability that the captured system survives disruption/merging and therefore becomes a binary can be significantly larger. It is found that in order to explain the observed rapid orbital decay of the PSR J0045-7319 binary system, retrograde rotation in the B-star is required.

The outermost layers of some neutron stars are likely to be dominated by hydrogen, as a result of fast gravitational settling of heavier elements. These layers directly mediate thermal radiation from the stars, and determine the characteristics of X-ray/EUV spectra. For a neutron star with surface temperature $T\lo 10^6$ K and magnetic field $B\go 10^{12}$ G, various forms of hydrogen can be present in the envelope, including atom, poly-molecules, and condensed metal. We study the physical properties of different hydrogen phases on the surface of a strongly magnetized neutron star for a wide range of field strength $B$ and surface temperature $T$. Depending on the values of $B$ and $T$, the outer envelope can be either in a nondegenerate gaseous phase or in a degenerate metallic phase. For $T\go 10^5$ K and moderately strong magnetic field, $B\lo 10^{13}$ G, the envelope is nondegenerate and the surface material gradually transforms into a degenerate Coulomb plasma as density increases. For higher field strength, $B>> 10^{13}$ G, there exists a first-order phase transition from the nondegenerate gaseous phase to the condensed metallic phase. The column density of saturated vapor above the metallic hydrogen decreases rapidly as the magnetic field increases or/and temperature decreases. Thus the thermal radiation can directly emerge from the degenerate metallic hydrogen surface. The characteristics of surface X-ray/EUV emission for different phases are discussed. A separate study concerning the possibility of magnetic field induced nuclear fusion of hydrogen on the neutron star surface is also presented.

We study how the neutron-star equation of state affects the onset of the dynamical instability in the equations of motion for inspiraling neutron-star binaries near coalescence. A combination of relativistic effects and Newtonian tidal effects cause the stars to begin their final, rapid, and dynamically-unstable plunge to merger when the stars are still well separated and the orbital frequency is $\approx$ 500 cycles/sec (i.e. the gravitational wave frequency is approximately 1000 Hz). The orbital frequency at which the dynamical instability occurs (i.e. the orbital frequency at the innermost stable circular orbit) shows modest sensitivity to the neutron-star equation of state (particularly the mass-radius ratio, $M/R_o$, of the stars). This suggests that information about the equation of state of nuclear matter is encoded in the gravitational waves emitted just prior to the merger.

What happens to a neutron star or white dwarf near its maximum mass limit when it is brought into a close binary orbit with a companion? Such situation may occur in the progenitors of Type Ia supernovae and in coalescing neutron star binaries. Using an energy variational principle, we show that tidal field reduces the central density of the compact object, making it more stable against radial collapse. For a cold white dwarf, the tidal field increases the maximum stable mass only slightly, but can actually lower the maximum central density by as much as $30\%$. Thus a white dwarf in a close binary may be more susceptible to general relativistic instability than the instability associated with electron capture and pycronuclear reaction (depending on the white dwarf composition). We analyse the radial stability of neutron star using post-Newtonian approximation with an ideal degenerate neutron gas equation of state. The tidal stablization effect implies that the neutron star in coalescing neutron star-neutron star or neutron star-black hole binaries does not collapse prior to merger or tidal disruption.

Recent timing observations of PSR J0045-7319 reveal that the neutron star/B star binary orbit is decaying on a time scale of $|\Porb/\dot\Porb|=0.5$ Myr, shorter than the characteristic age ($\tau_c=3$ Myr) of the pulsar (Kaspi et al.~1996a). We study mechanisms for the orbital decay. The standard weak friction theory based on static tide requires far too short a viscous time to explain the observed $\dot\Porb$. We show that dynamical tidal excitation of g-modes in the B star can be responsible for the orbital decay. However, to explain the observed short decay timescale, the B star must have some significant retrograde rotation with respect to the orbit --- The retrograde rotation brings lower-order g-modes, which couple much more strongly to the tidal potential, into closer ``resonances'' with the orbital motion, thus significantly enhancing the dynamical tide. A much less likely possibility is that the g-mode damping time is much shorter than the ordinary radiative damping time. The observed orbital decay timescale combined with a generic orbital evolution model based on dynamical tide can be used as a ``timer'', giving an upper limit of $1.4$ Myr for the age of the binary system since the neutron star formation. Thus the characteristic age of the pulsar is not a good age indicator. Assuming standard magnetic dipole braking for the pulsar and no significant magnetic field decay on a timescale $\lo 1$ Myr, the upper limit for the age implies that the initial spin of the neutron star at birth was close to its current value.

We study the effects of finite proton mass on the energy levels of hydrogen atoms moving transverse to a superstrong magnetic field $B$ with generalized pseudomomentum $K_\perp$. Field strengths of order $B\sim 10^{12}$ Gauss are typically found on the surfaces of neutron stars, but we also study the regime $B\go B_{crit}= 4.23\times 10^{13}$ Gauss, where the Landau excitation energy of the proton is large. We adopt two different approaches to the two-body problem in strong magnetic field, and obtain an approximate, but complete solution of the atomic energy as a function of $B$ and $K_\perp$. We show that, for $B>>B_{crit}$, there is an orthogonal set of bound states, which do not have any Landau excitation contribution in their energies. The states with very large $K_\perp$ have small binding energy and small transverse velocity, but are nevertheless distinct from the fully ionized states. The final results for the excitation energies are given in the form of analytical fitting formulae. The generalized Saha equation for the ionization-recombination equilibrium of hydrogen gas in the presence of a superstrong magnetic field is then derived. Although the maximum transverse velocity of a bound atom decreases as $B$ increases, the statistical weight due to transverse motion is actually increased by the strong magnetic field. For astrophysically interesting case of relatively low density and temperature, we obtain analytic approximations for the partition functions. The highly excited bound states have a smaller statistical weight than the fully ionized component.

We study the energy levels of H$_2$ molecules in a superstrong magnetic field ($B\go 10^{12}$ G), typically found on the surfaces of neutron stars. The interatomic interaction potentials are calculated by a Hartree-Fock method with multi-configurations assuming electrons are in the ground Landau state. Both the aligned configurations and arbitrary orientations of the molecular axis with respect to the magnetic field axis are considered. Different types of molecular excitations are then studied: electronic excitations, aligned (along the magnetic axis) vibrational excitations, transverse vibrational excitations (a constrained rotation of the molecular axis around the magnetic field line). Similar results for the molecular ion H$_2^+$ are also obtained and compared with previous variational calculations. Both numerical results and analytical fitting formulae are given for a wide range of field strengths. In contrast to the zero-field case, it is found that the transverse vibrational excitation energies can be larger than the aligned vibration excitation, and they both can be comparable or larger than the electronic excitations. For $B\go B_{crit}=4.23\times 10^{13}$ G, the Landau energy of proton is appreciable and there is some controversy regarding the dissociation energy of H$_2$. We show that H$_2$ is bound even for $B>>B_{crit}$ and that neither proton has a Landau excitation in the ground molecular state.

The recent claim that in a strong magnetic field hydrogenlike gas (e.g., excitons in certain semiconductors, neutron star surface layers) becomes superfluid is refuted. Molecules form by strong covalent bond along the magnetic field axis, which prohibits Bose-Einstein condensation.

The spin-induced quadrupole moment of a rapidly rotating star changes the orbital dynamics in a binary system, giving rise to advance (or regression) of periastron and precession of the orbital plane. We show that these effects are important in the recently discovered radio pulsar/main sequence star binary system PSR J0045$-$7319, and can reliably account for the observed peculiar timing residuals. Precise measurements of the apsidal motion and orbital plane precession can yield valuable information on the internal structure and rotation of the star. The detection of orbital precession implies that the spin of the companion star is not aligned with the orbital angular momentum, and suggests that the supernova gave the pulsar a kick out of the original orbital plane. Tidal excitation of g-mode oscillations in the PSR J0045$-$7319 system induces an orbital period change of order $|\Delta P_{\rm orb}/P_{\rm orb}|\sim 10^{-6}$ at each periastron passage, but the secular trend depends on the radiative damping time of the g-modes. We also discuss the spin-orbit coupling effects for the accreting X-ray pulsars and the other known radio pulsar/main sequence binary, PSR B1259$-$63.

We employ an approximate treatment of dissipative hydrodynamics in three dimensions to study the coalescence of binary neutron stars driven by the emission of gravitational waves. The stars are modeled as compressible ellipsoids obeying a polytropic equation of state; all internal fluid velocities are assumed to be linear functions of the coordinates. The hydrodynamic equations then reduce to a set of coupled ordinary differential equations for the evolution of the principal axes of the ellipsoids, the internal velocity parameters and the binary orbital parameters. Gravitational radiation reaction and viscous dissipation are both incorporated. We set up exact initial binary equilibrium configurations and follow the transition from the quasi-static, secular decay of the orbit at large separation to the rapid dynamical evolution of the configurations just prior to contact. A hydrodynamical instability resulting from tidal interactions significantly accelerates the coalescence at small separation, leading to appreciable radial infall velocity and tidal lag angles near contact. This behavior is reflected in the gravitational waveforms and may be observable by gravitational wave detectors under construction.

We study the secular evolution and gravitational wave signature of a newly-formed, rapidly rotating neutron star. The neutron star may arise from core collapse in a massive star or from the accretion-induced collapse of a white dwarf. After a brief dynamical phase, the nascent neutron star settles into an axisymmetric, secularly unstable equilibrium state. Gravitational radiation drives the star to a nonaxisymmetric, stationary equilibrium configuration via the bar-mode instability. The emitted quasi-periodic gravitational waves have a unique signature: the wave frequency sweeps downward from a few hundred Hertz to zero, while the wave amplitude increases from zero to a maximum and then decays back to zero. Such a wave signal could be detected by broad-band gravitational wave interferometers currently being constructed. We also characterize two other types of gravitational wave signals that could arise in principle from a rapidly rotating, secularly unstable neutron star: a high-frequency ($f\go 1000$ Hz) wave which increases the pattern-speed of the star, and a wave that actually increases the angular momentum of the star.

Tidal interaction in a coalescing neutron star binary can resonantly excite the g-mode oscillations of the neutron star when the frequency of the tidal driving force equals the intrinsic g-mode frequencies. We study the g-mode oscillations of cold neutron stars using recent microscopic nuclear equations of state, where we determine self-consistently the sound speed and Brunt-V\"ais\"al\"a frequency in the nuclear liquid core. The properties of the g-modes associated with the stable stratification of the core depend sensitively on the pressure-density relation as well as the symmetry energy of the dense nuclear matter. The frequencies of the first ten g-modes lie approximately in the range of $10-100$ Hz. Resonant excitations of these g-modes during the last few minutes of the binary coalescence result in energy transfer and angular momentum transfer from the binary orbit to the neutron star. The angular momentum transfer is possible because a dynamical tidal lag develops even in the absence of fluid viscosity. However, since the coupling between the g-mode and the tidal potential is rather weak, the amount of energy transfer during a resonance and the induced orbital phase error are very small. Resonant excitations of the g-modes play an important role in tidal heating of binary neutron stars. Without the resonances, viscous dissipation is effective only when the stars are close to contact. The resonant oscillations result in dissipation at much larger orbital separation. The actual amount of tidal heating depends on the viscosity of the neutron star. Using the microscopic viscosity, we find that the binary neutron stars are heated to a temperature $\sim 10^8$ K before they come into contact.

We develop a new formalism to study the dynamics of fluid polytropes in three dimensions. The stars are modeled as compressible ellipsoids and the hydrodynamic equations are reduced to a set of ordinary differential equations for the evolution of the principal axes and other global quantities. Both viscous dissipation and the gravitational radiation reaction are incorporated. We establish the validity of our approximations and demonstrate the simplicity and power of the method by rederiving a number of known results concerning the stability and dynamical oscillations of rapidly rotating polytropes. In particular, we present a generalization to compressible fluids of Chandrasekhar's classical results for the secular and dynamical instabilities of incompressible Maclaurin spheroids. We also present several applications of our method to astrophysical problems of great current interest such as the tidal disruption of a star by a massive black hole, the coalescence of compact binaries driven by the emission of gravitational waves, and the development of instabilities in close binary systems.

We present a new analytic study of the equilibrium and stability properties of close binary systems containing polytropic components. Our method is based on the use of ellipsoidal trial functions in an energy variational principle. We consider both synchronized and nonsynchronized systems, constructing the compressible generalizations of the classical Darwin and Darwin-Riemann configurations. Our method can be applied to a wide variety of binary models where the stellar masses, radii, spins, entropies, and polytropic indices are all allowed to vary over wide ranges and independently for each component. We find that both secular and dynamical instabilities can develop before a Roche limit or contact is reached along a sequence of models with decreasing binary separation. High incompressibility always makes a given binary system more susceptible to these instabilities, but the dependence on the mass ratio is more complicated. As simple applications, we construct models of double degenerate systems and of low-mass main-sequence-star binaries. We also discuss the orbital evolution of close binary systems under the combined influence of fluid viscosity and secular angular momentum losses from processes like gravitational radiation. We show that the existence of global fluid instabilities can have a profound effect on the terminal evolution of coalescing binaries. The validity of our analytic solutions is examined by means of detailed comparisons with the results of recent numerical fluid calculations in three dimensions.

We study the importance of hydrodynamic effects on the evolution of coalescing binary neutron stars. Using an approximate energy functional constructed from equilibrium solutions for polytropic binary configurations, we incorporate hydrodynamic effects into the calculation of the orbital decay driven by gravitational wave emission. In particular, we follow the transition between the quasi-static, secular decay of the orbit at large separation and the rapid dynamical evolution of configurations approaching contact. We show that a purely Newtonian hydrodynamic instability can significantly accelerate the coalescence at small separation. Such an instability occurs in all close binary configurations containing sufficiently incompressible stars. Calculations are performed for various neutron star masses, radii, spins, and effective polytropic indices. Typically, we find that the radial infall velocity just prior to contact is about 10\% of the tangential orbital velocity. Post-Newtonian effects can move the stability limit to a larger binary separation, and may induce an even larger radial velocity. We also consider the possibility of mass transfer from one neutron star to the other. We show that stable mass transfer is impossible except when the mass of one of the components is very small (less than about 0.4 solar mass) and the viscosity is high enough to maintain corotation.