Search

We report on a search for millicharged particles (mCPs) produced in cosmic ray proton atmospheric interactions using data collected during the first science run of the LUX-ZEPLIN experiment. The mCPs produced by two processes -- meson decay and proton bremsstrahlung -- are considered in this study. This search utilized a novel signature unique to liquid xenon (LXe) time projection chambers (TPCs), allowing sensitivity to mCPs with masses ranging from 10 to 1000 MeV/c$^2$ and fractional charges between 0.001 and 0.02 of the electron charge e. With an exposure of 60 live days and a 5.5 tonne fiducial mass, we observed no significant excess over background. This represents the first experimental search for atmospheric mCPs and the first search for mCPs using an underground LXe experiment.

We investigate the emergence of self-hybridized thermal magnetoplasmons in doped graphene nanodisks at finite temperatures when subjected to an external magnetic field. Using a semi-analytical approach, which fully describes the eigenmodes and polarizability of the graphene nanodisks, we show that the hybridization originates from the coupling of transitions between thermally populated Landau levels and localized magnetoplasmon resonances of the nanodisks. Owing to their origin, these modes combine the extraordinary magneto-optical response of graphene with the strong field enhancement of plasmons, making them an ideal tool for achieving strong chiral light-matter interactions, with the additional advantage of being tunable through carrier concentration, magnetic field, and temperature. As an illustration of their capabilities, we demonstrate that the thermal magnetoplasmons supported by an array of graphene nanodisks enable chiral perfect absorption and chiral thermal emission., Comment: 18 pages, 8 figures

In this paper, a recently developed snapshot hyperspectral imaging (HSI) system based on Computed Tomography Imaging Spectroscopy (CTIS) is utilized to determine Brix and pH values in Sheegene 20 table grapes through Partial Least Squares Regression (PLSR) modeling. The performance of the CTIS system is compared with that of a state-of-the-art line scan HSI system by imaging 100 grapes across both platforms. Reference measurements of Brix and pH values are obtained directly using a refractometer and a pH meter, as these parameters are essential for assessing the quality of table and wine grapes. The findings indicate that the spectra captured by the CTIS camera correlate well with the reference measurements, despite the system's narrower spectral range. The CTIS camera's advantages, including its lower cost, portability, and reduced susceptibility to motion errors, highlight its potential for promising in-field applications in grape quality assessment., Comment: 15 pages, 10 figures

This article introduces a groundbreaking analytical workflow designed for the holistic characterisation, modelling and physical simulation of device heterostructures. Our innovative workflow autonomously, comprehensively and locally characterises the crystallographic information and 3D orientation of the crystal phases, the elemental composition, and the strain maps of devices from (scanning) transmission electron microscopy data. It converts a manual characterisation process that traditionally takes days into an automatic routine completed in minutes. This is achieved through a physics-guided artificial intelligence model that combines unsupervised and supervised machine learning in a modular way to provide a representative 3D description of the devices, materials structures, or samples under analysis. To culminate the process, we integrate the extracted knowledge to automate the generation of both 3D finite element and atomic models of millions of atoms acting as digital twins, enabling simulations that yield essential physical and chemical insights crucial for understanding the device's behaviour in practical applications. We prove this end-to-end workflow with a state-of-the-art materials platform based on SiGe planar heterostructures for hosting coherent and scalable spin qubits. Our workflow connects representative digital twins of the experimental devices with their theoretical properties to reveal the true impact that every atom in the structure has on their electronic properties, and eventually, into their functional quantum performance. Notably, the versatility of our workflow is demonstrated through its successful application to a wide array of materials systems, device configurations and sample morphologies.

The XLZD collaboration is developing a two-phase xenon time projection chamber with an active mass of 60 to 80 t capable of probing the remaining WIMP-nucleon interaction parameter space down to the so-called neutrino fog. In this work we show that, based on the performance of currently operating detectors using the same technology and a realistic reduction of radioactivity in detector materials, such an experiment will also be able to competitively search for neutrinoless double beta decay in $^{136}$Xe using a natural-abundance xenon target. XLZD can reach a 3$\sigma$ discovery potential half-life of 5.7$\times$10$^{27}$ yr (and a 90% CL exclusion of 1.3$\times$10$^{28}$ yr) with 10 years of data taking, corresponding to a Majorana mass range of 7.3-31.3 meV (4.8-20.5 meV). XLZD will thus exclude the inverted neutrino mass ordering parameter space and will start to probe the normal ordering region for most of the nuclear matrix elements commonly considered by the community., Comment: 29 pages, 7 figures

This report describes the experimental strategy and technologies for a next-generation xenon observatory sensitive to dark matter and neutrino physics. The detector will have an active liquid xenon target mass of 60-80 tonnes and is proposed by the XENON-LUX-ZEPLIN-DARWIN (XLZD) collaboration. The design is based on the mature liquid xenon time projection chamber technology of the current-generation experiments, LZ and XENONnT. A baseline design and opportunities for further optimization of the individual detector components are discussed. The experiment envisaged here has the capability to explore parameter space for Weakly Interacting Massive Particle (WIMP) dark matter down to the neutrino fog, with a 3$\sigma$ evidence potential for the spin-independent WIMP-nucleon cross sections as low as $3\times10^{-49}\rm cm^2$ (at 40 GeV/c$^2$ WIMP mass). The observatory is also projected to have a 3$\sigma$ observation potential of neutrinoless double-beta decay of $^{136}$Xe at a half-life of up to $5.7\times 10^{27}$ years. Additionally, it is sensitive to astrophysical neutrinos from the atmosphere, sun, and galactic supernovae., Comment: 32 pages, 14 figures

We report results of a search for nuclear recoils induced by weakly interacting massive particle (WIMP) dark matter using the LUX-ZEPLIN (LZ) two-phase xenon time projection chamber. This analysis uses a total exposure of $4.2\pm0.1$ tonne-years from 280 live days of LZ operation, of which $3.3\pm0.1$ tonne-years and 220 live days are new. A technique to actively tag background electronic recoils from $^{214}$Pb $\beta$ decays is featured for the first time. Enhanced electron-ion recombination is observed in two-neutrino double electron capture decays of $^{124}$Xe, representing a noteworthy new background. After removal of artificial signal-like events injected into the data set to mitigate analyzer bias, we find no evidence for an excess over expected backgrounds. World-leading constraints are placed on spin-independent (SI) and spin-dependent WIMP-nucleon cross sections for masses $\geq$9 GeV/$c^2$. The strongest SI exclusion set is $2.1\times10^{-48}$ cm$^{2}$ at the 90% confidence level at a mass of 36 GeV/$c^2$, and the best SI median sensitivity achieved is $5.0\times10^{-48}$ cm$^{2}$ for a mass of 40 GeV/$c^2$., Comment: 9 pages, 7 figures. See https://www.hepdata.net/record/155182 for a data release related to this paper

Machine-learning models in chemistry - when based on descriptors of atoms embedded within molecules - face essential challenges in transferring the quality of predictions of local electronic structures and their associated properties across chemical compound space. In the present work, we make use of adversarial validation to elucidate certain intrinsic complications related to machine inferences of unseen chemistry. On this basis, we employ invariant and equivariant neural networks - both trained either exclusively on total molecular energies or a combination of these and data from atomic partitioning schemes - to evaluate how such models scale performance-wise between datasets of fundamentally different functionality and composition. We find the inference of local electronic properties to improve significantly when training models on augmented data that appropriately expose local functional features. However, molecular datasets for training purposes must themselves be sufficiently comprehensive and rich in composition to warrant any generalizations to larger systems, and even then, transferability can still only genuinely manifest if the body of atomic energies available for training purposes exposes the uniqueness of different functional moieties within molecules. We demonstrate this point by comparing machine models trained on atomic partitioning schemes based on the spatial locality of either native atomic or molecular orbitals., Comment: 23+7 pages, 3 tables, 6 figures. SI (as PDF) and training configurations (as accompanying YAML files) may be found as ancillary files

Social networks shape individuals' lives, influencing everything from career paths to health. This paper presents a registry-based, multi-layer and temporal network of the entire Danish population in the years 2008-2021 (roughly 7.2 mill. individuals). Our network maps the relationships formed through family, households, neighborhoods, colleagues and classmates. We outline key properties of this multiplex network, introducing both an individual-focused perspective as well as a bipartite representation. We show how to aggregate and combine the layers, and how to efficiently compute network measures such as shortest paths in large administrative networks. Our analysis reveals how past connections reappear later in other layers, that the number of relationships aggregated over time reflects the position in the income distribution, and that we can recover canonical shortest path length distributions when appropriately weighting connections. Along with the network data, we release a Python package that uses the bipartite network representation for efficient analysis.

ANNZ is a fast and simple algorithm which utilises artificial neural networks (ANNs), it was known as one of the pioneers of machine learning approaches to photometric redshift estimation decades ago. We enhanced the algorithm by introducing new activation functions like tanh, softplus, SiLU, Mish and ReLU variants; its new performance is then vigorously tested on legacy samples like the Luminous Red Galaxy (LRG) and Stripe-82 samples from SDSS, as well as modern galaxy samples like the Physics of the Accelerating Universe Survey (PAUS). This work focuses on testing the robustness of activation functions with respect to the choice of ANN architectures, particularly on its depth and width, in the context of galaxy photometric redshift estimation. Our upgraded algorithm, which we named ANNZ+, shows that the tanh and Leaky ReLU activation functions provide more consistent and stable results across deeper and wider architectures with > 1 per cent improvement in root-mean-square error ($\sigma_{\textrm{RMS}}$) and 68th percentile error ($\sigma_{68}$) when tested on SDSS data sets. While assessing its capabilities in handling high dimensional inputs, we achieved an improvement of 11 per cent in $\sigma_{\textrm{RMS}}$ and 6 per cent in $\sigma_{68}$ with the tanh activation function when tested on the 40-narrowband PAUS dataset; it even outperformed ANNZ2, its supposed successor, by 44 per cent in $\sigma_{\textrm{RMS}}$. This justifies the effort to upgrade the 20-year-old ANNZ, allowing it to remain viable and competitive within the photo-z community today. The updated algorithm ANNZ+ is publicly available at https://github.com/imdadmpt/ANNzPlus., Comment: 38 pages, 9 figures, revised version submitted to JCAP

The broad physics reach of the LUX-ZEPLIN (LZ) experiment covers rare phenomena beyond the direct detection of dark matter. We report precise measurements of the extremely rare decay of $^{124}$Xe through the process of two-neutrino double electron capture (2$\nu$2EC), utilizing a $1.39\,\mathrm{kg} \times \mathrm{yr}$ isotopic exposure from the first LZ science run. A half-life of $T_{1/2}^{2\nu2\mathrm{EC}} = (1.09 \pm 0.14_{\text{stat}} \pm 0.05_{\text{sys}}) \times 10^{22}\,\mathrm{yr}$ is observed with a statistical significance of $8.3\,\sigma$, in agreement with literature. First empirical measurements of the KK capture fraction relative to other K-shell modes were conducted, and demonstrate consistency with respect to recent signal models at the $1.4\,\sigma$ level., Comment: 15 pages, 3 figures

We present photometric redshifts for 1 341 559 galaxies from the Physics of the Accelerating Universe Survey (PAUS) over 50.38 ${\rm deg}^{2}$ of sky to $i_{\rm AB}=23$. Redshift estimation is performed using DEEPz, a deep-learning photometric redshift code. We analyse the photometric redshift precision when varying the photometric and spectroscopic samples. Furthermore, we examine observational and instrumental effects on the precision of the photometric redshifts, and we compare photometric redshift measurements with those obtained using a template method-fitting BCNz2. Finally, we examine the use of photometric redshifts in the identification of close galaxy pairs. We find that the combination of samples from W1+W3 in the training of DEEPz significantly enhances the precision of photometric redshifts. This also occurs when we recover narrow band fluxes using broad bands measurements. We show that DEEPz determines the redshifts of galaxies in the prevailing spectroscopic catalogue used in the training of DEEPz with greater precision. For the faintest galaxies ($i_{\rm AB}=21-23$), we find that DEEPz improves over BCNz2 both in terms of the precision (20-50 per cent smaller scatter) and in returning a smaller outlier fraction in two of the wide fields. The catalogues were tested for the identification of close galaxy pairs, showing that DEEPz is effective for the identification of close galaxy pairs for samples with $i_{\rm AB} < 22.5$ and redshift $0.2 < z < 0.6$. In addition, identifying close galaxy pairs common between DEEPz and BCNz2 is a promising approach to improving the purity of the catalogues of these systems., Comment: 21 pages, 13 figures, submitted to A&A, link to the catalogue

We present an improved zodiacal light model for COBE-DIRBE derived through global Bayesian analysis within the Cosmoglobe Data Release 2 framework. The parametric form of the ZL model is identical to that introduced by Kelsall et al. (1998), but the specific best-fit parameter values are re-derived using the combination of DIRBE Calibrated Individual Observations, Planck HFI sky maps, and WISE and Gaia compact object catalogs. Furthermore, the ZL parameters are fitted jointly with astrophysical parameters, such as thermal dust and starlight emission, and the new model takes into account excess radiation that appears stationary in solar-centric coordinates as reported in a companion paper. The relative differences between the predicted signals from K98 and our new model are $\lesssim 5\%$ in the 12 and 25 $\mu$m channels over the full sky. The zero-levels of the cleaned DR2 maps are lower than those of the K98 Zodiacal light Subtracted Mission Average maps by $\sim 10$ kJy/sr at 1.25-3.5 $\mu$m, which is larger than the entire predicted contribution from high-redshift galaxies to the Cosmic Infrared Background at the same wavelengths. The total RMS of each DR2 map at wavelengths up to and including 25 $\mu$m are $\sim 30$ $\%$ lower at high Galactic latitudes than the corresponding DIRBE ZSMA maps. Still, obvious ZL residuals can be seen in several of the DR2 maps, and further work is required to mitigate these. Joint analysis with existing and future high-resolution full-sky surveys such as AKARI, IRAS, Planck HFI, and SPHEREx will be essential both to break key degeneracies in the current model and to determine whether the reported solar-centric excess radiation has a ZL or instrumental origin. Thus, while the results presented in this paper do redefine the state-of-the-art for DIRBE modelling, it also only represents the first among many steps toward a future optimal Bayesian ZL model., Comment: 27 pages, 21 figures, submitted to A&A

We present the first global Bayesian analysis of the time-ordered DIRBE data within the Cosmoglobe framework, building on the same methodology that has previously been successfully applied to Planck LFI and WMAP. These data are analyzed jointly with COBE-FIRAS, Gaia, Planck HFI, and WISE, allowing for more accurate instrumental and astrophysical characterization than possible through single-experiment analysis. This paper provides an overview of the analysis pipeline and main results, and we present and characterize a new set of zodiacal light subtracted mission average (ZSMA) DIRBE maps spanning 1.25 to 240 $\mu$m. A novel aspect of this processing is the characterization and removal of excess radiation between 4.9 and 60$\,\mu$m that appears static in solar-centric coordinates. The DR2 ZSMA maps have several advantages with respect to the previously available maps, including 1) lower zodiacal light (and possibly straylight) residuals; 2) better determined zero-levels; 3) natively HEALPix tessellated maps with a $7'$ pixel size; 4) nearly white noise at pixel scales; and 5) a more complete and accurate noise characterization established through the combination of MCMC samples and half-mission maps. In addition, because the model has been simultaneously fitted with both DIRBE and HFI data, this is the first consistent unification of the infrared and CMB wavelength ranges into one global sky model covering 100 GHz to 1 $\mu$m. However, even though the new maps are improved with respect to the official maps, and should be preferred for most future analyses that require DIRBE sky maps, they still exhibit non-negligible zodiacal light residuals between 12 and 60$\,\mu$m. Further improvements should be made through joint analysis with complementary infrared experiments such IRAS, AKARI, WISE and SPHEREx, and releasing the full combined potential of all these powerful infrared observatories., Comment: 30 pages, 23 figures, submitted to A&A. Comments welcome! Data can be found at https://www.cosmoglobe.uio.no/products/cosmoglobe-dr2.html

SPIDER is a balloon-borne instrument designed to map the cosmic microwave background at degree-angular scales in the presence of Galactic foregrounds. SPIDER has mapped a large sky area in the Southern Hemisphere using more than 2000 transition-edge sensors (TESs) during two NASA Long Duration Balloon flights above the Antarctic continent. During its first flight in January 2015, SPIDER observed in the 95 GHz and 150 GHz frequency bands, setting constraints on the B-mode signature of primordial gravitational waves. Its second flight in the 2022-23 season added new receivers at 280 GHz, each using an array of TESs coupled to the sky through feedhorns formed from stacks of silicon wafers. These receivers are optimized to produce deep maps of polarized Galactic dust emission over a large sky area, providing a unique data set with lasting value to the field. In this work, we describe the instrument's performance during SPIDER's second flight., Comment: Submitted to SPIE Astronomical Telescopes + Instrumentation 2024, JATIS

Large angular scale surveys in the absence of atmosphere are essential for measuring the primordial $B$-mode power spectrum of the Cosmic Microwave Background (CMB). Since this proposed measurement is about three to four orders of magnitude fainter than the temperature anisotropies of the CMB, in-flight calibration of the instruments and active suppression of systematic effects are crucial. We investigate the effect of changing the parameters of the scanning strategy on the in-flight calibration effectiveness, the suppression of the systematic effects themselves, and the ability to distinguish systematic effects by null-tests. Next-generation missions such as LiteBIRD, modulated by a Half-Wave Plate (HWP), will be able to observe polarisation using a single detector, eliminating the need to combine several detectors to measure polarisation, as done in many previous experiments and hence avoiding the consequent systematic effects. While the HWP is expected to suppress many systematic effects, some of them will remain. We use an analytical approach to comprehensively address the mitigation of these systematic effects and identify the characteristics of scanning strategies that are the most effective for implementing a variety of calibration strategies in the multi-dimensional space of common spacecraft scan parameters. We also present Falcons, a fast spacecraft scanning simulator that we developed to investigate this scanning parameter space.

Local electronic-structure methods in quantum chemistry operate on the ability to compress electron correlations more efficiently in a basis of spatially localized molecular orbitals than in a parent set of canonical orbitals. However, many typical choices of localized orbitals tend to be related by selected, near-exact symmetry operations whenever a molecule belongs to a point group, a feature which remains largely unexploited in most local correlation methods. The present Letter demonstrates how to leverage a recent unitary protocol for enforcing symmetry properties among localized orbitals to yield a high-accuracy estimate of the exact ground-state correlation energy of benzene ($D_{6h}$) in correlation-consistent polarized basis sets of both double- and triple-$\zeta$ quality. Through an initial application to many-body expanded full configuration interaction (MBE-FCI) theory, we show how molecular point-group symmetry can lead to computational savings that are inversely proportional to the order of the point group in a manner generally applicable to the acceleration of modern local correlation methods., Comment: 13+5 pages, 5 figures. SI as an ancillary file

Using data from the first flight of SPIDER and from Planck HFI, we probe the properties of polarized emission from interstellar dust in the SPIDER observing region. Component separation algorithms operating in both the spatial and harmonic domains are applied to probe their consistency and to quantify modeling errors associated with their assumptions. Analyses spanning the full SPIDER region demonstrate that i) the spectral energy distribution of diffuse Galactic dust emission is broadly consistent with a modified-blackbody (MBB) model with a spectral index of $\beta_\mathrm{d}=1.45\pm0.05$ $(1.47\pm0.06)$ for $E$ ($B$)-mode polarization, slightly lower than that reported by Planck for the full sky; ii) its angular power spectrum is broadly consistent with a power law; and iii) there is no significant detection of line-of-sight decorrelation of the astrophysical polarization. The size of the SPIDER region further allows for a statistically meaningful analysis of the variation in foreground properties within it. Assuming a fixed dust temperature $T_\mathrm{d}=19.6$ K, an analysis of two independent sub-regions of that field results in inferred values of $\beta_\mathrm{d}=1.52\pm0.06$ and $\beta_\mathrm{d}=1.09\pm0.09$, which are inconsistent at the $3.9\,\sigma$ level. Furthermore, a joint analysis of SPIDER and Planck 217 and 353 GHz data within a subset of the SPIDER region is inconsistent with a simple MBB at more than $3\,\sigma$, assuming a common morphology of polarized dust emission over the full range of frequencies. These modeling uncertainties have a small--but non-negligible--impact on limits on the cosmological tensor-to-scalar ratio derived from the \spider dataset. The fidelity of the component separation approaches of future CMB polarization experiments may thus have a significant impact on their constraining power., Comment: 21 pages, 15 figures

In forensic genetics, short tandem repeats (STRs) are used for human identification (HID). Degraded biological trace samples with low amounts of short DNA fragments (low-quality DNA samples) pose a challenge for STR typing. Predefined single nucleotide polymorphisms (SNPs) can be amplified on short PCR fragments and used to generate SNP profiles from low-quality DNA samples. However, the stochastic results from low-quality DNA samples may result in frequent locus drop-outs and insufficient numbers of SNP genotypes for convincing identification of individuals. Shotgun DNA sequencing potentially analyses all DNA fragments in a sample in contrast to the targeted PCR-based sequencing methods and may be applied to DNA samples of very low quality, like heavily compromised crime-scene samples and ancient DNA samples. Here, we developed a statistical model for shotgun sequencing, sequence alignment, and genotype calling. Results from replicated shotgun sequencing of buccal swab (high-quality samples) and hair samples (low-quality samples) were arranged in a genotype-call confusion matrix to estimate the calling error probability by maximum likelihood and Bayesian inference. We developed formulas for calculating the evidential weight as a likelihood ratio (LR) based on data from dynamically selected SNPs from shotgun DNA sequencing. The method accounts for potential genotyping errors. Different genotype quality filters may be applied to account for genotyping errors. An error probability of zero resulted in the forensically commonly used LR formula. When considering a single SNP marker's contribution to the LR, error probabilities larger than zero reduced the LR contribution of matching genotypes and increased the LR in the case of a mismatch. We developed an open-source R package, wgsLR, which implements the method, including estimating the calling error probability and calculating LR values., Comment: 25 pages, 9 figures

We present the Lyman-$\alpha$ (Ly$\alpha$) and ultraviolet (UV) luminosity function (LF), in bins of redshift, of quasars selected in the Physics of the Accelerating Universe Survey (PAUS). A sample of 915 objects was selected at $2.710^{44}$ ${\rm erg\,s}$$^{-1}$. In order to obtain corrections for the LF estimation, and assess the accuracy of our selection method, we produced mock catalogs of $0

We assess the capabilities of the LiteBIRD mission to map the hot gas distribution in the Universe through the thermal Sunyaev-Zeldovich (SZ) effect. Our analysis relies on comprehensive simulations incorporating various sources of Galactic and extragalactic foreground emission, while accounting for specific instrumental characteristics of LiteBIRD, such as detector sensitivities, frequency-dependent beam convolution, inhomogeneous sky scanning, and $1/f$ noise. We implement a tailored component-separation pipeline to map the thermal SZ Compton $y$-parameter over 98% of the sky. Despite lower angular resolution for galaxy cluster science, LiteBIRD provides full-sky coverage and, compared to the Planck satellite, enhanced sensitivity, as well as more frequency bands to enable the construction of an all-sky $y$-map, with reduced foreground contamination at large and intermediate angular scales. By combining LiteBIRD and Planck channels in the component-separation pipeline, we obtain an optimal $y$-map that leverages the advantages of both experiments, with the higher angular resolution of the Planck channels enabling the recovery of compact clusters beyond the LiteBIRD beam limitations, and the numerous sensitive LiteBIRD channels further mitigating foregrounds. The added value of LiteBIRD is highlighted through the examination of maps, power spectra, and one-point statistics of the various sky components. After component separation, the $1/f$ noise from LiteBIRD is effectively mitigated below the thermal SZ signal at all multipoles. Cosmological constraints on $S_8=\sigma_8\left(\Omega_{\rm m}/0.3\right)^{0.5}$ obtained from the LiteBIRD-Planck combined $y$-map power spectrum exhibits a 15% reduction in uncertainty compared to constraints from Planck alone. This improvement can be attributed to the increased portion of uncontaminated sky available in the LiteBIRD-Planck combined $y$-map., Comment: 38 pages, 13 figures, abstract shortened. Updated to match version accepted by JCAP

We study a large family of generalized class groups of imaginary quadratic orders $O$ and prove that they act freely and (essentially) transitively on the set of primitively $O$-oriented elliptic curves over a field $k$ (assuming this set is non-empty) equipped with appropriate level structure. This extends, in several ways, a recent observation due to Galbraith, Perrin and Voloch for the ray class group. We show that this leads to a reinterpretation of the action of the class group of a suborder $O' \subseteq O$ on the set of $O'$-oriented elliptic curves, discuss several other examples, and briefly comment on the hardness of the corresponding vectorization problems., Comment: Paper accepted by the International Workshop on the Arithmetic of Finite Fields 2024. Comments welcome. 18 pages

We present a new theory for partitioning simulations of periodic and solid-state systems into physically sound atomic contributions at the level of Kohn-Sham density functional theory. Our theory is based on spatially localized linear combinations of crystalline Gaussian-type orbitals and, as such, capable of exposing local features within periodic electronic structures in a more intuitive and robust manner than alternatives based on the spatial distribution of atomic basis functions alone. Early decomposed cohesive energies of both molecular polymers and different crystalline polymorphs demonstrate how the atomic properties yielded by our theory convincingly align with the expected charge polarization in these systems, also whenever partial charges and Madelung energies may lend themselves somewhat ambiguous to interpretation., Comment: 24+8 pages, 4 figures. SI as an ancillary file

Wide-Area Linear Optical Polarimeter (WALOP)-South is the first wide-field and survey-capacity polarimeter in the optical wavelengths. On schedule for commissioning in 2024, it will be mounted on the 1 m SAAO telescope in Sutherland Observatory, South Africa to undertake the PASIPHAE sky survey. PASIPHAE program will create the first polarimetric sky map in the optical wavelengths, spanning more than 2000 square degrees of the southern Galactic region. The innovative design of WALOP-South will enable it to measure the linear polarization (Stokes parameters $q$ and $u$), in a single exposure, of all sources in a field of view (FoV) of $35\times35$ arcminutes-squared in the SDSS-r broadband and narrowband filters between 500-750 nm with 0.1 % polarization accuracy. The unique goals of the instrument place very stringent systems engineering goals, including on the performance of the optical, polarimetric, optomechanical, and electronic subsystems. All the subsystems have been designed carefully to meet the overall instrument performance goals. As of May 2024, all the instrument optical and mechanical subsystems have been assembled and are currently getting tested and integrated. The complete testing and characterization of the instrument in the lab is expected to be completed by August 2024. In this paper, we will present (a) the design and development of the entire instrument and its major subsystems, focusing on the opto-mechanical design which has not been reported before, and (b) assembly and integration of the instrument in the lab and early results from lab characterization of the instrument., Comment: 20 pages; presented at 2024 SPIE Astronomical Telescopes + Instrumentation conference

We present a wide-reaching revamp of the generalized many-body expanded full configuration interaction (MBE-FCI) method. First, we outline how to automatize the selection of reference active spaces whereby the inherent bias introduced through a manual identification is reduced, also within the context of traditional complete active space methods. Second, we allow for the use of compact orbital clusters as expansion objects, which works to circumvent the unfavorable scaling with the number of orbitals included in the space complementary to the reference orbitals. Finally, we present a new algorithm for ensuring that many-body expansions can be efficiently terminated while conservatively accounting for resulting errors. These developments are all tested on a variety of molecular systems and different orbital representations to illustrate the abilities of our algorithm to produce correlation energies within predetermined error bounds, significantly broadening the overall applicability of the MBE-FCI method., Comment: 29+7 pages, 8 figures, 1 table. SI as an ancillary file

The Carbon monOxide Mapping Array Project (COMAP) Pathfinder survey continues to demonstrate the feasibility of line-intensity mapping using high-redshift carbon monoxide (CO) line emission traced at cosmological scales. The latest COMAP Pathfinder power spectrum analysis is based on observations through the end of Season 2, covering the first three years of Pathfinder operations. We use our latest constraints on the CO(1-0) line-intensity power spectrum at $z\sim3$ to update corresponding constraints on the cosmological clustering of CO line emission and thus the cosmic molecular gas content at a key epoch of galaxy assembly. We first mirror the COMAP Early Science interpretation, considering how Season 2 results translate to limits on the shot noise power of CO fluctuations and the bias of CO emission as a tracer of the underlying dark matter distribution. The COMAP Season 2 results place the most stringent limits on the CO tracer bias to date, at $\langle{Tb}\rangle<4.8$ $\mu$K. These limits narrow the model space significantly compared to previous CO line-intensity mapping results while maintaining consistency with small-volume interferometric surveys of resolved line candidates. The results also express a weak preference for CO emission models used to guide fiducial forecasts from COMAP Early Science, including our data-driven priors. We also consider directly constraining a model of the halo-CO connection, and show qualitative hints of capturing the total contribution of faint CO emitters through the improved sensitivity of COMAP data. With continued observations and matching improvements in analysis, the COMAP Pathfinder remains on track for a detection of cosmological clustering of CO emission., Comment: 9 pages + bibliography and appendices (13 pages total); 9 figures, 1 table; v2 reflects minor changes made for version submitted to A&A, with no changes to top-line results

We present updated constraints on the cosmological 3D power spectrum of carbon monoxide CO(1-0) emission in the redshift range $2.4$-$3.4$. The constraints are derived from the two first seasons of Carbon monOxide Mapping Array Project (COMAP) Pathfinder line-intensity mapping observations aiming to trace star-formation during the Epoch of Galaxy Assembly. These results improve on the previous Early Science (ES) results through both increased data volume and improved data processing methodology. On the methodological side, we now perform cross-correlations between groups of detectors (''feed-groups''), as opposed to cross-correlations between single feeds, and this new feed-group pseudo power spectrum (FGPXS) is constructed to be more robust against systematic effects. In terms of data volume, the effective mapping speed is significantly increased due to an improved observational strategy as well as better data selection methodology. The updated spherically- and field-averaged FGPXS, $\tilde{C}(k)$, is consistent with zero, at a probability-to-exceed of around $34\,\%$, with an excess of $2.7\,\sigma$ in the most sensitive bin. Our power spectrum estimate is about an order of magnitude more sensitive in our six deepest bins across ${0.09\,\mathrm{Mpc}^{-1} < k < 0.73\,\mathrm{Mpc}^{-1}}$, as compared to the feed-feed pseudo power spectrum (FPXS) of COMAP ES. Each of these bins individually constrains the CO power spectrum to ${kP_\mathrm{CO}(k)< 2400-4900\,\mathrm{\mu K^2 Mpc^{2}}}$ at $95\,\%$ confidence. To monitor potential contamination from residual systematic effects, we analyze a set of 312 difference-map null tests and find that these are consistent with the instrumental noise prediction. In sum, these results provide the strongest direct constraints on the cosmological 3D CO(1-0) power spectrum published to date., Comment: 17 pages, 10 figures, v2 reflects changes made for version submitted to Astronomy and Astrophysics, addition of additional figure clarifying methodology, no change to final results

The CO Mapping Array Project (COMAP) Pathfinder is performing line intensity mapping of CO emission to trace the distribution of unresolved galaxies at redshift $z \sim 3$. We present an improved version of the COMAP data processing pipeline and apply this to the first two seasons of observations. This analysis improves on the COMAP Early Science (ES) results in several key aspects. On the observational side, all second season scans were made in constant-elevation mode, after noting that the previous Lissajous scans were associated with increased systematic errors; those scans accounted for 50% of the total Season 1 data volume. Secondly, all new observations were restricted to an elevation range of 35-65 degrees, to minimize sidelobe ground pickup. On the data processing side, more effective data cleaning in both the time- and map-domain has allowed us to eliminate all data-driven power spectrum-based cuts. This increases the overall data retention and reduces the risk of signal subtraction bias. On the other hand, due to the increased sensitivity, two new pointing-correlated systematic errors have emerged, and we introduce a new map-domain PCA filter to suppress these. Subtracting only 5 out of 256 PCA modes, we find that the standard deviation of the cleaned maps decreases by 67% on large angular scales, and after applying this filter, the maps appear consistent with instrumental noise. Combining all these improvements, we find that each hour of raw Season 2 observations yields on average 3.2 times more cleaned data compared to ES analysis. Combining this with the increase in raw observational hours, the effective amount of data available for high-level analysis is a factor of 8 higher than in ES. The resulting maps have reached an uncertainty of $25$-$50\,\mu K$ per voxel, providing by far the strongest constraints on cosmological CO line emission published to date., Comment: 23 pages, 22 figures, for submission to Astronomy and Astrophysics

The Physics of the Accelerating Universe (PAU) camera is an optical narrow band and broad band imaging instrument mounted at the prime focus of the William Herschel Telescope. We describe the image calibration procedure of the PAU Survey data. We rely on an external photometric catalogue to calibrate our narrow band data using stars that have been observed by both datasets. We fit stellar templates to the stellar broad band photometry of the Sloan Digital Sky Survey and synthesise narrow band photometry that we compare to the PAUS narrow band data to determine their calibration. Consequently, the PAUS data are in the AB system as inherited from its reference calibrator. We do several tests to check the performance of the calibration. We find it self-consistent when comparing repeated observations of the same objects, with a good overall accuracy to the AB system which we estimate to be at the 2\% precision level and no significant trends as a function of narrow band filter or wavelength. Repeated observations allow us to build a spatial map of the illumination pattern of the system. We also check the wavelength dependence of the calibration comparing to stellar spectra. We find that using only blue stars reduces the effects of variations in the stellar template fitting to broad-band colours, improving the overall precision of the calibration to around 1\% and its wavelength uniformity. The photometric redshift performance obtained with the PAUS data attests to the validity of our calibration to reach the PAUS science goals., Comment: 18 pages, 20 figures, accepted for publication in MNRAS

Weakly interacting massive particles (WIMPs) may interact with a virtual pion that is exchanged between nucleons. This interaction channel is important to consider in models where the spin-independent isoscalar channel is suppressed. Using data from the first science run of the LUX-ZEPLIN dark matter experiment, containing 60 live days of data in a 5.5~tonne fiducial mass of liquid xenon, we report the results on a search for WIMP-pion interactions. We observe no significant excess and set an upper limit of $1.5\times10^{-46}$~cm$^2$ at a 90\% confidence level for a WIMP mass of 33~GeV/c$^2$ for this interaction.

LiteBIRD, the next-generation cosmic microwave background (CMB) experiment, aims for a launch in Japan's fiscal year 2032, marking a major advancement in the exploration of primordial cosmology and fundamental physics. Orbiting the Sun-Earth Lagrangian point L2, this JAXA-led strategic L-class mission will conduct a comprehensive mapping of the CMB polarization across the entire sky. During its 3-year mission, LiteBIRD will employ three telescopes within 15 unique frequency bands (ranging from 34 through 448 GHz), targeting a sensitivity of 2.2\,$\mu$K-arcmin and a resolution of 0.5$^\circ$ at 100\,GHz. Its primary goal is to measure the tensor-to-scalar ratio $r$ with an uncertainty $\delta r = 0.001$, including systematic errors and margin. If $r \geq 0.01$, LiteBIRD expects to achieve a $>5\sigma$ detection in the $\ell=$2-10 and $\ell=$11-200 ranges separately, providing crucial insight into the early Universe. We describe LiteBIRD's scientific objectives, the application of systems engineering to mission requirements, the anticipated scientific impact, and the operations and scanning strategies vital to minimizing systematic effects. We will also highlight LiteBIRD's synergies with concurrent CMB projects., Comment: 23 pages, 9 figures, 1 table, SPIE Astronomical Telescopes + Instrumentation 2024

We derive new constraints on the CIB monopole spectrum from a set of reprocessed COBE-DIRBE sky maps that have lower instrumental and astrophysical contamination than the legacy DIRBE maps. These maps have been generated through a global Bayesian analysis framework that simultaneously fits cosmological, astrophysical, and instrumental parameters, as described in a series of papers collectively referred to as Cosmoglobe DR2. We have applied this method to the DIRBE TODs, complemented by selected HFI and FIRAS sky maps to break key astrophysical degeneracies, as well as WISE and Gaia compact object catalogs. In this paper, we focus on the CIB monopole constraints that result from this work. We report positive detections of an isotropic signal in six out of the ten DIRBE bands (1.25, 2.2, 3.5, 100, 140, and 240 $\mathrm{\mu m}$). For the 2.2 $\mu$m channel, we find an amplitude of $10.2\pm1.2\,\mathrm{nW\,m^{-2}\,sr^{-1}}$, 74 % lower than that reported from the official DIRBE maps. For the 240 $\mu\mathrm{m}$ channel, we find $6\pm3\mathrm{nW\,m^{-2}\,sr^{-1}}$, 56 % lower than the official DIRBE release. We interpret these lower values as resulting from improved zodiacal light and Galactic foreground modeling. For the bands between 4.9 and 60 $\mu\mathrm{m}$, the presence of excess radiation in solar-centric coordinates precludes the definition of lower limits. However, the analysis still provides well-defined upper limits. For the 12 $\mu\mathrm{m}$ channel, we find an upper 95 % confidence limit of 55 $\mathrm{nW\,m^{-2}\,sr^{-1}}$, more than a factor of eight lower than the corresponding legacy result of 468 $\mathrm{nW\,m^{-2}\,sr^{-1}}$. The results presented in this paper redefine the state-of-the-art CIB monopole constraints from COBE-DIRBE, and provide a real-world illustration of the power of global end-to-end analysis of multiple complementary datasets., Comment: 17 pages, 9 figures, submitted to A&A. Comments welcome! Data products can be found at https://www.cosmoglobe.uio.no/products/cosmoglobe-dr2.html

The radiative branching ratio of the Hoyle state is crucial to estimate the triple-$\alpha$ reaction rate in stellar environments at medium temperatures. Knowledge of the $\gamma$-decay channel is critical as this is the dominant radiative decay channel for the Hoyle state. A recent study by Kib\'edi et al. [Phys. Rev. Lett. 125, 182701 (2020)] has challenged our understanding of this astrophysically significant branching ratio and its constraints. The objective of this work was to perform a new measurement of the $\gamma$-decay branching ratio of the Hoyle state to deduce the radiative branching ratio of the Hoyle state. An additional objective was to independently verify aspects of the aforementioned measurement conducted by Kib\'edi et al. For the main experiment of this work, the Hoyle state was populated by the $^{12}\textrm{C}(p,p')$ reaction at 10.8 MeV at the Oslo Cyclotron Laboratory. The $\gamma$-decay branching ratio was deduced through triple-coincidence events, each consisting of a proton ejectile corresponding to the Hoyle state, and the subsequent $\gamma$-ray cascade. In the main experiment of this work, a $\gamma$-decay branching ratio of the Hoyle state of $\Gamma_{\gamma}/\Gamma=4.0(4)\times 10^{-4}$ was determined, yielding a corresponding radiative branching ratio of $\Gamma_{\textrm{rad}}/\Gamma=4.1(4) \times 10^{-4}$, which is in agreement with several recent studies, as well as the previously adopted ENSDF average of $\Gamma_{\textrm{rad}}/\Gamma=4.16(11)\times 10^{-4}$. Aspects of the analysis performed by Kib\'edi et al. were verified in this work and the source of discrepancy between the results of this work and that of Kib\'edi et al. could not be determined. Further independent and innovative studies for the radiative width of the Hoyle state will substantiate whether the discrepant result by Kib\'edi et al. should be excluded from future evaluations.

Narrowband galaxy surveys have recently gained interest as a promising method to achieve the necessary accuracy on the photometric redshift estimate of individual galaxies for stage-IV cosmological surveys. One key advantage is the ability to provide higher spectral resolution information about galaxies that should allow a more accurate and precise estimation of galaxy stellar population properties. However, the impact of adding narrow-band photometry on the stellar population properties estimate is largely unexplored. The scope of this work is two-fold: on one side, leveraging the predictive power of broad-band and narrow-band data to infer galaxy physical properties such as stellar masses, ages, star formation rates and metallicities. On the other hand, evaluating the improvement of performance in estimating galaxy properties when we use narrow-band data instead of broad-band. In this work we measure the stellar population properties of a sample of galaxies in the COSMOS field for which both narrowband and broadband data are available. In particular, we employ narrowband data from PAUS and broad-band data from CFHTLS. We use two different spectral energy distribution fitting codes to measure galaxy properties, namely CIGALE and Prospector. We find that the increased spectral resolution of narrow-band photometry does not yield a substantial improvement on constraining galaxy properties using spectral energy distribution fitting. Still we find that we obtain a more diverse distribution of metallicities and dust optical depths with cigale when employing the narrowband data. The effect is not as prominent as expected, which we relate this to the low narrowband SNR of a majority of the galaxies, the respective drawbacks of both codes as well as the coverage only in the optical regime. The measured properties are afterwards compared to the COSMOS2020 catalogue, showing good agreement., Comment: Published in A&A. 18 pages, 16 figures, 2 tables

The Data Acquisition System (DAQ) for the LUX-ZEPLIN (LZ) dark matter detector is described. The signals from 745 PMTs, distributed across three subsystems, are sampled with 100-MHz 32-channel digitizers (DDC-32s). A basic waveform analysis is carried out on the on-board Field Programmable Gate Arrays (FPGAs) to extract information about the observed scintillation and electroluminescence signals. This information is used to determine if the digitized waveforms should be preserved for offline analysis. The system is designed around the Kintex-7 FPGA. In addition to digitizing the PMT signals and providing basic event selection in real time, the flexibility provided by the use of FPGAs allows us to monitor the performance of the detector and the DAQ in parallel to normal data acquisition. The hardware and software/firmware of this FPGA-based Architecture for Data acquisition and Realtime monitoring (FADR) are discussed and performance measurements are described., Comment: 18 pages, 24 figures

We revisit the impact of finite time responses of bolometric detectors used for deep observations of the cosmic microwave background (CMB). Until now, bolometer transfer functions have been accounted for through a two-step procedure by first deconvolving an estimate of their Fourier-space representation from the raw time-ordered data (TOD), and then averaging the deconvolved TOD into pixelized maps. However, for many experiments, including the Planck High Frequency Instrument (HFI), it is necessary to apply an additional low-pass filter to avoid an excessive noise boost, which leads to an asymmetric effective beam. In this paper we demonstrate that this effect can be avoided if the transfer function deconvolution and pixelization operations are performed simultaneously through integrated maximum likelihood mapmaking. The resulting algorithm is structurally identical to the artDeco algorithm introduced by Keih\"anen & Reinecke (2012) for beam deconvolution. We illustrate the relevance of this method with simulated Planck HFI 143 GHz data, and find that the resulting effective beam is both more symmetric than with the two-step procedure, resulting in a sky-averaged ellipticity that is 64 % lower, and an effective beam full-width-at-half-maximum (FWHM) that is 2.3 % smaller. Similar improvements are expected for any other bolometer-based CMB experiments with long time constants., Comment: 7 pages, 11 figures

Anomalous Microwave Emission (AME) is a major component of Galactic emission in the frequency band 10 to 60 GHz and is commonly modelled as rapidly rotating spinning dust grains. The photodissociation region (PDR) at the boundary of the $\lambda$-Orionis Hii region has been identified by several recent analyses as one of the brightest spinning dust emitting sources in the sky. We investigate the Barnard 30 dark cloud, a dark cloud embedded within the $\lambda$-Orionis PDR. We use total-power observations of Barnard 30 from the CO Mapping Array Project (COMAP) pathfinder instrument at 26 to 34GHz with a resolution of 4.5 arcminutes alongside existing data from Planck, WISE, IRAS, ACT, and the 1.447GHz GALFACTS survey. We use aperture photometry and template fitting to measure the spectral energy distribution of Barnard 30. We find that the spinning dust is the dominant emission component in the 26 to 34GHz range at the $7 \sigma$ level ($S_{30GHz} = 2.85\pm0.43$Jy). We find no evidence that polycyclic aromatic hydrocarbons are the preferred carrier for the spinning dust emission, suggesting that the spinning dust carriers are due to a mixed population of very small grains. Finally, we find evidence for variations in spinning dust emissivity and peak frequency within Barnard 30, and that these variations are possibly driven by changes in dust grain population and the total radiation field. Confirming the origin of the variations in the spinning dust spectrum will require both future COMAP observations at 15GHz combined with spectroscopic mid-infrared data of Barnard 30., Comment: 21 pages, 13 figures. Submitted to MNRAS

LUX-ZEPLIN (LZ) is a tonne-scale experiment searching for direct dark matter interactions and other rare events. It is located at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. The core of the LZ detector is a dual-phase xenon time projection chamber (TPC), designed with the primary goal of detecting Weakly Interacting Massive Particles (WIMPs) via their induced low energy nuclear recoils. Surrounding the TPC, two veto detectors immersed in an ultra-pure water tank enable reducing background events to enhance the discovery potential. Intricate calibration systems are purposely designed to precisely understand the responses of these three detector volumes to various types of particle interactions and to demonstrate LZ's ability to discriminate between signals and backgrounds. In this paper, we present a comprehensive discussion of the key features, requirements, and performance of the LZ calibration systems, which play a crucial role in enabling LZ's WIMP-search and its broad science program. The thorough description of these calibration systems, with an emphasis on their novel aspects, is valuable for future calibration efforts in direct dark matter and other rare-event search experiments.

The first science run of the LUX-ZEPLIN (LZ) experiment, a dual-phase xenon time project chamber operating in the Sanford Underground Research Facility in South Dakota, USA, has reported leading limits on spin-independent WIMP-nucleon interactions and interactions described from a non-relativistic effective field theory (NREFT). Using the same 5.5~t fiducial mass and 60 live days of exposure we report on the results of a relativistic extension to the NREFT. We present constraints on couplings from covariant interactions arising from the coupling of vector, axial currents, and electric dipole moments of the nucleon to the magnetic and electric dipole moments of the WIMP which cannot be described by recasting previous results described by an NREFT. Using a profile-likelihood ratio analysis, in an energy region between 0~keV$_\text{nr}$ to 270~keV$_\text{nr}$, we report 90% confidence level exclusion limits on the coupling strength of five interactions in both the isoscalar and isovector bases., Comment: 7 pages, 4 figures

Corrosion, a naturally occurring process leading to the deterioration of metallic materials, demands diligent detection for quality control and the preservation of metal-based objects, especially within industrial contexts. Traditional techniques for corrosion identification, including ultrasonic testing, radio-graphic testing, and magnetic flux leakage, necessitate the deployment of expensive and bulky equipment on-site for effective data acquisition. An unexplored alternative involves employing lightweight, conventional camera systems, and state-of-the-art computer vision methods for its identification. In this work, we propose a complete system for semi-automated corrosion identification and mapping in industrial environments. We leverage recent advances in LiDAR-based methods for localization and mapping, with vision-based semantic segmentation deep learning techniques, in order to build semantic-geometric maps of industrial environments. Unlike previous corrosion identification systems available in the literature, our designed multi-modal system is low-cost, portable, semi-autonomous and allows collecting large datasets by untrained personnel. A set of experiments in an indoor laboratory environment, demonstrate quantitatively the high accuracy of the employed LiDAR based 3D mapping and localization system, with less then $0.05m$ and 0.02m average absolute and relative pose errors. Also, our data-driven semantic segmentation model, achieves around 70\% precision when trained with our pixel-wise manually annotated dataset.

We present the first degree-scale tomography map of the dusty magnetized interstellar medium (ISM) from stellar polarimetry and distance measurements. We used the RoboPol polarimeter at Skinakas Observatory to conduct a survey of starlight polarization in a region of the sky of 4 square degrees. We propose a Bayesian method to decompose the stellar-polarization source field along the distance to invert the 3D volume occupied by the observed stars. We used it to obtain the first 3D map of the dusty magnetized ISM. Specifically, we produced a tomography map of the orientation of the plane-of-sky (POS) component of the magnetic field threading the diffuse, dusty regions responsible for the stellar polarization. For the targeted region centered on Galactic coordinates $(l,b) \approx (103.3^\circ, 22.3^\circ)$, we identified several ISM clouds. Most of the lines of sight intersect more than one cloud. A very nearby component was detected in the foreground of a dominant component from which most of the polarization signal comes. Farther clouds, with a distance of up to 2~kpc, were similarly detected. Some of them likely correspond to intermediate-velocity clouds seen in HI spectra in this region of the sky. We found that the orientation of the POS component of the magnetic field changes along distance for most of the lines of sight. Our study demonstrates that starlight polarization data coupled to distance measures have the power to reveal the great complexity of the dusty magnetized ISM in 3D and, in particular, to provide local measurements of the POS component of the magnetic field. This demonstrates that the inversion of large data volumes, as expected from the PASIPHAE survey, will provide the necessary means to move forward in the modeling of the Galactic magnetic field and of the dusty magnetized ISM as a contaminant in observations of the cosmic microwave background polarization., Comment: Preprint of the accepted version at Astronomy & Astrophysics. The 3D map obtained in this paper can be visualized online at https://pasiphae.science/visualization and a video featuring it is accessible at https://www.youtube.com/watch?v=dB_6J1zhmPI

We present a novel implementation of the complete active space self-consistent field (CASSCF) method that makes use of the many-body expanded full configuration interaction (MBE-FCI) method to incrementally approximate electronic structures within large active spaces. On the basis of a hybrid first-order algorithm employing both Super-CI and quasi-Newton strategies for the optimization of molecular orbitals, we demonstrate both computational efficacy and high accuracy of the resulting MBE-CASSCF method. We assess the performance of our implementation on a set of established numerical tests before applying MBE-CASSCF in the investigation of the triplet-quintet spin gap of iron(II) porphyrin with active spaces as large as 50 electrons in 50 orbitals., Comment: 27+8 pages, 6 figures. SI as an ancillary file

We present detailed forecasts for the constraints on primordial magnetic fields (PMFs) that will be obtained with the LiteBIRD satellite. The constraints are driven by the effects of PMFs on the CMB anisotropies: the gravitational effects of magnetically-induced perturbations; the effects on the thermal and ionization history of the Universe; the Faraday rotation imprint on the CMB polarization; and the non-Gaussianities induced in polarization anisotropies. LiteBIRD represents a sensitive probe for PMFs and by exploiting all the physical effects, it will be able to improve the current limit coming from Planck. In particular, thanks to its accurate $B$-mode polarization measurement, LiteBIRD will improve the constraints on infrared configurations for the gravitational effect, giving $B_{\rm 1\,Mpc}^{n_{\rm B} =-2.9} < 0.8$ nG at 95% C.L., potentially opening the possibility to detect nanogauss fields with high significance. We also observe a significant improvement in the limits when marginalized over the spectral index, $B_{1\,{\rm Mpc}}^{\rm marg}< 2.2$ nG at 95% C.L. From the thermal history effect, which relies mainly on $E$-mode polarization data, we obtain a significant improvement for all PMF configurations, with the marginalized case, $\sqrt{\langle B^2\rangle}^{\rm marg}<0.50$ nG at 95% C.L. Faraday rotation constraints will take advantage of the wide frequency coverage of LiteBIRD and the high sensitivity in $B$ modes, improving the limits by orders of magnitude with respect to current results, $B_{1\,{\rm Mpc}}^{n_{\rm B} =-2.9} < 3.2$ nG at 95% C.L. Finally, non-Gaussianities of the $B$-mode polarization can probe PMFs at the level of 1 nG, again significantly improving the current bounds from Planck. Altogether our forecasts represent a broad collection of complementary probes, providing conservative limits on PMF characteristics that will be achieved with LiteBIRD., Comment: 51 pages, 24 figures, abstract shortened