High Energy Astrophysical Phenomena
High-energy radiation from the pulsar Equatorial Current Sheet
Ioannis Contopoulos, Jerome Petri and Ioannis Dimitropoulos
Pulsars emit beams of radiation that reveal the extreme physics of
neutron star magnetospheres. Yet, their understanding remains incomplete. Recent global Particle-in-
Cell (PIC) simulations have raised several questions that led us to question their validity and their extrapolation to realistic particle Lorentz factors, electric and
magnetic fields. We want to generate realistic sky maps of high-energy radiation from first principles. We propose a novel method to study the Equatorial Current Sheet (ECS) where most of the particle acceleration and the high-energy radiation is expected to originate. We first determine its shape and external
magnetic field with a steady-state ideal force-free solution. Then, we consider the extra electric and
magnetic field components that develop when dissipation is considered. Finally, we study the particle acceleration and radiation that is due to these extra field components for realistic field and particle parameters. We generate realistic sky maps of high-energy radiation and compare them with those obtained via PIC simulations. These sky maps may also be closely reproduced using the ECS of the split-monopole solution beyond the light cylinder. The ECS is probably stabilized by the normal
magnetic field component that is due to the global magnetospheric reconnection. Our method helps us better understand the origin of the pulsed high-energy radiation in the
pulsar magnetosphere.
Read the full abstract ▾
Earth and Planetary Astrophysics · Cross-listed from physics.hist-ph
The significance of Halley's Comet in the Bayeux Tapestry
Michael Lewis and Simon Portegies Zwart
A comet appears in the Bayeux Tapestry between the scene showing the death of the English king Edward the Confessor and the election of his successor, Harold Godwinson. The Tapestry's inscription only refers to this as a
star, though we can see from its depiction, shown with a hairy tail, that it is a comet, now known to us as Halleys Comet P1/Halley. Behind the exciting story of the Bayeux Tapestry, however, goes a rich mythological world of intrigue, deceit, the succession of kings, and earlier sightings of the same comet. In historical accounts of many cultures, comets are generally considered portents of change rather than disaster. Here we consider the significance of the Tapestry's comet in the context of the so-called English succession crisis of 1066 with reference to other contemporary accounts of comets. We conclude that although the tapestry's illustration is suggestive and unmistakably Halleys Comet, it is not a priority for its creators to give a precise account of its arrival in the sky, but rather connect it, likely for political reasons, albeit retrospectively, to the sequential events of Edwards death and Harolds coronation. In that, the tapestry's anonymous artists provide a unique telling of its arrival.
Read the full abstract ▾
Instrumentation and Methods for Astrophysics
Ice Giants Pre-Decadal Survey Mission Study Report
M. Hofstadter, A. Simon, K. Reh et al. (+26 more)
The Ice
Giants Study was commissioned by NASA to take a fresh look (as of 2017) at science priorities and concepts for missions to the Uranus and Neptune systems in preparation for the third Planetary Science Decadal Survey. This study was led by a Science Definition Team (SDT) and the Jet Propulsion Laboratory (JPL) with participation from Langley Research Center, Ames Research Center, The Aerospace Corporation, Purdue University, and the European Space Agency. The SDT was appointed by NASA and co-chaired by Mark Hofstadter (JPL) and Amy Simon (Goddard Space Flight Center). The Study Manager was Kim Reh (JPL) and the Study Lead was John Elliott (JPL). The study team assessed and prioritized science objectives taking into account advances since the last Decadal Survey, current and emerging technologies, mission implementation techniques and celestial mechanics. This study examined a wide range of mission architectures, flight elements, and instruments. Six of the prioritized concepts were studied via JPL's Team X process and resulting cost estimates were subjected to independent assessment by The Aerospace Corporation. Results presented herein show that high-value flagship-class mission concepts to either Uranus or Neptune are achievable within ground rule budgetary constraints.
Read the full abstract ▾
All authors M. Hofstadter, A. Simon, K. Reh, J. Elliott, S. Atreya, D. Banfield, J. Fortney, A. Hayes, M. Hedman, G. Hospodarsky, A. Masters, K. Mandt, M. Showalter, K. Soderlund, D. Turini, E. Turtle, P. Agrawal, T. Anderson, D. Atkinson, N. Arora, C. Borden, M. Brennan, J. Cutts, H. Hwang, M. Le, A. Petropoulos, S. Saikia, T. Spilker, W. Smythe
Cosmology and Nongalactic Astrophysics · Cross-listed from hep-th
Infrared Quantum Electrodynamics and the Rayleigh-Jeans Physics
Jorge Gamboa and Natalia Tapia Arellano
Infrared
quantum electrodynamics (IR-QED) acquires a natural
geometric interpretation once soft
photons are described as adiabatically transported electron-
photon clouds. Within this framework, the relevant infrared structure is encoded in a functional Berry phase associated with the space of gauge connections, and the corresponding Berry corrections modify the Rayleigh-Jeans spectrum. The infrared scaling
symmetry of the Rayleigh-Jeans law leads to a simple renormalization-group equation whose solution determines the frequency dependence of an effective factor $F_{\rm eff}(\omega)$ controlling the strength of the electron-
photon cloud dressing. As a result, the energy density of the
cosmic microwave background (CMB) receives a Berry-induced correction that scales as a power law and produces a frequency-dependent temperature excess in the
radio domain. Although the exponent $\gamma$ governing this scaling behaviour is not fixed internally by the present formulation of IR-QED and must instead be determined phenomenologically, the existence and structure of the excess are genuine predictions of the theory. Remarkably, the resulting expression is extremely simple and naturally aligns with the deviations suggested by the ARCADE 2 data. Taken together, these results indicate that Berry phases in IR-QED may lead to observable consequences in the low-frequency tail of the CMB spectrum.
Read the full abstract ▾
Cosmology and Nongalactic Astrophysics
Accessing the homogeneity scale with 21 cm intensity mapping surveys
Bruno B. Bizarria, Camila P. Novaes, Felipe Avila et al. (+4 more)
The homogeneity scale, $R_{\rm H}$, offers a fundamental test of the Cosmological Principle, yet it has not yet been measured with 21cm intensity mapping surveys. A key limitation for such a measurement is the
telescope beam, which artificially smooths the observed signal. We quantify this effect using the two-point correlation function and the correlation dimension, $\mathcal{D}_2(r)$, to model how beam convolution suppresses intrinsic clustering. For any given
redshift $z$, we identify a maximum beam width, $\sigma_{\rm max}(z)$, beyond which the homogeneity scale cannot be recovered. This limit defines an inaccessible region in the $\sigma \times z$ parameter space, where $R_{\rm H}$ is erased by beam smoothing. Applying this framework to several current and upcoming
radio telescopes, we assess their ability to probe $R_{\rm H}$. Our results provide the first quantitative forecast of the instrumental requirements for measuring the
cosmic homogeneity scale with 21cm IM, and establish a theoretical basis for future observational applications.
Read the full abstract ▾
All authors Bruno B. Bizarria, Camila P. Novaes, Felipe Avila, Rahima Mokeddem, Helissa H. da Costa, Carlos A. Wuensche, Gabriel A. S. Silva
Cosmology and Nongalactic Astrophysics
Thermodynamic Origin of the Tully-Fisher Relation in Dark Matter Dominated Galaxies: A Theoretical-Empirical Derivation
V.K. Oikonomou
In this work we introduce the concept of self-interacting dark matter with scale-dependent equation of state, in the context of which dark matter is collisional and its equation of state is radius-dependent and has the form $P(r)=K(r)\left(\frac{\rho(r)}{\rho_{\star}}\right)^{\gamma(r)}$. We confronted the effectively 2-parameter model with 174
galaxies from the SPARC data, and we found that the rotation curves of 100
galaxies can be perfectly fitted by the model. These
galaxies are dark matter dominated, mostly
dwarfs, low-luminosity and low-surface-brightness spiral
galaxies. We demonstrate that scale-dependent self-interacting dark matter solves the cusp-core issue for dark matter dominated
galaxies. More importantly, the structure of the scale-dependent SIDM model produces in a semi-theoretically and semi-empirically way the canonical Tully-Fisher and the baryonic Tully-Fisher relations when these 100 viable
dwarfs, low-surface-brightness and low-luminosity
galaxies are taken into account. The behavior of the
entropy function $K(r)$ is assumed to be $K(r)=K_0\times\left(1+\frac{r}{r_c}\right)^{-p}$. The perfect fits of the rotation curves come for a nearly isothermal and virialized dark matter
halo, which naturally predicts the correlation $K_0\sim V_{\mathrm{max}}^2$. This correlation holds true empirically as confirmed by the data and we also find empirically $L\sim K_0^2$ from the data, thus the canonical Tully-Fisher relation is reproduced semi-theoretically and semi-empirically. We perform the same task and we find theoretically, for dark matter dominated
galaxies, that $K_0\sim V_{\mathrm{flat}}^2$ which is also confirmed empirically from the data, along with the correlation $K_0\sim M_b^{0.5}$, hence the baryonic Tully-Fisher law naturally emerges in a semi-theoretical and semi-empirical manner.
Read the full abstract ▾
High Energy Astrophysical Phenomena
Multi-layered Fast Wind observed in XMM-Newton snapshot of Seyfert 1 Markarian 877
Xin Xiang, Jon M. Miller, Ehud Behar et al. (+4 more)
Ultra Fast
Outflows (UFOs) are powerful, highly ionized winds launched from the innermost regions of Active
Galactic Nuclei (AGNs), reaching velocities of 0.03 – 0.3 c and playing a key role in AGN feedback. We present a photoionization analysis of an 18 ks \xmm\ snapshot of the Seyfert 1 AGN Mrk 877, revealing three distinct UFO components with line-of-sight velocities of $0.10^{+0.005}_{-0.005}~c$ , $ 0.04^{+0.005}_{-0.004}~c$ , and $0.05^{+0.005}_{-0.004}~c$. These components span a broad range of ionization parameters and column densities, producing absorption features across both soft and hard X-ray bands. Even under the most conservative assumption for the volume filling factor, the fastest component exceeds $5\%$ of the Eddington luminosity, making it capable of driving strong
galaxy-scale feedback. The soft X-ray UFO component, despite its lower ionization, shares a similar velocity as a higher-ionization component, hinting at a two-phase medium likely shaped by clumpiness or interactions with ambient material. The density profile inferred from the Absorption Measurement Distributions (AMD) and the positive trend between
outflow momentum rate and radiation momentum flux suggest that wind is powered by a combination of radiative and
magnetic driving.
Read the full abstract ▾
All authors Xin Xiang, Jon M. Miller, Ehud Behar, W. N. Brandt, Luigi Gallo, Doyee Byun, Elena Gallo
Astrophysics of Galaxies
Mapping the Cosmic-Ray Ionization Rate in the Local Galaxy with H$_3^+$
Nick Indriolo, Alexei V. Ivlev, T. Pellegrin et al. (+6 more)
Chemistry in diffuse molecular clouds relies primarily on rapid ion-molecule reactions. Formation of the initial ions, H$^+$ and H$_2^+$, is dominated by
cosmic-ray ionization of H and H$_2$, making the
cosmic-ray ionization rate (denoted $\zeta({\rm X})$ for species X) an important parameter for chemical modeling. We have made observations targeting absorption lines of H$_3^+$, one of the most reliable tracers of $\zeta({\rm H_2})$, toward diffuse molecular cloud sight lines where the H$_2$ column density has been directly measured in the ultraviolet, detecting H$_3^+$ in 12 out of 27 sight lines. The 3D-PDR modeling method introduced by Obolentseva et al. (2024) was used to infer
cosmic-ray ionization rates in the clouds along these sight lines, and our combined sample has a mean ionization rate of $5.3\times10^{-17}$ s$^{-1}$ with standard deviation $2.5\times10^{-17}$ s$^{-1}$. By associating H$_3^+$ absorption with gas density peaks derived from the differential extinction maps of Edenhofer et al. (2024) we have constructed a sparsely sampled 3D map of the
cosmic-ray ionization rate in targeted regions within about 1~kpc of the Sun. Specific regions show reasonably uniform ionization rates over length scales of tens of parsecs, with the average ionization rate in each region being different. Large differences (factor of 5) in $\zeta({\rm H_2})$ are found over length scales of about 100 pc. This supports a picture where the
cosmic-ray ionization rate varies smoothly over small size scales, but is not uniform everywhere in the
Galactic disk, likely being controlled by proximity to particle acceleration sites.
Read the full abstract ▾
All authors Nick Indriolo, Alexei V. Ivlev, T. Pellegrin, M. Obolentseva, Paola Caselli, A. M. Jacob, D. A. Neufeld, Kedron Silsbee, M. G. Wolfire
High Energy Astrophysical Phenomena
Blazars as a Potential Origin of the KM3-230213A Event
The KM3NeT Collaboration
The KM3NeT collaboration has reported the detection of the highest energy
neutrino event observed to date. The energy of the event is of the order of 220 PeV hinting towards a
neutrino flux at the highest energies. In this article, the potential blazar origin for this event is explored. The publicly available Astro-Multimessenger Modeling software is used to model the blazar gamma-ray and
neutrino fluxes. It is concluded that a population of blazars could produce the diffuse flux compatible with the observation of the ultra-high energy event detected by the KM3NeT/ARCA detector. At the same time, the gamma-ray flux produced by such a population of blazars is consistent with the diffuse gamma-ray flux measured by the
Fermi Large Area
Telescope.
Read the full abstract ▾
Cosmology and Nongalactic Astrophysics
Large Scale White Noise and Cosmology
Gabriela Barenboim, Aurora Ireland and Albert Stebbins
The generation of white noise on large scales is a generic property of the dynamics of physical systems described by local non-linear partial differential equations. Non-linearities prevent the small scale dynamics from being erased by smoothing. Unresolved small scale dynamics act as an uncorrelated (white or Poissonian) noise (seemingly
stochastic but actually deterministic) contribution to large scale dynamics. This white noise exists even when the dynamics is very nearly linear. In cases where the power spectrum is sub-Poissonian on large scales, this noise will dominate on the largest scale power no matter the amplitude of the inhomogeneities. Such is the case in the standard model of cosmology, where the
primordial density power spectrum is expected to have an almost Harrison-Zel'dovich, $P[k]\sim k$, spectrum on a much broader range of scales than can be observed. Even though linear gravitational
evolution dominates non-linear corrections by a factor $\sim10^5$, the non-observation of white noise on the Hubble scale precludes the extrapolation of this power law below the comoving $1\,$pc scale. More generally, observation or non-observation of large scale white noise provides a powerful probe of the universe on very small scales in the early early universe. Gravitational radiation, phase transitions, vorticity, and running of the spectral index are all phenomena that can be probed with large scale white noise. Large scale white noise is a non-optional feature of all cosmological models but one which has not heretofore been appreciated.
Read the full abstract ▾
Solar and Stellar Astrophysics
Dual Origins of Rapid Flare Ribbon Downflows in an X9-class Solar Flare
Ryan J. French, William H. Ashfield IV, Cole A. Tamburri et al. (+4 more)
We detect rapid downflows of 150-217 km/s in IRIS Si IV 1402.77 nm measurements of an X9-class
solar flare on 2024 October 3rd. The fast
redshift values persist for over 15 minutes from
flare onset, and can be split into two distinct stages of behavior, suggesting that multiple mechanisms are responsible for the downwards acceleration of
flare ribbon
plasma. The first stage of rapid downflows are synchronized with peaks in emission from the Advanced Space-based
Solar Observatory Hard X-ray Imager (ASO-S/HXI) and Large Yield Radiometer (LYRA) Lyman-alpha measurements, indicative that the chromospheric downflows (with a maximum
redshift of 176 km/s) result from chromospheric condensations associated with impulsive energy release in the
solar flare. Later in the event, strong Si IV
flare ribbon downflows persist (to a maximum value of 217 km/s), despite the
magnetic flux rate falling to zero, and high-energy HXR and Lyman-alpha measurements returning to background levels. This is reflective of downflows in the
flare ribbon footpoints of
flare-induced coronal rain. Hard X-ray spectral analysis supports this scenario, revealing strong non-thermal emission during the initial downflow stage, falling near background levels by the second stage. Despite these distinct and contrasting stages of ribbon behavior, Si IV Doppler velocities exhibit quasi-periodic pulsations with a constant ~50 s period across the 15-minute
flare evolution (independently of loop length). We deduce that these pulsations are likely caused by MHD oscillations in the
magnetic arcade. Finally, we utilize machine learning K-means clustering methods to quantify line profile variations during the stages of rapid downflows.
Read the full abstract ▾
All authors Ryan J. French, William H. Ashfield IV, Cole A. Tamburri, Maria D. Kazachenko, Marie Dominique, Marcel Corchado Albelo, Amir Caspi
High Energy Astrophysical Phenomena
The Scoured Spike: Suppression of Indirect Dark Matter Signals by a Hidden Companion
Jaden Lopez and Stefano Profumo
A massive “dark companion” – such as an intermediate-mass black hole or other compact dark object – orbiting the supermassive black hole at the
Galactic Center can dynamically reshape the surrounding dark-matter spike. Through gravitational heating and angular-momentum exchange, the companion excavates a “scoured” region that lowers the inner density and suppresses the expected annihilation signal. We quantify this effect by computing the suppression of the dark-matter annihilation $J$-factor induced by such a companion, combining an analytic scouring-radius model with full numerical integrations of the modified density profile. We scan the parameter space of companion mass, orbital separation, system age, and spike slope, explicitly including the interplay with the annihilation plateau. \textcolor{black}{For canonical Gondolo–Silk spikes with $\gamma_{\rm sp} \gtrsim 2$, the scouring radius grows only weakly with injected energy because the binding energy is concentrated at small radii; however, once scouring occurs, the $J$-factor suppression can be substantial, for sufficiently massive or long-lived companions. In the shallower-spike regime with $\gamma_{\rm sp} \lesssim 2$, the same companion drives $r_{\rm scour}$ to much larger multiples of $R_{\rm core}$, and even a modest ${\sim}10^{4}\,M_\odot$ companion on a $\mathcal{O}(100)\,{\rm AU}$ orbit and $\mathcal{O}({\rm Gyr})$ age can suppress the annihilation flux by one to two orders of magnitude. Whether a steep spike is significantly suppressed therefore depends on the companion's mass and lifetime, not merely on the sign of $(2-\gamma_{\rm sp})$.} The numerical results are accurately captured (typically at the $\lesssim 10\%$ level) by a simple fitting formula in terms of a dimensionless scouring parameter that measures the ratio between the scoured region and the annihilation core. [...]
Read the full abstract ▾
High Energy Astrophysical Phenomena
Low-energy Radio Bursts from Magnetar XTE J1810$-$197: Implications for Fast Radio Bursts
Banshi Lal, Yogesh Maan, Moaz Abdelmaguid et al. (+3 more)
Magnetars are the leading candidate sources of fast
radio bursts (FRBs). However, the observational probes of the connections between magnetars and FRBs are severely limited by the paucity of detection of highly energetic
radio events from magnetars – to date, only one
radio burst as energetic as FRBs has been detected from a
Galactic magnetar. Here, we present a detailed analysis of a large sample of low-energy bursts detected from the magnetar XTE J1810$-$197, and probe their implications for FRB emission from magnetars. We report detection of over 97000 bright
radio pulses from 242 observations of the magnetar XTE J1810$-$197 over 4.5 years and two decades in frequency (300 MHz to 6.15 GHz), using the
Giant Meterwave
Radio Telescope and the Green Bank
Telescope, after its recent outburst onset in December 2018. We present detailed analysis of the burst fluence distributions and their trends with time as well as frequency, and the waiting time distribution. We show that XTE J1810$-$197 rapidly switches between
pulsar-like and
giant-pulse-like emission states, and magnetars like XTE J1810$-$197 remain viable and likely emitters of FRBs, in the form of
giant-pulses with energies comparable to FRBs. We also demonstrate that the lack of the detection of an underlying periodicity in the bursts from repeating FRBs might be caused by emission across a wide range of
spin phases.
Read the full abstract ▾
All authors Banshi Lal, Yogesh Maan, Moaz Abdelmaguid, Visweshwar Ram Marthi, Joseph D. Gelfand, Samayra Straal
Cosmology and Nongalactic Astrophysics
Exploring the Model Dependence of MCMC-Based 21 cm Power Spectrum Parameter Constraints
August Berklas and Jonathan Pober
Detection and analysis of the
cosmic 21 cm signal of neutral hydrogen has long been considered the most promising route towards exploration of the Epoch of Reionization (EoR). 21CMMC, a Markov Chain Monte Carlo sampler of the seminumerical simulation code 21cmFAST, has already been used in conjunction with published upper limits on the 21 cm signal from the Murchison Widefield Array, the Low Frequency Array, and the Hydrogen Epoch of Reionization Array to constrain the astrophysics of the EoR. Here, we investigate the extent to which analysis of the EoR performed using 21CMMC is dependent on the underlying seminumerical model. We used 21cmFAST to simulate two
datasets of 21 cm light-cones that differ only in the
algorithm used to identify ionized regions (the so-called "bubble-finding"
algorithm). We then tested 21CMMC's ability to return key astrophysical parameteres when using the different bubble-finding
algorithms. We find that the performance of 21CMMC depends sensitively on the agreement between the astrophysical model of our mock data and the model used for sampling. This result has important implications for the analysis of the 21 cm signal performed using 21CMMC and further motivates investigation into model-independent analysis techniques for 21 cm EoR data.
Read the full abstract ▾
Astrophysics of Galaxies
How accurately can obscured galaxy luminosities be measured using spectral energy distribution fitting of near- through far-infrared observations?
Duncan Farrah, Kiana Ejercito, Andreas Efstathiou et al. (+23 more)
Infrared-luminous
galaxies are important sites of
stellar and black hole mass assembly at most
redshifts. Their luminosities are often estimated by fitting spectral energy distribution (SED) models to near- to far-infrared data, but the dependence of these estimates on the data used is not well-understood. Here, using observations simulated from a well-studied local sample, we compare the effects of wavelength coverage, signal-to-noise (S/N), flux calibration, angular resolution, and
redshift on the recovery of starburst, AGN, and host luminosities. We show that the most important factors are wavelength coverage that spans the peak in a SED, with dense wavelength sampling. Such observations recover starburst and AGN infrared luminosities with systematic bias below $20\%$. Starburst luminosities are best recovered with far-infrared observations while AGN luminosities are best recovered with near- and mid-infrared observations, though the recovery of both are enhanced with near/mid-infrared, and far-infrared observations, respectively. Host luminosities are best recovered with near/far-infrared observations, but are usually biased low, by $\gtrsim20\%$. The recovery of starburst and AGN luminosity is enhanced by observing at high angular resolution. Starburst-dominated systems show more biased recovery of luminosities than do AGN-dominated systems. As
redshift increases, far-infrared observations become more capable, and mid-infrared observations less capable, at recovering luminosities. Our results highlight the transformative power of a far-infrared instrument with dense wavelength coverage from tens to hundreds of microns for studying infrared-luminous
galaxies. We tabulate estimates of systematic bias and random error for use with JWST and other observatories.
Read the full abstract ▾
All authors Duncan Farrah, Kiana Ejercito, Andreas Efstathiou, David Leisawitz, Athena Engholm, Irene Shivaei, Matteo Bonato, David L. Clements, Sara Petty, Lura. K. Pitchford, Charalambia Varnava, Jose Afonso, Carlotta Gruppioni, Evanthia Hatziminaoglou, Andrew Hoffman, Mark Lacy, Brenda C. Matthews, Conor Nixon, Chris Pearson, Berke Vow Ricketti, Dimitra Rigopoulou, Loren Robinson, Locke D. Spencer, Lingyu Wang, David B. Sanders, Gerard van Belle
Earth and Planetary Astrophysics
Inferring Planet and Disk Parameters from Protoplanetary Disk Images Using a Variational Autoencoder
Sayed Shafaat Mahmud, Sayantan Auddy, Neal Turner et al. (+1 more)
Dust-continuum observations of many protoplanetary disks reveal rings and gaps that are widely interpreted as evidence of ongoing
planet formation. Here we present the first framework for inferring
planet and disk parameters from such images using variational autoencoder (VAE) based generative machine learning (ML). The new framework is called VADER (Variational Autoencoder for Disks with Embedded Rings). We train VADER on synthetic images of
dust continuum emission, generated from FARGO3D hydrodynamic simulations post-processed with Monte Carlo radiative transfer calculations. VADER infers the masses of up to three embedded
planets as well as the disk parameters viscous $\alpha$,
dust-to-gas ratio, Stokes number, and flaring index. VADER returns a full posterior distribution for each of these quantities. We demonstrate that VADER reconstructs disk morphologies with high structural similarity (index $>$ 0.99), accurately recovers
planet parameters with $R^2 > 0.9$ across
planet masses, and reliably predicts disk parameters. Applied to ALMA
dust continuum images of 23 protoplanetary disks, our model returns mass estimates for embedded
planets of 0.3-2~$M_{\mathrm{Jup}}$ that agree to within $1\sigma$ of published values in most cases, and infers disk parameters consistent with current literature. Once trained, the VAE performs full posterior parameter inference in a matter of minutes, offering statistical rigor with enough computational speed for application to large-scale ALMA surveys. These results establish VAE-based models as powerful tools for inferring from disk structure the masses of embedded
planets and the global disk parameters, with their associated uncertainties.
Read the full abstract ▾
All authors Sayed Shafaat Mahmud, Sayantan Auddy, Neal Turner, Jeffrey S. Bary
High Energy Astrophysical Phenomena · Cross-listed from gr-qc
Characteristic Decomposition for Relativistic Numerical Simulations: II. Magnetohydrodynamics
Saul A. Teukolsky
The characteristic decomposition for GRMHD in the comoving frame of the fluid has been known for a long time. However, it has not been known in the coordinate frame of the simulation and in terms of the conserved variables evolved in typical numerical simulations. This paper applies the method of quasi-invertible transformations developed in Paper I to derive this decomposition. Among other benefits, this will allow us to use the most accurate known computational methods, such as full-wave Riemann solvers. The results turn out to be simpler than expected based on earlier attempts.
Read the full abstract ▾
High Energy Astrophysical Phenomena · Cross-listed from gr-qc
Characteristic Decomposition for Relativistic Numerical Simulations: I. Hydrodynamics
Saul A. Teukolsky
The characteristic decomposition for GRMHD is not known in a form useful for current numerical simulations. This prevents us from using the most accurate known computational methods, such as full-wave Riemann solvers. In this paper, we present a new method of finding decompositions. The method is based on transformations from the comoving frame, where the fluid flow is simplest and the decomposition has been known for a long time. The key innovation we introduce is that of quasi-invertible transformations. In this first paper, we introduce these transformations using the simpler example of relativistic hydrodynamics. We recover the known decomposition for relativistic hydrodynamics in somewhat simpler form than previously derived, and without the need for computer
algebra. A new result in this paper is the characteristic decomposition when the the
evolution tracks the composition of a fluid in
nuclear statistical
equilibrium. In Paper II of this series, we apply a quasi-invertible transformation to derive the complete characteristic decomposition for GRMHD in the conserved variables used in simulations.
Read the full abstract ▾
High Energy Astrophysical Phenomena
Estimation of Magnetic Fields from Synchrotron Emission: Numerical Tests
Nora B. Linzer, Lucia Armillotta, Eve C. Ostriker et al. (+1 more)
We use models of spectrally resolved
cosmic ray (CR) transport in TIGRESS MHD simulations of the local ISM to produce synthetic synchrotron emission and to test, on scales from a few kpc down to ~10 pc, the traditional estimate of
magnetic field strength based on the assumption of equipartition between the
magnetic and total CR energy densities. Our analysis shows that the traditional equipartition estimate works well at the kpc scale of the simulation box, but breaks down at smaller scales. We find that the predicted
magnetic field strength can be improved at small scales by assuming a constant CR energy density across each mock
radio observation. The large-scale mean CR energy density can be estimated by assuming equipartition with the large-scale mean
magnetic energy density, or as a function of additional observable quantities such as the
star formation rate surface density or gas weight. In addition to estimating the
magnetic field strength, we use synthetic polarized emission to create maps of the
magnetic field direction. We find that the true
magnetic field direction can be recovered well from the mock observations.
Read the full abstract ▾
All authors Nora B. Linzer, Lucia Armillotta, Eve C. Ostriker, Eliot Quataert