Showing 1–50 of 97 results for author: Lin, M

Rossby Wave Instability and Substructure Formation in 3D Non-Ideal MHD Wind-Launching Disks

Authors: Chun-Yen Hsu, Zhi-Yun Li, Yisheng Tu, Xiao Hu, Min-Kai Lin

Abstract: Rings and gaps are routinely observed in the dust continuum emission of protoplanetary discs (PPDs). How they form and evolve remains debated. Previous studies have demonstrated the possibility of spontaneous gas rings and gaps formation in wind-launching disks. Here, we show that such gas substructures are unstable to the Rossby Wave Instability (RWI) through numerical simulations. Specifically,… ▽ More Rings and gaps are routinely observed in the dust continuum emission of protoplanetary discs (PPDs). How they form and evolve remains debated. Previous studies have demonstrated the possibility of spontaneous gas rings and gaps formation in wind-launching disks. Here, we show that such gas substructures are unstable to the Rossby Wave Instability (RWI) through numerical simulations. Specifically, shorter wavelength azimuthal modes develop earlier, and longer wavelength ones dominate later, forming elongated (arc-like) anti-cyclonic vortices in the rings and (strongly magnetized) cyclonic vortices in the gaps that persist until the end of the simulation. Highly elongated vortices with aspect ratios of 10 or more are found to decay with time in our non-ideal MHD simulation, in contrast with the hydro case. This difference could be caused by magnetically induced motions, particularly strong meridional circulations with large values of the azimuthal component of the vorticity, which may be incompatible with the columnar structure preferred by vortices. The cyclonic and anti-cyclonic RWI vortices saturate at moderate levels, modifying but not destroying the rings and gaps in the radial gas distribution of the disk. In particular, they do not shut off the poloidal magnetic flux accumulation in low-density regions and the characteristic meridional flow patterns that are crucial to the ring and gap formation in wind-launching disks. Nevertheless, the RWI and their associated vortices open up the possibility of producing non-axisymmetric dust features observed in a small fraction of protoplanetary disks through non-ideal MHD, although detailed dust treatment is needed to explore this possibility. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Rossby wave instability in weakly ionized protoplanetary disks. I. azimuthal or vertical B-fields

Authors: Can Cui, Ashutosh Tripathi, Cong Yu, Min-Kai Lin, Andrew Youdin

Abstract: Rossby wave instability (RWI) is considered the underlying mechanism to crescent-shaped azimuthal asymmetries, discovered in (sub-)millimeter dust continuum of many protoplanetary disks. Previous works on linear theory were conducted in the hydrodynamic limit. Nevertheless, protoplanetary disks are likely magnetized and weakly ionized. We examine the influence of magnetic fields and non-ideal magn… ▽ More Rossby wave instability (RWI) is considered the underlying mechanism to crescent-shaped azimuthal asymmetries, discovered in (sub-)millimeter dust continuum of many protoplanetary disks. Previous works on linear theory were conducted in the hydrodynamic limit. Nevertheless, protoplanetary disks are likely magnetized and weakly ionized. We examine the influence of magnetic fields and non-ideal magnetohydrodynamic (MHD) effects - namely, Ohmic resistivity, Hall drift, and ambipolar diffusion - on the RWI unstable modes. We perform radially global linear analyses, employing constant azimuthal ($B_φ$) or vertical ($B_z$) background magnetic fields. It is found that, in the ideal MHD regime, magnetism can either enhance or diminish RWI growth. Strong non-ideal MHD effects cause RWI growth rates to recover hydrodynamic results. The sign of Hall Elsässer number subtly complicates the results, and vertical wavenumbers generically diminish growth rates. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 13 pages, 4 figures, submitted to MNRAS

Shoulder of Dust Rings Formed by Planet-disk Interactions

Authors: Jiaqing Bi, Min-Kai Lin

Abstract: Recent analyses of mm-wavelength protoplanetary disk observations have revealed several emission excesses on the previously identified dust rings, referred to as dust shoulders. The prevalence of dust shoulders suggests that they trace a common but unclear mechanism. In this work, we combine 3D, multifluid hydrodynamic simulations with radiative transfer calculations to explain the formation of du… ▽ More Recent analyses of mm-wavelength protoplanetary disk observations have revealed several emission excesses on the previously identified dust rings, referred to as dust shoulders. The prevalence of dust shoulders suggests that they trace a common but unclear mechanism. In this work, we combine 3D, multifluid hydrodynamic simulations with radiative transfer calculations to explain the formation of dust shoulders. We find that the ring-shoulder pairs can result from the 3D planet-disk interactions with massive, gap-opening planets. The key driver is the dust filtration effect at the local pressure maximum due to planet-driven outward gas flows. Our work provides a possible explanation for the outer dust shoulders in recent super-resolution analyses of ALMA observations. It also provides insights into the formation of the inner dust shoulder in the PDS 70 disk and highlights the role of 3D effects in planet-disk interaction studies. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: accepted to ApJ

Polar alignment of a dusty circumbinary disc -- I. Dust ring formation

Authors: Jeremy L. Smallwood, Min-Kai Lin, Hossam Aly, Rebecca Nealon, Cristiano Longarini

Abstract: We investigate the formation of dust traffic jams in polar-aligning circumbinary discs. We use 3D smoothed particle hydrodynamical simulations of both gas and dust to model an initially highly misaligned circumbinary disc around an eccentric binary. As the circumbinary disc evolves to a polar configuration (perpendicular to the binary orbital plane), the difference in the precession between the ga… ▽ More We investigate the formation of dust traffic jams in polar-aligning circumbinary discs. We use 3D smoothed particle hydrodynamical simulations of both gas and dust to model an initially highly misaligned circumbinary disc around an eccentric binary. As the circumbinary disc evolves to a polar configuration (perpendicular to the binary orbital plane), the difference in the precession between the gas and dust produces dust traffic jams, which become dense dust rings. We find the formation of dust rings exists for different Stokes number, binary eccentricity, and initial disc tilt. Dust rings are only produced while the circumbinary disc is misaligned to the binary orbital plane. When the disc becomes polar aligned, the dust rings are still present and long-lived. Once these dust rings are formed, they drift inward. The drift timescale depends on the Stokes number. The lower the Stokes number, the faster the dust ring drifts near the inner edge of the disc. The dust rings will have an increased midplane dust-to-go ratio, which may be a favourable environment for the steaming instability to operate. △ Less

Submitted 11 June, 2024; originally announced June 2024.

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

Convective overstability in radially global protoplanetary disks I -- Pure gas dynamics

Authors: Marius Lehmann, Min-Kai Lin

Abstract: Protoplanetary disks are prone to several hydrodynamic instabilities. One candidate, Convective Overstability (COS), can drive radial semi-convection that may influence dust dynamics and planetesimal formation. However, the COS has primarily been studied in local models. This paper investigates the COS near the mid-plane of radially global disk models. We first conduct a global linear stability an… ▽ More Protoplanetary disks are prone to several hydrodynamic instabilities. One candidate, Convective Overstability (COS), can drive radial semi-convection that may influence dust dynamics and planetesimal formation. However, the COS has primarily been studied in local models. This paper investigates the COS near the mid-plane of radially global disk models. We first conduct a global linear stability analysis, which shows that linear COS modes exist only radially inward of their Lindblad resonance (LR). The fastest-growing modes have LRs near the inner radial domain boundary with effective radial wavelengths that can be a substantial fraction of the disk radius. We then perform axisymmetric global simulations and find that the COS's nonlinear saturation is similar to previous incompressible shearing box simulations. In particular, we observe the onset of persistent zonal and elevator flows for sufficiently steep radial entropy gradients. In full 3D, non-axisymmetric global simulations, we find the COS produces large-scale, long-lived vortices, which induce outward radial transport of angular momentum via the excitation of spiral density waves. The corresponding $α$-viscosity values of order $10^{-3}$ agree well with those found in previous 3D compressible shearing box simulations. However, in global disks, significant modifications to their radial structure are found, including the formation of pressure bumps. Interestingly, the COS typically generates an outward radial mass transport, i.e. decretion. We briefly discuss the possible implications of our results for planetesimal formation and for interpreting dust rings and asymmetries observed in protoplanetary disks. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 30 pages, 24 figures

Analytical computation of bifurcation of orbits near collinear libration point in the restricted three-body problem

Authors: Mingpei Lin, Tong Luo, Hayato Chiba

Abstract: A unified analytical solution is presented for constructing the phase space near collinear libration points in the Circular Restricted Three-body Problem (CRTBP), encompassing Lissajous orbits and quasihalo orbits, their invariant manifolds, as well as transit and non-transit orbits. Traditional methods could only derive separate analytical solutions for the invariant manifolds of Lissajous orbits… ▽ More A unified analytical solution is presented for constructing the phase space near collinear libration points in the Circular Restricted Three-body Problem (CRTBP), encompassing Lissajous orbits and quasihalo orbits, their invariant manifolds, as well as transit and non-transit orbits. Traditional methods could only derive separate analytical solutions for the invariant manifolds of Lissajous orbits and halo orbits, falling short for the invariant manifolds of quasihalo orbits. By introducing a coupling coefficient η and a bifurcation equation, a unified series solution for these orbits is systematically developed using a coupling-induced bifurcation mechanism and Lindstedt-Poincaré method. Analyzing the third-order bifurcation equation reveals bifurcation conditions for halo orbits, quasihalo orbits, and their invariant manifolds. Furthermore, new families of periodic orbits similar to halo orbits are discovered, and non-periodic/quasi-periodic orbits, such as transit orbits and non-transit orbits, are found to undergo bifurcations. When η = 0, the series solution describes Lissajous orbits and their invariant manifolds, transit, and non-transit orbits. As η varies from zero to non-zero values, the solution seamlessly transitions to describe quasihalo orbits and their invariant manifolds, as well as newly bifurcated transit and non-transit orbits. This unified analytical framework provides a more comprehensive understanding of the complex phase space structures near collinear libration points in the CRTBP. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: Our unified analytical framework provides a holistic view of the phase space structures near collinear libration points in the CRTBP. It successfully addresses analytical challenges related to quasihalo orbits and their invariant manifolds. The proposed coupling-induced bifurcation mechanism and the corresponding coupling coefficients can also be applied to general dynamical systems

Bifurcation mechanism of quasihalo orbit from Lissajous orbit in the restricted three-body problem

Authors: Mingpei Lin, Hayato Chiba

Abstract: This paper presents a general analytical method to describe the center manifolds of collinear libration points in the Restricted Three-body Problem (RTBP). It is well-known that these center manifolds include Lissajous orbits, halo orbits, and quasihalo orbits. Previous studies have traditionally treated these orbits separately by iteratively constructing high-order series solutions using the Lind… ▽ More This paper presents a general analytical method to describe the center manifolds of collinear libration points in the Restricted Three-body Problem (RTBP). It is well-known that these center manifolds include Lissajous orbits, halo orbits, and quasihalo orbits. Previous studies have traditionally treated these orbits separately by iteratively constructing high-order series solutions using the Lindstedt-Poincaré method. Instead of relying on resonance between their frequencies, this study identifies that halo and quasihalo orbits arise due to intricate coupling interactions between in-plane and out-of-plane motions. To characterize this coupling effect, a novel concept, coupling coefficient $η$, is introduced in the RTBP, incorporating the coupling term $ηΔx$ in the $z$-direction dynamics equation, where $Δ$ represents a formal power series concerning the amplitudes. Subsequently, a uniform series solution for these orbits is constructed up to a specified order using the Lindstedt-Poincaré method. For any given paired in-plane and out-of-plane amplitudes, the coupling coefficient $η$ is determined by the bifurcation equation $Δ= 0$. When $η$ = 0, the proposed solution describes Lissajous orbits around libration points. As $η$ transitions from zero to non-zero values, the solution describes quasihalo orbits, which bifurcate from Lissajous orbits. Particularly, halo orbits bifurcate from planar Lyapunov orbits if the out-of-plane amplitude is zero. The proposed method provides a unified framework for understanding these intricate orbital behaviors in the RTBP. △ Less

Submitted 11 July, 2024; v1 submitted 26 January, 2024; originally announced January 2024.

Dust Dynamics in Hall-effected Protoplanetary Disks. I. Background Drift Hall Instability

Authors: Yinhao Wu, Min-Kai Lin, Can Cui, Leonardo Krapp, Yueh-Ning Lee, Andrew N. Youdin

Abstract: Recent studies have shown that the large-scale gas dynamics of protoplanetary disks (PPDs) are controlled by non-ideal magneto-hydrodynamics (MHD), but how this influences dust dynamics is not fully understood. To this end, we investigate the stability of dusty, magnetized disks subject to the Hall effect, which applies to planet-forming regions of PPDs. We find a novel Background Drift Hall Insta… ▽ More Recent studies have shown that the large-scale gas dynamics of protoplanetary disks (PPDs) are controlled by non-ideal magneto-hydrodynamics (MHD), but how this influences dust dynamics is not fully understood. To this end, we investigate the stability of dusty, magnetized disks subject to the Hall effect, which applies to planet-forming regions of PPDs. We find a novel Background Drift Hall Instability (BDHI) that may facilitate planetesimal formation in Hall-effected disk regions. Through a combination of linear analysis and nonlinear simulations, we demonstrate the viability and characteristics of BDHI. We find it can potentially dominate over the classical streaming instability (SI) and standard MHD instabilities at low dust-to-gas ratios and weak magnetic fields. We also identify magnetized versions of the classic SI, but these are usually subdominant. We highlight the complex interplay between magnetic fields and dust-gas dynamics in PPDs, underscoring the need to consider non-ideal MHD like the Hall effect in the broader narrative of planet formation. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: 23 pages, 12 figures, accepted by ApJ. Welcome any comments and suggestions!

Filament formation due to diffusive instabilities in dusty protoplanetary disks

Authors: Konstantin Gerbig, Min-Kai Lin, Marius Lehmann

Abstract: We report the finding of a new, local diffusion instability in a protoplanetary disk, which can operate in a dust fluid, subject to mass diffusion, shear viscosity, and dust-gas drag, provided diffusivity, viscosity, or both decrease sufficiently rapidly with increasing dust surface mass density. We devise a vertically averaged, axisymmetric hydrodynamic model to describe a dense, mid-plane dust l… ▽ More We report the finding of a new, local diffusion instability in a protoplanetary disk, which can operate in a dust fluid, subject to mass diffusion, shear viscosity, and dust-gas drag, provided diffusivity, viscosity, or both decrease sufficiently rapidly with increasing dust surface mass density. We devise a vertically averaged, axisymmetric hydrodynamic model to describe a dense, mid-plane dust layer in a protoplanetary disk. The gas is modeled as a passive component, imposing an effective, diffusion-dependent pressure, mass diffusivity, and viscosity onto the otherwise collisionless dust fluid, via turbulence excited by the gas alone, or dust and gas in combination. In particular, we argue that such conditions are met when the dust-gas mixture generates small-scale turbulence through the streaming instability, as supported by recent measurements of dust mass diffusion slopes in simulations. We hypothesize that the newly discovered instability may be the origin of filamentary features, almost ubiquitously found in simulations of the streaming instability. In addition, our model allows for growing oscillatory modes, which operate in a similar fashion as the axisymmetric viscous overstability in dense planetary rings. However, it remains speculative if the required conditions for such modes can be met in protoplanetary disks. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: Accepted for publication in ApJ, 23 pages, 13 figures

Cosmological shocks around galaxy clusters: A coherent investigation with DES, SPT & ACT

Authors: D. Anbajagane, C. Chang, E. J. Baxter, S. Charney, M. Lokken, M. Aguena, S. Allam, O. Alves, A. Amon, R. An, F. Andrade-Oliveira, D. Bacon, N. Battaglia, K. Bechtol, M. R. Becker, B. A. Benson, G. M. Bernstein, L. Bleem, S. Bocquet, J. R. Bond, D. Brooks, A. Carnero Rosell, M. Carrasco Kind, R. Chen, A. Choi , et al. (89 additional authors not shown)

Abstract: We search for signatures of cosmological shocks in gas pressure profiles of galaxy clusters using the cluster catalogs from three surveys: the Dark Energy Survey (DES) Year 3, the South Pole Telescope (SPT) SZ survey, and the Atacama Cosmology Telescope (ACT) data releases 4, 5, and 6, and using thermal Sunyaev-Zeldovich (SZ) maps from SPT and ACT. The combined cluster sample contains around… ▽ More We search for signatures of cosmological shocks in gas pressure profiles of galaxy clusters using the cluster catalogs from three surveys: the Dark Energy Survey (DES) Year 3, the South Pole Telescope (SPT) SZ survey, and the Atacama Cosmology Telescope (ACT) data releases 4, 5, and 6, and using thermal Sunyaev-Zeldovich (SZ) maps from SPT and ACT. The combined cluster sample contains around $10^5$ clusters with mass and redshift ranges $10^{13.7} < M_{\rm 200m}/M_\odot < 10^{15.5}$ and $0.1 < z < 2$, and the total sky coverage of the maps is $\approx 15,000 \,\,{\rm deg}^2$. We find a clear pressure deficit at $R/R_{\rm 200m}\approx 1.1$ in SZ profiles around both ACT and SPT clusters, estimated at $6σ$ significance, which is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions. The feature is not as clearly determined in profiles around DES clusters. We verify that measurements using SPT or ACT maps are consistent across all scales, including in the deficit feature. The SZ profiles of optically selected and SZ-selected clusters are also consistent for higher mass clusters. Those of less massive, optically selected clusters are suppressed on small scales by factors of 2-5 compared to predictions, and we discuss possible interpretations of this behavior. An oriented stacking of clusters -- where the orientation is inferred from the SZ image, the brightest cluster galaxy, or the surrounding large-scale structure measured using galaxy catalogs -- shows the normalization of the one-halo and two-halo terms vary with orientation. Finally, the location of the pressure deficit feature is statistically consistent with existing estimates of the splashback radius. △ Less

Submitted 12 December, 2023; v1 submitted 29 September, 2023; originally announced October 2023.

Comments: [v2]: Version accepted to MNRAS

Probing Parity Violation in the Stochastic Gravitational Wave Background with Astrometry

Authors: Qiuyue Liang, Meng-Xiang Lin, Mark Trodden, Sam S. C. Wong

Abstract: Astrometry holds the potential for testing fundamental physics through the effects of the Stochastic Gravitational Wave Background (SGWB) in the $\sim 1-100$ nHz frequency band on precision measurements of stellar positions. Such measurements are complementary to tests made possible by the detection of the SGWB using Pulsar Timing Arrays. Here, the feasibility of using astrometry for the identific… ▽ More Astrometry holds the potential for testing fundamental physics through the effects of the Stochastic Gravitational Wave Background (SGWB) in the $\sim 1-100$ nHz frequency band on precision measurements of stellar positions. Such measurements are complementary to tests made possible by the detection of the SGWB using Pulsar Timing Arrays. Here, the feasibility of using astrometry for the identification of parity-violating signals within the SGWB is investigated. This is achieved by defining and quantifying a non-vanishing $EB$ correlation function within astrometric correlation functions, and investigating how one might estimate the detectability of such signals. △ Less

Submitted 6 May, 2024; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: Match the published version, adding references. 8 pages, 2 figures

Late Time Modification of Structure Growth and the S8 Tension

Authors: Meng-Xiang Lin, Bhuvnesh Jain, Marco Raveri, Eric J. Baxter, Chihway Chang, Marco Gatti, Sujeong Lee, Jessica Muir

Abstract: The $S_8$ tension between low-redshift galaxy surveys and the primary CMB signals a possible breakdown of the $Λ$CDM model. Recently differing results have been obtained using low-redshift galaxy surveys and the higher redshifts probed by CMB lensing, motivating a possible time-dependent modification to the growth of structure. We investigate a simple phenomenological model in which the growth of… ▽ More The $S_8$ tension between low-redshift galaxy surveys and the primary CMB signals a possible breakdown of the $Λ$CDM model. Recently differing results have been obtained using low-redshift galaxy surveys and the higher redshifts probed by CMB lensing, motivating a possible time-dependent modification to the growth of structure. We investigate a simple phenomenological model in which the growth of structure deviates from the $Λ$CDM prediction at late times, in particular as a simple function of the dark energy density. Fitting to galaxy lensing, CMB lensing, BAO, and Supernovae datasets, we find significant evidence - 2.5 - 3$σ$, depending on analysis choices - for a non-zero value of the parameter quantifying a deviation from $Λ$CDM. The preferred model, which has a slower growth of structure below $z\sim 1$, improves the joint fit to the data over $Λ$CDM. While the overall fit is improved, there is weak evidence for galaxy and CMB lensing favoring different changes in the growth of structure. △ Less

Submitted 15 March, 2024; v1 submitted 30 August, 2023; originally announced August 2023.

Comments: Match the published version. Add more references. 6 pages, 5 figures, 1 table

Journal ref: Phys.Rev.D 109 (2024) 6, 063523

Distinguishing Magnetized Disc Winds from Turbulent Viscosity through Substructure Morphology in Planet-forming Discs

Authors: Yinhao Wu, Yi-Xian Chen, Haochang Jiang, Ruobing Dong, Enrique Macías, Min-Kai Lin, Giovanni P. Rosotti, Vardan Elbakyan

Abstract: The traditional paradigm of viscosity-dominated evolution of protoplanetary discs has been recently challenged by magnetized disc winds. However, distinguishing wind-driven and turbulence-driven accretion through observations has been difficult. In this study, we present a novel approach to identifying their separate contribution to angular momentum transport by studying the gap and ring morpholog… ▽ More The traditional paradigm of viscosity-dominated evolution of protoplanetary discs has been recently challenged by magnetized disc winds. However, distinguishing wind-driven and turbulence-driven accretion through observations has been difficult. In this study, we present a novel approach to identifying their separate contribution to angular momentum transport by studying the gap and ring morphology of planet-forming discs in the ALMA continuum. We model the gap-opening process of planets in discs with both viscous evolution and wind-driven accretion by 2D multi-fluid hydrodynamical simulations. Our results show that gap-opening planets in wind-driven accreting discs generate characteristic substructures that differ from those in purely viscous discs. Specifically, we demonstrate that discs, where wind-driven accretion dominates the production of substructures, exhibit significant asymmetries. Based on the diverse outputs of mock images in the ALMA continuum, we roughly divide the planet-induced features into four regimes (moderate-viscosity dominated, moderate-wind dominated, strong-wind dominated, inviscid). The classification of these regimes sets up a potential method to constrain the strength of magnetized disc wind and viscosity based on the observed gap and ring morphology. We discuss the asymmetry feature in our mock images and its potential manifestation in ALMA observations. △ Less

Submitted 27 May, 2023; v1 submitted 14 May, 2023; originally announced May 2023.

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

Instabilities in dusty non-isothermal proto-planetary discs

Authors: Marius Lehmann, Min-Kai Lin

Abstract: Protoplanetary discs (PPDs) can host a number of instabilities that may partake directly or indirectly in the process of planetesimal formation. These include the Vertical Shear Instability (VSI), Convective Overstability (COS), Streaming Instability (SI), and Dust Settling Instability (DSI), to name a few. Notably, the VSI and COS have mostly been studied in purely gaseous discs, while the SI and… ▽ More Protoplanetary discs (PPDs) can host a number of instabilities that may partake directly or indirectly in the process of planetesimal formation. These include the Vertical Shear Instability (VSI), Convective Overstability (COS), Streaming Instability (SI), and Dust Settling Instability (DSI), to name a few. Notably, the VSI and COS have mostly been studied in purely gaseous discs, while the SI and DSI have only been analyzed in isothermal discs. How these instabilities operate under more general conditions is therefore unclear. To this end, we devise a local model of a PPD describing a non-isothermal gas interacting with a single species of dust via drag forces. Using this, we find that dust drag sets minimum length scales below which the VSI and COS are suppressed. Similarly, we find that the SI can be suppressed on sufficiently small scales by the gas' radial buoyancy if it cools on roughly a dynamical timescale. We show that the DSI can be effectively stabilized by vertical buoyancy, except at special radial and vertical length scales. We also find novel instabilities unique to a dusty, non-isothermal gas. These result in a dusty analog of the COS that operates in slowly cooled discs, and a dusty version of the VSI that is strongly enhanced by dust settling. We briefly discuss the possible implications of our results on planetesimal formation. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Comments: 40 pages, 26 figures

A Test of Gravity with Pulsar Timing Arrays

Authors: Qiuyue Liang, Meng-Xiang Lin, Mark Trodden

Abstract: A successful measurement of the Stochastic Gravitational Wave Background (SGWB) in Pulsar Timing Arrays (PTAs) would open up a new window through which to test the predictions of General Relativity (GR). We consider how these measurements might reveal deviations from GR by studying the overlap reduction function -- the quantity that in GR is approximated by the Hellings-Downs curve -- in some samp… ▽ More A successful measurement of the Stochastic Gravitational Wave Background (SGWB) in Pulsar Timing Arrays (PTAs) would open up a new window through which to test the predictions of General Relativity (GR). We consider how these measurements might reveal deviations from GR by studying the overlap reduction function -- the quantity that in GR is approximated by the Hellings-Downs curve -- in some sample modifications of gravity, focusing on the generic prediction of a modified dispersion relation for gravitational waves. We find a distinct signature of such modifications to GR -- a shift in the minimum angle of the angular distribution -- and demonstrate that this shift is quantitatively sensitive to any change in the phase velocity. In a given modification of gravity, this result can be used, in some regions of parameter space, to distinguish the effect of a modified dispersion relation from that due to the presence of extra polarization modes. △ Less

Submitted 14 April, 2023; v1 submitted 5 April, 2023; originally announced April 2023.

Comments: v2: add references and discussions. 23 pages, 4 figures

How to form compact & other longer-lived planet-induced vortices: VSI, planet migration, or re-triggers, but not feedback

Authors: Michael Hammer, Min-Kai Lin

Abstract: Past computational studies of planet-induced vortices have shown that the dust asymmetries associated with these vortices can be long-lived enough that they should be much more common in mm/sub-mm observations of protoplanetary discs, even though they are quite rare. Observed asymmetries also have a range of azimuthal extents from compact to elongated even though computational studies have shown p… ▽ More Past computational studies of planet-induced vortices have shown that the dust asymmetries associated with these vortices can be long-lived enough that they should be much more common in mm/sub-mm observations of protoplanetary discs, even though they are quite rare. Observed asymmetries also have a range of azimuthal extents from compact to elongated even though computational studies have shown planet-induced vortices should be preferentially elongated. In this study, we use 2-D and 3-D hydrodynamic simulations to test whether those dust asymmetries should really be so long-lived or so elongated. With higher resolution (29 cells radially per scale height) than our previous work, we find that vortices can be more compact by developing compact cores when higher-mass planets cause them to re-form, or if they are seeded by tiny compact vortices from the vertical shear instability (VSI), but not through dust feedback in 3-D as was previously expected in general. Any case with a compact vortex or core(s) also has a longer lifetime. Even elongated vortices can have longer lifetimes with higher-mass planets or if the associated planet is allowed to migrate, the latter of which can cause the dust asymmetry to stop decaying as the planet migrates away from the vortex. These longer dust asymmetry lifetimes are even more inconsistent with observations, perhaps suggesting that discs still have an intermediate amount of effective viscosity. △ Less

Submitted 5 December, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Comments: 24 pages, 21 figures; Published in MNRAS

A Dark Matter Trigger for Early Dark Energy Coincidence

Authors: Meng-Xiang Lin, Evan McDonough, J. Colin Hill, Wayne Hu

Abstract: Early dark energy (EDE), whose cosmological role is localized in time around the epoch of matter-radiation equality in order to resolve the Hubble tension, introduces a new coincidence problem: why should the EDE dynamics occur near equality if EDE is decoupled from both matter and radiation? The resolution of this problem may lie in an {\it early dark sector} (EDS), wherein the dark matter mass i… ▽ More Early dark energy (EDE), whose cosmological role is localized in time around the epoch of matter-radiation equality in order to resolve the Hubble tension, introduces a new coincidence problem: why should the EDE dynamics occur near equality if EDE is decoupled from both matter and radiation? The resolution of this problem may lie in an {\it early dark sector} (EDS), wherein the dark matter mass is dependent on the EDE scalar field. Concretely, we consider a Planck-suppressed coupling of EDE to dark matter, as would naturally arise from breaking of the global $U(1)$ shift symmetry of the former by quantum gravity effects. With a sufficiently flat potential, the rise to dominance of dark matter at matter-radiation equality itself triggers the rolling and subsequent decay of the EDE. We show that this {\it trigger} EDS (tEDS) model can naturally resolve the EDE coincidence problem at the background level without any fine tuning of the coupling to dark matter or of the initial conditions. When fitting to current cosmological data, including that from the local distance ladder and the low-redshift amplitude of fluctuations, the tEDS maximum-likelihood model performs comparably to EDE for resolving the Hubble tension, achieving $H_0 =71.2$ km/s/Mpc. However, fitting the \emph{Planck} cosmic microwave background data requires a specific range of initial field positions to balance the scalar field fluctuations that drive acoustic oscillations, providing testable differences with other EDE models. △ Less

Submitted 15 December, 2022; originally announced December 2022.

Comments: 12 pages, 6 figures

Gap-opening Planets Make Dust Rings Wider

Authors: Jiaqing Bi, Min-Kai Lin, Ruobing Dong

Abstract: As one of the most commonly observed disk substructures, dust rings from high-resolution disk surveys appear to have different radial widths. Recent observations on PDS 70 and AB Aur reveal not only planets in the disk, but also the accompanying wide dust rings. We use three-dimensional dust-and-gas disk simulations to study whether gap-opening planets are responsible for the large ring width in d… ▽ More As one of the most commonly observed disk substructures, dust rings from high-resolution disk surveys appear to have different radial widths. Recent observations on PDS 70 and AB Aur reveal not only planets in the disk, but also the accompanying wide dust rings. We use three-dimensional dust-and-gas disk simulations to study whether gap-opening planets are responsible for the large ring width in disk observations. We find that gap-opening planets can widen rings of dust trapped at the pressure bump via planetary perturbations, even with the mid-plane dust-to-gas ratio approaching order unity and with the dust back-reaction accounted for. We show that the planet-related widening effect of dust rings can be quantified using diffusion-advection theory, and provide a generalized criterion for an equilibrated dust ring width in three-dimensional disk models. We also suggest that the ring width can be estimated using the gas turbulent viscosity $α_{\rm turb}$, but with cautions about the Schmidt number greater than order unity. △ Less

Submitted 20 October, 2022; originally announced October 2022.

Comments: 17 pages, 6+3 figures, 2 tables, accepted in ApJ

Nonlinear evolution of streaming instabilities in accreting protoplanetary disks

Authors: Chun-Yen Hsu, Min-Kai Lin

Abstract: The streaming instability (SI) is one of the most promising candidates for triggering planetesimal formation by producing dense dust clumps that undergo gravitational collapse. Understanding how the SI operates in realistic protoplanetary disks (PPDs) is therefore crucial to assess the efficiency of planetesimal formation. Modern models of PPDs show that large-scale magnetic torques or winds can d… ▽ More The streaming instability (SI) is one of the most promising candidates for triggering planetesimal formation by producing dense dust clumps that undergo gravitational collapse. Understanding how the SI operates in realistic protoplanetary disks (PPDs) is therefore crucial to assess the efficiency of planetesimal formation. Modern models of PPDs show that large-scale magnetic torques or winds can drive laminar gas accretion near the disk midplane. In a previous study, we identified a new linear dust-gas instability, the azimuthal drift SI (AdSI), applicable to such accreting disks and is powered by the relative azimuthal motion between dust and gas that results from the gas being torqued. In this work, we present the first nonlinear simulations of the AdSI. We show that it can destabilize an accreting, dusty disk even in the absence of a global radial pressure gradient, which is unlike the classic SI. We find the AdSI drives turbulence and the formation of vertically-extended dust filaments that undergo merging. In dust-rich disks, merged AdSI filaments reach maximum dust-to-gas ratios exceeding 100. Moreover, we find that even in dust-poor disks the AdSI can increase local dust densities by two orders of magnitude. We discuss the possible role of the AdSI in planetesimal formation, especially in regions of an accreting PPD with vanishing radial pressure gradients. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: Accepted by ApJ

Revealing the drag instability in one-fluid nonideal MHD simulations of a 1D isothermal C-shock

Authors: Pin-Gao Gu, Che-Yu Chen, Emma Shen, Chien-Chang Yen, Min-Kai Lin

Abstract: C-type shocks are believed to be ubiquitous in turbulent molecular clouds thanks to ambipolar diffusion. We investigate whether the drag instability in 1D isothermal C-shocks, inferred from the local linear theory of Gu & Chen, can appear in non-ideal magnetohydrodynamic simulations. Two C-shock models (with narrow and broad steady-state shock widths) are considered to represent the typical enviro… ▽ More C-type shocks are believed to be ubiquitous in turbulent molecular clouds thanks to ambipolar diffusion. We investigate whether the drag instability in 1D isothermal C-shocks, inferred from the local linear theory of Gu & Chen, can appear in non-ideal magnetohydrodynamic simulations. Two C-shock models (with narrow and broad steady-state shock widths) are considered to represent the typical environment of star-forming clouds. The ionization-recombination equilibrium is adopted for the one-fluid approach. In the 1D simulation, the inflow gas is continuously perturbed by a sinusoidal density fluctuation with a constant frequency. The perturbations clearly grow after entering the C-shock region until they start being damped at the transition to the postshock region. We show that the profiles of a predominant Fourier mode extracted locally from the simulated growing perturbation match those of the growing mode derived from the linear analysis. Moreover, the local growth rate and wave frequency derived from the predominant mode generally agree with those from the linear theory. Therefore, we confirm the presence of the drag instability in simulated 1D isothermal C-shocks. We also explore the nonlinear behavior of the instability by imposing larger-amplitude perturbations to the simulation. We find that the drag instability is subject to wave steepening, leading to saturated perturbation growth. Issues concerning local analysis, nonlinear effects, one-fluid approach, and astrophysical applications are discussed. △ Less

Submitted 19 August, 2022; v1 submitted 9 July, 2022; originally announced July 2022.

Comments: 19 pages, 10 figures

Journal ref: ApJ 2022, 935, 95

Modified gravitational wave propagation with higher modes and its degeneracies with lensing

Authors: Jose Maria Ezquiaga, Wayne Hu, Macarena Lagos, Meng-Xiang Lin, Fei Xu

Abstract: Low-energy alternatives to General Relativity (GR) generically modify the phase of gravitational waves (GWs) during their propagation. As detector sensitivities increase, it becomes key to understand how these modifications affect the GW higher modes and to disentangle possible degeneracies with astrophysical phenomena. We apply a general formalism -- the WKB approach -- for solving analytically w… ▽ More Low-energy alternatives to General Relativity (GR) generically modify the phase of gravitational waves (GWs) during their propagation. As detector sensitivities increase, it becomes key to understand how these modifications affect the GW higher modes and to disentangle possible degeneracies with astrophysical phenomena. We apply a general formalism -- the WKB approach -- for solving analytically wave propagation in the spatial domain with a modified dispersion relation (MDR). We compare this WKB approach to applying a stationary phase approximation (SPA) in the temporal domain with time delays associated to the group or particle velocity. To this end, we extend the SPA to generic signals with higher modes, keeping careful track of reference phases and arrival times. We find that the WKB approach coincides with the SPA using the group velocity, in agreement with the principles of wave propagation. We then explore the degeneracies between a GW propagation with an MDR and a strongly-lensed GW in GR, since the latter can introduce a frequency-independent phase shift which is not degenerate with source parameters in the presence of higher modes. We find that for a particular MDR there is an exact degeneracy for wave propagation, unlike with the SPA for particle propagation. For the other cases, we search for the values of the MDR parameters that minimize the $χ^2$ and conclude that strongly-lensed GR GWs could be misinterpreted as GWs in modified gravity. Future MDR constraints with higher mode GWs should include the possibility of frequency-independent phase shifts, allowing for the identification of modified gravity and strong lensing distortions at the same time. △ Less

Submitted 15 August, 2022; v1 submitted 24 March, 2022; originally announced March 2022.

Comments: 32 pages, 6 figures. Match the published version

Journal ref: Journal of Cosmology and Astroparticle Physics, Volume 2022, Issue 08, id.016, 33 pp

Hydro-, Magnetohydro-, and Dust-Gas Dynamics of Protoplanetary Disks

Authors: G. Lesur, B. Ercolano, M. Flock, M. -K. Lin, C. -C. Yang, J. A. Barranco, P. Benitez-Llambay, J. Goodman, A. Johansen, H. Klahr, G. Laibe, W. Lyra, P. Marcus, R. P. Nelson, J. Squire, J. B. Simon, N. Turner, O. M. Umurhan, A. N. Youdin

Abstract: The building of planetary systems is controlled by the gas and dust dynamics of protoplanetary disks. While the gas is simultaneously accreted onto the central star and dissipated away by winds, dust grains aggregate and collapse to form planetesimals and eventually planets. This dust and gas dynamics involves instabilities, turbulence and complex non-linear interactions which ultimately control t… ▽ More The building of planetary systems is controlled by the gas and dust dynamics of protoplanetary disks. While the gas is simultaneously accreted onto the central star and dissipated away by winds, dust grains aggregate and collapse to form planetesimals and eventually planets. This dust and gas dynamics involves instabilities, turbulence and complex non-linear interactions which ultimately control the observational appearance and the secular evolution of these disks. This chapter is dedicated to the most recent developments in our understanding of the dynamics of gaseous and dusty disks, covering hydrodynamic and magnetohydrodynamic turbulence, gas-dust instabilities, dust clumping and disk winds. We show how these physical processes have been tested from observations and highlight standing questions that should be addressed in the future. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: 32 pages, 7 figures, to appear in Protostars and Planets VII, eds: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

Using Bayesian Deep Learning to infer Planet Mass from Gaps in Protoplanetary Disks

Authors: Sayantan Auddy, Ramit Dey, Min-Kai Lin, Daniel Carrera, Jacob B. Simon

Abstract: Planet induced sub-structures, like annular gaps, observed in dust emission from protoplanetary disks provide a unique probe to characterize unseen young planets. While deep learning based model has an edge in characterizing the planet's properties over traditional methods, like customized simulations and empirical relations, it lacks in its ability to quantify the uncertainty associated with its… ▽ More Planet induced sub-structures, like annular gaps, observed in dust emission from protoplanetary disks provide a unique probe to characterize unseen young planets. While deep learning based model has an edge in characterizing the planet's properties over traditional methods, like customized simulations and empirical relations, it lacks in its ability to quantify the uncertainty associated with its predictions. In this paper, we introduce a Bayesian deep learning network "DPNNet-Bayesian" that can predict planet mass from disk gaps and provides uncertainties associated with the prediction. A unique feature of our approach is that it can distinguish between the uncertainty associated with the deep learning architecture and uncertainty inherent in the input data due to measurement noise. The model is trained on a data set generated from disk-planet simulations using the \textsc{fargo3d} hydrodynamics code with a newly implemented fixed grain size module and improved initial conditions. The Bayesian framework enables estimating a gauge/confidence interval over the validity of the prediction when applied to unknown observations. As a proof-of-concept, we apply DPNNet-Bayesian to dust gaps observed in HL Tau. The network predicts masses of $ 86.0 \pm 5.5 M_{\Earth} $, $ 43.8 \pm 3.3 M_{\Earth} $, and $ 92.2 \pm 5.1 M_{\Earth} $ respectively, which are comparable to other studies based on specialized simulations. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: 14 pages, 6 figures, submitted to ApJ

The Early Dark Sector, the Hubble Tension, and the Swampland

Authors: Evan McDonough, Meng-Xiang Lin, J. Colin Hill, Wayne Hu, Shengjia Zhou

Abstract: We consider the interplay of the Early Dark Energy (EDE) model, the Swampland Distance Conjecture (SDC), and cosmological parameter tensions. EDE is a proposed resolution of the Hubble tension relying upon a near-Planckian scalar field excursion, while the SDC predicts an exponential sensitivity of masses of other fields to such an excursion, $m\propto e^{-c|Δφ|/M_{\rm pl}}$ with… ▽ More We consider the interplay of the Early Dark Energy (EDE) model, the Swampland Distance Conjecture (SDC), and cosmological parameter tensions. EDE is a proposed resolution of the Hubble tension relying upon a near-Planckian scalar field excursion, while the SDC predicts an exponential sensitivity of masses of other fields to such an excursion, $m\propto e^{-c|Δφ|/M_{\rm pl}}$ with $c\sim{\cal O}(1)$. Meanwhile, EDE is in tension with large-scale structure (LSS) data, due to shifts in the standard $Λ$CDM parameters necessary to fit the cosmic microwave background (CMB). One might hope that a proper treatment of the model, e.g., accounting for the SDC, may ameliorate the tension with LSS. Motivated by these considerations, we introduce the Early Dark Sector (EDS) model, wherein the mass of dark matter is exponentially sensitive to super-Planckian field excursions of the EDE scalar. The EDS model exhibits new phenomenology in both the early and late universe, the latter due to an EDE-mediated dark matter self-interaction. This dark matter-philic "fifth force", while constrained to be small, remains active in the late universe and is not screened in virialized halos. We find that the new interaction with dark matter partially resolves the LSS tension. However, the marginalized posteriors are nonetheless consistent with $f_{\rm EDE}=0$ at 95$\%$ CL once the Dark Energy Survey Year 3 measurement of $S_8$ is included. We study constraints on the model from Atacama Cosmology Telescope data, and find a factor of two improvement on the error bar on the SDC parameter $c$, along with an increased preference for the EDE component. We discuss the implications of these constraints for the SDC, and find the tightest observational constraints to date on a swampland parameter, suggesting that an EDE description of cosmological data is in tension with the SDC. △ Less

Submitted 26 August, 2022; v1 submitted 16 December, 2021; originally announced December 2021.

Comments: v2: version published in Phys. Rev. D

Journal ref: Phys. Rev. D 106, 043525 2022

Impact of Local Pressure Enhancements on Dust Concentration inTurbulent Protoplanetary Discs

Authors: Marius Lehmann, Min-Kai Lin

Abstract: We investigate the evolution of dust and gas in the vicinity of local pressure enhancements (pressure bumps) in a protoplanetary disc (PPD) with turbulence due to the Vertical Shear Instability (VSI). We perform global 2D axisymmetric and 3D simulations of dust and gas for a range of values for Z (ratio of dust-to-gas (d/g) surface mass densities or metallicity), particle Stokes numbers tau, and p… ▽ More We investigate the evolution of dust and gas in the vicinity of local pressure enhancements (pressure bumps) in a protoplanetary disc (PPD) with turbulence due to the Vertical Shear Instability (VSI). We perform global 2D axisymmetric and 3D simulations of dust and gas for a range of values for Z (ratio of dust-to-gas (d/g) surface mass densities or metallicity), particle Stokes numbers tau, and pressure bump amplitude A. Dust feedback onto the gas is included. For the first time we demonstrate in global 3D simulations the collection of dust in long-lived vortices induced by the VSI. Without a pressure bump and for Z~0.01 and tau~0.01 we find that such vortices reach d/g density ratios slightly below unity in the discs mid-plane, while for Z>0.05 long-lived vortices are largely absent. In presence of a pressure bump, for Z~0.01 & tau~0.01 a dusty vortex forms reaching d/g ratios of a few times unity, such that the SI is expected to develop, before it eventually shears out into a turbulent dust ring. For Z~0.03 this occurs for tau~0.005, with a weaker, more short-lived vortex, while for larger tau only a turbulent dust ring forms. For Z>0.03 we find that the dust ring becomes increasingly axisymmetric for increasing tau and d/g ratios reach ~1 for tau>0.005. Furthermore, the disc's vertical mass flow profile is strongly affected by dust for Z>0.03, such that gas is transported inward near the mid-plane and outward at larger heights, i.e. the reversed situation compared to simulations with zero or small amounts of dust. Viscous $α$-values drop moderately as 0.001-0.0001 for increasing Z=0-0.05. Our results suggest that the VSI can play an active role in planetesimal formation through the formation of vortices for plausible values of Z and tau. Also it may provide a natural explanation for the presence/absence of asymmetries of observed dust rings in PPDs, depending on the value of Z. △ Less

Submitted 12 December, 2021; originally announced December 2021.

Journal ref: A&A 658, A156 (2022)

Streaming instabilities in accreting and magnetized laminar protoplanetary disks

Authors: Min-Kai Lin, Chun-Yen Hsu

Abstract: The streaming instability is one of the most promising pathways to the formation of planetesimals from pebbles. Understanding how this instability operates under realistic conditions expected in protoplanetary disks is therefore crucial to assess the efficiency of planet formation. Contemporary models of protoplanetary disks show that magnetic fields are key to driving gas accretion through large-… ▽ More The streaming instability is one of the most promising pathways to the formation of planetesimals from pebbles. Understanding how this instability operates under realistic conditions expected in protoplanetary disks is therefore crucial to assess the efficiency of planet formation. Contemporary models of protoplanetary disks show that magnetic fields are key to driving gas accretion through large-scale, laminar magnetic stresses. However, the effect of such magnetic fields on the streaming instability has not been examined in detail. To this end, we study the stability of dusty, magnetized gas in a protoplanetary disk. We find the streaming instability can be enhanced by passive magnetic torques and even persist in the absence of a global radial pressure gradient. In this case, instability is attributed to the azimuthal drift between dust and gas, unlike the classical streaming instability, which is driven by radial drift. This suggests that the streaming instability can remain effective inside dust-trapping pressure bumps in accreting disks. When a live vertical field is considered, we find the magneto-rotational instability can be damped by dust feedback, while the classic streaming instability can be stabilized by magnetic perturbations. We also find that Alfvén waves can be destabilized by dust-gas drift, but this instability requires nearly ideal conditions. We discuss the possible implications of these results for dust dynamics and planetesimal formation in protoplanetary disks. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Comments: Accepted by ApJ

Gravitational wave propagation beyond general relativity: waveform distortions and echoes

Authors: Jose Maria Ezquiaga, Wayne Hu, Macarena Lagos, Meng-Xiang Lin

Abstract: We study the cosmological propagation of gravitational waves (GWs) beyond general relativity (GR) across homogeneous and isotropic backgrounds. We consider scenarios in which GWs interact with an additional tensor field and use a parametrized phenomenological approach that generically describes their coupled equations of motion. We analyze four distinct classes of derivative and non-derivative int… ▽ More We study the cosmological propagation of gravitational waves (GWs) beyond general relativity (GR) across homogeneous and isotropic backgrounds. We consider scenarios in which GWs interact with an additional tensor field and use a parametrized phenomenological approach that generically describes their coupled equations of motion. We analyze four distinct classes of derivative and non-derivative interactions: mass, friction, velocity, and chiral. We apply the WKB formalism to account for the cosmological evolution and obtain analytical solutions to these equations. We corroborate these results by analyzing numerically the propagation of a toy GW signal. We then proceed to use the analytical results to study the modified propagation of realistic GWs from merging compact binaries, assuming that the GW signal emitted is the same as in GR. We generically find that tensor interactions lead to copies of the originally emitted GW signal, each one with its own possibly modified dispersion relation. These copies can travel coherently and interfere with each other leading to a scrambled GW signal, or propagate decoherently and lead to echoes arriving at different times at the observer that could be misidentified as independent GW events. Depending on the type of tensor interaction, the detected GW signal may exhibit amplitude and phase distortions with respect to a GW waveform in GR, as well as birefringence effects. We discuss observational probes of these tensor interactions with both individual GW events, as well as population studies for both ground- and space-based detectors. △ Less

Submitted 11 November, 2021; v1 submitted 24 August, 2021; originally announced August 2021.

Comments: Updated to journal version

DPNNet-2.0 Part I: Finding hidden planets from simulated images of protoplanetary disk gaps

Authors: Sayantan Auddy, Ramit Dey, Min-Kai Lin, Cassandra Hall

Abstract: The observed sub-structures, like annular gaps, in dust emissions from protoplanetary disk, are often interpreted as signatures of embedded planets. Fitting a model of planetary gaps to these observed features using customized simulations or empirical relations can reveal the characteristics of the hidden planets. However, customized fitting is often impractical owing to the increasing sample size… ▽ More The observed sub-structures, like annular gaps, in dust emissions from protoplanetary disk, are often interpreted as signatures of embedded planets. Fitting a model of planetary gaps to these observed features using customized simulations or empirical relations can reveal the characteristics of the hidden planets. However, customized fitting is often impractical owing to the increasing sample size and the complexity of disk-planet interaction. In this paper we introduce the architecture of DPNNet-2.0, second in the series after DPNNet \citep{aud20}, designed using a Convolutional Neural Network ( CNN, here specifically ResNet50) for predicting exoplanet masses directly from simulated images of protoplanetary disks hosting a single planet. DPNNet-2.0 additionally consists of a multi-input framework that uses both a CNN and multi-layer perceptron (a class of artificial neural network) for processing image and disk parameters simultaneously. This enables DPNNet-2.0 to be trained using images directly, with the added option of considering disk parameters (disk viscosities, disk temperatures, disk surface density profiles, dust abundances, and particle Stokes numbers) generated from disk-planet hydrodynamic simulations as inputs. This work provides the required framework and is the first step towards the use of computer vision (implementing CNN) to directly extract mass of an exoplanet from planetary gaps observed in dust-surface density maps by telescopes such as the Atacama Large (sub-)Millimeter Array. △ Less

Submitted 19 July, 2021; originally announced July 2021.

Comments: 15 pages, 10 figures, to appear in ApJ

On the Vertical Shear Instability in Magnetized Protoplanetary Disks

Authors: Can Cui, Min-Kai Lin

Abstract: The vertical shear instability (VSI) is a robust phenomenon in irradiated protoplanetary disks (PPDs). While there is extensive literature on the VSI in the hydrodynamic limit, PPDs are expected to be magnetized and their extremely low ionization fractions imply that non-ideal magneto-hydrodynamic (MHD) effects should be properly considered. To this end, we present linear analyses of the VSI in ma… ▽ More The vertical shear instability (VSI) is a robust phenomenon in irradiated protoplanetary disks (PPDs). While there is extensive literature on the VSI in the hydrodynamic limit, PPDs are expected to be magnetized and their extremely low ionization fractions imply that non-ideal magneto-hydrodynamic (MHD) effects should be properly considered. To this end, we present linear analyses of the VSI in magnetized disks with Ohmic resistivity. We primarily consider toroidal magnetic fields, which are likely to dominate the field geometry in PPDs. We perform vertically global and radially local analyses to capture characteristic VSI modes with extended vertical structures. To focus on the effect of magnetism, we use a locally isothermal equation of state. We find that magnetism provides a stabilizing effect to dampen the VSI, with surface modes, rather than body modes, being the first to vanish with increasing magnetization. Subdued VSI modes can be revived by Ohmic resistivity, where sufficient magnetic diffusion overcome magnetic stabilization, and hydrodynamic results are recovered. We also briefly consider poloidal fields to account for the magnetorotational instability (MRI), which may develop towards surface layers in the outer parts of PPDs. The MRI grows efficiently at small radial wavenumbers, in contrast to the VSI. When resistivity is considered, we find the VSI dominates over the MRI for Ohmic Elsässer numbers $\lesssim 0.09$ at plasma beta parameter $β_Z \sim 10^4$. △ Less

Submitted 24 May, 2021; originally announced May 2021.

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

LLAMA: Stellar populations in the nuclei of ultra hard X-ray selected AGN and matched inactive galaxies

Authors: L. Burtscher, R. I. Davies, T. T. Shimizu, R. Riffel, D. J. Rosario, E. K. S. Hicks, M. -Y. Lin, R. A. Riffel, M. Schartmann, A. Schnorr-Müller, T. Storchi-Bergmann, G. Orban de Xivry, S. Veilleux

Abstract: The relation between nuclear ($\lesssim$ 50 pc) star formation and nuclear galactic activity is still elusive: theoretical models predict a link between the two, but it is unclear whether active galactic nuclei (AGNs) should appear at the same time, before or after nuclear star formation activity is ongoing. We present a study of this relation in a complete, volume-limited sample of nine of the mo… ▽ More The relation between nuclear ($\lesssim$ 50 pc) star formation and nuclear galactic activity is still elusive: theoretical models predict a link between the two, but it is unclear whether active galactic nuclei (AGNs) should appear at the same time, before or after nuclear star formation activity is ongoing. We present a study of this relation in a complete, volume-limited sample of nine of the most luminous ($\log L_{\rm 14-195 keV} > 10^{42.5}$ erg/s) local AGNs (the LLAMA sample), including a sample of 18 inactive control galaxies (6 star-forming; 12 passive) that are matched by Hubble type, stellar mass (9.5 $\lesssim$ log M_star/M_sun $\lesssim$ 10.5), inclination and distance. This allows us to calibrate our methods on the control sample and perform a differential analysis between the AGN and control samples. We perform stellar population synthesis on VLT/X-SHOOTER spectra in an aperture corresponding to a physical radius of $\approx$ 150 pc. We find young ($\lesssim$ 30 Myr) stellar populations in seven out of nine AGNs and in four out of six star-forming control galaxies. In the non-star-forming control population, in contrast, only two out of twelve galaxies show such a population. We further show that these young populations are not indicative of ongoing star-formation, providing evidence for models that see AGN activity as a consequence of nuclear star formation. Based on the similar nuclear star-formation histories of AGNs and star-forming control galaxies, we speculate that the latter may turn into the former for some fraction of their time. Under this assumption, and making use of the volume-completeness of our sample, we infer that the AGN phase lasts for about 5 % of the nuclear starburst phase. △ Less

Submitted 11 May, 2021; originally announced May 2021.

Comments: 53 pages, 44 figures, accepted for publication by A&A

Journal ref: A&A 654, A132 (2021)

Which planets trigger longer-lived vortices: low-mass or high-mass?

Authors: Michael Hammer, Min-Kai Lin, Kaitlin M. Kratter, Paola Pinilla

Abstract: Recent ALMA observations have found many protoplanetary discs with rings that can be explained by gap-opening planets less massive than Jupiter. Meanwhile, recent studies have suggested that protoplanetary discs should have low levels of turbulence. Past computational work on low-viscosity discs has hinted that these two developments might not be self-consistent because even low-mass planets can b… ▽ More Recent ALMA observations have found many protoplanetary discs with rings that can be explained by gap-opening planets less massive than Jupiter. Meanwhile, recent studies have suggested that protoplanetary discs should have low levels of turbulence. Past computational work on low-viscosity discs has hinted that these two developments might not be self-consistent because even low-mass planets can be accompanied by vortices instead of conventional double rings. We investigate this potential discrepancy by conducting hydrodynamic simulations of growing planetary cores in discs with various aspect ratios ($H/r=0.04$, 0.06, 0.08) and viscosities ($1.5 \times 10^{-5} \lesssim α\lesssim 3 \times 10^{-4}$), having these cores accrete their gas mass directly from the disc. With $α< 10^{-4}$, we find that sub-Saturn-mass planets in discs with $H/r \le 0.06$ are more likely to be accompanied by dust asymmetries compared to Jupiter-mass planets because they can trigger several generations of vortices in succession. We also find that vortices with $H/r = 0.08$ survive $>6000$ planet orbits regardless of the planet mass or disc mass because they are less affected by the planet's spiral waves. We connect our results to observations and find that the outward migration of vortices with $H/r \ge 0.08$ may be able to explain the cavity in Oph IRS 48 or the two clumps in MWC 758. Lastly, we show that the lack of observed asymmetries in the disc population in Taurus is unexpected given the long asymmetry lifetimes in our low viscosity simulations ($α\sim 2 \times 10^{-5}$), a discrepancy we suggest is due to these discs having higher viscosities. △ Less

Submitted 9 May, 2021; v1 submitted 6 April, 2021; originally announced April 2021.

Comments: 24 pages, 20 figures, 2 tables; Accepted to MNRAS

Puffed up Edges of Planet-opened Gaps in Protoplanetary Disks. I. hydrodynamic simulations

Authors: Jiaqing Bi, Min-Kai Lin, Ruobing Dong

Abstract: Dust gaps and rings appear ubiquitous in bright protoplanetary disks. Disk-planet interaction with dust-trapping at the edges of planet-induced gaps is one plausible explanation. However, the sharpness of some observed dust rings indicate that sub-mm-sized dust grains have settled to a thin layer in some systems. We test whether or not such dust around gas gaps opened by planets can remain settled… ▽ More Dust gaps and rings appear ubiquitous in bright protoplanetary disks. Disk-planet interaction with dust-trapping at the edges of planet-induced gaps is one plausible explanation. However, the sharpness of some observed dust rings indicate that sub-mm-sized dust grains have settled to a thin layer in some systems. We test whether or not such dust around gas gaps opened by planets can remain settled by performing three-dimensional, dust-plus-gas simulations of protoplanetary disks with an embedded planet. We find planets massive enough to open gas gaps stir small, sub-mm-sized dust grains to high disk elevations at the gap edges, where the dust scale-height can reach ~70% of the gas scale-height. We attribute this dust 'puff-up' to the planet-induced meridional gas flows previously identified by Fung & Chiang and others. We thus emphasize the importance of explicit 3D simulations to obtain the vertical distribution of sub-mm-sized grains around gas gaps opened by massive planets. We caution that the gas-gap-opening planet interpretation of well-defined dust rings is only self-consistent with large grains exceeding mm in size. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: 15 pages, 7+2 figures, accepted in ApJ

Stratified and vertically-shearing streaming instabilities in protoplanetary disks

Authors: Min-Kai Lin

Abstract: Under the right conditions, the streaming instability between imperfectly coupled dust and gas is a powerful mechanism for planetesimal formation as it can concentrate dust grains to the point of gravitational collapse. In its simplest form, the streaming instability can be captured by analyzing the linear stability of unstratified disk models, which represent the midplane of protoplanetary disks.… ▽ More Under the right conditions, the streaming instability between imperfectly coupled dust and gas is a powerful mechanism for planetesimal formation as it can concentrate dust grains to the point of gravitational collapse. In its simplest form, the streaming instability can be captured by analyzing the linear stability of unstratified disk models, which represent the midplane of protoplanetary disks. We extend such studies by carrying out vertically-global linear stability analyses of dust layers in protoplanetary disks. We find the dominant form of instability in stratified dust layers is one driven by the vertical gradient in the rotation velocity of the dust-gas mixture, but also requires partial dust-gas coupling. These vertically-shearing streaming instabilities grow on orbital timescales and occur on radial length scales $\sim10^{-3}H_\mathrm{g}$, where $H_\mathrm{g}$ is the local pressure scale height. The classic streaming instability, associated with the relative radial drift between dust and gas, occur on radial length scales $\sim10^{-2}H_\mathrm{g}$, but have much smaller growth rates than vertically-shearing streaming instabilities. Including gas viscosity is strongly stabilizing and leads to vertically-elongated disturbances. We briefly discuss the potential effects of vertically-shearing streaming instabilities on planetesimal formation. △ Less

Submitted 24 November, 2020; originally announced November 2020.

Comments: Accepted by ApJ; project source codes are hosted at https://github.com/minkailin/stratsi; for a lay summary, see https://minkailin.wixsite.com/minkailin/post/streaming-instabilities-in-stratified-protoplanetary-disks

Testing H0 in Acoustic Dark Energy Models with Planck and ACT Polarization

Authors: Meng-Xiang Lin, Wayne Hu, Marco Raveri

Abstract: The canonical acoustic dark energy model (cADE), which is based on a scalar field with a canonical kinetic term that rapidly converts potential to kinetic energy around matter radiation equality, alleviates the Hubble tension found in $Λ$CDM. We show that it successfully passes new consistency tests in the CMB damping tail provided by the ACT data, while being increasingly constrained and distingu… ▽ More The canonical acoustic dark energy model (cADE), which is based on a scalar field with a canonical kinetic term that rapidly converts potential to kinetic energy around matter radiation equality, alleviates the Hubble tension found in $Λ$CDM. We show that it successfully passes new consistency tests in the CMB damping tail provided by the ACT data, while being increasingly constrained and distinguished from alternate mechanisms by the improved CMB acoustic polarization data from Planck. The best fit cADE model to a suite of cosmological observations, including the SH0ES $H_0$ measurement, has $H_0=70.25$ compared with $68.23$ (km s$^{-1}$ Mpc$^{-1}$) in $Λ$CDM and a finite cADE component is preferred at the $2.8σ$ level. The ability to raise $H_0$ is now mainly constrained by the improved Planck acoustic polarization data, which also plays a crucial role in distinguishing cADE from the wider class of early dark energy models. ACT and Planck TE polarization data are currently mildly discrepant in normalization and drive correspondingly different preferences in parameters. Improved constraints on intermediate scale polarization approaching the cosmic variance limit will be an incisive test of the acoustic dynamics of these models and their alternatives. △ Less

Submitted 14 December, 2020; v1 submitted 18 September, 2020; originally announced September 2020.

Comments: v2: Added more references; matched the published version. 10 pages, 8 figures

Journal ref: Phys. Rev. D 102, 123523 (2020)

A Machine Learning model to infer planet masses from gaps observed in protoplanetary disks

Authors: Sayantan Auddy, Min-Kai Lin

Abstract: Observations of bright protoplanetary disks often show annular gaps in their dust emission. One interpretation of these gaps is disk-planet interaction. If so, fitting models of planetary gaps to observed protoplanetary disk gaps can reveal the presence of hidden planets. However, future surveys are expected to produce an ever-increasing number of protoplanetary disks with gaps. In this case, perf… ▽ More Observations of bright protoplanetary disks often show annular gaps in their dust emission. One interpretation of these gaps is disk-planet interaction. If so, fitting models of planetary gaps to observed protoplanetary disk gaps can reveal the presence of hidden planets. However, future surveys are expected to produce an ever-increasing number of protoplanetary disks with gaps. In this case, performing a customized fitting for each target becomes impractical owing to the complexity of disk-planet interaction. To this end, we introduce DPNNet (Disk Planet Neural Network), an efficient model of planetary gaps by exploiting the power of machine learning. We train a deep neural network with a large number of dusty disk-planet hydrodynamic simulations across a range of planet masses, disk temperatures, disk viscosities, disk surface density profiles, particle Stokes numbers, and dust abundances. The network can then be deployed to extract the planet mass for a given gap morphology. In this work, first in a series, we focus on the basic concepts of our machine learning framework. We demonstrate its utility by applying it to the dust gaps observed in the protoplanetary disk around HL Tau at $10$ au, $30$ au, and $80$ au. Our network predict planet masses of $80 \, M_{\rm Earth}$, $63 \, M_{\rm Earth}$, and $70 \, M_{\rm Earth}$, respectively, which are comparable to other studies based on specialized simulations. We discuss the key advantages of our DPNNet in its flexibility to incorporate new physics, any number of parameters and predictions, and its potential to ultimately replace hydrodynamical simulations for disk observers and modelers. △ Less

Submitted 27 July, 2020; originally announced July 2020.

Comments: 16 Pages, 8 Figures, to appear in ApJ

Migrating Low-Mass Planets in Inviscid Dusty Protoplanetary Discs

Authors: He-Feng Hsieh, Min-Kai Lin

Abstract: Disc-driven planet migration is integral to the formation of planetary systems. In standard, gas-dominated protoplanetary discs, low-mass planets or planetary cores undergo rapid inwards migration and are lost to the central star. However, several recent studies indicate that the solid component in protoplanetary discs can have a significant dynamical effect on disc-planet interaction, especially… ▽ More Disc-driven planet migration is integral to the formation of planetary systems. In standard, gas-dominated protoplanetary discs, low-mass planets or planetary cores undergo rapid inwards migration and are lost to the central star. However, several recent studies indicate that the solid component in protoplanetary discs can have a significant dynamical effect on disc-planet interaction, especially when the solid-to-gas mass ratio approaches unity or larger and the dust-on-gas drag forces become significant. As there are several ways to raise the solid abundance in protoplanetary discs, for example through disc winds and dust-trapping in pressure bumps, it is important to understand how planets migrate through a dusty environment. To this end, we study planet migration in dust-rich discs via a systematic set of high-resolution, two-dimensional numerical simulations. We show that the inwards migration of low-mass planets can be slowed down by dusty dynamical corotation torques. We also identify a new regime of stochastic migration applicable to discs with dust-to-gas mass ratios $\gtrsim 0.3$ and particle Stokes numbers $\gtrsim 0.03$. In these cases, disc-planet interaction leads to the continuous development of small-scale, intense dust vortices that scatter the planet, which can potentially halt or even reverse the inwards planet migration. We briefly discuss the observational implications of our results and highlight directions for future work. △ Less

Submitted 16 July, 2020; originally announced July 2020.

Comments: Accepted by MNRAS

Ionized outflows in local luminous AGN: what are the real densities and outflow rates?

Authors: R. Davies, D. Baron, T. Shimizu, H. Netzer, L. Burtscher, P. T. de Zeeuw, R. Genzel, E. K. S. Hicks, M. Koss, M. -Y. Lin, D. Lutz, W. Maciejewski, F. Müller-Sánchez, G. Orban de Xivry, C. Ricci, R. Riffel, R. A. Riffel, D. Rosario, M. Schartmann, A. Schnorr-Müller, J. Shangguan, A. Sternberg, E. Sturm, T. Storchi-Bergmann, L. Tacconi , et al. (1 additional authors not shown)

Abstract: We report on the determination of electron densities, and their impact on the outflow masses and rates, measured in the central few hundred parsecs of 11 local luminous active galaxies. We show that the peak of the integrated line emission in the AGN is significantly offset from the systemic velocity as traced by the stellar absorption features, indicating that the profiles are dominated by outflo… ▽ More We report on the determination of electron densities, and their impact on the outflow masses and rates, measured in the central few hundred parsecs of 11 local luminous active galaxies. We show that the peak of the integrated line emission in the AGN is significantly offset from the systemic velocity as traced by the stellar absorption features, indicating that the profiles are dominated by outflow. In contrast, matched inactive galaxies are characterised by a systemic peak and weaker outflow wing. We present three independent estimates of the electron density in these AGN, discussing the merits of the different methods. The electron density derived from the [SII] doublet is significantly lower than than that found with a method developed in the last decade using auroral and transauroral lines, as well as a recently introduced method based on the ionization parameter. The reason is that, for gas photoionized by an AGN, much of the [SII] emission arises in an extended partially ionized zone where the implicit assumption that the electron density traces the hydrogen density is invalid. We propose ways to deal with this situation and we derive the associated outflow rates for ionized gas, which are in the range 0.001--0.5 M$_{\odot}$yr$^{-1}$ for our AGN sample. We compare these outflow rates to the relation between $\dot{M}_{out}$ and $L_{AGN}$ in the literature, and argue that it may need to be modified and rescaled towards lower mass outflow rates. △ Less

Submitted 11 August, 2020; v1 submitted 13 March, 2020; originally announced March 2020.

Comments: accepted by MNRAS; 30 pages (the last 8 of which are extra figures in an appendix). This is an update of the original version posted

How efficient is the streaming instability in viscous protoplanetary disks?

Authors: Kan Chen, Min-Kai Lin

Abstract: The streaming instability is a popular candidate for planetesimal formation by concentrating dust particles to trigger gravitational collapse. However, its robustness against physical conditions expected in protoplanetary disks is unclear. In particular, particle stirring by turbulence may impede the instability. To quantify this effect, we develop the linear theory of the streaming instability wi… ▽ More The streaming instability is a popular candidate for planetesimal formation by concentrating dust particles to trigger gravitational collapse. However, its robustness against physical conditions expected in protoplanetary disks is unclear. In particular, particle stirring by turbulence may impede the instability. To quantify this effect, we develop the linear theory of the streaming instability with external turbulence modelled by gas viscosity and particle diffusion. We find the streaming instability is sensitive to turbulence, with growth rates becoming negligible for alpha-viscosity parameters $α\gtrsim \mathrm{St} ^{1.5}$, where $\mathrm{St}$ is the particle Stokes number. We explore the effect of non-linear drag laws, which may be applicable to porous dust particles, and find growth rates are modestly reduced. We also find that gas compressibility increase growth rates by reducing the effect of diffusion. We then apply linear theory to global models of viscous protoplanetary disks. For minimum-mass Solar nebula disk models, we find the streaming instability only grows within disk lifetimes beyond $\sim 10$s of AU, even for cm-sized particles and weak turbulence ($α\sim 10^{-4}$). Our results suggest it is rather difficult to trigger the streaming instability in non-laminar protoplanetary disks, especially for small particles. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: Accepted by ApJ, see https://minkailin.wixsite.com/minkailin/post/making-planets-in-turbulent-disks-is-not-so-easy for a lay summary

LLAMA: The $M_{BH}$ - $σ_{\star}$ Relation of the most luminous local AGNs

Authors: Turgay Caglar, Leonard Burtscher, Bernhard Brand, Jarle Brinchmann, Richard I. Davies, Erin K. S. Hicks, Michael Koss, Ming-Yi Lin, Witold Maciejewski, Francisco Müller-Sánchez, Rogemar A. Riffel, Rogério Riffel, David J. Rosario, Marc Schartmann, Allan Schnorr-Müller, T. Taro Shimizu, Thaisa Storchi-Bergmann, Sylvain Veilleux, Gilles O. de Xivry, Vardha N. Bennert

Abstract: The $M_{BH}$ - $σ_{\star}$ relation is considered a result of co-evolution between the host galaxies and their super-massive black holes. For elliptical-bulge hosting inactive galaxies, this relation is well established, but there is still a debate whether active galaxies follow the same relation. In this paper, we estimate black hole masses for a sample of 19 local luminous AGNs (LLAMA) in order… ▽ More The $M_{BH}$ - $σ_{\star}$ relation is considered a result of co-evolution between the host galaxies and their super-massive black holes. For elliptical-bulge hosting inactive galaxies, this relation is well established, but there is still a debate whether active galaxies follow the same relation. In this paper, we estimate black hole masses for a sample of 19 local luminous AGNs (LLAMA) in order to test their location on the $M_{BH}$ - $σ_{\star}$ relation. Super-massive black hole masses ($M_{BH}$) were derived from the broad-line based relations for H$α$, H$β$ and Pa$β$ emission line profiles for the Type 1 AGNs. We compare the bulge stellar velocity dispersion ($σ_{\star}$) as determined from the Ca II triplet (CaT) with the dispersion measured from the near-infrared CO (2-0) absorption features for each AGN and find them to be consistent with each other. We apply an extinction correction to the observed broad line fluxes and we correct the stellar velocity dispersion by an average rotation contribution as determined from spatially resolved stellar kinematic maps. The H$α$-based black hole masses of our sample of AGNs were estimated in the range 6.34 $\leq$ $\log{M_{BH}}$ $\leq$ 7.75 M$_\odot$ and the $σ_{\star CaT}$ estimates range between 73 $\leq$ $σ_{\star CaT}$ $\leq$ 227 km s$^{-1}$. From the so-constructed $M_{BH}$ - $σ_{\star}$ relation for our Type 1 AGNs, we estimate the black hole masses for the Type 2 AGNs and the inactive galaxies in our sample. In conclusion, we find that our sample of local luminous AGNs is consistent with the $M_{BH}$ - $σ_{\star}$ relation of lower luminosity AGNs and inactive galaxies, after correcting for dust extinction and the rotational contribution to the stellar velocity dispersion. △ Less

Submitted 16 December, 2019; originally announced December 2019.

Comments: Accepted by A&A, 30 pages, 19 figures

Journal ref: A&A 634, A114 (2020)

SCUBA-2 Ultra Deep Imaging Eao Survey (Studies) III: Multi-wavelength properties, luminosity functions and preliminary source catalog of 450-$μ$m-selected galaxies

Authors: Chen-Fatt Lim, Wei-Hao Wang, Ian Smail, Douglas Scott, Chian-Chou Chen, Yu-Yen Chang, James M. Simpson, Yoshiki Toba, Xinwen Shu, Dave Clements, Josh Greenslade, YiPing Ao, Arif Babul, Jack Birkin, Scott C. Chapman, Tai-An Cheng, Brian S. Cho, Helmut Dannerbauer, Ugnė Dudzevičiūtė, James Dunlop, Yu Gao, Tomotsugu Goto, Luis C. Ho, Li-Ting Hsu, Ho Seong Hwang , et al. (13 additional authors not shown)

Abstract: We construct a SCUBA-2 450-$μ$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$σ$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $μ$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of… ▽ More We construct a SCUBA-2 450-$μ$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$σ$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $μ$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of $z = 1.79^{+0.03}_{-0.15}$. About $35^{+32}_{-25}$% of our sources are classified as starburst galaxies based on their total star-formation rates (SFRs) and stellar masses ($M_{\ast}$). By fitting the far-infrared spectral energy distributions, we find that our 450-$μ$m-selected sample has a wide range of dust temperatures (20 K $ \lesssim T_{\rm d} \lesssim$ 60 K), with a median of ${T}_{\rm d} = 38.3^{+0.4}_{-0.9}$ K. We do not find a redshift evolution in dust temperature for sources with $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ at $z<3$. However, we find a moderate correlation where dust temperature increases with the deviation from the SFR-$M_{\ast}$ relation. The increase in dust temperature also correlates with optical morphology, which is consistent with merger-triggered starbursts in sub-millimeter galaxies. Our galaxies do not show the tight IRX-$β_{\rm UV}$ correlation that has been observed in the local Universe. We construct the infrared luminosity functions of our 450-$μ$m sources and measure their comoving SFR densities. The contribution of the $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ population to the SFR density rises dramatically from $z$ = 0 to 2 ($\propto$ ($1+z$)$^{3.9\pm1.1}$) and dominates the total SFR density at $z \gtrsim 2$. △ Less

Submitted 8 March, 2020; v1 submitted 8 December, 2019; originally announced December 2019.

Comments: Published in the ApJ

Journal ref: ApJ 889 (2020) 80

Curvature Perturbations in the Effective Field Theory of Inflation

Authors: Macarena Lagos, Meng-Xiang Lin, Wayne Hu

Abstract: We discuss the difference between various gauge-invariant quantities typically used in single-field inflation, namely synchronous $ζ_s$, comoving $ζ_c$, and unitary $ζ_u$ curvatures. We show that conservation of $ζ_c$ outside the horizon is quite restrictive on models as it leads to conservation of $ζ_s$ and $ζ_u$, whereas the reverse does not hold. We illustrate the consequence of these differenc… ▽ More We discuss the difference between various gauge-invariant quantities typically used in single-field inflation, namely synchronous $ζ_s$, comoving $ζ_c$, and unitary $ζ_u$ curvatures. We show that conservation of $ζ_c$ outside the horizon is quite restrictive on models as it leads to conservation of $ζ_s$ and $ζ_u$, whereas the reverse does not hold. We illustrate the consequence of these differences with two inflationary models: ultra-slow-roll (USR) and braiding-ultra-slow-roll (BUSR). In USR, we show that out of the three curvatures, only $ζ_s$ is conserved outside the horizon, and we connect this result to the concepts of separate universe and the usage of the $δN$ formalism. We find that even though $ζ_s$ is conserved, there is still a mild violation of the separate universe approximation in the continuity equation. Nevertheless, the $δN$ formalism can still be applied to calculate the primordial power spectrum of some gauge-invariant quantities such as $ζ_u$, although it breaks down for others such as the uniform-density curvature. In BUSR, we show that both $ζ_u$ and $ζ_s$ are conserved outside the horizon, but take different values. Additionally, since $ζ_u\not=ζ_c$ we find that the prediction for observable curvature fluctuations after inflation does not reflect $ζ_c$ at horizon crossing during inflation and moreover involves not just $ζ_u$ at that epoch but also the manner in which the braiding phase ends. △ Less

Submitted 22 August, 2019; originally announced August 2019.

Journal ref: Phys. Rev. D 100, 123507 (2019)

The multiphase gas structure and kinematics in the circumnuclear region of NGC 5728

Authors: T. Taro Shimizu, R. I. Davies, D. Lutz, L. Burtscher, M. Lin, D. Baron, R. L. Davies, R. Genzel, E. K. S. Hicks, M. Koss, W. Maciejewski, F. Mu ller-Sanchez, G. Orban de Xivry, S. H. Price, C. Ricci, R. Riffel, R. A. Riffel, D. Rosario, M. Schartmann, A. Schnorr-Muller, A. Sternberg, E. Sturm, T. Storchi-Bergmann, L. Tacconi, S. Veilleux

Abstract: We report on our combined analysis of HST, VLT/MUSE, VLT/SINFONI, and ALMA observations of the local Seyfert 2 galaxy, NGC 5728 to investigate in detail the feeding and feedback of the AGN. The datasets simultaneously probe the morphology, excitation, and kinematics of the stars, ionized gas, and molecular gas over a large range of spatial scales (10 pc--10 kpc). NGC 5728 contains a large stellar… ▽ More We report on our combined analysis of HST, VLT/MUSE, VLT/SINFONI, and ALMA observations of the local Seyfert 2 galaxy, NGC 5728 to investigate in detail the feeding and feedback of the AGN. The datasets simultaneously probe the morphology, excitation, and kinematics of the stars, ionized gas, and molecular gas over a large range of spatial scales (10 pc--10 kpc). NGC 5728 contains a large stellar bar which is driving gas along prominent dust lanes to the inner 1 kpc where the gas settles into a circumnuclear ring. The ring is strongly star forming and contains a substantial population of young stars as indicated by the lowered stellar velocity dispersion and gas excitation consistent with HII regions. We model the kinematics of the ring using the velocity field of the CO (2--1) emission and stars and find it is consistent with a rotating disk. The outer regions of the disk, where the dust lanes meet the ring, show signatures of inflow at a rate of 1 M$_{\sun}$ yr$^{-1}$. Inside the ring, we observe three molecular gas components corresponding to the circular rotation of the outer ring, a warped disk, and the nuclear stellar bar. The AGN is driving an ionized gas outflow that reaches a radius of 250 pc with a mass outflow rate of 0.08 M$_{\sun}$ yr$^{-1}$ consistent with its luminosity and scaling relations from previous studies. While we observe distinct holes in CO emission which could be signs of molecular gas removal, we find that largely the AGN is not disrupting the structure of the circumnuclear region. △ Less

Submitted 2 October, 2019; v1 submitted 8 July, 2019; originally announced July 2019.

Comments: accepted for publication in MNRAS; 22 pages + Appendix, 16 figures

Vortex instabilities triggered by low-mass planets in pebble-rich, inviscid protoplanetary discs

Authors: Arnaud Pierens, Min-Kai Lin, Sean Raymond

Abstract: In the innermost regions of protoplanerary discs, the solid-to-gas ratio can be increased considerably by a number of processes, including photoevaporative and particle drift. MHD disc models also suggest the existence of a dead-zone at $R\lesssim 10$ AU, where the regions close to the midplane remain laminar. In this context, we use two-fluid hydrodynamical simulations to study the interaction be… ▽ More In the innermost regions of protoplanerary discs, the solid-to-gas ratio can be increased considerably by a number of processes, including photoevaporative and particle drift. MHD disc models also suggest the existence of a dead-zone at $R\lesssim 10$ AU, where the regions close to the midplane remain laminar. In this context, we use two-fluid hydrodynamical simulations to study the interaction between a low-mass planet ($\sim 1.7 \;{\rm M_\oplus}$) on a fixed orbit and an inviscid pebble-rich disc with solid-to-gas ratio $ε\ge 0.5$. For pebbles with Stokes numbers St=0.1, 0.5, multiple dusty vortices are formed through the Rossby Wave Instability at the planet separatrix. Effects due to gas drag then lead to a strong enhancement in the solid-to-gas ratio, which can increase by a factor of $\sim 10^3$ for marginally coupled particles with St=0.5. As in streaming instabilities, pebble clumps reorganize into filaments that may plausibly collapse to form planetesimals. When the planet is allowed to migrate in a MMSN disc, the vortex instability is delayed due to migration but sets in once inward migration stops due a strong positive pebble torque. Again, particle filaments evolving in a gap are formed in the disc while the planet undergoes an episode of outward migration. Our results suggest that vortex instabilities triggered by low-mass planets could play an important role in forming planetesimals in pebble-rich, inviscid discs, and may significantly modify the migration of low-mass planets. They also imply that planetary dust gaps may not necessarily contain planets if these migrated away. △ Less

Submitted 19 June, 2019; originally announced June 2019.

Comments: Accepted in MNRAS

Acoustic Dark Energy: Potential Conversion of the Hubble Tension

Authors: Meng-Xiang Lin, Giampaolo Benevento, Wayne Hu, Marco Raveri

Abstract: We discuss the ability of a dark fluid becoming relevant around the time of matter radiation equality to significantly relieve the tension between local measurements of the Hubble constant and CMB inference, within the $Λ$CDM model. We show that the gravitational impact of acoustic oscillations in the dark fluid balance the effects on the CMB and result in an improved fit to CMB measurements thems… ▽ More We discuss the ability of a dark fluid becoming relevant around the time of matter radiation equality to significantly relieve the tension between local measurements of the Hubble constant and CMB inference, within the $Λ$CDM model. We show that the gravitational impact of acoustic oscillations in the dark fluid balance the effects on the CMB and result in an improved fit to CMB measurements themselves while simultaneously raising the Hubble constant. The required balance favors a model where the fluid is a scalar field that converts its potential to kinetic energy around matter radiation equality which then quickly redshifts away. We derive the requirements on the potential for this conversion mechanism and find that a simple canonical scalar with two free parameters for its local slope and amplitude robustly improves the fit to the combined data by $Δχ^2 \approx 12.7$ over $Λ$CDM. We uncover the CMB polarization signatures that can definitively test this scenario with future data. △ Less

Submitted 29 October, 2019; v1 submitted 29 May, 2019; originally announced May 2019.

Comments: v3: matched the published version

Journal ref: Phys. Rev. D 100, 063542 (2019)

Dust settling against hydrodynamic turbulence in protoplanetary discs

Authors: Min-Kai Lin

Abstract: Enhancing the local dust-to-gas ratio in protoplanetary discs is a necessary first step to planetesimal formation. In laminar discs, dust settling is an efficient mechanism to raise the dust-to-gas ratio at the disc midplane. However, turbulence, if present, can stir and lift dust particles, which ultimately hinders planetesimal formation. In this work, we study dust settling in protoplanetary dis… ▽ More Enhancing the local dust-to-gas ratio in protoplanetary discs is a necessary first step to planetesimal formation. In laminar discs, dust settling is an efficient mechanism to raise the dust-to-gas ratio at the disc midplane. However, turbulence, if present, can stir and lift dust particles, which ultimately hinders planetesimal formation. In this work, we study dust settling in protoplanetary discs with hydrodynamic turbulence sustained by the vertical shear instability. We perform axisymmetric numerical simulations to investigate the effect of turbulence, particle size, and solid abundance or metallicity on dust settling. We highlight the positive role of drag forces exerted onto the gas by the dust for settling to overcome the vertical shear instability. In typical disc models we find particles with a Stokes number $\sim 10^{-3}$ can sediment to $\lesssim 10\%$ of the gas scale-height, provided that $Σ_\mathrm{d}/Σ_\mathrm{g}\gtrsim 0.02$-$0.05$, where $Σ_\mathrm{d,g}$ are the surface densities in dust and gas, respectively. This coincides with the metallicity condition for small particles to undergo clumping via the streaming instability. Super-solar metallicities, at least locally, are thus required for a self-consistent picture of planetesimal formation. Our results also imply that dust rings observed in protoplanetary discs should have smaller scale-heights than dust gaps, provided that the metallicity contrast between rings and gaps exceed the corresponding contrast in gas density. △ Less

Submitted 8 March, 2019; originally announced March 2019.

Comments: Accepted by MNRAS; see https://minkailin.wixsite.com/minkailin/blog/dust-settling-in-turbulent-discs for a summary

Phenomenology of Modified Gravity at Recombination

Authors: Meng-Xiang Lin, Marco Raveri, Wayne Hu

Abstract: We discuss the phenomenological imprints of modifications to gravity in the early universe with a specific focus on the time of recombination. We derive several interesting results regarding the effect that such modifications have on cosmological observables, especially on the driving and phasing of acoustic oscillations, observed in the CMB and BAO, as well as the weak gravitational lensing of th… ▽ More We discuss the phenomenological imprints of modifications to gravity in the early universe with a specific focus on the time of recombination. We derive several interesting results regarding the effect that such modifications have on cosmological observables, especially on the driving and phasing of acoustic oscillations, observed in the CMB and BAO, as well as the weak gravitational lensing of the CMB and of galaxy shapes. This widens the pool of measurements that can be used to test gravity with present and future surveys, in particular realizing the full constraining power of the structure of the primary peaks of the CMB spectrum. We investigate whether such a phenomenology can relax tensions between cosmological measurements and find that a modification of the gravitational constant at recombination would help in reconciling measurements of the CMB with local measurements of the Hubble constant. △ Less

Submitted 4 October, 2018; originally announced October 2018.

Comments: 22 pages, 21 figures

Journal ref: Phys. Rev. D 99, 043514 (2019)

Observational diagnostics of elongated planet-induced vortices with realistic planet formation timescales

Authors: Michael Hammer, Paola Pinilla, Kaitlin M. Kratter, Min-Kai Lin

Abstract: Gap-opening planets can generate dust-trapping vortices that may explain some of the latest discoveries of high-contrast crescent-shaped dust asymmetries in transition discs. While planet-induced vortices were previously thought to have concentrated shapes, recent computational work has shown that these features naturally become much more elongated in the gas when simulations account for the relat… ▽ More Gap-opening planets can generate dust-trapping vortices that may explain some of the latest discoveries of high-contrast crescent-shaped dust asymmetries in transition discs. While planet-induced vortices were previously thought to have concentrated shapes, recent computational work has shown that these features naturally become much more elongated in the gas when simulations account for the relatively long timescale over which planets accrete their mass. In this work, we conduct two-fluid hydrodynamical simulations of vortices induced by slowly-growing Jupiter-mass planets in discs with very low viscosity ($α= 3 \times 10^{-5}$). We simulate the dust dynamics for four particle sizes spanning 0.3 mm to 1 cm in order to produce synthetic ALMA images. In our simulations, we find that an elongated vortex still traps dust, but not directly at its center. With a flatter pressure bump and disruptions from the planet's overlapping spiral density waves, the dust instead circulates around the vortex. This motion (1) typically carries the peak off-center, (2) spreads the dust out over a wider azimuthal extent $\geq 180^{\circ}$, (3) skews the azimuthal profile towards the front of the vortex, and (4) can also create double peaks in newly-formed vortices. In particular, we expect that the most defining observational signature, a peak offset of more than $30^{\circ}$, should be detectable $>30\%$ of the time in observations with a beam diameter of at most the planet's separation from its star. △ Less

Submitted 29 October, 2018; v1 submitted 19 September, 2018; originally announced September 2018.

Comments: Accepted to MNRAS. 13 pages, 8 figures. Movies available at: https://lavinia.as.arizona.edu/~mhammer/vortex_signatures.html

Vortex survival in 3D self-gravitating accretion discs

Authors: Min-Kai Lin, Arnaud Pierens

Abstract: Large-scale, dust-trapping vortices may account for observations of asymmetric protoplanetary discs. Disc vortices are also potential sites for accelerated planetesimal formation by concentrating dust grains. However, in 3D discs vortices are subject to destructive `elliptic instabilities', which reduces their viability as dust traps. The survival of vortices in 3D accretion discs is thus an impor… ▽ More Large-scale, dust-trapping vortices may account for observations of asymmetric protoplanetary discs. Disc vortices are also potential sites for accelerated planetesimal formation by concentrating dust grains. However, in 3D discs vortices are subject to destructive `elliptic instabilities', which reduces their viability as dust traps. The survival of vortices in 3D accretion discs is thus an important issue to address. In this work, we perform shearing box simulations to show that disc self-gravity enhances the survival of 3D vortices, even when self-gravity is weak in the classic sense (e.g. with a Toomre $Q\simeq5$). We find a 3D, self-gravitating vortex can grow on secular timescales in spite of the elliptic instability. The vortex aspect-ratio decreases as it strengthens, which feeds the elliptic instability. The result is a 3D vortex with a turbulent core that persists for $\sim 10^{3}$ orbits. We find when gravitational and hydrodynamic stresses become comparable, the vortex may undergo episodic bursts, which we interpret as interaction between elliptic and gravitational instabilities. We estimate the distribution of dust particles in self-gravitating, turbulent vortices. Our results suggest large-scale vortices in protoplanetary discs are more easily observed at large radii. △ Less

Submitted 15 May, 2018; originally announced May 2018.

Comments: Accepted by MNRAS; full resolution article at https://docs.wixstatic.com/ugd/f3eb2c_500e5bc1e9354ba295e7b954b8f5247e.pdf

On the evolution of vortices in massive protoplanetary discs

Authors: Arnaud Pierens, Min-Kai Lin

Abstract: It is expected that a pressure bump can be formed at the inner edge of a dead-zone, and where vortices can develop through the Rossby Wave Instability (RWI). It has been suggested that self-gravity can significantly affect the evolution of such vortices. We present the results of 2D hydrodynamical simulations of the evolution of vortices forming at a pressure bump in self-gravitating discs with To… ▽ More It is expected that a pressure bump can be formed at the inner edge of a dead-zone, and where vortices can develop through the Rossby Wave Instability (RWI). It has been suggested that self-gravity can significantly affect the evolution of such vortices. We present the results of 2D hydrodynamical simulations of the evolution of vortices forming at a pressure bump in self-gravitating discs with Toomre parameter in the range $4-30$. We consider isothermal plus non-isothermal disc models that employ either the classical $β$ prescription or a more realistic treatment for cooling. The main aim is to investigate whether the condensating effect of self-gravity can stabilize vortices in sufficiently massive discs. We confirm that in isothermal disc models with ${\cal Q} \gtrsim 15$, vortex decay occurs due to the vortex self-gravitational torque. For discs with $3\lesssim {\cal Q} \lesssim 7$, the vortex develops gravitational instabilities within its core and undergoes gravitational collapse, whereas more massive discs give rise to the formation of global eccentric modes. In non-isothermal discs with $β$ cooling, the vortex maintains a turbulent core prior to undergoing gravitational collapse for $β\lesssim 0.1$, whereas it decays if $β\ge 1$. In models that incorpore both self-gravity and a better treatment for cooling, however, a stable vortex is formed with aspect ratio $χ\sim 3-4$. Our results indicate that self-gravity significantly impacts the evolution of vortices forming in protoplanetary discs, although the thermodynamical structure of the vortex is equally important for determining its long-term dynamics. △ Less

Submitted 15 May, 2018; originally announced May 2018.

Comments: Accepted in MNRAS

Dusty disc-planet interaction with dust-free simulations

Authors: Jhih-Wei Chen, Min-Kai Lin

Abstract: Protoplanets may be born into dust-rich environments if planetesimals formed through streaming or gravitational instabilities, or if the protoplanetary disc is undergoing mass loss due to disc winds or photoevaporation. Motivated by this possibility, we explore the interaction between low mass planets and dusty protoplanetary discs with focus on disc-planet torques. We implement Lin & Youdin's new… ▽ More Protoplanets may be born into dust-rich environments if planetesimals formed through streaming or gravitational instabilities, or if the protoplanetary disc is undergoing mass loss due to disc winds or photoevaporation. Motivated by this possibility, we explore the interaction between low mass planets and dusty protoplanetary discs with focus on disc-planet torques. We implement Lin & Youdin's newly developed, purely hydrodynamic model of dusty gas into the PLUTO code to simulate dusty protoplanetary discs with an embedded planet. We find that for imperfectly coupled dust and high metallicity, e.g. Stokes number $10^{-3}$ and dust-to-gas ratio $Σ_\mathrm{d}/Σ_\mathrm{g}=0.5$ , a `bubble' develops inside the planet's co-orbital region, which introduces unsteadiness in the flow. The resulting disc-planet torques sustain large amplitude oscillations that persists well beyond that in simulations with perfectly coupled dust or low dust-loading, where co-rotation torques are always damped. We show that the desaturation of the co-rotation torques by finite-sized particles is related to potential vorticity generation from the misalignment of dust and gas densities. We briefly discuss possible implications for the orbital evolution of protoplanets in dust-rich discs. We also demonstrate Lin & Youdin's dust-free framework reproduces previous results pertaining to dusty protoplanetary discs, including dust-trapping by pressure bumps, dust settling, and the streaming instability. △ Less

Submitted 2 May, 2018; originally announced May 2018.

Comments: Accepted by MNRAS