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 469 records · Page 26

NuSTAR broad-band X-ray observation campaign of energetic pulsar wind nebulae in synergy with VERITAS, HAWC and Fermi gamma-ray telescopes

We report recent progress on the on-going NuSTAR observation campaign of 8 TeV-detected middle-aged pulsar wind nebulae (PWNe). This campaign constitutes a major part of our NuSTAR study of some of the most energetic TeV sources in our Galaxy detected by VERITAS and HAWC. NuSTAR is the first focusing X-ray telescope operating above 10 keV in space with sub-arcminute angular resolution. Broad-band X-ray imaging and spectroscopy data, obtained by NuSTAR, allow us to probe sub-PeV electron populations through detecting synchrotron X-ray radiation. Our targets include PeVatron candidates detected by HAWC, the Boomerang nebula, PWNe crushed by supernova remnant shocks (or relic PWNe) and G0.9‚0.1 in the Galactic Center. Combined with Fermi-LAT data and available TeV data, we aim to provide a complete, multi-wavelength view of a diverse class of middle-aged (~10–100 kyrs old) PWNe. Our NuSTAR analysis detected hard X-ray emission from the Eel and Boomerang PWNe and characterized their broad-band X-ray spectra most accurately. We plan to apply both time-evolution and multi-zone PWNe models to multi-wavelength spectral energy distribution (SED) data over the radio, X-ray, GeV and TeV bands. In this proceeding, we will review our observation campaign and discuss the preliminary results for some PWNe.

79 ASTRONOMY AND ASTROPHYSICS↗

Large-scale Structure of the Heliospheric Current Sheets within the Heliosheath Inferred from Voyager 2 Observations

In this paper, we conducted an analysis of the heliospheric current sheets (HCSs) in the heliosheath (HS), utilizing observations by Voyager 2 between 2008 and 2018. Employing rigorous criteria, we identified a total of 34 HCSs that indicate significant changes in magnetic polarity. These occurrences were more prevalent during solar maximum periods when the HCS expanded to higher latitudes, coinciding with the spacecraft positioned at an average latitude of –31° from the solar equator. We determined certain features of the large-scale structures of the identified HCSs. Most importantly, employing two distinct methods indicates that the thickness of the HCSs within the HS ranges from ~0.003 to ~0.4 au with an average thickness of ~0.03 to ~0.1 au, depending on methods of event selection and fitting techniques. This thickness surpasses that known near 1 au or other heliospheric distances. It is also notably thicker than the typical proton inertial length, implying unfavorable conditions for magnetic reconnection. Additionally, our analysis reveals a frequent tilt of HCS planes relative to the solar equatorial plane by a varying angle up to several tens of degrees, likely implying a common occurrence of a warped structure of the HCS within the HS. Longitudinally, the HCS planes closely align with the Parker spiral field direction expected in the HS. Finally, for a large fraction of the identified HCS events, the HCS planes are likely characterized by a rotational discontinuity. These findings are valid within the limits of the 1 hr resolution data used in this study.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of BOSS galaxy data with weighted skew-spectra

We present the first application of the weighted skew-spectra to analyze non-Gaussian information in galaxy survey data. Using the tree-level galaxy skew-spectra together with the one-loop power spectrum multipoles, we analyze the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) galaxy clustering data, and target our search towards the equilateral bispectrum shape of primordial non-Gaussianity. We use the Effective Field Theory model for the galaxy power spectrum and bispectrum, and account for systematic effects, such as the survey geometry. From our likelihood analysis, we find f NL equil = -34 +296 -334 at 68% CL, consistent with previous works, while systematic errors from our treatment of the survey geometry lead to an unreliable estimation of f NL ortho . As a result, we further constrain the bias and counterterm parameters, while keeping the cosmology fixed to Planck 2018 values. As a check, we also validate our analysis pipeline using theNseries simulation suite.

79 ASTRONOMY AND ASTROPHYSICS↗

TARDIS. II. Synergistic Density Reconstruction from Lyα Forest and Spectroscopic Galaxy Surveys with Applications to Protoclusters and the Cosmic Web

Looking at this work, we expand upon the Tomographic Absorption Reconstruction and Density Inference Scheme (TARDIS) in order to include multiple tracers while reconstructing matter density fields at Cosmic Noon (z ~ 2–3). In particular, we jointly reconstruct the underlying density field from simulated Lyα forest observations at z ~ 2.5 and an overlapping galaxy survey. We find that these data are synergistic, with the Lyα forest providing reconstruction of low-density regions and galaxy surveys tracing the density peaks. We find a more accurate power spectra reconstruction going to higher scales when fitting these two data sets simultaneously than when using either one individually. When applied to cosmic web analysis, we find that performing the joint analysis is equivalent to an Lγα survey with significantly increased sight-line spacing. Because we reconstruct the velocity field and matter field jointly, we demonstrate the ability to evolve the mock observed volume further to z = 0, allowing us to create a rigorous definition of a “protocluster” as regions that will evolve into clusters. We apply our reconstructions to study protocluster structure and evolution, finding for realistic survey parameters that we can provide accurate mass estimates of the z ≈ 2 structures and their z = 0 fate.

79 ASTRONOMY AND ASTROPHYSICS↗

Cosmological constraints from H ii starburst galaxy apparent magnitude and other cosmological measurements

ABSTRACT We use H ii starburst galaxy apparent magnitude measurements to constrain cosmological parameters in six cosmological models. A joint analysis of H ii galaxy, quasar angular size, baryon acoustic oscillations peak length scale, and Hubble parameter measurements result in relatively model-independent and restrictive estimates of the current values of the non-relativistic matter density parameter $\Omega _{\rm m_0}$ and the Hubble constant H0. These estimates favour a 2.0–3.4σ (depending on cosmological model) lower H0 than what is measured from the local expansion rate. The combined data are consistent with dark energy being a cosmological constant and with flat spatial hypersurfaces, but do not strongly rule out mild dark energy dynamics or slightly non-flat spatial geometries.

79 ASTRONOMY AND ASTROPHYSICS↗

Tree-based algorithms for weakly supervised anomaly detection

Weakly supervised methods have emerged as a powerful tool for model-agnostic anomaly detection at the Large Hadron Collider (LHC). While these methods have shown remarkable performance on specific signatures such as dijet resonances, their application in a more model-agnostic manner requires dealing with a larger number of potentially noisy input features. In this paper, we show that using boosted decision trees as classifiers in weakly supervised anomaly detection gives superior performance compared to deep neural networks. Boosted decision trees are well known for their effectiveness in tabular data analysis. Our results show that they not only offer significantly faster training and evaluation times, but they are also robust to a large number of noisy input features. By using advanced gradient boosted decision trees in combination with ensembling techniques and an extended set of features, we significantly improve the performance of weakly supervised methods for anomaly detection at the LHC. This advance is a crucial step toward a more model-agnostic search for new physics. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Singling out modified gravity parameters and data sets reveals a dichotomy between Planck and lensing

ABSTRACT An important route to testing general relativity (GR) at cosmological scales is usually done by constraining modified gravity (MG) parameters added to the Einstein perturbed equations. Most studies have analysed so far constraints on pairs of MG parameters, but here, we explore constraints on one parameter at a time while fixing the other at its GR value. This allows us to analyse various models while benefiting from a stronger constraining power from the data. We also explore which specific data sets are in tension with GR. We find that models with (μ = 1, η) and (μ, η = 1) exhibit a 3.9σ and 3.8σ departure from GR when using Planck18 + Supernovae type Ia (SNe) + Baryon acoustic oscillation (BAO), while (μ, η) shows a tension of 3.4σ. We find no tension with GR for models with the MG parameter Σ fixed to its GR value. Using a Bayesian model selection analysis, we find that some one-parameter MG models are moderately favoured over ΛCDM when using all data set combinations except Planck cosmic microwave background lensing and dark energy survey data. Namely, Planck18 shows a moderate tension with GR that only increases when adding any combination of redshift space distortion, SNe, or BAO. However, adding lensing diminishes or removes these tensions, which can be attributed to the ability of lensing in constraining the MG parameter Σ. The two overall groups of data sets are found to have a dichotomy when performing consistency tests with GR, which may be due to systematic effects, lack of constraining power, or modelling. These findings warrant further investigation using more precise data from ongoing and future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

guppy i : a code for reducing the storage requirements of cosmological simulations

ABSTRACT As cosmological simulations have grown in size, the permanent storage requirements of their particle data have also grown. Even modest simulations present a major logistical challenge for the groups which run these boxes and researchers without access to high performance computing facilities often need to restrict their analysis to lower quality data. In this paper, we present guppy, a compression algorithm and code base tailored to reduce the sizes of dark matter-only cosmological simulations by approximately an order of magnitude. guppy is a ‘lossy’ algorithm, meaning that it injects a small amount of controlled and uncorrelated noise into particle properties. We perform extensive tests on the impact that this noise has on the internal structure of dark matter haloes, and identify conservative accuracy limits which ensure that compression has no practical impact on single-snapshot halo properties, profiles, and abundances. We also release functional prototype libraries in C, Python, and Go for reading and creating guppy data.

79 ASTRONOMY AND ASTROPHYSICS↗

A Joint Roman Space Telescope and Rubin Observatory Synthetic Wide-Field Imaging Survey

We present and validate 20 deg(2) of overlapping synthetic imaging surveys representing the full depth of the Nancy Grace Roman Space Telescope High-Latitude Imaging Survey (HLIS) and 5 yr of observations of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The two synthetic surveys are summarized, with reference to the existing 300 deg(2) of LSST simulated imaging produced as part of Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2). Both synthetic surveys observe the same simulated DESC DC2 universe. For the synthetic Roman survey, we simulate for the first time fully chromatic images along with the detailed physics of the Sensor Chip Assemblies derived from lab measurements using the flight detectors. The simulated imaging and resulting pixel-level measurements of photometric properties of objects span a wavelength range of similar to 0.3-2.0 mu m. We also describe updates to the Roman simulation pipeline, changes in how astrophysical objects are simulated relative to the original DC2 simulations, and the resulting simulated Roman data products. We use these simulations to explore the relative fraction of unrecognized blends in LSST images, finding that 20-30 per cent of objects identified in LSST images with i-band magnitudes brighter than 25 can be identified as multiple objects in Roman images. These simulations provide a unique testing ground for the development and validation of joint pixel-level analysis techniques of ground- and space-based imaging data sets in the second half of the 2020s - in particular the case of joint Roman-LSST analyses.

gravitational lensing↗

TOI-1235 b: A Keystone Super-Earth for Testing Radius Valley Emergence Models around Early M Dwarfs

Small planets on close-in orbits tend to exhibit envelope mass fractions of either effectively zero or up to a few percent depending on their size and orbital period. Models of thermally driven atmospheric mass loss and of terrestrial planet formation in a gas-poor environment make distinct predictions regarding the location of this rocky/nonrocky transition in period–radius space. Here we present the confirmation of TOI-1235 b (P = 3.44 days, r{sub p}=1.738{sub −0.076}{sup +0.087} R{sub ⊕}), a planet whose size and period are intermediate between the competing model predictions, thus making the system an important test case for emergence models of the rocky/nonrocky transition around early M dwarfs (R {sub s} = 0.630 ± 0.015 R{sub ⊙}, M {sub s} = 0.640 ± 0.016 M{sub ⊙}). We confirm the TESS planet discovery using reconnaissance spectroscopy, ground-based photometry, high-resolution imaging, and a set of 38 precise radial velocities (RVs) from HARPS-N and HIRES. We measure a planet mass of 6.91{sub −0.85}{sup +0.75} M{sub ⊕}, which implies an iron core mass fraction of 20{sub −12}{sup +15}% in the absence of a gaseous envelope. The bulk composition of TOI-1235 b is therefore consistent with being Earth-like, and we constrain an H/He envelope mass fraction to be <0.5% at 90% confidence. Our results are consistent with model predictions from thermally driven atmospheric mass loss but not with gas-poor formation, suggesting that the former class of processes remains efficient at sculpting close-in planets around early M dwarfs. Our RV analysis also reveals a strong periodicity close to the first harmonic of the photometrically determined stellar rotation period that we treat as stellar activity, despite other lines of evidence favoring a planetary origin (P=21.8{sub −0.8}{sup +0.9} days, m{sub p}sin i=13.0{sub −5.3}{sup +3.8} M{sub ⊕}) that cannot be firmly ruled out by our data.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiscale Flow for robust and optimal cosmological analysis

We propose Multiscale Flow, a generative Normalizing Flow that creates samples and models the field-level likelihood of two-dimensional cosmological data such as weak lensing. Multiscale Flow uses hierarchical decomposition of cosmological fields via a wavelet basis and then models different wavelet components separately as Normalizing Flows. The log-likelihood of the original cosmological field can be recovered by summing over the log-likelihood of each wavelet term. This decomposition allows us to separate the information from different scales and identify distribution shifts in the data such as unknown scale-dependent systematics. The resulting likelihood analysis can not only identify these types of systematics, but can also be made optimal, in the sense that the Multiscale Flow can learn the full likelihood at the field without any dimensionality reduction. We apply Multiscale Flow to weak lensing mock datasets for cosmological inference and show that it significantly outperforms traditional summary statistics such as power spectrum and peak counts, as well as machine learning–based summary statistics such as scattering transform and convolutional neural networks. We further show that Multiscale Flow is able to identify distribution shifts not in the training data such as baryonic effects. Finally, we demonstrate that Multiscale Flow can be used to generate realistic samples of weak lensing data.

79 ASTRONOMY AND ASTROPHYSICS↗

TIC 278956474: Two Close Binaries in One Young Quadruple System Identified by TESS

We have identified a quadruple system with two close eclipsing binaries in Transiting Exoplanet Survey Satellite (TESS) data. The object is unresolved in Gaia and appears as a single source at parallax 1.08 ± 0.01 mas. Both binaries have observable primary and secondary eclipses and were monitored throughout TESS Cycle 1 (sectors 1–13), falling within the TESS Continuous Viewing Zone. In one eclipsing binary (P = 5.488 days), the smaller star is completely occluded by the larger star during the secondary eclipse; in the other (P = 5.674 days) both eclipses are grazing. Using these data, spectroscopy, speckle photometry, spectral energy distribution analysis, and evolutionary stellar tracks, we have constrained the masses and radii of the four stars in the two eclipsing binaries. The Li i equivalent width indicates an age of 10–50 Myr and, with an outer period of 858{sub −5}{sup +7} days, our analysis indicates this is one of the most compact young 2 + 2 quadruple systems known.

79 ASTRONOMY AND ASTROPHYSICS↗

High-velocity CP2 stars in the Galactic halo

The second subclass of chemically peculiar stars, the CP2 stars, are early-type stars exhibiting anomalous abundances with variable line strengths possibly also accompanied by photometric variability that typically belong to the Galactic disk. However, a small fraction of these objects were recently found to be located far from the Galactic plane and are thought to belong to the Galactic halo, which is unexpected for this class of objects. Our study investigates the halo membership of the proposed halo CP2 star candidates based on their velocities and Galactic orbits, to determine their points of origin. In addition, we searched for further halo CP2 star candidates using an as-yet-unpublished sample of CP stars. Our analysis relied on the astrometry from the early third Gaia data release and radial velocities based on the spectroscopy from LAMOST and SDSS, which was also employed in spectral classification. The light variability of the confirmed CP2 stars in our sample was analyzed using data from the ZTF and ATLAS surveys. After filtering our initial sample using kinematic and spectroscopic criteria, we identified six CP2 stars with kinematical properties consistent with a halo membership. The orbits of these stars are in agreement with an origin in the Galactic disk where they were probably ejected through dynamical interactions or in the binary supernova scenario, making them the first runaway CP2 stars known.

79 ASTRONOMY AND ASTROPHYSICS↗

Exploring modified gravity: constraints on the μ and Σ parametrization with WMAP , ACT, and SPT

ABSTRACT The cosmic acceleration problem remains one of the most significant challenges in cosmology. One of the proposed solutions to this problem is the modification of gravity on large scales. In this paper, we explore the well-known μ–Σ parametrization scenarios and confront them with observational data, including the cosmic microwave background (CMB) radiation from the Wilkinson Microwave Anisotropy Probe (WMAP), Atacama Cosmology Telescope (ACT), and South Pole Telescope (SPT), as well as large-scale structure data from the Sloan Digital Sky Survey (SDSS; baryon acoustic oscillation + redshift-space distortion) and Pantheon supernova (SN) catalogue. We employ a Bayesian framework to constrain the model parameters and discuss the implications of our results on the viability of modified gravity theories. Our analysis reveals the strengths and limitations of the μ–Σ parametrization and provides valuable insights into the nature of gravity on cosmological scales. From the joint analysis of the ACT + WMAP + SDSS + SN, we find μ0 − 1 = 0.02 ± 0.19 and Σ0 − 1 = 0.021 ± 0.068 at 68 per cent confidence level (CL). In light of the SPT + WMAP + SDSS + SN, we find μ0 − 1 = 0.07 ± 0.18 and $\Sigma _0 -1 = -0.009^{+0.078}_{-0.11}$ at 68 per cent CL. In all the analyses carried out, we do not find any deviations from the theory of general relativity. Our results represent an observational update on the well-known μ–Σ parametrization in view of current CMB data, independent of and competitive with the constraints obtained with the Planck data.

79 ASTRONOMY AND ASTROPHYSICS↗

Summary of NDC Capacity Building Workshop and Regional Seismic Travel Time (RSTT) in combination with Data Sharing and Integration Training

The NNSA Seismic Cooperation Program (SCP) sponsored Stephen Myers (LLNL), Michael Begnaud (LANL), Brian Young (SNL) and Istvan Bondar (Research Center for Astronomy and Earth Sciences, Hungary) to serve as a presenters/trainers at the “NDC Capacity Building Workshop and Regional Seismic Travel Time (RSTT) in combination with Data Sharing and Integration Training” September 4-8 2022 in Muscat, Oman (See Appendix A for the agenda). The workshop and training (workshop from here forward) was organized by the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO) Provisional Technical Secretariat (PTS). The first half of the week was devoted to NDC workshop activities, and the second half was devoted to RSTT training. Fifty-five participants from 27 countries and the CTBTO-PTS attended the 5-day workshop (See Appendix B for list of participants and countries of origin). Presentations from the PTS described the International Monitoring System (IMS), International Data Centre (IDC) products, and metrics of regional data utilization. Contributed presentations from each country’s scientists included descriptions of regional and national networks, methods of data analysis, and needs for material and technical assistance. Training included an overview of the RSTT method and instruction on how to locate seismic events with the iLoc program, which utilizes RSTT travel times to reduce bias in event location estimates. Methods of seismic tomography and the need for a high-quality tomographic set, including seismological “ground truth”, were emphasized. Seismological “ground truth” or “GT” is a term that has come to mean both events with known location and events with well-characterized locations that are estimated using seismological data, typically with epicenter accuracy of 5 km or better. Notably, the instructional platform has migrated from UNIX shell scripts to Jupyter Notebooks. Jupyter Notebooks have the advantage being more visually intuitive, including display of graphics within the notebook. Each notebook includes every processing step that participants need to reproduce the entire exercise.

58 GEOSCIENCES↗

DMTN-107: Options for Alert Production in LSST Operations Year 1

This document reviews five options for alert production in LSST Operations Year 1 (LOY1), taking into account any implications on LSST formal requirements including up-scopes, down-scopes or explicit violations. The Data Management System Science Team's preferred option for maximizing LSST science is to generate template images from as much of the data from the commissioning and/or science verification phases as possible and use them to run Difference Image Analysis and alert production during LOY1. A proposal to increase the sky area covered during commissioning via a "filler" scheduler program is also presented. As a potential moderate up-scope, this study presents an option to build interim templates on a $\sim$monthly basis during LOY1, which could increase the accessible sky area by ~1000-2000 deg**2 per month, and should be reconsidered closer to the start of Operations.

79 ASTRONOMY AND ASTROPHYSICS↗

Milestone 7714 National Opacity Program: High-Fidelity Iron Data on NIF

Establish the level of reproducibility of NIF iron opacity shots, quantify to-date improvements in data quality through, e.g. background reduction, improved signal, and improved analysis methods, and formalize the approach that will be taken in the future to compare NIF data, Z data, and theory.

74 ATOMIC AND MOLECULAR PHYSICS↗

Prospects for future studies using deep imaging: analysis of individual Galactic cirrus filaments

The presence of Galactic cirrus is an obstacle for studying both faint objects in our Galaxy and low surface brightness extragalactic structures. With the aim of studying individual cirrus filaments in Sloan Digital Sky Survey (SDSS) Stripe 82 data, we develop techniques based on machine learning and neural networks that allow one to isolate filaments from foreground and background sources in the entirety of Stripe 82 with a precision similar to that of the human expert. Our photometric study of individual filaments indicates that only those brighter than 26 mag arcsec −2 in the SDSS r band are likely to be identified in SDSS Stripe 82 data by their distinctive colours in the optical bands. We also show a significant impact of data processing (e.g. flat-fielding, masking of bright stars, and sky subtraction) on colour estimation. Analysing the distribution of filaments’ colours with the help of mock simulations, we conclude that most filaments have colours in the following ranges: 0.55 ≤g − r ≤ 0.73 and 0.01 ≤ r − i ≤ 0.33. Our work provides a useful framework for an analysis of all types of low surface brightness features (cirri, tidal tails, stellar streams, etc.) in existing and future deep optical surveys. For practical purposes, we provide the catalogue of dust filaments.

79 ASTRONOMY AND ASTROPHYSICS↗