Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Astronomy data analysis”

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 577 records · Page 32

Constraining Hamiltonians from chiral effective field theory with neutron-star data

Multi-messenger observations of neutron stars (NSs) and their mergers have placed strong constraints on the dense-matter equation of state (EOS). The EOS, in turn, depends on microscopic nuclear interactions that are described by nuclear Hamiltonians. These Hamiltonians are commonly derived within chiral effective field theory (EFT). Ideally, multi-messenger observations of NSs could be used to directly inform our understanding of EFT interactions, but such a direct inference necessitates millions of model evaluations. This is computationally prohibitive because each evaluation requires us to calculate the EOS from a Hamiltonian by solving the quantum many-body problem with methods such as auxiliary-field diffusion Monte Carlo (AFDMC), which provides very accurate and precise solutions but at a significant computational cost. Additionally, we need to solve the stellar structure equations for each EOS which further slows down each model evaluation by a few seconds. In this work, we combine emulators for AFDMC calculations of neutron matter, built using parametric matrix models, and for the stellar structure equations, built using multilayer perceptron neural networks, with the PyCBC data-analysis framework to enable a direct inference of coupling constants in an EFT Hamiltonian using multi-messenger observations of NSs. We find that astrophysical data can provide informative constraints on two-nucleon couplings despite the high densities probed in NS interiors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Search for topological defect dark matter with a global network of optical magnetometers

Ultralight bosons such as axion-like particles are viable candidates for dark matter. They can form stable, macroscopic field configurations in the form of topological defects that could concentrate the dark matter density into many distinct, compact spatial regions that are small compared with the Galaxy but much larger than the Earth. Here we report the results of the search for transient signals from the domain walls of axion-like particles by using the global network of optical magnetometers for exotic (GNOME) physics searches. We search the data, consisting of correlated measurements from optical atomic magnetometers located in laboratories all over the world, for patterns of signals propagating through the network consistent with domain walls. The analysis of these data from a continuous month-long operation of GNOME finds no statistically significant signals, thus placing experimental constraints on such dark matter scenarios.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Interacting dark sector within ETHOS: Cosmological constraints from SPT cluster abundance with DES and HST weak lensing data

We use galaxy cluster abundance measurements from the South Pole Telescope enhanced by multicomponent matched filter confirmation and complemented with mass information obtained using weak-lensing data from Dark Energy Survey Year 3 (DES Y3) and targeted Hubble Space Telescope observations for probing deviations from the cold dark matter paradigm. Concretely, we consider a class of dark sector models featuring interactions between dark matter (DM) and a dark radiation (DR) component within the framework of the effective theory of structure formation (ETHOS). We focus on scenarios that lead to power suppression over a wide range of scales, and thus can be tested with data sensitive to large scales, as realized, for example, for DM–DR interactions following from an unbroken non-Abelian 𝑆⁢𝑈⁡(𝑁) gauge theory (interaction rate with power-law index 𝑛 = 0 within the ETHOS parametrization). Cluster abundance measurements are mostly sensitive to the amount of DR interacting with DM, parametrized by the ratio of DR temperature to the cosmic microwave background (CMB) temperature, 𝜉 DR = 𝑇 DR /𝑇 CMB . We find an upper limit 𝜉 DR < 17% at 95% credibility. When the cluster data are combined with Planck 2018 CMB data along with baryon acoustic oscillation (BAO) measurements we find 𝜉 DR < 10%, corresponding to a limit on the abundance of interacting DR that is around 3 times tighter than that from CMB + BAO data alone. We also discuss the complementarity of weak lensing informed cluster abundance studies with probes sensitive to smaller scales, explore the impact on our analysis of massive neutrinos, and comment on a slight preference for the presence of a nonzero interacting DR abundance, which enables a physical solution to the 𝑆 8 tension.

79 ASTRONOMY AND ASTROPHYSICS↗

First Results on Dark Matter Substructure from Astrometric Weak Lensing

Low-mass structures of dark matter (DM) are expected to be entirely devoid of light-emitting regions and baryons. Precisely because of this lack of baryonic feedback, small-scale substructures of the Milky Way are a relatively pristine testing ground for discovering aspects of DM microphysics and primordial fluctuations on subgalactic scales. In this work, we report results from the first search for Galactic DM subhalos with time-domain astrometric weak gravitational lensing. The analysis is based on a matched-filter template of local lensing corrections to the proper motion of stars in the Magellanic Clouds. We describe a data analysis pipeline detailing sample selection, background subtraction, and handling outliers and other systematics. For tentative candidate lenses, we identify a signature based on an anomalous parallax template that can unequivocally confirm the presence of a DM lens, opening up prospects for robust discovery potential with full time-series data. We present our constraints on substructure fraction f l ≲5 at 90% CL (and f l ≲2 at 50% CL) for compact lenses with radii r l <1pc, with best sensitivity reached for lens masses M l around 10 7 -10 8 M ⊙ . Parametric improvements are expected with future astrometric data sets; by end of mission, Gaia could reach f l ≲10 -3 for these massive point-like objects, and be sensitive to lighter and/or more extended subhalos for O(1) substructure fractions.

79 ASTRONOMY AND ASTROPHYSICS↗

The Milky Way Revealed by Variable Stars. I. Sample Selection of RR Lyrae Stars and Evidence for Merger History

In order to study the Milky Way, RR Lyrae (RRL) variable stars identified by Gaia, ASAS-SN, and ZTF sky survey projects have been analyzed as tracers in this work. Photometric and spectroscopic information of 3417 RRLs including proper motions, radial velocity, and metallcity are obtained from observational data of Gaia, LAMOST, GALAH, APOGEE, and RAVE. Precise distances of RRLs with typical uncertainties less than 3% are derived by using a recent comprehensive period–luminosity–metallicity relation. Our results from kinematical and chemical analysis provide important clues for the assembly history of the Milky Way, especially for the Gaia–Sausage ancient merger. The kinematical and chemical trends found in this work are consistent with those of recent simulations that indicated that the Gaia–Sausage merger had a dual origin in the Galactic thick disk and halo. As recent similar works have found, the halo RRL sample in this work contains a subset of radially biased orbits besides a more isotropic component. This higher orbital anisotropy component amounts to β ≃ 0.8, and it contributes between 42% and 83% of the halo RRLs at 4 < R( kpc) < 20.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling Dust Production, Growth, and Destruction in Reionization-era Galaxies with the CROC Simulations: Methods and Parameter Exploration

We introduce a model for the explicit evolution of interstellar dust in a cosmological galaxy formation simulation. We post-process a simulation from the Cosmic Reionization on Computers project, integrating an ordinary differential equation for the evolution of the dust-to-gas ratio along pathlines in the simulation sampled with a tracer particle technique. This model incorporates the effects of dust grain production in asymptotic giant branch star winds and supernovae, grain growth due to the accretion of heavy elements from the gas phase of the interstellar medium, and grain destruction due to thermal sputtering in the high-temperature gas of supernova remnants due to thermal sputtering. A main conclusion of our analysis is the importance of a carefully chosen dust destruction model, for which different reasonable parameterizations can predict very different values at the ∼100 pc resolution of the interstellar medium (ISM) in our simulations. We run this dust model on the single most massive galaxy in a 10 h$^{−1}$ comoving megaparsec box, which attains a stellar mass of ∼2 × 10$^{9}$ M $_{⊙}$ by z = 5. We find that the model is capable of reproducing dust masses and dust-sensitive observable quantities broadly consistent with existing data from high-redshift galaxies. The total dust mass in the simulated galaxy is somewhat sensitive to parameter choices for the dust model, especially the timescale for grain growth due to accretion in the ISM. Consequently, observable quantity observations that can constrain galaxy dust masses at these epochs are potentially useful for placing constraints on dust physics.

79 ASTRONOMY AND ASTROPHYSICS↗

The COMBS Survey - II. Distinguishing the metal-poor bulge from the halo interlopers

The metal-poor stars in the bulge are important relics of the Milky Way’s formation history, as simulations predict that they are some of the oldest stars in the Galaxy. In order to determine if they are truly ancient stars, we must understand their origins. Currently, it is unclear if the metal-poor stars in the bulge ([Fe/H] < -1 dex) are merely halo interlopers, a unique accreted population, part of the boxy/peanut-shaped bulge, or a classical bulge population. In this work, we use spectra from the VLT/FLAMES spectrograph to obtain metallicity estimates using the Ca-II triplet of 473 bulge stars (187 of which have [Fe/H] < -1 dex), targeted using SkyMapper photometry. We also use Gaia DR2 data to infer the Galactic positions and velocities along with orbital properties for 523 stars. We employ a probabilistic orbit analysis and find that about half of our sample has a >50 percent probability of being bound to the bulge, and half are halo interlopers. We also see that the occurrence rate of halo interlopers increases steadily with decreasing metallicity across the full range of our sample (-3 < [Fe/H] < 0.5). Our examination of the kinematics of the confined compared to the unbound stars indicates the metal-poor bulge comprises at least two populations; those confined to the boxy/peanut bulge and halo stars passing through the inner galaxy. We conclude that an orbital analysis approach, as we have employed, is important to understand the composite nature of the metal-poor stars in the inner region.

79 ASTRONOMY AND ASTROPHYSICS↗

Reconstruction with velocities

ABSTRACT Reconstruction is becoming a crucial procedure of galaxy clustering analysis for future spectroscopic redshift surveys to obtain subper cent level measurement of the baryon acoustic oscillation scale. Most reconstruction algorithms rely on an estimation of the displacement field from the observed galaxy distribution. However, the displacement reconstruction degrades near the survey boundary due to incomplete data and the boundary effects extend to ${\sim}100\, \mathrm{Mpc}/h$ within the interior of the survey volume. We study the possibility of using radial velocities measured from the cosmic microwave background observation through the kinematic Sunyaev–Zeldovich effect to improve performance near the boundary. We find that the boundary effect can be reduced to ${\sim}30-40\, \mathrm{Mpc}/h$ with the velocity information from Simons Observatory. This is especially helpful for dense low redshift surveys where the volume is relatively small and a large fraction of total volume is affected by the boundary.

79 ASTRONOMY AND ASTROPHYSICS↗

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential for specialized applications in segmented inverse beta decay (IBD) neutrino detectors, astronomy, machine learning, and more. Although this method may easily extend to 3D scalar fields, our focus here is on 2D real-valued fields as it directly applies to directionality. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Data Analysis, Statistics and Probability (physics↗

Impact of heavy quark and quarkonium data on nuclear gluon PDFs

A clear understanding of nuclear parton distribution functions (nPDFs) plays a crucial role in the interpretation of collider data taken at the Relativistic Heavy Ion Collider, the Large Hadron Collider (LHC), and in the near future at the Electron-Ion Collider. Even with the recent inclusions of vector boson and light meson production data, the uncertainty of the gluon PDF remains substantial and limits the interpretation of heavy ion collision data. To obtain new constraints on the nuclear gluon PDF, we extend our recent nCTEQ15WZ+SIH analysis to inclusive quarkonium and open heavy-flavor meson production data from the LHC. This vast new data set covers a wide kinematic range and puts strong constraints on the nuclear gluon PDF down to x ≲ 10 -5 . The theoretical predictions for these data sets are obtained from a data-driven approach, where proton-proton data are used to determine effective scattering matrix elements. This approach is validated with detailed comparisons to existing next-to-leading order calculations in nonrelativistic QCD for quarkonia and in the general-mass variable-flavor-number scheme for the open heavy-flavored mesons. In addition, the uncertainties from the data-driven approach are determined using the Hessian method and accounted for in the PDF fits. This extension of our previous analyses represents an important step toward the next generation of PDFs not only by including new data sets, but also by exploring new methods for future analyses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The luminosity function of ringed galaxies

ABSTRACT We perform an analysis of the luminosity functions (LFs) of two types of ringed galaxies – polar-ring galaxies and collisional ring galaxies – using data from the Sloan Digital Sky Survey (SDSS). Both classes of galaxies were formed as a result of interaction with their environment and they are very rare objects. We constructed LFs of galaxies by different methods and found their approximations by the Schechter function. The luminosity functions of both types of galaxies show a systematic fall-off at low luminosities. The polar structures around bright (Mr ≤ −20m) and red (g − r > +0.8) galaxies are about twice as common as around blue ones. The LF of collisional rings is shifted towards brighter luminosities compared to polar-ring galaxies. We analysed the published data on the ringed galaxies in several deep fields and confirmed the increase in their volume density with redshift: up to z ∼ 1 their density grows as (1 + z)m, where m ≳ 5.

79 ASTRONOMY AND ASTROPHYSICS↗

Systematic study of projection biases in the weak lensing analysis of cosmic shear and the combination of galaxy clustering and galaxy-galaxy lensing

This paper presents the results of a systematic study of projection biases in the weak lensing analysis of cosmic shear and the combination of galaxy clustering and galaxy-galaxy lensing using data collected during the first year of running the Dark Energy Survey experiment. The study uses Lambda cold dark matter ( Λ CDM ) as the cosmological model and two-point correlation functions for the weak lensing (WL) analysis. The results in this paper show that, independent of the WL analysis, projection biases of more than 1 σ exist and are a function of the position of the true values of the parameters h , n s , Ω b , and Ω ν h 2 with respect to their prior probabilities. For cosmic shear, and the combination of galaxy clustering and galaxy-galaxy lensing, this study shows that the coverage probability of the 68.27% credible intervals ranges from as high as 93% to as low as 16% and that these credible intervals are inflated, on average, by 29% for cosmic shear and 20% for the combination of galaxy clustering and galaxy-galaxy lensing. The results of the study also show that, in six out of nine tested cases, the reduction in error bars obtained by transforming credible intervals into confidence intervals is equivalent to an increase in the amount of data by a factor of 3.

79 ASTRONOMY AND ASTROPHYSICS↗

New insights into supradense matter from dissecting scaled stellar structure equations

The strong-field gravity in general relativity (GR) realized in neutron stars (NSs) renders the equation of state (EOS) P(ε) of supradense neutron star matter to be essentially nonlinear and refines the upper bound for Φ ≡ P/ε to be much smaller than the special relativity (SR) requirement with linear EOSs, where P and ε are respectively the pressure and energy density of the system considered. Specifically, a tight bound Φ ≲ 0.374 is obtained by perturbatively anatomizing the intrinsic structures of the scaled Tolman–Oppenheimer–Volkoff (TOV) equations without using any input nuclear EOS. New insights gained from this novel analysis provide EOS-model-independent constraints on the properties (e.g., density profiles of the sound speed squared s 2 = dP/dε and trace anomaly Δ = 1/3 – Φ) of cold supradense matter in NS cores. Using the gravity-matter duality in theories describing NSs, we investigate the impact of gravity on supradense matter EOS in NSs. In particular, we show that the NS mass M NS , radius R, and compactness ξ ≡ M NS /R scale with certain combinations of its central pressure and energy density (encapsulating its central EOS). Thus, observational data on these properties of NSs can straightforwardly constrain NS central EOSs without relying on any specific nuclear EOS model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Deep Search for Gamma-Ray Emission from the Accreting X-Ray Pulsar 1A 0535+262

Binary systems are a well-established subclass of gamma-ray sources. The high mass X-ray binary pulsar 1A 0535+262 has been considered to be a possible gamma-ray emitter for a long time, although former gamma-ray searches using the Fermi Large Area Telescope (LAT) and VERITAS data resulted in upper limits only. We aim at a deep search for gamma-ray emission and pulsations from 1A 0535+262 using more than 13 yr of Fermi-LAT data. The analysis was performed for both the whole Fermi-LAT data set, as well as for the X-ray outbursts that 1A 0535+262 has experienced since the launch of Fermi. Various X-ray observations have been used to generate the ephemeris for the pulsation search. We also investigate the long-term gamma-ray flux variability and perform orbital phase-resolved analysis for the outbursts. We did not detect any steady or pulsed gamma-ray emission from 1A 0535+262 during the whole Fermi-LAT mission span or its X-ray outbursts. We thus derived the deepest gamma-ray luminosity upper limits to date at the 95% confidence level to be around (2.3–4.7) × 10 32 erg s -1 depending on different spectral indices assumed, which results in a ratio of L γ to L X (2–150 keV) being (1.9–3.9) × 10 -6 .

79 ASTRONOMY AND ASTROPHYSICS↗

New dark matter analysis of milky way dwarf satellite galaxies with madhatv2

We obtain bounds on dark matter annihilation using 14 years of publicly available Fermi-LAT data from a set of 54 dwarf spheroidal galaxies, using spectral information from 16 energy bins. We perform this analysis using our updated and publicly available code , which can be used to test a variety of models for dark matter particle physics and astrophysics in an accessible manner. In particular, we note that including Carina III in the analysis strengthens constraints on s -wave annihilation into two-body Standard Model final states by a factor of ∼ 3 but broadens the error on the constraint due to the large uncertainty of its J -factor. Our findings illustrate the importance of verifying if Carina III is in fact a dwarf spheroidal galaxy and measuring more precisely its J -factor. More generally, they highlight the significance of forthcoming discoveries of nearby ultrafaint dwarfs for dark matter indirect detection. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Identification of Blue Horizontal Branch Stars with Multimodal Fusion

Blue Horizontal Branch stars (BHBs) are ideal tracers to probe the global structure of the milky Way (MW), and the increased size of the BHB star sample could be helpful to accurately calculate the MW's enclosed mass and kinematics. Large survey telescopes have produced an increasing number of astronomical images and spectra. However, traditional methods of identifying BHBs are limited in dealing with the large scale of astronomical data. A fast and efficient way of identifying BHBs can provide a more significant sample for further analysis and research. Therefore, in order to fully use the various data observed and further improve the identification accuracy of BHBs, we have innovatively proposed and implemented a Bi-level attention mechanism-based Transformer multimodal fusion model, called Bi-level Attention in the Transformer with Multimodality (BATMM). The model consists of a spectrum encoder, an image encoder, and a Transformer multimodal fusion module. The Transformer enables the effective fusion of data from two modalities, namely image and spectrum, by using the proposed Bi-level attention mechanism, including cross-attention and self-attention. As a result, the information from the different modalities complements each other, thus improving the accuracy of the identification of BHBs. The experimental results show that the F1 score of the proposed BATMM is 94.78%, which is 21.77% and 2.76% higher than the image and spectral unimodality, respectively. It is therefore demonstrated that higher identification accuracy of BHBs can be achieved by means of using data from multiple modalities and employing an efficient data fusion strategy.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct nonparametric multimessenger constraints on the equation of state of cold dense nuclear matter

We utilize the now substantial amount of astrophysical observations of neutron stars (NSs), along with perturbative quantum chromodynamics (pQCD) calculations at high density, to directly constrain the NS equation of state (EOS). To this end, we construct nonparametric EOS priors by using Gaussian processes trained on 75 EOSs, which include models with either hadrons, hyperons, or quarks at high densities. We create a prior using the full EOS sample (model agnostic), and one prior for each EOS family to test model discrimination. We introduce a novel inference approach, which allows the simultaneous sampling of intrinsic and extrinsic parameters of binary NS mergers, as well as a nonparametric equation of state. We showcase this method in a Bayesian updating scheme by first performing a complete analysis of the binary NS merger event GW170817 with minimal assumptions, and sequentially adding information from x-ray and radio NS observations, along with pQCD calculations. Besides providing standard constraints, such as the pressure at twice nuclear saturation density 𝑝⁡(2⁢𝜌 sat ) = 4.3$^{+0.6}_{−0.6}$ × 10 34 dyne/cm 2 , at 95% confidence level, for the model agnostic prior, our methodology shows how the choice of EOS families used in conditioning changes the inferred astrophysical properties of the EOS, namely tidal deformability and maximum supported NS mass. We find hyperonic priors predicting higher tidal deformabilities for a 1.4⁢𝑀 ⊙ NS, and hadronic priors being preferred by the considered astrophysical data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Deciphering Fingerprints of Stellar Nucleosynthesis Through Nuclear Reaction Rate Measurements and Isotopic Analyses of Stardust

Microscopic stardust grains, found in trace amounts in primitive meteorites, afford a unique opportunity to study stellar nucleosynthesis in the laboratory by measuring their heavy element isotopic compositions. Applying LLNL’s state-of-the art material characterization capabilities to these micrometer-sized grains yields large new datasets with unprecedented precision. However, the astrophysical models necessary to interpret the isotopic record require improved nuclear physics data to constrain stellar nucleosynthesis conditions. For example, the branch point in the s-process path at 95 Zr controls the relative amounts of heavier nuclides but is poorly understood. We conducted a coordinated effort across disciplines to combine a new determination of the neutron capture cross section of 95 Zr with the laboratory analysis of hundreds of stardust grains to yield a view into stellar interiors. We developed new methods to quantify multiple isotopic systems simultaneously (e.g., Zr, Ba, and W) in stardust grains and improved the determination of the neutron capture cross section of unstable 95 Zr via the surrogate reaction method. This interdisciplinary approach is allowing us to fill gaps in our understanding of the environments in which elements are produced while enhancing LLNL’s experimental and theoretical capabilities relevant to nuclear threat reduction missions.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗