$N$-body simulations of non-resonant tightly-packed planetary systems have found that their survival time (i.e. time to first close encounter) grows exponentially with their interplanetary spacing and planetary masses. Although this result has important consequences for the assembly of planetary systems by giants collisions and their long-term evolution, this underlying exponential dependence is not understood from first principles, and previous attempts based on orbital diffusion have only yielded power-law scalings. We propose a different picture, where large deviations of the system from its initial conditions is due to few slowly developing high-order resonances. Thus, we show that the survival time of the system $T$ can be estimated using a heuristic motivated by Nekhoroshev's theorem, and obtain a formula for systems away from overlapping two-body mean-motion resonances as: $T/P=c_1 \frac{a}{\Delta a} \exp \left(c_2 \frac{\Delta a}{a} /\mu^{1/4}\right)$, where $P$ is the average Keplerian period, $a$ is the average semi major axis, $\Delta a\ll a$ is the difference between the semi major axes of neighbouring planets, $\mu$ is the planet to star mass ratio, and $c_1$ and $c_2$ are dimensionless constants. We show that this formula is in good agreement with numerical N-body experiments for $c_1=5 \cdot 10^{-4}$ and $c_2=8$.

Exoplanetary systems host giant planets on substantially non-circular, close-in orbits. We propose that these eccentricities arise in a phase of giant impacts, analogous to the final stage of Solar System assembly that formed Earth's Moon. In this scenario, the planets scatter each other and collide, with corresponding mass growth as they merge. We numerically integrate an ensemble of systems with varying total planet mass, allowing for collisional growth, to show that (1) the high-eccentricity giants observed today may have formed preferentially in systems of higher initial total planet mass, and (2) the upper bound on the observed giant planet eccentricity distribution is consistent with planet-planet scattering. We predict that mergers will produce a population of high-mass giant planets between 1 and 8 au from their stars.

There is significant observational evidence that a large fraction of galactic centers, including those in the Milky Way and M31, host a supermassive black hole (SMBH) embedded in a triaxial nuclear star cluster. In this work, we study the secular orbital evolution of binaries in these environments, and characterize the regions and morphological properties of nuclear star clusters that lead to gravitational wave mergers and/or tidal captures. We show that even a modest level of triaxiality in the density distribution of a cluster (an ellipsoid with axis ratios of 0.7 and 0.95) dramatically enhances the merger rates in the central parsecs of the Galaxy by a factor of up to ~10-30 relative to a spherical density distribution. Moreover, we show that the merger fraction of binaries with semi-major axes in the range 10-100 AU remains above 10% for the entire central parsec of the cluster, reaching values close to unity at a distance of ~0.2-0.4 pc from the SMBH. We understand this large merger efficiency in terms of two distinct mechanisms: i) eccentricity oscillations driven by the dominant axisymmetric part of the cluster potential that are enhanced by the slow modulation of a binary's angular momentum from the triaxial contribution, similar to the well-known octupole-level dynamics in three-body systems; ii) chaotic diffusion of eccentricities arising when the nodal precession timescale of a binary's orbit about the SMBH becomes comparable to its characteristic secular timescale. Overall, our results indicate that galactic centers are significantly more collisional than previously thought, with mergers taking place up to the effective radii of their nuclear star clusters.

Due to fortuitous circumstances, the two giant planets around Kepler-419 have well characterized 3-D orbits. They are nearly coplanar to each other; the inner one has a large eccentricity (~0.82); and the apses of the two orbits librate around anti-alignment. Such a state defies available proposals for large eccentricities. We argue that it is instead uniquely produced by a decaying protoplanetary disk. When the disk was massive, its precessional effect on the planets forced the two apses to center around an anti-aligned state. And as the disk is gradually eroded, the pair of planets are adiabatically transported to a new state where most of the eccentricity (or rather, the angular momentum deficit) is transferred to the inner planet, and the two apses are largely anti-aligned. During this transport, any initial mutual inclination may be reduced or enhanced; either may be compatible with the current constraints. So a primordial disk can drive up planet eccentricities both in resonant planet pairs (as has been shown for GJ 876) and in secularly interacting, non-resonant pairs. The mechanism discussed here may be relevant for forming hot Jupiters and for explaining the observed eccentricities of warm and cold giant planets.

We report the discovery of a transiting planet first identified as a candidate in Sector 1 of the Transiting Exoplanet Survey Satellite (TESS), and then confirmed with precision radial velocities. HD1397b has a mass of ${\rm M_P}$ = $0.335_{-0.018}^{+0.018}$ ${\rm M_J}$, a radius of ${\rm R_P}$ = $1.021_{-0.014}^{+0.015}$ ${\rm M_J}$, and orbits its bright host star ($V = 7.8$ mag) with an orbital period of $11.53508 \pm 0.00057 $ d, on a moderately eccentric orbit ($e$ = 0.210 $\pm$ 0.038). With a mass of ${\rm M}_{\star}$ = $1.284_{-0.016}^{+0.020}$ ${\rm M}_{J}$, a radius of ${\rm R}_{\star}$ = 2.314$_{-0.042}^{+0.049}$ ${\rm R_J}$, and an age of 4.7 $\pm$ 0.2 Gyr, the solar metallicity host star has already departed from the main sequence. We find evidence in the radial velocity measurements for a long term acceleration, and a $P \approx 18$ d periodic signal that we attribute to rotational modulation by stellar activity. The HD1397 system is among the brightest systems currently known to host a transiting planet, which will make it possible to perform detailed follow-up observations in order to characterize the properties of giant planets orbiting evolved stars.

Binaries within the sphere of influence of a massive black hole (MBH) in galactic nuclei are susceptible to the Lidov-Kozai (LK) mechanism, which can drive orbits to high eccentricities and trigger strong interactions within the binary such as the emission of gravitational waves (GWs), and mergers of compact objects. These events are potential sources for GW detectors such as Advanced LIGO and VIRGO. The LK mechanism is only effective if the binary is highly inclined with respect to its orbit around the MBH (within a few degrees of 90 deg), implying low rates. However, close to an MBH, torques from the stellar cluster give rise to the process of vector resonant relaxation (VRR). VRR can bring a low-inclination binary into an `active' LK regime in which high eccentricities and strong interactions are triggered in the binary. Here, we study the coupled LK-VRR dynamics, with implications for LIGO and VIRGO GW sources. We carry out Monte Carlo simulations and find that the merger fraction enhancement due to LK-VRR dynamics is up to a factor of ~10 for the lower end of assumed MBH masses (M_MBH = 1e4 MSun), and decreases sharply with increasing M_MBH. We find that, even in our most optimistic scenario, the baseline BH-BH merger rate is small, and the enhancement by LK-VRR coupling is not large enough to increase the rate to well above the LIGO/VIRGO lower limit, 12 Gpc^{-3} yr^{-1}. For the Galactic Center, the LK-VRR-enhanced rate is ~100 times lower than the LIGO/VIRGO limit, and for M_MBH = 1e4 MSun, the rate barely reaches 12 Gpc^{-3} yr^{-1}.

The binary black hole mergers detected by Advanced LIGO/Virgo have shown no evidence of large black hole spins. However, because LIGO/Virgo best measures the effective combination of the two spins along the orbital angular momentum ($\chi_{\rm eff}$), it is difficult to distinguish between binaries with slowly-spinning black holes and binaries with spins lying in the orbital plane. Here, we study the spin dynamics for binaries with a distant black hole companion. For spins initially aligned with the orbital angular momentum of the binary, we find that $\chi_{\rm eff}$ "freezes" near zero as the orbit decays through the emission of gravitational waves. Through a population study, we show that this process predominantly leads to merging black hole binaries with near-zero $\chi_{\rm eff}$. We conclude that if the detected black hole binaries were formed in triples, then this would explain their low $\chi_{\rm eff}$ without the need to invoke near-zero spins or initially large spin-orbit angles.

We constrain the intrinsic architecture of Kepler planetary systems by modeling the observed multiplicities of the transiting planets (tranets) and their transit timing variations (TTVs). We robustly determine that the fraction of Sun-like stars with Kepler-like planets, $\eta_{\rm Kepler}$, is $30\pm3\%$. Here Kepler-like planets are planets that have radii $R_{\rm p} \gtrsim R_\oplus$ and orbital periods $P<400$~days. Our result thus significantly revises previous claims that more than 50\% of Sun-like stars have such planets. Combining with the average number of Kepler planets per star ($\sim0.9$), we obtain that on average each planetary system has $3.0\pm0.3$ planets within 400 days. We also find that the dispersion in orbital inclinations of planets within a given planetary system, $\sigma_{i,k}$, is a steep function of its number of planets, $k$. This can be parameterized as $\sigma_{i,k}\propto k^\alpha$ and we find that $-4<\alpha<-2$ at 2-$\sigma$ level. Such a distribution well describes the observed multiplicities of both tranets and TTVs with no excess of single tranets. Therefore we do not find evidence supporting the so-called "Kepler dichotomy." Together with a previous study on orbital eccentricities, we now have a consistent picture: the fewer planets in a system, the hotter it is dynamically. We discuss briefly possible scenarios that lead to such a trend. Despite our Solar system not belonging to the Kepler club, it is interesting to notice that the Solar system also has three planets within 400 days and that the inclination dispersion is similar to Kepler systems of the same multiplicity.

Over a hundred rocky planets orbiting Sun-like stars in very short orbital periods (<1 day) have been discovered by the Kepler mission. These planets, known as ultra-short-period (USP) planets, are unlikely to have attained their orbits in situ. Instead, they must have migrated in. Here we propose that these planets reach their current orbits by high-eccentricity migration. In a scaled-down version of the dynamics that may have been experienced by their high mass analog, the hot Jupiters, these planets reach high eccentricities via chaotic secular interactions with their companion planets and then undergo orbital circularization due to dissipation of tides raised on the planet. This proposal is motivated by the following observations: planetary systems observed by Kepler often contain several super-Earths with non-negligible eccentricities and inclinations, and possibly extending beyond AU distances; while only a small fraction of USP planets have known transiting companions, and none closely spaced, we argue that most of them should have companions at periods of ~10-50 days. The outer sibling planets, through secular chaos, can remove angular momentum from the inner most planet, originally at periods of ~5-10 days. When the latter reaches an eccentricity higher than 0.8, it is tidally captured by the central star and becomes an USP planet. This scenario naturally explains the observation that most USP planets have significantly more distant transiting companions compared to their counterparts at slightly longer periods (1-3 days), a feature un-accounted for in other proposed scenarios. Our model also predicts that USP planets should have: (i) spin-orbit angles, and inclinations relative to outer planets, in the range of ~10-50 degrees; (ii) several outer planetary companions extending to beyond AU distances, both of which may be tested by TESS and its follow-up observations.

A significant fraction of main sequence stars are part of a triple system. We study the long-term stability and dynamical outcomes of triple stellar systems using a large number of long-term direct N-body integrations with relativistic precession. We find that the previously proposed stability criteria by Eggleton & Kiseleva 1995 and Mardling & Aarseth 2001 predict the stability against ejections reasonably well for a wide range of parameters. Assuming that the triple stellar systems follow orbital and mass distributions from FGK binary stars in the field, we find that in ~1% and ~0.5% of the triple systems lead to a direct head-on collision (impact velocity ~ escape velocity) between main sequence (MS) stars and between a MS star and a stellar-mass compact object, respectively. We conclude that triple interactions are the dominant channel for direct collisions involving a MS star in the field with a rate of one event every ~100 years in the Milky Way. We estimate that the fraction of triple systems that forms short-period binaries is up to ~23% with only up to ~13% being the result of three-body interactions with tidal dissipation, which is consistent with previous work using a secular code.

The Milky Way and a significant fraction of galaxies are observed to host a central Massive Black Hole (MBH) embedded in a non-spherical nuclear star cluster. We study the secular orbital evolution of compact-object binaries in these environments and characterize the excitation of extremely large eccentricities that can lead to mergers by gravitational radiation. We find that the eccentricity excitation occurs most efficiently when the nodal precession timescale of the binary's orbit around the MBH due to the non-spherical cluster becomes comparable (within a factor of ~10) to the timescale on which the binary is torqued by the MBH due to the Lidov-Kozai (LK) mechanism. We show that in this regime the perturbations due to the cluster increase the fraction of systems that reach extreme eccentricities ($1-e\sim10^{-4}-10^{-6}$) by a factor of ~10-100 compared to the idealized case of a spherical cluster, increasing the merger rates of compact objects by a similar factor. We identify two main channels that lead to this extreme eccentricity excitation: (i) chaotic diffusion of the eccentricities due to resonance overlap; (ii) cluster-driven variations of the mutual inclinations between the binary orbit and its center-of-mass orbit around the MBH, which can intensify the LK oscillations. We estimate that our mechanism can produce black hole-black hole and black hole-neutron star binary merger rates of up to $\approx 15$ $\rm{Gpc}^{-3}yr^{-1}$ and $\approx 0.4$ $\rm{Gpc}^{-3}yr^{-1}$, respectively. Thus, we propose the cluster-enhanced Lidov-Kozai mechanism as a new channel for the merger of compact-object binaries, competing with scenarios that invoke isolated binary evolution or dynamical formation in globular clusters.

TRAPPIST-1 is a late M-dwarf orbited by seven Earth-sized planets with orbital period ratios near a chain of mean motion resonances. Due to uncertain system parameters, most orbital configurations drawn from the inferred posterior distribution are unstable on short timescales of $\sim$ 0.5 Myr, even when including the eccentricity damping effect of tides. By contrast, we show that most physically plausible resonant configurations generated through disk migration are stable even without tidal dissipation on timescales of at least 50 Myr ($10^{10}$ orbits), an increase of at least two orders of magnitude. This result, together with the remarkable chain of period ratios in the system, provide strong evidence for convergent migration naturally emplacing the system near an equilibrium configuration forced by the resonant chain. We provide an openly available database of physically plausible initial conditions for TRAPPIST-1 generated through parametrized planet-disk interactions, as well as bit-by-bit reproducible N-body integrations over $10^9-10^{10}$ orbits.

The hundreds of multiple planetary systems discovered by the \textit{Kepler} mission are typically observed to reside in close-in ($\lesssim0.5$ AU), low-eccentricity, and low-inclination orbits. We run N-body experiments to study the effect that unstable outer ($\gtrsim1$ AU) giant planets, whose end orbital configurations resemble those in the Radial Velocity population, have on these close-in multiple super-Earth systems. Our experiments show that the giant planets greatly reduce the multiplicity of the inner super-Earths and the surviving population can have large eccentricities ($e\gtrsim0.3$) and inclinations ($i\gtrsim20^\circ$) at levels that anti-correlate with multiplicity. Consequently, this model predicts the existence of a population of dynamically hot single-transiting planets with typical eccentricities and inclinations %in the ranges of $\sim 0.1-0.5$ and $\sim 10^\circ-40^\circ$. We show that these results can explain the following observations: (i) the recent eccentricity measurements of \textit{Kepler} super-Earths from transit durations; (ii) the tentative observation that single-transiting systems have a wider distribution of stellar obliquity angles compared to the multiple-transiting systems; (iii) the architecture of some eccentric super-Earths discovered by Radial Velocity surveys such as HD\,125612c. Future observations from \textit{TESS} will reveal many more dynamically hot single transiting planets, for which follow up Radial Velocity studies will be able to test our models and see whether they have outer giant planets.

The Atacama Large Millimeter/submilimeter Array (ALMA) recently revealed a set of nearly concentric gaps in the protoplanetary disk surrounding the young star HL Tau. If these are carved by forming gas giants, this provides the first set of orbital initial conditions for planets as they emerge from their birth disks. Using N-body integrations, we have followed the evolution of the system for 5 Gyr to explore the possible outcomes. We find that HL Tau initial conditions scaled down to the size of typically observed exoplanet orbits naturally produce several populations in the observed exoplanet sample. First, for a plausible range of planetary masses, we can match the observed eccentricity distribution of dynamically excited radial velocity giant planets with eccentricities $>$ 0.2. Second, we roughly obtain the observed rate of hot Jupiters around FGK stars. Finally, we obtain a large efficiency of planetary ejections of $\approx 2$ per HL Tau-like system, but the small fraction of stars observed to host giant planets makes it hard to match the rate of free-floating planets inferred from microlensing observations. In view of upcoming GAIA results, we also provide predictions for the expected mutual inclination distribution, which is significantly broader than the absolute inclination distributions typically considered by previous studies.

The requirement that planetary systems be dynamically stable is often used to vet new discoveries or set limits on unconstrained masses or orbital elements. This is typically carried out via computationally expensive N-body simulations. We show that characterizing the complicated and multi-dimensional stability boundary of tightly packed systems is amenable to machine learning methods. We find that training an XGBoost machine learning algorithm on physically motivated features yields an accurate classifier of stability in packed systems. On the stability timescale investigated ($10^7$ orbits), it is 3 orders of magnitude faster than direct N-body simulations. Optimized machine learning classifiers for dynamical stability may thus prove useful across the discipline, e.g., to characterize the exoplanet sample discovered by the upcoming Transiting Exoplanet Survey Satellite (TESS). This proof of concept motivates investing computational resources to train algorithms capable of predicting stability over longer timescales and over broader regions of phase space.

The presence of a planetary system can shield a planetesimal disk from the secular gravitational perturbations due to distant outer massive objects (planets or stellar companions). As the host star evolves off the main sequence to become a white dwarf, these planets can be engulfed, triggering secular instabilities and leading to the tidal disruptions of small rocky bodies. These disrupted bodies can feed the white dwarfs with rocky material and possibly explain the high-metallicity material in their atmospheres. We illustrate how this mechanism can operate when the gravitational perturbations are due to the Kozai-Lidov mechanism from a stellar binary companion. We show that this mechanism can explain the observed levels of accretion if: (1) the planetary engulfment happens fast compared to the secular timescale, which is generally the case for wide binaries ($>100$ AU) and planetary engulfment during the Asymptotic Giant Branch; (2) the planetesimal disk has a total mass of $\sim10^{-4}-10^{-2}M_\oplus$. We show that this new mechanism can provide a steady supply of material throughout the entire life of the white dwarfs for all cooling ages and can account for a large fraction (up to nearly half) of the observed polluted WDs.

Most warm Jupiters (gas-giant planets with $0.1~{\rm AU}\lesssim a \lesssim1$ AU) have pericenter distances that are too large for significant orbital migration by tidal friction. We study the possibility that the warm Jupiters are undergoing secular eccentricity oscillations excited by an outer companion (a planet or star) in an eccentric and/or mutually inclined orbit. In this model the warm Jupiters migrate periodically, in the high-eccentricity phase of the oscillation when the pericenter distance is small, but are typically observed at much lower eccentricities. We show that the steady-state eccentricity distribution of the warm Jupiters migrating by this mechanism is approximately flat, which is consistent with the observed distribution if and only if we restrict the sample to warm Jupiters that have outer companions detected by radial-velocity surveys. The eccentricity distribution of warm Jupiters without companions exhibits a peak at low eccentricities ($e\lesssim 0.2$) that must be explained by a different formation mechanism. Based on a population-synthesis study we find that high-eccentricity migration excited by an outer planetary companion (i) can account for $\sim20\%$ of the warm Jupiters and most of the warm Jupiters with $e\gtrsim 0.4$, a fraction that is consistent with the observed population of warm Jupiters with outer companions; (ii) can produce most of the observed population of hot Jupiters, with a semimajor axis distribution that matches the observations, but fails to account adequately for $\sim 60\%$ of hot Jupiters with projected obliquities $\lesssim 20^\circ$. Thus $\sim 20\%$ of the warm Jupiters and $\sim 60\%$ of the hot Jupiters can be produced by high-eccentricity migration. We also provide predictions for the expected mutual inclinations and spin-orbit angles of the planetary systems with hot and warm Jupiters produced by high-eccentricity migration.

Aims. For more than 1.5 years we monitored spectroscopically the star KELT-6 (BD+312447), known to host the transiting hot Saturn KELT-6b, because a previously observed long-term trend in radial velocity time series suggested the existence of an outer companion. Methods. We collected a total of 93 new spectra with the HARPS-N and TRES spectrographs. A spectroscopic transit of KELT-6b was observed with HARPS-N, and simultaneous photometry was obtained with the IAC-80 telescope. Results. We proved the existence of an outer planet with a mininum mass M$_{\rm p}$sini=3.71$\pm$0.21 M$_{\rm Jup}$ and a moderately eccentric orbit ($e=0.21_{-0.036}^{+0.039}$) of period P$\sim$3.5 years. We improved the orbital solution of KELT-6b and obtained the first measurement of the Rossiter-McLaughlin effect, showing that the planet has a likely circular, prograde, and slightly misaligned orbit, with a projected spin-orbit angle $\lambda$=$-$36$\pm$11 degrees. We improved the KELT-6b transit ephemeris from photometry, and we provided new measurements of the stellar parameters. KELT-6 appears as an interesting case to study the formation and evolution of multi-planet systems.

We study the possibility that hot Jupiters are formed through the secular gravitational interactions between two planets in eccentric orbits with relatively low mutual inclinations ($\lesssim20^\circ$) and friction due to tides raised on the planet by the host star. We term this migration mechanism Coplanar High-eccentricity Migration because, like disk migration, it allows for migration to occur on the same plane in which the planets formed. Coplanar High-eccentricity Migration can operate from the following typical initial configurations: (i) inner planet in a circular orbit and the outer planet with an eccentricity $\gtrsim0.67$ for $m_{\rm in}/m_{\rm out}(a_{\rm in}/a_{\rm out})^{1/2}\lesssim0.3$; (ii) two eccentric ($\gtrsim0.5$) orbits for $m_{\rm in}/m_{\rm out}(a_{\rm in}/a_{\rm out})^{1/2}\lesssim0.16$. A population synthesis study of hierarchical systems of two giant planets using the observed eccentricity distribution of giant planets shows that Coplanar High-eccentricity Migration produces hot Jupiters with low stellar obliquities ($\lesssim30^\circ$), with a semi-major axis distribution that matches the observations, and at a rate that can account for their observed occurrence. A different mechanism is needed to create large obliquity hot Jupiters, either a different migration channel or a mechanism that tilts the star or the proto-planetary disk. Coplanar High-eccentricity Migration predicts that hot Jupiters should have distant ($a\gtrsim5$ AU) and massive (most likely $\sim1-3$ more massive than the hot Jupiter) companions with relatively low mutual inclinations ($\lesssim 20^\circ$) and moderately high eccentricities ($e\sim0.2-0.5$)

We study the dynamical stability and fates of hierarchical (in semi-major axis) two-planet systems with arbitrary eccentricities and mutual inclinations. We run a large number of long-term numerical integrations and use the Support Vector Machine algorithm to search for an empirical boundary that best separates stable systems from systems experiencing either ejections or collisions with the star. We propose the following new criterion for dynamical stability: $a_{\rm out}(1-e_{\rm out})/[a_{\rm in}(1+e_{\rm in})]>2.4\left[\max(\mu_{\rm in},\mu_{\rm out})\right]^{1/3}(a_{\rm out}/a_{\rm in})^{1/2}+1.15$, which should be applicable to planet-star mass ratios $\mu_{\rm in},\mu_{\rm out}=10^{-4}-10^{-2}$, integration times up to $10^8$ orbits of the inner planet, and mutual inclinations $\lesssim40^\circ$. Systems that do not satisfy this condition by a margin of $\gtrsim0.5$ are expected to be unstable, mostly leading to planet ejections if $\mu_{\rm in}>\mu_{\rm out}$, while slightly favoring collisions with the star for $\mu_{\rm in}<\mu_{\rm out}$. We use our numerical integrations to test other stability criteria that have been proposed in the literature and show that our stability criterion performs significantly better for the range of system parameters that we have explored.

We study the steady-state orbital distributions of giant planets migrating through the combination of the Kozai-Lidov (KL) mechanism due to a stellar companion and friction due to tides raised on the planet by the host star. We run a large set of Monte Carlo simulations that describe the secular evolution of a star-planet-star triple system including the effects from general relativistic precession, stellar and planetary spin evolution, and tides. Our simulations show that KL migration produces Hot Jupiters (HJs) with semi-major axes that are generally smaller than in the observations and they can only explain the observations if the following are both true: (i) tidal dissipation at high eccentricities is at least $\sim 150$ times more efficient than the upper limit inferred from the Jupiter-Io interaction; (ii) highly eccentric planets get tidally disrupted at distances $\gtrsim 0.015$ AU. Based on the occurrence rate and semi-major axis distribution of HJs, we find that KL migration in stellar binaries can produce at most $\sim 20\%$ of the observed HJs. Almost no intermediate-period (semi-major axis $\sim0.1-2$ AU) planets are formed by this mechanism - migrating planets spend most of their lifetimes undergoing KL oscillations at large orbital separations ($>2$ AU) or as Hot Jupiters.

When a rotating neutron star loses angular momentum, the progressive reduction of the centrifugal force makes it contract. This perturbs each fluid element, raising the local pressure and originating deviations from beta equilibrium, inducing reactions that release heat (rotochemical heating). This effect has previously been studied by Fern\'andez & Reisenegger (2005) for non-superfluid neutron stars and by Petrovich & Reisenegger (2010) for superfluid millisecond pulsars. Both studies found that pulsars reach a quasi-steady state in which the compression driving the matter out of beta equilibrium is balanced by the reactions trying to restore the equilibrium. We extend previous studies by considering the effect of density-dependence and anisotropy of the superfluid energy gaps, for the case in which the dominant reactions are the modified Urca processes, the protons are non-superconducting, and the neutron superfluidity is parametrized by models proposed in the literature. By comparing our predictions with the surface temperature of the millisecond pulsar PSR J0437-4715 and upper limits for twenty-one classical pulsars, we find the millisecond pulsar can be only explained by the models with the effectively largest energy gaps (type B models), the classical pulsars require with the gap models that vanish for some angle (type C) and two different envelope compositions. Thus, no single model for neutron superfluidity can simultaneously account for the thermal emission of all available observations of non-accreting neutron stars, possibly due to our neglect of proton superconductivity.

Many exoplanets in close-in orbits are observed to have relatively high eccentricities and large stellar obliquities. We explore the possibility that these result from planet-planet scattering by studying the dynamical outcomes from a large number of orbit integrations in systems with two and three gas-giant planets in close-in orbits (0.05 AU < a < 0.15 AU). We find that at these orbital separations, unstable systems starting with low eccentricities and mutual inclinations ($e\lesssim0.1$, $i\lesssim0.1$) generally lead to planet-planet collisions in which the collision product is a planet on a low-eccentricity, low-inclination orbit. This result is inconsistent with the observations. We conclude that eccentricity and inclination excitation from planet-planet scattering must precede migration of planets into short-period orbits. This result constrains theories of planet migration: the semi-major axis must shrink by 1-2 orders of magnitude without damping the eccentricity and inclination.

The multiple-planet systems discovered by the Kepler mission exhibit the following feature: planet pairs near first-order mean-motion resonances prefer orbits just outside the nominal resonance, while avoiding those just inside the resonance. We explore an extremely simple dynamical model for planet formation, in which planets grow in mass at a prescribed rate without orbital migration or dissipation. We develop an analytic version of this model for two-planet systems in two limiting cases: the planet mass grows quickly or slowly relative to the characteristic resonant libration time. In both cases the distribution of systems in period ratio develops a characteristic asymmetric peak-trough structure around the resonance, qualitatively similar to that observed in the Kepler sample. We verify this result with numerical integrations of the restricted three-body problem. We show that for the 3:2 resonance, where the observed peak-trough structure is strongest, our simple model is consistent with the observations for a range of mean planet masses 20-100M_{\oplus}. This mass range is higher than expected, by at least a factor of three, from the few Kepler planets with measured masses, but part of this discrepancy could be due to oversimplifications in the dynamical model or uncertainties in the planetary mass-radius relation.

Gravitational coupling between a gaseous disk and an orbiting perturber leads to angular momentum exchange between them which can result in gap opening by planets in protoplanetary disks and clearing of gas by binary supermassive black holes (SMBHs) embedded in accretion disks. Understanding the co-evolution of the disk and the orbit of the perturber in these circumstances requires knowledge of the spatial distribution of the torque exerted by the latter on a highly nonuniform disk. Here we explore disk-satellite interaction in disks with gaps in linear approximation both in Fourier and in physical space, explicitly incorporating the disk non-uniformity in the fluid equations. Density gradients strongly displace the positions of Lindblad resonances in the disk (which often occur at multiple locations), and the waveforms of modes excited close to the gap edge get modified compared to the uniform disk case. The spatial distribution of the excitation torque density is found to be quite different from the existing prescriptions: most of the torque is exerted in a rather narrow region near the gap edge where Lindblad resonances accumulate, followed by an exponential fall-off with the distance from the perturber. Despite these differences, for a given gap profile the full integrated torque exerted on the disk agrees with the conventional uniform disk theory prediction at the level of ~10%. The nonlinearity of the density wave excited by the perturber is shown to decrease as the wave travels out of the gap, slowing down its nonlinear evolution and damping. Our results suggest that gap opening in protoplanetary disks and gas clearing around SMBH binaries can be more efficient than the existing theories predict. They pave the way for self-consistent calculations of the gap structure and the orbital evolution of the perturber using accurate prescription for the torque density behavior.

Tidal interaction between a gaseous disk and a massive orbiting perturber is known to result in angular momentum exchange between them. Understanding astrophysical manifestations of this coupling such as gap opening by planets in protoplanetary disks or clearing of gas by binary supermassive black holes (SMBHs) embedded in accretion disks requires knowledge of the spatial distribution of the torque exerted on the disk by a perturber. Recent hydrodynamical simulations by Dong et al (2011) have shown evidence for the tidal torque density produced in a uniform disk to change sign at the radial separation of $\approx 3.2$ scale heights from the perturber's orbit, in clear conflict with the previous studies. To clarify this issue we carry out a linear calculation of the disk-satellite interaction putting special emphasis on understanding the behavior of the perturbed fluid variables in physical space. Using analytical as well as numerical methods we confirm the reality of the negative torque density phenomenon and trace its origin to the overlap of Lindblad resonances in the vicinity of the perturber's orbit - an effect not accounted for in previous studies. These results suggest that calculations of the gap and cavity opening in disks by planets and binary SMBHs should rely on more realistic torque density prescriptions than the ones used at present.

Density waves excited by planets embedded in protoplanetary disks play a central role in planetary migration and gap opening processes. We carry out 2D shearing sheet simulations to study the linear regime of wave evolution with the grid-based code Athena, and provide detailed comparisons with the theoretical predictions. Low mass planets (down to ~0.03 Earth mass at 1 AU) and high spatial resolution (256 grid points per scale height) are chosen to mitigate the effects of wave nonlinearity. To complement the existing numerical studies, we focus on the primary physical variables such as the spatial profile of the wave, torque density, and the angular momentum flux carried by the wave, instead of secondary quantities such as the planetary migration rate. Our results show percent level agreement with theory in both physical and Fourier space. New phenomena such as the change of the toque density sign far from the planet are discovered and discussed. Also, we explore the effect of the numerical algorithms, and find that a high order of accuracy, high resolution, and an accurate planetary potential are crucial to achieve good agreement with the theory. We find that the use of a too large time-step without properly resolving the dynamical time scale around the planet produces incorrect results, and may lead to spurious gap opening. Global simulations of planet migration and gap opening violating this requirement may be affected by spurious effects resulting in e.g. the incorrect planetary migration rate and gap opening mass.

In previous papers, we have shown that, as the rotation of a neutron star slows down, it will be internally heated as a consequence of the progressively changing mix of particles (rotochemical heating). In previously studied cases (non-superfluid neutron stars or superfluid stars with only modified Urca reactions), this leads to a quasi-steady state in which the star radiates thermal photons for a long time, possibly accounting for the ultraviolet radiation observed from the millisecond pulsar J0437-4715. For the first time, we explore the phenomenology of rotochemical heating with direct Urca reactions and uniform and isotropic superfluid energy gaps of different sizes. We first do exploratory work by integrating the thermal and chemical evolution equations numerically for different energy gaps, which uncovers a rich phenomenology of stable and unstable solutions. To understand these, we perform a stability analysis around the quasi-steady state, identifying the characteristic times of growing, decaying, and oscillating solutions. For small gaps, the phenomenology is similar to the previously studied cases, in the sense that the solutions quickly converge to a quasi-steady state. For large gaps ($\gtrsim 0.05$ MeV), these solutions become unstable, leading to a limit-cycle behavior of periodicity $\sim 10^{6-7} yr, in which the star is hot ($T_s\gtrsim 10^5$ K) for a small fraction of the cycle ($\sim 5- 20 %$), and cold for a longer time.

Context: When a rotating neutron star loses angular momentum, the reduction in the centrifugal force makes it contract. This perturbs each fluid element, raising the local pressure and originating deviations from beta equilibrium that enhance the neutrino emissivity and produce thermal energy. This mechanism is named rotochemical heating and has previously been studied for neutron stars of nonsuperfluid matter, finding that they reach a quasi-steady configuration in which the rate at which the spin-down modifies the equilibrium concentrations is the same at which neutrino reactions restore the equilibrium. Aims: We describe the thermal effects of Cooper pairing with spatially uniform energy gaps of neutrons \Delta_n and protons \Delta_p on the rotochemical heating in millisecond pulsars (MSPs) when only modified Urca reactions are allowed. By this, we may determine the amplitude of the superfluid energy gaps for the neutron and protons needed to produce different thermal evolution of MSPs. Results: We find that the chemical imbalances in the star grow up to the threshold value \Delta_{thr}= min(\Delta_n+ 3\Delta_p, 3\Delta_n+\Delta_p), which is higher than the quasi-steady state achieved in absence of superfluidity. Therefore, the superfluid MSPs will take longer to reach the quasi-steady state than their nonsuperfluid counterparts, and they will have a higher a luminosity in this state, given by L_\gamma ~ (1-4) 10^{32}\Delta_{thr}/MeV \dot{P}_{-20}/P_{ms}^3 erg s^-1. We can explain the UV emission of the PSR J0437-4715 for 0.05 MeV<\Delta_{thr}<0.45 MeV.

When a rotating neutron star loses angular momentum, the reduction of the centrifugal force makes it contract. This perturbs each fluid element, raising the local pressure and originating deviations from beta equilibrium, inducing reactions that release heat (rotochemical heating). This effect has previously been studied by Fern\'andez and Reisenegger for neutron stars of non-superfluid matter and by Petrovich and Reisenegger for superfluid matter, finding that the system in both cases reaches a quasi-steady state, corresponding to a partial equilibration between compression, due to the loss of angular momentum, and reactions that try to restore the equilibrium. However, Petrovich and Reisenegger assumes a constant value of the superfluid energy gap, whereas theoretical models predict density-dependent gap amplitudes, and therefore gaps that depend on the location in the star. In this work, we try to discriminate between several proposed gap models, comparing predicted surface temperatures to the value measured for the nearest millisecond pulsar, J0437-4715.

When a rotating neutron star loses angular momentum, the reduction in the centrifugal force makes it contract. This perturbs each fluid element, raising the local pressure and originating deviations from beta equilibrium that enhance the neutrino emissivity and produce thermal energy. This mechanism is named rotochemical heating and has previously been studied for neutron stars of non-superfluid matter, finding that they reach a quasi-steady state in which the rate that the spin-down modifies the equilibrium concentrations is the same to that of the neutrino reactions restoring the equilibrium. On the other hand, the neutron star interior is believed to contain superfluid nucleons, which affect the thermal evolution of the star by suppressing the neutrino reactions and the specific heat, and opening new Cooper pairing reactions. In this work we describe the thermal effects of Cooper pairing with spatially uniform energy gaps of neutrons and protons on rotochemical heating in millisecond pulsars (MSPs) when only modified Urca reactions are allowed. We find that the chemical imbalances grow up to a value close to the energy gaps, which is higher than the one of the nonsuperfluid case. Therefore, the surface temperatures predicted with Cooper pairing are higher and explain the recent measurement of MSP J0437-4715.