Aeolian-erosion is a destructive process which can erode small-size planetary objects through their interaction with a gaseous environment. Aeolian-erosion operates in a wide range of environments and under various conditions. Aeolian-erosion has been extensively explored in the context of geophysics in terrestrial planets. Here we show that aeolian-erosion of cobbles, boulders and small planetesimals in protoplanetary-discs can constitute a significant barrier for the early stages of planet formation. We use analytic calculations to show that under the conditions prevailing in protoplanetary-discs small bodies ($10-10^4 \rm{m}$) are highly susceptible to gas-drag aeolian-erosion. At this size-range aeolian-erosion can efficiently erode the planetesimals down to tens-cm size and quench any further growth of such small bodies. It thereby raises potential difficulties for channels suggested to alleviate the metre-size barrier. Nevertheless, the population of $\sim$decimetre-size cobbles resulting from aeolian-erosion might boost the growth of larger (>km size) planetesimals and planetary embryos through increasing the efficiency of pebble-accretion, once/if such large planetesimals and planetary embryos exist in the disc.

The formation of the first planetesimals and the final growth of planetary cores likely relies on the abundance of small pebbles. The efficiencies of both the streaming instability (SI) process, suggested to catalyze the early growth of planetesimals, and the pebble-accretion process, suggested to accelerate the growth of planetary cores, depend on the sizes of solids residing in the disk. In particular, these processes were found to be sensitive to size-distribution of solids, and efficient planetesimal formation and growth through these channels require a limited pebble size-distribution. Here we show that aeolian-erosion, a process that efficiently grinds down boulders into a mono-sized distribution of pebbles, provides a natural upper cut-off for the maximal pebble sizes. We find the dependence of the critical solids size cut-off on the radial separation, disk age, turbulence strength, and the grain-size composition of the boulders in the disk. This size cut-off shapes the size distribution of small pebbles and thereby helps catalyze the early onset of SI formation of planetesimals, and the later core accretion growth through pebble accretion.

We investigate the extent to which resonances between an oscillating background of ultra-light axion and a binary Keplerian system can affect the motion of the latter. These resonances lead to perturbations in the instantaneous time-of-arrivals, and to secular variations in the period of the binary. While the secular changes at exact resonance have recently been explored, the instantaneous effects have been overlooked. In this paper, we examine the latter using N-body simulations including the external oscillatory forcing induced by the axion background. While the secular effects are restricted to a narrow width near the resonance, the instantaneous changes, albeit strongest close to resonances, are apparent for wide range of configurations. We compute the signal-to-noise ratio (SNR) as a function of semi-major axis for a detection of axion oscillations through the R\{o} mer delay. The latter can be extracted from the time-of-arrivals of binary pulsars. The SNR broadly increases with increasing binary eccentricity in agreement with the secular expectation. However, we find that it differs significantly from the scaling a^{5/2} around the lowest orders of resonance. Future observations could probe these effects away from resonances and, therefore, constrain a much broader range of axion masses provided that binary pulsar systems are found near the central region of our Galaxy, and that the time-or-arrival measurement accuracy reaches < 10 ns

Scalar, tensor waves induce oscillatory perturbations in Keplerian systems which can be probed with measurements of pulsar timing residuals. In this paper, we consider the imprint of coherent oscillations produced by ultra-light axion dark matter on the Roemer time delay. We use the angle-action formalism to calculate the time evolution of the observed signal and its dependence on the orbital parameters and the axion phase. We emphasize the similarity of the expected signal-to-noise with the response of an harmonic oscillator to an external oscillatory driving. We validate our theoretical predictions with numerical simulations. Our results furnish a useful benchmark for numerical codes and analysis procedures and, hopefully, will motivate the search for such imprints in real data.

Following its flyby and first imaging the Pluto-Charon binary, the New Horizons spacecraft visited the Kuiper-Belt-Object (KBO) (486958) 2014 MU69 (Arrokoth). Imaging showed MU69 to be a contact-binary, made of two individual lobes connected by a narrow neck, rotating at low spin period (15.92 h), and having high obliquity (~98 deg), similar to other KBO contact-binaries inferred through photometric observations. The origin of such peculiar configurations is puzzling, and all scenarios suggested for the origins of contact-binaries fail to reproduce such properties and their likely high frequency. Here we show that semi-secular perturbations operating only on ultra-wide (~0.1-0.4 Hill-radius) KBO-binaries can robustly lead to gentle, slow-speed binary mergers at arbitrarily high obliquities, but low rotational velocities, that can reproduce MU69's (and similar oblique contact binaries) characteristics. Using N-body simulations, we find that ~15% of all ultra-wide binaries with cosine-uniform inclination distribution are likely to merge through this process. Moreover, we find that such mergers are sufficiently gentle as to only slightly deform the KBO shape, and can produce the measured rotation speed of MU69. The semi-secular contact-binary formation channel not only explains the observed properties of MU69, but could also apply for other Kuiper/asteroid belt binaries, and for Solar/extra-solar moon systems.

Context. GJ 1148 is an M-dwarf star hosting a planetary system composed of two Saturn-mass planets in eccentric orbits with periods of 41.38 and 532.02 days. Aims. We reanalyze the orbital configuration and dynamics of the GJ 1148 multi-planetary system based on new precise radial velocity (RV) measurements taken with CARMENES. Methods. We combined new and archival precise Doppler measurements from CARMENES with those available from HIRES for GJ 1148 and modeled these data with a self-consistent dynamical model. We studied the orbital dynamics of the system using the secular theory and direct N-body integrations. The prospects of potentially habitable moons around GJ 1148 b were examined. Results. The refined dynamical analyses show that the GJ 1148 system is long-term stable in a large phase-space of orbital parameters with an orbital configuration suggesting apsidal alignment, but not in any particular high-order mean-motion resonant commensurability. GJ 1148 b orbits inside the optimistic habitable zone (HZ). We find only a narrow stability region around the planet where exomoons can exist. However, in this stable region exomoons exhibit quick orbital decay due to tidal interaction with the planet. Conclusions. The GJ 1148 planetary system is a very rare M-dwarf planetary system consisting of a pair of gas giants, the inner of which resides in the HZ. We conclude that habitable exomoons around GJ 1148 b are very unlikely to exist.

Detection of gravitational-wave (GW) sources enables the characterisation of binary compact objects and of their in-spiral. However, other dissipative processes can affect the in-spiral. Here we show that the in-spiral of compact objects through a gaseous common-envelope (CE) arising from an evolved stellar companion produces a novel type of GW-sources, whose evolution is dominated by the dissipative gas dynamical friction effects from the CE, rather than the GW-emission itself. The evolution and properties of the GW-signals differ from those of isolated gas-poor mergers significantly. We find characteristic strains of $\sim10^{-23}$-$10^{-21}$ ($10{\rm kpc}/{D}$) for such sources -- observable by next-generation space-based GW-detectors. The evolution of the GW-signal can serve as a probe of the interior regions of the evolved star, and the final stages of CE-evolution, otherwise inaccessible through other observational means. Moreover, such CE-mergers are frequently followed by observable explosive electromagnetic counterparts and/or the formation of exotic stars.

Tertiary tides (TTs), or the continuous tidal distortion of the tertiary in a hierarchical triple system, can extract energy from the inner binary, inducing within it a proclivity to merge. Despite previous work on the subject, which established that it is significant for certain close triple systems, it is still not a well-understood process. A portion of our ignorance in this regard stems from our inability to integrate a simulation of this phenomenon into conventional stellar evolution codes, since full calculations of these tidal interactions are computationally expensive on stellar evolution timescales. Thus, to attain a better understanding of how these TTs act on longer timescales, an empirical expression of its effects as a function of parameters of the triple system involved is required. In our work, we evaluate the rate at which TTs extract energy from the inner binary within a series of constructed hierarchical triple systems under varying parameters, and study the rate at which the inner binary orbital separation shrinks as a function of those parameters. We find that this rate varies little with the absolute values of the masses of the three component objects, but is very sensitive to the mass ratio of the inner binary $q$, the tertiary radius $R_{\rm 3}$, the inner binary orbital separation $a_{\rm 1}$, the outer orbital separation $a_{\rm 2}$, and the viscoelastic relaxation time of the tertiary $\tau$. More specifically, we find that the percentage by which $a_{\rm 1}$ shrinks per unit time can be reasonably approximated by (1/$a_{\rm 1}$)(d$a_{\rm 1}$/d$t$)=$\left(2.22{\times}10^{-8}{\rm yrs}^{-1}\right)4q\left(1+q\right)^{-2}(R_{\rm 3}/100{\rm R}_{\odot})^{5.2}(a_{\rm 1}/0.2{\rm AU})^{4.8}(a_{\rm 2}/2{\rm AU})^{-10.2}$ $({\tau}/0.534{\rm yrs})^{-1.0}$. We also provide tests of how precise this fitting function is.

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.

Nuclear star clusters surrounding supermassive black holes (SMBHs) in galactic nuclei contain large numbers of stars, black holes (BHs) and neutron stars (NSs), a fraction of which are likely to form binaries. These binaries were suggested to form a triple system with the SMBH, which acts as a perturber and may enhance BH and NS mergers via the Lidov-Kozai mechanism. We follow-up previous studies, but for the first time perform an extensive statistical study of BH-BH, NS-NS and BH-NS binary mergers by means of direct high-precision regularized $N$-body simulations, including Post-Newtonian (PN) terms up to order PN2.5. We consider different SMBH masses, slopes for the BH mass function, binary semi-major axis and eccentricity distributions, and different spatial distributions for the binaries. We find that the merger rates are a decreasing function of the SMBH mass and are in the ranges $\sim 0.17$-$0.52 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$, $\sim 0.06$-$0.10 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$ and $\sim 0.04$-$0.16 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$ for BH-BH, BH-NS and NS-NS binaries, respectively. However, the rate estimate from this channel remains highly uncertain and depends on the specific assumptions regarding the star-formation history in galactic nuclei and the supply rate of compact objects. We find that $\sim 10\%$--$20\%$ of the mergers enter the LIGO band with eccentricities $\gtrsim 0.1$. We also compare our results to the secular approximation, and show that $N$-body simulations generally predict a larger number of mergers. Finally, these events can also be observable via their electromagnetic counterparts, thus making these compact object mergers especially valuable for cosmological and astrophysical purposes.

Planet formation begins with collisional growth of small planetesimals accumulating into larger ones. Such growth occurs while planetesimals are embedded in a gaseous protoplanetary disc. However, small-planetesimals experience collisions and gas-drag that lead to their destruction on short timescales, not allowing, or requiring fine tuned conditions for the efficient growth of $\sim$ metre-size objects. Here we show that $\sim10^{4}$ interstellar objects such as the recently detected 1I/2017-U1 ('Oumuamua) could have been captured, and become part of the young Solar System, together with up to hundreds of $\sim$ km sized ones. The capture rates are robust even for conservative assumptions on the protoplanetary disc structure, local stellar environment and planetesimal ISM density. 'Seeding' of such planetesimals then catalyze further planetary growth into planetary embryos, and potentially alleviate the main-challenges with the meter-size growth-"barrier". The capture model is in synergy with the current leading planet formation theories, providing the missing link to the first planetesimals. Moreover, planetesimal capture provides a far more efficient route for lithopanspermia than previously thought.

In hierarchical triple systems, the inner binary is perturbed by a distant companion. For large mutual inclinations, the Lidov-Kozai mechanism secularly excites large eccentricity and inclination oscillations of the inner binary. The maximal eccentricity attained, $e_{\rm max}$ is simply derived and widely used. However, for mildly hierarchical systems (i.e. the companion is relatively close and massive), non-secular perturbations affect the evolution. Here we account for fast non-secular variations and find new analytic formula for $e_{{\rm max}}$, in terms of the system's hierarchy level, correcting previous work and reproducing the orbital flip criteria. We find that $e_{{\rm max}}$ is generally enhanced, allowing closer encounters between the inner binary components, thus significantly changing their interaction and its final outcome. We then extend our approach to include additional relativistic and tidal forces. Using our results, we show that the merger time of gravitational-wave (GW) sources orbiting massive black-holes in galactic nuclei is enhanced compared with previous analysis accounting only for the secular regime. Consequently, this affects the distribution and rates of such GW sources in the relevant mild-hierarchy regime. We test and confirm our predictions with direct N-body and 2.5-level Post-Newtonian codes. Finally, we calculate the formation and disruption rates of hot-Jupiters (HJ) in planetary systems using a statistical approach, which incorporates our novel results for $e_{{\rm max}}$. We find that more HJ migrate from further out, but they are also tidally disrupted more frequently. Remarkably, the overall formation rate of HJs remains similar to that found in previous studies. Nevertheless, the different rates could manifest in different underlying distribution of observed warm-Jupiters.

We study the chaotic and secular evolution of hierarchical quadruple systems in the $3+1$ configuration, focusing on the evolution of mutual inclination of the inner binaries as the system undergoes coupled Lidov-Kozai (LK) oscillations. We include short-range forces (SRF; such as those due to tidal and rotational distortions) that control the eccentricity excitation of the inner binary. The evolution of mutual inclination is described, a priori, by two dimensionless parameters, $\pazocal{R}_0$, the ratio between the inner and outer LK time-scales and $\epsilon_{SRF}$, the ratio between the SRF precession and the inner LK precession rates. We find that the chaotic zones for the mutual inclination depend mainly on $\pazocal{R}_0$, while $\epsilon_{SRF}$ controls mainly the range of eccentricity excitation. The mutual inclination evolves chaotically for $1\lesssim \pazocal{R}_0\lesssim 10$, leading to large misalignments. For $0.4 \lesssim \pazocal{R}_0 \lesssim 0.8$, the system could be weakly excited and produce bimodal distribution of mutual inclination angles. Our results can be applied to exomoons-planets in stellar binaries and Warm/Hot Jupiters in stellar triples. Such systems could develop large mutual inclination angles if the inner binary is tight enough, and also high eccentricities, depending of the strength of the short-range forces. Future detections of tilted Warm/Hot Jupiters and exomoons could put our mechanism under observational tests.

A fundamental aspect of the three-body problem is its stability. Most stability studies have focused on the co-planar three-body problem, deriving analytic criteria for the dynamical stability of such pro/retrograde systems. Numerical studies of inclined systems phenomenologically mapped their stability regions, but neither complement it by theoretical framework, nor provided satisfactory fit for their dependence on mutual inclinations. Here we present a novel approach to study the stability of hierarchical three-body systems at arbitrary inclinations, which accounts not only for the instantaneous stability of such systems, but also for the secular stability and evolution through Lidov-Kozai cycles and evection. We generalize the Hill-stability criteria to arbitrarily inclined triple systems, explain the existence of quasi-stable regimes and characterize the inclination dependence of their stability. We complement the analytic treatment with an extensive numerical study, to test our analytic results. We find excellent correspondence up to high inclinations $(\sim120^{\circ}$), beyond which the agreement is marginal. At such high inclinations the stability radius is larger, the ratio between the outer and inner periods becomes comparable, and our secular averaging approach is no longer strictly valid. We therefore combine our analytic results with polynomial fits to the numerical results to obtain a generalized stability formula for triple systems at arbitrary inclinations. Besides providing a generalized secular-based physical explanation for the stability of non co-planar systems, our results have direct implications for any triple systems, and in particular binary planets and moon/satellite systems; we briefly discuss the latter as a test case for our models.

One of first the stages of planet formation is the growth of small planetesimals and their accumulation into large planetesimals and planetary embryos. This early stage occurs much before the dispersal of most of the gas from the protoplanetary disk. At this stage gas-planetesimal interactions play a key role in the dynamical evolution of \emph{single} intermediate-mass planetesimals ($m_{p}\sim10^{21}-10^{25}g$) \emph{through gas dynamical friction} (GDF). A significant fraction of all Solar system planetesimals (asteroids and Kuiper-belt objects) are known to be binary planetesimals (BPs). Here, we explore the effects of GDF on the evolution of \emph{binary} planetesimals embedded in a gaseous disk using an N-body code with a fiducial external force accounting for GDF. We find that GDF can induce binary mergers on timescales shorter than the disk lifetime for masses above $m_{p}\gtrsim10^{22}g$ at 1AU, independent of the binary initial separation and eccentricity. Such mergers can affect the structure of merger-formed planetesimals, and the GDF-induced binary inspiral can play a role in the evolution of the planetesimal disk. In addition, binaries on eccentric orbits around the star may evolve in the supersonic regime, where the torque reverses and the binary expands, which would enhance the cross section for planetesimal encounters with the binary. Highly inclined binaries with small mass ratios, evolve due to the combined effects of Kozai-Lidov cycles with GDF which lead to chaotic evolution. Prograde binaries go through semi-regular Kozai-Lidov evolution, while retrograde binaries frequently flip their inclination and $\sim50\%$ of them are destroyed.

The growth of small planetesimals into large planetary embryos occurs much before the dispersal of the gas from the protoplanetary disk. The planetesimal - gaseous-disk interactions give rise to migration and orbital evolution of the planetesimals/planets. Small planetesimals are dominated by aerodynamic gas drag. Large protoplanets, $m\sim0.1M_{\oplus}$, are dominated by type I migration \emph{differential} torque. There is an additional mass range, $m\sim10^{21}-10^{25}g$ of \emph{intermediate mass} planetesimals (IMPs), where gravitational interactions with the disk dominate over aerodynamic gas drag, but for which such interactions were typically neglected. Here we model these interactions using the \emph{gas dynamical friction} (GDF) approach, previously used to study the disk-planet interactions at the planetary mass range. We find the critical size where GDF dominates over gas drag, and then study the implications of GDF on single IMPs. We find that planetesimals with small inclinations rapidly become co-planar. Eccentric orbits circularize within a few Myrs, provided the the planetesimal mass is large, $m\gtrsim10^{23}g$ and that the initial eccentricity is low, $e\lesssim0.1$. Planetesimals of higher masses, $m\sim10^{24}-10^{25}g$ inspiral on a time-scale of a few Myrs, leading to \emph{an embryonic migration} to the inner disk. This can lead to an over-abundance of rocky material (in the form of IMPs) in the inner protoplanetary disk ($<1$AU) and induce rapid planetary growth. This can explain the origin of super-Earth planets. In addition, GDF damps the velocities of IMPs, thereby cooling the planetesimal disk and affecting its collisional evolution through quenching the effects of viscous stirring by the large bodies.

The dwarf planet Pluto is known to host an extended system of five co-planar satellites. Previous studies have explored the formation and evolution of the system in isolation, neglecting perturbative effects by the Sun. Here we show that secular evolution due to the Sun can strongly affect the evolution of outer satellites and rings in the system, if such exist. Although precession due to extended gravitational potential from the inner Pluto-Charon binary quench such secular evolution up to $a_{crit}\sim0.0035$ AU ($\sim0.09$ $R_{Hill}$ the Hill radius; including all of the currently known satellites), outer orbits can be significantly altered. In particular, we find that \emph{co-planar} rings and satellites should not exist beyond $a_{crit}$; rather, satellites and dust particles in these regions secularly evolve on timescales ranging between $10^{4}-10^{6}$ yrs, and quasi-periodically change their inclinations and eccentricities through secular evolution (Lidov-Kozai oscillations). Such oscillations can lead to high inclinations and eccentricities, constraining the range where such satellites (and dust particles) can exist without crossing the orbits of the inner satellites, or crossing the outer Hill stability range. Outer satellites, if such exist are therefore likely to be \emph{irregular} satellites, with orbits limited to be non-circular and/or highly inclined. These could be potentially detected and probed by the New-Horizon mission, possibly providing direct evidence for the secular evolution of the Pluto satellite system, and shedding new light on its origins.