BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Cristobal Petrovich

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

ResearcherUnclaimed source record
astro-ph.EPastro-ph.HEastro-ph.SRastro-ph.GAcond-mat.supr-congr-qc

Catalog footprint

What is connected

20works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

20 published item(s)

preprint2022arXiv

Eccentric Mergers of Intermediate-Mass Black Holes from Evection Resonances in AGN Disks

We apply the theory of nonlinear resonance capture to the problem of a black hole binary (BHB) orbiting a supermassive black hole (SMBH) while embedded in the accretion disk of an active galactic nucleus (AGN). If successful, resonance capture can trigger dramatic growth in the BHB eccentricity, with important consequences for the BHB merger timescale, as well as for the gravitational wave (GW) signature of such an eccentric merger. This resonance capture may occur when the orbital period around the SMBH (the "outer binary") and the apsidal precession of the BHB (the "inner binary") are in a 1:1 commensurability. This effect is analogous to the phenomenon of lunar evection resonance in the early Sun-Earth-Moon system, with the distinction that in the present case, the BHB apsidal precession is due to general relativity, rather than rotationally-induced distortion. In contrast to the case of lunar evection, however, the inner binary undergoes orbital decay driven by GW emission, rather than orbital expansion driven by tidal dissipation. This distinction fundamentally alters the three-body dynamics, forbidding resonance capture, and limiting eccentricity growth. However, if the BHB migrates through of a gaseous AGN disk, the change in the outer binary can counterbalance the suppressing effect of BHB decay, permitting evection resonance capture and the production of eccentric BHB mergers. We compute the likelihood of resonance capture assuming an agnostic distribution of parameters for the three bodies involved and for the properties of the AGN disk. We find that intermediate-mass ratio BHBs (involving an intermediate-mass black hole and a stellar-mass black hole) are the most likely to be captured into evection resonance and thus undergo an eccentric merger. We also compute the GW signature of these mergers, showing that they can enter the LISA band while eccentric.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mutual Inclination of Ultra-Short-Period Planets with Time Varying Stellar J2-moment

Systems with ultra-short-period planets (USPs) tend to possess larger mutual inclinations compared to those with planets located farther from their host stars. This could be explained due to precession caused by stellar oblateness at early times when the host star was rapidly spinning. However, stellar oblateness reduces over time due to the decrease in the stellar rotation rate, and this may further shape the planetary mutual inclinations. In this work, we investigate in detail how the final mutual inclination varies under the effect of a decreasing $J_2$. We find that different initial parameters (e.g., the magnitude of $J_2$ and planetary inclinations) will contribute to different final mutual inclinations, providing a constraint on the formation mechanisms of USPs. In general, if the inner planets start in the same plane as the stellar equator (or co-planar while misaligned with the stellar spin-axis), the mutual inclination decreases (or increases then decreases) over time due to the decay of the $J_2$ moment. This is because the inner orbit typically possesses less orbital angular momentum than the outer ones. However, if the outer planet is initially aligned with the stellar spin while the inner one is misaligned, the mutual inclination nearly stays the same. Overall, our results suggest that either the USP planets formed early and acquired significant inclinations (e.g., $\gtrsim30^\circ$ with its companion or $\gtrsim10^\circ$ with its host star spin-axis for Kepler-653c) or they formed late ($\gtrsim$Gyr) when their host stars rotate slower.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Substructures in protoplanetary disks imprinted by compact planetary systems

The substructures observed in protoplanetary disks may be the signposts of embedded planets carving gaps or creating vortices. The inferred masses of these planets often fall in the Jovian regime despite their low abundance compared to lower-mass planets, partly because previous works often assume that a single substructure (a gap or vortex) is caused by a single planet. In this work, we study the possible imprints of compact systems composed of Neptune-like planets ($\sim10-30\;M_\oplus$) and show that long-standing vortices are a prevalent outcome when their inter-planetary separation ($Δa$) falls below $\sim8$ times $H_{\rm p}$ -- the average disk's scale height at the planets locations. In simulations where a single planet is unable to produce long-lived vortices, two-planet systems can preserve them for at least $5,000$ orbits in two regimes: i) fully-shared density gaps with elongated vortices around the stable Lagrange points $L_4$ and $L_5$ for the most compact planet pairs ($Δa \lesssim 4.6\; H_{\rm p}$); ii) partially-shared gaps for more widely spaced planets ($Δa \sim 4.6 - 8\;H_{\rm p}$) forming vortices in a density ring between the planets through the Rossby wave instability. The latter case can produce vortices with a wide range of aspect ratios down to $\sim3$ and can occur for planets captured into the 3:2 (2:1) mean-motion resonances for disk's aspects ratios of $h\gtrsim 0.033$ ($h\gtrsim 0.057$). We suggest that their long lifetimes are sustained by the interaction of spiral density waves launched by the neighboring planets. Overall, our results show that distinguishing imprint of compact systems with Neptune-mass planets are long-lived vortices inside the density gaps, which in turn are shallower than single-planet gaps for a fixed gap width.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Warm Neptune GJ 3470b has a Polar Orbit

The warm Neptune GJ 3470b transits a nearby ($d=29$pc) bright slowly rotating M1.5-dwarf star. Using spectroscopic observations during two transits with the newly commissioned NEID spectrometer on the WIYN 3.5m Telescope at Kitt Peak Observatory, we model the classical Rossiter-Mclaughlin effect yielding a sky-projected obliquity of $λ=98_{-12}^{+15\:\circ}$ and a $v \sin i = 0.85_{-0.33}^{+0.27}$km/s. Leveraging information about the rotation period and size of the host star, our analysis yields a true obliquity of $ψ=95_{-8}^{+9\:\circ}$, revealing that GJ 3470b is on a polar orbit. Using radial velocities from HIRES, HARPS and the Habitable-zone Planet Finder, we show that the data are compatible with a long-term RV slope of $\dotγ = -0.0022 \pm 0.0011$m/s/day over a baseline of 12.9 years. If the RV slope is due to acceleration from another companion in the system, we show that such a companion is capable of explaining the polar and mildly eccentric orbit of GJ 3470b using two different secular excitation models. The existence of an outer companion can be further constrained with additional RV observations, Gaia astrometry, and future high-contrast imaging observations. Lastly, we show that tidal heating from GJ 3470b's mild eccentricity has most likely inflated the radius of GJ 3470b by a factor of $\sim$1.5-1.7, which could help account for its evaporating atmosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Formation of counter-rotating and highly eccentric massive black hole binaries in galaxy mergers

Supermassive black hole (SMBH) binaries represent the main target for missions such as the Laser Interferometer Space Antenna and Pulsar Timing Arrays. The understanding of their dynamical evolution prior to coalescence is therefore crucial to improving detection strategies and for the astrophysical interpretation of the gravitational wave data. In this paper, we use high-resolution $N$-body simulations to model the merger of two equal-mass galaxies hosting a central SMBH. In our models, all binaries are initially prograde with respect to the galaxy sense of rotation. But, binaries that form with a high eccentricity, $e\gtrsim 0.7$, quickly reverse their sense of rotation and become almost perfectly retrograde at the moment of binary formation. The evolution of these binaries proceeds towards larger eccentricities, as expected for a binary hardening in a counter-rotating stellar distribution. Binaries that form with lower eccentricities remain prograde and at comparatively low eccentricities. We study the origin of the orbital flip by using an analytical model that describes the early stages of binary evolution. This model indicates that the orbital plane flip is due to the torque from the triaxial background mass distribution that naturally arises from the galactic merger process. Our results imply the existence of a population of SMBH binaries with a high eccentricity and could have significant implications for the detection of the gravitational wave signal emitted by these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nekhoroshev Estimates for the Survival Time of Tightly Packed Planetary Systems

$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}{Δa} \exp \left(c_2 \frac{Δa}{a} /μ^{1/4}\right)$, where $P$ is the average Keplerian period, $a$ is the average semi major axis, $Δa\ll a$ is the difference between the semi major axes of neighbouring planets, $μ$ 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$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Compact-Object Mergers in the Galactic Center: Evolution in Triaxial Clusters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Warm Jupiters from secular planet-planet interactions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hot Jupiters from Coplanar High-eccentricity Migration

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The stability and fates of hierarchical two-planet systems

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(μ_{\rm in},μ_{\rm out})\right]^{1/3}(a_{\rm out}/a_{\rm in})^{1/2}+1.15$, which should be applicable to planet-star mass ratios $μ_{\rm in},μ_{\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 $μ_{\rm in}>μ_{\rm out}$, while slightly favoring collisions with the star for $μ_{\rm in}<μ_{\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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rotochemical heating of millisecond and classical pulsars with anisotropic and density-dependent superfluid gap models

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ández & 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scattering outcomes of close-in planets: constraints on planet migration

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Steady-state planet migration by the Kozai-Lidov mechanism in stellar binaries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Disk-satellite interaction in disks with density gaps

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planets Near Mean-Motion Resonances

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Density Waves Excited by Low-Mass Planets in Protoplanetary Disks I: Linear Regime

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Long-period thermal oscillations in superfluid millisecond pulsars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The origin of the negative torque density in disk-satellite interaction

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Rotochemical heating in millisecond pulsars with Cooper pairing

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Rotochemical heating in millisecond pulsars: modified Urca reactions with uniform Cooper pairing gaps

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 Δ_n and protons Δ_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 Δ_{thr}= min(Δ_n+ 3Δ_p, 3Δ_n+Δ_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_γ~ (1-4) 10^{32}Δ_{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<Δ_{thr}<0.45 MeV.

0 citations0 reviewsNo code linkedOpen access