Showing 1–50 of 173 results for author: Lin, D N C

Radiation Hydrodynamic Simulations of Massive Stars in Gas-rich Environments: Accretion of AGN Stars Suppressed By Thermal Feedback

Authors: Yi-Xian Chen, Yan-Fei Jiang, Jeremy Goodman, Douglas N. C. Lin

Abstract: 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 proce… ▽ More 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/τ$ 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/τ<c_s$, due to adiabatic accretion, the stellar envelope expands quickly to become hydrostatic and further net accretion occurs on thermal timescales in the absence of self-gravity. When the radiation entropy of the medium is less than that of the star, however, this hydrostatic envelope can become more massive than the star itself. Within this sub-regime, self-gravity of the envelope excites runaway growth. Applying our results to realistic environments, moderately massive stars ($\lesssim 100M_\odot$) embedded in AGN disks typically accrete in the fast-diffusion regime, leading to reduction of steady-state accretion rate 1-2 orders of magnitudes lower than expected by previous 1D calculations and $R_{\rm crit}$ smaller than the disk scale height, except in the opacity window at temperature $T\sim 2000$K. Accretion in slow diffusion regime occurs in regions with very high density $ρ\gtrsim 10^{-9}$g/cm$^3$, and needs to be treated with caution in 1D long-term calculations. △ Less

Submitted 21 August, 2024; originally announced August 2024.

Comments: Accepted to ApJ, 20 pages, 16 figures

Evolution of Extremely Soft Binaries in Dense Star Clusters: On the Jupiter Mass Binary Objects

Authors: Yihan Wang, Rosalba Perna, Zhaohuan Zhu, Douglas N. C. Lin

Abstract: 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 dy… ▽ More 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. △ Less

Submitted 29 July, 2024; originally announced July 2024.

Comments: Comments are welcome

Concurrent Accretion and Migration of Giant Planets in their Natal Disks with Consistent Accretion Torque

Authors: Ya-Ping Li, Yi-Xian Chen, Douglas N. C. Lin

Abstract: 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 hydrodyn… ▽ More 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 $α\sim 0.003$ for Jupiter-mass planets. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 17 pages, 9 figures, Accepted by ApJ

White dwarf magnetospheres: Shielding volatile content of icy objects and implications for volatile pollution scarcity

Authors: Wen-Han Zhou, Shang-Fei Liu, Douglas N. C. Lin

Abstract: 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 inve… ▽ More 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. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 13 pages, 11 figures. Accepted for publication in A&A

The Population of Massive Stars in AGN Disks

Authors: Yi-Xian Chen, Douglas N. C. Lin

Abstract: 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… ▽ More 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. △ Less

Submitted 30 May, 2024; v1 submitted 12 April, 2024; originally announced April 2024.

Comments: Published in ApJ

Episodic eruptions of young accreting stars: the key role of disc thermal instability due to Hydrogen ionisation

Authors: Sergei Nayakshin, Fernando Cruz Saenz de Miera, Agnes Kospal, Aleksandra Calovic, Jochen Eisloffel, Douglas N. C. Lin

Abstract: 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 interm… ▽ More 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. △ Less

Submitted 25 March, 2024; v1 submitted 7 March, 2024; originally announced March 2024.

Comments: accepted to MNRAS, 18 pages

Dust Accumulation near the Magnetospheric Truncation of Protoplanetary Discs. II. The Effects of Opacity and Thermal Evolution

Authors: Rixin Li, Yi-Xian Chen, Douglas N. C. Lin

Abstract: 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… ▽ More 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. △ Less

Submitted 22 February, 2024; originally announced February 2024.

Comments: Accepted by MNRAS and in production; 11 pages, 8 figures

Magnetic field of gas giant exoplanets and its influence on the retention of their exomoons

Authors: Xing Wei, D. N. C. Lin

Abstract: We study the magnetic and tidal interactions of a gas-giant exoplanet with its host star and with its exomoons, and focus on their retention. We briefly revisit the scaling law for planetary dynamo in terms of its mass, radius and luminosity. Based on the virial theorem, we construct an evolution law for planetary magnetic field and find that its initial entropy is important for the field evolutio… ▽ More We study the magnetic and tidal interactions of a gas-giant exoplanet with its host star and with its exomoons, and focus on their retention. We briefly revisit the scaling law for planetary dynamo in terms of its mass, radius and luminosity. Based on the virial theorem, we construct an evolution law for planetary magnetic field and find that its initial entropy is important for the field evolution of a high-mass planet. We estimate the magnetic torques on orbit arising from the star-planet and planet-moon magnetic interactions, and find that it can compensate tidal torques and bypass frequency valleys where dynamical-tide response is ineffective. For exomoon's retention we consider two situations. In the presence of a circumplanetary disk (CPD), by comparison between CPD's inner and outer radii, we find that planets with too strong magnetic fields or too small distance from its host star tend not to host exomoons. During the subsequent CPD-free evolution, we find, by comparison between planet's spindown and moon's migration timescales, that hot Jupiters with periods of several days are unlikely to retain large exomoons, albeit they could be surrounded by rings from the debris of tidally disrupted moons. In contrast, moons, if formed around warm or cold Jupiters, can be preserved. Finally, we estimate the radio power and flux density due to the star-planet and planet-moon magnetic interactions and give the upper limit of detection distance by FAST. △ Less

Submitted 11 March, 2024; v1 submitted 11 February, 2024; originally announced February 2024.

Chaotic Type I Migration in Turbulent Discs

Authors: Yinhao Wu, Yi-Xian Chen, Douglas N. C. Lin

Abstract: By performing global hydrodynamical simulations of accretion discs with driven turbulence models, we demonstrate that elevated levels of turbulence induce highly stochastic migration torques on low-mass companions embedded in these discs. This scenario applies to planets migrating within gravito-turbulent regions of protoplanetary discs as well as stars and black holes embedded in the outskirts of… ▽ More By performing global hydrodynamical simulations of accretion discs with driven turbulence models, we demonstrate that elevated levels of turbulence induce highly stochastic migration torques on low-mass companions embedded in these discs. This scenario applies to planets migrating within gravito-turbulent regions of protoplanetary discs as well as stars and black holes embedded in the outskirts of active galactic nuclei (AGN) accretion discs. When the turbulence level is low, linear Lindblad torques persists in the background of stochastic forces and its accumulative effect can still dominate over relatively long timescales. However, in the presence of very stronger turbulence, classical flow patterns around the companion embedded in the disc are disrupted, leading to significant deviations from the expectations of classical Type I migration theory over arbitrarily long timescales. Our findings suggest that the stochastic nature of turbulent migration can prevent low-mass companions from monotonically settling into universal migration traps within the traditional laminar disc framework, thus reducing the frequency of three-body interactions and hierarchical mergers compared to previously expected. We propose a scaling for the transition mass ratio from classical to chaotic migration $q\propto α_R$, where $α_R$ is the Reynolds viscosity stress parameter, which can be further tested and refined by conducting extensive simulations over the relevant parameter space. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 6 pages, 7 figures, accepted by MNRAS Letters. Welcome any comments and suggestions!

3D Global Simulations of Accretion onto Gap-opening Planets: Implications for Circumplanetary Disc Structures and Accretion Rates

Authors: Ya-Ping Li, Yi-Xian Chen, Douglas N. C. Lin

Abstract: 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 inflo… ▽ More 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. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: 19 pages, 17 figures, Accepted for publication in MNRAS

Changing-Look AGN Behaviour Induced by Disk-Captured Tidal Disruption Events

Authors: Yihan Wang, Douglas N. C. Lin, Bing Zhang, Zhaohuan Zhu

Abstract: 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 pro… ▽ More 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. △ Less

Submitted 29 September, 2023; originally announced October 2023.

Comments: Submitted to ApJL. Comments are welcome

The impermanent fate of massive stars in AGN disks

Authors: Mohamad Ali-Dib, Douglas N. C. Lin

Abstract: 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… ▽ More 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 α 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-α and α-chain reactions generate α 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 △ Less

Submitted 18 September, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 17 pages, 7 figures, Accepted for publication in MNRAS

Metal Enrichment due to Embedded Stars in AGN Discs

Authors: Jiamu Huang, Douglas N. C. Lin, Gregory Shields

Abstract: 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 α abundance, though to a lesser degree than previously reported. We als… ▽ More 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 α 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. △ Less

Submitted 30 August, 2023; originally announced August 2023.

Comments: 17 pages, 8 figures, accepted for publication in MNRAS

Stellar/BH Population in AGN Disks: Direct Binary Formation from Capture Objects in Nuclei Clusters

Authors: Yihan Wang, Zhaohuan Zhu, Douglas N. C. Lin

Abstract: 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 numer… ▽ More 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. △ Less

Submitted 26 August, 2024; v1 submitted 17 August, 2023; originally announced August 2023.

Comments: Typo corrections

A massive hot Jupiter orbiting a metal-rich early-M star discovered in the TESS full frame images

Authors: Tianjun Gan, Charles Cadieux, Farbod Jahandar, Allona Vazan, Sharon X. Wang, Shude Mao, Jaime A. Alvarado-Montes, D. N. C. Lin, Étienne Artigau, Neil J. Cook, René Doyon, Andrew W. Mann, Keivan G. Stassun, Adam J. Burgasser, Benjamin V. Rackham, Steve B. Howell, Karen A. Collins, Khalid Barkaoui, Avi Shporer, Jerome de Leon, Luc Arnold, George R. Ricker, Roland Vanderspek, David W. Latham, Sara Seager , et al. (19 additional authors not shown)

Abstract: Observations and statistical studies have shown that giant planets are rare around M dwarfs compared with Sun-like stars. The formation mechanism of these extreme systems remains under debate for decades. With the help of the TESS mission and ground based follow-up observations, we report the discovery of TOI-4201b, the most massive and densest hot Jupiter around an M dwarf known so far with a rad… ▽ More Observations and statistical studies have shown that giant planets are rare around M dwarfs compared with Sun-like stars. The formation mechanism of these extreme systems remains under debate for decades. With the help of the TESS mission and ground based follow-up observations, we report the discovery of TOI-4201b, the most massive and densest hot Jupiter around an M dwarf known so far with a radius of $1.22\pm 0.04\ R_J$ and a mass of $2.48\pm0.09\ M_J$, about 5 times heavier than most other giant planets around M dwarfs. It also has the highest planet-to-star mass ratio ($q\sim 4\times 10^{-3}$) among such systems. The host star is an early-M dwarf with a mass of $0.61\pm0.02\ M_{\odot}$ and a radius of $0.63\pm0.02\ R_{\odot}$. It has significant super-solar iron abundance ([Fe/H]=$0.52\pm 0.08$ dex). However, interior structure modeling suggests that its planet TOI-4201b is metal-poor, which challenges the classical core-accretion correlation of stellar-planet metallicity, unless the planet is inflated by additional energy sources. Building on the detection of this planet, we compare the stellar metallicity distribution of four planetary groups: hot/warm Jupiters around G/M dwarfs. We find that hot/warm Jupiters show a similar metallicity dependence around G-type stars. For M dwarf host stars, the occurrence of hot Jupiters shows a much stronger correlation with iron abundance, while warm Jupiters display a weaker preference, indicating possible different formation histories. △ Less

Submitted 13 September, 2023; v1 submitted 13 July, 2023; originally announced July 2023.

Comments: 24 pages, 14 figures, 4 tables, accepted to AJ

The Orbital Structure and Selection Effects of the Galactic Center S-Star Cluster

Authors: Andreas Burkert, Stefan Gillessen, Douglas N. C. Lin, Xiaochen Zheng, Philipp Schoeller, Frank Eisenhauer, Reinhard Genzel

Abstract: 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 fou… ▽ More 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}$. △ Less

Submitted 3 June, 2023; originally announced June 2023.

Comments: 15 pages, 4 figures, submitted to ApJ, comments very welcome

Chaotic Gas Accretion by Black Holes Embedded in AGN Discs as Cause of Low-spin Signatures in Gravitational Wave Events

Authors: Yi-Xian Chen, Douglas N. C. Lin

Abstract: 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 large… ▽ More 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. △ Less

Submitted 29 March, 2023; originally announced March 2023.

Comments: accepted by MNRAS, 11 pages, 8 figures

Dynamical Evolution of Closely Packed Multiple Planetary Systems Subject to Atmospheric Mass-Loss

Authors: S. Wang, D. N. C. Lin

Abstract: A gap in exoplanets' radius distribution has been widely attributed to the photo-evaporation threshold of their progenitors' gaseous envelope. Giant impacts can also lead to substantial mass-loss. The outflowing gas endures tidal torque from the planets and their host stars. Alongside the planet-star tidal and magnetic interaction, this effect leads to planets' orbital evolution. In multiple super… ▽ More A gap in exoplanets' radius distribution has been widely attributed to the photo-evaporation threshold of their progenitors' gaseous envelope. Giant impacts can also lead to substantial mass-loss. The outflowing gas endures tidal torque from the planets and their host stars. Alongside the planet-star tidal and magnetic interaction, this effect leads to planets' orbital evolution. In multiple super-Earth systems, especially in those which are closely spaced and/or contain planets locked in mean motion resonances (MMRs), modest mass-loss can lead to dynamical instabilities. In order to place some constraints on the extent of planets' mass-loss, we study the evolution of a series of idealized systems of multiple planets with equal masses and a general scaled separation. We consider mass-loss from one or more planets either in the conservative limit or with angular momentum loss from the system. We show that the stable preservation of idealized multiple planetary systems requires either a wide initial separation or a modest upper limit in the amount of mass-loss. This constraint is stringent for the multiple planetary systems in compact and resonant chains. Perturbation due to either impulsive giant impacts between super-Earths or greater than a few percent mass-loss can lead to dynamical instabilities. △ Less

Submitted 1 March, 2023; originally announced March 2023.

Comments: 21 pages, 13 figures, accepted for publication in AJ

Tidal disruption of stellar clusters and their remnants' spatial distribution near the galactic center

Authors: Long Wang, D. N. C. Lin

Abstract: The accretion of massive star clusters via dynamical friction has previously been established to be a likely scenario for the build up of nuclear stellar clusters (NSCs). A remaining issue is whether strong external tidal perturbation may lead to the severe disruption of loosely-bound clusters well before they sink deeply into the center of their host galaxies. We carry out a series of $N$-body si… ▽ More The accretion of massive star clusters via dynamical friction has previously been established to be a likely scenario for the build up of nuclear stellar clusters (NSCs). A remaining issue is whether strong external tidal perturbation may lead to the severe disruption of loosely-bound clusters well before they sink deeply into the center of their host galaxies. We carry out a series of $N$-body simulations and verify our early idealized analytic models. We show if the density profile of the host galaxies can be described by a power-law distribution with an index, $α<1$, the cluster would be compressed in the radial direction by the external galactic tidal field. In contrast, the galactic tidal perturbation is disruptive in regions with a steep, $α>1$, density fall-off or in the very center where gravity is dominated by the point-mass potential of super-massive black holes (SMBHs). This sufficient criterion supplements the conventional necessary Roche-lobe-filling condition in determining the preservation versus disintegration of satellite stellar systems. We simulate the disruption of stellar clusters which venture on nearly-circular, modestly- or highly-eccentric orbits into the center of galaxies with a range of background density profiles and SMBHs. We obtain the spatial distribution of the stellar-cluster remnants. We apply these results to the NSC within a few parsecs from SMBH Sgr A$^\ast$ at the Galactic Center. Recent observations indicate the coexistence of two populations of stars with distinctively separate ages and metallicities. We verify that the subsolar-metalicity population can be the debris of disrupted stellar clusters. △ Less

Submitted 18 December, 2022; originally announced December 2022.

Comments: 21 pages, 18 figures, accepted for publication in ApJ

Hydrodynamical Simulations of Circumbinary Accretion: Balance between Heating and Cooling

Authors: Hai-Yang Wang, Xue-Ning Bai, Dong Lai, Douglas N. C. Lin

Abstract: 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-dimensio… ▽ More 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$-β$ cooling ansatz, where $β$ 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 ($β\gtrsim 0.1$) leads to binary inspiral as opposed to outspiral and the CBD cavity becomes more symmetric. When adopting a dynamically varying $α-$viscosity, binary inspiral only occurs within a narrow range of cooling time (corresponding to $β$ around 0.5). △ Less

Submitted 19 September, 2023; v1 submitted 14 December, 2022; originally announced December 2022.

Comments: 20 pages, 20 figures, Accepted for publication in MNRAS

Rimmed and Rippled Accretion Disc Models to Explain AGN Continuum Lags

Authors: D. A. Starkey, Jiamu Huang, Keith Horne, Douglas N. C. Lin

Abstract: 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… ▽ More 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. △ Less

Submitted 2 December, 2022; originally announced December 2022.

Comments: 17 pages, 7 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

The evolution of circumstellar discs in the Galactic Centre: an application to the G-clouds

Authors: James E. Owen, Douglas N. C. Lin

Abstract: 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 com… ▽ More 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. △ Less

Submitted 30 November, 2022; originally announced December 2022.

Comments: Accepted for publication in MNRAS

Turbulent Transport of Dust Particles in Protostellar Disks: The Effect of Upstream Diffusion

Authors: Tingtao Zhou, Hongping Deng, Yi-Xian Chen, Douglas N. C. Lin

Abstract: 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… ▽ More 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. △ Less

Submitted 14 November, 2022; v1 submitted 19 October, 2022; originally announced October 2022.

Comments: ApJ published, 14 pages, 12 figures

Observable tests for the light-sail scenario of interstellar objects

Authors: Wen-Han Zhou, Shang-Fei Liu, Yun Zhang, Douglas N. C. Lin

Abstract: 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… ▽ More 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σ. 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. △ Less

Submitted 5 September, 2022; v1 submitted 29 August, 2022; originally announced August 2022.

Comments: 17 pages, 13 figures Comments are welcome

Journal ref: A&A 667, A108 (2022)

Blackhole Mergers Through Evection Resonances

Authors: Hareesh Gautham Bhaskar, Gongjie Li, Douglas N. C. Lin

Abstract: 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, includin… ▽ More 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}$. △ Less

Submitted 18 October, 2022; v1 submitted 14 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ApJ

Spin Evolution of Stellar-mass Black Holes Embedded in AGN disks: Orbital Eccentricity Produces Retrograde Circumstellar Flows

Authors: Ya-Ping Li, Yi-Xian Chen, Douglas N. C. Lin, Zhuoxiao Wang

Abstract: 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… ▽ More 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. △ Less

Submitted 25 March, 2022; v1 submitted 10 March, 2022; originally announced March 2022.

Comments: ApJL published

Secular Spin-orbit Resonances of Black Hole Binaries in AGN Disks

Authors: Gongjie Li, Hareesh Gautham Bhaskar, Bence Kocsis, Douglas N. C. Lin

Abstract: 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 warpi… ▽ More 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 $α$-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. △ Less

Submitted 12 April, 2023; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: 11 pages, 13 figures. Accepted for publication in ApJ

Dust Accumulation near the Magnetospheric Truncation of Protoplanetary Discs around T Tauri Stars

Authors: Rixin Li, Yi-Xian Chen, Douglas N. C. Lin

Abstract: 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 establis… ▽ More 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. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: Accepted by MNRAS and in production; 20 pages, 14 figures

The effect of late giant collisions on the atmospheres of protoplanets and the formation of cold sub-Saturns

Authors: Mohamad Ali-Dib, Andrew Cumming, Douglas N. C. Lin

Abstract: 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 gia… ▽ More 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. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Comments: Accepted for publication in MNRAS, 17 pages, 9 figures

Wide Dust Gaps in Protoplanetary Disks Induced by Eccentric Planets: A Mass-Eccentricity Degeneracy

Authors: Yi-Xian Chen, Zhuoxiao Wang, Ya-Ping Li, Clément Baruteau, Douglas N. C. Lin

Abstract: 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 sa… ▽ More 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. △ Less

Submitted 14 October, 2021; v1 submitted 7 September, 2021; originally announced September 2021.

Comments: 13 pages, 7 figures, ApJ in press

A Solar System formation analogue in the Ophiuchus star-forming complex

Authors: John C. Forbes, João Alves, Douglas N. C. Lin

Abstract: 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. $γ$-rays from one such living SLR, $^{26}$Al, are detected in only a few nearby star-forming regions. Here we employ multi-waveleng… ▽ More 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. $γ$-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. △ Less

Submitted 20 August, 2021; originally announced August 2021.

Comments: Published in Nature Astronomy, August 16, 2021

Journal ref: Nature Astronomy (2021)

Orbital migration and circularization of tidal debris by Alfvén-wave drag: circumstellar debris and pollution around white dwarfs

Authors: Yun Zhang, Shang-Fei Liu, Douglas N. C. Lin

Abstract: 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… ▽ More 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én-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. △ Less

Submitted 13 May, 2021; originally announced May 2021.

Comments: Accepted for publication in ApJ. 16 pages, 4 figures and 1 table

The influence of the secular perturbation of an intermediate-mass companion: II. Ejection of hypervelocity stars from the Galactic Center

Authors: Xiaochen Zheng, Douglas N. C. Lin, Shude Mao

Abstract: 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 hypothetic… ▽ More 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. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 17 pages, 10 figures, to be published in ApJ

The TESS Objects of Interest Catalog from the TESS Prime Mission

Authors: Natalia M. Guerrero, S. Seager, Chelsea X. Huang, Andrew Vanderburg, Aylin Garcia Soto, Ismael Mireles, Katharine Hesse, William Fong, Ana Glidden, Avi Shporer, David W. Latham, Karen A. Collins, Samuel N. Quinn, Jennifer Burt, Diana Dragomir, Ian Crossfield, Roland Vanderspek, Michael Fausnaugh, Christopher J. Burke, George Ricker, Tansu Daylan, Zahra Essack, Maximilian N. Günther, Hugh P. Osborn, Joshua Pepper , et al. (80 additional authors not shown)

Abstract: 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 t… ▽ More 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. △ Less

Submitted 24 March, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

Comments: 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/

Departure from the Exact Location of Mean Motion Resonances Induced by the Gas Disk in the Systems Observed by Kepler

Authors: Su Wang, D. N. C. Lin, Xiaochen Zheng, Jianghui Ji

Abstract: The statistical results of transiting planets show that there are two peaks around 1.5 and 2.0 in the distribution of orbital period ratios. A large number of planet pairs are found near the exact location of mean motion resonances (MMRs). In this work, we find out that the depletion and structures of gas disk play crucial roles in driving planet pairs out of exact location of MMRs. Under such sce… ▽ More The statistical results of transiting planets show that there are two peaks around 1.5 and 2.0 in the distribution of orbital period ratios. A large number of planet pairs are found near the exact location of mean motion resonances (MMRs). In this work, we find out that the depletion and structures of gas disk play crucial roles in driving planet pairs out of exact location of MMRs. Under such scenario, planet pairs are trapped into exact MMRs during orbital migration firstly and keep migrating in a same pace. The eccentricities can be excited. Due to the existence of gas disk, eccentricities can be damped leading to the change of orbital period. It will make planet pairs depart from the exact location of MMRs. With depletion timescales larger than 1 Myr, near MMRs configurations are formed easily. Planet pairs have higher possibilities to escape from MMRs with higher disk aspect ratio. Additionally, with weaker corotation torque, planet pairs can depart farther from exact location of MMRs. The final location of the innermost planets in systems are directly related to the transition radius from optically thick region to inner optically thin disk. While the transition radius is smaller than 0.2 AU at the late stage of star evolution process, the innermost planets can reach around 10 days. Our formation scenario is a possible mechanism to explain the formation of near MMRs configuration with the innermost planet farther than 0.1 AU. △ Less

Submitted 11 December, 2020; originally announced December 2020.

Comments: 18 pages, 7 figures, accepted for publication in AJ

The influence of the secular perturbation of an intermediate-mass companion: I. Eccentricity excitation of disk stars at the Galactic center

Authors: Xiaochen Zheng, Douglas N. C. Lin, Shude Mao

Abstract: 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 lifes… ▽ More 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. △ Less

Submitted 9 November, 2020; originally announced November 2020.

Comments: 21 pages, 7 figures, accepted by ApJ

Accretion of Gas Giants Constrained by the Tidal Barrier

Authors: Ya-Ping Li, Yi-Xian Chen, Douglas N. C. Lin, Xiaojia Zhang

Abstract: 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 barri… ▽ More 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. △ Less

Submitted 5 November, 2020; originally announced November 2020.

Comments: Accepted by ApJ

The secular evolution of a uniform density star cluster immersed in a compressible galactic tidal field

Authors: P. B. Ivanov, D. N. C. Lin

Abstract: Nuclear stellar clusters are common in the center of galaxies. We consider the possibility that their progenitors assumed to be globular clusters may have formed elsewhere, migrated to and assembled near their present location. The main challenge for this scenario is whether globular clusters can withstand the tidal field of their host galaxies. Our analysis suggests that provided the mass-density… ▽ More Nuclear stellar clusters are common in the center of galaxies. We consider the possibility that their progenitors assumed to be globular clusters may have formed elsewhere, migrated to and assembled near their present location. The main challenge for this scenario is whether globular clusters can withstand the tidal field of their host galaxies. Our analysis suggests that provided the mass-density distribution of background potential is relatively shallow, as in some galaxies with relatively flat surface brightness profiles, the tidal field near the center of galaxies may be shown to be able to compress rather than disrupt a globular cluster at a distance from the center much smaller than the conventionally defined `tidal disruption radius', $r_t$. To do so, we adopt a previously constructed formalism and consider the secular evolution of star clusters with a homogeneous mass density distribution. We analytically solve the secular equations in the limit that the mass density of stars in the galactic center approaches a uniform distribution. Our model indicates that a star cluster could travel to distances much smaller than $r_t$ without disruption, thus potentially contributing to the formation of the nuclear cluster. However, appropriate numerical N-body simulations are needed to confirm our analytic findings. △ Less

Submitted 28 September, 2020; originally announced September 2020.

Comments: Accepted for publication in The Astrophysical Journal

Stellar Evolution in AGN Disks

Authors: Matteo Cantiello, Adam S. Jermyn, Douglas N. C. Lin

Abstract: 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 strongl… ▽ More 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. △ Less

Submitted 7 March, 2021; v1 submitted 8 September, 2020; originally announced September 2020.

Comments: 26 pages, 16 figures. Accepted by ApJ

Making massive stars in the Galactic Centre via accretion onto low-mass stars within an accretion disc

Authors: Melvyn B. Davies, Doug N. C. Lin

Abstract: The origin of the population of very massive stars observed within $\sim 0.4$ pc of the supermassive black hole in the Galactic Centre is a mystery. Tidal forces from the black hole would likely inhibit {\it in situ} star formation whilst the youth of the massive stars would seem to exclude formation elsewhere followed by transportation (somehow) into the Galactic centre. Here we consider a third… ▽ More The origin of the population of very massive stars observed within $\sim 0.4$ pc of the supermassive black hole in the Galactic Centre is a mystery. Tidal forces from the black hole would likely inhibit {\it in situ} star formation whilst the youth of the massive stars would seem to exclude formation elsewhere followed by transportation (somehow) into the Galactic centre. Here we consider a third way to produce these massive stars from the lower-mass stars contained in the nuclear stellar cluster which surrounds the supermassive black hole. A passing gas cloud can be tidally shredded by the supermassive black hole forming an accretion disc around the black hole. Stars embedded within this accretion disc will accrete gas from the disc via Bondi-Hoyle accretion, where the accretion rate onto a star, $\dot{M}_\star \propto M_\star^2$. This super-exponential growth of accretion can lead to a steep increase in stellar masses, reaching the required 40-50 M$_\odot$ in some cases. The mass growth rate depends sensitively on the stellar orbital eccentricities and their inclinations. The evolution of the orbital inclinations and/or their eccentricities as stars are trapped by the disc, and their orbits are circularised, will increase the number of massive stars produced. Thus accretion onto low-mass stars can lead to a top heavy stellar mass function in the Galactic Centre and other galactic nuclei. The massive stars produced will pollute the environment via supernova explosions and potentially produce compact binaries whose mergers may be detectable by the LIGO-VIRGO gravitational waves observatories. △ Less

Submitted 23 August, 2020; originally announced August 2020.

Comments: Accepted for publication in MNRAS

Retention of Long-Period Gas Giant Planets: Type II Migration Revisited

Authors: Yi-Xian Chen, Xiaojia Zhang, Ya-Ping Li, Hui Li, Douglas N. C. Lin

Abstract: 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… ▽ More 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. △ Less

Submitted 29 July, 2020; originally announced July 2020.

Comments: 21 pages, 12 figures, accepted by ApJ

The Preservation of Super Earths and the Emergence of Gas Giants after Their Progenitor Cores have Entered the Pebble Isolation Phase

Authors: Yi-Xian Chen, Ya-Ping Li, Hui Li, Douglas N. C. Lin

Abstract: 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… ▽ More 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. △ Less

Submitted 21 May, 2020; originally announced May 2020.

Comments: 23 pages, 10 figures, accepted by ApJ

Tidal fragmentation as the origin of 1I/2017 U1 ('Oumuamua)

Authors: Yun Zhang, Douglas N. C. Lin

Abstract: 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… ▽ More 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. △ Less

Submitted 15 April, 2020; originally announced April 2020.

Comments: 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

A Shining Death of Unequal Supermassive Black Hole Binaries

Authors: Xian Chen, D. N. C. Lin, Xiaojia Zhang

Abstract: In the $Λ$CDM scenario, small galaxies merge to produce larger entities. Since supermassive black holes (SMBHs) are found in galaxies of all sizes, SMBH binaries (SMBHBs) are generally expected to form during the amalgamation of galaxies. It is unclear what fraction of these binaries could eventually merge, but a general consensus is that initially the orbital decay is mediated by the surrounding… ▽ More In the $Λ$CDM scenario, small galaxies merge to produce larger entities. Since supermassive black holes (SMBHs) are found in galaxies of all sizes, SMBH binaries (SMBHBs) are generally expected to form during the amalgamation of galaxies. It is unclear what fraction of these binaries could eventually merge, but a general consensus is that initially the orbital decay is mediated by the surrounding gas and stars. In this Letter, we show that in active galactic nulcei (AGNs) the radiation field also causes the orbits of the accreting SMBHs to shrink. The corresponding mechanism, known as the "Poynting-Robertson drag" (PR drag), takes effect on a well-defined timescale $CT_{\rm Sal}$, where $T_{\rm Sal}$ is the Salpeter timescale of the AGN, presumably coincide with the primary SMBH, and $C=4ξ^{-1}ε^{-1}q^{1/3}(1+q)^{2/3}(1-ε)$ is a constant determined by the radiative efficiency $ε$, the mass ratio $q$ of the two black holes, and a parameter $ξ$ characterizing the size of the circum-secondary accretion disk. We find that when $q\lesssim$a few$\times10^{-5}$, the PR drag is more efficient in shrinking the binary than many other mechanisms, such as dynamical friction and type-I migration. Our finding points to a possible new channel for the coalescence of unequal SMBHBs and the clearing of intermediate-massive black holes in AGNs. △ Less

Submitted 26 March, 2020; originally announced March 2020.

Comments: 6 pages, 2 figures, accepted for publication in ApJ Letters

The Imprint of the Protoplanetary Disk in the Accretion of Super-Earth Envelopes

Authors: Mohamad Ali-Dib, Andrew Cumming, Douglas N. C. Lin

Abstract: 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 the… ▽ More 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. △ Less

Submitted 25 March, 2020; v1 submitted 3 March, 2020; originally announced March 2020.

Comments: 10 pages, 6 figures, accepted for publication in MNRAS

On the Dust Signatures Induced by Eccentric Super-Earths in Protoplanetary Disks

Authors: Ya-Ping Li, Hui Li, Shengtai Li, Douglas N. C. Lin

Abstract: 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 cir… ▽ More 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 ($α\lesssim2\times10^{-4}$) and the circularization is faster than the dust radial drift. △ Less

Submitted 7 October, 2019; originally announced October 2019.

Comments: 11 pages, accepted for publication in ApJ

The Effect of Star--Disk Interactions on Highly Eccentric Stellar Orbits in Active Galactic Nuclei: A Disk Loss Cone and Implications for Stellar Tidal Disruption Events

Authors: Morgan MacLeod, Douglas N. C. Lin

Abstract: 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 diss… ▽ More 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. △ Less

Submitted 15 January, 2020; v1 submitted 20 September, 2019; originally announced September 2019.

Comments: Submitted to AAS Journals; comments are welcome. Associated software at: https://github.com/morganemacleod/NSC_dynamics

Chondrule Formation by the Jovian Sweeping Secular Resonance

Authors: Munan Gong, Xiaochen Zheng, Douglas N. C. Lin, Kedron Silsbee, Clement Baruteau, Shude Mao

Abstract: 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 protoplanet… ▽ More 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 $τ_\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. △ Less

Submitted 26 August, 2019; originally announced August 2019.

Comments: Accepted for publication by ApJ

Journal ref: The Astrophysical Journal, Volume 883, Issue 2, article id. 164, 12 pp. (2019)

A 9-Hr CV With One Outburst in 4 Years of Kepler Data

Authors: Zhifei Yu, John Thorstensen, Saul Rappaport, Andrew Mann, Thomas Jacobs, Lorne Nelson, Boris T. Gaensicke, Daryll LaCourse, Tamás Borkovits, Joshua Aiken, Daniel Steeghs, Odette Toloza, Andrew Vanderburg, Douglas N. C. Lin

Abstract: 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.… ▽ More 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$α$ 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. △ Less

Submitted 5 August, 2019; originally announced August 2019.

Comments: 15 pages, 13 figures, 4 tables; accepted for publication in MNRAS

Globular Cluster Formation from Colliding Substructure

Authors: Piero Madau, Alessandro Lupi, Juerg Diemand, Andreas Burkert, Douglas N. C. Lin

Abstract: 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,… ▽ More 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. △ Less

Submitted 15 January, 2020; v1 submitted 22 May, 2019; originally announced May 2019.

Comments: 15 pages, 5 figures, accepted for publication by the Astrophysical Journal