Orbital Eccentricity Research Articles

Merging hierarchical triple black hole systems with intermediate-mass black holes in population III star clusters

ABSTRACT Theoretical predictions suggest that very massive stars have the potential to form through multiple collisions and eventually evolve into intermediate-mass black holes (IMBHs) within Population III star clusters embedded in mini dark matter haloes. In this study, we investigate the long-term evolution of Population III star clusters, including models with a primordial binary fraction of $f_{\rm b}=0$ and 1, using the N-body simulation code petar. We comprehensively examine the phenomenon of hierarchical triple black holes in the clusters, specifically focusing on their merging inner binary black holes (BBHs), with post-Newtonian correction, by using the tsunami code. Our findings suggest a high likelihood of the inner BBHs containing IMBHs with masses on the order of $\mathcal {O}(100)\,{\rm M}_{\odot }$, and as a result, their merger rate could be up to $0.1{\rm Gpc}^{-3}{\rm yr}^{-3}$. The orbital eccentricities of some merging inner BBHs oscillate over time periodically, known as the Kozai–Lidov oscillation, due to dynamical perturbations. Detectable merging inner BBHs for mHz GW detectors LISA/TianQin/Taiji concentrate within $z\lt 3$. More distant sources would be detectable for CE/ET/LIGO/KAGRA/DECIGO, which are sensitive from $\mathcal {O}(0.1)$Hz to $\mathcal {O}(100)$ Hz. Furthermore, compared with merging isolated BBHs, merging inner BBHs affected by dynamical perturbations from tertiary BHs tend to have higher eccentricities, with a significant fraction of sources with eccentricities closing to 1 at mHz bands. GW observations would help constrain the formation channels of merging BBHs, whether through isolated evolution or dynamical interaction, by examining eccentricities.

Capture of field stars by dark substructures

ABSTRACT We use analytical and N-body methods to study the capture of field stars by gravitating substructures moving across a galactic environment. The majority of stars captured by a substructure move on temporarily bound orbits that are lost to galactic tides after a few orbital revolutions. In numerical experiments where a substructure model is immersed into a sea of field particles on a circular orbit, we find a population of particles that remain bound to the substructure potential for indefinitely long times. This population is absent from substructure models, initially placed outside the galaxy on an eccentric orbit. We show that gravitational capture is most efficient in dwarf spheroidal galaxies (dSphs) on account of their low velocity dispersions and high stellar phase-space densities. In these galaxies, ‘dark’ sub-subhaloes, which do not experience in situ star formation, may capture field stars and become visible as stellar overdensities with unusual properties: (i) they would have a large size for their luminosity, (ii) contain stellar populations indistinguishable from the host galaxy, and (iii) exhibit dark matter (DM)-dominated mass-to-light ratios. We discuss the nature of several ‘anomalous’ stellar systems reported as star clusters in the Fornax and Eridanus II dSphs that exhibit some of these characteristics. DM sub-subhaloes with a mass function ${\rm d}N/{\rm d}M_\bullet \sim M_\bullet ^{-\alpha }$ are expected to generate stellar systems with a luminosity function, ${\rm d}N/{\rm d}M_\star \sim M_\star ^{-\beta }$, where $\beta =(2\alpha +1)/3=1.6$ for $\alpha =1.9$. Detecting and characterizing these objects in dSphs would provide unprecedented constraints on the particle mass and cross-section of a large range of DM particle candidates.

Different Planetary Eccentricity-period (PEP) Distributions of Small and Giant Planets

We used the database of 1040 short-period (1 ≤ P < 200 days) exoplanets radial-velocity orbits to study the planetary eccentricity-period (PEP) distribution. We first divided the sample into low- and high-mass exoplanet subsamples based on the distribution of the (minimum) planetary masses, which displays a clear two-Gaussian distribution, separated at 0.165M J. We then selected 216 orbits, low- and high-mass alike, with eccentricities significantly distinct from circular orbits. The 131 giant-planet eccentric orbits display a clear upper envelope, which we model quantitatively, rises monotonically from zero eccentricity and reaches an eccentricity of 0.8 at P ∼ 100 days. Conversely, the 85 low-mass planetary orbits display a flat eccentricity distribution between 0.1 and 0.5, with almost no dependence on the orbital period. We show that the striking difference between the two PEP distributions is not a result of the detection technique used. The upper envelope of the high-mass planets, also seen in short-period binary stars, is a clear signature of tidal circularization, which probably took place inside the planets, while the small-planet PEP distribution suggests that the circularization was not effective, probably due to dynamical interactions with neighboring planets.

Counter-rotation and slow precession in aligned eccentric nuclear discs due to gravitational wave recoil kicks

ABSTRACT The M31 nucleus contains a supermassive black hole embedded in a massive stellar disc of apsidally aligned eccentric orbits. It has recently been shown that this disc is slowly precessing at a rate consistent with zero. Here, we demonstrate using N-body methods that apsidally aligned eccentric discs can form with a significant ($\sim$0.5) fraction of orbits counter-rotating as the result of a gravitational wave recoil kick of merging supermassive black holes. Higher amplitude kicks map to a larger retrograde fraction in the surrounding stellar population, which in turn results in slow precession. We furthermore show that discs with significant counter-rotation are more stable (i.e. apsidal alignment is most pronounced and long lasting), more eccentric, and have the highest rates of stars entering the black hole’s tidal disruption radius.

Exocomet orbital distribution around beta Pictoris

The $ sim 23$\,Myr young star beta \,Pictoris (beta \,Pic) is a laboratory for planet formation studies because of its observed debris disk, its directly imaged super-Jovian planets beta \,Pic\,b and c, and the evidence of extrasolar comets that regularly transit in front of the star. The most recent evidence of exocometary transits around beta \,Pic came from stellar photometric time series obtained with the TESS space mission. Previous analyses of these transits constrained the orbital distribution of the underlying exocomet population to a range between about 0.03 and 1.3\,AU assuming a fixed transit impact parameter. We examine the distribution of the observed transit durations Delta t$) to infer the orbital surface density distribution (delta ) of the underlying exocomet sample. The effect of the geometric transit probability for circular orbits was properly taken into account, but we assumed that the radius of the transiting comets and their possible clouds of evaporating material are much smaller than the stellar radius. We show that a narrow belt of exocomets around beta \,Pic, in which the transit impact parameters are randomized but the orbital semimajor axes are equal, results in a pile-up of long transit durations. This is in contrast to observations, which reveal a pile-up of short transit durations Delta t 0.1$\,d) and a tail of only a few transits with Delta t &gt; 0.4$\,d. A flat density distribution of exocomets between about 0.03 and 2.5\,AU results in a better match between the resulting Delta t$ distribution and the observations, but the slope of the predicted $ Delta t$ histogram is not sufficiently steep. An even better match to the observations can be produced with a $ beta $ power law. Our modeling reveals a best fit between the observed and predicted $ Delta t$ distribution for $ A more reasonable scenario in which the exocometary trajectories are modeled as hyperbolic orbits can also reproduce the observed $ Delta t$ distribution to some extent. Future studies might reproduce the observed $ Delta t$ distribution with a full exploration of the four-dimensional parameter space of highly eccentric orbits, and they might need to relax our assumption that the transiting objects are smaller than the stellar disk. The number of observed exocometary transits around beta \,Pic is currently too small to validate the previously reported distinction of two distinct exocomet families, but this might be possible with future TESS observations of this star. Our results nevertheless imply that cometary material exists on highly eccentric orbits with a more extended range of semimajor axes than suggested by previous spectroscopic observations.

Primeval very low-mass stars and brown dwarfs – VIII. The first age benchmark L subdwarf, a wide companion to a halo white dwarf

Abstract We report the discovery of five white dwarf + ultracool dwarf systems identified as common proper motion wide binaries in the Gaia Catalogue of Nearby Stars. The discoveries include a white dwarf + L subdwarf binary, VVV 1256−62AB, a gravitationally-bound system located 75.6$^{+1.9}_{-1.8}$ pc away with a projected separation of 1375$^{+35}_{-33}$ au. The primary is a cool DC white dwarf with a hydrogen dominated atmosphere, and has a total age of $10.5^{+3.3}_{-2.1}$ Gyr, based on white dwarf model fitting. The secondary is an L subdwarf with a metallicity of [M/H] = $-0.72^{+0.08}_{-0.10}$ (i.e. [Fe/H] = −0.81 ± 0.10) and Teff = 2298$^{+45}_{-43}$ K based on atmospheric model fitting of its optical to near infrared spectrum, and likely has a mass just above the stellar/substellar boundary. The subsolar metallicity of the L subdwarf and the system’s total space velocity of 406 km/s indicates membership in the Galactic halo, and it has a flat eccentric Galactic orbit passing within 1 kpc of the centre of the Milky Way every ∼0.4 Gyr and extending to 15-31 kpc at apogal. VVV 1256−62B is the first L subdwarf to have a well-constrained age, making it an ideal benchmark of metal-poor ultracool dwarf atmospheres and evolution.

Dust rings trap protoplanets on eccentric orbits and get consumed by them

ABSTRACT We study the orbital evolution and mass growth of protoplanets with masses $M \in [0.1{\!-\!}8]$ M$_{{\oplus }}$ in the vicinity of a dusty ring, using three-dimensional numerical simulations with a two-fluid model and nested-meshes. We find two stable, eccentric orbits that lock the planet in the ring vicinity, thereby inhibiting its migration and allowing it to accrete dust from the ring. One of these orbits has an eccentricity comparable to the aspect ratio of the gaseous disc and has its periastron within the ring, enabling intermittent accretion during each pass. The other orbit has a smaller eccentricity and an apoastron slightly inside the ring. A planet locked at the outer orbit efficiently accretes from the ring and can reach the critical mass for runaway gas accretion on time-scales $\gtrsim 10^5$ yr (for a 10 M$_{{\oplus }}$ dust ring at 10 au), while a planet locked at the inner orbit has a slower growth and might not supersede the super-Earth stage over the disc lifetime. While in our runs a low-mass embryo forming within the ring eventually joins the outer orbit, it is likely that the path taken depends on the specific details of the ring. The trapping on the outer orbit arises from an intermittent, strong thermal force at each passage through the ring, where the accretion rate spikes. It is insensitive to uncertainties that plague models considering planets trapped on circular orbits in rings. It is highly robust and could allow a growing planet to follow an expanding ring over large distances.

Tracking Advanced Planetary Systems (TAPAS) with HARPS-N

Aims. The star HD 118203, classified as a K0 subgiant, was known to harbour a transiting hot Jupiter planet on a 6.1-day eccentric orbit. Previous studies also revealed a linear trend in the radial velocity (RV) domain, indicative of a companion on a wide orbit. Such a hierarchical orbital architecture could be helpful in studies of the origins of hot Jupiters. Methods. We acquired precise RV measurements over 17 yr using the 9.2 m Hobby-Eberly Telescope and the 3.6 m Telescopio Nazionale Galileo. Combining these observations with space-born photometric time series from the Transiting Exoplanet Survey Satellite, we constructed a two-planetary model for the system. Astrometric observations from HIPPARCOS and Gaia were used to constrain the orbital inclination of the wide-orbit companion and its mass. Numerical simulations were used to investigate the dynamics of the system. The photometric data were searched for additional transit-like flux drops. Results. We found that the additional companion is an 11-Jupiter mass planet orbiting HD 118203 on a 14-yr moderately eccentric orbit, constituting a hierarchical planetary system with the hot Jupiter. Both planets were found to be dynamically decoupled mainly due to the general relativistic apsidal precession of the inner planet, marginalising secular interactions. The orbits of both planets might have a relatively low mutual inclination unless the longitudes of the ascending node differ substantially. This configuration favours the coplanar high-eccentricity migration as a path to the present-day orbital configuration. No other transiting planets with radii down to 2 Earth radii and orbital periods less than 100 days were found in the system.

To Be or Not To Be: The Role of Rotation in Modeling Galactic Be X-Ray Binaries

Be X-ray binaries (Be-XRBs) are one of the largest subclasses of high-mass X-ray binaries, comprised of a rapidly rotating Be star and neutron star companion in an eccentric orbit, intermittently accreting material from a decretion disk around the donor. Originating from binary stellar evolution, Be-XRBs are of significant interest to binary population synthesis (BPS) studies, encapsulating the physics of supernovae, common envelope, and mass transfer (MT). Using the state-of-the-art BPS code, POSYDON, which relies on precomputed grids of detailed binary stellar evolution models, we investigate the Galactic Be-XRB population. POSYDON incorporates stellar rotation self-consistently during MT phases, enabling detailed examination of the rotational distribution of Be stars in multiple phases of evolution. Our fiducial BPS and Be-XRB model aligns well with the orbital properties of Galactic Be-XRBs, emphasizing the role of rotational constraints. Our modeling reveals a rapidly rotating population (ω/ω crit ≳ 0.3) of Be-XRB-like systems with a strong peak at intermediate rotation rates (ω/ω crit ≃ 0.6) in close alignment with observations. All Be-XRBs undergo an MT phase before the first compact object forms, with over half experiencing a second MT phase from a stripped helium companion (Case BB). Computing rotationally limited MT efficiencies and applying them to our population, we derive a physically motivated MT efficiency distribution, finding that most Be-XRBs have undergone highly nonconservative MT ( β¯rot≃0.15 ). Our study underscores the importance of detailed angular momentum modeling during MT in interpreting Be-XRB populations, emphasizing this population as a key probe for the stability and efficiency of MT in interacting binaries.

Revisiting the Tertiary-induced Binary Black Hole Mergers: The Role of Superthermal Wide Tertiary Eccentricity Distributions

Recent studies show that the eccentricity distribution of wide binaries (semimajor axis ≳103 au) observed by Gaia tends to favor large eccentricities more strongly than the canonical thermal distribution (P(e) ∝ e)—such distributions are termed “superthermal.” Motivated by this observation, we revisit the formation channel of black hole (BH) binary mergers in triple stellar systems and study the impact of superthermal eccentricity distributions in the outer binaries. We explore the persistence of the highly eccentric outer orbits after each component in a stellar triple has undergone mass loss due to supernova explosions. We find that the outer eccentricity distribution can remain significantly superthermal for modestly hierarchical BH triples satisfying a in/a out ≳ 0.005 (where a in and a out are the semimajor axes of the inner and outer orbits), and are otherwise shaped by mass-loss induced kicks and dynamical instability. We then study the impact of these different outer eccentricity distributions of the remaining BH triples on mergers via the tertiary-induced channel. Of interest, we find that mergers can sometimes be produced even when the initial stellar orbits are near alignment (not subject to the von-Zeipel–Lidov–Kozai effect; ZLK effect) as long as the system is sufficiently hierarchical. On the other hand, although the impact of the octupole-order ZLK effect is much greater when the outer binary is more eccentric, we find that the merger fraction only changes modestly for extreme outer eccentricity distributions because the largest eccentricities tend to lead to dynamical instability.

The Orbit and Dynamical Mass of Polaris: Observations with the CHARA Array

The 30 yr orbit of the Cepheid Polaris has been followed with observations by the Center for High Angular Resolution Astronomy (CHARA) Array from 2016 through 2021. An additional measurement has been made with speckle interferometry at the Apache Point Observatory. Detection of the companion is complicated by its comparative faintness—an extreme flux ratio. Angular diameter measurements appear to show some variation with pulsation phase. Astrometric positions of the companion were measured with a custom grid-based model-fitting procedure and confirmed with the CANDID software. These positions were combined with the extensive radial velocities (RVs) discussed by Torres to fit an orbit. Because of the imbalance of the sizes of the astrometry and RV data sets, several methods of weighting are discussed. The resulting mass of the Cepheid is 5.13 ± 0.28 M ⊙. Because of the comparatively large eccentricity of the orbit (0.63), the mass derived is sensitive to the value found for the eccentricity. The mass combined with the distance shows that the Cepheid is more luminous than predicted for this mass from evolutionary tracks. The identification of surface spots is discussed. This would give credence to the identification of a radial velocity variation with a period of approximately 120 days as a rotation period. Polaris has some unusual properties (rapid period change, a phase jump, variable amplitude, and unusual polarization). However, a pulsation scenario involving pulsation mode, orbital periastron passage, and low pulsation amplitude can explain these characteristics within the framework of pulsation seen in Cepheids.

Obliquity Constraints for the Extremely Eccentric Sub-Saturn Kepler-1656 b

The orbits of close-in exoplanets provide clues to their formation and evolutionary history. Many close-in exoplanets likely formed far out in their protoplanetary disks and migrated to their current orbits, perhaps via high-eccentricity migration (HEM), a process that can also excite obliquities. A handful of known exoplanets are perhaps caught in the act of HEM, as they are observed on highly eccentric orbits with tidal circularization timescales shorter than their ages. One such exoplanet is Kepler-1656 b, which is also the only known nongiant exoplanet (<100 M ⊕) with an extreme eccentricity (e = 0.84). We measured the sky-projected obliquity of Kepler-1656 b by observing the Rossiter–McLaughlin effect during a transit with the Keck Planet Finder. Our data are consistent with an aligned orbit but are also consistent with moderate misalignment with λ < 50° at 95% confidence, with the most likely solution of 35.0−21.6+14.9 deg. A low obliquity would be an unlikely outcome of most eccentricity-exciting scenarios, but we show that the properties of the outer companion in the system are consistent with the coplanar HEM mechanism. Alternatively, if the system is not relatively coplanar (≲20° mutual inclination), Kepler-1656 b may be presently at a rare snapshot of long-lived eccentricity oscillations that do not induce migration. Kepler-1656 b is only the fourth exoplanet with e > 0.8 to have its obliquity constrained; expanding this population will help establish the degree to which orbital misalignment accompanies migration. Future work that constrains the mutual inclinations of outer perturbers will be key for distinguishing plausible mechanisms.

Reducing reentry casualty risk for spacecraft on long-term reentering eccentric inclined geosynchronous (Tundra) disposal orbits

This paper considers disposal of spacecraft and upper stages on eccentric inclined geosynchronous orbits (eccentric IGSOs), specifically the class known as Tundra orbits. The U.S. Government Orbital Debris Mitigation Standard Practices (ODMSP) released in December 2019 include a disposal option to use orbital eccentricity growth for long-term reentry within 200 years that could be used for Tundra orbits. A condition of this disposal option is that reentry casualty risk be limited. An analysis was performed to assess reentry casualty risk for a generic Tundra spacecraft and typical upper stages on Tundra disposal orbits. From propagation sweeps accounting for eccentricity growth of Tundra disposal orbits, several Tundra disposal orbit cases were selected for reentry analysis. A reentry risk analysis for these cases assuming reentry from near-circular orbits was performed. Results show that predicted casualty risk for the generic Tundra spacecraft and typical upper stages well exceed allowable risk limits in the ODMSP. An analysis of reentry from the selected Tundra disposal orbits accounting for the high eccentricity due to eccentricity growth just before reentry was then performed. Results show that the distribution of reentry points on the Earth can be concentrated over ocean in the southern hemisphere where there is less human population. The generic Tundra spacecraft and one of the upper stages considered were then found to be compliant with the limit of 1 in 10,000 expected casualties in the ODMSP. These results indicate promise for wider usage of the long-term re-entry option in the ODMSP.

Astronomically forcing salinity variations in a marginal-marine environment, Bohai Bay Basin, NE China

Exploring the interdependence between watermass conditions and climate change in marginal marine regions is essential for gaining insights into contemporary environmental challenges arising from global warming. Salinity, a key parameter of watermass conditions in coastal areas, plays a crucial role in influencing nutrient availability, stratification, and the spatial-temporal distribution of redox conditions. Here, we focus on the cyclostratigraphic analysis of gamma-ray (GR) data obtained from Eocene strata in the Luo-69 drillcore from the Bohai Bay Basin (NE China). This analysis establishes a high-resolution astronomical time scale (ATS) for the lower third member of the Shahejie Formation (Es3l), spanning 39.04 Ma to 41.46 million years ago (Ma). Various proxies for watermass conditions, such as salinity, sedimentary structures, mineral content, total organic carbon (TOC), organic carbon isotopes and nitrogen isotopes (δ13Corg and δ15N), exhibit a strong correlation with 405-kyr long orbital eccentricity climate forcing. Notably, spectral analysis of salinity data (represented by the Sr/Ba proxy) reveals a clear astronomical control, with 405-kyr eccentricity cycles that align closely with variations in the GR data. Peaks in paleosalinity, as indicated by Sr/Ba (as well as B/Ga and S/TOC), correspond with maxima in 405-kyr long orbital eccentricity cycles in the GR data. As such, we infer that astronomical forcing regulated long-term changes in watermass conditions, likely through mediation of the East Asia monsoon system. Specifically, during eccentricity maxima, the basin would have experienced a warm and humid climate with pronounced seasonality and a robust summer monsoon. During such periods, increased continental runoff would lead to a rise in lake-levels, facilitating better connections with the ocean, and hence increased salinity. By contrast, during eccentricity minima, cooler and more arid conditions with weaker seasonality would result in lower lake levels which together with a global eustatic lowstand would result in as weakened connection with the open ocean restricting the ingress of seawater into the basin lowering salinity. Our dataset spans the Middle Eocene Climatic Optimum (MECO, ∼40.60–40.10 Ma), identified by elevated GR values and a conspicuous negative excursion in oxygen isotopes. This interval is characterized by a significant decline in salinity, likely attributed to enhanced freshwater runoff driven by an extremely warm and humid climate. The paleosalinity fluctuation recorded by the proxies B/Ga, Sr/Ba, and S/TOC in Bohai Bay Basin demonstrate the sensitivity of marginal-marine environment to subtle changes in climate driven by astronomical forcing. Our results provide new insights into how elemental salinity proxies can be utilized in paleoenvironmental reconstruction studies.

Chaotic libration and control of a space tether-sail system in Earth polar orbits with J2 perturbation

The recently proposed space tether-sail system could enable various interesting and important potential space missions thanks to its thin and long connecting tether and solar radiation pressure (SRP) force without consuming propellants. This paper focuses on nonlinear libration and its control of a tether-sail system in Earth polar orbits. Firstly, nonlinear coupled in-plane and out-of-plane dynamics of the system with a rigid mass-distributed tether and two point masses (a chief satellite and sail respectively) is established using the Lagrangian equation method. Secondly, the predicted occurrence of chaotic in- and out-of-plane librational motions is verified utilizing numerical tools such as presenting the time history of libration, the phase plane, the Poincaré section and power spectral density(PSD). Specifically, the paper investigates the sensitivity of nonlinear dynamical responses to initial values, focuses on the chaotic motion caused by the coupling between the in- and out-of-plane librations, the J2 perturbation of the Earth gravitational field and the small eccentricity of the orbits, also studies the nonlinear librational characteristics influenced by the initial mechanical energy (the Hamiltonian). Thirdly, chaotic librational motions will be controlled merely using modulatable SRP force by the sail propellantlessly. A sliding mode controller based on exponential approaching law is developed. To mitigate the effects of control torque saturation, an input compensator based on Radial Basis Function (RBF) neural network is designed. In addition, to address the cumbersome and inefficient manual tuning of control parameters for the sliding mode controller, a parameter tuning algorithm based on Learning Automata is designed. This algorithm is capable of tuning a high-quality set of controller parameters within 100 iterations. The effectiveness of the propellantless actuators for chaotic libration motion control and developed control algorithm is supported by numerical simulations.

Fate of a remnant solid disk around an eccentric giant planet

The composition of giant planets' atmospheres is an important tracer of their formation history. While many theoretical studies investigate the heavy-element accretion within a gaseous protoplanetary disk, the possibility of solid accretion after disk dissipation has not been explored. Here, we focus on the case of a gas giant planet excited to an eccentric orbit and assess the likelihood of solid accretion after disk dissipation. We follow the orbital evolution of the surrounding solid materials and investigate the scattering and accretion of heavy elements in the remnant solid disks. We perform N-body simulations of planetesimals and embryos around an eccentric giant planet. We consider various sizes and orbits for the eccentric planet and determine the fate of planetesimals and embryos. We find that the orbital evolution of solids, such as planetesimals and embryos, is regulated by weak encounters with the eccentric planet rather than strong close encounters. Even in the region where the Safronov number is smaller than unity, most solid materials fall onto the central star or are ejected from the planetary system. We also develop an analytical model of the solid accretion along the orbital evolution of a giant planet, where the accretion probability is obtained as a function of the planetary mass, radius, semi-major axis, eccentricity, inclination, and solid disk thickness. Our model predicts that sim 0.01-0.1 $M_ of solids is accreted onto an eccentric planet orbiting in the outer disk ($ The accreted heavy-element mass increases (decreases) with the eccentricity (inclination) of the planet. We also discuss the possibility of collisions of terrestrial planets and find that $ of the hot Jupiters formed via high-eccentric migration collide with a planet of $10M_ However, we find that solid accretion and collisions with terrestrial planets are minor events for planets in the inner orbit, and a different accretion process is required to enrich eccentric giant planets with heavy elements.

Runaway Eccentricity Growth: A Pathway for Binary Black Hole Mergers in AGN Disks

Binary black holes (BBHs) embedded within the accretion disks that fuel active galactic nuclei (AGN) are promising progenitors for the source of gravitational wave (GW) events detected by LIGO/VIRGO. Several recent studies have shown that when these binaries form, they are likely to be highly eccentric and retrograde. However, many uncertainties remain concerning the orbital evolution of these binaries as they either inspiral toward merger or disassociate. Previous hydrodynamical simulations exploring their orbital evolution have been predominantly two-dimensional or have been restricted to binaries on nearly circular orbits. We present the first high-resolution, three-dimensional local shearing-box simulations of both prograde and retrograde eccentric BBHs embedded in AGN disks. We find that retrograde binaries shrink several times faster than their prograde counterparts and exhibit significant orbital eccentricity growth, the rate of which monotonically increases with binary eccentricity. Our results suggest that retrograde binaries may experience runaway orbital eccentricity growth, which may bring them close enough together at pericenter for GW emission to drive them to coalescence. Although their eccentricity is damped, prograde binaries shrink much faster than their orbital eccentricity decays, suggesting they should remain modestly eccentric as they contract toward merger. Finally, binary precession driven by the AGN disk may dominate over precession induced by the supermassive black hole depending on the binary accretion rate and its location in the AGN disk, which can subdue the evection resonance and von Ziepel–Lidov–Kozai cycles.

Eccentric orbits may enhance the habitability of Earth-like exoplanets

ABSTRACT The detection and characterization of Earth-like planets around Sun-like stars is an important goal of exoplanetary research, given their promise for hosting potentially habitable conditions. Key orbital parameters, such as eccentricity, can influence a planet’s climate response and, as a consequence, affect its potential habitability. Utilizing the Earth System Model – the Whole Atmosphere Community Climate Model (WACCM6), we simulated Earth-like exoplanets with two different orbital parameters: one circular ($e = 0$) and another highly eccentric ($e = 0.4$), both with zero obliquity but fixing the annual mean insolation. The highly eccentric case exhibits a 1.9 K warmer surface temperature due to lower surface and cloud albedo and a weaker longwave cloud forcing. Exploring the annual global mean climate difference, we analysed latitudinal and seasonal variations in hydrological cycle variables, such as sea ice, land snow, and clouds. Land habitability metrics based on temperature and precipitation reveal that the $e=0.4$ case has over 25 per cent more habitable land area for more than 80 per cent of its orbit, compared with the $e=0$ case. Additionally, the global circulation pattern shifts from a three-cell to a two-cell system in the $e=0.4$ case, expanding the Hadley cell to higher latitudes, enhancing meridional latent heat transport, and improving land habitability at higher latitudes. Our study suggests that Earth-like exoplanets with high eccentricity orbiting Sun-like stars may have greater land habitability than their circular counterparts, due to seasonally warmer surface temperatures and more evenly distributed precipitation over land.

Measuring eccentricity and gas-induced perturbation from gravitational waves of LISA massive black hole binaries

ABSTRACT We assess the possibility of detecting both eccentricity and gas effects (migration and accretion) in the gravitational wave (GW) signal from LISA massive black hole binaries at redshift $z=1$. Gas induces a phase correction to the GW signal with an effective amplitude ($C_{\rm g}$) and a semimajor axis dependence (assumed to follow a power-law with slope $n_{\rm g}$). We use a complete model of the LISA response and employ a gas-corrected post-Newtonian inspiral-only waveform model TaylorF2Ecc. By using the Fisher formalism and Bayesian inference, we constrain $C_{\rm g}$ together with the initial eccentricity $e_0$, the total redshifted mass $M_z$, the primary-to-secondary mass ratio q, the dimensionless spins $\chi _{1,2}$ of both component BHs, and the time of coalescence $t_c$. We find that simultaneously constraining $C_{\rm g}$ and $e_0$ leads to worse constraints on both parameters with respect to when considered individually. For a standard thin viscous accretion disc around $M_z=10^5~{\rm M}_{\odot }$, $q=8$, $\chi _{1,2}=0.9$, and $t_c=4$ years MBHB, we can confidently measure (with a relative error of $\lt 50$ per cent) an Eddington ratio ${\rm f}_{\rm Edd}\sim 0.1$ for a circular binary and ${\rm f}_{\rm Edd}\sim 1$ for an eccentric system assuming $\mathcal {O}(10)$ stronger gas torque near-merger than at the currently explored much-wider binary separations. The minimum measurable eccentricity is $e_0\gtrsim 10^{-2.75}$ in vacuum and $e_0\gtrsim 10^{-2}$ in gas. A weak environmental perturbation (${\rm f}_{\rm Edd}\lesssim 1$) to a circular binary can be mimicked by an orbital eccentricity during inspiral, implying that an electromagnetic counterpart would be required to confirm the presence of an accretion disc.

A hot-Jupiter progenitor on a super-eccentric retrograde orbit.

Giant exoplanets orbiting close to their host stars are unlikely to have formed in their present configurations1. These 'hotJupiter' planets are instead thought to have migrated inward from beyond the ice line and several viable migration channels have been proposed, including eccentricity excitation through angular-momentum exchange with a third body followed by tidally driven orbital circularization2,3. The discovery of the extremely eccentric (e = 0.93) giant exoplanet HD 80606 b (ref. 4) provided observational evidence that hot Jupiters may have formed through this high-eccentricity tidal-migration pathway5. However, no similar hot-Jupiter progenitors have been found and simulations predict that one factor affecting the efficacy of this mechanism is exoplanet mass, as low-mass planets are more likely to be tidally disrupted during periastron passage6-8. Here we present spectroscopic and photometric observations of TIC 241249530 b, a high-mass, transiting warm Jupiter with an extreme orbital eccentricity of e = 0.94. The orbit of TIC 241249530 b is consistent with a history of eccentricity oscillations and a future tidal circularization trajectory. Our analysis of the mass and eccentricity distributions of the transiting-warm-Jupiter population further reveals a correlation between high mass and high eccentricity.