Search NASA⌕ Search

SEARCH · Search NASA

Results for “time domain astronomy”

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 127 records · Page 7

The large-scale environment of thermonuclear and core-collapse supernovae

ABSTRACT The new generation of wide-field time-domain surveys has made it feasible to study the clustering of supernova (SN) host galaxies in the large-scale structure (LSS) for the first time. We investigate the LSS environment of SN populations, using 106 dark matter density realisations with a resolution of ∼3.8 Mpc, constrained by the 2M+ + galaxy survey. We limit our analysis to redshift z < 0.036, using samples of 498 thermonuclear and 782 core-collapse SNe from the Zwicky Transient Facility’s Bright Transient Survey and Census of the Local Universe catalogues. We detect clustering of SNe with high significance; the observed clustering of the two SNe populations is consistent with each other. Further, the clustering of SN hosts is consistent with that of the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey DR12 spectroscopic galaxy sample in the same redshift range. Using a tidal shear classifier, we classify the LSS into voids, sheets, filaments, and knots. We find that both SNe and SDSS galaxies are predominantly found in sheets and filaments. SNe are significantly under-represented in voids and over-represented in knots compared to the volume fraction in these structures. This work opens the potential for using forthcoming wide-field deep SN surveys as a complementary LSS probe.

Astronomy & Astrophysics↗

The Time Domain Spectroscopic Survey: Changing-look Quasar Candidates from Multi-epoch Spectroscopy in SDSS-IV

Active galactic nuclei (AGN) can vary significantly in their rest-frame optical/UV continuum emission, and with strong associated changes in broad line emission, on much shorter timescales than predicted by standard models of accretion disks around supermassive black holes. Most such changing-look or changing-state AGN—and at higher luminosities, changing-look quasars (CLQs)—have been found via spectroscopic follow-up of known quasars showing strong photometric variability. The Time Domain Spectroscopic Survey of the Sloan Digital Sky Survey IV (SDSS-IV) includes repeat spectroscopy of large numbers of previously known quasars, many selected irrespective of photometric variability, and with spectral epochs separated by months to decades. Our visual examination of these repeat spectra for strong broad line variability yielded 61 newly discovered CLQ candidates. We quantitatively compare spectral epochs to measure changes in continuum and Hβ broad line emission, finding 19 CLQs, of which 15 are newly recognized. The parent sample includes only broad line quasars, so our study tends to find objects that have dimmed, i.e., turn-off CLQs. However, we nevertheless find four turn-on CLQs that meet our criteria, albeit with broad lines in both dim and bright states. We study the response of Hβ and Mg ii emission lines to continuum changes. The Eddington ratios of CLQs are low, and/or their Hβ broad line width is large relative to the overall quasar population. Repeat quasar spectroscopy in the upcoming SDSS-V black hole Mapper program will reveal significant numbers of CLQs, enhancing our understanding of the frequency and duty cycle of such strong variability, and the physics and dynamics of the phenomenon.

79 ASTRONOMY AND ASTROPHYSICS↗

Extreme Variability Quasars in Their Various States. I. The Sample Selection and Composite SDSS Spectra

Extremely variable quasars (EVQs) are a population of sources showing large optical photometric variability revealed by time-domain surveys. The physical origin of such extreme variability is yet unclear. In this first paper of a series, we construct the largest-ever sample of 14,012 EVQs using more than 15 yr of photometric data from Sloan Digital Sky Survey (SDSS) and Pan-STARRS1. We divide the EVQs into five subsamples according to the relative brightness of each EVQ during SDSS spectroscopic observation compared with the mean brightness from photometric observations. Corresponding control samples of normal quasars are built with matched redshift, bolometric luminosity, and supermassive black hole mass. We obtain the composite SDSS spectra of EVQs in various states and their corresponding control samples. We find EVQs exhibit clearly bluer SDSS spectra during bright states and clearly redder spectra during dim states, consistent with the "bluer-when-brighter" trend widely seen in normal quasars. We further find that the line equivalent widths (EWs) of broad Mg II, C IV and [O III] (but not broad Hβ, which is yet puzzling) gradually decreases from the dim state to the bright state, similar to the so-called intrinsic Baldwin effect commonly seen in normal active galactic nuclei. In addition, EVQs have consistently larger line EWs compared with the control samples. We also see that EVQs show slight excess in the very broad line component compared with control samples. Possible explanations for the discoveries are discussed. Our findings support the hypothesis that EVQs are in the tail of a broad distribution of quasar properties but are not a distinct population.

79 ASTRONOMY AND ASTROPHYSICS↗

Octofitter: Fast, Flexible, and Accurate Orbit Modeling to Detect Exoplanets

As next-generation imaging instruments and interferometers search for planets closer to their stars, they must contend with increasing orbital motion and longer integration times. These compounding effects make it difficult to detect faint planets but also present an opportunity. Increased orbital motion makes it possible to move the search for planets into the orbital domain, where direct images can be freely combined with the radial velocity and proper motion anomaly, even without a confirmed detection in any single epoch. In this paper, we present a fast and differentiable multimethod orbit-modeling and planet detection code called Octofitter. This code is designed to be highly modular and allows users to easily adjust priors, change parameterizations, and specify arbitrary function relations between the parameters of one or more planets. Octofitter further supplies tools for examining model outputs including prior and posterior predictive checks and simulation-based calibration. We demonstrate the capabilities of Octofitter on real and simulated data from different instruments and methods, including HD 91312, simulated JWST/NIRISS aperture masking interferometry observations, radial velocity curves, and grids of images from the Gemini Planet Imager. We show that Octofitter can reliably recover faint planets in long sequences of images with arbitrary orbital motion. This publicly available tool will enable the broad application of multiepoch and multimethod exoplanet detection, which could improve how future targeted ground- and space-based surveys are performed. Finally, its rapid convergence makes it a useful addition to the existing ecosystem of tools for modeling the orbits of directly imaged planets.

79 ASTRONOMY AND ASTROPHYSICS↗

Universal CMB 𝐵-mode spectrum from early causal tensor sources

Many early Universe scenarios predict postinflationary tensor perturbations from causality-limited, subhorizon sources. While the microphysical details may be different, as long as these sources are bounded in duration and correlation length, their tensor power spectra exhibit a universal scaling behavior at small wave number: 𝒫 ℎ ⁡(𝑘)∝𝑘 3 , corresponding to white noise on superhorizon scales at the time of production. If these early causal tensor sources (ECTs) exclusively produce gravitational waves before redshift 𝑧 ∼10 5 , this scaling is realized on all of the scales observed in the CMB, and thus yields a universal multipole distribution for the 𝐵-mode angular power spectrum. Unlike the scale-invariant distributions of inflationary 𝐵 modes, ECTs generically predict enhanced power on small scales and suppressed power on large scales, which allows these source classes to be distinguished given measurements over a sufficient range of angular scales. In this paper, we introduce a unified framework for characterizing ECTs and demonstrate how their universal infrared scaling manifests in low-frequency observables, including CMB 𝐵 modes and stochastic gravitational wave spectral densities. We illustrate this mapping with representative case studies of this universality class involving first-order phase transitions, topological defects, and enhanced scalar perturbations, which source tensor modes at second order in perturbation theory.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraining Type Ia supernova explosions and early flux excesses with the Zwicky Transient Factory

In the new era of time-domain surveys, Type Ia supernovae are being caught sooner after explosion, which has exposed significant variation in their early light curves. Two driving factors for early-time evolution are the distribution of 56 Ni in the ejecta and the presence of flux excesses of various causes. We perform an analysis of the largest young SN Ia sample to date. We compare 115 SN Ia light curves from the Zwicky Transient Facility to the turtls model grid containing light curves of Chandrasekhar mass explosions with a range of 56 Ni masses, 56 Ni distributions, and explosion energies. We find that the majority of our observed light curves are well reproduced by Chandrasekhar mass explosion models with a preference for highly extended 56 Ni distributions. We identify six SNe Ia with an early-time flux excess in our gr-band data (four ‘blue’ and two ‘red’ flux excesses). We find an intrinsic rate of 18 ± 11 percent of early flux excesses in SNe Ia at z < 0.07, based on three detected flux excesses out of 30 (10 percent) observed SNe Ia with a simulated efficiency of 57 percent. This is comparable to rates of flux excesses in the literature but also accounts for detection efficiencies. Two of these events are mostly consistent with circumstellar material interaction, while the other four have longer lifetimes in agreement with companion interaction and 56 Ni-clump models. We find a higher frequency of flux excesses in 91T/99aa-like events (44 ± 13 percent).

79 ASTRONOMY AND ASTROPHYSICS↗

A Simulation-based Method for Correcting Mode Coupling in CMB Angular Power Spectra

Modern cosmic microwave background (CMB) analysis pipelines regularly employ complex time-domain filters, beam models, masking, and other techniques during the production of sky maps and their corresponding angular power spectra. However, these processes can generate couplings between multipoles from the same spectrum and from different spectra, in addition to the typical power attenuation. Within the context of pseudo-C ℓ based, MASTER-style analyses, the net effect of the time-domain filtering is commonly approximated by a multiplicative transfer function, F ℓ , that can fail to capture mode mixing and is dependent on the spectrum of the signal. To address these shortcomings, we have developed a simulation-based spectral correction approach that constructs a two-dimensional transfer matrix, ${J}_{{\ell }{\ell }^{\prime} }$, which contains information about mode mixing in addition to mode attenuation. We demonstrate the application of this approach on data from the first flight of the Spider balloon-borne CMB experiment.

79 ASTRONOMY AND ASTROPHYSICS↗

Data Quality Monitoring for the Hadron Calorimeters Using Transfer Learning for Anomaly Detection

The proliferation of sensors brings an immense volume of spatio-temporal (ST) data in many domains, including monitoring, diagnostics, and prognostics applications. Data curation is a time-consuming process for a large volume of data, making it challenging and expensive to deploy data analytics platforms in new environments. Transfer learning (TL) mechanisms promise to mitigate data sparsity and model complexity by utilizing pre-trained models for a new task. Despite the triumph of TL in fields like computer vision and natural language processing, efforts on complex ST models for anomaly detection (AD) applications are limited. In this study, we present the potential of TL within the context of high-dimensional ST AD with a hybrid autoencoder architecture, incorporating convolutional, graph, and recurrent neural networks. Motivated by the need for improved model accuracy and robustness, particularly in scenarios with limited training data on systems with thousands of sensors, this research investigates the transferability of models trained on different sections of the Hadron Calorimeter of the Compact Muon Solenoid experiment at CERN. The key contributions of the study include exploring TL’s potential and limitations within the context of encoder and decoder networks, revealing insights into model initialization and training configurations that enhance performance while substantially reducing trainable parameters and mitigating data contamination effects.

47 OTHER INSTRUMENTATION↗

The Simons Observatory: Beam Characterization for the Small Aperture Telescopes

Abstract We use time-domain simulations of Jupiter observations to test and develop a beam reconstruction pipeline for the Simons Observatory Small Aperture Telescopes. The method relies on a mapmaker that estimates and subtracts correlated atmospheric noise and a beam fitting code designed to compensate for the bias caused by the mapmaker. We test our reconstruction performance for four different frequency bands against various algorithmic parameters, atmospheric conditions, and input beams. We additionally show the reconstruction quality as a function of the number of available observations and investigate how different calibration strategies affect the beam uncertainty. For all of the cases considered, we find good agreement between the fitted results and the input beam model within an ∼1.5% error for a multipole range ℓ = 30–700 and an ∼0.5% error for a multipole range ℓ = 50–200. We conclude by using a harmonic-domain component separation algorithm to verify that the beam reconstruction errors and biases observed in our analysis do not significantly bias the Simons Observatory r -measurement

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Brighter-fatter Effect in Near-infrared Detectors—III. Fourier-domain Treatment of Flat Field Correlations and Application to WFIRST

Weak gravitational lensing has emerged as a leading probe of the growth of cosmic structure. However, the shear signal is very small and accurate measurement depends critically on our ability to understand how non-ideal instrumental effects affect astronomical images. The Wide-Field Infrared Survey Telescope (WFIRST) will fly a focal plane containing 18 Teledyne H4RG-10 near-infrared detector arrays, which present different instrument calibration challenges from previous weak lensing observations. Previous work [Paper I: Hirata & Choi, PASP, 132, 014501 (2020); and Paper II: Choi & Hirata, PASP, 132, 014502 (2020)] has shown that correlation functions of flat field images, including cross-correlations between different time slices that are enabled by the non-destructive read capability of the infrared detectors, are effective tools for disentangling linear and nonlinear inter-pixel capacitance (IPC) and the brighter-fatter effect (BFE). Here we present a Fourier-domain treatment of the flat field correlations, which allows us to expand the previous formalism to all orders in IPC, BFE, and classical nonlinearity. We show that biases in simulated flat field analyses in Paper I are greatly reduced through the use of this formalism. We then apply this updated formalism to flat field data from three WFIRST flight candidate detectors, and explore the robustness to variations in the analysis. We find that the BFE is present in all three detectors, and that its contribution to the flat field correlations dominates over the nonlinear IPC, in accordance with the results from Paper II on a development detector. The magnitude of the BFE is such that the effective area of a pixel is increased by (3.54 ± 0.03) × 10 -7 for every electron deposited in a neighboring pixel (sensor chip assembly [SCA] 20829, statistical error, not IPC-deconvolved). We compare IPC maps from flat field autocorrelation measurements to those obtained from the single pixel reset method and find a median difference of 0.113% for SCA 20829. After further diagnosis of this difference, we ascribe it largely to an additional source of cross-talk, the vertical trailing pixel effect, and recommend further work to develop a model for this effect. These results represent a significant step toward calibration of the non-ideal effects in WFIRST detectors.

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↗

Deep drilling in the time domain with DECam: survey characterization

ABSTRACT This paper presents a new optical imaging survey of four deep drilling fields (DDFs), two Galactic and two extragalactic, with the Dark Energy Camera (DECam) on the 4-m Blanco telescope at the Cerro Tololo Inter-American Observatory (CTIO). During the first year of observations in 2021, >4000 images covering 21 deg2 (seven DECam pointings), with ∼40 epochs (nights) per field and 5 to 6 images per night per filter in g, r, i, and/or z have become publicly available (the proprietary period for this program is waived). We describe the real-time difference-image pipeline and how alerts are distributed to brokers via the same distribution system as the Zwicky Transient Facility (ZTF). In this paper, we focus on the two extragalactic deep fields (COSMOS and ELAIS-S1) characterizing the detected sources, and demonstrating that the survey design is effective for probing the discovery space of faint and fast variable and transient sources. We describe and make publicly available 4413 calibrated light curves based on difference-image detection photometry of transients and variables in the extragalactic fields. We also present preliminary scientific analysis regarding the Solar system small bodies, stellar flares and variables, Galactic anomaly detection, fast-rising transients and variables, supernovae, and active Galactic nuclei.

79 ASTRONOMY AND ASTROPHYSICS↗

Astrometric microlensing probes of the isolated neutron star population with Roman

Notoriously hard to detect and study, isolated neutron stars (NSs) might provide valuable answers to fundamental questions about stellar evolution and explosion physics. With the upcoming Roman Space Telescope, scheduled for launch in 2026, a new and powerful channel for their detection will become available: astrometric microlensing. We set out to create a realistic sample of simulated gravitational microlensing events as observed by Roman with the Galactic Bulge Time Domain Survey. We focus in particular on the population of NS lenses, which has until now been largely understudied. We used dedicated Galactic models tailored for application to microlensing by compact objects. In addition to populations of stars, white dwarfs, and black holes, we simulated four different NS populations with Maxwellian natal kick distributions: v = (150, 250, 350, 450) km/s. For each simulation, we applied projected Roman precision, cadence, and detectability criteria. We found that the parameter space log10 t E –log 10 θ E , which will be accessible to Roman observations, is efficient for the classification of stellar remnants. We found a feature in this space that is characteristic of NSs; using this feature, optimal samples of NS candidates can be constructed from Roman -like datasets. We describe the dependence of the observable parameter distributions on the assumed mean kick velocities. As the effects of natal kicks are very complex and mutually counteracting, we suggest that more detailed studies focused on the dynamics of NSs are needed in anticipation of Roman and future surveys. We estimate that Roman will observe approximately 11 000 microlensing events, including ~100 with NS lenses, whose photometric and astrometric signals are detectable; the event yield decreases by 38% when gap-filling low-cadence observations are not included. We make all simulated microlensing event datasets publicly available in preparation for Roman data.

79 ASTRONOMY AND ASTROPHYSICS↗

LAMOST Time-Domain survey: first results of four K2 plates

From Oct. 2019 to Apr. 2020, LAMOST performed a time-domain (TD) spectroscopic survey of four K2 plates with both low- and medium-resolution observations. The low-resolution spectroscopic survey acquired 282 exposures ( ≈ 46.6 h) over 25 nights, yielding a total of about 767 000 spectra, and the medium-resolution survey took 177 exposures ( ≈ 49.1 h) over 27 nights, collecting about 478 000 spectra. More than 70%/50% of low-resolution/medium-resolution spectra have signal-to-noise ratio higher than 10. We determine stellar parameters (e.g., T eff , log g, [Fe/H]) and radial velocity (RV) with different methods, including LASP, DD-Payne and SLAM. In general, these parameter estimations from different methods show good agreement, and the stellar parameter values are consistent with those of APOGEE. We use the Gaia DR2 RV values to calculate a median RV zero point (RVZP) for each spectrograph exposure by exposure, and the RVZP-corrected RVs agree well with the APOGEE data. The stellar evolutionary and spectroscopic masses are estimated based on the stellar parameters, multi-band magnitudes, distances and extinction values. Finally, we construct a binary catalog including about 2700 candidates by analyzing their light curves, fitting the RV data, calculating the binarity parameters from medium-resolution spectra and cross-matching the spatially resolved binary catalog from Gaia EDR3. In conclusion, the LAMOST TD survey is expected to represent a breakthrough in various scientific topics, such as binary systems, stellar activity, stellar pulsation, etc.

79 ASTRONOMY AND ASTROPHYSICS↗

The LSST DESC DC2 Simulated Sky Survey

Here, we describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep area of approximately 300 deg 2 , as well as a deep drilling field of approximately 1 deg 2 . We simulate 5 yr of the planned 10 yr survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the data set to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic test bed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time domain cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

J-PLUS: Support vector regression to measure stellar parameters

Stellar parameters are among the most important characteristics in studies of stars which, in traditional methods, are based on atmosphere models. However, time, cost, and brightness limits restrain the efficiency of spectral observations. The Javalambre Photometric Local Universe Survey (J-PLUS) is an observational campaign that aims to obtain photometry in 12 bands. Owing to its characteristics, J-PLUS data have become a valuable resource for studies of stars. Machine learning provides powerful tools for efficiently analyzing large data sets, such as the one from J-PLUS, and enables us to expand the research domain to stellar parameters. The main goal of this study is to construct a support vector regression (SVR) algorithm to estimate stellar parameters of the stars in the first data release of the J-PLUS observational campaign. The training data for the parameter's regressions are featured with 12-waveband photometry from J-PLUS and are crossidentified with spectrum-based catalogs. These catalogs are from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, the Apache Point Observatory Galactic Evolution Experiment, and the Sloan Extension for Galactic Understanding and Exploration. We then label them with the stellar effective temperature, the surface gravity, and the metallicity. Ten percent of the sample is held out to apply a blind test. We develop a new method, a multi-model approach, in order to fully take into account, the uncertainties of both the magnitudes and the stellar parameters. The method utilizes more than 200 models to apply the uncertainty analysis. We present a catalog of 2 493 424 stars with the root mean square error of 160 K in the effective temperature regression, 0.35 in the surface gravity regression, and 0.25 in the metallicity regression. We also discuss the advantages of this multi-model approach and compare it to other machine-learning methods.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimization of the Observing Cadence for the Rubin Observatory Legacy Survey of Space and Time: A Pioneering Process of Community-focused Experimental Design

Vera C. Rubin Observatory is a ground-based astronomical facility under construction, a joint project of the National Science Foundation and the U.S. Department of Energy, designed to conduct a multipurpose 10 yr optical survey of the Southern Hemisphere sky: the Legacy Survey of Space and Time. Significant flexibility in survey strategy remains within the constraints imposed by the core science goals of probing dark energy and dark matter, cataloging the solar system, exploring the transient optical sky, and mapping the Milky Way. The survey’s massive data throughput will be transformational for many other astrophysics domains and Rubin’s data access policy sets the stage for a huge community of potential users. To ensure that the survey science potential is maximized while serving as broad a community as possible, Rubin Observatory has involved the scientific community at large in the process of setting and refining the details of the observing strategy. The motivation, history, and decision-making process of this strategy optimization are detailed in this paper, giving context to the science-driven proposals and recommendations for the survey strategy included in this Focus Issue.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

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↗