SEARCH · Engineering Papers

Results for “variable stars”

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.

SourceResearch connectionFrom yearThrough year

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

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Wavelength Dependence of Activity-induced Photometric Variations for Young Cool Stars in Hyades

We investigate photometric variations due to stellar activity that induce systematic radial-velocity errors (so-called “jitter”) for the four targets in the Hyades open cluster observed by the K2 mission (EPIC 210721261, EPIC 210923016, EPIC 247122957, and EPIC 247783757). Applying Gaussian process regressions to the K2 light curves and the near-infrared (NIR) light curves observed with the IRSF 1.4 m telescope, we derive the wavelength dependences of the photometric signals due to stellar activity. To estimate the temporal variations in the photometric variability amplitudes between the two observation periods of K2 and IRSF, separated by more than 2 yr, we analyze a number of K2 targets in Hyades that have also been observed in Campaigns 4 and 13 and find a representative variation rate over 2 yr of 38% ± 71%. Taking this temporal variation into account, we constrain projected sizes and temperature contrast properties of the starspots in the stellar photosphere to be approximately 10% and 0.95%, respectively. These starspot properties can induce relatively large differences in the variability amplitude over different observational passbands, and we find that radial-velocity jitter may be more suppressed in the NIR than previously expected. Our result supports profits of ongoing exoplanet search projects that are attempting to detect or confirm young planets in open clusters via radial-velocity measurements in the NIR.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

Spatially resolved properties of early-type group-dominant galaxies with MUSE: gas content, ionization mechanisms, and metallicity gradients

ABSTRACT With the goal of a thorough investigation of the ionized gas and its origin in early-type group-dominant galaxies, we present archival MUSE data for 18 galaxies from the Complete Local-Volume Groups Sample (CLoGS). These data allowed us to study the spatially resolved warm gas properties, including the morphology of the ionized gas, EW(H α), and kinematics as well as the gas-phase metallicity (12 + log(O/H)) of these systems. In order to distinguish between different ionization mechanisms, we used the emission-line ratios [O iii]/H β and [N ii]/H α in the BPT diagrams and EW(H α). We find that the ionization sources in our sample have variable impacts at different radii; central regions are more influenced by low-luminosity active galactic nuclei, while extended regions of low-ionization nuclear emission-line region-like emission are ionized by other mechanisms with post-asymptotic giant branch stars photoionization likely contributing significantly. We classified our sample into three H α +[N ii] emission morphology types. We calculate the gas-phase metallicity assuming several methods and ionization sources. In general, 12 + log(O/H) decreases with radius from the centre for all galaxies, independently of nebular morphology type, indicating a metallicity gradient in the abundance profiles. Interestingly, the more extended filamentary structures and all extranuclear star-forming regions present shallow metallicity gradients. Within the uncertainties these extended structures can be considered chemically homogeneous. We suggest that group-dominant galaxies in our sample likely acquired their cold gas in the past as a consequence of one or more mechanisms, e.g. gas-clouds or satellite mergers/accretion and/or cooling flows that contribute to the growth of the ionized gas structures.

Lagos, P. (ORCID:0000000223218657)↗

DOE OSTI2022

Quasar Identification Using Multivariate Probability Density Estimated from Nonparametric Conditional Probabilities

Nonparametric estimation for a probability density function that describes multivariate data has typically been addressed by kernel density estimation (KDE). A novel density estimator recently developed by Farmer and Jacobs offers an alternative high-throughput automated approach to univariate nonparametric density estimation based on maximum entropy and order statistics, improving accuracy over univariate KDE. This article presents an extension of the single variable case to multiple variables. The univariate estimator is used to recursively calculate a product array of one-dimensional conditional probabilities. In combination with interpolation methods, a complete joint probability density estimate is generated for multiple variables. Good accuracy and speed performance in synthetic data are demonstrated by a numerical study using known distributions over a range of sample sizes from 100 to 10 6 for two to six variables. Performance in terms of speed and accuracy is compared to KDE. The multivariate density estimate developed here tends to perform better as the number of samples and/or variables increases. As an example application, measurements are analyzed over five filters of photometric data from the Sloan Digital Sky Survey Data Release 17. The multivariate estimation is used to form the basis for a binary classifier that distinguishes quasars from galaxies and stars with up to 94% accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Image Deconvolution and Point-spread Function Reconstruction with STARRED: A Wavelet-based Two-channel Method Optimized for Light-curve Extraction

We present starred, a point-spread function (PSF) reconstruction, two-channel deconvolution, and light-curve extraction method designed for high-precision photometric measurements in imaging time series. An improved resolution of the data is targeted rather than an infinite one, thereby minimizing deconvolution artifacts. In addition, starred performs a joint deconvolution of all available data, accounting for epoch-to-epoch variations of the PSF and decomposing the resulting deconvolved image into a point source and an extended source channel. The output is a high-signal-to-noise-ratio, high-resolution frame combining all data and the photometry of all point sources in the field of view as a function of time. Of note, starred also provides exquisite PSF models for each data frame. We showcase three applications of starred in the context of the imminent LSST survey and of JWST imaging: (i) the extraction of supernovae light curves and the scene representation of their host galaxy; (ii) the extraction of lensed quasar light curves for time-delay cosmography; and (iii) the measurement of the spectral energy distribution of globular clusters in the "Sparkler," a galaxy at redshift z = 1.378 strongly lensed by the galaxy cluster SMACS J0723.3-7327. starred is implemented in jax, leveraging automatic differentiation and graphics processing unit acceleration. This enables the rapid processing of large time-domain data sets, positioning the method as a powerful tool for extracting light curves from the multitude of lensed or unlensed variable and transient objects in the Rubin-LSST data, even when blended with intervening objects.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

The Three Hundred Project: Modeling baryon and hot-gas fraction evolution in simulated clusters

The baryon fraction of galaxy clusters, expressed as the ratio between the mass in baryons (including both stars and cold or hot gas) and the total mass, is a powerful tool to provide information on the cosmological parameters, while the hot-gas fraction provides indications on the physics of the intracluster plasma and its interplay with the processes that drive galaxy formation. Using cosmological hydrodynamical simulations of about 300 simulated massive galaxy clusters with a median mass M 500 ≈ 7 × 10 14 M ⊙ at z = 0, we model the relations between total mass and either baryon fraction or the hot gas fractions at overdensities Δ = 2500, 500, and 200 with respect to the cosmic critical density, and their evolution from z ∼ 0 to z ∼ 1.3. We utilized the simulated galaxy clusters from the Three Hundred project, which include star formation and feedback from both supernovae and active galactic nuclei. We fit the simulation results for such scaling relations against three analytic forms (linear, quadratic, and logarithmic in a logarithmic plane) and three forms for the redshift dependence, and we considered as a variable both the inverse of the cosmic scale factor, (1 + z), and the Hubble expansion rate, E(z). We show that power-law dependencies on cluster mass poorly describe the investigated relations. A power law fails to simultaneously capture the flattening of the total baryon and gas fractions at high masses, their drop at low masses, and the transition between these two regimes. The other two functional forms provide a more accurate description of the curvature in mass scaling. The fractions measured within smaller radii exhibit a stronger evolution than those measured within larger radii. From the analysis of these simulations, we evince that as long as we include systems in the mass range herein investigated, the baryon or gas fraction can be accurately related to the total mass through either a parabola or a logarithm in the logarithmic plane. The trends are common to all modern hydro simulations, although the amplitude of the drop at low masses might differ. Being able to observationally determine the gas fraction in groups will thus provide constraints on the baryonic physics.

galaxy clusters↗

DOE OSTI2023

Empirical radius formulas for canonical neutron stars from bidirectionally selecting features of equations of state in extended Bayesian analyses of observational data

Significant advancement in Bayesian inference of nuclear equation of state (EOS) from gravitational wave and x-ray observations of neutron stars (NSs) has been made by the nuclear astrophysics community especially since GW170817. By extending the traditional Bayesian analysis which normally ends at presenting the marginalized posterior probability distribution functions (PDFs) of individual EOS parameters and their correlations (or sometimes only the Pearson correlation coefficients which are only reliably useful when the variables are linearly correlated while they are actually often not), we search for a data-driven and robust empirical formula for the radius 𝑅 1.4 of canonical NSs in terms of the characteristic EOS parameters (features). We also identify the single most important but currently poorly known EOS parameter for determining the 𝑅 1.4 . Using three regression-model-building methodologies: bidirectional stepwise feature selection, least absolute shrinkage selection operator (LASSO) regression, and neural network regression on a large set of posterior EOSs and the corresponding 𝑅 1.4 values inferred from earlier comprehensive Bayesian analyses of NS observational data, we systematically and rigorously develop the most probable 𝑅 1.4 formulas with varying statistical accuracy and technical complexity. Here, the most important EOS parameters for determining 𝑅 1.4 are found consistently in each of the feature selection processes to be (in order of decreasing importance): curvature 𝐾 sym , slope 𝐿, skewness 𝐽 sym of nuclear symmetry energy, skewness 𝐽 0 , incompressibility 𝐾 0 of symmetric nuclear matter, and the magnitude 𝐸 sym ⁡(𝜌 0 ) of symmetry energy at the saturation density 𝜌 0 of nuclear matter.

Bayesian methods↗

DOE OSTI2023

Neutron star mass estimates from gamma-ray eclipses in spider millisecond pulsar binaries

Reliable neutron star mass measurements are key to determining the equation of state of cold nuclear matter, but such measurements are rare. Black widows and redbacks are compact binaries consisting of millisecond pulsars and semi-degenerate companion stars. Spectroscopy of the optically bright companions can determine their radial velocities, providing inclination-dependent pulsar mass estimates. Although inclinations can be inferred from subtle features in optical light curves, such estimates may be systematically biased due to incomplete heating models and poorly understood variability. Using data from the Fermi Large Area Telescope, we have searched for gamma-ray eclipses from 49 spider systems, discovering significant eclipses in 7 systems, including the prototypical black widow PSR B1957+20. Gamma-ray eclipses require direct occultation of the pulsar by the companion, and so the detection, or significant exclusion, of a gamma-ray eclipse strictly limits the binary inclination angle, providing new robust, model-independent pulsar mass constraints. For PSR B1957+20, the eclipse implies a much lighter pulsar (1.81 ± 0.07 solar masses) than inferred from optical light curve modelling.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

A Deep, High-angular-resolution 3D Dust Map of the Southern Galactic Plane

We present a deep, high-angular-resolution 3D dust map of the southern Galactic plane over 239° < l < 6° and ∣ b∣ < 10° built on photometry from the DECaPS2 survey, in combination with photometry from VISTA Variables in the Via Lactea, the Two Micron All Sky Survey, and “Unofficial” Wide-field Infrared Survey Explorer and parallaxes from Gaia Data Release 3 where available. To construct the map, we first infer the distance, extinction, and stellar types of over 700 million stars using the brutus stellar inference framework with a set of theoretical MESA Isochrone and Stellar Tracks (MIST) stellar models. Our resultant 3D dust map has an angular resolution of 1′ , roughly an order of magnitude finer than existing 3D dust maps and comparable to the angular resolution of the Herschel 2D dust emission maps. We detect complexes at the range of distances associated with the Sagittarius-Carina and Scutum-Centaurus arms in the fourth quadrant, as well as more distant structures out to a maximum reliable distance of d ≈ 10 kpc from the Sun. The map is sensitive up to a maximum extinction of roughly A V ≈ 12 mag. We publicly release both the stellar catalog and the 3D dust map, the latter of which can easily be queried via the Python package dustmaps. When combined with the existing Bayestar19 3D dust map of the northern sky, the DECaPS 3D dust map fills in the missing piece of the Galactic plane, enabling extinction corrections over the entire disk ∣b∣ < 10°. Our map serves as a pathfinder for the future of 3D dust mapping in the era of LSST and Roman, targeting regimes accessible with deep optical and near-infrared photometry but often inaccessible with Gaia.

Milky Way galaxy↗

DOE OSTI2026

Robust Measurement of Stellar Streams around the Milky Way: Correcting Spatially Variable Observational Selection Effects in Optical Imaging Surveys

Observations of density variations in stellar streams are a promising probe of low-mass dark matter substructure in the Milky Way. However, survey systematics such as variations in seeing and sky brightness can also induce artificial fluctuations in the observed densities of known stellar streams. These variations arise because survey conditions affect both object detection and star–galaxy misclassification rates. To mitigate these effects, we use Balrog synthetic source injections in the Dark Energy Survey (DES) Y3 data to calculate detection rate variations and classification rates as functions of survey properties. We show that these rates are nearly separable with respect to survey properties and can be estimated with sufficient statistics from the synthetic catalogs. Applying these corrections reduces the standard deviation of relative detection rates across the DES footprint by a factor of 5, and our corrections significantly change the inferred linear density of the Phoenix stream when including faint objects. Additionally, for artificial streams with DES-like survey properties we are able to recover density power spectra with reduced bias. We also find that uncorrected power-spectrum results for Legacy Survey of Space and Time (LSST)-like data can be around 5 times more biased, highlighting the need for such corrections in future ground-based surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Magnetic Navigation for GPS-Denied Airborne Applications

Most current flight systems are dependent on GPS for navigation. Recently, however, navigation in GPS-denied environments has become an area of intensive research. Additional navigation sensor data can be obtained from visual observations (stars or terrain), inertial measurement units, radar, measurements of the local magnetic field, or perhaps even gravity. Absolute and relative positioning via magnetic field measurements have been shown to be viable in many applications including ground navigation, low altitude aircraft flight, and spaceflight. There is greater variability in the magnetic field over shorter distances when flying at low altitude and in ground applications, leading to more accurate positioning. However, ground-based magnetic navigation is often heavily influenced by man-made structures, especially in urban environments. This is not the case for airborne magnetic navigation since the influence of buildings, roads, etc. is negligible for typical aircraft altitudes. For absolute magnetic navigation, the positioning accuracy decreases as altitude increases for a given vehicle velocity, but the observed time variability in the field can be reclaimed by traveling faster through the field. Thus, navigation accuracy becomes a balance of speed and altitude since the higher altitude can be counterbalanced by higher velocity. To understand these effects quantitatively, we explored various techniques to aid a simulated inertial measurement unit with magnetic information. Using a technique known as two-dimensional magnetic map matching, we simulated the performance of airborne magnetic navigation at fixed speed while varying the altitude, flight direction, magnetometer data collection time, reference magnetic map bias error, and type of trajectory (over land or over ocean).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Orbital Motion, Variability, and Masses in the T Tauri Triple System

We present results from adaptive optics imaging of the T Tauri triple system obtained at the Keck and Gemini Observatories in 2015−2019. We fit the orbital motion of T Tau Sb relative to Sa and model the astrometric motion of their center of mass relative to T Tau N. Using the distance measured by Gaia, we derived dynamical masses of M{sub Sa}=2.05±0.14 M {sub ⊙} and M {sub Sb} = 0.43 ± 0.06 M{sub ⊙}. The precision in the masses is expected to improve with continued observations that map the motion through a complete orbital period; this is particularly important as the system approaches periastron passage in 2023. Based on published properties and recent evolutionary tracks, we estimate a mass of ∼2 M {sub ⊙} for T Tau N, suggesting that T Tau N is similar in mass to T Tau Sa. Narrowband infrared photometry shows that T Tau N remained relatively constant between late 2017 and early 2019 with an average value of K = 5.54 ± 0.07 mag. Using T Tau N to calibrate relative flux measurements since 2015, we found that T Tau Sa varied dramatically between 7.0 and 8.8 mag in the K band over timescales of a few months, while T Tau Sb faded steadily from 8.5 to 11.1 mag in the K band. Over the 27 yr orbital period of the T Tau S binary, both components have shown 3–4 mag of variability in the K band, relative to T Tau N.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Red supergiant candidates for multimessenger monitoring of the next Galactic supernova

ABSTRACT We compile a catalogue of 578 highly probable and 62 likely red supergiants (RSGs) of the Milky Way, which represents the largest list of Galactic RSG candidates designed for continuous follow-up efforts to date. We match distances measured by Gaia DR3, 2MASS photometry, and a 3D Galactic dust map to obtain luminous bright late-type stars. Determining the stars’ bolometric luminosities and effective temperatures, we compare to Geneva stellar evolution tracks to determine likely RSG candidates, and quantify contamination using a catalogue of Galactic AGB in the same luminosity-temperature space. We add details for common or interesting characteristics of RSG, such as multistar system membership, variability, and classification as a runaway. As potential future core-collapse supernova progenitors, we study the ability of the catalogue to inform the Supernova Early Warning System (SNEWS) coincidence network made to automate pointing, and show that for 3D position estimates made possible by neutrinos, the number of progenitor candidates can be significantly reduced, improving our ability to observe the progenitor pre-explosion and the early phases of core-collapse supernovae.

Astronomy & Astrophysics↗

DOE OSTI2023

Evidence for Black Holes in Green Peas from WISE Colors and Variability

We explore the presence of active galactic nuclei (AGNs)/black holes in Green Pea galaxies (GPs), motivated by the presence of high-ionization emission lines such as He ii and [Ne iii] in their optical spectra. In order to identify AGN candidates, we used mid-infrared (MIR) photometric observations from the all-sky Wide-field Infrared Survey Explorer (WISE) mission for a sample of 1004 GPs. Considering only >5σ detections with no contamination from neighboring sources in AllWISE, we select 31 GPs out of 134 as candidate AGNs based on a stringent three-band WISE color diagnostic. Using multi-epoch photometry in W1 and W2 bands based on time-resolved unWISE coadd images, we find two sources exhibiting variability in both the WISE bands among 112 GPs with W1 ≤16 mag and no contamination from neighboring sources in unWISE. These two variable sources were selected as AGNs by the WISE three-band color diagnostic as well. Compared to variable AGN fractions observed among low-mass galaxy samples in previous studies, we find a higher fraction (~1.8%) of MIR variable sources among GPs, which demonstrates the uniqueness and importance of studying these extreme objects. Through this work, we demonstrate that MIR diagnostics are promising tools to select AGNs that may be missed by other selection techniques (including optical emission-line ratios and X-ray emission) in star formation-dominated, low-mass, low-metallicity galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Magnetic braking saturates: evidence from the orbital period distribution of low-mass detached eclipsing binaries from ZTF

ABSTRACT We constrain the orbital period (Porb) distribution of low-mass detached main-sequence eclipsing binaries (EBs) with light-curves from the Zwicky Transient Facility (ZTF), which provides a well-understood selection function and sensitivity to faint stars. At short periods (Porb ≲ 2 d), binaries are predicted to evolve significantly due to magnetic braking (MB), which shrinks orbits and ultimately brings detached binaries into contact. The period distribution is thus a sensitive probe of MB. We find that the intrinsic period distribution of low-mass (0.1 ≲ M1/M⊙ < 0.9) binaries is basically flat (${\rm d}N/{\rm d}P_{\rm orb} \propto P_{\rm orb}^0$) from Porb = 10 d down to the contact limit. This is strongly inconsistent with predictions of classical MB models based on the Skumanich relation, which are widely used in binary evolution calculations and predict ${\rm d}N/{\rm d}P_{\rm orb} \propto P_{\rm orb}^{7/3}$ at short periods. The observed distributions are best reproduced by models in which the magnetic field saturates at short periods with a MB torque that scales roughly as $\dot{J}\propto P_{\rm orb}^{-1}$, as opposed to $\dot{J} \propto P_{\rm orb}^{-3}$ in the standard Skumanich law. We also find no significant difference between the period distributions of binaries containing fully and partially convective stars. Our results confirm that a saturated MB law, which was previously found to describe the spin-down of rapidly rotating isolated M dwarfs, also operates in tidally locked binaries. We advocate using saturated MB models in binary evolution calculations. Our work supports previous suggestions that MB in cataclysmic variables (CVs) is much weaker than assumed in the standard evolutionary model, unless mass transfer leads to significant additional angular momentum loss in CVs.

El-Badry, Kareem (ORCID:0000000268711752)↗

DOE OSTI2026

Chasing Gamma-Ray Signals from Binary Neutron Star Coalescences with the Cherenkov Telescope Array: Prospects and Observing Strategies

The detection of gravitational waves (GWs) from a binary neutron star (BNS) merger by Advanced LIGO and Advanced Virgo (GW170817), together with its electromagnetic counterpart, the short gamma-ray burst GRB 170817A, heralded the birth of multimessenger astronomy. The detection of TeV emission from GRBs motivates follow-up observations with the Cherenkov Telescope Array Observatory (CTAO), which is ideal for detecting such signals due to its unprecedented sensitivity, rapid response, and wide-field survey capabilities. The aim of this work is to evaluate GeV–TeV GW follow-up strategies for CTAO using a multistep simulation pipeline and to estimate the expected rate of joint GW–GRB detections during observing run O5. Using a simulated sample of BNS systems with corresponding GW detections, gamma-ray emission is simulated through phenomenological prescriptions based on the observed population of short GRBs, including off-axis jet scenarios. CTAO observations are simulated to account for instrument response, sky tiling strategies, integration times, and varying observing conditions. Strategies with variable and constant integration times are investigated. We find that, via an optimized follow-up strategy, about 5% of simulated GW-associated short GRBs produce GeV–TeV radiation detectable by CTAO. Detectability is strongly influenced by the jet opening angle and viewing angle, suggesting that even rough estimates of the viewing angle in GW alerts could enhance targeting. This framework motivates future follow-ups of GW-detectable events, including neutron star–black hole mergers, and further supports the development of advanced strategies incorporating galaxy distributions and synergies with future detectors such as the Einstein Telescope.

Abe, S. [University of Tokyo] (ORCID:0000000172503↗

DOE OSTI2025

Do We Know How to Model Reionization?

I compare the power spectra of the radiation fields from two recent sets of fully-coupled simulations that model cosmic reionization: “Cosmic Reionization On Computers” (CROC) and “Thesan”. While both simulations have similar power spectra of the radiation sources, the power spectra of the photoionization rate are significantly different at the same values of cosmic time or the same values of the mean neutral hydrogen fraction. However, the power spectra of the photoionization rate can be matched at large scales for the two simulations when the matching snapshots are allowed to vary independently. I.e., on large scales, the clustering of the radiation field in two simulations evolves similarly, but the exact timing of this evolution is different in different simulations and is not parameterized by an easily interpretable physical quantity like the mean neutral fraction or the mean free path. On small scales, large differences are present and remain partially unexplained. Both CROC and Thesan use the Variable Eddington Tensor approximation for modeling radiative transfer, but adopt different closure relations (optically thin OTVET versus M1). The role of this key difference is tested by using smaller simulations with a new cosmological simulation code that implements both closure relations in a controlled environment (the same hydro, cooling, and gravity solvers and the star formation recipe). In these controlled tests, both the M1 closure and the OTVET ansatz follow the expected behavior from a simple analytical approximation, demonstrating that the differences in the 2-point function of the radiation field induced by the choice of the Eddington tensor are not dominant.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Deep X-Ray Observation of NGC 3221: Everything Everywhere All at Once

We present a comprehensive analysis of 475 ks (438 ks unpublished and 37 ks archival) XMM-Newton/EPIC-pn observations of a nearby, highly inclined, star-forming, luminous infrared galaxy NGC 3221 through spatial, temporal, and spectral information. We confirm the presence of a low-luminosity (presumably Compton-thick) active galactic nucleus (AGN). The 0.4–12 keV luminosity and the hardness ratio of the six ultraluminous X-ray sources previously identified in Chandra data exhibit diverse variability on day scales. The collective emission from unresolved sources exhibits a different day-scale variability. We have also discovered two new predominantly soft (<1 keV) sources. One of these has an enigmatic spectral shape featuring a soft component, which we interpret as a superbubble in NGC 3221, and a variable hard component from a compact object, unresolved from the superbubble. We do not confidently detect any X-ray emission from SN 1961L. The hot gas in the interstellar medium (ISM, out to ±6 kpc from the disk plane) and that in the extraplanar region (6–12 kpc) both require two thermal phases at ∼0.15 keV and ∼0.55 keV. The ∼0.55 keV component is fainter in the ISM than the ∼0.15 keV component, but the emission from the latter falls off more steeply with disk height than the emission from the former. This makes the extraplanar region hotter and less dense than the ISM. The proximity of NGC 3221 and the occurrence of the underluminous AGN offer a unique observing opportunity to study the hot diffuse medium along with nuclear and diskwide point sources.

Das, Sanskriti [Stanford Univ., CA (United States)↗

DOE OSTI2026

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗