Search NASA⌕ Search

SEARCH · Search NASA

Results for “Methods: 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 289 records · Page 16

Discovering Ca II absorption lines with a neural network

Quasar absorption line analysis is critical for studying gas and dust components and their physical and chemical properties as well as the evolution and formation of galaxies in the early universe. Calcium II (Ca II ) absorbers, which are one of the dustiest absorbers and are located at lower redshifts than most other absorbers, are especially valuable when studying physical processes and conditions in recent galaxies. However, the number of known quasar Ca II absorbers is relatively low due to the difficulty of detecting them with traditional methods. In this work, we developed an accurate and quick approach to search for Ca II absorption lines using deep learning. In our deep learning model, a convolutional neural network, tuned using simulated data, is used for the classification task. The simulated training data are generated by inserting artificial Ca II absorption lines into original quasar spectra from the Sloan Digital Sky Survey (SDSS), while an existing Ca II catalogue is adopted as the test set. The resulting model achieves an accuracy of 96 per cent on the real data in the test set. Our solution runs thousands of times faster than traditional methods, taking a fraction of a second to analyse thousands of quasars, while traditional methods may take days to weeks. The trained neural network is applied to quasar spectra from SDSS’s DR7 and DR12 and discovered 399 new quasar Ca II absorbers. In addition, we confirmed 409 known quasar Ca II absorbers identified previously by other research groups through traditional methods.

79 ASTRONOMY AND ASTROPHYSICS↗

J-PLUS: Support vector machine applied to STAR-GALAXY-QSO classification

Context. In modern astronomy, machine learning has proved to be efficient and effective in mining big data from the newest telescopes. Aims. In this study, we construct a supervised machine-learning algorithm to classify the objects in the Javalambre Photometric Local Universe Survey first data release (J-PLUS DR1). Methods. The sample set is featured with 12-waveband photometry and labeled with spectrum-based catalogs, including Sloan Digital Sky Survey spectroscopic data, the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, and VERONCAT – the Veron Catalog of Quasars & AGN. The performance of the classifier is presented with the applications of blind test validations based on RAdial Velocity Extension, the Kepler Input Catalog, the Two Micron All Sky Survey Redshift Survey, and the UV-bright Quasar Survey. A new algorithm was applied to constrain the potential extrapolation that could decrease the performance of the machine-learning classifier. Results. The accuracies of the classifier are 96.5% in the blind test and 97.0% in training cross-validation. The F1-scores for each class are presented to show the balance between the precision and the recall of the classifier. We also discuss different methods to constrain the potential extrapolation.

79 ASTRONOMY AND ASTROPHYSICS↗

Exploring Telematics Big Data for Truck Platooning Opportunities

NREL completed a temporal and geospatial analysis of telematics data to estimate the fraction of platoonable miles traveled by class 8 tractor trailers currently in operation. This paper discusses the value and limitations of very large but low time-resolution data sets, and the fuel consumption reduction opportunities from large scale adoption of platooning technology for class 8 highway vehicles in the US based on telematics data. The telematics data set consist of about 57,000 unique vehicles traveling over 210 million miles combined during a two-week period. 75% of the total fuel consumption result from vehicles operating in top gear, suggesting heavy highway utilization. The data is at a one-hour resolution, resulting in a significant fraction of data be uncategorizable, yet significant value can still be extracted from the remaining data. Multiple analysis methods to estimate platoonable miles are discussed. Results indicate that 63% of total miles driven at known hourly-average speeds happens at speeds amenable to platooning. When also considering availability of nearby partner vehicles, results indicate 55.7% of all classifiable miles driven were platoonable. Analysis also address the availability of numerous partners enabling platoons greater than 2 trucks and the percentage of trucks that would be required to be equipped with platooning equipment to realize more than 50% of the possible savings.

33 ADVANCED PROPULSION SYSTEMS↗

Impact of point spread function higher moments error on weak gravitational lensing

ABSTRACT Weak gravitational lensing is one of the most powerful tools for cosmology, while subject to challenges in quantifying subtle systematic biases. The point spread function (PSF) can cause biases in weak lensing shear inference when the PSF model does not match the true PSF that is convolved with the galaxy light profile. Although the effect of PSF size and shape errors – i.e. errors in second moments – is well studied, weak lensing systematics associated with errors in higher moments of the PSF model require further investigation. The goal of our study is to estimate their potential impact for LSST weak lensing analysis. We go beyond second moments of the PSF by using image simulations to relate multiplicative bias in shear to errors in the higher moments of the PSF model. We find that the current level of errors in higher moments of the PSF model in data from the Hyper Suprime-Cam survey can induce a ∼0.05 per cent shear bias, making this effect unimportant for ongoing surveys but relevant at the precision of upcoming surveys such as LSST.

79 ASTRONOMY AND ASTROPHYSICS↗

GAMA/DEVILS: constraining the cosmic star formation history from improved measurements of the 0.3–2.2 μm extragalactic background light

We present a revised measurement of the optical extragalactic background light (EBL), based on the contribution of resolved galaxies to the integrated galaxy light (IGL). The cosmic optical background radiation (COB), encodes the light generated by star formation, and provides a wealth of information about the cosmic star formation history (CSFH). We combine wide and deep galaxy number counts from the Galaxy And Mass Assembly survey (GAMA) and Deep Extragalactic VIsible Legacy Survey (DEVILS), along with the Hubble Space Telescope (HST) archive and other deep survey data sets, in nine multiwavelength filters to measure the COB in the range from 0.35 μm to 2.2 μm. We derive the luminosity density in each band independently and show good agreement with recent and complementary estimates of the optical-EBL from very high-energy (VHE) experiments. Our error analysis suggests that the IGL and γ-ray measurements are now fully consistent to within $\sim 10{{\ \rm per\ cent}}$, suggesting little need for any additional source of diffuse light beyond the known galaxy population. We use our revised IGL measurements to constrain the CSFH, and place amplitude constraints on a number of recent estimates. As a consistency check, we can now demonstrate convincingly, that the CSFH, stellar mass growth, and the optical-EBL provide a fully consistent picture of galaxy evolution. Here, we conclude that the peak of star formation rate lies in the range 0.066–0.076 M ⊙ yr -1 Mpc -3 at a lookback time of 9.1 to 10.9 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗

Reconstructing probability distributions with Gaussian processes

Modern cosmological analyses constrain physical parameters using Markov Chain Monte Carlo (MCMC) or similar sampling techniques. Oftentimes, these techniques are computationally expensive to run and require up to thousands of CPU hours to complete. Here we present a method for reconstructing the log-probability distributions of completed experiments from an existing chain (or any set of posterior samples). Here, the reconstruction is performed using Gaussian process regression for interpolating the log-probability. This allows for easy resampling, importance sampling, marginalization, testing different samplers, investigating chain convergence, and other operations. As an example use case, we reconstruct the posterior distribution of the most recent Planck 2018 analysis. We then resample the posterior, and generate a new chain with 40 times as many points in only 30 min. Our likelihood reconstruction tool is made publicly available online.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A general framework for removing point-spread function additive systematics in cosmological weak lensing analysis

ABSTRACT Cosmological weak lensing measurements rely on a precise measurement of the shear two-point correlation function (2PCF) along with a deep understanding of systematics that affect it. In this work, we demonstrate a general framework for detecting and modelling the impact of PSF systematics on the cosmic shear 2PCF and mitigating its impact on cosmological analysis. Our framework can detect PSF leakage and modelling error from all spin-2 quantities contributed by the PSF second and higher moments, rather than just the second moments, using the cross-correlations between galaxy shapes and PSF moments. We interpret null tests using the HSC Year 3 (Y3) catalogs with this formalism and find that leakage from the spin-2 combination of PSF fourth moments is the leading contributor to additive shear systematics, with total contamination that is an order-of-magnitude higher than that contributed by PSF second moments alone. We conducted a mock cosmic shear analysis for HSC Y3 and find that, if uncorrected, PSF systematics can bias the cosmological parameters Ωm and S8 by ∼0.3σ. The traditional second moment-based model can only correct for a 0.1σ bias, leaving the contamination largely uncorrected. We conclude it is necessary to model both PSF second and fourth moment contaminations for HSC Y3 cosmic shear analysis. We also reanalyse the HSC Y1 cosmic shear analysis with our updated systematics model and identify a 0.07σ bias on Ωm when using the more restricted second moment model from the original analysis. We demonstrate how to self-consistently use the method in both real space and Fourier space, assess shear systematics in tomographic bins, and test for PSF model overfitting.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Median Statistics Analysis of Deuterium Abundance Measurements and Spatial Curvature Constraints

Zavarygin et al. compiled a list of 15 deuterium abundance measurements; two were discarded because the remaining 13 measurements are then consistent with Gaussianity. Here, those authors found that the weighted-mean baryon density (Ω b h 2) determined from the 13 measurements is mildly discrepant (1.6σ) with that determined from the Planck 2015 cosmic microwave background anisotropy data in a flat cosmogony. We find that a median statistic central estimate of Ω b h 2 from all 15 deuterium abundance measurements is a more accurate estimate; is very consistent with Ω b h 2 estimated from Planck 2015 data in a flat cosmogony, but is about 2σ lower than that found in a closed cosmogonical model from the Planck 2015 data.

79 ASTRONOMY AND ASTROPHYSICS↗

Planck constraints on the tensor-to-scalar ratio

In this work, we present constraints on the tensor-to-scalar ratio r using Planck data. We use the latest release of Planck maps (PR4), processed with the NPIPE code, which produces calibrated frequency maps in temperature and polarization for all Planck channels from 30 GHz to 857 GHz using the same pipeline. We computed constraints on r using the BB angular power spectrum, and we also discuss constraints coming from the TT spectrum. Given Planck's noise level, the TT spectrum gives constraints on r that are cosmic-variance limited (with $\sigma$ r = 0.093), but we show that the marginalized posterior peaks towards negative values of r at about the 1.2$\sigma$ level. We derived Planck constraints using the BB power spectrum at both large angular scales (the 'reionization bump') and intermediate angular scales (the 'recombination bump') from $\ell$=2 to 150, and find a stronger constraint than that from TT, with $\sigma$ r = 0.069. The Planck BB spectrum shows no systematic bias, and is compatible with zero, given both the statistical noise and the systematic uncertainties. The likelihood analysis using B modes yields the constraint r < 0.158 at 95% confidence using more than 50% of the sky. This upper limit tightens to r < 0.069 when Planck EE, BB, and EB power spectra are combined consistently, and it tightens further to r < 0.056 when the Planck TT power spectrum is included in the combination. Finally, combining Planck with BICEP2/Keck 2015 data yields an upper limit of r < 0.044.

79 ASTRONOMY AND ASTROPHYSICS↗

SALT2 versus SALT3: updated model surfaces and their impacts on type Ia supernova cosmology

ABSTRACT For the past decade, SALT2 has been the most common model used to fit Type Ia supernova (SN Ia) light curves for dark energy analyses. Recently, the SALT3 model was released, which upgraded a number of model features but has not yet been used for measurements of dark energy. Here, we evaluate the impact of switching from SALT2 to SALT3 for a SN cosmology analysis. We train SALT2 and SALT3 on an identical training sample of 1083 well-calibrated Type Ia supernovae, ensuring that any differences found come from the underlying model framework. We publicly release the results of this training (the SALT ‘surfaces’). We then run a cosmology analysis on the public Dark Energy Survey 3-Yr Supernova data sample (DES-SN3YR), and on realistic simulations of those data. We provide the first estimate of the SN + CMB systematic uncertainty arising from the choice of SALT model framework (i.e. SALT2 versus SALT3), Δw = + 0.001 ± 0.005 – a negligible effect at the current level of dark energy analyses. We also find that the updated surfaces are less sensitive to photometric calibration uncertainties than previous SALT2 surfaces, with the average spectral energy density dispersion reduced by a factor of two over optical wavelengths. This offers an opportunity to reduce the contribution of calibration errors to SN cosmology uncertainty budgets.

79 ASTRONOMY AND ASTROPHYSICS↗

21 cm Power Spectrum Analysis of North Celestial Pole Observations with the Tianlai Dish Pathfinder Array

The Tianlai Dish Pathfinder Array (TDPA) is a radio interferometer designed to test techniques for 21 cm intensity mapping in the post-reionization Universe as a means of measuring large-scale cosmic structure. Using nine nights of observations targeting the North Celestial Pole field, totaling approximately 107 hr of integration time, we analyze data in the frequency range 700–800 MHz (corresponding to redshift z ∼ 0.9). We do the data format conversion, radio frequency interference flagging, calibration, imaging and point source subtraction, and foreground removal via Singular Value Decomposition. The spherically averaged power spectrum Δ 2 (k) is obtained. Furthermore, this work successfully establishes and validates a comprehensive data analysis framework for the TDPA. We identify key improvements including sky model refinement, increased integration time, and pipeline optimization that will enable future detection of the 21 cm signal through auto-correlation and cross-correlation with optical galaxy surveys.

cosmology: large-scale structure of universe↗

Constraining Temporal Oscillations of Cosmological Parameters Using SNe Ia

The existing set of SNe Ia is now sufficient to detect oscillatory deviations from the canonical ΛCDM cosmology. In this work we determine that the Fourier spectrum of the Pantheon data set of spectroscopically well-observed SNe Ia is consistent with the predictions of ΛCDM. We also develop and describe two complementary techniques for using SNe Ia to constrain those alternate cosmological models that predict deviations from ΛCDM that are oscillatory in conformal time. The first technique uses the reduced $χ^2$ statistic to determine the likelihood that the observed data would result from a given model. The second technique uses bootstrap analysis to determine the likelihood that the Fourier spectrum of a proposed model could result from statistical fluctuations around ΛCDM. We constrain three oscillatory alternate cosmological models: one in which the dark energy equation of state parameter oscillates around the canonical value of $w_Λ = -1$, one in which the energy density of dark energy oscillates around its ΛCDM value, and one in which gravity derives from a scalar field evolving under an oscillatory coupling. We further determine that any alternate cosmological model that produces distance modulus residuals with a Fourier amplitude of $\simeq$36 millimags is strongly ruled out, given the existing data, for frequencies between $\simeq$0.08 Gyr$^{-1} h_{100}$ and $\simeq$80 Gyr$^{-1} h_{100}$.

79 ASTRONOMY AND ASTROPHYSICS↗

A computationally efficient approach for calculating galaxytwo-point correlations

We developed a modification to the calculation of the two-point correlation function commonly used in the analysis of large-scale structure in cosmology. An estimator of the two-point correlation function is constructed by contrasting the observed distribution of galaxies with that of a uniformly populated random catalogue. The assumption that the distribution of random galaxies in redshift is independent of angular position allows us to replace pair-wise combinatorics with fast integration over probability maps. The new method significantly reduces the computation time while simultaneously increasing the precision of the calculation. Also, it only requires cosmological parameters at the last and least computationally expensive stage, which allows for faster exploration of the parameter space.

79 ASTRONOMY AND ASTROPHYSICS↗

The Effective Number Density of Galaxies for Weak Lensing Measurements in the LSST Project

Future weak lensing surveys potentially hold the highest statistical power for constrain ing cosmological parameters compared to other cosmological probes. The statistical power of a weak lensing survey is determined by the sky coverage, the inverse of the noise in shear measurements, and the galaxy number density. The combination of the latter two factors is often expressed in terms of n eff – the “effective number density of galaxies used for weak lensing measurements”. In this work, we estimate n eff for the Large Synoptic Survey Telescope (LSST) project, the most powerful ground-based lensing survey planned for the next two decades. We investigate how the following factors affect the resulting n eff of the survey with detailed simulations: (1) survey time, (2) shear measurement algorithm, (3) algorithm for combining multiple exposures, (4) inclusion of data from multiple filter bands, (5) redshift distribution of the galaxies, and (6) masking and blending. For the first time, we quantify in a general weak lensing analysis pipeline the sensitivity of n eff to the above factors. We find that with current weak lensing algorithms, expected distributions of observing parameters, and all lensing data (r- and i-band, covering 18,000 degree 2 of sky) for LSST, n eff ≈ 37 arcmin -2 before considering blending and masking, n eff ≈ 31 arcmin -2 when rejecting seriously blended galaxies and n eff ≈ 26 arcmin -2 when considering an additional 15% loss of galaxies due to masking. With future improvements in weak lensing algorithms, these values could be expected to increase by up to 20%. Throughout the paper, we also stress the ways in which n eff depends on our ability to understand and control systematic effects in the measurements.

Astrophysics↗

Realistic galaxy image simulation via score-based generative models

ABSTRACT We show that a denoising diffusion probabilistic model (DDPM), a class of score-based generative model, can be used to produce realistic mock images that mimic observations of galaxies. Our method is tested with Dark Energy Spectroscopic Instrument (DESI) grz imaging of galaxies from the Photometry and Rotation curve OBservations from Extragalactic Surveys (PROBES) sample and galaxies selected from the Sloan Digital Sky Survey. Subjectively, the generated galaxies are highly realistic when compared with samples from the real data set. We quantify the similarity by borrowing from the deep generative learning literature, using the ‘Fréchet inception distance’ to test for subjective and morphological similarity. We also introduce the ‘synthetic galaxy distance’ metric to compare the emergent physical properties (such as total magnitude, colour, and half-light radius) of a ground truth parent and synthesized child data set. We argue that the DDPM approach produces sharper and more realistic images than other generative methods such as adversarial networks (with the downside of more costly inference), and could be used to produce large samples of synthetic observations tailored to a specific imaging survey. We demonstrate two potential uses of the DDPM: (1) accurate inpainting of occluded data, such as satellite trails, and (2) domain transfer, where new input images can be processed to mimic the properties of the DDPM training set. Here we ‘DESI-fy’ cartoon images as a proof of concept for domain transfer. Finally, we suggest potential applications for score-based approaches that could motivate further research on this topic within the astronomical community.

79 ASTRONOMY AND ASTROPHYSICS↗

Planck intermediate results

We present the NPIPE processing pipeline, which produces calibrated frequency maps in temperature and polarization from data from the Planck Low Frequency Instrument (LFI) and High Frequency Instrument (HFI) using high-performance computers. NPIPE represents a natural evolution of previous Planck analysis efforts, and combines some of the most powerful features of the separate LFI and HFI analysis pipelines. For example, following the LFI 2018 processing procedure, NPIPE uses foreground polarization priors during the calibration stage in order to break scanning-induced degeneracies. Similarly, NPIPE employs the HFI 2018 time-domain processing methodology to correct for bandpass mismatch at all frequencies. In addition, NPIPE introduces several improvements, including, but not limited to: inclusion of the 8% of data collected during repointing manoeuvres; smoothing of the LFI reference load data streams; in-flight estimation of detector polarization parameters; and construction of maximally independent detector-set split maps. For component-separation purposes, important improvements include: maps that retain the CMB Solar dipole, allowing for high-precision relative calibration in higher-level analyses; well-defined single-detector maps, allowing for robust CO extraction; and HFI temperature maps between 217 and 857 GHz that are binned into 0'.9 pixels (Nside = 4096), ensuring that the full angular information in the data is represented in the maps even at the highest Planck resolutions. The net effect of these improvements is lower levels of noise and systematics in both frequency and component maps at essentially all angular scales, as well as notably improved internal consistency between the various frequency channels. Based on the NPIPE maps, we present the first estimate of the Solar dipole determined through component separation across all nine Planck frequencies. The amplitude is (3366.6 ± 2.7) μK, consistent with, albeit slightly higher than, earlier estimates. From the large-scale polarization data, we derive an updated estimate of the optical depth of reionization of τ = 0.051 ± 0.006, which appears robust with respect to data and sky cuts. There are 600 complete signal, noise and systematics simulations of the full-frequency and detector-set maps. As a Planck first, these simulations include full time-domain processing of the beam-convolved CMB anisotropies. The release of NPIPE maps and simulations is accompanied with a complete suite of raw and processed time-ordered data and the software, scripts, auxiliary data, and parameter files needed to improve further on the analysis and to run matching simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Galaxy Zoo DECaLS: Detailed visual morphology measurements from volunteers and deep learning for 314,000 galaxies

We present Galaxy Zoo DECaLS: detailed visual morphological classifications for Dark Energy Camera Legacy Survey images of galaxies within the SDSS DR8 footprint. Deeper DECaLS images (r = 23.6 versus r = 22.2 from SDSS) reveal spiral arms, weak bars, and tidal features not previously visible in SDSS imaging. To best exploit the greater depth of DECaLS images, volunteers select from a new set of answers designed to improve our sensitivity to mergers and bars. Galaxy Zoo volunteers provide 7.5 million individual classifications over 314 000 galaxies. 140 000 galaxies receive at least 30 classifications, sufficient to accurately measure detailed morphology like bars, and the remainder receive approximately 5. All classifications are used to train an ensemble of Bayesian convolutional neural networks (a state-of-the-art deep learning method) to predict posteriors for the detailed morphology of all 314 000 galaxies. We use active learning to focus our volunteer effort on the galaxies which, if labelled, would be most informative for training our ensemble. When measured against confident volunteer classifications, the trained networks are approximately 99 per cent accurate on every question. Morphology is a fundamental feature of every galaxy; our human and machine classifications are an accurate and detailed resource for understanding how galaxies evolve.

79 ASTRONOMY AND ASTROPHYSICS↗

The first full orbit of η Carinae seen by Fermi

The binary system η Carinae has completed its first 5.54 y orbit since the beginning of science operation of the Fermi Large Area Telescope (LAT). We are now able to investigate the high-energy γ-ray source at the position of η Carinae over its full orbital period. By this, we can address and confirm earlier predictions for temporal and spectral variability. Here, newer versions of the LAT datasets, instrument response functions and background models allow for a more accurate analysis. Therefore it is important to re-evaluate the previously analyzed time period along with the new data to further constrain location, spectral shape, and flux time history of the γ-ray source. As a result, we confirm earlier predictions of increasing flux values above 10 GeV toward the next periastron passage. For the most recent part of the data sample, flux values as high as those before the first periastron passage in 2008 are recorded. A comparison of spectral energy distributions around periastron and apastron passages reveals strong variation in the high-energy band. This is due to a second spectral component that is present only around periastron. In conclusion, improved spatial consistency with the γ-ray source at the position of η Carinae along with the confirmation of temporal variability above 10 GeV in conjunction with the orbital period strengthens the argument for unambiguous source identification. Spectral variability provides additional constraints for future modeling of the particle acceleration and γ-ray emission in colliding-wind binary systems.

79 ASTRONOMY AND ASTROPHYSICS↗