More Like this

1. Length and Velocity Scales in Protoplanetary Disk Turbulence

   https://doi.org/10.3847/1538-4357/ad2c89  

   Sengupta, Debanjan; Cuzzi, Jeffrey N; Umurhan, Orkan M; Lyra, Wladimir (April 2024, The Astrophysical Journal)

   Abstract In the theory of protoplanetary disk turbulence, a widely adopted ansatz, or assumption, is that the turnover frequency of the largest turbulent eddy, Ω_L, is the local Keplerian frequency Ω_K. In terms of the standard dimensionless Shakura–Sunyaevαparameter that quantifies turbulent viscosity or diffusivity, this assumption leads to characteristic length and velocity scales given respectively by $\sqrt{\alpha }H$and $\sqrt{\alpha }c$, in whichHandcare the local gas scale height and sound speed. However, this assumption is not applicable in cases when turbulence is forced numerically or driven by some natural processes such as vertical shear instability. Here, we explore the more general case where Ω_L≥ Ω_Kand show that, under these conditions, the characteristic length and velocity scales are respectively $\sqrt{\alpha /R\prime }H$and $\sqrt{\alpha R\prime }c$, where $R\prime \equiv {\mathrm{\Omega }}_L/{\mathrm{\Omega }}_K$is twice the Rossby number. It follows that $\alpha =\tilde{\alpha }/R\prime$, where $\sqrt{\tilde{\alpha }}c$is the root-mean-square average of the turbulent velocities. Properly allowing for this effect naturally explains the reduced particle scale heights produced in shearing box simulations of particles in forced turbulence, and it may help with interpreting recent edge-on disk observations; more general implications for observations are also presented. For $R\prime >1$, the effective particle Stokes numbers are increased, which has implications for particle collision dynamics and growth, as well as for planetesimal formation. 

   more » « less
2. Divergence beneath the Brillouin sphere and the phenomenology of prediction error in spherical harmonic series approximations of the gravitational field

   https://doi.org/10.1088/1361-6633/ad44d5  

   Bevis, M; Ogle, C; Costin, O; Jekeli, C; Costin, R D; Guo, J; Fowler, J; Dunne, G V; Shum, C K; Snow, K (June 2024, Reports on Progress in Physics)

   The Brillouin sphere is defined as the smallest sphere, centered at the origin of the geocentric coordinate system, that incorporates all the condensed matter composing the planet. The Brillouin sphere touches the Earth at a single point, and the radial line that begins at the origin and passes through that point is called the singular radial line. For about 60 years there has been a persistent anxiety about whether or not a spherical harmonic (SH) expansion of the external gravitational potential,V, will converge beneath the Brillouin sphere. Recently, it was proven that the probability of such convergence is zero. One of these proofs provided an asymptotic relation, called Costin’s formula, for the upper bound,E_N, on the absolute value of the prediction error,e_N, of a SH series model, $V_N(\theta ,\lambda ,r)$, truncated at some maximum degree, $N=n_{max}$. When the SH series is restricted to (or projected onto) a particular radial line, it reduces to a Taylor series (TS) in $1/r$. Costin’s formula is $E_N\simeq B{N}^{-b}(R/r{)}^N$, whereRis the radius of the Brillouin sphere. This formula depends on two positive parameters:b, which controls the decay of error amplitude as a function ofNwhenris fixed, and a scale factorB. We show here that Costin’s formula derives from a similar asymptotic relation for the upper bound,A_non the absolute value of the TS coefficients,a_n, for the same radial line. This formula, $A_n\simeq K{n}^{-k}$, depends on degree,n, and two positive parameters,kandK, that are analogous tobandB. We use synthetic planets, for which we can compute the potential,V, and also the radial component of gravitational acceleration, $g_r=\partial V/\partial r$, to hundreds of significant digits, to validate both of these asymptotic formulas. Let superscriptVrefer to asymptotic parameters associated with the coefficients and prediction errors for gravitational potential, and superscriptgto the coefficients and predictions errors associated withg_r. For polyhedral planets of uniform density we show that $b^V=k^V=7/2$and $b^g=k^g=5/2$almost everywhere. We show that the frequency of oscillation (around zero) of the TS coefficients and the series prediction errors, for a given radial line, is controlled by the geocentric angle,α, between that radial line and the singular radial line. We also derive useful identities connecting $K^V,B^V,K^g$, andB^g. These identities are expressed in terms of quotients of the various scale factors. The only other quantities involved in these identities areαandR. The phenomenology of ‘series divergence’ and prediction error (whenr < R) can be described as a function of the truncation degree,N, or the depth,d, beneath the Brillouin sphere. For a fixed $r\text{⩽}R$, asNincreases from very low values, the upper error boundE_Nshrinks until it reaches its minimum (best) value whenNreaches some particular or optimum value, $N_{\mathrm{opt}}$. When $N>N_{\mathrm{opt}}$, prediction error grows asNcontinues to increase. Eventually, when $N\gg N_{\mathrm{opt}}$, prediction errors increase exponentially with risingN. If we fix the value ofNand allow $R/r$to vary, then we find that prediction error in free space beneath the Brillouin sphere increases exponentially with depth,d, beneath the Brillouin sphere. Because $b^g=b^V-1$everywhere, divergence driven prediction error intensifies more rapidly forg_rthan forV, both in terms of its dependence onNandd. If we fix bothNandd, and focus on the ‘lateral’ variations in prediction error, we observe that divergence and prediction error tend to increase (as doesB) as we approach high-amplitude topography. 

   more » « less
3. Acceleration and Spectral Redistribution of Cosmic Rays in Radio-jet Shear Flows

   https://doi.org/10.3847/1538-4357/acfda9  

   Webb, G. M.; Xu, Y.; Biermann, P. L.; Al-Nussirat, S.; Mostafavi, P.; Li, G.; Barghouty, A. F.; Zank, G. P. (November 2023, The Astrophysical Journal)

   Abstract A steady-state, semi-analytical model of energetic particle acceleration in radio-jet shear flows due to cosmic-ray viscosity obtained by Webb et al. is generalized to take into account more general cosmic-ray boundary spectra. This involves solving a mixed Dirichlet–Von Neumann boundary value problem at the edge of the jet. The energetic particle distribution functionf₀(r,p) at cylindrical radiusrfrom the jet axis (assumed to lie along thez-axis) is given by convolving the particle momentum spectrum $f_0(\infty ,p\prime )$with the Green’s function $G(r,p;p\prime )$, which describes the monoenergetic spectrum solution in which $f_0\to \delta (p-p\prime )$asr→ ∞ . Previous work by Webb et al. studied only the Green’s function solution for $G(r,p;p\prime )$. In this paper, we explore for the first time, solutions for more general and realistic forms for $f_0(\infty ,p\prime )$. The flow velocityu=u(r)e_zis along the axis of the jet (thez-axis).uis independent ofz, andu(r) is a monotonic decreasing function ofr. The scattering time $\tau {(r,p)={\tau }_0\left(p/p_0\right)}^{\alpha }$in the shear flow region 0 <r<r₂, and $\tau {(r,p)={\tau }_0\left(p/p_0\right)}^{\alpha }{\left(r/r_2\right)}^s$, wheres> 0 in the regionr>r₂is outside the jet. Other original aspects of the analysis are (i) the use of cosmic ray flow lines in (r,p) space to clarify the particle spatial transport and momentum changes and (ii) the determination of the probability distribution ${\psi }_p(r,p;p\prime )$that particles observed at (r,p) originated fromr→ ∞ with momentum $p\prime$. The acceleration of ultrahigh-energy cosmic rays in active galactic nuclei jet sources is discussed. Leaky box models for electron acceleration are described. 

   more » « less
4. The Bimodal Absorption System Imaging Campaign (BASIC). I. A Dual Population of Low-metallicity Absorbers at z < 1

   https://doi.org/10.3847/1538-4357/acb047  

   Berg, Michelle_A; Lehner, Nicolas; Howk, J_Christopher; O’Meara, John_M; Schaye, Joop; Straka, Lorrie_A; Cooksey, Kathy_L; Tripp, Todd_M; Prochaska, J_Xavier; Oppenheimer, Benjamin_D; et al (February 2023, The Astrophysical Journal)

   Abstract The bimodal absorption system imaging campaign (BASIC) aims to characterize the galaxy environments of a sample of 36 Hi-selected partial Lyman limit systems (pLLSs) and Lyman limit systems (LLSs) in 23 QSO fields atz≲ 1. These pLLSs/LLSs provide a unique sample of absorbers with unbiased and well-constrained metallicities, allowing us to explore the origins of metal-rich and low-metallicity circumgalactic medium (CGM) atz< 1. Here we present Keck/KCWI and Very Large Telescope/MUSE observations of 11 of these QSO fields (19 pLLSs) that we combine with Hubble Space Telescope/Advanced Camera for Surveys imaging to identify and characterize the absorber-associated galaxies at 0.16 ≲z≲ 0.84. We find 23 unique absorber-associated galaxies, with an average of one associated galaxy per absorber. For seven absorbers, all with <10% solar metallicities, we find no associated galaxies with $\log M_{\star }\gtrsim 9.0$withinρ/R_virand ∣Δv∣/v_esc≤ 1.5 with respect to the absorber. We do not find any strong correlations between the metallicities or Hicolumn densities of the gas and most of the galaxy properties, except for the stellar mass of the galaxies: the low-metallicity ([X/H] ≤ −1.4) systems have a probability of ${0.39}_{-0.15}^{+0.16}$for having a host galaxy with $\log M_{\star }\ge 9.0$withinρ/R_vir≤ 1.5, while the higher metallicity absorbers have a probability of ${0.78}_{-0.13}^{+0.10}$. This implies metal-enriched pLLSs/LLSs atz< 1 are typically associated with the CGM of galaxies with $\log M_{\star }>9.0$, whereas low-metallicity pLLSs/LLSs are found in more diverse locations, with one population arising in the CGM of galaxies and another more broadly distributed in overdense regions of the universe. Using absorbers not associated with galaxies, we estimate the unweighted geometric mean metallicity of the intergalactic medium to be [X/H] ≲ −2.1 atz< 1, which is lower than previously estimated. 

   more » « less
5. The Stellar Mass–Black Hole Mass Relation at z ∼ 2 down to ${}_{\mathrm{BH}}\sim {10}^7{M}_{\odot }$ Determined by HETDEX

   https://doi.org/10.3847/1538-4357/acc2c2  

   Zhang, Yechi; Ouchi, Masami; Gebhardt, Karl; Liu, Chenxu; Harikane, Yuichi; Cooper, Erin_Mentuch; Davis, Dustin; Farrow, Daniel_J; Gawiser, Eric; Hill, Gary_J; et al (May 2023, The Astrophysical Journal)

   Abstract We investigate the stellar mass–black hole mass ( ${\mathit{}}_{*}–{\mathit{}}_{\mathrm{BH}}$) relation with type 1 active galactic nuclei (AGNs) down to ${\mathit{}}_{\mathrm{BH}}={10}^7{M}_{\odot }$, corresponding to a ≃ −21 absolute magnitude in rest-frame ultraviolet, atz= 2–2.5. Exploiting the deep and large-area spectroscopic survey of the Hobby–Eberly Telescope Dark Energy Experiment (HETDEX), we identify 66 type 1 AGNs with ${\mathit{}}_{\mathrm{BH}}$ranging from 10⁷–10¹⁰M_⊙that are measured with single-epoch virial method using Civemission lines detected in the HETDEX spectra. ${\mathit{}}_{*}$of the host galaxies are estimated from optical to near-infrared photometric data taken with Spitzer, the Wide-field Infrared Survey Explorer, and ground-based 4–8 m class telescopes byCIGALEspectral energy distribution (SED) fitting. We further assess the validity of SED fitting in two cases by host-nuclear decomposition performed through surface brightness profile fitting on spatially resolved host galaxies with the James Webb Space Telescope/NIRCam CEERS data. We obtain the ${\mathit{}}_{*}–{\mathit{}}_{\mathrm{BH}}$relation covering the unexplored low-mass ranges of ${\mathit{}}_{\mathrm{BH}}\sim {10}^7–{10}^8{M}_{\odot }$, and conduct forward modeling to fully account for the selection biases and observational uncertainties. The intrinsic ${\mathit{}}_{*}–{\mathit{}}_{\mathrm{BH}}$relation atz∼ 2 has a moderate positive offset of 0.52 ± 0.14 dex from the local relation, suggestive of more efficient black hole growth at higher redshift even in the low-mass regime of ${\mathit{}}_{\mathrm{BH}}\sim {10}^7–{10}^8{M}_{\odot }$. Our ${\mathit{}}_{*}–{\mathit{}}_{\mathrm{BH}}$relation is inconsistent with the ${\mathit{}}_{\mathrm{BH}}$suppression at the low- ${\mathit{}}_{*}$regime predicted by recent hydrodynamic simulations at a 98% confidence level, suggesting that feedback in the low-mass systems may be weaker than those produced in hydrodynamic simulations. 

   more » « less