Search NASA⌕ Search

SEARCH · Search NASA

Results for “Inclination”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evidence of a Dwarf Galaxy Stream Populating the Inner Milky Way Halo

Stellar streams produced from dwarf galaxies provide direct evidence of the hierarchical formation of the Milky Way. Here, we present the first comprehensive study of the LMS-1 stellar stream, that we detect by searching for wide streams in the Gaia EDR3 data set using the STREAMFINDER algorithm. This stream was recently discovered by Yuan et al. We detect LMS-1 as a 60° long stream to the north of the galactic bulge, at a distance of ~20 kpc from the Sun, together with additional components that suggest that the overall stream is completely wrapped around the inner Galaxy. Using spectroscopic measurements from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, the Sloan Digital Sky Survey, and the Apache Point Observatory Galactic Evolution Experiment, we infer that the stream is very metal-poor (<[Fe/H]> = –2.1) with a significant metallicity dispersion (σ [Fe/H] = 0.4), and it possesses a large radial velocity dispersion (σ v = 20 ± 4 km s–1 ). These estimates together imply that LMS-1 is a dwarf galaxy stream. The orbit of LMS-1 is close to polar, with an inclination of 75° to the galactic plane. Both the orbit and metallicity of LMS-1 are remarkably similar to the globular clusters NGC 5053, NGC 5024, and the stellar stream Indus. As a result, these findings make LMS-1 an important contributor to the stellar population of the inner Milky Way halo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

SDSS-IV MaNGA: Understanding Ionized Gas Turbulence Using Integral Field Spectroscopy of 4500 Star-forming Disk Galaxies

The Sloan Digital Sky Survey MaNGA program has now obtained integral field spectroscopy for over 10,000 galaxies in the nearby universe. We use the final MaNGA data release DR17 to study the correlation between ionized gas velocity dispersion and galactic star formation rate, finding a tight correlation in which σ Hα from galactic H ιι regions increases significantly from ~18–30 km s –1 , broadly in keeping with previous studies. In contrast, σ Hα from diffuse ionized gas increases more rapidly from 20–60 km s –1 . Using the statistical power of MaNGA, we investigate these correlations in greater detail using multiple emission lines and determine that the observed correlation of σ Hα with local star formation rate surface density is driven primarily by the global relation of increasing velocity dispersion at higher total star formation rate, as are apparent correlations with stellar mass. Assuming H ιι region models consistent with our finding that σ [O III] < σ Hα < σ [O I] , we estimate the velocity dispersion of the molecular gas in which the individual H ιι regions are embedded, finding values σ Mol = 5–30 km s –1 consistent with ALMA observations in a similar mass range. Finally, we use variations in the relation with inclination and disk azimuthal angle to constrain the velocity dispersion ellipsoid of the ionized gas σ z /σ r = 0.84 ± 0.03 and σ Φ /σ r = 0.91 ± 0.03, similar to that of young stars in the Galactic disk. Our results are most consistent with the theoretical models in which turbulence in modern galactic disks is driven primarily by star formation feedback.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Lick AGN Monitoring Project 2016: Dynamical Modeling of Velocity-resolved H β Lags in Luminous Seyfert Galaxies

We have modeled the velocity-resolved reverberation response of the Hβ broad emission line in nine Seyfert 1 galaxies from the Lick Active Galactic Nucleus (AGN) Monitoring Project 2016 sample, drawing inferences on the geometry and structure of the low-ionization broad-line region (BLR) and the mass of the central supermassive black hole. Overall, we find that the Hβ BLR is generally a thick disk viewed at low to moderate inclination angles. We combine our sample with prior studies and investigate line-profile shape dependence, such as ${\mathrm{log}}_{10}(\mathrm{FWHM}/\sigma )$, on BLR structure and kinematics and search for any BLR luminosity-dependent trends. We find marginal evidence for an anticorrelation between the profile shape of the broad Hβ emission line and the Eddington ratio, when using the rms spectrum. However, we do not find any luminosity-dependent trends, and conclude that AGNs have diverse BLR structure and kinematics, consistent with the hypothesis of transient AGN/BLR conditions rather than systematic trends.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Long-term Evolution of Tightly Packed Stellar Black Holes in AGN Disks: Formation of Merging Black Hole Binaries via Close Encounters

We study the long-term evolution of two or more stellar black holes (BHs) on initially separated but unstable circular orbits around a supermassive BH (SMBH). Such a close-packed orbital configuration can naturally arise from BH migrations in the AGN disk. Dynamical instability of the orbits leads to recurring close encounters between two BHs, during which the BH separation r p becomes less than the Hill radius R H . In rare very close encounters, a tight merging BH binary can form with the help of gravitational wave emission. We use N-body simulations to study the time evolution of close encounters of various degrees of closeness. For a typical "SMBH+2BH" system, the averaged cumulative number of close encounters (with r p ≲ R H ) scales approximately as ∝ t 0.5 . The minimum encounter separation r p follows a cumulative distribution P( p ) ∝ r p for r p « R H . We obtain a semi-analytical expression for the averaged rate of binary captures that lead to BH mergers. Our results suggest that close-packed BHs in AGN disks may take a long time (≳10 7 orbits around the SMBH) to experience a sufficiently close encounter and form a bound binary. This time can be shorter if the initial BH orbits are highly aligned. The BH binary mergers produced in this scenario have high eccentricities when entering the LIGO band and broad distribution of orbital inclinations relative to the original AGN disk. We explore the effects of the gas disk and find that simple gas drags on the BHs do not necessarily lead to an enhanced BH binary capture rate.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Physical Properties of Low-redshift FeLoBAL Quasars. II. The Rest-frame Optical Emission Line Properties

We report the results of an analysis of the Hβ emission line region of a sample of 30 low-redshift (z < 1) iron low-ionization broad absorption line quasars (FeLoBALQs). Eleven of these objects are newly classified as FeLoBALQs. A matched sample of 132 unabsorbed quasars was analyzed in parallel. The emission lines showed the well-known anticorrelation between the [O III ] and Fe II emission. Using a summary statistic called E1 to quantify this anticorrelation, we found that while the distribution of E1 for the unabsorbed quasars has a single peak, the FeLoBALQs have a bimodal shape in this parameter. Previous studies have shown that the line emission properties of BAL and non-BALQs are consistent; therefore, the difference in the Hβ region emission between FeLoBALQs and unabsorbed quasars is a new result. The two populations of FeLoBALQs are characterized by low and high bolometric luminosities and Eddington ratios. Some previous studies have suggested that BALQs are high accretion rate objects and therefore the discovery of the low accretion rate branch of FeLoBAL quasars was unexpected. We also found that the Hβ FWHM is systematically broader among the FeLoBALQs, implying a higher inclination viewing angle or a dearth of low velocity line emitting gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

BASS XXXII: Studying the Nuclear Millimeter-wave Continuum Emission of AGNs with ALMA at Scales ≲100–200 pc

To understand the origin of nuclear (≲100 pc) millimeter-wave (mm-wave) continuum emission in active galactic nuclei (AGNs), we systematically analyzed subarcsecond resolution Band-6 (211–275 GHz) Atacama Large Millimeter/submillimeter Array data of 98 nearby AGNs (z < 0.05) from the 70 month Swift/BAT catalog. The sample, almost unbiased for obscured systems, provides the largest number of AGNs to date with high mm-wave spatial resolution sampling (~1–200 pc), and spans broad ranges of 14–150 keV luminosity {$40\lt \mathrm{log}[{L}_{14-150}/(\mathrm{erg}\,{{\rm{s}}}^{-1})]\lt 45$}, black hole mass $[5\lt \mathrm{log}({M}_{\mathrm{BH}}/{M}_{\odot })\lt 10$], and Eddington ratio ($-4\lt \mathrm{log}{\lambda }_{\mathrm{Edd}}\lt 2$). We find a significant correlation between 1.3 mm (230 GHz) and 14–150 keV luminosities. Its scatter is ≈0.36 dex, and the mm-wave emission may serve as a good proxy of the AGN luminosity, free of dust extinction up to N H ~ 10 26 cm –2 . While the mm-wave emission could be self-absorbed synchrotron radiation around the X-ray corona according to past works, we also discuss different possible origins of the mm-wave emission: AGN-related dust emission, outflow-driven shocks, and a small-scale (<200 pc) jet. The dust emission is unlikely to be dominant, as the mm-wave slope is generally flatter than expected. Also, due to no increase in the mm-wave luminosity with the Eddington ratio, a radiation-driven outflow model is possibly not the common mechanism. Furthermore, we find independence of the mm-wave luminosity on indicators of the inclination angle from the polar axis of the nuclear structure, which is inconsistent with a jet model whose luminosity depends only on the angle.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Magnetic Energy Dissipation and γ-Ray Emission in Energetic Pulsars

Abstract Some of the most energetic pulsars exhibit rotation-modulated γ -ray emission in the 0.1–100 GeV band. The luminosity of this emission is typically 0.1%–10% of the pulsar spin-down power ( γ -ray efficiency), implying that a significant fraction of the available electromagnetic energy is dissipated in the magnetosphere and reradiated as high-energy photons. To investigate this phenomenon we model a pulsar magnetosphere using 3D particle-in-cell simulations with strong synchrotron cooling. We particularly focus on the dynamics of the equatorial current sheet where magnetic reconnection and energy dissipation take place. Our simulations demonstrate that a fraction of the spin-down power dissipated in the magnetospheric current sheet is controlled by the rate of magnetic reconnection at microphysical plasma scales and only depends on the pulsar inclination angle. We demonstrate that the maximum energy and the distribution function of accelerated pairs is controlled by the available magnetic energy per particle near the current sheet, the magnetization parameter. The shape and the extent of the plasma distribution is imprinted in the observed synchrotron emission, in particular, in the peak and the cutoff of the observed spectrum. We study how the strength of synchrotron cooling affects the observed variety of spectral shapes. Our conclusions naturally explain why pulsars with higher spin-down power have wider spectral shapes and, as a result, lower γ -ray efficiency.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Chemodynamical Analysis of Metal-rich High-eccentricity Stars in the Milky Way's Disk

We present a chemodynamical analysis of 11,562 metal-rich, high-eccentricity halo-like main-sequence stars, which have been referred to as the Splash or Splashed Disk, selected from the Sloan Digital Sky Survey and Large Sky Area Multi-Object Fiber Spectroscopic Telescope. When divided into two groups, a low-[α/Fe] population (LAP) and a high-[α/Fe] population (HAP), based on kinematics and chemistry, we find that they exhibit very distinct properties, indicative of different origins. From a detailed analysis of their orbital inclinations, we suggest that the HAP arises from a large fraction (~90%) of heated disk stars and a small fraction (~10%) of in situ stars from a starburst population, likely induced by interaction of the Milky Way with the Gaia-Sausage/Enceladus (GSE) or another early merger. The LAP comprises about half accreted stars from the GSE and half formed by the GSE-induced starburst. Our findings further imply that the Splash stars in our sample originated from at least three different mechanisms: accretion, disk heating, and a merger-induced starburst.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Unraveling Joint Evolution of Bars, Star Formation, and Active Galactic Nuclei of Disk Galaxies

We aim to unravel the interplay between bars, star formation (SF), and active galactic nuclei (AGNs) in barred galaxies. To this end, we utilize the SDSS DR12 to select a sample of nearby (0.02 < z < 0.06) disk galaxies that are suitable for bar examination (M r < –20.12 and inclination ≲53°). We identify 3662 barred galaxies and measure the length and axis ratio of each bar. We invent new bar parameters that mitigate the stellar and bulge mass biases and show, for the first time, that the evolution of non-AGN and AGN-hosting barred galaxies should be tracked using different bar parameters: the bar length for non-AGN galaxies and the bar axis ratio for AGN-hosting galaxies. Our analysis confirms that barred galaxies have a higher specific SF rate than unbarred control galaxies. Moreover, we find a positive correlation of bar length with both the SF enhancement and the centrally star-forming galaxy fraction, indicating the interconnectivity of bars and SF through the bar-driven gas inflow. We also find that, while the AGN fraction of barred galaxies is the same as that of the unbarred control sample, galaxies hosting more massive black holes (BHs) have rounder (i.e., higher axis ratio) bars, implying that the bar is not a cause of AGN activity; rather, AGNs appear to regulate bars. Our findings corroborate theoretical predictions that bars in non-AGN galaxies grow in length, and bars in AGN-hosting galaxies become rounder as BHs grow and eventually get destroyed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Distant Echoes of the Milky Way’s Last Major Merger

The majority of the Milky Way’s stellar halo consists of debris from our galaxy’s last major merger, the Gaia-Sausage-Enceladus (GSE). In the past few years, stars from the GSE have been kinematically and chemically studied in the inner 30 kpc of our galaxy. However, simulations predict that accreted debris could lie at greater distances, forming substructures in the outer halo. Here we derive metallicities and distances using Gaia DR3 XP spectra for an all-sky sample of luminous red giant stars, and map the outer halo with kinematics and metallicities out to 100 kpc. We obtain follow-up spectra of stars in two strong overdensities—including the previously identified outer Virgo Overdensity—and find them to be relatively metal rich and on predominantly retrograde orbits, matching predictions from simulations of the GSE merger. We argue that these are apocentric shells of GSE debris, forming 60–90 kpc counterparts to the 15–20 kpc shells that are known to dominate the inner stellar halo. Extending our search across the sky with literature radial velocities, we find evidence for a coherent stream of retrograde stars encircling the Milky Way from 50 to 100 kpc, in the same plane as the Sagittarius Stream but moving in the opposite direction. These are the first discoveries of distant and structured imprints from the GSE merger, cementing the picture of an inclined and retrograde collision that built up our galaxy’s stellar halo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Latent Stochastic Differential Equations for Modeling Quasar Variability and Inferring Black Hole Properties

Quasars are bright and unobscured active galactic nuclei (AGN) thought to be powered by the accretion of matter around supermassive black holes at the centers of galaxies. The temporal variability of a quasar’s brightness contains valuable information about its physical properties. The UV/optical variability is thought to be a stochastic process, often represented as a damped random walk described by a stochastic differential equation (SDE). Upcoming wide-field telescopes such as the Rubin Observatory Legacy Survey of Space and Time (LSST) are expected to observe tens of millions of AGN in multiple filters over a ten year period, so there is a need for efficient and automated modeling techniques that can handle the large volume of data. Latent SDEs are machine learning models well suited for modeling quasar variability, as they can explicitly capture the underlying stochastic dynamics. In this work, we adapt latent SDEs to jointly reconstruct multivariate quasar light curves and infer their physical properties such as the black hole mass, inclination angle, and temperature slope. Our model is trained on realistic simulations of LSST ten year quasar light curves, and we demonstrate its ability to reconstruct quasar light curves even in the presence of long seasonal gaps and irregular sampling across different bands, outperforming a multioutput Gaussian process regression baseline. Our method has the potential to provide a deeper understanding of the physical properties of quasars and is applicable to a wide range of other multivariate time series with missing data and irregular sampling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Study of Magnetic Field and Turbulence in the TeV Halo around the Monogem Pulsar

Magnetic fields are ubiquitous in the interstellar medium, including extended objects such as supernova remnants and diffuse halos around pulsars. Its turbulent characteristics govern the diffusion of cosmic rays and the multiwavelength emission from pulsar wind nebulae (PWNe). However, the geometry and turbulence nature of the magnetic fields in the ambient region of PWN is still unknown. Recent gamma-ray observations from HAWC and synchrotron observations suggest a highly suppressed diffusion coefficient compared to the mean interstellar value. In this study, we present the first direct observational evidence of the orientation of the mean magnetic field and turbulent characteristics by employing a recently developed statistical parameter "Y turb " in the extended halo around the Monogem pulsar. Our study points to two possible scenarios: nearly aligned toward the line of sight (LOS) with compressible modes dominance or high inclination angle toward the LOS and characterized by Alfvénic turbulence. The first scenario appears consistent with other observational signatures. Furthermore, we report that the magnetic field has an observed correlation length of approximately 3 ± 0.6 pc in the Monogem halo. Our study highlights the pivotal role of magnetic field and turbulence in unraveling the physical processes in TeV halos and cosmic-ray transport.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Slant, Fan, and Narrow: The Response of Stellar Streams to a Tilting Galactic Disk

Stellar streams are sensitive tracers of the gravitational potential, which is typically assumed to be static in the inner Galaxy. However, massive mergers like Gaia-Sausage-Enceladus can impart torques on the stellar disk of the Milky Way that result in the disk tilting at rates of up to 10°–20° Gyr –1 . Here, we demonstrate the effects of disk tilting on the morphology and kinematics of stellar streams. Through a series of numerical experiments, we find that streams with nearby apocenters (r apo ≲ 20 kpc) are sensitive to disk tilting, with the primary effect being changes to the stream's on-sky track and width. Interestingly, disk tilting can produce both more diffuse streams and more narrow streams, depending on the orbital inclination of the progenitor and the direction in which the disk is tilting. Our model of Pal 5's tidal tails for a tilting rate of 15° Gyr –1 is in excellent agreement with the observed stream's track and width, and reproduces the extreme narrowing of the trailing tail. We also find that failure to account for a tilting disk can bias constraints on shape parameters of the Milky Way's local dark matter distribution at the level of 5%–10%, with the direction of the bias changing for different streams. Disk tilting could therefore explain discrepancies in the Milky Way's dark matter halo shape inferred using different streams.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Distance Estimate Method for Asymptotic Giant Branch Stars Using Infrared Spectral Energy Distributions

We present a method to estimate distances to asymptotic giant branch (AGB) stars in the Galaxy, using spectral energy distributions (SEDs) in the near- and mid-infrared. By assuming that a given set of source properties (initial mass, stellar temperature, composition, and evolutionary stage) will provide a typical SED shape and brightness, sources are color matched to a distance-calibrated template and thereafter scaled to extract the distance. The method is tested by comparing the distances obtained to those estimated from very long baseline interferometry or Gaia parallax measurements, yielding a strong correlation in both cases. Additional templates are formed by constructing a source sample likely to be close to the Galactic center, and thus with a common, typical distance for calibration of the templates. These first results provide statistical distance estimates to a set of almost 15,000 Milky Way AGB stars belonging to the Bulge Asymmetries and Dynamical Evolution (BAaDE) survey, with typical distance errors of ±35%. With these statistical distances, a map of the intermediate-age population of stars traced by AGBs is formed, and a clear bar structure can be discerned, consistent with the previously reported inclination angle of 30° to the GC–Sun direction vector. These results motivate deeper studies of the AGB population to tease out the intermediate-age stellar distribution throughout the Galaxy, as well as determining statistical properties of the AGB population luminosity and mass-loss-rate distributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Multifrequency Analysis of Favored Models for the Messier 87* Accretion Flow

Abstract The polarized images of the supermassive black hole Messier 87* (M87*) produced by the Event Horizon Telescope (EHT) provide a direct view of the near-horizon emission from a black hole accretion and jet system. The EHT theoretical analysis of the polarized M87* images compared thousands of snapshots from numerical models with a variety of spins, magnetization states, viewing inclinations, and electron energy distributions, and found a small subset consistent with the observed image. In this article, we examine two models favored by EHT analyses: a magnetically arrested disk with moderate retrograde spin and a magnetically arrested disk with high prograde spin. Both have electron distribution functions that lead to strong depolarization by cold electrons. We ray trace five snapshots from each model at 22, 43, 86, 230, 345, and 690 GHz to forecast future very long baseline interferometry (VLBI) observations and examine limitations in numerical models. We find that even at low frequencies where optical and Faraday rotation depths are large, approximately rotationally symmetric polarization persists, suggesting that shallow depths dominate the polarization signal. However, morphology and spectra suggest that the assumed thermal electron distribution is not adequate to describe emission from the jet. We find 86 GHz images show a ringlike shape determined by a combination of plasma and spacetime imprints, smaller in diameter than recent results from the Global mm-VLBI Array. We find that the photon ring becomes more apparent with increasing frequency, and is more apparent in the retrograde model, leading to large differences between models in asymmetry and polarization structure.

Astronomy & Astrophysics↗
DOE OSTI2024

Spindown of Pulsars Interacting with Companion Winds: The Impact of Magnetospheric Compression

Pulsars in binary systems with strong companion winds can have the magnetopause separating their magnetosphere from the wind located well within their light cylinder. This bow-like enclosure effectively creates a waveguide that confines the pulsar's electromagnetic fields and can significantly alter its spindown. In this paper, we study the spindown of compressed pulsar magnetospheres in such systems. We parameterize the confinement as the ratio between the equatorial position of the magnetopause (or standoff distance) R m and the pulsar's light cylinder R LC . Using particle-in-cell simulations, we quantify the pulsar spindown for a range of compressions, R m /R LC = 1/3–1, and inclination angles, χ = 0°...90°, between magnetic and rotation axes. Our strongly confined models (R m /R LC = 1/3) show two distinct limits. For χ = 0°, the spindown of a compressed pulsar magnetosphere is enhanced by approximately a factor of three compared to an isolated pulsar due to the increased number of open magnetic field lines. Conversely, for χ = 90°, the compressed pulsar spins down at less than 40% of the rate of an isolated reference pulsar due to the mismatch between the pulsar wind stripe wavelength and the waveguide size. We apply our analysis to the 2.77 s oblique rotator (χ = 60°) in the double-pulsar system PSR J0737-3039. With the numerically derived spindown estimate, we constrain its surface magnetic field to B* ≈ (7.3 ± 0.2) × 10 11 G. We discuss the time modulation of its period derivative, the effects of compression on its braking index, and implications for the radio eclipse in PSR J0737-3039.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Applying the Velocity Gradient Technique in NGC 1333: Comparison with Dust Polarization Observations

Magnetic fields (B-fields) are ubiquitous in the interstellar medium (ISM), and they play an essential role in the formation of molecular clouds and subsequent star formation. However, B-fields in interstellar environments remain challenging to measure, and their properties typically need to be inferred from dust polarization observations over multiple physical scales. In this work, we seek to use a recently proposed approach called the velocity gradient technique (VGT) to study B-fields in star-forming regions and compare the results with dust polarization observations in different wavelengths. The VGT is based on the anisotropic properties of eddies in magnetized turbulence to derive B-field properties in the ISM. We investigate that this technique is synergistic with dust polarimetry when applied to a turbulent diffused medium for the purpose of measuring its magnetization. Specifically, we use the VGT on molecular line data toward the NGC 1333 star-forming region ( 12 CO, 13 CO, C 18 O, and N 2 H + ), and we compare the derived B-field properties with those inferred from 214 and 850 μm dust polarization observations of the region using Stratospheric Observatory for Infrared Astronomy/High-Resolution Airborne Wide-band Camera Plus and James Clerk Maxwell Telescope/POL-2, respectively. We estimate both the inclination angle and the 3D Alfvénic Mach number M A from the molecular line gradients. Crucially, testing this technique on gravitationally bound, dynamic, and turbulent regions, and comparing the results with those obtained from polarization observations at different wavelengths, such as the plane-of-sky field orientation, is an important test on the applicability of the VGT in various density regimes of the ISM. We in general do not find a close correlation between the velocity gradient inferred orientations and the dust inferred magnetic field orientations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Jipole: A Differentiable ipole-based Code for Radiative Transfer in Curved Spacetimes

Recent imaging of supermassive black holes by the Event Horizon Telescope has relied on exhaustive parameter-space searches, matching observations to large, precomputed libraries of theoretical models. As observational data become increasingly precise, the limitations of this computationally expensive approach grow more acute, creating a pressing need for more efficient methods. In this work, we present Jipole, an automatically differentiable (AD), ipole-based code for radiative transfer in curved spacetimes, designed to compute image gradients with respect to underlying model parameters. These gradients quantify how parameter changes—such as the black hole’s spin or the observer’s inclination—affect the image, enabling more efficient parameter estimation and reducing the number of required images. We validate Jipole against ipole in two analytical tests and then compare pixelwise intensity derivatives from AD with those from finite-difference methods. We then demonstrate the utility of these gradients by performing parameter recovery for an analytical model in three increasingly complex cases for the injected image: ideal, blurred, and blurred with added noise. In most cases, high-accuracy fits are obtained in only a few optimization steps, failing only in cases with extremely low signal-to-noise ratios. These results highlight the potential of AD-based methods to accelerate robust, high-fidelity model-data comparisons in current and future black hole imaging efforts.

79 ASTRONOMY AND ASTROPHYSICS↗