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

Planck 2015 results: XII. Full focal plane simulations

In this paper, we present the 8th full focal plane simulation set (FFP8), deployed in support of the Planck 2015 results. FFP8 consists of 10 fiducial mission realizations reduced to 18 144 maps, together with the most massive suite of Monte Carlo realizations of instrument noise and CMB ever generated, comprising 10 4 mission realizations reduced to about 10 6 maps. The resulting maps incorporate the dominant instrumental, scanning, and data analysis effects, and the remaining subdominant effects will be included in future updates. Finally, generated at a cost of some 25 million CPU-hours spread across multiple high-performance-computing (HPC) platforms, FFP8 is used to validate and verify analysis algorithms and their implementations, and to remove biases from and quantify uncertainties in the results of analyses of the real data.

79 ASTRONOMY AND ASTROPHYSICS↗

Planck intermediate results: XLVIII. Disentangling Galactic dust emission and cosmic infrared background anisotropies

Using the Planck 2015 data release (PR2) temperature maps, we separate Galactic thermal dust emission from cosmic infrared background (CIB) anisotropies. For this purpose, we implement a specifically tailored component-separation method, the so-called generalized needlet internal linear combination (GNILC) method, which uses spatial information (the angular powerspectra) to disentangle the Galactic dust emission and CIB anisotropies. We produce significantly improved all-sky maps of Planck thermal dust emission, with reduced CIB contamination, at 353, 545, and 857 GHz. By reducing the CIB contamination of the thermal dust maps, we provide more accurate estimates of the local dust temperature and dust spectral index over the sky with reduced dispersion, especially at high Galactic latitudes above b = ±20°. We find that the dust temperature is T = (19.4 ± 1.3) K and the dust spectral index is β = 1.6 ± 0.1 averaged over the whole sky, while T = (19.4 ± 1.5) K and β = 1.6 ± 0.2 on 21% of the sky at high latitudes. Moreover, subtracting the new CIB-removed thermal dust maps from the CMB-removed Planck maps gives access to the CIB anisotropies over 60% of the sky at Galactic latitudes |b| > 20°. Because they are a significant improvement over previous Planck products, the GNILC maps are recommended for thermal dust science. The new CIB maps can be regarded as indirect tracers of the dark matter and they are recommended for exploring cross-correlations with lensing and large-scale structure optical surveys. The reconstructed GNILC thermal dust and CIB maps are delivered as Planck products.

79 ASTRONOMY AND ASTROPHYSICS↗

SCALA: In situ calibration for integral field spectrographs

The scientific yield of current and future optical surveys is increasingly limited by systematic uncertainties in the flux calibration. This is the case for type Ia supernova (SN Ia) cosmology programs, where an improved calibration directly translates into improved cosmological constraints. Current methodology rests on models of stars. Here we aim to obtain flux calibration that is traceable to state-of-the-art detector-based calibration. Methods. We present the SNIFS Calibration Apparatus (SCALA), a color (relative) flux calibration system developed for the SuperNova integral field spectrograph (SNIFS), operating at the University of Hawaii 2.2 m (UH 88) telescope. Results. By comparing the color trend of the illumination generated by SCALA during two commissioning runs, and to previous laboratory measurements, we show that we can determine the light emitted by SCALA with a long-term repeatability better than 1%. We describe the calibration procedure necessary to control for system aging. We present measurements of the SNIFS throughput as estimated by SCALA observations. Conclusions. The SCALA calibration unit is now fully deployed at the UH 88 telescope, and with it color-calibration between 4000 Å and 9000 Å is stable at the percent level over a one-year baseline.

79 ASTRONOMY AND ASTROPHYSICS↗

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses: XVI. Time delays for the quadruply imaged quasar DES J0408–5354 with high-cadence photometric monitoring* [COSMOGRAIL XVI: Time delays for the quadruply imaged quasar DES J0408–5354 with high-cadence photometric monitoring]

Here, we present time-delay measurements for the new quadruple imaged quasar DES J0408–5354, the first quadruple imaged quasar found in the Dark Energy Survey (DES). Our result is made possible by implementing a new observational strategy using almost daily observations with the MPIA 2.2 m telescope at La Silla observatory and deep exposures reaching a signal-to-noise ratio of about 1000 per quasar image. This data qualityallows us to catch small photometric variations (a few mmag rms) of the quasar, acting on temporal scales much shorter than microlensing, and hence making the time delay measurement very robust against microlensing. In only seven months we very accurately measured one of the time delays in DES J0408–5354: Δt(AB) = –112.1 ± 2.1 days (1.8%) using only the MPIA 2.2 m data. In combination with data taken with the 1.2 m Euler Swiss telescope, we also measured two delays involving the D component of the system Δt(AD) = –155.5 ± 12.8 days (8.2%) and Δt(BD) = –42.4 ± 17.6 days (41%), where all the error bars include systematics. Turning these time delays into cosmological constraints will require deep Hubble Space Telescope (HST) imaging or ground-based adaptive optics (AO), and information on the velocity field of the lensing galaxy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Planck intermediate results: LIII. Detection of velocity dispersion from the kinetic Sunyaev-Zeldovich effect

Using the Planck full-mission data, we present a detection of the temperature (and therefore velocity) dispersion due to the kinetic Sunyaev-Zeldovich (kSZ) effect from clusters of galaxies. To suppress the primary CMB and instrumental noise we derive a matched filter and then convolve it with the Planck foreground-cleaned `2D-ILC' maps. By using the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC), we determine the normalized rms dispersion of the temperature fluctuations at the positions of clusters, finding that this shows excess variance compared with the noise expectation. We then build an unbiased statistical estimator of the signal, determining that the normalized mean temperature dispersion of 1526 clusters is ((ΔT/T) 2 ) = (1.64±0.48)×10 -11 . However, comparison with analytic calculations and simulations suggest that around 0.7σ of this result is due to cluster lensing rather than the kSZ effect. By correcting this, the temperature dispersion is measured to be ((ΔT/T) 2 ) = (1.35±0.48)×10 -11 , which gives a detection at the 2.8σ level. We further convert uniform-weight temperature dispersion into a measurement of the line-of-sight velocity dispersion, by using estimates of the optical depth of each cluster (which introduces additional uncertainty into the estimate). We find that the velocity dispersion is (v 2 )=(123000±71000)(kms -1 ) 2 , which is consistent with findings from other large-scale structure studies, and provides direct evidence of statistical homogeneity on scales of 600 h -1 Mpc. Our study shows the promise of using cross-correlations of the kSZ effect with large-scale structure in order to constrain the growth of structure.

79 ASTRONOMY AND ASTROPHYSICS↗

COSMOGRAIL: XVII. Time delays for the quadruply imaged quasar PG 1115+080

Here, we present time-delay estimates for the quadruply imaged quasar PG 1115+080. Our results are based on almost daily observations for seven months at the ESO MPIA 2.2m telescope at La Silla Observatory, reaching a signal-to-noise ratio of about 1000 per quasar image. In addition, we re-analyze existing light curves from the literature that we complete with an additional three seasons of monitoring with the Mercator telescope at La Palma Observatory. When exploring the possible source of bias we consider the so-called microlensing time delay, a potential source of systematic error so far never directly accounted for in previous time-delay publications. In fifteen years of data on PG 1115+080, we find no strong evidence of microlensing time delay. Therefore not accounting for this effect, our time-delay estimates on the individual data sets are in good agreement with each other and with the literature. Combining the data sets, we obtain the most precise time-delay estimates to date on PG 1115+080, with Δt(AB) = 8.3 +1.5 –1.6 days (18.7% precision), Δt(AC) = 9.9 +1.1 –1.1 days (11.1%) and Δt(BC) = 18.8 +1.6 –1.6 days (8.5%). Turning these time delays into cosmological constraints is done in a companion paper that makes use of ground-based Adaptive Optics (AO) with the Keck telescope.

79 ASTRONOMY AND ASTROPHYSICS↗

TDCOSMO: II. Six new time delays in lensed quasars from high-cadence monitoring at the MPIA 2.2m telescope

In this paper we present six new time-delay measurements obtained from $R_c$-band monitoring data acquired at the Max Planck Institute for Astrophysics (MPIA) 2.2 m telescope at La Silla observatory between October 2016 and February 2020. The lensed quasars HE 0047–1756, WG 0214–2105, DES 0407–5006, 2M 1134–2103, PSJ 1606–2333, and DES 2325–5229 were observed almost daily at high signal-to-noise ratio to obtain high-quality light curves where we can record fast and small-amplitude variations of the quasars. We measured time delays between all pairs of multiple images with only one or two seasons of monitoring with the exception of the time delays relative to image D of PSJ 1606–2333. The most precise estimate was obtained for the delay between image A and image B of DES 0407–5006, where $τ_{AB} = –128.4_{–3.8}^{+3.5}$ d (2.8% precision) including systematics due to extrinsic variability in the light curves. For HE 0047–1756, we combined our high-cadence data with measurements from decade-long light curves from previous COSMOGRAIL campaigns, and reach a precision of 0.9 d on the final measurement. The present work demonstrates the feasibility of measuring time delays in lensed quasars in only one or two seasons, provided high signal-to-noise ratio data are obtained at a cadence close to daily.

79 ASTRONOMY AND ASTROPHYSICS↗

Image segmentation for analyzing galaxy-galaxy strong lensing systems

The goal of this Letter is to develop a machine learning model to analyze the main gravitational lens and detect dark substructure (subhalos) within simulated images of strongly lensed galaxies. Using the technique of image segmentation, we turn the task of identifying subhalos into a classification problem, where we label each pixel in an image as coming from the main lens, a subhalo within a binned mass range, or neither. Our network is only trained on images with a single smooth lens and either zero or one subhalo near the Einstein ring. On an independent test set with lenses with large ellipticities, quadrupole and octopole moments, and for source apparent magnitudes between 17–25, the area of the main lens is recovered accurately. On average, only 1.3% of the true area is missed and 1.2% of the true area is added to another part of the lens. In addition, subhalos as light as 10 8.5 M ⊙ can be detected if they lie in bright pixels along the Einstein ring. Furthermore, the model is able to generalize to new contexts it has not been trained on, such as locating multiple subhalos with varying masses or more than one large smooth lens.

79 ASTRONOMY AND ASTROPHYSICS↗

Substructure in the stellar halo near the Sun: I. Data-driven clustering in integrals-of-motion space

Context. Merger debris is expected to populate the stellar haloes of galaxies. In the case of the Milky Way, this debris should be apparent as clumps in a space defined by the orbital integrals of motion of the stars. Aims. Our aim is to develop a data-driven and statistics-based method for finding these clumps in integrals-of-motion space for nearby halo stars and to evaluate their significance robustly. Methods. We used data from Gaia EDR3, extended with radial velocities from ground-based spectroscopic surveys, to construct a sample of halo stars within 2.5 kpc from the Sun. We applied a hierarchical clustering method that makes exhaustive use of the single linkage algorithm in three-dimensional space defined by the commonly used integrals of motion energy E, together with two components of the angular momentum, L z and L ⊥ . To evaluate the statistical significance of the clusters, we compared the density within an ellipsoidal region centred on the cluster to that of random sets with similar global dynamical properties. By selecting the signal at the location of their maximum statistical significance in the hierarchical tree, we extracted a set of significant unique clusters. By describing these clusters with ellipsoids, we estimated the proximity of a star to the cluster centre using the Mahalanobis distance. Additionally, we applied the HDBSCAN clustering algorithm in velocity space to each cluster to extract subgroups representing debris with different orbital phases. Results. Our procedure identifies 67 highly significant clusters (> 3σ), containing 12% of the sources in our halo set, and 232 subgroups or individual streams in velocity space. In total, 13.8% of the stars in our data set can be confidently associated with a significant cluster based on their Mahalanobis distance. Inspection of the hierarchical tree describing our data set reveals a complex web of relations between the significant clusters, suggesting that they can be tentatively grouped into at least six main large structures, many of which can be associated with previously identified halo substructures, and a number of independent substructures. This preliminary conclusion is further explored in a companion paper, in which we also characterise the substructures in terms of their stellar populations. Conclusions. Our method allows us to systematically detect kinematic substructures in the Galactic stellar halo with a data-driven and interpretable algorithm. The list of the clusters and the associated star catalogue are provided in two tables available at the CDS.

79 ASTRONOMY AND ASTROPHYSICS↗

Absolute colors and phase coefficients of asteroids

Context: We use the phase curves of small bodies to measure absolute magnitudes and, together with complementary theoretical and laboratory results, to understand the micro- and macroscopic properties of their surfaces. Although we can observe asteroids up to phase angles of about 30 deg, the range of phase angles covered by outer Solar System objects usually does not go further than 7-10 deg for centaurs and 2 deg for trans-Neptunian objects, and a linear relation between magnitude and phase angle may be assumed. Aims: We aim to directly compare data taken for objects in the inner Solar System (inside the orbit of Jupiter) with data for centaurs and trans-Neptunian objects. Methods: We use the SLOAN Moving Objects Catalog (MOC) data to construct phase curves restricted to phase angles of less than or equal to 7.5 deg, which is compatible with the angles observed for the trans-Neptunian and Centaur population. We assume a linear model for the photometric behavior to obtain absolute magnitudes and phase coefficients in the ugirz, V, and R filters. Results: We obtain absolute magnitudes in seven filters for over 4000 objects. Our comparison with outer Solar System objects suggests the surfaces have a common property: intrinsically redder objects become blue with increasing phase angle, while the opposite happens for intrinsically bluer objects.

79 ASTRONOMY AND ASTROPHYSICS↗

BEYONDPLANCK VI. Noise characterization and modeling

We present a Bayesian method for estimating instrumental noise parameters and propagating noise uncertainties within the global BeyondPlanck Gibbs sampling framework, and apply this to Planck Low Frequency Instrument (LFI) time-ordered data. Following previous literature, we initially adopt a 1/ f model for the noise power spectral density (PSD), but find the need for an additional lognormal component in the noise model for the 30 and 44 GHz bands. We implement an optimal Wiener-filter (or constrained realization) gap-filling procedure to account for masked data. We then use this procedure to both estimate the gapless correlated noise in the time-domain, n corr , and to sample the noise PSD parameters, ξ n = {σ 0 , f knee , α, A p }. In contrast to previous Planck analyses, we assume piecewise stationary noise only within each pointing period (PID), not throughout the full mission, but we adopt the LFI Data Processing Center (DPC) results as priors on α and f knee . On average, we find best-fit correlated noise parameters that are mostly consistent with previous results, with a few notable exceptions. However, a detailed inspection of the time-dependent results reveals many important findings. First and foremost, we find strong evidence for statistically significant temporal variations in all noise PSD parameters, many of which are directly correlated with satellite housekeeping data. Second, while the simple 1/ f model appears to be an excellent fit for the LFI 70 GHz channel, there is evidence for additional correlated noise not described by a 1/ f model in the 30 and 44 GHz channels, including within the primary science frequency range of 0.1–1 Hz. In general, most 30 and 44 GHz channels exhibit deviations from 1/ f at the 2–3σ level in each one hour pointing period, motivating the addition of the lognormal noise component for these bands. For some periods of time, we also find evidence of strong common mode noise fluctuations across the entire focal plane. Overall, we conclude that a simple 1/ f profile is not adequate to fully characterize the Planck LFI noise, even when fitted hour-by-hour, and a more general model is required. These findings have important implications for large-scale CMB polarization reconstruction with the Planck LFI data, and the current work is a first attempt at understanding and mitigating these issues.

79 ASTRONOMY AND ASTROPHYSICS↗

Discovering strongly lensed quasar candidates with catalogue-based methods from DESI Legacy Surveys

The Hubble tension, revealed by a ~5σ discrepancy between measurements of the Hubble-Lemaitre constant among observations of the early and local Universe, is one of the most significant problems in modern cosmology. In order to better understand the origin of this mismatch, independent techniques to measure H 0 , such as strong lensing time delays, are required. Notably, the sample size of such systems is key to minimising the statistical uncertainties and cosmic variance, which can be improved by exploring the datasets of large-scale sky surveys such as Dark Energy Spectroscopic Instrument (DESI). We identify possible strong lensing time-delay systems within DESI by selecting candidate multiply imaged lensed quasars from a catalogue of 24 440 816 candidate QSOs contained in the ninth data release of the DESI Legacy Imaging Surveys (DESI-LS). Using a friend-of-friends-like algorithm on spatial co-ordinates, our method generates an initial list of compact quasar groups. This list is subsequently filtered using a measure of the similarity of colours among a group’s members and the likelihood that they are quasars. A visual inspection finally selects candidate strong lensing systems based on the spatial configuration of the group members. We identified 620 new candidate multiply imaged lensed quasars (101 grade-A, 214 grade-B, 305 grade-C). This number excludes 53 known spectroscopically confirmed systems and existing candidate systems identified in other similar catalogues. When available, these new candidates will be further checked by combining the spectroscopic and photometric data from DESI.

79 ASTRONOMY AND ASTROPHYSICS↗

Bicep2. III. INSTRUMENTAL SYSTEMATICS

In a companion paper, we have reported a >5σ detection of degree scale B-mode polarization at 150 GHz by the Bicep2 experiment. Here we provide a detailed study of potential instrumental systematic contamination to that measurement. We focus extensively on spurious polarization that can potentially arise from beam imperfections. We present a heuristic classification of beam imperfections according to their symmetries and uniformities, and discuss how resulting contamination adds or cancels in maps that combine observations made at multiple orientations of the telescope about its boresight axis. We introduce a technique, which we call “deprojection,” for filtering the leading order beam-induced contamination from time-ordered data, and show that it reduces power in Bicep2's actual and null-test BB spectra consistent with predictions using high signal-to-noise beam shape measurements. We detail the simulation pipeline that we use to directly simulate instrumental systematics and the calibration data used as input to that pipeline. Finally, we present the constraints on BB contamination from individual sources of potential systematics. We find that systematics contribute BB power that is a factor of ∼10× below Bicep2's three-year statistical uncertainty, and negligible compared to the observed BB signal. The contribution to the best-fit tensor/scalar ratio is at a level equivalent to r = (3–6) × 10{sup −3}.

79 ASTRONOMY AND ASTROPHYSICS↗

Instrumental Response Model and Detrending for the Dark Energy Camera

We describe the model for mapping from sky brightness to the digital output of the Dark Energy Camera (DECam) and the algorithms adopted by the Dark Energy Survey (DES) for inverting this model to obtain photometric measures of celestial objects from the raw camera output. This calibration aims for fluxes that are uniform across the camera field of view and across the full angular and temporal span of the DES observations, approaching the accuracy limits set by shot noise for the full dynamic range of DES observations. The DES pipeline incorporates several substantive advances over standard detrending techniques, including principal-components-based sky and fringe subtraction; correction of the "brighter-fatter" nonlinearity; use of internal consistency in on-sky observations to disentangle the influences of quantum efficiency, pixel-size variations, and scattered light in the dome flats; and pixel-by-pixel characterization of instrument spectral response, through combination of internal-consistency constraints with auxiliary calibration data. This article provides conceptual derivations of the detrending/calibration steps, and the procedures for obtaining the necessary calibration data. Other publications will describe the implementation of these concepts for the DES operational pipeline, the detailed methods, and the validation that the techniques can bring DECam photometry and astrometry within $\approx 2$ mmag and $\approx 3$ mas, respectively, of fundamental atmospheric and statistical limits. In conclusion, the DES techniques should be broadly applicable to wide-field imagers.

79 ASTRONOMY AND ASTROPHYSICS↗

Photometric Characterization of the Dark Energy Camera

We characterize the variation in photometric response of the Dark Energy Camera (DECam) across its 520 Mpix science array during 4 years of operation. These variations are measured using high signal-to-noise aperture photometry of >10 7 stellar images in thousands of exposures of a few selected fields, with the telescope dithered to move the sources around the array. A calibration procedure based on these results brings the rms variation in aperture magnitudes of bright stars on cloudless nights down to 2–3 mmag, with <1 mmag of correlated photometric errors for stars separated by ≥20''. On cloudless nights, any departures of the exposure zeropoints from a secant airmass law exceeding 1 mmag are plausibly attributable to spatial/temporal variations in aperture corrections. These variations can be inferred and corrected by measuring the fraction of stellar light in an annulus between 6'' and 8'' diameter. Key elements of this calibration include: correction of amplifier nonlinearities; distinguishing pixel-area variations and stray light from quantum-efficiency variations in the flat fields; field-dependent color corrections; and the use of an aperture-correction proxy. The DECam response pattern across the 2° field drifts over months by up to ±9 mmag, in a nearly wavelength-independent low-order pattern. Here, we find no fundamental barriers to pushing global photometric calibrations toward mmag accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗

Stellar loci IV. red giant stars

In the fourth paper of this series, we present the metallicity-dependent Sloan Digital Sky Survey (SDSS) stellar color loci of red giant stars, using a spectroscopic sample of red giants in the SDSS Stripe 82 region. The stars span a range of 0.55 – 1.2 mag in color g – i, –0.3 – –2.5 in metallicity [Fe/H], and have values of surface gravity log g smaller than 3.5 dex. As in the case of main-sequence (MS) stars, the intrinsic widths of loci of red giants are also found to be quite narrow, a few mmag at maximum. There are however systematic differences between the metallicity-dependent stellar loci of red giants and MS stars. The colors of red giants are less sensitive to metallicity than those of MS stars. With good photometry, photometric metallicities of red giants can be reliably determined by fitting the u – g, g – r, r – i, and i – z colors simultaneously to an accuracy of 0.2 – 0.25 dex, comparable to the precision achievable with low-resolution spectroscopy for a signal-to-noise ratio of 10. By comparing fitting results to the stellar loci of red giants and MS stars, we propose a new technique to discriminate between red giants and MS stars based on the SDSS photometry. The technique achieves completeness of ~70 per cent and efficiency of ~80 per cent in selecting metal-poor red giant stars of [Fe/H] ≤ –1.2. It thus provides an important tool to probe the structure and assemblage history of the Galactic halo using red giant stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Identifications of RR Lyrae Stars and Quasars from the Simulated Data of Mephisto-W Survey

We have investigated the feasibilities and accuracies of the identifications of RR Lyrae stars and quasars from the simulated data of the Multi-channel Photometric Survey Telescope (Mephisto) W Survey. Based on the variable sources light curve libraries from the Sloan Digital Sky Survey (SDSS) Stripe 82 data and the observation history simulation from the Mephisto-W Survey Scheduler, we have simulated the uvgriz multi-band light curves of RR Lyrae stars, quasars and other variable sources for the first-year observation of Mephisto W Survey. We have applied the ensemble machine learning algorithm Random Forest Classifier (RFC) to identify RR Lyrae stars and quasars, respectively. We build training and test samples and extract ~150 features from the simulated light curves and train two RFCs respectively for the RR Lyrae star and quasar classification. We find that, our RFCs are able to select the RR Lyrae stars and quasars with remarkably high precision and completeness, with purity = 95.4% and completeness = 96.9% for the RR Lyrae RFC and purity = 91.4% and completeness = 90.2% for the quasar RFC. In conclusion, we have also derived relative importances of the extracted features utilized to classify RR Lyrae stars and quasars.

(galaxies:) quasars: general↗

Elliptic multipoles and the modeling of narrow-gap bend magnets in accelerators

We highlight the virtues of 2D elliptic-multipole field expansions in modeling the magnetic fields of narrow-aperture, straight-axis bending magnets with parallel faces, addressing the limitations of the conventional circular multipole series when the beam-orbit sagitta exceeds the magnet's vertical half-gap. The elliptic multipoles provide a convenient way to represent the field in all aspects of the magnet development (design, particle-tracking simulations, measurements). We propose a numerically robust method of data analysis to determine the elliptic (or circular) multipoles from stretched-wire measurements with the wire moving on an arbitrary path.

Venturini, Marco↗