skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Length and Velocity Scales in Protoplanetary Disk Turbulence
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 α H and α 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 α / R H and α R c , where R Ω L / Ω K is twice the Rossby number. It follows that α = α ˜ / R , where α ˜ 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 > 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
Award ID(s):
2007422
PAR ID:
10554145
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
The Astrophysical Journal
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
966
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
90
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We measure the CO-to-H2conversion factor (αCO) in 37 galaxies at 2 kpc resolution, using the dust surface density inferred from far-infrared emission as a tracer of the gas surface density and assuming a constant dust-to-metal ratio. In total, we have ∼790 and ∼610 independent measurements ofαCOfor CO (2–1) and (1–0), respectively. The mean values forαCO (2–1)andαCO (1–0)are 9.3 5.4 + 4.6 and 4.2 2.0 + 1.9 M pc 2 ( K km s 1 ) 1 , respectively. The CO-intensity-weighted mean is 5.69 forαCO (2–1)and 3.33 forαCO (1–0). We examine howαCOscales with several physical quantities, e.g., the star formation rate (SFR), stellar mass, and dust-mass-weighted average interstellar radiation field strength ( U ¯ ). Among them, U ¯ , ΣSFR, and the integrated CO intensity (WCO) have the strongest anticorrelation with spatially resolvedαCO. We provide linear regression results toαCOfor all quantities tested. At galaxy-integrated scales, we observe significant correlations betweenαCOandWCO, metallicity, U ¯ , and ΣSFR. We also find thatαCOin each galaxy decreases with the stellar mass surface density (Σ) in high-surface-density regions (Σ≥ 100Mpc−2), following the power-law relations α CO ( 2 1 ) Σ 0.5 and α CO ( 1 0 ) Σ 0.2 . The power-law index is insensitive to the assumed dust-to-metal ratio. We interpret the decrease inαCOwith increasing Σas a result of higher velocity dispersion compared to isolated, self-gravitating clouds due to the additional gravitational force from stellar sources, which leads to the reduction inαCO. The decrease inαCOat high Σis important for accurately assessing molecular gas content and star formation efficiency in the centers of galaxies, which bridge “Milky Way–like” to “starburst-like” conversion factors. 
    more » « less
  2. Abstract Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral atz> 7 and largely ionized byz∼ 5.5. However, most methods to measure the IGM neutral fraction are highly model dependent and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low ( x ¯ H I 10 3 ) or close to unity ( x ¯ H I 1 ). In particular, the neutral fraction evolution of the IGM at the critical redshift range ofz= 6–7 is poorly constrained. We present new constraints on x ¯ H I atz∼ 5.1–6.8 by analyzing deep optical spectra of 53 quasars at 5.73 <z< 7.09. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of “dark” pixels identified in the Lyαand Lyβforests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction atz∼ 6.2–6.8 using quasar absorption measurements. Our results give upper limits of x ¯ H I ( z = 6.3 ) < 0.79 ± 0.04 (1σ), x ¯ H I ( z = 6.5 ) < 0.87 ± 0.03 (1σ), and x ¯ H I ( z = 6.7 ) < 0.94 0.09 + 0.06 (1σ). The dark pixel fractions atz> 6.1 are consistent with the redshift evolution of the neutral fraction of the IGM derived from Planck 2018. 
    more » « less
  3. Abstract State transitions in black hole X-ray binaries are likely caused by gas evaporation from a thin accretion disk into a hot corona. We present a height-integrated version of this process, which is suitable for analytical and numerical studies. With radiusrscaled to Schwarzschild units and coronal mass accretion rate m ̇ c to Eddington units, the results of the model are independent of black hole mass. State transitions should thus be similar in X-ray binaries and an active galactic nucleus. The corona solution consists of two power-law segments separated at a break radiusrb∼ 103(α/0.3)−2, whereαis the viscosity parameter. Gas evaporates from the disk to the corona forr>rb, and condenses back forr<rb. Atrb, m ̇ c reaches its maximum, m ̇ c , max 0.02 ( α / 0.3 ) 3 . If atr≫rbthe thin disk accretes with m ̇ d < m ̇ c , max , then the disk evaporates fully before reachingrb, giving the hard state. Otherwise, the disk survives at all radii, giving the thermal state. While the basic model considers only bremsstrahlung cooling and viscous heating, we also discuss a more realistic model that includes Compton cooling and direct coronal heating by energy transport from the disk. Solutions are again independent of black hole mass, andrbremains unchanged. This model predicts strong coronal winds forr>rb, and aT∼ 5 × 108K Compton-cooled corona forr<rb. Two-temperature effects are ignored, but may be important at small radii. 
    more » « less
  4. Abstract The growth of supermassive black holes is strongly linked to their galaxies. It has been shown that the population mean black hole accretion rate ( BHAR ¯ ) primarily correlates with the galaxy stellar mass (M) and redshift for the general galaxy population. This work aims to provide the best measurements of BHAR ¯ as a function ofMand redshift over ranges of 109.5<M< 1012Mandz< 4. We compile an unprecedentedly large sample with 8000 active galactic nuclei (AGNs) and 1.3 million normal galaxies from nine high-quality survey fields following a wedding cake design. We further develop a semiparametric Bayesian method that can reasonably estimate BHAR ¯ and the corresponding uncertainties, even for sparsely populated regions in the parameter space. BHAR ¯ is constrained by X-ray surveys sampling the AGN accretion power and UV-to-infrared multiwavelength surveys sampling the galaxy population. Our results can independently predict the X-ray luminosity function (XLF) from the galaxy stellar mass function (SMF), and the prediction is consistent with the observed XLF. We also try adding external constraints from the observed SMF and XLF. We further measure BHAR ¯ for star-forming and quiescent galaxies and show that star-forming BHAR ¯ is generally larger than or at least comparable to the quiescent BHAR ¯
    more » « less
  5. 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 functionf0(r,p) at cylindrical radiusrfrom the jet axis (assumed to lie along thez-axis) is given by convolving the particle momentum spectrum f 0 ( , p ) with the Green’s function G ( r , p ; p ) , which describes the monoenergetic spectrum solution in which f 0 δ ( p p ) asr→ ∞ . Previous work by Webb et al. studied only the Green’s function solution for G ( r , p ; p ) . In this paper, we explore for the first time, solutions for more general and realistic forms for f 0 ( , p ) . The flow velocityu=u(r)ezis along the axis of the jet (thez-axis).uis independent ofz, andu(r) is a monotonic decreasing function ofr. The scattering time τ ( r , p ) = τ 0 ( p / p 0 ) α in the shear flow region 0 <r<r2, and τ ( r , p ) = τ 0 ( p / p 0 ) α ( r / r 2 ) s , wheres> 0 in the regionr>r2is 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 ψ p ( r , p ; p ) that particles observed at (r,p) originated fromr→ ∞ with momentum p . 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