Search

A large fraction of white dwarfs are accreting or have recently accreted rocky material from their planetary systems, thereby polluting their atmospheres with elements heavier than helium. In recent years, the quest for mechanisms that can deliver planetesimals to the immediate vicinity of their central white dwarfs has stimulated a flurry of modelling efforts. The observed time evolution of the accretion rates of white dwarfs through their multi-Gyr lifetime is a crucial test for dynamical models of evolved planetary systems. Recent studies of cool white dwarfs samples have identified a significant decrease of the mass accretion rates of cool, old white dwarfs over Gyr timescales. Here, we revisit those results using updated white dwarf models and larger samples of old polluted H- and He-atmosphere white dwarfs. We find no compelling evidence for a strong decrease of their time-averaged mass accretion rates for cooling times between 1 and 8 Gyrs. Over this period, the mass accretion rates decrease by no more than a factor of the order of 10, which is one order of magnitude smaller than the decay rate found in recent works. Our results require mechanisms that can efficiently and consistently deliver planetesimals inside the Roche radius of white dwarfs over at least 8 Gyrs., 9 pages, 6 figures, 1 table. Accepted for publication in MNRAS

White dwarfs containing orbiting planetesimals or their debris represent crucial benchmarks by which theoretical investigations of post-main-sequence planetary systems may be calibrated. The photometric transit signatures of likely planetary debris in the ZTF J0139+5245 white dwarf system has an orbital period of about 110 days. An asteroid which breaks up to produce this debris may spin itself to destruction through repeated close encounters with the star without entering its Roche radius and without influence from the white dwarf's luminosity. Here, we place coupled constraints on the orbital pericentre ($q$) and the ratio ($\beta$) of the middle to longest semiaxes of a triaxial asteroid which disrupts outside of this white dwarf's Roche radius ($r_{\rm Roche}$) soon after attaining its 110-day orbit. We find that disruption within tens of years is likely when $\beta \lesssim 0.6$ and $q\approx 1.0-2.0r_{\rm Roche}$, and when $\beta \lesssim 0.2$ out to $q\approx 2.5r_{\rm Roche}$. Analysing the longer-timescale disruption of triaxial asteroids around ZTF J0139+5245 is desirable but may require either an analytical approach relying on ergodic theory or novel numerical techniques., Comment: Accepted for publication in MNRAS

We perform simulations of giant impacts onto the young Uranus using smoothed particle hydrodynamics (SPH) with over 100 million particles. This 100--1000$\times$ improvement in particle number reveals that simulations with below 10^7 particles fail to converge on even bulk properties like the post-impact rotation period, or on the detailed erosion of the atmosphere. Higher resolutions appear to determine these large-scale results reliably, but even 10^8 particles may not be sufficient to study the detailed composition of the debris -- finding that almost an order of magnitude more rock is ejected beyond the Roche radius than with 10^5 particles. We present two software developments that enable this increase in the feasible number of particles. First, we present an algorithm to place any number of particles in a spherical shell such that they all have an SPH density within 1% of the desired value. Particles in model planets built from these nested shells have a root-mean-squared velocity below 1% of the escape speed, which avoids the need for long precursor simulations to produce relaxed initial conditions. Second, we develop the hydrodynamics code SWIFT for planetary simulations. SWIFT uses task-based parallelism and other modern algorithmic approaches to take full advantage of contemporary supercomputer architectures. Both the particle placement code and SWIFT are publicly released., Fixed typo in Appx B equation and updated urls

The origin of the Martian moons, Phobos and Deimos, remains elusive. While the morphology and their cratered surfaces suggest an asteroidal origin1–3, capture has been questioned because of potential dynamical difficulties in achieving the current near-circular, near-equatorial orbits4,5. To circumvent this, in situ formation models have been proposed as alternatives6–9. Yet, explaining the present location of the moons on opposite sides of the synchronous radius, their small sizes and apparent compositional differences with Mars2 has proved challenging. Here, we combine geophysical and tidal-evolution modelling of a Mars–satellite system to propose that Phobos and Deimos originated from disintegration of a common progenitor that was possibly formed in situ. We show that tidal dissipation within a Mars–satellite system, enhanced by the physical libration of the satellite, circularizes the post-disrupted eccentric orbits in

We consider the geodesic deviation equation, describing the relative accelerations of nearby particles, and the Raychaudhuri equation, giving the evolution of the kinematical quantities associated with deformations (expansion, shear and rotation) in the Weyl type $f(Q,T)$ gravity, in which the nonmetricity $Q$ is represented in the standard Weyl form, fully determined by the Weyl vector, while $T$ represents the trace of the matter energy-momentum tensor. The effects of the Weyl geometry and of the extra force induced by the nonmetricity-matter coupling are explicitly taken into account. The Newtonian limit of the theory is investigated, and the generalized Poisson equation, containing correction terms coming from the Weyl geometry, and from the geometry matter coupling, is derived. As a physical application of the geodesic deviation equation the modifications of the tidal forces, due to the nonmetricity-matter coupling, are obtained in the weak field approximation. The tidal motion of test particles is directly influenced by the gradients of the extra force, and of the Weyl vector. As a concrete astrophysical example we obtain the expression of the Roche limit (the orbital distance at which a satellite begins to be tidally torn apart by the body it orbits) in the Weyl type $f(Q,T)$ gravity., Comment: 18 pages, no figures, accepted for publication in EPJC; minor corrections; references added

We numerically investigate how an asteroid's elongation controls the sensitivity of its surface to tidal effects during a distant planetary encounter beyond the Roche limit. We analyze the surface slope and its variation by considering the shape elongation, as well as the spin period and orbital conditions. A more elongated asteroid tends to have a higher slope variation, while there may not be a monotonic increase of the total area having such a variation, 12 pages, 1 figures, to be published in Icarus

Tidal dissipation due to turbulent viscosity in the convective regions of giant stars plays an important role in shaping the orbits of pre-common envelope systems. Such systems are possible sources of transients and close compact binary systems that will eventually merge and produce detectable gravitational wave signals. Most previous studies of the onset of common envelope episodes have focused on circular orbits and synchronously rotating donor stars under the assumption that tidal dissipation can quickly spin up the primary and circularize the orbit before the binary reaches Roche-lobe overflow (RLO). We test this assumption by coupling numerical models of the post main sequence stellar evolution of massive-stars with the model for tidal dissipation in convective envelopes developed in Vick & Lai (2020) $-$ a tidal model that is accurate even for highly eccentric orbits with small pericentre distances. We find that, in many cases, tidal dissipation does not circularize the orbit before RLO. For a $10~M_\odot$ ($15~M_\odot$) primary star interacting with a $1.4~M_\odot$ companion, systems with pericentre distances within 3 AU (6 AU) when the primary leaves the main sequence will retain the initial orbital eccentricity when the primary grows to the Roche radius. Even in systems that tidally circularize before RLO, the donor star may be rotating subsynchronously at the onset of mass transfer. Our results demonstrate that some possible precursors to double neutron star systems are likely eccentric at the Roche radius. The effects of pre-common envelope eccentricity on the resulting compact binary merit further study., 15 pages, 9 figures, Accepted to MNRAS

The Moon migrated to $r_{\leftmoon}\simeq3.8\times10^{10}$ cm over a characteristic time $r/v=10^{10}$ Gyr by tidal interaction with the Earth's oceans at a present velocity of $v=3.8$ cm yr$^{-1}$. We derive scaling of global dissipation that covers the entire history over the past 4.52 Gyr. Off-resonance tidal interactions at relatively short tidal periods in the past reveal the need for scaling {with amplitude}. The global properties of the complex spatio-temporal dynamics and dissipation in broad spectrum ocean waves is modeled by damping $\epsilon = h F/(2Q_0)$, where $h$ is the tidal wave amplitude, $F$ is the tidal frequency, and $Q_0$ is the $Q$-factor at the present time. It satisfies $Q_0\simeq 14$ for consistency of migration time and age of the Moon consistent with observations for a near-resonance state today. It shows a startingly fast eviction of the Moon from an unstable near-synchronous orbit close to the Roche limit, probably in a protolunar disk. Rapid spin down of the Earth from an intial $\sim30\%$ of break-up by the Moon favored early formation of a clement global climate. Our theory suggests moons may be similarly advantageous to potentially habitable exoplanets., Comment: 8 pages, 5 figures

The class of ultra-hot Jupiters comprises giant exoplanets undergoing intense irradiation from their host stars. They have proved to be a particularly interesting population for their orbital and atmospheric properties. One such planet, WASP-121~b, is in a highly misaligned orbit close to its Roche limit, and its atmosphere exhibits a thermal inversion. These properties make WASP-121~b an interesting target for additional atmospheric characterization. In this paper, we present analysis of archival high-resolution optical spectra obtained during transits of WASP-121~b. Artifacts from the Rossiter-McLaughlin effect and Center-to-Limb Variation are deemed negligible. However, we discuss scenarios where these effects warrant more careful treatment by modeling the WASP-121 system and varying its properties. We report a new detection of atmospheric absorption from H$\alpha$ in the planet with a transit depth of $1.87\pm0.11\%$. We further confirm a previous detection of the Na I doublet, and report a new detection of Fe I via cross-correlation with a model template. We attribute the H$\alpha$ absorption to an extended Hydrogen atmosphere, potentially undergoing escape, and the Fe I to equilibrium chemistry at the planetary photosphere. These detections help to constrain the composition and chemical processes in the atmosphere of WASP-121~b., Comment: 15 pages, 14 figures, submitted to MNRAS after addressing reviewer comments

Since 2013, dense and narrow rings are known around the small Centaur object Chariklo and the dwarf planet Haumea. Dense material has also been detected around the Centaur Chiron, although its nature is debated. This is the first time ever that rings are observed elsewhere than around the giant planets, suggesting that those features are more common than previously thought. The origins of those rings remain unclear. In particular, it is not known if the same generic process can explain the presence of material around Chariklo, Chiron, Haumea, or if each object has a very different history. Nonetheless, a specific aspect of small bodies is that they may possess a nonaxisymmetric shape (topographic features and/or elongation) that are essentially absent in giant planets. This creates strong resonances between the spin rate of the object and the mean motion of ring particles. In particular, Lindblad-type resonances tend to clear the region around the corotation (or synchronous) orbit, where the particles orbital period matches that of the body. Whatever the origin of the ring is, modest topographic features or elongations of Chariklo and Haumea explain why their rings should be found beyond the outermost 1/2 resonance, where the particles complete one revolution, while the body completes two rotations. Comparison of the resonant locations relative to the Roche limit of the body shows that fast rotators are favored for being surrounded by rings. We discuss in more detail the phase portraits of the 1/2 and 1/3 resonances, and the consequences of a ring presence on satellite formation.

The early thermal evolution of Moon has been numerically simulated to understand the magnitude of the impact induced heating and the initially stored thermal energy of the accreting Moonlets. The main objective of the present study is to understand the nature of processes leading to core-mantle differentiation and the production and cooling of the initial convective magma ocean. The accretion of Moon was commenced over a timescale of 100 years after the giant impact event around 30-100 million years in the early solar system. We studied the dependence of the planetary processes on the impact scenarios, the initial average temperature of the accreting moonlets and the size of the protoMoon that accreted rapidly beyond the Roche limit within the initial one year after the giant impact. The simulations indicate that the accreting Moonlets should have a minimum initial averaged temperature around 1600 K. The impacts would provide additional thermal energy. The initial thermal state of the moonlets depends upon the environment prevailing within the Roche limit that experienced episodes of extensive vaporization and re-condensation of silicates. The initial convective magma ocean of depth more than 1000 km is produced in the majority of simulations along with the global core-mantle differentiation in case the melt percolation of the molten metal through porous flow from bulk silicates was not the major mode of core-mantle differentiation. The possibility of shallow magma oceans cannot be ruled out in the presence of the porous flow. Our simulations indicate the core-mantle differentiation within the initial $10^2$-$10^3$ years of the Moon accretion. The majority of the convective magma ocean cooled down for crystallization within the initial $10^3$-$10^4$ years., Comment: Accepted in Meteoritics & Planetary Science, 39 pages, 5 figures

Binary systems undergoing unstable Roche Lobe overflow spill gas into their circumbinary environment as their orbits decay toward coalescence. In this paper, we use a suite of hydrodynamic models of coalescing binaries involving an extended donor and a more compact accretor. We focus on the period of unstable Roche Lobe overflow that ends as the accretor plunges within the envelope of the donor at the onset of a common envelope phase. During this stage, mass is removed from the donor and flung into the circumbinary environment. Across a wide range of binary mass ratios, we find that the mass expelled as the separation decreases from the Roche limit to the donor's original radius is of the order of 25% of the accretor's mass. We study the kinematics of this ejecta and its dependencies on binary properties and find that it assembles into a toroidal circumbinary distribution. These circumbinary tori have approximately constant specific angular momentum due to momentum transport by spiral shocks launched from the orbiting binary. We show that an analytic model with these torus properties captures many of the main features of the azimuthally-averaged profiles of our hydrodynamic simulations. Our results, in particular the simple relationship between accretor mass and expelled mass and its spatial distribution, may be useful in interpreting stellar coalescence transients like luminous red novae, and in initializing hydrodynamic simulations of the subsequent common envelope phase., Comment: The Astrophysical Journal, Volume 895, Issue 1, id.29. Software to reproduce these results availible at: https://github.com/morganemacleod/PreCEMassLoss, software to model pre-CE unstable mass transfer and circumbinary material based on simulation results availible at: https://github.com/morganemacleod/RLOF

Original content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI., From CCD observations carried out with different telescopes, we present short-term photometric measurements of the large trans-Neptunian object Varuna in 10 epochs, spanning around 19 years. We observe that the amplitude of the rotational light curve has changed considerably during this period of time from 0.41 to 0.55 mag. In order to explain this variation, we constructed a model in which Varuna has a simple triaxial shape, assuming that the main effect comes from the change of the aspect angle as seen from Earth, due to Varuna's orbital motion in the 19 year time span. The best fits to the data correspond to a family of solutions with axial ratios b/a between 0.56 and 0.60. This constrains the pole orientation in two different ranges of solutions presented here as maps. Apart from the remarkable variation of the amplitude, we have detected changes in the overall shape of the rotational light curve over shorter timescales. After the analysis of the periodogram of the residuals to a 6.343572 hr double-peaked rotational light-curve fit, we find a clear additional periodicity. We propose that these changes in the rotational light-curve shape are due to a large and close-in satellite whose rotation induces the additional periodicity. The peak-to-valley amplitude of this oscillation is in the order of 0.04 mag. We estimate that the satellite orbits Varuna with a period of 11.9819 hr (or 23.9638 hr), assuming that the satellite is tidally locked, at a distance of ∼1300 km (or ∼2000 km) from Varuna, outside the Roche limit.© 2019. The American Astronomical Society.., E.F.V. acknowledges UCF 2017 Preeminent Postdoctoral Program (P3). Part of the research leading to these results has received funding from the European Unions Horizon 2020 Research and Innovation Programme, under Grant Agreement No. 687378 (SBNAF). We would like to acknowledge financial support by the Spanish grant AYA-2017-84637-R and the financial support from the State Agency for Research of the Spanish MCIU through the >Center of Excellence Severo Ochoa> award for the Instituto de Astrofisica de Andalucia (SEV-2017-0709). This work is partially based on observations made with the Italian Telescopio Nazionale Galileo (TNG) operated on the island of La Palma by the Fundacion Galileo Galilei of the INAF (Istituto Nazionale di Astrofisica) at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias. We thank the director of the Telescopio Nazionale Galileo for allocation of Director's Discretionary Time. We acknowledge Raul Carballo-Rubio for providing useful comments that helped improve this manuscript and Pedro Gutierrez for helpful discussions.

The discovery of the intact minor planet embedded in the debris disc orbiting SDSS J1228+1040 raises questions about the dynamical history of the system. Further, the recent passage of the potentially interstellar object 1I/'Oumuamua within the solar system has re-ignited interest in minor body flux through exoplanetary systems. Here, we utilize the new analytical formalism from Grishin, Perets and Avni (2019) to estimate the rate at which the gaseous components of typical white dwarf discs trap an exo-planetesimal. We compare the types of captured orbits which arise from planetesimals originating from the interstellar medium, exo-Kuiper belts and exo-Oort clouds. We find that the rate of interstellar medium injection is negligible, whereas capture of both exo-Kuiper and exo-Oort cloud planetesimals is viable, but strongly size-dependent. For a gaseous disc which extends much beyond its Roche limit, capture is more probable than disruption at the Roche limit. We find that the capture probability linearly increases with the radial extent of the disc. Even in systems without minor planets, capture of smaller bodies will change the disc size distribution and potentially its temporal variability. Our formalism is general enough to be applied to future discoveries of embedded planetesimals in white dwarf debris discs., Accepted for publication in MNRAS

Observational evidence of white dwarf planetary systems is dominated by the remains of exo-asteroids through accreted metals, debris discs, and orbiting planetesimals. However, exo-planets in these systems play crucial roles as perturbing agents, and can themselves be perturbed close to the white dwarf Roche radius. Here, we illustrate a procedure for computing the tidal interaction between a white dwarf and a near-spherical solid planet. This method determines the planet's inward and/or outward drift, and whether the planet will reach the Roche radius and be destroyed. We avoid constant tidal lag formulations and instead employ the self-consistent secular Darwin-Kaula expansions from Bou�� & Efroimsky (2019), which feature an arbitrary frequency dependence on the quality functions. We adopt wide ranges of dynamic viscosities and spin rates for the planet in order to straddle many possible outcomes, and provide a foundation for the future study of individual systems with known or assumed rheologies. We find that: (i) massive Super-Earths are destroyed more readily than minor planets (such as the ones orbiting WD 1145+017 and SDSS J1228+1040), (ii) low-viscosity planets are destroyed more easily than high-viscosity planets, and (iii) the boundary between survival and destruction is likely to be fractal and chaotic., Accepted for publication in MNRAS

We present follow-up photometry and spectroscopy of ZTF J0328$-$1219 strengthening its status as a white dwarf exhibiting transiting planetary debris. Using TESS and Zwicky Transient Facility photometry, along with follow-up high speed photometry from various observatories, we find evidence for two significant periods of variability at 9.937 and 11.2 hr. We interpret these as most likely the orbital periods of different debris clumps. Changes in the detailed dip structures within the light curves are observed on nightly, weekly, and monthly timescales, reminiscent of the dynamic behavior observed in the first white dwarf discovered to harbor a disintegrating asteroid, WD 1145+017. We fit previously published spectroscopy along with broadband photometry to obtain new atmospheric parameters for the white dwarf, with $M_{\star} = 0.731 \pm 0.023\,M_{\odot}$, $T_{\mathrm{eff}} = 7630 \pm 140\,$K, and $\mathrm{[Ca/He]}=-9.55\pm0.12$. With new high-resolution spectroscopy, we detect prominent and narrow Na D absorption features likely of circumstellar origin, with velocities $21.4\pm1.0$ km s$^{-1}$ blue-shifted relative to atmospheric lines. We attribute the periodically modulated photometric signal to dusty effluents from small orbiting bodies such as asteroids or comets, but are unable to identify the most likely material that is being sublimated, or otherwise ejected, as the environmental temperatures range from roughly 400K to 600K., 24 pages, 10 figures, 9 tables, accepted to ApJ

A model of the viscous evolution of a two-phase, vapor/melt protolunar disk is described. Droplets condense from the vapor and "rain-out", forming a stratified structure with a mid-plane magma layer surrounded by a vapor reservoir. The magma layer is gravitationally unstable, but material interior to the Roche distance cannot fragment, and instead develops an effective viscosity. However, magma flowing across the Roche limit can fragment and accrete into moonlets, while magma spreading inward is accreted by Earth. As mass leaves the melt layer, it is replenished by vapor condensation, leading to a quasi steady-state (QSS). The layer's mass is maintained at ~14% of a lunar mass, and replaced every ~ 3 yrs. The vapor atmosphere steadily decreases, and once exhausted, the disk would cool below condensation temperature and spread as a time-dependent disk until it too is exhausted. The timescale of the QSS is regulated by the disk's ability to radiate all the released latent heat plus viscous dissipation energy to space at the photospheric temperature, i.e., ~2000 K for silicon phase-equilibrium. For the protolunar disk, latent heat dominates viscous heating, and the QSS for a ~2 lunar mass disk lasts for ~ 50 yrs. For comparison, a hypothetical water/steam disk orbiting an ice giant planet in which viscous heating dominates is also modeled. The photospheric temperature is closer to ~373 K, and the replacement time of the water layer is ~125 yrs. Finally, disks confined by external torques at either the planet-disk boundary or at its outer edge are also briefly examined.

Post-main-sequence planetary science has been galvanised by the striking variability, depth and shape of the photometric transit curves due to objects orbiting white dwarf WD 1145+017, a star which also hosts a dusty debris disc and circumstellar gas, and displays strong metal atmospheric pollution. However, the physical properties of the likely asteroid which is discharging disintegrating fragments remain largely unconstrained from the observations. This process has not yet been modelled numerically. Here, we use the N-body code PKDGRAV to compute dissipation properties for asteroids of different spins, densities, masses, and eccentricities. We simulate both homogeneous and differentiated asteroids, for up to two years, and find that the disruption timescale is strongly dependent on density and eccentricity, but weakly dependent on mass and spin. We find that primarily rocky differentiated bodies with moderate (~3-4 g/cm^3) bulk densities on near-circular (e, Comment: Accepted in MNRAS. Movies here!: http://www.star.bris.ac.uk/pcarter/WD1145_asteroid_disruption/

The Martian satellite Phobos is criss-crossed by linear grooves and crater chains whose origin is unexplained. Anomalous grooves are relatively young, and crosscut tidally predicted stress fields as Phobos spirals towards Mars. Here we report strong correspondence between these anomalous features and reaccretion patterns of sesquinary ejecta from impacts on Phobos. Escaping ejecta persistently imprint Phobos with linear, low-velocity crater chains (catenae) that match the geometry and morphology of prominent features that do not fit the tidal model. We prove that these cannot be older than Phobos' current orbit inside Mars' Roche limit. Distinctive reimpact patterns allow sesquinary craters to be traced back to their source, for the first time across any planetary body, creating a novel way to probe planetary surface characteristics. For example, we show that catena-producing craters likely formed in the gravity regime, providing constraints on the ejecta velocity field and knowledge of source crater material properties., The origin of the numerous linear grooves and craters that litter the Martian moon Phobos' surface remains enigmatic. Here, by modelling low-velocity escaping ejecta from impacts to Phobos, the authors show that several of these chains can be explained by reimpacting sesquinary ejecta shortly after ejection.

All four giant planets are encircled by distinctive systems of rings and small, inner satellites. These all reside within or near their central planet's Roche limit, the rough boundary within which bodies held together by self-gravity will be disrupted by tidal forces. However, the similarities of the four ring-moon systems end here; in most other regards, they are remarkably diverse. We study these systems for three key reasons: (1) for the information they reveal about the properties, history and ongoing evolution of the planetary systems of which they are a part; (2) as dynamical analogues for other astrophysical systems such as protoplanetary disks; and (3) for the wealth of fascinating properties and origin scenarios that make them worthy of study in their own right. The inner Uranus system is characterized by 10 narrow rings, some quite dense, as well as a variety of more tenuous structures. These are accompanied by 13 known moons all orbiting interior to Miranda. Nine of these, Bianca through Perdita, comprise the most densely packed set of moons in the solar system, with orbits so close that their interactions appear to drive chaos over time scales approximately 106years. Neptune has five named rings, all optically thin, interleaved with seven inner moons. The most notable feature is a set of arcs embedded within the Adams ring; two of these arcs have been stable for time scales of decades.This article is part of a discussion meeting issue ‘Future exploration of ice giant systems’.

Regions of Hill stability were determined using the full three-body theory for the planetary candidates that were found by the Kepler mission to lie in the Habitable Zone of the host star with radii larger than 3 Earth radii. The values obtained give a more accurate assessment of the sphere of influence of a planet than does the Hill radius using the restricted three-body expression. The values given by the two estimates for circular orbits are compared. The variation in the full three-body Hill regions were examined for a variety of possible moon masses and orbital eccentricities and inclinations. It was found that the critical moon planet separation decreased with increased moon orbital eccentricity and inclination but to a far greater extent when the planetary orbital eccentricity increased suggesting they would lie inside the Roche limit and are unlikely to have formed or exist now. The angular separation of the moon and planet and the radial velocity amplitudes of the planets likely to be observed are also affected. This is important for future observations and missions. It was also found that the giant planets in the habitable zone lie in a narrow band that can be fitted by a linear regression relation between log10(Planet luminosity) and log10 (Planet temperature) with a slope of 3.34.

We report the discovery of a white dwarf exhibiting deep, irregularly shaped transits, indicative of circumstellar planetary debris. Using Zwicky Transient Facility DR2 photometry of ZTF$\,$J013906.17+524536.89 and follow-up observations from the Las Cumbres Observatory, we identify multiple transit events that recur every ${\approx}\,107.2\,$d, much longer than the $4.5{-}4.9\,$h orbital periods observed in WD$\,$1145+017, the only other white dwarf known with transiting planetary debris. The transits vary in both depth and duration, lasting $15{-}25\,$d and reaching $20{-}45\,\%$ dips in flux. Optical spectra reveal strong Balmer lines, identifying the white dwarf as a DA with $T_{\mathrm{eff}}=10{,}530\pm140\,\mathrm{K}$ and $\log(g)=7.86\pm0.06$. A $\mathrm{Ca\,II\,K}$ absorption feature is present in all spectra both in and out of transit. Spectra obtained during one night at roughly $15\,\%$ transit depth show increased $\mathrm{Ca\,II\,K}$ absorption with a model atmospheric fit suggesting $[\mathrm{Ca/H}]=-4.6\pm0.3$, whereas spectra taken on three nights out of transit have $[\mathrm{Ca/H}]$ of -5.5, -5.3, and -4.9 with similar uncertainties. While the $\mathrm{Ca\,II\,K}$ line strength varies by only 2-sigma, we consider a predominantly interstellar origin for Ca absorption unlikely. We suggest a larger column density of circumstellar metallic gas along the line of site or increased accretion of material onto the white dwarf's surface are responsible for the Ca absorption, but further spectroscopic studies are required. In addition, high-speed time series photometry out of transit reveals variability with periods of 900 and 1030$\,$s, consistent with ZZ Ceti pulsations., Submitted to ApJ, Accepted May 26, 2020. 5 Figures, 1 Table

We consider the impact of electromagnetic induction and Ohmic heating on a conducting planetary object that orbits a magnetic star. Power dissipated as heat saps orbital energy. If this heat is trapped by an insulating crust or mantle, interior temperatures increase substantially. We provide a quantitative description of this behavior and discuss the astrophysical scenarios in which it might occur. Magnetic fields around some main-sequence stars and white dwarfs are strong enough to cause the decay of close-in orbits of asteroids and dwarf planets, drawing them through the Roche limit on Myr time scales. We confirm that Ohmic heating around neutron stars is driven by the rotation of the stellar magnetic dipole, not orbital dynamics. In any case, heating can raise interior temperatures of asteroids or dwarf planets on close-in orbits to well above liquidus. Hot material escaping to the surface may lead to volcanic ejections that can obscure the host star (as in the light curve of KIC 8462852) and pollute its atmosphere (as observed with metal-rich white dwarfs). We speculate that mixing of a volatile-rich mantle or crust with material from an induction-heated core may lead to an explosion that could destroy the asteroid prior to tidal break-up., 36 pages, 5 figures, 2 tables, ApJ, accepted

Most massive stars, if not all, are in binary configuration or higher multiples. These massive stars undergo supernova explosions and end their lives as either black holes or neutron stars. Recent observations have suggested that neutron stars and perhaps even black holes receive large velocity kicks at birth. Such natal kicks and the sudden mass loss can significantly alter the orbital configuration of the system. Here we derive general analytical expressions that describe the effects of natal kicks in binaries on hierarchical triple systems. We explore several proof-of-concept applications such as black hole and neutron stars binaries and X-ray binaries with either stellar or Supermassive Black Hole (SMBH) companions on a wide orbit. Kicks can disrupt the hierarchical configuration, although it is harder to escape the potential well of an SMBH. Some binary systems do escape the SMBH system resulting in hyper-velocity binary system. Furthermore, kicks can result in increasing or decreasing the orbital separations. Decreasing the orbital separation may have significant consequences in these astrophysical systems. For example, shrinking the separation post-supernova kick can lead to the shrinking of an inner compact binary that then may merge via Gravitational Wave (GW) emission. This process yields a supernova that is shortly followed by a possible GW-LIGO event. Interestingly, we find that the natal kick can result in shrinking the outer orbit, and the binary may cross the tertiary Roche limit, breaking up the inner binary. Thus, in the case of SMBH companion, this process can lead to either a tidal disruption event or a GW-LISA detection event with a supernova precursor.

Observed planetary debris in white dwarf atmospheres predominately originate from the destruction of small bodies on highly eccentric ($>0.99$) orbits. Despite their importance, these minor planets have coupled physical and orbital evolution which has remained largely unexplored. Here, we present a novel approach for estimating the influence of fast chaotic rotation on the orbital evolution of high-eccentricity triaxial asteroids, and formally characterize the propagation of their angular rotation velocities and orbital elements as random time processes. By employing the impulse approximation, we demonstrate that the violent gravitational interactions during periastron passages transfer energy between the orbit and asteroid's rotation. If the distribution of spin impulses were symmetric around zero, then the net result would be a secular decrease of the semimajor axis and a further increase of the eccentricity. We find evidence, however, that the chaotic rotation may be self-regulated in such a manner that these effects are reduced or nullified. We discover that asteroids on highly eccentric orbits can break themselves apart --- in a type of YORP-less rotational fission --- without actually entering the Roche radius, with potentially significant consequences for the distribution of debris and energy requirements for gravitational scattering in metal-polluted white dwarf planetary systems. This mechanism provides a steady stream of material impacting a white dwarf without rapidly depleting the number of small bodies in the stellar system., Revised manuscript resubmitted to ApJ after referee report

We analyze the orbital motion of two natural satellites initially in co-orbital configuration with a third (guiding) satellite embedded into a circumplanetary gas disc and undergoing tidal interactions with the central planet. By solving the exact equations of motion, including the dissipative effects and the mutual gravitational perturbations, we investigate the configuration of the system soon after the guiding satellite is disrupted when crossing the Roche limit with the central planet during its orbital decay. The application is done for the system composed of Saturn and the two small co-orbital satellites Janus and Epimetheus. We perform a large number of numerical simulations, varying the mass of the hypothetical guiding satellite, Saturn's tidal quality factor and the initial configuration and the masses of the small satellites. Analyzing the orbital configurations of Janus and Epimetheus soon after the disruption of the guiding satellite, we obtain that the mutual horseshoe co-orbital orbits appear as a natural outcome of the numerical simulations., Comment: Accepted for publication in MNRAS

Infrared excesses due to dusty disks have been observed orbiting white dwarfs with effective temperatures between 7200 K and 25000 K, suggesting that the rate of tidal disruption of minor bodies massive enough to create a coherent disk declines sharply beyond 1~Gyr after white dwarf formation. We report the discovery that the candidate white dwarf LSPM J0207+3331, via the Backyard Worlds: Planet 9 citizen science project and Keck Observatory follow-up spectroscopy, is hydrogen-dominated with a luminous compact disk (L$_{\rm IR}$/L$_{\star}$=14%) and an effective temperature nearly 1000K cooler than any known white dwarf with an infrared excess. The discovery of this object places the latest time for large scale tidal disruption events to occur at $\sim$3 Gyr past the formation of the host white dwarf, making new demands of dynamical models for planetesimal perturbation and disruption around post main sequence planetary systems. Curiously, the mid-IR photometry of the disk cannot be fully explained by a geometrically thin, optically thick dust disk as seen for other dusty white dwarfs, but requires a second ring of dust near the white dwarf's Roche radius. In the process of confirming this discovery, we found that careful measurements of WISE source positions can reveal when infrared excesses for white dwarfs are co-moving with their hosts, helping distinguish them from confusion noise., Comment: 7 pages, Published in ApJ Letters

A gas giant planet which survives the giant branch stages of evolution at a distance of many au and then is subsequently perturbed sufficiently close to a white dwarf will experience orbital shrinkage and circularization due to star-planet tides. The circularization timescale, when combined with a known white dwarf cooling age, can place coupled constraints on the scattering epoch as well as the active tidal mechanisms. Here, we explore this coupling across the entire plausible parameter phase space by computing orbit shrinkage and potential self-disruption due to chaotic f-mode excitation and heating in planets on orbits with eccentricities near unity, followed by weakly dissipative equilibrium tides. We find that chaotic f-mode evolution activates only for orbital pericentres which are within twice the white dwarf Roche radius, and easily restructures or destroys ice giants but not gas giants. This type of internal thermal destruction provides an additional potential source of white dwarf metal pollution. Subsequent tidal evolution for the surviving planets is dominated by non-chaotic equilibrium and dynamical tides which may be well-constrained by observations of giant planets around white dwarfs at early cooling ages., Comment: Accepted for publication in MNRAS

Dense and narrow rings have been discovered recently around the small Centaur object Chariklo1 and the dwarf planet Haumea2, while being suspected around the Centaur Chiron3, although this point is debated4. They are the first rings observed in the Solar System elsewhere than around giant planets. In contrast to giant planets, gravitational fields of small bodies may exhibit large non-axisymmetric terms that create strong resonances between the spin of the object and the mean motion of ring particles. Here we show that modest topographic features or elongations of Chariklo and Haumea explain why their rings are relatively far away from the central body, when scaled to those of the giant planets5. Resonances actually clear on decadal timescales an initial collisional disk that straddles the corotation resonance (where the particles' mean motion matches the spin rate of the body). Quite generically, the disk material inside the corotation radius migrates onto the body, while the material outside the corotation radius is pushed outside the 1/2 resonance, where the particles complete one revolution while the body completes two rotations. Consequently, the existence of rings around non-axisymmetric bodies requires that the 1/2 resonance resides inside the Roche limit of the body, favouring faster rotators for being surrounded by rings. © 2018, The Author(s), under exclusive licence to Springer Nature Limited., The work leading to these results has received funding from the European Research Council under the European Community's H2020 2014-2020 ERC Grant Agreement No. 669416 'Lucky Star'. P.S.-S. acknowledges financial support by the European Union's Horizon 2020 Research and Innovation Programme, under Grant Agreement No. 687378 ('SBNAF').

We present 1"-resolution ALMA observations of the circumnuclear disk (CND) and the environment around SgrA*. The images unveil the presence of small spatial scale CO (J=3-2) molecular "cloudlets" within the central pc of the Milky Way, moving at high speeds, up to 300 km/s along the line-of-sight. The CO-emitting structures show intricate morphologies: extended and filamentary at high negative-velocities (v_LSR < -150 km/s), more localized and clumpy at extreme positive-velocities (v_LSR > +200 km/s). Based on the pencil-beam CO absorption spectrum toward SgrA* synchrotron emission, we also present evidence for a diffuse gas component producing absorption features at more extreme negative-velocities (v_LSR < -200 km/s). The CND shows a clumpy spatial distribution. Its motion requires a bundle of non-uniformly rotating streams of slightly different inclinations. The inferred gas density peaks are lower than the local Roche limit. This supports that CND molecular cores are transient. We apply the two standard orbit models, spirals vs. ellipses, invoked to explain the kinematics of the ionized gas streamers around SgrA*. The location and velocities of the CO cloudlets are inconsistent with the spiral model, and only two of them are consistent with the Keplerian ellipse model. Most cloudlets, however, show similar velocities that are incompatible with the motions of the ionized streamers or with gas bounded to the central gravity. We speculate that they are leftovers of more massive, tidally disrupted, clouds that fall into the cavity, or that they originate from instabilities in the inner rim of the CND and infall from there. Molecular cloudlets, all together with a mass of several 10 M_Sun, exist around SgrA*. Most of them must be short-lived: photoevaporated by the intense stellar radiation field, blown away by winds from massive stars, or disrupted by strong gravitational shears., Accepted for publication in A&A. 21 pages, 19 figures (includes corrections by language editor)

Using numerical simulations, we studied several coupled translational and rotational solutions of the two-finite-body problem with one spherical and one triaxial body. The aim was to investigate which types of orbits and planetary bodies could produce spin-induced orbital perturbations relevant enough to add to models dealing with other perturbations. To fully assess the strengths and consequences of this perturbation, we did not include any other perturbation even when a more realistic scenario would have required it. Interesting results concern planet–star mass ratios like a hot Jupiter or a super-Jupiter around a star like the Sun or the red dwarf Proxima Centauri. The short-period chaotic effect of the gravitational spin–orbit perturbation on highly eccentric orbits in the vicinity of the Roche limit can be a prominent feature. It should be taken into account when studying the tidal evolution of such a planet or its interactions with any companion in the neighborhood of the star.

Recent observations have revealed that the remnants of stellar-coalescence transients are bipolar. This raises the questions of how these bipolar morphologies arise and what they teach us about the mechanisms of mass ejection during stellar mergers and common-envelope phases. In this paper, we analyze hydrodynamic simulations of the lead-in to binary coalescence, a phase of unstable Roche lobe overflow that takes the binary from the Roche limit separation to the engulfment of the more compact accretor within the envelope of the extended donor. As mass transfer runs away at increasing rates, gas trails away from the binary. Contrary to previous expectations, early mass loss from the system remains bound to the binary and forms a circumbinary torus. Later ejecta, generated as the accretor grazes the surface of the donor, have very different morphologies and are unbound. These two components of mass loss from the binary interact as later, higher-velocity ejecta collide with the circumbinary torus formed by earlier mass loss. Unbound ejecta are redirected toward the poles, and escaping material creates a bipolar outflow. Our findings show that the transition from bound to unbound ejecta from coalescing binaries can explain the bipolar nature of their remnants, with implications for our understanding of the origin of bipolar remnants of stellar-coalescence transients and, perhaps, some preplanetary nebulae., ApJ accepted

Saturns F-ring is a unique, narrow ring that lies (radially) close to the tidally disruptive Roche limit of water ice for Saturn. Significant work has been done that shows it to be one of the most dynamic places in the Solar System. Aggregates that are fortunate enough to form constantly battle against the strong tidal forces of Saturn and the nearby moons Prometheus and Pandora, which act to gravitationally stir up ring material. Planetary rings are also known to radially spread. Therefore, as the F ring lies at the edge of the main rings, we investigate the effect of an outwardly migrated F ring and its interaction with Prometheus. An increase in the maximum number density of particles at the channel edges is observed with decreasing local tidal environment. Radial velocity dispersions are also observed to fall below the typical escape velocity of a 150m icy moonlet (, 25 pages, 8 figures

We perform a suite of smoothed particle hydrodynamics simulations to investigate in detail the results of a giant impact on the young Uranus. We study the internal structure, rotation rate, and atmospheric retention of the post-impact planet, as well as the composition of material ejected into orbit. Most of the material from the impactor's rocky core falls in to the core of the target. However, for higher angular momentum impacts, significant amounts become embedded anisotropically as lumps in the ice layer. Furthermore, most of the impactor's ice and energy is deposited in a hot, high-entropy shell at a radius of ~3 Earth radii. This could explain Uranus' observed lack of heat flow from the interior and be relevant for understanding its asymmetric magnetic field. We verify the results from the single previous study of lower resolution simulations that an impactor with a mass of at least 2 Earth masses can produce sufficiently rapid rotation in the post-impact Uranus for a range of angular momenta. At least 90% of the atmosphere remains bound to the final planet after the collision, but over half can be ejected beyond the Roche radius by a 2 or 3 Earth mass impactor. This atmospheric erosion peaks for intermediate impactor angular momenta (~3*10^36 kg m^2 s^-1). Rock is more efficiently placed into orbit and made available for satellite formation by 2 Earth mass impactors than 3 Earth mass ones, because it requires tidal disruption that is suppressed by the more massive impactors., Comment: Published in ApJ. 14 pages, 11 figures

We address the thermal history of the Earth after the Moon-forming impact, taking tidal heating and thermal blanketing by the atmosphere into account. The atmosphere sets an upper bound of ~100 W/m^2 on how quickly the Earth can cool. The liquid magma ocean cools over 2-10 Myrs, with longer times corresponding to high angular-momentum events. Tidal heating is focused mostly in mantle materials that are just beginning to freeze. The atmosphere's control over cooling sets up a negative feedback between viscosity-dependent tidal heating and temperature-dependent viscosity of the magma ocean. While the feedback holds, evolution of the Moon's orbit is limited by the modest radiative cooling rate of Earth's atmosphere. Orbital evolution is orders of magnitude slower than in conventional constant Q models, which promotes capture by resonances. The evection resonance is encountered early, when the Earth is molten. Capture by the evection resonance appears certain but unlikely to generate much eccentricity because it is encountered early when the Earth is molten and Q_Earth >> Q_Moon. Tidal dissipation in the Earth becomes more efficient (Q_Earth << Q_Moon) later when the Moon is between ~20 R_Earth and ~40 R_Earth. If lunar eccentricity grew great, this was when it did so, perhaps setting the table for some other process to leave its mark on the inclination of the Moon., 34 pages, 10 figures, preprint format