Source author record

Dong Lai

Dong Lai appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

85works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

85 published item(s)

preprint2022arXiv

Dynamical Instability in Multi-Orbiter Systems with Gas Friction

Closely-packed multi-planet systems are known to experience dynamical instability if the spacings between the planets are too small. Such instability can be tempered by the frictional forces acting on the planets from gaseous discs. A similar situation applies to stellar-mass black holes embedded in AGN discs around supermassive black holes. In this paper, we use $N$-body integrations to evaluate how the frictional damping of orbital eccentricity affects the growth of dynamical instability for a wide range of planetary spacing and planet-to-star mass ratios. We find that the stability of a system depends on the damping timescale $τ$ relative to the zero-friction instability growth timescale $t_{\rm inst}$. In a two-planet system, the frictional damping can stabilise the dynamical evolution if $t_{\rm inst}\gtrsimτ$. With three planets, $t_{\rm inst} \gtrsim 10τ- 100τ$ is needed for stabilisation. When the separations between the planetary orbits are sufficiently small, $t_{\rm inst}$ can be less than the synodic period between the planets, which makes frictional stabilisation unlikely to occur. As the orbital spacing increases, the instability timescale tends to grow exponentially on average, but it can vary by a few orders of magnitude depending on the initial orbital phases of the planets. In general, the stable region (at large orbital spacings) and unstable region (at small orbital spacings) are separated by a transition zone, in which the (in)stability of the system is not guaranteed. We also devise a linear map to analyse the dynamical instability of the "planet + test-mass" system, and we find qualitatively similar results to the $N$-body simulations.

preprint2022arXiv

Evolution of Stellar Orbits Around Merging Massive Black-Hole Binary

We study the long-term orbital evolution of stars around a merging massive or supermassive black-hole (BH) binary, taking into account the general relativistic effect induced by the BH spin. When the BH spin is significant compared to and misaligned with the binary orbital angular momentum, the orbital axis ($\hat{\mathbf{l}}$) of the circumbinary star can undergo significant evolution during the binary orbital decay driven by gravitational radiation. Including the spin effect of the primary (more massive) BH, we find that starting from nearly coplanar orbital orientations, the orbital axes $\hat{\mathbf{l}}$ of circumbinary stars preferentially evolve towards the spin direction after the merger of the BH binary, regardless of the initial BH spin orientation. Such alignment phenomenon, i.e., small final misalignment angle between $\hat{\mathbf{l}}$ and the spin axis of the remanent BH $\hat{\mathbf{S}}$, can be understood analytically using the principle of adiabatic invariance. For the BH binaries with extremely mass ratio ($m_2/m_1\lesssim0.01$), $\hat{\mathbf{l}}$ may experience more complicated evolution as adiabatic invariance breaks down, but the trend of alignment still works reasonably well when the initial binary spin-orbit angle is relatively small. Our result suggests that the correlation between the orientations of stellar orbits and the spin axis of the central BH could provide a potential signature of the merger history of the massive BH.

preprint2022arXiv

Long-term Evolution of Tightly-Packed Stellar Black Holes in AGN Disks: Formation of Merging Black-Hole Binaries via Close Encounters

We study the long-term evolution of two or more stellar black holes (BHs) on initially separated but unstable circular orbits around a supermassive BH (SMBH). Such a close-packed orbital configuration can naturally arise from BH migrations in the AGN disk. Dynamical instability of the orbits leads to recurring close encounters between two BHs, during which the BH separation $r_{\rm p}$ becomes less than the Hill radius $R_{\rm H}$. In the rare very close encounters (with $r_{\rm p}$ several orders of magnitude less than $R_{\rm H}$), a tight merging BH binary can form with the help of gravitational wave emission. We use $N$-body simulations to study the time evolution of close encounters of various degrees of "closeness" and the property of the resulting binary BH mergers. For a typical "SMBH + 2 BHs" system, the averaged cumulative number of close encounters (with $r_{\rm p} \lesssim R_{\rm H}$) scales approximately as $\propto t^{0.5}$. The minimum encounter separation $r_{\rm p}$ follows a linear cumulative distribution $P(<r_{\rm p}) \propto r_{\rm p}$ for $r_{\rm p} \ll R_{\rm H}$. From these, we obtain a semi-analytical expression for the averaged rate of binary captures that lead to BH mergers. Our results suggest that close-packed BHs in AGN disks may take a long time ($\gtrsim 10^7$ orbits around the SMBH) to experience a sufficiently close encounter and form a bound binary, although this time can be shorter if the initial BH orbits are highly aligned. The BH binary mergers produced in this scenario always have high eccentricities when entering the LIGO band, and have a broad distribution of orbital inclinations relative to the original AGN disk. We also explore the effects of the gas disk and find that simple gas drags on the BHs do not necessarily lead to an enhanced BH binary capture rate.

preprint2022arXiv

Probing the Spins of Supermassive Black Holes with Gravitational Waves from Surrounding Compact Binaries

Merging compact black-hole (BH) binaries are likely to exist in the nuclear star clusters around supermassive BHs (SMBHs), such as Sgr A$^\ast$. They may also form in the accretion disks of active galactic nuclei. Such compact binaries can emit gravitational waves (GWs) in the low-frequency band (0.001-1 Hz) that are detectable by several planned space-borne GW observatories. We show that the orbital axis of the compact binary may experience significant variation due to the frame-dragging effect associated with the spin of the SMBH. The dynamical behavior of the orbital axis can be understood analytically as a resonance phenomenon. We show that rate of change of the binary orbital axis encodes the information on the spin of the SMBH. Therefore detecting GWs from compact binaries around SMBHs, particularly the modulation of the waveform associated with the variation of the binary orbital axis, can provide a new probe on the spins of SMBHs.

preprint2022arXiv

Secular chaos in white-dwarf planetary systems: origins of metal pollution and short-period planetary companions

Secular oscillations in multi-planet systems can drive chaotic evolution of a small inner body through non-linear resonant perturbations. This "secular chaos" readily pushes the inner body to an extreme eccentricity, triggering tidal interactions or collision with the central star. We present a numerical study of secular chaos in systems with two planets and test particles using the ring-averaging method, with emphasis on the relationship between the planets' properties and the time-scale and efficiency of chaotic diffusion. We find that secular chaos can excite extreme eccentricities on time-scales spanning several orders of magnitude in a given system. We apply our results to the evolution of planetary systems around white dwarfs (WDs), specifically the tidal disruption and high-eccentricity migration of planetesimals and planets. We find that secular chaos in a planetesimal belt driven by large ($\gtrsim 10 M_{\oplus}$), distant ($\gtrsim 10 \, {\rm au}$) planets can sustain metal accretion onto a WD over Gyr time-scales. We constrain the total mass of planetesimals initially present within the chaotic zone by requiring that the predicted mass delivery rate to the Roche limit be consistent with the observed metal accretion rates of WDs with atmospheric pollution throughout the cooling sequence. Based on the occurrence of long-period exoplanets and exo-asteroid belts, we conclude that secular chaos can be a significant (perhaps dominant) channel for polluting solitary WDs. Secular chaos can also produce short-period planets and planetesimals around WDs in concert with various circularization mechanisms. We discuss prospects for detecting exoplanets driving secular chaos around WDs using direct imaging and microlensing.

preprint2021arXiv

Dynamical Tides in Eccentric Binaries Containing Massive Main-Sequence Stars: Analytical Expressions

Tidal evolution of eccentric binary systems containing at least one massive main-sequence (MS) star plays an important role in the formation scenarios of merging compact-object binaries. The dominant dissipation mechanism in such systems involves tidal excitation of outgoing internal gravity waves at the convective-radiative boundary and dissipation of the waves at the stellar envelope/surface. We have derived analytical expressions for the tidal torque and tidal energy transfer rate in such binaries for arbitrary orbital eccentricities and stellar rotation rates. These expressions can be used to study the spin and orbital evolution of eccentric binaries containing massive MS stars, such as the progenitors of merging neutron star binaries. Applying our results to the PSR J0045-7319 system, which has a massive B-star companion and an observed, rapidly decaying orbit, we find that for the standard radius of convective core based on non-rotating stellar models, the B-star must have a significant retrograde and differential rotation in order to explain the observed orbital decay rate. Alternatively, we suggest that the convective core may be larger as a result of rapid stellar rotation and/or mass transfer to the B-star in the recent past during the post-MS evolution of the pulsar progenitor.

preprint2021arXiv

Enhanced Lidov-Kozai migration and the formation of the transiting giant planet WD1856+534b

We investigate the possible origin of the transiting giant planet WD1856+534b, the first strong exoplanet candidate orbiting a white dwarf, through high-eccentricity migration (HEM) driven by the Lidov-Kozai (LK) effect. The host system's overall architecture is an hierarchical quadruple in the '2+2' configuration, owing to the presence of a tertiary companion system of two M-dwarfs. We show that a secular inclination resonance in 2+2 systems can significantly broaden the LK window for extreme eccentricity excitation ($e \gtrsim 0.999$), allowing the giant planet to migrate for a wide range of initial orbital inclinations. Octupole effects can also contribute to the broadening of this 'extreme' LK window. By requiring that perturbations from the companion stars be able to overcome short-range forces and excite the planet's eccentricity to $e \simeq 1$, we obtain an absolute limit of $a_{1} \gtrsim 8 \, {\rm AU} \, (a_{3} / 1500 \, {\rm AU})^{6/7}$ for the planet's semi-major axis just before migration (where $a_{3}$ is the semi-major axis of the 'outer' orbit). We suggest that, to achieve a wide LK window through the 2+2 resonance, WD1856b likely migrated from $30 \, {\rm AU} \lesssim a_{1} \lesssim 60 \, {\rm AU}$, corresponding to $\sim 10$--$20 \, {\rm AU}$ during the host's main-sequence phase. We discuss possible difficulties of all flavours of HEM affecting the occurrence rate of short-period giant planets around white dwarfs.

preprint2021arXiv

Hierarchical Black-Hole Mergers in Multiple Systems: Constrain the Formation of GW190412, GW190814 and GW190521-like events

The merging black-hole (BH) binaries GW190412, GW190814 and GW190521 from the third LIGO/VIRGO observing run exhibit some extraordinary properties, including highly asymmetric masses, significant spin, and component mass in the "mass gap". These features can be explained if one or both components of the binary are the remnants of previous mergers. In this paper, we explore hierarchical mergers in multiple stellar systems, taking into account the natal kick and mass loss due to the supernova explosion (SN) on each component, as well as the merger kick received by the merger remnant. The binaries that have survived the SNe and kicks generally have too wide orbital separations to merge by themselves, but can merge with the aid of an external companion that gives rise to Lidov-Kozai oscillations. The BH binaries that consist of second-generation BHs can also be assembled in dense star clusters through binary interactions. We characterize the parameter space of these BH binaries by merger fractions in an analytical approach. Combining the distributions of the survived binaries, we further constrain the parameters of the external companion, using the analytically formulated tertiary perturbation strength. We find that to produce the three LIGO/VIRGO O3 events, the external companions must be at least a few hundreds $M_\odot$, and fall in the intermediate-mass BH and supermassive BH range. We suggest that GW190412, GW190814 and GW190521 could all be produced via hierarchical mergers in multiples, likely in a nuclear star cluster, with the final merger induced by a massive BH.

preprint2021arXiv

Inclination Dynamics of Resonant Planets under the Influence of an Inclined External Companion

Recent observations suggest that a large fraction of Kepler super-Earth systems have external giant planet companions (cold Jupiters), which can shape the architecture of the inner planets, in particular their mutual inclinations. The dynamical perturbation from cold Jupiters may account for the population of misaligned planets in the Kepler data. The effectiveness of this mechanism can be hindered by a strong planet--planet coupling in the inner system. In this paper, we study how mean-motion resonances (MMRs) affect this coupling and the expected misalignment. We derive approximate analytical expressions for the mutual inclination excitations in the inner planet system induced by an inclined companion, for various period ratios and perturber properties. In most cases, the mutual inclination is proportional to a dimensionless parameter that characterizes the strength of the perturber relative to the coupling in the inner system. We show that the MMR strengthens the inner coupling, reducing the mutual inclination induced by the perturber by a factor of a few. We find that the resonance is resilient to the perturbation, and derive a criterion for the libration of the resonant angle. Our results have applications for constraining unseen planetary perturbers, and for understanding the architecture of multiplanet systems.

preprint2020arXiv

Excitation of Spin-Orbit Misalignments in Stellar Binaries with Circumbinary Disks: Application to DI Herculis

The large spin-orbit misalignments in the DI Herculis stellar binary system have resolved the decades-long puzzle of the anomalously slow apsidal precession rate, but raise new questions regarding the origin of the obliquities. This paper investigates obliquity evolution in stellar binaries hosting modestly-inclined circumbinary disks. As the disk and binary axes undergo mutual precession, each oblate star experiences a torque from its companion star, so that the spin and orbital axes undergo mutual precession. As the disk loses mass through a combination of winds and accretion, the system may be captured into a high-obliquity Cassini state (a spin-orbit resonance). The final obliquity depends on the details of the disk dispersal. We construct a simple disk model to emulate disk dispersal due to viscous accretion and photoevaporation, and identify the necessary disk properties for producing the observed obliquities in DI Herculis. The disk must be massive (at least $10 \%$ of the binary mass). If accretion onto the binary is suppressed, the observed high stellar obliquities are reproduced with a binary-disk inclination of $\sim 5^\circ - 10^\circ$, but if substantial accretion occurs, the inclination must be larger, $\sim 20^\circ - 30^\circ$. If moderate accretion occurs, initially the disk must lose its mass slowly, but eventually lose its remaining mass abruptly, analogous to the observed two-timescale behavior for disks around T-Tauri stars. The spin feedback on the binary orbit causes the binary-disk inclination to decay as the obliquity evolves, a feature that is absent from the standard Cassini state treatment.

preprint2020arXiv

Fast radio burst counterparts and their implications for the central engine

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.

preprint2020arXiv

Merging Compact Binaries Near a Rotating Supermassive Black Hole: Eccentricity Excitation due to Apsidal Precession Resonance

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.

preprint2020arXiv

Periodic Fast Radio Bursts with Neutron Star Free/Radiative Precession

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.

preprint2020arXiv

Planetary Spin and Obliquity from Mergers

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 $Θ_{\rm SL}$ of the merger products have simple analytical forms: $f_{S} \propto S$ and $f_{\cosΘ_{\rm SL}} \propto (1-\cos^2Θ_{\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.

preprint2020arXiv

Tidal Evolution of Eccentric Binaries Driven by Convective Turbulent Viscosity

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 $ω_{\rm tide}$; when $ω_{\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 giant branch (GB) 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 GB stars on eccentric, small pericentre orbits, but is suppressed by a few orders of magnitude in solar-type stars due to viscosity reduction.

preprint2019arXiv

Circumbinary Accretion from Finite and Infinite Disks

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Ω_{\rm{b}}$ for $q_{\rm{b}}>0.5$ and $Ω_{\rm{b}}$ for $q_{\rm{b}}<0.5$, ($Ω_{\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Ω_{\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.

preprint2019arXiv

In-Situ Scattering of Warm Jupiters and Implications for Dynamical Histories

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.

preprint2019arXiv

Pulsed Disc Accretion Driven by Hot Jupiters

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.

preprint2016arXiv

Axion-photon Propagation in Magnetized Universe

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

preprint2016arXiv

Disruption of Planetary Orbits Through Evection Resonance with an External Companion: Circumbinary Planets and Multiplanet Systems

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.

preprint2016arXiv

Extended Transiting Disks and Rings Around Planets and Brown Dwarfs: Theoretical Constraints

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.

preprint2016arXiv

Formation and Stellar Spin-Orbit Misalignment of Hot Jupiters from Lidov-Kozai Oscillations in Stellar Binaries

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.

preprint2016arXiv

Pulsed Accretion onto Eccentric and Circular Binaries

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.

preprint2016arXiv

Solar Obliquity Induced by Planet Nine: Simple Calculation

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.

preprint2016arXiv

The formation efficiency of close-in planets via Lidov-Kozai migration: analytic calculations

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.

preprint2016arXiv

Tidal Dissipation and Evolution of White Dwarfs Around Massive Black Holes: An Eccentric Path to Tidal Disruption

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.

preprint2016arXiv

Viscous Hydrodynamics Simulations of Circumbinary Accretion Discs: Variability, Quasi-Steady State, and Angular Momentum Transfer

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.

preprint2015arXiv

Analytical Model of Tidal Distortion and Dissipation for a Giant Planet with a Viscoelastic Core

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.

preprint2015arXiv

Electromagnetic Torques, Precession and Evolution of Magnetic Inclination of Pulsars

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.

preprint2015arXiv

Inertial-Acoustic Oscillations of Black-Hole Accretion Discs with Large-Scale Poloidal Magnetic Fields

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.

preprint2015arXiv

Rossby Wave Instability and Long-Term Evolution of Dead Zones in Protoplanetary Discs

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

preprint2015arXiv

Suppression of extreme orbital evolution in triple systems with short range forces

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.

preprint2015arXiv

Survival of Planets Around Shrinking Stellar Binaries

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.

preprint2015arXiv

Tidal Truncation of Inclined Circumstellar and Circumbinary Discs in Young Stellar Binaries

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.

preprint2015arXiv

Viscous Driving of Global Oscillations in Accretion Discs Around Black Holes

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 $α$-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 $κ_\mathrm{max}$, non-axisymmetric with frequencies close to multiples of of the innermost stable orbit frequency $Ω_\mathrm{ISCO}$, or hybrid modes whose frequencies are linear combinations of these two frequencies. Small values of the viscosity parameter $α$ primarily produce non-axisymmetric modes, while axisymmetric modes become dominant for large $α$. 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.

preprint2014arXiv

Chaotic Dynamics of Stellar Spin Driven by Planets Undergoing Lidov-Kozai Oscillations: Resonances and Origin of Chaos

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 $ε$, 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 $ε$) 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.

preprint2014arXiv

Chaotic Dynamics of Stellar Spin in Binaries and the Production of Misaligned Hot Jupiters

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.

preprint2014arXiv

Dynamical Tides in Compact White Dwarf Binaries: Influence of Rotation

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.

preprint2014arXiv

Evolution of linear warps in accretion discs and applications to protoplanetary discs in binaries

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 $α$, 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.

preprint2014arXiv

Physics in Very Strong Magnetic Fields: Introduction and Overview

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.

preprint2014arXiv

Star-Disc-Binary Interactions in Protoplanetary Disc Systems and Primordial Spin-Orbit Misalignments

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.

preprint2014arXiv

The Light Curve and Internal Magnetic Field of the Mode-Switching Pulsar PSR B0943+10

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.

preprint2014arXiv

Theory of Disk Accretion onto Magnetic Stars

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.

preprint2013arXiv

Assembly of Protoplanetary Disks and Inclinations of Circumbinary Planets

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.

preprint2013arXiv

Corotation resonance and overstable oscillations in black-hole accretion discs: general-relativistic calculations

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.

preprint2013arXiv

Non-radial Oscillations in Rotating Giant Planets with Solid Cores: Application to Saturn and its Rings

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.

preprint2013arXiv

Viscoelastic Tidal Dissipation in Giant Planets and Formation of Hot Jupiters Through High-Eccentricity Migration

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(ω)]$, as a function of the forcing frequency $ω$. 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(ω)]$, 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.

preprint2013arXiv

White Paper for Blazar Observations with a GEMS-like X-ray Polarimetry Mission

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.

preprint2013arXiv

White Paper on GEMS Study of Polarized X-rays from Neutron Stars

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.

preprint2013arXiv

X-ray Polarization from Black Holes: GEMS Scientific White Paper

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.

preprint2013arXiv

XIPE: the X-ray Imaging Polarimetry Explorer

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).

preprint2012arXiv

DC Circuit Powered by Orbital Motion: Magnetic Interactions in Compact Object Binaries and Exoplanetary Systems

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.

preprint2012arXiv

Dynamical Tides in Compact White Dwarf Binaries: Helium Core White Dwarfs, Tidal Heating, and Observational Signatures

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 ω. Like carbon-oxygen (CO) WDs studied in our previous paper, we find that the dimensionless tidal torque F(ω) (inversely proportional to the effective tidal quality factor) has an erratic dependence on ω. On average, F(ω) scales approximately as ω^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.

preprint2012arXiv

Dynamical Tides in Compact White Dwarf Binaries: Tidal Synchronization and Dissipation

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 ω=2(Ω-Ω_s) (where Ωis the orbital frequency and Ω_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 ωin an erratic way. On average, \dot E_{\rm tide} scales approximately as Ω^5ω^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 Ω_c, corresponding to an orbital period of order one hour (depending on WD models), dynamical tides efficiently drive Ω_s toward Ω, although a small, almost constant degree of asynchronization (Ω-Ω_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.

preprint2012arXiv

Dynamics of the Innermost Accretion Flows Around Compact Objects: Magnetosphere-Disc Interface, Global Oscillations and Instabilities

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.

preprint2012arXiv

High-Frequency QPOs and Overstable Oscillations of Black-Hole Accretion Disks

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.

preprint2012arXiv

Merging neutron star binaries: equation of state and electrodynamics

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.

preprint2012arXiv

Simulations of Overstable Inertial-acoustic Modes in Black-Hole Accretion Discs

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.

preprint2012arXiv

The Hidden Magnetic Field of The Young Neutron Star in Kesteven 79

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.

preprint2012arXiv

Tidal Dissipation in Planet-Hosting Stars: Damping of Spin-Orbit Misalignment and Survival of Hot Jupiters

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.

preprint2012arXiv

Tidal Excitation of Oscillation Modes in Compact White Dwarf Binaries: I. Linear Theory

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.

preprint2012arXiv

Tidal Novae in Compact Binary White Dwarfs

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.

preprint2011arXiv

Distorted, non-spherical transiting planets: impact on the transit depth and on the radius determination

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.

preprint2011arXiv

Dynamical Tides in Eccentric Binaries and Tidally-Excited Stellar Pulsations in KEPLER KOI-54

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.

preprint2011arXiv

Faraday Rotation in Pulsar Magnetosphere

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 $λ^{0.5}$ instead of the usual $λ^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.

preprint2010arXiv

Corotational Instability of Inertial-Acoustic Modes in Black-Hole Accretion Discs: Non-Barotropic Flows

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, $ζ=κ^2/(2ΩΣ)$ (where $Ω, κ, Σ$ 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äsälä 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 $α$-viscosity parameterizations. We find that for disc fluids with constant adiabatic index $Γ$, super-reflection and mode growth depend on the gradient of the effective vortensity, $ζ_{\rm eff} = ζ/S^{2/Γ}$ (where $S \equiv P/Σ^Γ$ measures the entropy): when $dζ_{\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.

preprint2010arXiv

Corotational Instability, Magnetic Resonances and Global Inertial-Acoustic Oscillations in Magnetized Black-Hole Accretion Discs

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 $ω/m$ equals the disc rotation rate $Ω$), when the fluid vortensity, $ζ=κ^2/(2ΩΣ)$ (where $κ$ and $Σ$ 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=ω-mΩ$, 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.

preprint2010arXiv

Evolution of Spin Direction of Accreting Magnetic Protostars and Spin-Orbit Misalignment in Exoplanetary Systems

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.

preprint2010arXiv

Evolution of Spin Direction of Accreting Magnetic Protostars and Spin-Orbit Misalignment in Exoplanetary Systems: II. Warped Discs

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.

preprint2010arXiv

Low-$T/|W|$ instabilities in differentially rotating proto-neutron stars with magnetic fields

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.

preprint2010arXiv

Mass Transfer, Transiting Stream and Magnetopause in Close-in Exoplanetary Systems with Applications to WASP-12

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.

preprint2010arXiv

Pair cascades in the magnetospheres of strongly-magnetized neutron stars

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.

preprint2010arXiv

Papaloizou-Pringle Instability of Magnetized Accretion Tori

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ϕ}\go 0.2rΩ$ (where $Ω$ 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ϕ}\go 0.2 rΩ$). 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.

preprint2009arXiv

Polarization Changes of Pulsars due to Wave Propagation Through Magnetospheres

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.

preprint2009arXiv

Polarization Evolution in A Strongly Magnetized Vacuum: QED Effect and Polarized X-ray Emission from Magnetized Neutron Stars

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.

preprint2008arXiv

Condensed Surfaces of Magnetic Neutron Stars, Thermal Surface Emission, and Particle Acceleration Above Pulsar Polar Caps

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.

preprint2008arXiv

Corotational Damping of Diskoseismic C-modes in Black Hole Accretion Discs

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.

preprint2008arXiv

Effects of Magnetic Fields on the Diskoseismic Modes of Accreting Black Holes

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.

preprint2008arXiv

Interface Modes and Their Instabilities in Accretion Disc Boundary Layers

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.

preprint2008arXiv

Radiative transitions of the helium atom in highly magnetized neutron star atmospheres

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.

preprint1999arXiv

Tidal Interaction between a Fluid Star and a Kerr Black Hole in Circular Orbit

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.

preprint1996arXiv

Innermost Stable Circular Orbit of Inspiraling Neutron-Star Binaries: Tidal Effects, Post-Newtonian Effects and the Neutron-Star Equation of State

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.

preprint1994arXiv

Resonant Oscillations and Tidal Heating in Coalescing Binary Neutron Stars

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äisälä 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.