Search

Massive stars may form in or be captured into AGN disks. Recent 1D studies employing stellar-evolution codes have demonstrated the potential for rapid growth of such stars through accretion up to a few hundred $M_\odot$. We perform 3D radiation hydrodynamic simulations of moderately massive stars' envelopes, in order to determine the rate and critical radius $R_{\rm crit}$ of their accretion process in an isotropic gas-rich environment in the absence of luminosity-driven mass loss. We find that in the ``fast-diffusion" regime where characteristic radiative diffusion speed $c/\tau$ is faster than the gas sound speed $c_s$, the accretion rate is suppressed by feedback from gravitational and radiative advection energy flux, in addition to the stellar luminosity. Alternatively, in the ``slow-diffusion" regime where $c/\tau

Star-forming regions, characterized by dense environments, experience frequent encounters that significantly influence binary systems, leading to their hardening, softening, or ionization. We extend the Hut \& Bahcall formalism to derive an analytical expression for the ionization cross-section in extreme mass ratio binary systems, allowing us to investigate the orbital evolution and population dynamics of binary planets and binary brown dwarfs in star clusters, while considering ongoing binary system formation. Our findings reveal that for low-mass soft binaries, the semi-major axis distribution asymptotes to a universal power law between $\propto a^{-8/3}$ and $\propto a^{-5/3}$ over the derived ionization timescale. We also discuss the implications of our results for the candidate Jupiter-mass binary objects putatively reported in the Trapezium cluster. We demonstrate that if their existence is verified, they likely form continuously with a spectrum proportional to $a^{1}$, aligning better with the ejection mechanism than with the in-situ formation mechanism, which predicts a distribution roughly proportional to $a^{-1}$. However, this implies an impractically high ejection formation rate. Alternatively, if these objects are binary brown dwarfs, continuous in-situ formation ($\propto a^{-1}$) with an initial minimal semi-major axis around 20 AU and a formation rate of 100 Myr$^{-1}$ plausibly matches the observed number of single objects, binary number, binary fraction, and semi-major axis distribution., Comment: Comments are welcome

Migration commonly occurs during the epoch of planet formation. For emerging gas giant planets, it proceeds concurrently with their growth through the accretion of gas from their natal protoplanetary disks. Similar migration process should also be applied to the stellar-mass black holes embedded in active galactic nucleus disks. In this work, we perform high resolution 3D and 2D numerical hydrodynamical simulations to study the migration dynamics for accreting embedded objects over the disk viscous timescales in a self-consistent manner. We find that an accreting planet embedded in a predominantly viscous disk has a tendency to migrate outward, in contrast to the inward orbital decay of non-accreting planets. 3D and 2D simulations find the consistent outward migration results for the accreting planets. Under this circumstance, the accreting planet's outward migration is mainly due to the asymmetric spiral arms feeding from the global disk into the Hill radius. This is analogous to the unsaturated corotation torque although the imbalance is due to material accretion within the libration timescale rather than diffusion onto the inner disk. In a disk with a relatively small viscosity, the accreting planets clear deep gaps near their orbits. The tendency of inward migration is recovered, albeit with suppressed rates. By performing a parameter survey with a range of disks' viscosity, we find that the transition from outward to inward migration occurs with the effective viscous efficiency factor $\alpha\sim 0.003$ for Jupiter-mass planets., Comment: 17 pages, 9 figures, Accepted by ApJ

Context. About 25% -- 50% of white dwarfs are found to be contaminated by heavy elements, which are believed to originate from external sources such as planetary materials. Elemental abundances suggest that most of the pollutants are rocky objects and only a small fraction of white dwarfs bear traces of volatile accretion. Aims. In order to account for the scarcity of volatile pollution, we investigate the role of the white dwarfs' magnetospheres in shielding the volatile content of icy objects. Methods. We estimated the volatile sublimation of inward-drifting exocomets. We assume the orbits of the exocomets are circularized by the Alfven wing drag that is effective for long-period comets. Results. Volatile material can sublimate outside the corotation radius and be shielded by the magnetic field. {The two conditions for this volatile-shielded mechanism are that the magnetosphere radius must be larger than the corotation radius and that the volatiles are depleted outside the corotation radius, which requires a sufficiently slow orbital circularization process.} We applied our model to nine white dwarfs with known rotational periods, magnetic fields, and atmosphere compositions. Our volatile-shielded model may explain the excess of volatile elements such as C and S in the disk relative to the white dwarf atmosphere in WD2326+049 (G29-38). Nevertheless, given the sensitivity of our model to the circularization process and material properties of icy objects, there remains considerable uncertainty in our results. Conclusions. Our work suggests a possible explanation for the scarcity of volatile-accretion signatures among white dwarfs. We also identify a correlation between the magnetic field strength, the spin period, and the composition of pollutants in white dwarf atmospheres., Comment: 13 pages, 11 figures. Accepted for publication in A&A

Gravitational instability in the outskirts of Active Galactic Nuclei (AGN) disks lead to disk fragmentation and formation of super-massive (several 10^2Msun) stars with potentially long lifetimes. Alternatively, stars can be captured ex-situ and grow from gas accretion in the AGN disk. However, the number density distribution throughout the disk is limited by thermal feedback as their luminosities provide the dominant heating source. We derive equilibrium stellar surface density profiles under two limiting contexts: in the case where the stellar lifetimes are prolonged due to recycling of hydrogen rich disk gas, only the fraction of gas converted into heat is removed from the disk accretion flow. Alternatively, if stellar composition recycling is inefficient and stars can evolve off the main sequence, the disk accretion rate is quenched towards smaller radii resembling a classical star-burst disk, albeit the effective removal rate depends not only on the stellar lifetime, but also the mass of stellar remnants. For AGNs with central Supermassive Black Hole (SMBH) masses of \sim 10^6 to 10^8Msun accreting at \sim 0.1 Eddington efficiency, we estimate a total number of 10^3 to 10^5 coexisting massive stars and the rate of stellar mergers to be 10^-3 to 1 per year. We motivate the detailed study of interaction between a swarm of massive stars through hydro and N body simulations to provide better prescriptions of dynamical processes in AGN disks, and to constrain more accurate estimates of the stellar population., Comment: Published in ApJ

In the classical grouping of large magnitude episodic variability of young accreting stars, FUORs outshine their stars by a factor of $\sim$ 100, and can last for up to centuries; EXORs are dimmer, and last months to a year. A disc Hydrogen ionisation Thermal Instability (TI) scenario was previously proposed for FUORs but required unrealistically low disc viscosity. In the last decade, many intermediate type objects, e.g., FUOR-like in luminosity and spectra but EXOR-like in duration were found. Here we show that the intermediate type bursters Gaia20eae, PTF14jg, Gaia19bey and Gaia21bty may be naturally explained by the TI scenario with realistic viscosity values. We argue that TI predicts a dearth (desert) of bursts with peak accretion rates between $\dot M \sim 10^{-6} M_\odot$/yr and $\dot M \sim 10^{-5} M_\odot$/yr, and that this desert is seen in the sample of all the bursters with previously determined $\dot M$ burst. Most classic EXORs (FUORs) appear to be on the cold (hot) branch of the S-curve during the peak light of their eruptions; thus TI may play a role in this class differentiation. At the same time, TI is unable to explain how classic FUORs can last for up to centuries, and over-predicts the occurrence rate of short FUORs by at least an order of magnitude. We conclude that TI is a required ingredient of episodic accretion operating at R < 0.1 au, but additional physics must play a role at larger scales. Knowledge of TI inner workings from related disciplines may enable its use as a tool to constrain the nature of this additional physics., Comment: accepted to MNRAS, 18 pages

Dust trapping in the global pressure bump induced by magnetospheric truncation offers a promising formation mechanism for close-in super-Earths/sub-Neptunes. These planets likely form in evolved protoplanetary discs, where the gas temperature at the expanding truncation radius become amiable to refractory solids. However, dust accumulation may alter the disc opacity such that thermal evolution is inevitable. To better understand how thermodynamics affects this planet formation pathway, we conduct a suite of local dust evolution simulations in an idealized inner disc model. Our calculations take into account self-consistent opacity-dependent temperature changes as well as dust evaporation and vapour condensation. We find that disc thermal evolution regulates dust growth and evolution, discouraging any accumulation of small particles that drives the increase of opacity and temperature. Significant retention of dust mass takes place when the disc environments allow runaway growth of large solids beyond the fragmentation barrier, where small particles are then swept up and preserved. Our results further validate dust accumulation near disc truncation as a promising mechanism to form close-in planets., Comment: Accepted by MNRAS and in production; 11 pages, 8 figures

We perform a series of 3D simulations to study the accretion of giant planet embedded in protoplanetary discs (PPDs) over gap-opening timescales. We find that the accretion mass flux mainly comes from the intermediate latitude above the disc midplane. The circumplanetary disc (CPD) for a super-thermal planet is rotation-supported up to $\sim$20-30\% of the planet Hill radius. While both mass inflow and outflow exists in the CPD midplane, the overall trend is an outflow that forms a meridional circulation with high-latitude inflows. We confirm the absence of accretion outburst from disc eccentricity excited by massive planets in our 3D simulations, contrary to the consensus of previous 2D simulations. This suggests the necessity of 3D simulations of accretion even for super-Jupiters. The accretion rates of planets measured in steady-state can be decomposed into the ``geometric" and ``density depletion" factors. Through extensive parameter survey, we identify a power-law scaling for the geometric factor $\propto q_{\rm th}^{2/3}$ for super-thermal planets ($q_{\rm th}$ being the thermal mass ratio), which transforms to $\propto q_{\rm th}^{2}$ for less massive cases. The density depletion factor is limited by the disc accretion rate for mildly super-thermal planets, and by gap-opening for highly super-thermal ones. Moderate planetary eccentricities can enhance the accretion rates by a factor of $2-3$ through making the gap shallower, but does not impact the flow geometry. We have applied our simulations results to accreting protoplanet system PDS 70 and can satisfactorily explain the accretion rate and CPD size in observations., Comment: 19 pages, 17 figures, Accepted for publication in MNRAS

Recent observations of changing-look active galactic nuclei (AGN) hint at a frequency of accretion activity not fully explained by tidal disruption events (TDEs) stemming from relaxation processes in nucleus star clusters (NSCs), traditionally estimated to occur at rates of $10^{-4}$ to $10^{-5}$ yr$^{-1}$ per galaxy. In this letter, we propose an enhanced TDE rate through the AGN disk capture process, presenting a viable explanation for the frequent transitions observed in changing-look AGN. Specifically, we investigate the interaction between the accretion disk and retrograde stars within NSCs, resulting in the rapid occurrence of TDEs within a condensed time frame. Through detailed calculations, we derive the time-dependent TDE rates for both relaxation-induced TDE and disk-captured TDE. Our analysis reveals that TDEs triggered by the disk capture process can notably amplify the TDE rate by several orders of magnitude during the AGN phase. This mechanism offers a potential explanation for the enhanced high-energy variability characteristic of changing-look AGNs., Comment: Submitted to ApJL. Comments are welcome

Stars are likely to form or to be captured in AGN disks. Their mass reaches an equilibrium when their rate of accretion is balanced by that of wind. If the exchanged gas is well mixed with the stellar core, this metabolic process would indefinitely sustain an "immortal" state on the main sequence (MS) and pollute the disk with He byproducts. This theoretical extrapolation is inconsistent with the super-solar {\alpha} element and Fe abundances inferred from the broad emission lines in active AGNs with modest He concentration. We show this paradox can be resolved with a highly-efficient retention of the He ashes or the suppression of chemical blending. The latter mechanism is robust in the geometrically-thin, dense, sub-pc regions of the disk where the embedded-stars' mass is limited by the gap-formation condition. These stars contain a radiative zone between their mass-exchange stellar surface and the nuclear-burning core. Insulation of the core lead to the gradual decrease of its H fuel and the stars' equilibrium masses. These stars transition to their post-main-sequence (PostMS) tracks on a chemical evolution time scale of a few Myr. Subsequently, the triple-{\alpha} and {\alpha}-chain reactions generate {\alpha} and Fe byproducts which are released into their natal disks. These PostMS stars also undergo core collapse, set off type II supernova, and leave behind a few solar-mass residual black holes or neutron stars, Comment: 17 pages, 7 figures, Accepted for publication in MNRAS

We separately assess elemental abundances in AGNs' broad and narrow emission line regions (BLR and NLR), based on a critical assessment of published results together with new photoionization models. We find 1) He/H enhancements in some AGN, exceeding what can be explained by normal chemical evolution and confirm 2) super-solar {\alpha} abundance, though to a lesser degree than previously reported. We also reaffirm 3) a N/O ratio consistent with secondary production; 4) solar or slightly sub-solar Fe abundance; and 5) red-shift independent metallicity, in contrast with galactic chemical evolution. We interpret 6) the larger metallicity in the BLR than NRL in terms of an in situ stellar evolution and pollution in AGN discs (SEPAD) model. We attribute: a) the redshift independence to the heavy element pollutants being disposed into the disc and accreted onto the central supermassive black hole (SMBH); b) the limited He excess to the accretion-wind metabolism of a top-heavy population of evolving massive main sequence stars; c) the super-solar CNO enrichment to the nuclear synthesis during their post-main-sequence evolution; d) the large N/O to the byproduct of multiple stellar generations; and e) the Mg, Si, and Fe to the ejecta of type II supernovae in the disc. These results provide supporting evidence for f) ongoing self-regulated star formation, g) adequate stellar luminosity to maintain marginal gravitational stability, h) prolific production of seeds and i) dense coexistence of subsequently-grown residual black hole populations in AGN discs., Comment: 17 pages, 8 figures, accepted for publication in MNRAS

The Active Galatic Nuclei(AGN) disk has been proposed as a potential channel for the merger of binary black holes. The population of massive stars and black holes in AGN disks captured from the nuclei cluster plays a crucial role in determining the efficiency of binary formation and final merger rate within the AGN disks. In this paper, we investigate the capture process using analytical and numerical approaches. We discover a new constant integral of motion for one object's capture process. Applying this result to the whole population of the nuclei cluster captured by the AGN disk, we find that the population of captured objects depends on the angular density and eccentricity distribution of the nuclei clusters and is effectively independent of the radial density profile of the nuclei cluster and disk models. An isotropic nuclei cluster with thermal eccentricity distribution predicts a captured profile $d N/d r \propto r^{-1/4}$. The captured objects are found to be dynamically crowded within the disk. Direct binary formation right after the capture would be promising, especially for stars. The conventional migration traps that help pile up single objects in AGN disks for black hole mergers might not be required., Comment: Typo corrections

The orbital distribution of the S-star cluster surrounding the supermassive black hole in the center of the Milky Way is analyzed. A tight, roughly exponential dependence of the pericenter distance r$_{p}$ on orbital eccentricity e$_{\star}$ is found, $\log ($r$_p)\sim$(1-e$_{\star}$), which cannot be explained simply by a random distribution of semi-major axes and eccentricities. No stars are found in the region with high e$_{\star}$ and large log r$_{p}$ or in the region with low e$_{\star}$ and small log r$_{p}$. G-clouds follow the same correlation. The likelihood P(log r$_p$,(1-e$_{\star}$)) to determine the orbital parameters of S-stars is determined. P is very small for stars with large e$_{\star}$ and large log r$_{p}$. S-stars might exist in this region. To determine their orbital parameters, one however needs observations over a longer time period. On the other hand, if stars would exist in the region of low log r$_{p}$ and small e$_{\star}$, their orbital parameters should by now have been determined. That this region is unpopulated therefore indicates that no S-stars exist with these orbital characteristics, providing constraints for their formation. We call this region, defined by $\log$ (r$_p$/AU) $<$ 1.57+2.6(1-e$_{\star})$, the zone of avoidance. Finally, it is shown that the observed frequency of eccentricities and pericenter distances is consistent with a random sampling of log r$_{p}$ and e$_{\star}$. However, only if one takes into account that no stars exist in the zone of avoidance and that orbital parameters cannot yet be determined for stars with large r$_{p}$ and large e$_{\star}$., Comment: 15 pages, 4 figures, submitted to ApJ, comments very welcome

Accretion discs around super-massive black holes (SMBH) not only power active galactic nuclei (AGNs), but also host single and binary embedded stellar-mass black holes (EBHs) that grow rapidly from gas accretion. The merger of these EBHs provides a promising mechanism for the excitation of some gravitational wave events observed by LIGO-Virgo, especially those with source masses considerably larger than isolated stellar-mass black hole binaries. In addition to their mass and mass-ratio distribution, their hitherto enigmatic small spin-parameters chi_effective carry important clues and stringent constraints on their formation channels and evolutionary pathways. Here we show that, between each coalescence, the typical rapid spin of the merged EBHs is suppressed by their subsequent accretion of gas from a turbulent environment, due to its ability to randomize the flow's spin orientation with respect to that of the EBHs on an eddy-turnover timescale. This theory provides supporting evidence for the prolificacy of EBH mergers and suggests that their mass growth may be dominated by gas accretion rather than their coalescence in AGN discs., Comment: accepted by MNRAS, 11 pages, 8 figures

Hydrodynamical interaction in circumbinary discs (CBDs) plays a crucial role in various astrophysical systems, ranging from young stellar binaries to supermassive black hole binaries in galactic centers. Most previous simulations of binary-disc systems have adopted locally isothermal equation of state. In this study, we use the grid-based code $\texttt{Athena++}$ to conduct a suite of two-dimensional viscous hydrodynamical simulations of circumbinary accretion on a cartesian grid, resolving the central cavity of the binary. The gas thermodynamics is treated by thermal relaxation towards an equilibrium temperature (based on the constant$-\beta$ cooling ansatz, where $\beta$ is the cooling time in units of the local Keplerian time). Focusing on equal mass, circular binaries in CBDs with (equilibrium) disc aspect ratio $H/R=0.1$, we find that the cooling of the disc gas significantly influences the binary orbital evolution, accretion variability, and CBD morphology, and the effect depends sensitively on the disc viscosity prescriptions. When adopting a constant kinematic viscosity, a finite cooling time ($\beta \gtrsim 0.1$) leads to binary inspiral as opposed to outspiral and the CBD cavity becomes more symmetric. When adopting a dynamically varying $\alpha-$viscosity, binary inspiral only occurs within a narrow range of cooling time (corresponding to $\beta$ around 0.5)., Comment: 20 pages, 20 figures, Accepted for publication in MNRAS

We propose a solution to the problem of accretion disc sizes in active galactic nuclei being larger when measured by reverberation mapping than predicted by theory. Considering blackbody reprocessing on a disc with thickness profile $H(r)$, our solution invokes a steep rim or rippled structures irradiated by the central lamp-post. We model the continuum lags and the faint and bright disc spectral energy distribution (SED) in the best-studied case NGC 5548 (black hole mass $M = 7\times10^{7} M_\odot$, disc inclination $i=45^\circ$). With the lamp-post off, the observed disc SED requires a low accretion rate ($\dot{M} \sim 0.0014 M_\odot$/yr) and high prograde black hole spin ($a \sim 0.93$). Reprocessing on the thin disc gives time lags increasing with wavelength but 3 times smaller than observed. Introducing a steep $H(r)$ rim, or multiple crests, near $r = 5$ light days, reprocessing on their steep centre-facing slopes increases temperatures from $\sim1500$ K to $\sim6000$ K and this increases optical lags to match the lag data. Most of the disc surface maintains the cooler $T\propto r^{-3/4}$ temperature profile that matches the SED. The bright lamp-post may be powered by magnetic links tapping the black hole spin. The steep rim occurs near the sublimation radius for dust in the disc, as in the "failed disc wind model" for broad-line clouds. Lens-Thirring torques aligning the disc and black hole spin may also raise a warp and associated waves. In both scenarios, the small density scale height implied by the inferred value of $H(r)$ suggests possible marginal gravitational instability in the disc., Comment: 17 pages, 7 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

The Galactic Centre is known to have undergone a recent star formation episode a few Myrs ago, which likely produced many T Tauri stars hosting circumstellar discs. It has been suggested that these discs may be the compact and dusty ionized sources identified as ``G-clouds''. Given the Galactic Centre's hostile environment, we study the possible evolutionary pathways these discs experience. We compute new external photoevaporation models applicable to discs in the Galactic Centre that account for the sub-sonic launching of the wind and absorption of UV photons by dust. Using evolutionary disc calculations, we find that photoevaporation's rapid truncation of the disc causes them to accrete onto the central star rapidly. Ultimately, an accreting circumstellar disc has a lifetime $\lesssim1~$Myr, which would fail to live long enough to explain the G-clouds. However, we identify a new evolutionary pathway for circumstellar discs in the Galactic Centre. Removal of disc material by photoevaporation prevents the young star from spinning down due to magnetic braking, ultimately causing the rapidly spinning young star to torque the disc into a ``decretion disc'' state which prevents accretion. At the same time, any planetary companion in the disc will trap dust outside its orbit, shutting down photoevaporation. The disc can survive for up to $\sim$10 Myr in this state. Encounters with other stars are likely to remove the planet on Myr timescales, causing photoevaporation to restart, giving rise to a G-cloud signature. A giant planet fraction of $\sim10\%$ can explain the number of observed G-clouds., Comment: Accepted for publication in MNRAS

We study the long-term radial transport of micron to mm-size grain in protostellar disks (PSDs) based on diffusion and viscosity coefficients measured from 3D global stratified-disk simulations with a Lagrangian hydrodynamic method. While gas-drag tend to transport dust species radially inwards, stochastic diffusion can spread a considerable fraction of dust radially outwards (upstream) depending on the nature of turbulence. In gravitationally unstable disks, we measure a high radial diffusion coefficient Dr with little dependence on altitude. This leads to strong and vertically homogeneous upstream diffusion in early PSDs. In the solar nebula, the robust upstream diffusion of micron to mm size grains not only efficiently transports highly refractory mocron-size grains (such as those identified in the samples of comet 81P/Wild 2) from their regions of formation inside the snow line out to the Kuiper Belt, but can also spread mm-size CAI formed in the stellar proximity to distances where they can be assimilated into chondritic meteorites. In disks dominated by magnetorotational instability (MRI), the upstream diffusion effect is generally milder, with a separating feature due to diffusion being stronger in the surface layer than the midplane. This variation becomes much more pronounced if we additionally consider a quiescent midplane with lower turbulence and larger characteristic dust size due to non-ideal MHD effects. This segregation scenario helps to account for dichotomy of two dust populations' spatial distribution as observed in scattered light and ALMA images., Comment: ApJ published, 14 pages, 12 figures

We scrutinize the light sail scenario of the first interstellar object (ISO) 1I/2017 U1 (`Oumuamua) by making comparisons between physical models and observational constraints. These analyses can be generalized for future surveys of `Oumuamua-like objects. The light sail goes through a drift in the interstellar space due to the magnetic field and gas atoms, which poses challenges to the navigation system. When the light sail enters the inner solar system, the sideways radiation pressure leads to a considerable non-radial displacement. The immensely high dimensional ratio and the tumbling motion could cause a light curve with an extremely large amplitude and could even make the light sail invisible from time to time. These observational features allow us to examine the light sail scenario of interstellar objects. Our results show that the drift of the freely rotating light sail in the interstellar medium is $\sim 100\,$au even if the travel distance is only 1 pc. The probability is < 1.5\% for the expected brightness modulation of the light sail to match with `Oumuamua's observed variation amplitude ($\sim$ 2.5 -- 3). In addition, the probability is 0.4% for the tumbling light-sail to be visible (brighter than V=27) in all 55 observations spread over two months after discovery. Radiation pressure could cause a larger displacement that is normal to the orbital plane for a lightsail than that for `Oumuamua. Also, the ratio of anti-solar to sideways acceleration of `Oumuamua deviates from that of the light sail by ~ 1.5{\sigma}. We suggest that `Oumuamua is unlikely a light sail. The dynamics of an intruding light sail, if exist, has distinct observational signatures, which can be quantitatively identified and analyzed with our methods in future surveys., Comment: 17 pages, 13 figures Comments are welcome

Mechanisms have been proposed to enhance the merger rate of stellar mass black hole binaries, such as the Von Zeipel-Lidov-Kozai mechanism (vZLK). However, high inclinations are required in order to greatly excite the eccentricity and to reduce the merger time through vZLK. Here, we propose a novel pathway through which compact binaries could merge due to eccentricity increase in general, including in a near coplanar configuration. Specifically, a compact binary migrating in an AGN disk could be captured in an evection resonance, when the precession rate of the binary equals their orbital period around the supermassive black hole. In our study we include precession to due first-order post Newtonian precession as well as that due to disk around one or both components of the binary. Eccentricity is excited when the binary sweeps through the resonance which happens only when it migrates on a timescale 10-100 times the libration timescale of the resonance. Libration timescale decreases as the mass of the disk increases. The eccentricity excitation of the binary can reduce the merger timescale by a factor up to $\sim 10^{3-5}$., Comment: Accepted for publication in ApJ

Spin evolution of stellar-mass Black Holes (sBHs) embedded in AGN accretion disks is an important process relevant to production of gravitaional waves from binary Black Hole (BBH) merger events through the AGN channel. Since embedded sBHs are surrounded by circum-stellar disks (CSDs), the rotation of CSD gas flows determine the direction of the angular momentum it accretes. In this Letter, we use global 2D hydrodynamic simulations to show that while a disk-embedded sBH on a circular orbit transforms the initial retrograde Keplerian shear of the background accretion disk into a prograde CSD flow, as in the classical picture of companion-disk interaction theory, moderate orbital eccentricity could disrupt the steady-state tidal perturbation and preserve a retrograde CSD flow around the sBH. This switch of CSD orientation occurs at a transition eccentricity that scales nearly proportional with local sound speed. This bifurcation in the CSD flow and thereafter spin-up direction of SBHs leads to formation of a population of nearly anti-aligned sBHs and should be incorporated in future population models of sBH and BBH evolutions., Comment: ApJL published

The spin-orbit misalignment of stellar-mass black hole (sBH) binaries provide important constraints on the formation channels of merging sBHs. Here, we study the role of secular spin-orbit resonance in the evolution of a sBH binary component around a supermassive BH (SMBH) in an AGN disk. We consider the sBH's spin-precession due to the $J_2$ moment introduced by a circum-sBH disk within the warping/breaking radius of the disk. We find that the sBH's spin-orbit misalignment (obliquity) can be excited via spin-orbit resonance between the sBH binary's orbital nodal precession and the sBH spin-precession driven by a massive circum-sBH disk. Using an $\alpha$-disk model with Bondi-Hoyle-Lyttleton accretion, the resonances typically occur for sBH binaries with semi-major axis of $1$AU, and at a distance of $\sim 1000$AU around a $10^7$\msun SMBH. The spin-orbit resonances can lead to high sBH obliquities, and a broad distribution of sBH binary spin-spin misalignments. However, we note that the Bondi-Hoyle-Lyttleton accretion is much higher than that of Eddington accretion, which typically results in spin precession being too low to trigger spin-orbit resonances. Thus, the secular spin-orbit resonances can be quite rare for sBHs in AGN disks., Comment: 11 pages, 13 figures. Accepted for publication in ApJ

The prevalence of short-period super-Earths that are independent of host metallicity challenges the theoretical construction of their origin. We propose that dust trapping in the global pressure bump induced by magnetospheric truncation in evolved protoplanetary discs (PPDs) around T Tauri stars offers a promising formation mechanism for super-Earths, where the host metallicity is already established. To better understand this planet forming scenario, we construct a toy inner disc model and focus on the evolution of dust trapped in the bump, taking into account the supply from drifting pebbles and loss due to funnel flows. We develop an implicit coagulation-fragmentation code, $\mathtt{Rubble}$, and perform a suite of simulations to evolve the local dust size distributions. Our study for the first time considers dust feedback effect on turbulent diffusion in this kind of model. We report that efficient dust growth and significant accumulation of dust mass is possible in less turbulent disc with sturdier solids and with faster external supply, laying out a solid foundation for further growth towards planetesimals and planetary embryos. We further find that, depending on the dominant process, solid mass may predominantly accumulate in cm-sized grains or particles in runaway growth, indicating different ways of forming planetesimals. Furthermore, these various outcomes show different efficiencies in saving dust from funnel flows, suggesting that they may be distinguishable by constraining the opacity of funnel flows. Also, these diverse dust behaviours may help explain the observed dipper stars and rapidly varying shadows in PPDs., Comment: Accepted by MNRAS and in production; 20 pages, 14 figures

We investigate the origins of cold sub-Saturns (CSS), an exoplanetary population inferred from microlensing surveys. If confirmed, these planets would rebut a theorised gap in planets' mass distribution between those of Neptune and Jupiter caused by the rapid runaway accretion of super-critical cores. In an attempt to resolve this theoretical-observational disparity, we examine the outcomes of giant collisions between sub-critical protoplanets. Due to the secular interaction among protoplanets, these events may occur in rapidly depleting discs. We show that impactors ~ 5% the mass of near-runaway envelopes around massive cores can efficiently remove these envelopes entirely via a thermally-driven super-Eddington wind emanating from the core itself, in contrast with the stellar Parker winds usually considered. After a brief cooling phase, the merged cores resume accretion. But, the evolution timescale of transitional discs is too brief for the cores to acquire sufficiently massive envelopes to undergo runaway accretion despite their large combined masses. Consequently, these events lead to the emergence of CSS without their transformation into gas giants. We show that these results are robust for a wide range of disc densities, grain opacities and silicate abundance in the envelope. Our fiducial case reproduces CSS with heavy (>= 30 M_Earth) cores and less massive (a few M_Earth) sub-critical envelopes. We also investigate the other limiting cases, where continuous mergers of comparable-mass cores yield CSS with wider ranges of core-to-envelope mass ratios and envelope opacities. Our results indicate that it is possible for CSS and Uranus and Neptune to emerge within the framework of well studied processes and they may be more common than previously postulated., Comment: Accepted for publication in MNRAS, 17 pages, 9 figures

The tidal perturbation of embedded protoplanets on their natal disks has been widely attributed to be the cause of gap-ring structures in sub-mm images of protoplanetary disks around T Tauri stars. Numerical simulations of this process have been used to propose scalings of characteristic dust gap width/gap-ring distance with respect to planet mass. Applying such scalings to analyze observed gap samples yields a continuous mass distribution for a rich population of hypothetical planets in the range of several Earth to Jupiter masses. In contrast, the conventional core-accretion scenario of planet formation predicts a bi-modal mass function due to 1) the onset of runaway gas accretion above \sim20 Earth masses and 2) suppression of accretion induced by gap opening. Here we examine the dust disk response to the tidal perturbation of eccentric planets as a possible resolution of this paradox. Based on simulated gas and dust distributions, we show the gap-ring separation of Neptune-mass planets with small eccentricities might become comparable to that induced by Saturn-mass planets on circular orbits. This degeneracy may obliterate the discrepancy between the theoretical bi-modal mass distribution and the observed continuous gap width distribution. Despite damping due to planet-disk interaction, modest eccentricity may be sustained either in the outer regions of relatively thick disks or through resonant excitation among multiple super Earths. Moreover, the ring-like dust distribution induced by planets with small eccentricities is axisymmetric even in low viscosity environments, consistent with the paucity of vortices in ALMA images., Comment: 13 pages, 7 figures, ApJ in press

Anomalies among the daughter nuclei of the extinct short-lived radionuclides (SLRs) in the calcium-aluminum-rich inclusions (CAIs) indicate that the Solar System must have been born near a source of the SLRs so that they could be incorporated before they decayed away. $\gamma$-rays from one such living SLR, $^{26}$Al, are detected in only a few nearby star-forming regions. Here we employ multi-wavelength observations to demonstrate that one such region, Ophiuchus, containing many pre-stellar cores that may serve as analogs for the emerging Solar System, is inundated with $^{26}$Al from the neighboring Upper-Scorpius association, and so may provide concrete guidance for how SLR enrichment proceeded in the Solar System complementary to the meteoritics. We demonstrate via Bayesian forward modeling drawing on a wide range of observational and theoretical results that this $^{26}$Al likely 1) arises from supernova explosions, 2) arises from multiple stars, 3) has enriched the gas prior to the formation of the cores, and 4) gives rise to a broad distribution of core enrichment spanning about two orders of magnitude. This means that if the spread in CAI ages is small, as it is in the Solar System, protoplanetary disks must suffer a global heating event., Comment: Published in Nature Astronomy, August 16, 2021

A significant fraction of white dwarfs (WDs) exhibit signs of ongoing accretion of refractory elements at rates $\sim10^3$--$10^7$ kg s$^{-1}$, among which, 37 WDs were detected to harbor dusty debris disks. Such a concurrence requires not only fertile reservoirs of planetary material, but also a high duty cycle of metal delivery. It has been commonly suggested that this material could be supplied by Solar System analogs of Main Belt asteroids or Kuiper Belt objects. Here we consider the primary progenitors of WD pollutants as a population of residual high-eccentricity planetesimals, de-volatilized during the stellar giant phases. Equivalent to the Solar System's long-period comets, they are scattered to the proximity of WDs by perturbations from remaining planets, Galactic tides, passing molecular clouds, and nearby stars. These objects undergo downsizing when they venture within the tidal disruption limit. We show quantitatively how the breakup condition and fragment sizes are determined by material strength and gravity. Thereafter, the fragments' semi-major axes need to decay by at least $\sim$6 orders of magnitude before their constituents are eventually accreted onto the surface of WDs. We investigate the orbital evolution of these fragments around WDs and show that WDs' magnetic fields induce an Alfv\'en-wave drag during their periastron passages and rapidly circularize their orbits. This process could be responsible for the observed accretion rates of heavy-elements and the generation of circum-WD debris disks. A speculative implication is that giant planets may be common around WDs' progenitors and they may still be bound to some WDs today., Comment: Accepted for publication in ApJ. 16 pages, 4 figures and 1 table

There is a population of stars with velocities in excess of 500 km s$^{-1}$ relative to the Galactic center. Many, perhaps most, of these hyper-velocity stars (HVSs) are B stars, similar to the disk and S stars in a nuclear cluster around a super-massive black hole (SMBH) near $\rm Sgr~A^{\star}$. In the paper I of this series, we showed that the eccentricity of the stars emerged from a hypothetical disk around the SMBH can be rapidly excited by the secular perturbation of its intermediate-mass companion (IMC), and we suggested IRS 13E as a potential candidate for the IMC. Here, we show that this process leads to an influx of stars on parabolic orbits to the proximity of $\rm Sgr~A^{\star}$ on a secular timescale of a few Myr. This timescale is much shorter than the diffusion timescale into the lost cone through either the classical or the resonant relaxation. Precession of the highly-eccentric stars' longitude of periastron, relative to that of the IMC, brings them to its proximity within a few Myr. The IMC's gravitational perturbation scatters a fraction of the stars from nearly parabolic to hyperbolic orbits, with respect to the SMBH. Their follow-up close encounters with the SMBH induce them to escape with hyper-velocity. This scenario is a variant of the hypothesis proposed by Hills based on the anticipated breakup of some progenitor binary stars in the proximity of the SMBH, and its main objective is to account for the limited lifespan of the known HVSs. We generalize our previous numerical simulations of this process with a much wider range of orbital configurations. We demonstrate the robustness and evaluate the efficiency of this channel of HVS formation. From these numerical simulations, we infer observable kinematic properties for the HVSs., Comment: 17 pages, 10 figures, to be published in ApJ

We present 2,241 exoplanet candidates identified with data from the Transiting Exoplanet Survey Satellite (TESS) during its two-year prime mission. We list these candidates in the TESS Objects of Interest (TOI) Catalog, which includes both new planet candidates found by TESS and previously-known planets recovered by TESS observations. We describe the process used to identify TOIs and investigate the characteristics of the new planet candidates, and discuss some notable TESS planet discoveries. The TOI Catalog includes an unprecedented number of small planet candidates around nearby bright stars, which are well-suited for detailed follow-up observations. The TESS data products for the Prime Mission (Sectors 1-26), including the TOI Catalog, light curves, full-frame images, and target pixel files, are publicly available on the Mikulski Archive for Space Telescopes., Comment: 39 pages, 16 figures. The Prime Mission TOI Catalog is included in the ancillary data as a CSV. For the most up-to-date catalog, refer to https://tess.mit.edu/toi-releases/

There is a dense group of OB and Wolf-Rayet stars within a fraction of a parsec from the super-massive black hole (SMBH) at the Galactic Center. These stars appear to be coeval and relatively massive. A subgroup of these stars orbits on the same plane. If they emerged with low to modest eccentricity orbits from a common gaseous disk around the central super-massive black hole, their inferred lifespan would not be sufficiently long to account for the excitation of their high orbital eccentricity through dynamical relaxation. Here we analyze the secular perturbation on Galactic Center stars by an intermediate-mass companion (IMC) as a potential mechanism to account for these young disk stars' high eccentricity. This IMC may be either an intermediate-mass black hole (IMBH) or a compact cluster such as IRS-13E. If its orbital angular momentum vector is anti-parallel to that of the disk stars, this perturbation would be effective in exciting the eccentricity of stars with orbital precession rates which resonate with IMC's precession rate. If it orbits around the SMBH in the same direction as the disk stars, the eccentricity of the young stars can still be highly excited by the IMC during the depletion of their natal disk, possible associated with the launch of the Fermi bubble. In this scenario, IMC's precession rate decreases and its secular resonance sweeps through the proximity of the young stars. We carry out numerical simulations with various inclination angles between the orbits of IMC and the disk stars and show this secular interaction is a robust mechanism to excite the eccentricity and inclination of some disk stars., Comment: 21 pages, 7 figures, accepted by ApJ

After protoplanets have acquired sufficient mass to open partial gaps in their natal protostellar disks, residual gas continues to diffuse onto horseshoe streamlines under effect of viscous dissipation, and meander in and out of the planets' Hill sphere. Within the Hill sphere, the horseshoe streamlines intercept gas flow in circumplanetary disks. The host stars' tidal perturbation induces a barrier across the converging streamlines' interface. Viscous transfer of angular momentum across this tidal barrier determines the rate of mass diffusion from the horseshoe streamlines onto the circumplanetary disks, and eventually the accretion rate onto the protoplanets. We carry out a series of numerical simulations to test the influence of this tidal barrier on super thermal planets. In weakly viscous disks, protoplanets' accretion rate steeply decreases with their masses above the thermal limit. As their growth timescale exceeds the gas depletion time scale, their masses reach asymptotic values comparable to that of Jupiter. In relatively thick and strongly viscous disks, protoplanets' asymptotic masses exceed several times that of Jupiter. Two dimensional numerical simulations show that such massive protoplanets strongly excite the eccentricity of nearby horseshoe streamlines, destabilize orderly flow, substantially enhance the diffusion rate across the tidal barrier, and elevate their growth rate until their natal disk is severely depleted. In contrast, eccentric streamlines remain stable in three dimensional simulations. Based on the upper falloff in the observe mass distribution of known exoplanets, we suggest their natal disks had relatively low viscosity alpha sim 0.001, modest thickness H/R sim 0.03 to 0.05, and limited masses comparable to that of minimum mass solar nebula model., Comment: Accepted by ApJ

Active Galactic Nuclei are powered by geometrically-thin accretion disks surrounding a central supermassive black hole. Here we explore the evolution of stars embedded in these extreme astrophysical environments (AGN stars). Because AGN disks are much hotter and denser than the interstellar medium, AGN stars are subject to very different boundary conditions than normal stars. They are also strongly affected by both mass accretion, which can runaway given the vast mass of the disk, and mass loss due to super-Eddington winds. Moreover, chemical mixing plays a critical role in the evolution of these stars by allowing fresh hydrogen accreted from the disk to mix into their cores. We find that, depending on the local AGN density and sound speed and the duration of the AGN phase, AGN stars can rapidly become very massive (M > 100 M$_\odot$). These stars undergo core-collapse, leave behind compact remnants and contribute to polluting the disk with heavy elements. We show that the evolution of AGN stars can have a profound impact on the evolution of AGN metallicities, as well as the production of gravitational waves sources observed by LIGO-Virgo. We point to our galactic center as a region well-suited to test some of our predictions of this exotic stellar evolutionary channel., Comment: 26 pages, 16 figures. Accepted by ApJ

During their formation, emerging protoplanets tidally interact with their natal disks. Proto-gas-giant planets, with Hills radius larger than the disk thickness, open gaps and quench gas flow in the vicinity of their orbits. It is usually assumed that their type II migration is coupled to the viscous evolution of the disk. Although this hypothesis provides an explanation for the origin of close-in planets, it also encounter predicament on the retention of long-period orbits for most gas giant planets. Moreover, numerical simulations indicate that planets migrations are not solely determined by the viscous diffusion of their natal disk. Here we carry out a series of hydrodynamic simulations combined with analytic studies to examine the transition between different paradigms of type II migration. We find a range of planetary mass for which gas continues to flow through a severely depleted gap so that the surface density distribution in the disk region beyond the gap is maintained in a quasi-steady state. The associated gap profile modifies the location of corotation \& Lindblad resonances. In the proximity of the planet's orbit, high-order Lindblad \& corotation torque are weakened by the gas depletion in the gap while low-order Lindblad torques near the gap walls preserves their magnitude. Consequently, the intrinsic surface density distribution of the disk determines delicately both pace and direction of planets' type II migration. We show that this effect might stall the inward migration of giant planets and preserve them in disk regions where the surface density is steep., Comment: 21 pages, 12 figures, accepted by ApJ

The omnipresence of super-Earths suggests that they are able to be retained in natal disks around low-mass stars, whereas exoplanets' mass distribution indicates that some cores have transformed into gas giants through runaway gas accretion at 1AU from solar-type stars. In this paper, we show that transition to runaway gas accretion by cores may be self-impeded by an increase of the grain opacity in their envelope after they have acquired sufficient mass (typically 10Mearth) to enter a pebble-isolation phase. The accumulation of mm-m size pebbles in their migration barriers enhances their local fragmentation rates. The freshly produced sub-mm grains pass through the barrier, elevate the effective dust opacity and reduce the radiative flux in the cores envelope. These effects alone are adequate to suppress the transition to runaway accretion and preserve super-Earths in the stellar proximity (0.1 AU), albeit entropy advection between the envelope and the disk can further reduce the accretion rate. At intermediate distance (1AU) from their host stars, the escalation in the dust opacity dominates over entropy advection in stalling the transition to runaway accretion for marginally pebble-isolated cores. Beyond a few AU, the transformation of more massive cores to gas giants is reachable before severe depletion of disk gas. This requirement can be satisfied either in extended disks with large scale height via orderly accretion of migrating pebbles or through the mergers of oligarchic protoplanetary embryos, and can account for the correlated occurrence of long-period gas giants and close-in super-Earths., Comment: 23 pages, 10 figures, accepted by ApJ

The first discovered interstellar object (ISO), `Oumuamua (1I/2017 U1) shows a dry and rocky surface, an unusually elongated short-to-long axis ratio $c/a \lesssim 1/6$, a low velocity relative to the local standard of rest ($\sim 10$ km s$^{-1}$), non-gravitational accelerations, and tumbles on a few hours timescale. The inferred number density ($\sim 3.5 \times 10^{13} - 2 \times 10^{15}$ pc$^{-3}$) for a population of asteroidal ISOs outnumbers cometary ISOs by $\geq 10^3$, in contrast to the much lower ratio ($\lesssim 10^{-2}$) of rocky/icy Kuiper belt objects. Although some scenarios can cause the ejection of asteroidal ISOs, a unified formation theory has yet to comprehensively link all `Oumuamua's puzzling characteristics and to account for the population. Here we show by numerical simulations that `Oumuamua-like ISOs can be prolifically produced through extensive tidal fragmentation and ejected during close encounters of their volatile-rich parent bodies with their host stars. Material strength enhanced by the intensive heating during periastron passages enables the emergence of extremely elongated triaxial ISOs with shape $c/a \lesssim 1/10$, sizes $a \sim 100$ m, and rocky surfaces. Although volatiles with low sublimation temperature (such as CO) are concurrently depleted, H$_2$O buried under surfaces is preserved in these ISOs, providing an outgassing source without measurable cometary activities for `Oumuamua's non-gravitational accelerations during its passage through the inner Solar System. We infer that the progenitors of `Oumuamua-like ISOs may be km-sized long-period comets from Oort clouds, km-sized residual planetesimals from debris disks, or planet-size bodies at a few AU, orbiting around low-mass main-sequence stars or white dwarfs. These provide abundant reservoirs to account for `Oumuamua's occurrence rate., Comment: 33 pages, 3 figures, Published in Nature Astronomy on 13 April 2020 for associated supplementary files, see http://www.nature.com/articles/s41550-020-1065-8

Super-Earths are by far the most dominant type of exoplanet, yet their formation is still not well understood. In particular, planet formation models predict that many of them should have accreted enough gas to become gas giants. Here we examine the role of the protoplanetary disk in the cooling and contraction of the protoplanetary envelope. In particular, we investigate the effects of 1) the thermal state of the disk as set by the relative size of heating by accretion or irradiation, and whether its energy is transported by radiation or convection, and 2) advection of entropy into the outer envelope by disk flows that penetrate the Hill sphere, as found in 3D global simulations. We find that, at 5 and 1 AU, this flow at the level reported in the non-isothermal simulations where it penetrates only to ~ 0.3 times the Hill radius has little effect on the cooling rate since most of the envelope mass is concentrated close to the core, and far from the flow. On the other hand, at 0.1 AU, the envelope quickly becomes fully-radiative, nearly isothermal, and thus cannot cool down, stalling gas accretion. This effect is significantly more pronounced in convective disks, leading to envelope mass orders of magnitude lower. Entropy advection at 0.1 AU in either radiative or convective disks could therefore explain why super-Earths failed to undergo runaway accretion. These results highlight the importance of the conditions and energy transport in the protoplanetary disk for the accretion of planetary envelopes., Comment: 10 pages, 6 figures, accepted for publication in MNRAS

We investigate the impact of a highly eccentric 10 $M_{\rm \oplus}$ (where $M_{\rm \oplus}$ is the Earth mass) planet embedded in a dusty protoplanetary disk on the dust dynamics and its observational implications. By carrying out high-resolution 2D gas and dust two-fluid hydrodynamical simulations, we find that the planet's orbit can be circularized at large radii. After the planet's orbit is circularized, partial gap opening and dust ring formation happen close to the planet's circularization radius, which can explain the observed gaps/rings at the outer region of disks. When the disk mass and viscosity become low, we find that an eccentric planet can even open gaps and produce dust rings close to the pericenter and apocenter radii before its circularization. This offers alternative scenarios for explaining the observed dust rings and gaps in protoplanetary disks. A lower disk viscosity is favored to produce brighter rings in observations. An eccentric planet can also potentially slow down the dust radial drift in the outer region of the disk when the disk viscosity is low ($\alpha \lesssim2\times10^{-4}$) and the circularization is faster than the dust radial drift., Comment: 11 pages, accepted for publication in ApJ

Galactic center black holes appear to be nearly universally surrounded by dense stellar clusters. When these black holes go through an active accretion phase, the multiple components of the accretion disk, stellar cluster, and black hole system all coexist. We analyze the effect of drag forces on highly eccentric stellar orbits incurred as stars puncture through the disk plane. Disk crossings dissipate orbital energy, drawing eccentric stars into more circular orbits. For high surface density disks, such as those found around black holes accreting near the Eddington mass accretion limit, the magnitude of this energy dissipation can be larger than the mean scatterings that stars receive by two body relaxation. One implication of this is the presence of a disk "loss cone" for highly eccentric stellar orbits where the dissipation from disk interaction outweighs scatter via two body relaxation. The disk loss cone is larger than the tidal disruption loss cone for near-Eddington black hole accretion rates. Stellar orbits within the disk loss cone are lost from the overall cluster as stellar orbits are circularized and stars are potentially ablated by their high-velocity impacts with the disk. We find, however, that the presence of the disk loss cone has a minimal effect on the overall rate of stellar tidal disruptions. Stars are still efficiently fed to the black hole from more-distant stellar orbits that receive large-enough per orbit scatter to jump over the disk loss cone and end up tidally disrupted., Comment: Submitted to AAS Journals; comments are welcome. Associated software at: https://github.com/morganemacleod/NSC_dynamics

Chondrules are silicate spheroids found in meteorites, serving as important fossil records of the early solar system. In order to form chondrules, chondrule precursors must be heated to temperatures much higher than the typical conditions in the current asteroid belt. One proposed mechanism for chondrule heating is the passage through bow shocks of highly eccentric planetesimals in the protoplanetary disk in the early solar system. However, it is difficult for planetesimals to gain and maintain such high eccentricities. In this paper, we present a new scenario in which planetesimals in the asteroid belt region are excited to high eccentricities by the Jovian sweeping secular resonance in a depleting disk, leading to efficient formation of chondrules. We study the orbital evolution of planetesimals in the disk using semi-analytic models and numerical simulations. We investigate the dependence of eccentricity excitation on the planetesimal's size as well as the physical environment, and calculate the probability for chondrule formation. We find that 50 - 2000 km planetesimals can obtain eccentricities larger than 0.6 and cause effective chondrule heating. Most chondrules form in high velocity shocks, in low density gas, and in the inner disk. The fraction of chondrule precursors which become chondrules is about 4 - 9 % between 1.5 - 3 AU. Our model implies that the disk depletion timescale is $\tau_\mathrm{dep}\approx 1~\mathrm{Myr}$, comparable to the age spread of chondrules; and that Jupiter formed before chondrules, no more than 0.7 Myr after the formation of the CAIs., Comment: Accepted for publication by ApJ

During a visual search through the Kepler main-field lightcurves, we have discovered a cataclysmic variable (CV) that experienced only a single 4-day long outburst over four years, rising to three times the quiescent flux. During the four years of non-outburst data the Kepler photometry of KIC 5608384 exhibits ellipsoidal light variations (`ELV') with a $\sim$12% amplitude and period of 8.7 hours. Follow-up ground-based spectral observations have yielded a high-quality radial velocity curve and the associated mass function. Additionally, H$\alpha$ emission lines were present in the spectra even though these were taken while the source was presumably in quiescence. These emission lines are at least partially eclipsed by the companion K star. We utilize the available constraints of the mass function, the ELV amplitude, Roche-lobe filling condition, and inferred radius of the K star to derive the system masses and orbital inclination angle: $M_{\rm wd} \simeq 0.46 \pm 0.02 \, M_\odot$, $M_{\rm K} \simeq 0.41 \pm 0.03 \, M_\odot$, and $i \gtrsim 70^\circ$. The value of $M_{\rm wd}$ is the lowest reported for any accreting WD in a cataclysmic variable. We have also run binary evolution models using MESA to infer the most likely parameters of the pre-cataclysmic binary. Using the mass-transfer rates from the model evolution tracks we conclude that although the rates are close to the critical value for accretion disk stability, we expect KIC 5608384 to exhibit dwarf nova outbursts. We also conclude that the accreting white dwarf most likely descended from a hot subdwarf and, most notably, that this binary is one of the first bona fide examples of a progenitor of AM CVn binaries to have evolved through the CV channel., Comment: 15 pages, 13 figures, 4 tables; accepted for publication in MNRAS

We investigate a scenario where the formation of Globular Clusters (GCs) is triggered by high-speed collisions between infalling atomic-cooling subhalos during the assembly of the main galaxy host, a special dynamical mode of star formation that operates at high gas pressures and is intimately tied to LCDM hierarchical galaxy assembly. The proposed mechanism would give origin to "naked" globulars, as colliding dark matter subhalos and their stars will simply pass through one another while the warm gas within them clashes at highly supersonic speed and decouples from the collisionless component, in a process reminiscent of the Bullet galaxy cluster. We find that the resulting shock-compressed layer cools on a timescale that is typically shorter than the crossing time, first by atomic line emission and then via fine-structure metal-line emission, and is subject to gravitational instability and fragmentation. Through a combination of kinetic theory approximation and high-resolution $N$-body simulations, we show that this model may produce: (a) a GC number-halo mass relation that is linear down to dwarf galaxy scales and agrees with the trend observed over five orders of magnitude in galaxy mass; (b) a population of old globulars with a median age of 12 Gyr and an age spread similar to that observed; (c) a spatial distribution that is biased relative to the overall mass profile of the host; and (d) a bimodal metallicity distribution with a spread similar to that observed in massive galaxies., Comment: 15 pages, 5 figures, accepted for publication by the Astrophysical Journal

We compare the planet-to-star mass-ratio distribution measured by gravitational microlensing to core accretion theory predictions from population synthesis models. The core accretion theory's runaway gas accretion process predicts a dearth of intermediate-mass giant planets that is not seen in the microlensing results. In particular, the models predict $\sim10\,\times$ fewer planets at mass ratios of $10^{-4} \leq q \leq 4 \times 10^{-4}$ than inferred from microlensing observations. This tension implies that gas giant formation may involve processes that have hitherto been overlooked by existing core accretion models or that the planet-forming environment varies considerably as a function of host-star mass. Variation from the usual assumptions for the protoplanetary disk viscosity and thickness could reduce this discrepancy, but such changes might conflict with microlensing results at larger or smaller mass ratios, or with other observations. The resolution of this discrepancy may have important implications for planetary habitability because it has been suggested that the runaway gas accretion process may have triggered the delivery of water to our inner solar system. So, an understanding of giant planet formation may help us to determine the occurrence rate of habitable planets., Comment: 12 pages, 2 figures, 1 table, accepted for publication in ApJL

Cold clouds in hot media are quickly crushed, shredded, and then accelerated as a result of their interaction with the background gas. The persistence of cold clouds moving at substantial velocities in harsh environments is a common yet puzzling feature of many astrophysical systems, from quasar absorption lines probing galactic halos to clouds of dust passing near Sgr $A^*$. Here we run a set of idealized numerical experiments, subjecting a line of cold clouds at a series of mutual separations to a hot background wind. We find that this stream of clouds is able to shield itself from hydrodynamic destruction by accelerating the hot background material, creating a protective layer of co-moving gas. We write down a simple diffusion equation that reproduces the behavior of the simulations, and we discuss the implications for cosmological gas accretion and G2., Comment: Submitted to AAS Journals, comments welcome

Extended gaps in the debris disks of both Vega and Fomalhaut have been observed. These structures have been attributed to tidal perturbations by multiple super-Jupiter gas giant planets. Within the current observational limits, however, no such massive planets have been detected. Here we propose a less stringent `lone-planet' scenario to account for the observed structure with a single eccentric gas giant and suggest that clearing of these wide gaps is induced by its sweeping secular resonance. During the depletion of the disk gas, the planet's secular resonance propagates inward and clears a wide gap over an extended region of the disk. Although some residual intermediate-size planetesimals may remain in the gap, their surface density is too low to either produce super-Earths or lead to sufficiently frequent disruptive collisions to generate any observable dusty signatures. The main advantage of this lone-planet sweeping-secular-resonance model over the previous multiple gas giant tidal truncation scenario is the relaxed requirement on the number of gas giants. The observationally inferred upper mass limit can also be satisfied provided the hypothetical planet has a significant eccentricity. A significant fraction of solar or more massive stars bear gas giant planets with significant eccentricities. If these planets acquired their present-day kinematic properties prior to the depletion of their natal disks, their sweeping secular resonance would effectively impede the retention of neighboring planets and planetesimals over a wide range of orbital semi-major axes., Comment: 20 pages, 12 figures. Accepted for publication in ApJ

Despite more than 20 years since the discovery of the first gas giant planet with an anomalously large radius, the mechanism for planet inflation remains unknown. Here, we report the discovery of EPIC228754001.01, an inflated gas giant planet found with the NASA K2 Mission, and a revised mass for another inflated planet, K2-97b. These planets reside on ~9 day orbits around host stars which recently evolved into red giants. We constrain the irradiation history of these planets using models constrained by asteroseismology and Keck/HIRES spectroscopy and radial velocity measurements. We measure planet radii of 1.31 +\- 0.11 Rjup and and 1.30 +\- 0.07 Rjup, respectively. These radii are typical for planets receiving the current irradiation, but not the former, zero age main sequence irradiation of these planets. This suggests that the current sizes of these planets are directly correlated to their current irradiation. Our precise constraints of the masses and radii of the stars and planets in these systems allow us to constrain the planetary heating efficiency of both systems as 0.03% +0.03%/-0.02%. These results are consistent with a planet re-inflation scenario, but suggest the efficiency of planet re-inflation may be lower than previously theorized. Finally, we discuss the agreement within 10% of stellar masses and radii, and planet masses, radii, and orbital periods of both systems and speculate that this may be due to selection bias in searching for planets around evolved stars., Comment: 18 pages, 15 figures, accepted to AJ. Figures 11, 12, and 13 are the key figures of the paper

The Kepler mission has discovered that multiple close-in super-Earth planets are common around solar-type stars, but their period ratios do not show strong pile-ups near mean motion resonances (MMRs). One scenario is that super-Earths form in a gas-rich disk, and they interact gravitationally with the surrounding gas, inducing their orbital migration. Disk migration theory predicts, however, that planets would end up at resonant orbits due to their differential migration speed. Motivated by the discrepancy between observation and theory, we seek for a mechanism that moves planets out of resonances. We examine the orbital evolution of planet pairs near the magnetospheric cavity during the gas disk dispersal phase. Our study determines the conditions under which planets can escape resonances. We perform two-planet N-body simulations, varying the planet masses, stellar magnetic field strengths, disk accretion rates and gas disk depletion timescales. As planets migrate outward with the expanding magnetospheric cavity, their dynamical configurations can be rearranged. Migration of planets is substantial (minor) in a massive (light) disk. When the outer planet is more massive than the inner planet, the period ratio of two planets increases through outward migration. On the other hand, when the inner planet is more massive, the final period ratio tends to remain similar to the initial one. Larger stellar magnetic field strengths result in planets stopping their migration at longer periods. We highlight \textit{magnetospheric rebound} as an important ingredient able to reconcile disk migration theory with observations. Even when planets are trapped into MMR during the early gas-rich stage, subsequent cavity expansion would induce substantial changes to their orbits, moving them out of resonance., Comment: 10 pages, 5 figures, accepted for publication in A&A

The distribution of heavy elements is anomalously low in the asteroid main belt region compared with elsewhere in the solar system. Observational surveys also indicate a deficit in the number of small ($ \le 50$~km size) asteroids that is two orders of magnitude lower than what is expected from the single power-law distribution that results from a collisional coagulation and fragmentation equilibrium. Here, we consider the possibility that a major fraction of the original asteroid population may have been cleared out by Jupiter's secular resonance, as it swept through the main asteroid belt during the depletion of the solar nebula. This effect leads to the excitation of the asteroids' orbital eccentricities. Concurrently, hydrodynamic drag and planet-disk tidal interaction effectively damp the eccentricities of sub-100 km-size and of super-lunar-size planetesimals, respectively. These combined effects lead to the asteroids' orbital decay and clearing from the present-day main belt region ($\sim 2.1-3.3$~AU). The intermediate-size (50 to several hundreds of kilometers) planetesimals therefore preferentially remain as main belt asteroids near their birthplaces, with modest asymptotic eccentricities. The smaller asteroids are the fragments of subsequent disruptive collisions at later times as suggested by the present-day asteroid families. This scenario provides a natural explanation for both the observed low surface density and the size distribution of asteroids in the main belt. It also offers an explanation for the confined spatial extent of the terrestrial planet building blocks without the requirement of extensive migration of Jupiter., Comment: 26 pages, 18figures, accepted for publication in the ApJ