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.


This content will become publicly available on January 27, 2026

Title: Extreme Scale Height Variations and Nozzle Shocks in Warped Disks
Abstract Accretion disks around both stellar-mass and supermassive black holes (BHs) are likely often warped. Whenever a disk is warped, its scale height varies with azimuth. Sufficiently strong warps cause extreme compressions of the scale height, which fluid parcels “bounce” off of twice per orbit to high latitudes. We study the dynamics of strong warps using (i) the nearly analytic “ring theory” of Fairbairn & Ogilvie, which we generalize to the Kerr metric, and (ii) three-dimensional general-relativistic hydrodynamic simulations of tori (“rings”) around BHs, using theH-AMRcode. We initialize a ring with a warp and study its evolution on tens of orbital periods. The simulations agree excellently with the ring theory until the warp amplitude,ψ, reaches a critical valueψc. Whenψ > ψc, the rings enter the bouncing regime. We analytically derive (and numerically validate) that ψ c ( r / r g ) 1 / 2 in the non-Keplerian regime, whererg = GM/c2is the gravitational radius, andMis the mass of the central object. Whenever the scale height bounces, the vertical velocity becomes supersonic, leading to “nozzle shocks” as gas collides at the scale height minima. Nozzle shocks damp the warp within ≈10–20 orbits, which is not captured by the ring theory. Nozzle shock dissipation leads to inflow timescales 1–2 orders of magnitude shorter than unwarpedαdisks, which may result in rapid variability, such as in changing-look active galactic nuclei or in the soft state of X-ray binaries. We propose that steady disks with strong warps may self-regulate to have amplitudes nearψc more » « less
Award ID(s):
2407475
PAR ID:
10589503
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
IOP Publishing
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
979
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
192
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract It remains unclear what mechanism is driving the evolution of protoplanetary disks. Direct detection of the main candidates, either turbulence driven by magnetorotational instabilities or magnetohydrodynamical disk winds, has proven difficult, leaving the time evolution of the disk size as one of the most promising observables able to differentiate between these two mechanisms. But to do so successfully, we need to understand what the observed gas disk size actually traces. We studied the relation betweenRCO,90%, the radius that encloses 90% of the12CO flux, andRc, the radius that encodes the physical disk size, in order to provide simple prescriptions for conversions between these two sizes. For an extensive grid of thermochemical models, we calculateRCO,90%from synthetic observations and relate properties measured at this radius, such as the gas column density, to bulk disk properties, such asRcand the disk massMdisk. We found an empirical correlation between the gas column density atRCO,90%and disk mass: N gas ( R CO , 90 % ) 3.73 × 10 21 ( M disk / M ) 0.34 cm 2 . Using this correlation we derive an analytical prescription ofRCO,90%that only depends onRcandMdisk. We deriveRcfor disks in Lupus, Upper Sco, Taurus, and the DSHARP sample, finding that disks in the older Upper Sco region are significantly smaller (〈Rc〉 = 4.8 au) than disks in the younger Lupus and Taurus regions (〈Rc〉 = 19.8 and 20.9 au, respectively). This temporal decrease inRcgoes against predictions of both viscous and wind-driven evolution, but could be a sign of significant external photoevaporation truncating disks in Upper Sco. 
    more » « less
  2. 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
  3. A measurement of the ratio of branching fractions R ( J / ψ ) = B ( B c + J / ψ τ + ν τ ) / B ( B c + J / ψ μ + ν μ ) in the J / ψ μ + μ , τ + μ + ν μ ν ¯ τ decay channel is presented. This measurement uses a sample of proton-proton collision data collected at a center-of-mass energy of 13 TeV by the CMS experiment in 2018, corresponding to an integrated luminosity of 59.7 fb 1 . The measured ratio, R ( J / ψ ) = 0.1 7 0.17 + 0.18 ( stat ) 0.22 + 0.21 ( syst ) 0.18 + 0.19 ( theo ) = 0.17 ± 0.33 , agrees with the value of 0.2582 ± 0.0038 predicted by the standard model, which assumes lepton flavor universality. By testing lepton flavor universality, this measurement is a probe of new physics using B c + mesons, which are currently only produced at the LHC. 
    more » « less
  4. Abstract We study the linear stability of a planar interface separating two fluids in relative motion, focusing on conditions appropriate for the boundaries of relativistic jets. The jet is magnetically dominated, whereas the ambient wind is gas-pressure-dominated. We derive the most general form of the dispersion relation and provide an analytical approximation of its solution for an ambient sound speed much smaller than the jet Alfvén speedvA, as appropriate for realistic systems. The stability properties are chiefly determined by the angleψbetween the wavevector and the jet magnetic field. Forψ=π/2, magnetic tension plays no role, and our solution resembles the one of a gas-pressure-dominated jet. Here, only sub-Alfvénic jets are unstable ( 0 < M e ( v / v A ) cos θ < 1 , wherevis the shear velocity andθthe angle between the velocity and the wavevector). Forψ= 0, the free energy in the velocity shear needs to overcome the magnetic tension, and only super-Alfvénic jets are unstable ( 1 < M e < ( 1 + Γ w 2 ) / [ 1 + ( v A / c ) 2 Γ w 2 ] , with Γwthe wind adiabatic index). Our results have important implications for the propagation and emission of relativistic magnetized jets. 
    more » « less
  5. The production of the ψ ( 2 S ) charmonium state was measured with ALICE in Pb-Pb collisions at s NN = 5.02 TeV , in the dimuon decay channel. A significant signal was observed for the first time at LHC energies down to zero transverse momentum, at forward rapidity ( 2.5 < y < 4 ). The measurement of the ratio of the inclusive production cross sections of the ψ ( 2 S ) and J / ψ resonances is reported as a function of the centrality of the collisions and of transverse momentum, in the region p T < 12 GeV / c . The results are compared with the corresponding measurements in p p collisions, by forming the double ratio [ σ ψ ( 2 S ) / σ J / ψ ] Pb Pb / [ σ ψ ( 2 S ) / σ J / ψ ] p p . It is found that in Pb-Pb collisions the ψ ( 2 S ) is suppressed by a factor of 2 with respect to the J / ψ . The ψ ( 2 S ) nuclear modification factor R AA was also obtained as a function of both centrality and p T . The results show that the ψ ( 2 S ) resonance yield is strongly suppressed in Pb-Pb collisions, by a factor of up to 3 with respect to p p . Comparisons of cross section ratios with previous Super Proton Synchrotron findings by the NA50 experiment and of R AA with higher- p T results at LHC energy are also reported. These results and the corresponding comparisons with calculations of transport and statistical models address questions on the presence and properties of charmonium states in the quark-gluon plasma formed in nuclear collisions at the LHC. © 2024 CERN, for the ALICE Collaboration2024CERN 
    more » « less