Showing 1–50 of 324 results for author: Lin, D

Constraints on Relativistic Jets from the Fast X-ray Transient 210423 using Prompt Radio Follow-Up Observations

Authors: Dina Ibrahimzade, R. Margutti, J. S. Bright, P. Blanchard, K. Paterson, D. Lin, H. Sears, A. Polzin, I. Andreoni, G. Schroeder, K. D. Alexander, E. Berger, D. L. Coppejans, A. Hajela, J. Irwin, T. Laskar, B. D. Metzger, J. C. Rastinejad, L. Rhodes

Abstract: Fast X-ray Transients (FXTs) are a new observational class of phenomena with no clear physical origin. This is at least partially a consequence of limited multi-wavelength follow up of this class of transients in real time. Here we present deep optical ($g-$ and $i-$ band) photometry with Keck, and prompt radio observations with the VLA of FXT 210423 obtained at ${δt \approx 14-36}$ days since the… ▽ More Fast X-ray Transients (FXTs) are a new observational class of phenomena with no clear physical origin. This is at least partially a consequence of limited multi-wavelength follow up of this class of transients in real time. Here we present deep optical ($g-$ and $i-$ band) photometry with Keck, and prompt radio observations with the VLA of FXT 210423 obtained at ${δt \approx 14-36}$ days since the X-ray trigger. We use these multi-band observations, combined with publicly available data sets, to constrain the presence and physical properties of on-axis and off-axis relativistic jets such as those that can be launched by neutron-star mergers and tidal disruption events, which are among the proposed theoretical scenarios of FXTs. Considering a wide range of possible redshifts $z\le3.5$, circumstellar medium (CSM) density $n={10^{-6}-10^{-1}\,\rm{cm^{-3}}}$, isotropic-equivalent jet kinetic energy $E_{k,iso}={10^{48}-10^{55}\,\rm{erg}}$, we find that we can rule out wide jets with opening angle ${θ_{j}=15^{\circ}}$ viewed within ${10^{\circ}}$ off-axis. For more collimated jets (${θ_{j}=3^{\circ}}$) we can only rule out on-axis (${θ_{obs}=0^{\circ}}$) orientations. This study highlights the constraining power of prompt multi-wavelength observations of FXTs discovered in real time by current (e.g., Einstein Probe) and future facilities. △ Less

Submitted 11 July, 2024; v1 submitted 9 July, 2024; originally announced July 2024.

Comments: 14 pages, 6 figures

GRB 221009A/SN 2022xiw: A Supernova Obscured by a Gamma-Ray Burst Afterglow?

Authors: De-Feng Kong, Xiang-Gao Wang, WeiKang Zheng, Hou-Jun Lü, L. P. Xin, Da-Bin Lin, Jia-Xin Cao, Ming-Xuan Lu, B. Ren, Edgar P. Vidal, J. Y. Wei, En-Wei Liang, Alexei V. Filippenko

Abstract: We present optical photometry for the afterglow of GRB 221009A, in some respects the most extraordinary gamma-ray burst (GRB) ever observed. Good quality in the R-band light curve is obtained, covering 0.32-19.57 days since the Fermi-GBM trigger. We find that a weak bump emerges fromthe declining afterglow at $t \approx 11$ days; a supernova (SN) may be responsible. We use a smooth broken power-la… ▽ More We present optical photometry for the afterglow of GRB 221009A, in some respects the most extraordinary gamma-ray burst (GRB) ever observed. Good quality in the R-band light curve is obtained, covering 0.32-19.57 days since the Fermi-GBM trigger. We find that a weak bump emerges fromthe declining afterglow at $t \approx 11$ days; a supernova (SN) may be responsible. We use a smooth broken power-law and $^{56}\mathrm{Ni}$ model to fit the light curve. The best-fitting results reveal that the SN ejected a total mass of $M_\mathrm{ej} = 3.70 M_\odot$, a $^{56}\mathrm{Ni}$ mass of $M_\mathrm{Ni} = 0.23 M_\odot$, and a kinetic energy of $E_\mathrm{SN,K} = 2.35 \times 10^{52} \mathrm{erg}$. We also compare GRB 221009A with other GRB-SN events based on a GRB-associated SN sample, and find that only SN 2003lw and SN 2011kl can be obviously revealed in the afterglow of GRB 221009A by setting these objects at its distance. This suggests that a supernova (SN 2022xiw) is possibly obscured by the brighter afterglow emission from GRB 221009A. △ Less

Submitted 30 June, 2024; originally announced July 2024.

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

Revisiting the Constraint on Equation of State of Neutron Star based on the Binary Neutron Star Mergers

Authors: Yun-Peng Li, Zhi-Lin Chen, Da-Bin Lin, En-Wei Liang

Abstract: The merger of neutron star (NS)-NS binary can form different production of the compact remnant, among which the supramassive NS (SMNS) could create an internal plateau and the followed steep decay marks the collapse of the SMNS. The proportion of SMNS and the corresponding collapse-time are often used to constrain the NS equation of state (EoS). This paper revisits this topic by considering the ef… ▽ More The merger of neutron star (NS)-NS binary can form different production of the compact remnant, among which the supramassive NS (SMNS) could create an internal plateau and the followed steep decay marks the collapse of the SMNS. The proportion of SMNS and the corresponding collapse-time are often used to constrain the NS equation of state (EoS). This paper revisits this topic by considering the effect of an accretion disk on the compact remnant, which is not considered in previous works. Compared with previous works, the collapse-time distribution (peaks $\sim$100 s) of the SMNSs formed from NS-NS merger is almost unaffected by the initial surface magnetic ($B_{{\rm s},i}$) of NS, but the total energy output of the magnetic dipole radiation from the SMNSs depends on $B_{{\rm s},i}$ significantly. Coupling the constraints from the SMNS fraction, we exclude some EoSs and obtain three candidate EoSs, i.e., DD2, ENG, and MPA1. By comparing the distributions of the collapse-time and the luminosity of the internal plateau (in the short gamma-ray bursts) for those from observations with those obtained based on the three candidate EoSs, it is shown that only the EoS of ENG is favored. Our sample based on the ENG EOS and a mass distribution motivated by Galactic systems suggests that approximately $99\%$ of NS-NS mergers collapse to form a black hole within $10^7$s. This includes scenarios forming a BH promptly ($36.5\%$), a SMNS ($60.7\%$), or a stable NS that transitions into a BH or a SMNS following accretion ($2.1\%$). It also indicates that the remnants for GW170817 and GW190425, and the second object of GW190814 are more likely to be BHs. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Protoplanetary Disk Polarization at Multiple Wavelengths: Are Dust Populations Diverse?

Authors: Rachel E. Harrison, Zhe-Yu Daniel Lin, Leslie W. Looney, Zhi-Yun Li, Haifeng Yang, Ian Stephens, Manuel Fernández-López

Abstract: Millimeter and sub-millimeter observations of continuum linear dust polarization provide insight into dust grain growth in protoplanetary disks, which are the progenitors of planetary systems. We present the results of the first survey of dust polarization in protoplanetary disks at 870 $μ$m and 3 mm. We find that protoplanetary disks in the same molecular cloud at similar evolutionary stages can… ▽ More Millimeter and sub-millimeter observations of continuum linear dust polarization provide insight into dust grain growth in protoplanetary disks, which are the progenitors of planetary systems. We present the results of the first survey of dust polarization in protoplanetary disks at 870 $μ$m and 3 mm. We find that protoplanetary disks in the same molecular cloud at similar evolutionary stages can exhibit different correlations between observing wavelength and polarization morphology and fraction. We explore possible origins for these differences in polarization, including differences in dust populations and protostar properties. For RY Tau and MWC 480, which are consistent with scattering at both wavelengths, we present models of the scattering polarization from several dust grain size distributions. These models aim to reproduce two features of the observational results for these disks: (1) both disks have an observable degree of polarization at both wavelengths and (2) the polarization fraction is higher at 3 mm than at 870 $μ$m in the centers of the disks. For both disks, these features can be reproduced by a power-law distribution of spherical dust grains with a maximum radius of 200 $μ$m and high optical depth. In MWC 480, we can also reproduce features (1) and (2) with a model containing large grains ($a_{max}$ = 490 $μ$m ) near the disk midplane and small grains ($a_{max}$ = 140 $μ$m) above and below the midplane. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 22 pages, 12 figures

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

Early Planet Formation in Embedded Disks (eDisk) XIII: Aligned Disks with Non-Settled Dust Around the Newly Resolved Class 0 Protobinary R CrA IRAS 32

Authors: Frankie J. Encalada, Leslie W. Looney, Shigehisa Takakuwa, John J. Tobin, Nagayoshi Ohashi, Jes K. Jørgensen, Zhi-Yun Li, Yuri Aikawa, Yusuke Aso, Patrick M. Koch, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Zhe-Yu Daniel Lin, Alejandro Santamarıa-Miranda, Itziar de Gregorio-Monsalvo, Nguyen Thi Phuong, Adele Plunkett, Jinshi Sai, Rajeeb Sharma, Hsi-Wei Yen, Ilseung Han

Abstract: Young protostellar binary systems, with expected ages less than $\sim$10$^5$ years, are little modified since birth, providing key clues to binary formation and evolution. We present a first look at the young, Class 0 binary protostellar system R CrA IRAS 32 from the Early Planet Formation in Embedded Disks (eDisk) ALMA large program, which observed the system in the 1.3 mm continuum emission,… ▽ More Young protostellar binary systems, with expected ages less than $\sim$10$^5$ years, are little modified since birth, providing key clues to binary formation and evolution. We present a first look at the young, Class 0 binary protostellar system R CrA IRAS 32 from the Early Planet Formation in Embedded Disks (eDisk) ALMA large program, which observed the system in the 1.3 mm continuum emission, $^{12}$CO (2-1), $^{13}$CO (2-1), C$^{18}$O (2-1), SO (6$_5$-5$_4$), and nine other molecular lines that trace disk, envelope, shocks, and outflows. With a continuum resolution of $\sim$0.03$^{\prime\prime}$ ($\sim$5 au, at a distance of 150 pc), we characterize the newly discovered binary system with a separation of 207 au, their circumstellar disks, and a circumbinary disk-like structure. The circumstellar disk radii are 26.9$\pm$0.3 and 22.8$\pm$0.3 au for sources A and B, respectively, and their circumstellar disk dust masses are estimated as 22.5$\pm$1.1 and 12.4$\pm$0.6 M$_{\Earth}$. The circumstellar disks and the circumbinary structure have well aligned position angles and inclinations, indicating formation in a smooth, ordered process such as disk fragmentation. In addition, the circumstellar disks have a near/far-side asymmetry in the continuum emission suggesting that the dust has yet to settle into a thin layer near the midplane. Spectral analysis of CO isotopologues reveals outflows that originate from both of the sources and possibly from the circumbinary disk-like structure. Furthermore, we detect Keplerian rotation in the $^{13}$CO isotopologues toward both circumstellar disks and likely Keplerian rotation in the circumbinary structure; the latter suggests that it is probably a circumbinary disk. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Comments: 33 pages, 29 figures, 6 tables

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

Turbulent vortex with moderate dust settling probed by scattering-induced polarization in the IRS 48 system

Authors: Haifeng Yang, Manuel Fernández-López, Zhi-Yun Li, Ian W. Stephens, Leslie W. Looney, Zhe-Yu Daniel Lin, Rachel Harrison

Abstract: We investigate the crescent-shaped dust trap in the transition disk, Oph IRS 48, using well-resolved (sub)millimeter polarimetric observations at ALMA Band 7 (870 $μ$m). The dust polarization map reveals patterns consistent with dust scattering-induced polarization. There is a relative displacement between the polarized flux and the total flux, which holds the key to understanding the dust scale h… ▽ More We investigate the crescent-shaped dust trap in the transition disk, Oph IRS 48, using well-resolved (sub)millimeter polarimetric observations at ALMA Band 7 (870 $μ$m). The dust polarization map reveals patterns consistent with dust scattering-induced polarization. There is a relative displacement between the polarized flux and the total flux, which holds the key to understanding the dust scale heights in this system. We model the polarization observations, focusing on the effects of dust scale heights. We find that the interplay between the inclination-induced polarization and the polarization arising from radiation anisotropy in the crescent determines the observed polarization; the anisotropy is controlled by the dust optical depth along the midplane, which is, in turn, determined by the dust scale height in the vertical direction. We find that the dust grains can neither be completely settled nor well mixed with the gas. The completely settled case produces little radial displacement between the total and polarized flux, while the well-mixed case produces an azimuthal pattern in the outer (radial) edge of the crescent that is not observed. Our best model has a gas-to-dust scale height ratio of 2, and can reproduce both the radial displacement and the azimuthal displacement between the total and polarized flux. We infer an effective turbulence $α$ parameter of approximately $0.0001-0.005$. The scattering-induced polarization provides insight into a turbulent vortex with a moderate level of dust settling in the IRS 48 system, which is hard to achieve otherwise. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 9 pages in the main text, 5 figures, accepted by ApJ

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.

Revisit the Emission Polarization of the Internal-shock for the blazars' Jet

Authors: Hao-Qiang Zhang, Da-Bin Lin, Kuan Liu, En-Wei Liang

Abstract: Recent Imaging X-ray Polarimetry Explorer (IXPE) observations of blazars tend to support the shock model for the X-ray emission, but report a low polarization degree ($Π\sim 10\%$) in X-rays compared with the previous theoretical expectations in the shock model. In order to reconcile the theoretical expectations with observations, we revisit the polarization of the shock emission by considering di… ▽ More Recent Imaging X-ray Polarimetry Explorer (IXPE) observations of blazars tend to support the shock model for the X-ray emission, but report a low polarization degree ($Π\sim 10\%$) in X-rays compared with the previous theoretical expectations in the shock model. In order to reconcile the theoretical expectations with observations, we revisit the polarization of the shock emission by considering different kind of direction distribution for the shock-generated magnetic fields (sgMFs). Here, $w'_{\rm sg}\propto(\sinθ')^{ζ_{\rm sg}}$ with $θ'=0$ along the shock normal direction is used to describe the direction distribution of the sgMFs in the shock co-moving frame. It is found that the polarization in the X-ray and radio emission for a general jet in blazars can be described as $Π\sim 44.5[1-\exp(-ζ_{\rm sg}/2.6)]\%$ and $Π\sim 20[1-\exp(-ζ_{\rm sg}/2.4)]\%$, respectively. Correspondingly, one can have $ζ_{\rm sg}\sim 1-1.5$ according to the IXPE observations. Besides the sgMFs, the magnetic fields generated by the Richmyer-Meshkov instability (rmMFs) is supposed to present in the jets. The direction of the rmMFs is mainly distributed along the shock normal in the simulations and thus $w'_{\rm rm}\propto(\cos θ')^{ζ_{\rm rm}}$ is adopted to describe the direction distribution of rmMFs. We find that the rmMFs is likely to significantly affect the polarization properties at the low-frequency emission, especially when the sgMFs decay rapidly. Based on the contemporaneous radio and X-ray observations, we find the the emission of the electrons in the rmMFs make a significant contribution in the low-frequency emission and the the ordered background magnetic fields (obMFs) can be neglected. △ Less

Submitted 10 February, 2024; originally announced February 2024.

Early Planet Formation in Embedded Disks (eDisk) XIV: Flared Dust Distribution and Viscous Accretion Heating of the Disk around R CrA IRS 7B-a

Authors: Shigehisa Takakuwa, Kazuya Saigo, Miyu Kido, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Yuri Aikawa, Yusuke Aso, Sacha Gavino, Ilseung Han, Patrick M. Koch, Woojin Kwon, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Shoji Mori, Jinshi Sai, Rajeeb Sharma, Patrick Sheehan, Kengo Tomida, Jonathan P. Williams, Yoshihide Yamato, Hsi-Wei Yen

Abstract: We performed radiative transfer calculations and observing simulations to reproduce the 1.3-mm dust-continuum and C$^{18}$O (2-1) images in the Class I protostar R CrA IRS7B-a, observed with the ALMA Large Program ``Early Planet Formation in Embedded Disks (eDisk)". We found that the dust disk model passively heated by the central protostar cannot reproduce the observed peak brightness temperature… ▽ More We performed radiative transfer calculations and observing simulations to reproduce the 1.3-mm dust-continuum and C$^{18}$O (2-1) images in the Class I protostar R CrA IRS7B-a, observed with the ALMA Large Program ``Early Planet Formation in Embedded Disks (eDisk)". We found that the dust disk model passively heated by the central protostar cannot reproduce the observed peak brightness temperature of the 1.3-mm continuum emission ($\sim$195 K), regardless of the assumptions about the dust opacity. Our calculation suggests that viscous accretion heating in the disk is required to reproduce the observed high brightness temperature. The observed intensity profile of the 1.3-mm dust-continuum emission along the disk minor axis is skewed toward the disk far side. Our modeling reveals that such an asymmetric intensity distribution requires flaring of the dust along the disk's vertical direction with the scale-height following $h/r \sim r^{0.3}$ as function of radius. These results are in sharp contrast to those of Class II disks, which show geometrically flat dust distributions and lower dust temperatures. From our modeling of the C$^{18}$O (2-1) emission, the outermost radius of the gas disk is estimated to be $\sim$80 au, larger than that of the dust disk ($\sim$62 au), to reproduce the observed distribution of the C$^{18}$O (2-1) emission in IRS 7B-a. Our modeling unveils a hot and thick dust disk plus a larger gas disk around one of the eDisk targets, which could be applicable to other protostellar sources in contrast to more evolved sources. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 26 pages, 13 figures

Rapid dimming followed by a state transition: a study of the highly variable nuclear transient AT 2019avd over 1000+ days

Authors: Yanan Wang, Dheeraj R. Pasham, Diego Altamirano, Andres Gurpide, Noel Castro Segura, Matthew Middleton, Long Ji, Santiago del Palacio, Muryel Guolo, Poshak Gandhi, Shuang-Nan Zhang, Ronald Remillard, Dacheng Lin, Megan Masterson, Ranieri D. Baldi, Francesco Tombesi, Jon M. Miller, Wenda Zhang, Andrea Sanna

Abstract: The tidal disruption of a star around a supermassive black hole (SMBH) offers a unique opportunity to study accretion onto a SMBH on a human-timescale. We present results from our 1000+ days NICER, Swift and Chandra monitoring campaign of AT 2019avd, a nuclear transient with TDE-like properties. Our primary finding is that approximately 225 days following the peak of X-ray emission, there is a rap… ▽ More The tidal disruption of a star around a supermassive black hole (SMBH) offers a unique opportunity to study accretion onto a SMBH on a human-timescale. We present results from our 1000+ days NICER, Swift and Chandra monitoring campaign of AT 2019avd, a nuclear transient with TDE-like properties. Our primary finding is that approximately 225 days following the peak of X-ray emission, there is a rapid drop in luminosity exceeding two orders of magnitude. This X-ray drop-off is accompanied by X-ray spectral hardening, followed by a 740-day plateau phase. During this phase, the spectral index decreases from 6.2+-1.1 to 2.3+-0.4, while the disk temperature remains constant. Additionally, we detect pronounced X-ray variability, with an average fractional root mean squared amplitude of 47%, manifesting over timescales of a few dozen minutes. We propose that this phenomenon may be attributed to intervening clumpy outflows. The overall properties of AT 2019avd suggest that the accretion disk evolves from a super-Eddington to a sub-Eddington luminosity state, possibly associated with a compact jet. This evolution follows a pattern in the hardness-intensity diagram similar to that observed in stellar-mass black holes, supporting the mass invariance of accretion-ejection processes around black holes. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: 21 pages, 11 figures, 1 table; accepted for publication in ApJ

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!

Aligned Grains and Scattered Light Found in Gaps of Planet-Forming Disk

Authors: Ian W. Stephens, Zhe-Yu Daniel Lin, Manuel Fernandez-Lopez, Zhi-Yun Li, Leslie W. Looney, Haifeng Yang, Rachel Harrison, Akimasa Kataoka, Carlos Carrasco-Gonzalez, Satoshi Okuzumi, Ryo Tazaki

Abstract: Polarized (sub)millimeter emission from dust grains in circumstellar disks was initially thought to be due to grains aligned with the magnetic field. However, higher resolution multi-wavelength observations along with improved models found that this polarization is dominated by self-scattering at shorter wavelengths (e.g., 870 $μ$m) and by grains aligned with something other than magnetic fields a… ▽ More Polarized (sub)millimeter emission from dust grains in circumstellar disks was initially thought to be due to grains aligned with the magnetic field. However, higher resolution multi-wavelength observations along with improved models found that this polarization is dominated by self-scattering at shorter wavelengths (e.g., 870 $μ$m) and by grains aligned with something other than magnetic fields at longer wavelengths (e.g., 3 mm). Nevertheless, the polarization signal is expected to depend on the underlying substructure, and observations hitherto have been unable to resolve polarization in multiple rings and gaps. HL Tau, a protoplanetary disk located 147.3 $\pm$ 0.5 pc away, is the brightest Class I or Class II disk at millimeter/submillimeter wavelengths. Here we show deep, high-resolution 870 $μ$m polarization observations of HL Tau, resolving polarization in both the rings and gaps. We find that the gaps have polarization angles with a significant azimuthal component and a higher polarization fraction than the rings. Our models show that the disk polarization is due to both scattering and emission from aligned effectively prolate grains. The intrinsic polarization of aligned dust grains is likely over 10%, which is much higher than what was expected in low resolution observations (~1%). Asymmetries and dust features are revealed in the polarization observations that are not seen in non-polarimetric observations. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: Published in Nature

Searching for Associations Between Short Gamma-ray Bursts and Fast Radio Burst

Authors: Ming-Xuan Lu, Long Li, Xiang-Gao Wang, Can-Min Deng, Yun-Feng Liang, Da-Bin Lin, En-Wei Liang

Abstract: The physical origin of fast radio bursts (FRBs) is still unclear. However, young magnetars associated with short-duration gamma-ray bursts (SGRBs) have been thought to be possible central engines for some FRBs. In this paper, we perform a systematic search for SGRBs that are associated with FRBs in a sample including 623 FRBs (601 one-off bursts and 22 repeaters) and 168 SGRBs with precise localiz… ▽ More The physical origin of fast radio bursts (FRBs) is still unclear. However, young magnetars associated with short-duration gamma-ray bursts (SGRBs) have been thought to be possible central engines for some FRBs. In this paper, we perform a systematic search for SGRBs that are associated with FRBs in a sample including 623 FRBs (601 one-off bursts and 22 repeaters) and 168 SGRBs with precise localizations. We find that FRB 190309A is spatially associated with GRB 060502B, with a chance probability of 0.05 when temporal and redshift information is taken into account. Considering the high chance probability (the statistical significance is < 3σ), we examine other observational properties such as the host galaxy, the dispersion measure, and the energy budget of the central engine to check the possibility of their association. Although the available observational information is insufficient to determine whether they are physically associated, it does not rule out such a possibility. As the only pair of FRB and GRB that are spatially associated, it remains an interesting case worthy of further attention △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 19 pages, 4 figures, accepted to ApJ

Early Planet Formation in Embedded Disks (eDisk) X: Compact Disks, Extended Infall, and a Fossil Outburst in the Class I Oph IRS43 Binary

Authors: Suchitra Narayanan, Jonathan P. Williams, John J. Tobin, Jes K. Jorgensen, Nagayoshi Ohashi, Zhe-Yu Daniel Lin, Merel L. R. van't Hoff, Zhi-Yun Li, Adele L. Plunkett, Leslie W. Looney, Shigehisa Takakuwa, Hsi-Wei Yen, Yusuke Aso, Christian Flores, Jeong-Eun Lee, Shih-Ping Lai, Woojin Kwon, Itziar de Gregorio-Monsalvo, Rajeeb Sharma, Chang Won Lee

Abstract: We present the first results from the Early Planet Formation in Embedded Disks (eDisk) ALMA Large Program toward Oph IRS43, a binary system of solar mass protostars. The 1.3 mm dust continuum observations resolve a compact disk, ~6au radius, around the northern component and show that the disk around the southern component is even smaller, <~3 au. CO, 13CO, and C18O maps reveal a large cavity in a… ▽ More We present the first results from the Early Planet Formation in Embedded Disks (eDisk) ALMA Large Program toward Oph IRS43, a binary system of solar mass protostars. The 1.3 mm dust continuum observations resolve a compact disk, ~6au radius, around the northern component and show that the disk around the southern component is even smaller, <~3 au. CO, 13CO, and C18O maps reveal a large cavity in a low mass envelope that shows kinematic signatures of rotation and infall extending out to ~ 2000au. An expanding CO bubble centered on the extrapolated location of the source ~130 years ago suggests a recent outburst. Despite the small size of the disks, the overall picture is of a remarkably large and dynamically active region. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: Paper 10 of the ALMA eDisk Large Program. Accepted for publication in ApJ

Early Planet Formation in Embedded Disks (eDisk) XII: Accretion streamers, protoplanetary disk, and outflow in the Class I source Oph IRS63

Authors: Christian Flores, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Shigehisa Takakuwa, Zhi-Yun Li, Zhe-Yu Daniel Lin, Merel L. R. van 't Hoff, Adele L. Plunkett, Yoshihide Yamato, Jinshi Sai, Patrick M. Koch, Hsi-Wei Yen, Yuri Aikawa, Yusuke Aso, Itziar de Gregorio-Monsalvo, Miyu Kido, Woojin Kwon, Jeong-Eun Lee, Chang Won Lee, Leslie W. Looney, Alejandro Santamaría-Miranda, Rajeeb Sharma, Travis J. Thieme, Jonathan P. Williams , et al. (3 additional authors not shown)

Abstract: We present ALMA observations of the Class I source Oph IRS63 in the context of the Early Planet Formation in Embedded Disks (eDisk) large program. Our ALMA observations of Oph IRS63 show a myriad of protostellar features, such as a shell-like bipolar outflow (in $^{12}$CO), an extended rotating envelope structure (in $^{13}$CO), a streamer connecting the envelope to the disk (in C$^{18}$O), and se… ▽ More We present ALMA observations of the Class I source Oph IRS63 in the context of the Early Planet Formation in Embedded Disks (eDisk) large program. Our ALMA observations of Oph IRS63 show a myriad of protostellar features, such as a shell-like bipolar outflow (in $^{12}$CO), an extended rotating envelope structure (in $^{13}$CO), a streamer connecting the envelope to the disk (in C$^{18}$O), and several small-scale spiral structures seen towards the edge of the dust continuum (in SO). By analyzing the velocity pattern of $^{13}$CO and C$^{18}$O, we measure a protostellar mass of $\rm M_\star = 0.5 \pm 0.2 $~$\rm M_\odot$ and confirm the presence of a disk rotating at almost Keplerian velocity that extends up to $\sim260$ au. These calculations also show that the gaseous disk is about four times larger than the dust disk, which could indicate dust evolution and radial drift. Furthermore, we model the C$^{18}$O streamer and SO spiral structures as features originating from an infalling rotating structure that continuously feeds the young protostellar disk. We compute an envelope-to-disk mass infall rate of $\sim 10^{-6}$~$\rm M_\odot \, yr^{-1}$ and compare it to the disk-to-star mass accretion rate of $\sim 10^{-8}$~$\rm M_\odot \, yr^{-1}$, from which we infer that the protostellar disk is in a mass build-up phase. At the current mass infall rate, we speculate that soon the disk will become too massive to be gravitationally stable. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 26 pages and 17 figures

Early Planet Formation in Embedded Disks (eDisk). VIII. A Small Protostellar Disk around the Extremely Low-Mass and Young Class 0 Protostar, IRAS 15398-3359

Authors: Travis J. Thieme, Shih-Ping Lai, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Jinshi Sai, Yusuke Aso, Jonathan P. Williams, Yoshihide Yamato, Yuri Aikawa, Itziar de Gregorio-Monsalvo, Ilseung Han, Woojin Kwon, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Suchitra Narayanan, Nguyen Thi Phuong, Adele L. Plunkett, Alejandro Santamaría-Miranda, Rajeeb Sharma, Shigehisa Takakuwa, Hsi-Wei Yen

Abstract: Protostellar disks are a ubiquitous part of the star formation process and the future sites of planet formation. As part of the Early Planet Formation in Embedded Disks (eDisk) large program, we present high-angular resolution dust continuum ($\sim40\,$mas) and molecular line ($\sim150\,$mas) observations of the Class 0 protostar, IRAS 15398-3359. The dust continuum is small, compact, and centrall… ▽ More Protostellar disks are a ubiquitous part of the star formation process and the future sites of planet formation. As part of the Early Planet Formation in Embedded Disks (eDisk) large program, we present high-angular resolution dust continuum ($\sim40\,$mas) and molecular line ($\sim150\,$mas) observations of the Class 0 protostar, IRAS 15398-3359. The dust continuum is small, compact, and centrally peaked, while more extended dust structures are found in the outflow directions. We perform a 2D Gaussian fitting to find the deconvolved size and $2σ$ radius of the dust disk to be $4.5\times2.8\,\mathrm{au}$ and $3.8\,\mathrm{au}$, respectively. We estimate the gas+dust disk mass assuming optically thin continuum emission to be $0.6-1.8\,M_\mathrm{jup}$, indicating a very low-mass disk. The CO isotopologues trace components of the outflows and inner envelope, while SO traces a compact, rotating disk-like component. Using several rotation curve fittings on the PV diagram of the SO emission, the lower limits of the protostellar mass and gas disk radius are $0.022\,M_\odot$ and $31.2\,\mathrm{au}$ from our Modified 2 single power-law fitting. A conservative upper limit of the protostellar mass is inferred to be $0.1\,M_\odot$. The protostellar mass-accretion rate and the specific angular momentum at the protostellar disk edge are found to be between $1.3-6.1\times10^{-6}\,M_\odot\,\mathrm{yr^{-1}}$ and $1.2-3.8\times10^{-4}\,\mathrm{km\,s^{-1}\,pc}$, respectively, with an age estimated between $0.4-7.5\times10^{4}\,$yr. At this young age with no clear substructures in the disk, planet formation would likely not yet have started. This study highlights the importance of high-resolution observations and systematic fitting procedures when deriving dynamical properties of deeply embedded Class 0 protostars. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 28 pages, 16 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

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

Panchromatic (Sub)millimeter Polarization Observations of HL Tau Unveil Aligned Scattering Grains

Authors: Zhe-Yu Daniel Lin, Zhi-Yun Li, Ian W. Stephens, Manuel Fernández-López, Carlos Carrasco-González, Claire J. Chandler, Alice Pasetto, Leslie W. Looney, Haifeng Yang, Rachel E. Harrison, Sarah I. Sadavoy, Thomas Henning, A. Meredith Hughes, Akimasa Kataoka, Woojin Kwon, Takayuki Muto, Dominique Segura-Cox

Abstract: Polarization is a unique tool to study the properties of dust grains of protoplanetary disks and detail the initial conditions of planet formation. Polarization around HL Tau was previously imaged using the Atacama Large Millimeter/submillimeter Array (ALMA) at Bands 3 (3.1 mm), 6 (1.3 mm), and 7 (0.87 mm), showing that the polarization orientation changes across wavelength $λ$. The polarization m… ▽ More Polarization is a unique tool to study the properties of dust grains of protoplanetary disks and detail the initial conditions of planet formation. Polarization around HL Tau was previously imaged using the Atacama Large Millimeter/submillimeter Array (ALMA) at Bands 3 (3.1 mm), 6 (1.3 mm), and 7 (0.87 mm), showing that the polarization orientation changes across wavelength $λ$. The polarization morphology at Band 7 is predominantly parallel to the disk minor axis but appears azimuthally oriented at Band 3, with the morphology at Band 6 in between the two. We present new ~0.2" (29 au) polarization observations at Q-Band (7.0 mm) using the Karl G. Jansky Very Large Array (VLA) and at Bands 4 (2.1 mm), 5 (1.5 mm), and 7 using ALMA, consolidating HL Tau's position as the protoplanetary disk with the most complete wavelength coverage in dust polarization. The polarization patterns at Bands 4 and 5 continue to follow the morphological transition with wavelength previously identified in Bands 3, 6, and 7. Based on the azimuthal variation, we decompose the polarization into contributions from scattering ($s$) and thermal emission ($t$). We find that $s$ decreases slowly with increasing $λ$, and $t$ increases more rapidly with $λ$ which are expected from optical depth effects of toroidally aligned, scattering prolate grains. The relatively weak $λ$ dependence of $s$ is consistent with large, porous grains. The sparse polarization detections from the Q-band image are also consistent with toroidally aligned prolate grains. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 22 pages, 13 figures, submitted to MNRAS

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

Early Planet Formation in Embedded Disks (eDisk) VI: Kinematic Structures around the Very Low Mass Protostar IRAS 16253-2429

Authors: Yusuke Aso, Woojin Kwon, Nagayoshi Ohashi, Jes K. Jorgensen, John J. Tobin, Yuri Aikawa, Itziar de Gregorio-Monsalvo, Ilseung Han, Miyu Kido, Patrick M. Koch, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Suchitra Narayanan, Nguyen Thi Phuong, Jinshi Sai, Kazuya Saigo, Alejandro Santamaria-Miranda, Rajeeb Sharma, Shigehisa Takakuwa, Travis J. Thieme, Kengo Tomida , et al. (2 additional authors not shown)

Abstract: Precise estimates of protostellar masses are crucial to characterize the formation of stars of low masses down to brown-dwarfs (BDs; M* < 0.08 Msun). The most accurate estimation of protostellar mass uses the Keplerian rotation in the circumstellar disk around the protostar. To apply the Keplerian rotation method to a protostar at the low-mass end, we have observed the Class 0 protostar IRAS 16253… ▽ More Precise estimates of protostellar masses are crucial to characterize the formation of stars of low masses down to brown-dwarfs (BDs; M* < 0.08 Msun). The most accurate estimation of protostellar mass uses the Keplerian rotation in the circumstellar disk around the protostar. To apply the Keplerian rotation method to a protostar at the low-mass end, we have observed the Class 0 protostar IRAS 16253-2429 using the Atacama Large Millimeter/submillimeter Array (ALMA) in the 1.3 mm continuum at an angular resolution of 0.07" (10 au), and in the 12CO, C18O, 13CO (J=2-1), and SO (J_N = 6_5-5_4) molecular lines, as part of the ALMA Large Program Early Planet Formation in Embedded Disks (eDisk). The continuum emission traces a non-axisymmetric, disk-like structure perpendicular to the associated 12CO outflow. The position-velocity (PV) diagrams in the C18O and 13CO lines can be interpreted as infalling and rotating motions. In contrast, the PV diagram along the major axis of the disk-like structure in the 12CO line allows us to identify Keplerian rotation. The central stellar mass and the disk radius are estimated to be ~0.12-0.17 Msun and ~13-19 au, respectively. The SO line suggests the existence of an accretion shock at a ring (r~28 au) surrounding the disk and a streamer from the eastern side of the envelope. IRAS 16253-2429 is not a proto-BD but has a central stellar mass close to the BD mass regime, and our results provide a typical picture of such very low-mass protostars. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: 41 pages, 14 figures

Journal ref: ApJ 954 (2023) 101-117

Early Planet Formation in Embedded Disks (eDisk) IX: High-resolution ALMA Observations of the Class 0 Protostar R CrA IRS5N and its surrounding

Authors: Rajeeb Sharma, Jes K. Jørgensen, Sacha Gavino, Nagayoshi Ohashi, John J. Tobin, Zhe-Yu Daniel Lin, Zhi-Yun Li, Shigehisa Takakuwa, Chang Won Lee, Jinshi Sai, Woojin Kwon, Itziar de Gregorio-Monsalvo, Alejandro Santamaría-Miranda, Hsi-Wei Yen, Yuri Aikawa, Yusuke Aso, Shih-Ping Lai, Jeong-Eun Lee, Leslie W. Looney, Nguyen Thi Phuong, Travis J. Thieme, Jonathan P. Williams

Abstract: We present high-resolution, high-sensitivity observations of the Class 0 protostar RCrA IRS5N as part of the Atacama Large Milimeter/submilimeter Array (ALMA) large program Early Planet Formation in Embedded Disks (eDisk). The 1.3 mm continuum emission reveals a flattened continuum structure around IRS5N, consistent with a protostellar disk in the early phases of evolution. The continuum emission… ▽ More We present high-resolution, high-sensitivity observations of the Class 0 protostar RCrA IRS5N as part of the Atacama Large Milimeter/submilimeter Array (ALMA) large program Early Planet Formation in Embedded Disks (eDisk). The 1.3 mm continuum emission reveals a flattened continuum structure around IRS5N, consistent with a protostellar disk in the early phases of evolution. The continuum emission appears smooth and shows no substructures. However, a brightness asymmetry is observed along the minor axis of the disk, suggesting the disk is optically and geometrically thick. We estimate the disk mass to be between 0.007 and 0.02 M$_{\odot}$. Furthermore, molecular emission has been detected from various species, including C$^{18}$O (2$-$1), $^{12}$CO (2$-$1), $^{13}$CO (2$-$1), and H$_2$CO (3$_{0,3}-2_{0,2}$, 3$_{2,1}-2_{2,0}$, and 3$_{2,2}-2_{2,1}$). By conducting a position-velocity analysis of the C$^{18}$O (2$-$1) emission, we find that the disk of IRS5N exhibits characteristics consistent with Keplerian rotation around a central protostar with a mass of approximately 0.3 M$_{\odot}$. Additionally, we observe dust continuum emission from the nearby binary source, IRS5a/b. The emission in $^{12}$CO toward IRS5a/b seems to emanate from IRS5b and flow into IRS5a, suggesting material transport between their mutual orbits. The lack of a detected outflow and large-scale negatives in \tlvco~observed toward IRS5N suggests that much of the flux from IRS5N is being resolved out. Due to this substantial surrounding envelope, the central IRS5N protostar is expected to be significantly more massive in the future. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 25 pages, 24 figures. Accepted for publication in ApJ

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 June, 2024; v1 submitted 17 August, 2023; originally announced August 2023.

Comments: Typo corrections

Early Planet Formation in Embedded Disks (eDisk) V: Possible Annular Substructure in a Circumstellar Disk in the Ced110 IRS4 System

Authors: Jinshi Sai, Hsi-Wei Yen, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Shigehisa Takakuwa, Kazuya Saigo, Yusuke Aso, Zhe-Yu Daniel Lin, Patrick M. Koch, Yuri Aikawa, Christian Flores, Itziar de Gregorio-Monsalvo, Ilseung Han, Miyu Kido, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Leslie W. Looney, Shoji Mori, Nguyen Thi Phuong, Alejandro Santamaría-Miranda, Rajeeb Sharma , et al. (3 additional authors not shown)

Abstract: We have observed the Class 0/I protostellar system Ced110 IRS4 at an angular resolution of $0.05''$ ($\sim$10 au) as a part of the ALMA large program; Early Planet Formation in the Embedded Disks (eDisk). The 1.3 mm dust continuum emission reveals that Ced110 IRS4 is a binary system with a projected separation of $\sim$250 au. The continuum emissions associated with the main source and its compani… ▽ More We have observed the Class 0/I protostellar system Ced110 IRS4 at an angular resolution of $0.05''$ ($\sim$10 au) as a part of the ALMA large program; Early Planet Formation in the Embedded Disks (eDisk). The 1.3 mm dust continuum emission reveals that Ced110 IRS4 is a binary system with a projected separation of $\sim$250 au. The continuum emissions associated with the main source and its companion, named Ced110 IRS4A and IRS4B respectively, exhibit disk-like shapes and likely arise from dust disks around the protostars. The continuum emission of Ced110 IRS4A has a radius of $\sim$91.7 au ($\sim0.485''$), and shows bumps along its major axis with an asymmetry. The bumps can be interpreted as an shallow, ring-like structure at a radius of $\sim$40 au ($\sim0.2''$) in the continuum emission, as demonstrated from two-dimensional intensity distribution models. A rotation curve analysis on the C$^{18}$O and $^{13}$CO $J=2$-1 lines reveals the presence of a Keplerian disk within a radius of 120 au around Ced110 IRS4A, which supports the interpretation that the dust continuum emission arises from a disk. The ring-like structure in the dust continuum emission might indicate a possible, annular substructure in the surface density of the embedded disk, although the possibility that it is an apparent structure due to the optically thick continuum emission cannot be ruled out. △ Less

Submitted 31 August, 2023; v1 submitted 17 July, 2023; originally announced July 2023.

Comments: 32 pages, 23 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

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

Event Rate of Fast Radio Burst from Binary Neutron-star Mergers

Authors: Zhi-Lin Chen, Rui-Chong Hu, Da-Bin Lin, En-Wei Liang

Abstract: It is proposed that one-off fast radio burst (FRB) with periodic structures may be produced during the inspiral phase of a binary neutron-star (BNS) merger. In this paper, we study the event rate of such kind of FRB. We first investigate the properties of two one-off FRBs with periodic structures (i.e., FRB~20191221A and FRB~20210213A) in this scenario, by assuming the fast magnetosonic wave is re… ▽ More It is proposed that one-off fast radio burst (FRB) with periodic structures may be produced during the inspiral phase of a binary neutron-star (BNS) merger. In this paper, we study the event rate of such kind of FRB. We first investigate the properties of two one-off FRBs with periodic structures (i.e., FRB~20191221A and FRB~20210213A) in this scenario, by assuming the fast magnetosonic wave is responsible for their radio emission. For the luminosities and periods of these bursts, it is found that the pre-merger BNS with magnetic field strength $B\gtrsim 10^{12}\,{\rm Gs}$ is required. This is relatively high compared with that of the most of the BNSs observed in our Galaxy, of which the magnetic field is around $10^{9}\,{\rm Gs}$. Since the observed BNSs in our Galaxy are the binaries without suffering merger, a credited event rate of BNS-merger originated FRBs should be estimated by considering the evolution of both the BNS systems and their magnetic fields. Based on the population synthesis and adopting a decaying magnetic field of NSs, we estimate the event rate of BNS-mergers relative to their final magnetic fields. We find that the rapid merged BNSs tend to merge with high magnetization, and the event rate of BNS-merger originated FRBs, i.e., the BNS-mergers with both NSs' magnetic field being higher than $10^{12}\,{\rm Gs}$ is $\sim8\times10^{4}\,\rm{yr}^{-1}$ ($19 \%$ of the total BNS-mergers) in redshift $z<1$. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 11 pages, 2 figures, accepted for publication in Astrophysical Journal

Early Planet Formation in Embedded Disks (eDisk). VII. Keplerian Disk, Disk Substructure, and Accretion Streamers in the Class 0 Protostar IRAS 16544-1604 in CB 68

Authors: Miyu Kido, Shigehisa Takakuwa, Kazuya Saigo, Nagayoshi Ohashi, John J. Tobin, Jes K, Jørgensen, Yuri Aikawa, Yusuke Aso, Frankie J. Encalada, Christian Flores, Sacha Gavino, Itziar de Gregorio-Monsalvo, Ilseung Han, Shingo Hirano, Patrick M. Koch, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Shoji Mori, Suchitra Narayanan , et al. (12 additional authors not shown)

Abstract: We present observations of the Class 0 protostar IRAS 16544-1604 in CB 68 from the ''Early Planet Formation in Embedded Disks (eDisk)'' ALMA Large program. The ALMA observations target continuum and lines at 1.3-mm with an angular resolution of $\sim$5 au. The continuum image reveals a dusty protostellar disk with a radius of $\sim$30 au seen close to edge-on, and asymmetric structures both along… ▽ More We present observations of the Class 0 protostar IRAS 16544-1604 in CB 68 from the ''Early Planet Formation in Embedded Disks (eDisk)'' ALMA Large program. The ALMA observations target continuum and lines at 1.3-mm with an angular resolution of $\sim$5 au. The continuum image reveals a dusty protostellar disk with a radius of $\sim$30 au seen close to edge-on, and asymmetric structures both along the major and minor axes. While the asymmetry along the minor axis can be interpreted as the effect of the dust flaring, the asymmetry along the major axis comes from a real non-axisymmetric structure. The C$^{18}$O image cubes clearly show the gas in the disk that follows a Keplerian rotation pattern around a $\sim$0.14 $M_{\odot}$ central protostar. Furthermore, there are $\sim$1500 au-scale streamer-like features of gas connecting from North-East, North-North-West, and North-West to the disk, as well as the bending outflow as seen in the $^{12}$CO (2-1) emission. At the apparent landing point of NE streamer, there are SO (6$_5$-5$_4$) and SiO (5-4) emission detected. The spatial and velocity structure of NE streamer can be interpreted as a free-falling gas with a conserved specific angular momentum, and the detection of the SO and SiO emission at the tip of the streamer implies presence of accretion shocks. Our eDisk observations have unveiled that the Class 0 protostar in CB 68 has a Keplerian rotating disk with flaring and non-axisymmetric structure associated with accretion streamers and outflows. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 30 pages, 24 figures, accepted for publication in The Astrophysical Journal as one of the first-look papers of the eDisk ALMA Large Program

Early Planet Formation in Embedded Disks (eDisk). II. Limited Dust Settling and Prominent Snow Surfaces in the Edge-on Class I Disk IRAS 04302+2247

Authors: Zhe-Yu Daniel Lin, Zhi-Yun Li, John J. Tobin, Nagayoshi Ohashi, Jes Kristian Jørgensen, Leslie W. Looney, Yusuke Aso, Shigehisa Takakuwa, Yuri Aikawa, Merel L. R. van 't Hoff, Itziar de Gregorio-Monsalvo, Frankie J. Encalada, Christian Flores, Sacha Gavino, Ilseung Han, Miyu Kido, Patrick M. Koch, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Nguyen Thi Phuong, Jinshi Sai, Rajeeb Sharma, Patrick Sheehan , et al. (4 additional authors not shown)

Abstract: While dust disks around optically visible, Class II protostars are found to be vertically thin, when and how dust settles to the midplane are unclear. As part of the Atacama Large Millimeter/submillimeter Array (ALMA) large program, Early Planet Formation in Embedded Disks, we analyze the edge-on, embedded, Class I protostar IRAS 04302+2247, also nicknamed the ``Butterfly Star." With a resolution… ▽ More While dust disks around optically visible, Class II protostars are found to be vertically thin, when and how dust settles to the midplane are unclear. As part of the Atacama Large Millimeter/submillimeter Array (ALMA) large program, Early Planet Formation in Embedded Disks, we analyze the edge-on, embedded, Class I protostar IRAS 04302+2247, also nicknamed the ``Butterfly Star." With a resolution of 0.05" (8~au), the 1.3 mm continuum shows an asymmetry along the minor axis which is evidence of an optically thick and geometrically thick disk viewed nearly edge-on. There is no evidence of rings and gaps, which could be due to the lack of radial substructure or the highly inclined and optically thick view. With 0.1" (16~au) resolution, we resolve the 2D snow surfaces, i.e., the boundary region between freeze-out and sublimation, for $^{12}$CO $J$=2--1, $^{13}$CO $J$=2--1, C$^{18}$O $J$=2--1, $H_{2}$CO $J$=$3_{0,3}$--$2_{0,2}$, and SO $J$=$6_{5}$--$5_{4}$, and constrain the CO midplane snow line to $\sim 130$ au. We find Keplerian rotation around a protostar of $1.6 \pm 0.4 M_{\odot}$ using C$^{18}$O. Through forward ray-tracing using RADMC-3D, we find that the dust scale height is $\sim 6$ au at a radius of 100~au from the central star and is comparable to the gas pressure scale height. The results suggest that the dust of this Class~I source has yet to vertically settle significantly. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 33 pages, 21 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

Early Planet Formation in Embedded Disks (eDisk). IV. The Ringed and Warped Structure of the Disk around the Class I Protostar L1489 IRS

Authors: Yoshihide Yamato, Yuri Aikawa, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Shigehisa Takakuwa, Yusuke Aso, Jinshi Sai, Christian Flores, Itziar de Gregorio-Monsalvo, Shingo Hirano, Ilseung Han, Miyu Kido, Patrick M. Koch, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Shoji Mori, Suchitra Narayanan, Nguyen Thi Phuong, Kazuya Saigo , et al. (6 additional authors not shown)

Abstract: Constraining the physical and chemical structure of young embedded disks is crucial to understanding the earliest stages of planet formation. As part of the Early Planet Formation in Embedded Disks Atacama Large Millimeter/submillimeter Array Large Program, we present high spatial resolution ($\sim$0$.\!\!^{\prime\prime}$1 or $\sim$15 au) observations of the 1.3 mm continuum and $^{13}$CO $J=$ 2-1… ▽ More Constraining the physical and chemical structure of young embedded disks is crucial to understanding the earliest stages of planet formation. As part of the Early Planet Formation in Embedded Disks Atacama Large Millimeter/submillimeter Array Large Program, we present high spatial resolution ($\sim$0$.\!\!^{\prime\prime}$1 or $\sim$15 au) observations of the 1.3 mm continuum and $^{13}$CO $J=$ 2-1, C$^{18}$O $J=$ 2-1, and SO $J_N=$ $6_5$-$5_4$ molecular lines toward the disk around the Class I protostar L1489 IRS. The continuum emission shows a ring-like structure at 56 au from the central protostar and a tenuous, optically thin emission extending beyond $\sim$300 au. The $^{13}$CO emission traces the warm disk surface, while the C$^{18}$O emission originates from near the disk midplane. The coincidence of the radial emission peak of C$^{18}$O with the dust ring may indicate a gap-ring structure in the gaseous disk as well. The SO emission shows a highly complex distribution, including a compact, prominent component at $\lesssim$30 au, which is likely to originate from thermally sublimated SO molecules. The compact SO emission also shows a velocity gradient along a slightly ($\sim15^\circ$) tilted direction with respect to the major axis of the dust disk, which we interpret as an inner warped disk in addition to the warp around $\sim$200 au suggested by previous work. These warped structures may be formed by a planet or companion with an inclined orbit, or by a gradual change in the angular momentum axis during gas infall. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 24 pages, 12 figures. Accepted for publication in The Astrophysical Journal as one of the first-look papers of the eDisk ALMA Large Program

Early Planet Formation in Embedded Disks (eDisk) III: A first high-resolution view of sub-mm continuum and molecular line emission toward the Class 0 protostar L1527 IRS

Authors: Merel L. R. van 't Hoff, John J. Tobin, Zhi-Yun Li, Nagayoshi Ohashi, Jes K. Jørgensen, Zhe-Yu Daniel Lin, Yuri Aikawa, Yusuke Aso, Itziar de Gregorio-Monsalvo, Sacha Gavino, Ilseung Han, Patrick M. Koch, Woojin Kwon, Chang Won Lee, Jeong-Eun Lee, Leslie W. Looney, Suchitra Narayanan, Adele Plunkett, Jinshi Sai, Alejandro Santamaría-Miranda, Rajeeb Sharma, Patrick D. Sheehan, Shigehisa Takakuwa, Travis J. Thieme, Jonathan P. Williams , et al. (3 additional authors not shown)

Abstract: Studying the physical and chemical conditions of young embedded disks is crucial to constrain the initial conditions for planet formation. Here, we present Atacama Large Millimeter/submillimeter Array (ALMA) observations of dust continuum at $\sim$0.06" (8 au) resolution and molecular line emission at $\sim$0.17" (24 au) resolution toward the Class 0 protostar L1527 IRS from the Large Program eDis… ▽ More Studying the physical and chemical conditions of young embedded disks is crucial to constrain the initial conditions for planet formation. Here, we present Atacama Large Millimeter/submillimeter Array (ALMA) observations of dust continuum at $\sim$0.06" (8 au) resolution and molecular line emission at $\sim$0.17" (24 au) resolution toward the Class 0 protostar L1527 IRS from the Large Program eDisk (Early Planet Formation in Embedded Disks). The continuum emission is smooth without substructures, but asymmetric along both the major and minor axes of the disk as previously observed. The detected lines of $^{12}$CO, $^{13}$CO, C$^{18}$O, H$_2$CO, c-C$_3$H$_2$, SO, SiO, and DCN trace different components of the protostellar system, with a disk wind potentially visible in $^{12}$CO. The $^{13}$CO brightness temperature and the H$_2$CO line ratio confirm that the disk is too warm for CO freeze out, with the snowline located at $\sim$350 au in the envelope. Both molecules show potential evidence of a temperature increase around the disk-envelope interface. SO seems to originate predominantly in UV-irradiated regions such as the disk surface and the outflow cavity walls rather than at the disk-envelope interface as previously suggested. Finally, the continuum asymmetry along the minor axis is consistent with the inclination derived from the large-scale (100" or 14,000 au) outflow, but opposite to that based on the molecular jet and envelope emission, suggesting a misalignment in the system. Overall, these results highlight the importance of observing multiple molecular species in multiple transitions to characterize the physical and chemical environment of young disks. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 27 pages, 16 figures, 2 tables, 10 pages appendix with 12 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

Early Planet Formation in Embedded Disks (eDisk). I. Overview of the Program and First Results

Authors: Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Shigehisa Takakuwa, Patrick Sheehan, Yuri Aikawa, Zhi-Yun Li, Leslie W. Looney, Jonathan P. Willians, Yusuke Aso, Rajeeb Sharma, Jinshi Sai, Yoshihide Yamato, Jeong-Eun Lee, Kengo Tomida, Hsi-Wei Yen, Frankie J Encalada, Christian Flores, Sacha Gavino, Miyu Kido, Ilseung Han, Zhe-Yu Daniel Lin, Suchitra Narayanan, Nguyen Thi Phuong, Alejandro Santamaría-Miranda , et al. (12 additional authors not shown)

Abstract: We present an overview of the Large Program, ``Early Planet Formation in Embedded Disks (eDisk)'', conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The ubiquitous detections of substructures, particularly rings and gaps, in protoplanetary disks around T Tauri stars raise the possibility that at least some planet formation may have already started during the embedded stages o… ▽ More We present an overview of the Large Program, ``Early Planet Formation in Embedded Disks (eDisk)'', conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The ubiquitous detections of substructures, particularly rings and gaps, in protoplanetary disks around T Tauri stars raise the possibility that at least some planet formation may have already started during the embedded stages of star formation. In order to address exactly how and when planet formation is initiated, the program focuses on searching for substructures in disks around 12 Class 0 and 7 Class I protostars in nearby ($< $200 pc) star-forming regions through 1.3 mm continuum observations at a resolution of $\sim7$ au (0.04"). The initial results show that the continuum emission, mostly arising from dust disks around the sample protostars, has relatively few distinctive substructures, such as rings and spirals, in marked contrast to Class II disks. The dramatic difference may suggest that substructures quickly develop in disks when the systems evolve from protostars to Class II sources or alternatively that high optical depth of the continuum emission could obscure internal structures. Kinematic information obtained through CO isotopologue lines and other lines reveals the presence of Keplerian disks around protostars, providing us with crucial physical parameters, in particular, the dynamical mass of the central protostars. We describe the background of the eDisk program, the sample selection and their ALMA observations, the data reduction, and also highlight representative first-look results. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: This is a publication of a series of eDisk ALMA large program first-look papers

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

Tidal disruption events and quasi periodic eruptions

Authors: Natalie A. Webb, Didier Barret, Olivier Godet, Maitrayee Gupta, Dacheng Lin, Erwan Quintin, Hugo Tranin

Abstract: Tidal disruption events (TDEs) occur when a star passes close to a massive black hole, so that the tidal forces of the black hole exceed the binding energy of a star and cause it to be ripped apart. Part of the matter will fall onto the black hole, causing a strong increase in the luminosity. Such events are often seen in the optical or the X-ray (or both) or even at other wavelengths such as in t… ▽ More Tidal disruption events (TDEs) occur when a star passes close to a massive black hole, so that the tidal forces of the black hole exceed the binding energy of a star and cause it to be ripped apart. Part of the matter will fall onto the black hole, causing a strong increase in the luminosity. Such events are often seen in the optical or the X-ray (or both) or even at other wavelengths such as in the radio, where the diversity of observed emission is still poorly understood. The XMM-Newton catalogue of approximately a million X-ray detections covering 1283$^2$ degrees of sky contains a number of these events. Here I will show the diverse nature of a number of TDEs discovered in the catalogue and discuss their relationship with quasi periodic eruptions. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: 7 pages, 1 figure, accepted version for the proceedings of the 'Black Hole Accretion Under the X-ray Microscope' Meeting held at ESAC in June 2022. Publisher : Astronomische Nachrichten

Eccentric Dust Ring in the IRS 48 Transition Disk

Authors: Haifeng Yang, Manuel Fernandez-Lopez, Zhi-Yun Li, Ian W. Stephens, Leslie W. Looney, Zhe-Yu Daniel Lin, Rachel Harrison

Abstract: Crescent-shaped structures in transition disks hold the key to studying the putative companions to the central stars. The dust dynamics, especially that of different grain sizes, is important to understanding the role of pressure bumps in planet formation. In this work, we present deep dust continuum observation with high resolution towards the Oph IRS 48 system. For the first time, we are able to… ▽ More Crescent-shaped structures in transition disks hold the key to studying the putative companions to the central stars. The dust dynamics, especially that of different grain sizes, is important to understanding the role of pressure bumps in planet formation. In this work, we present deep dust continuum observation with high resolution towards the Oph IRS 48 system. For the first time, we are able to significantly trace and detect emission along $95\%$ of the ring crossing the crescent-shaped structure. The ring is highly eccentric with an eccentricity of $0.27$. The flux density contrast between the peak of the flux and its counter part along the ring is $\sim 270$. In addition, we detect a compact emission toward the central star. If the emission is an inner circumstellar disk inside the cavity, it has a radius of at most a couple of astronomical units with a dust mass of $1.5\times 10^{-8}\rm\, M_\odot$, or $0.005\rm\, M_\oplus$. We also discuss the implications of the potential eccentric orbit on the proper motion of the crescent, the putative secondary companion, and the asymmetry in velocity maps. △ Less

Submitted 6 April, 2023; originally announced April 2023.

Comments: 9 pages, 5 figures, accepted for publication on AAS Journals

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

Enhanced Blackhole mergers in AGN discs due to Precession induced resonances

Authors: Hareesh Gautham Bhaskar, Gongjie Li, Doug Lin

Abstract: Recent studies have shown that AGN discs can host sources of gravitational waves. Compact binaries can form and merge in AGN discs through their interactions with the gas and other compact objects in the disc. It is also possible for the binaries to shorten the merging timescale due to eccentricity excitation caused by perturbations from the supermassive blackhole (SMBH). In this paper we focus on… ▽ More Recent studies have shown that AGN discs can host sources of gravitational waves. Compact binaries can form and merge in AGN discs through their interactions with the gas and other compact objects in the disc. It is also possible for the binaries to shorten the merging timescale due to eccentricity excitation caused by perturbations from the supermassive blackhole (SMBH). In this paper we focus on effects due to precession-induced (eviction-like) resonances, where nodal and apsidal precession rates of the binary is commensurable with the mean motion of the binary around the SMBH. We focus on intermediate mass black hole (IMBH)-stellar mass black hole (SBH) binaries, and consider binary orbit inclined from the circum-IMBH disk which leads to the orbital $J_2$ precession. We show that if a binary is captured in these resonances and is migrating towards the companion, it can undergo large eccentricity and inclination variations. We derive analytical expressions for the location of fixed points, libration timescale and width for these resonances, and identified two resonances in the near coplanar regime (the evection and eviction resonances) as well as two resonances in the near polar regime that can lead to mergers. We also derive analytical expressions for the maximum eccentricity that a migrating binary can achieve for given initial conditions. Specifically, the maximum eccentricity can reach 0.9 when captured in these resonances before orbital decay due to gravitational wave emission dominates, and the capture is only possible for slow migration ($\sim 10$ Myr) 2-3 order of magnitude longer than the resonance libration timescale. We also show that capture into multiple resonances is possible, and can further excite eccentricities. △ Less

Submitted 16 August, 2023; v1 submitted 20 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

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

Detection of Quasi-periodic Oscillations in SGR 150228213

Authors: Run-Chao Chen, Can-Min Deng, Xiang-Gao Wang, Zi-Min Zhou, Xing Yang, Da-Bin Lin, Qi Wang, En-Wei Liang

Abstract: The detection of quasi-periodic oscillations (QPOs) in magnetar giant flares (GFs) has brought a new perspective to study the mechanism of magnetar bursts. Due to the scarcity of GFs, searching QPOs from magnetar short bursts is reasonable. Here we report the detection of a high frequency QPO at approximately 110 Hz and a wide QPO at approximately 60 Hz in a short magnetar burst SGR 150228213, wit… ▽ More The detection of quasi-periodic oscillations (QPOs) in magnetar giant flares (GFs) has brought a new perspective to study the mechanism of magnetar bursts. Due to the scarcity of GFs, searching QPOs from magnetar short bursts is reasonable. Here we report the detection of a high frequency QPO at approximately 110 Hz and a wide QPO at approximately 60 Hz in a short magnetar burst SGR 150228213, with a confidence level of 3.35$σ$. This burst was initially attributed to 4U 0142+61 by $Fermi$/GBM on location, but we haven't detected such QPOs in other bursts from this magnetar. We also found that there was a repeating fast radio burst associated with SGR 150228213 on location. Finally, we discuss the possible origins of SGR 150228213. △ Less

Submitted 29 May, 2023; v1 submitted 21 February, 2023; originally announced February 2023.

The star formation history of the Sco-Cen association. Coherent star formation patterns in space and time

Authors: Sebastian Ratzenböck, Josefa E. Großschedl, João Alves, Núria Miret-Roig, Immanuel Bomze, John Forbes, Alyssa Goodman, Alvaro Hacar, Doug Lin, Stefan Meingast, Torsten Möller, Martin Piecka, Laura Posch, Alena Rottensteiner, Cameren Swiggum, Catherine Zucker

Abstract: We reconstruct the star formation history of the Sco-Cen OB association using a novel high-resolution age map of the region. We develop an approach to produce robust ages for Sco-Cen's recently identified 37 stellar clusters using the SigMA algorithm. The Sco-Cen star formation timeline reveals four periods of enhanced star formation activity, or bursts, remarkably separated by about 5 Myr. Of the… ▽ More We reconstruct the star formation history of the Sco-Cen OB association using a novel high-resolution age map of the region. We develop an approach to produce robust ages for Sco-Cen's recently identified 37 stellar clusters using the SigMA algorithm. The Sco-Cen star formation timeline reveals four periods of enhanced star formation activity, or bursts, remarkably separated by about 5 Myr. Of these, the second burst, which occurred about 15 million years ago, is by far the dominant, and most of Sco-Cen's stars and clusters were in place by the end of this burst. The formation of stars and clusters in Sco-Cen is correlated but not linearly, implying that more stars were formed per cluster during the peak of the star formation rate. Most of the clusters that are large enough to have supernova precursors were formed during the 15 Myr period. Star and cluster formation activity has been continuously declining since then. We have clear evidence that Sco-Cen formed from the inside out and contains 100-pc long chains of contiguous clusters exhibiting well-defined age gradients, from massive older clusters to smaller young clusters. These observables suggest an important role for feedback in forming about half of Sco-Cen stars, although follow-up work is needed to quantify this statement. Finally, we confirm that the Upper-Sco age controversy discussed in the literature during the last decades is solved: the nine clusters previously lumped together as Upper-Sco, a benchmark region for planet formation studies, exhibit a wide range of ages from 3 to 19 Myr. △ Less

Submitted 14 May, 2023; v1 submitted 15 February, 2023; originally announced February 2023.

Comments: 20 pages, 15 figures, 3 tables, submitted to A&A

Journal ref: A&A 678, A71 (2023)

GRB 200829A: External Shock Origin of the Very Early Prompt Emission?

Authors: Jing Li, Da-Bin Lin, Rui-Jing Lu, Lu-Yao Jiang, Wen-Qiang Liang, Zhi-Lin Chen, Xiao-Yan Li, Xiang-Gao Wang, En-Wei Liang

Abstract: Long-duration GRB~200829A was detected by Fermi-GBM and Swift-BAT/XRT, and then rapidly observed by other ground-based telescopes. It has a weak $γ$-ray emission in the very early phase and followed by a bright spiky $γ$-ray emission pulse. The radiation spectrum of the very early emission is best fitted by a power-law function with index $\sim -1.7$. However, the bright spiky $γ$-ray pulse, espec… ▽ More Long-duration GRB~200829A was detected by Fermi-GBM and Swift-BAT/XRT, and then rapidly observed by other ground-based telescopes. It has a weak $γ$-ray emission in the very early phase and followed by a bright spiky $γ$-ray emission pulse. The radiation spectrum of the very early emission is best fitted by a power-law function with index $\sim -1.7$. However, the bright spiky $γ$-ray pulse, especially the time around the peak, exhibits a distinct two-component radiation spectra, i.e., Band function combined with a blackbody radiation spectrum. We infer the photospheric properties and reveal a medium magnetization at photospheric position by adopting the initial size of the outflow as $r_0=10^9$~cm. It implies that Band component in this pulse may be formed during the dissipation of magnetic field. The power-law radiation spectra found in the very early prompt emission may imply the external-shock origination of this phase. Then, we perform Markov Chain Monte Carlo method fitting on the light-curves of this burst, where the jet corresponding to the $γ$-ray pulses at around $20$~s is used to refresh the external-shock. It is shown that the light-curves of very early phase and X-ray afterglow after $40$~s, involving the X-ray bump at around $100$~s, can be well modelled in the external-shock scenario. For the obtained initial outflow, we estimate the minimum magnetization factor of the jet based on the fact that the photospheric emission of this jet is missed in the very early phase. △ Less

Submitted 2 January, 2023; originally announced January 2023.

Comments: Accepted for publication in Astrophysical Journal

Jets in a Gamma-Ray Burst During its Prompt Emission: Evolution of Lorentz Factor

Authors: Jing Li, Da-Bin Lin, Rui-Jing Lu, Yun Wang, Lu-Yao Jiang, Shen-Shi Du, Wen-Qiang Liang, Xiang-Gao Wang, En-Wei Liang

Abstract: Knowledge about the Lorentz factor and its evolution of relativistic jets in gamma-ray bursts (GRBs) is crucial to understand their physics. An exact value of bulk Lorentz factor can be estimated based on a high-energy spectral cutoff, which may appear in GRBs' prompt emission owing to the absorption of photon-photon pair production. In this work, we focus on the investigation of the bulk Lorentz… ▽ More Knowledge about the Lorentz factor and its evolution of relativistic jets in gamma-ray bursts (GRBs) is crucial to understand their physics. An exact value of bulk Lorentz factor can be estimated based on a high-energy spectral cutoff, which may appear in GRBs' prompt emission owing to the absorption of photon-photon pair production. In this work, we focus on the investigation of the bulk Lorentz factor evolution of jets in an individual burst. Based on \textsl{Fermi} observations, we search for the bursts with multiple $γ$-ray pulses characterized by a high-energy spectral cutoff, and nine GRBs are obtained. Together with the estimation of the pulse duration and radiation spectrum, the Lorentz factor of jets corresponding to different pulses in an individual GRB are estimated. It is shown that the Lorentz factor of jets in an individual GRB fluctuates within a certain range and without a general trend in these nine GRBs. In addition, the Lorentz factors of the jets in GRBs~130821A, 160509A and 160625B seem to increase with time. We also study the relations among $L_{\rm iso }$, $E_{\rm p,z}$, and $Γ$ for the pulses in our sample, which is found to be consistent with that found in previous works. △ Less

Submitted 26 December, 2022; originally announced December 2022.

Comments: Accepted for publication in ApJ

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