We present the development of two complementary amplifier architectures for axion haloscope experiments, based on two types of Josephson Parametric Amplifiers (JPAs). The first employs a multi-chip module of flux-driven JPAs in a parallel–series configuration, enabling near quantum-limited amplification over an extended tunable range of between 1.2 and 1.5 GHz. The second design features a lumped-element JPA, offering continuous tunability across a wide frequency range from 2.4 to 4 GHz. Both approaches demonstrate near-quantum-limited noise performance and are compatible with operation in cryogenic environments. These amplifiers significantly enhance the sensitivity and frequency coverage of axion search experiments, and also provide new opportunities for broadband quantum sensing applications.
Superconducting Transition Edge Sensors (TESs) are a promising technology for fundamental physics applications due to their low dark count rates, excellent energy resolution, and high detection efficiency. On the DESY campus, we have been developing a program to characterize cryogenic quantum sensors for fundamental physics applications, particularly focused on TESs. We currently have two fully equipped dilution refrigerators that enable simultaneous TES characterization and fundamental physics searches. In this paper, we summarize the current status of our TES characterization, including recent calibration efforts and efficiency measurements, as well as simulations to better understand TES behavior and its backgrounds. Additionally, we summarize some physics applications that we are already exploring or planning to explore. We will give preliminary projections on a direct dark matter search with our TES, where exploiting low-threshold electron scattering in superconducting materials allows us to search for sub-MeV-scale dark matter. We are also working toward performing a measurement of the even-number photon distribution (beyond one pair) of a quantum-squeezed light source. Finally, if it proves to meet the requirements, our TES detector may be used as a second, independent detection system to search for an axion signal at the ALPS II experiment.
Of the three major axion search experimental strategies, light-shining-through-wall experiments, haloscopes, and helioscopes, this paper focuses on the latter. IAXO, the International AXion Observatory, will be a next-generation helioscope following in the footsteps of previous experiments like SUMICO and CAST. Helioscopes aim to detect axions produced in the Sun, utilizing a magnetic field to couple them to X-ray photons. BabyIAXO represents a near-term step toward IAXO, designed to test custom components while delivering competitive results in axion searches. The experimental components are currently under development and construction. Further research into the applications of BabyIAXO beyond baseline axion searches is being conducted.
Consider the dynamical system x ˙ = f ( x ) , where x ∈ R n is the state vector, x ˙ is the time or spatial derivative, and f is the system model. We wish to identify unknown f from its time-series or spatial data. For this, we propose a Bayesian framework based on the maximum a posteriori (MAP) point estimate, to give a generalized Tikhonov regularization method with the residual and regularization terms identified, respectively, with the negative logarithms of the likelihood and prior distributions. As well as estimates of the model coefficients, the Bayesian interpretation provides access to the full Bayesian apparatus, including the ranking of models, the quantification of model uncertainties, and the estimation of unknown (nuisance) hyperparameters. For multivariate Gaussian likelihood and prior distributions, the Bayesian formulation gives a Gaussian posterior distribution, in which the numerator contains a Mahalanobis distance or “Gaussian norm”. In this study, two Bayesian algorithms for the estimation of hyperparameters—the joint maximum a posteriori (JMAP) and variational Bayesian approximation (VBA)—are compared to the popular SINDy, LASSO, and ridge regression algorithms for the analysis of several dynamical systems with additive noise. We consider two dynamical systems, the Lorenz convection system and the Shil’nikov cubic system, with four choices of noise model: symmetric Gaussian or Laplace noise and skewed Rayleigh or Erlang noise, with different magnitudes. The posterior Gaussian norm is found to provide a robust metric for quantitative model selection—with quantification of the model uncertainties—across all dynamical systems and noise models examined.
The modelling of plasma–wall interactions (PWIs) depends on distributions describing the angle and energy distribution of particles scattered at the first wall of fusion devices. Most PWI codes rely on extensive tables based on data from reflection simulations, employing a Monte Carlo method. At first glance, the uncertainty distribution of the data should be assumed Gaussian. However, in order to obtain the resulting particle distribution, the reflected ions are counted within angle sections of the upper hemisphere, which hints to a Poisson uncertainty distribution. In this paper, we let Bayesian model comparison decide which uncertainty model should be taken.
The static topology of surface characteristics and active sites in catalysis overlooks a crucial element: the dynamic processes of optimal pattern formation over time and the creation of intermediate structures that enhance reactions. Nature’s principle of coupling reaction and motion in catalytic processes by enzymes or higher organisms offers a new perspective. This work explores a novel theoretical approach by adding the time dimension to optimise topological variations using the Maximum Entropy Production (MEP) assumption. This approach recognises that the catalyst surface is not an unchanging energy landscape but can change dynamically. The time-dependent transport problem of molecules is here interpreted by a non-equilibrium model used for modelling and predicting dynamic pattern formation in excitable media, a class of active matter requiring an activation threshold. We present a nonlocal reaction–cross-diffusion (RXD) formulation of catalytic reactions that can capture the catalyst’s interaction with the target molecule in space and time. The approach provides a theoretical basis for future deep learning models and multiphysics upscaling of catalysts and their support structures across multiphysics fields. The particular advantage of the RXD approach is that it allows each scale to investigate dynamic pattern-forming processes using linear and nonlinear stability analysis, thus establishing a rule base for developing new catalysts.
Groundwater vulnerability assessment (GVA) is critical for understanding contaminant migration into groundwater systems, yet conventional methods often overlook its probabilistic nature. Bayesian inference offers a robust framework using Bayes’ rule to enhance decision-making through posterior probability calculations. This study introduces inverse Bayesian methods for GVA using spatial-series data, focusing on nitrate concentrations in groundwater as an indicator of groundwater vulnerability in agricultural catchments. Using the joint maximum a-posteriori (JMAP) and variational Bayesian approximation (VBA) algorithms, the advantages of the Bayesian framework over traditional index-based methods are demonstrated for GVA of the Burdekin Basin, Queensland, Australia. This provides an evidence-based methodology for GVA which enables model ranking, parameter estimation, and uncertainty quantification.
Building appropriate models is crucial for imaging tasks in many fields but often challenging due to the richness of the systems. In radio astronomy, for example, wide-field observations can contain various and superposed structures that require different descriptions, such as filaments, point sources or compact objects. This work presents an automatic pipeline that iteratively adapts probabilistic models for such complex systems in order to improve the reconstructed images. It uses the Bayesian imaging library NIFTy, which is formulated in the language of information field theory. Starting with a preliminary reconstruction using a simple and flexible model, the pipeline employs deep learning and clustering methods to identify and separate different objects. In a further step, these objects are described by adding new building blocks to the model, allowing for a component separation in the next reconstruction step. This procedure can be repeated several times for refinement to iteratively improve the overall reconstruction. In addition, the individual components can be modeled at different resolutions allowing us to focus on important parts of the emission field without getting computationally too expensive.
Magnetic confinement nuclear fusion offers a promising solution to the world’s growing energy demands. The DEMO reactor presented here aims to bridge the gap between laboratory fusion experiments and practical electricity generation, posing unique challenges for magnetic plasma diagnostics due to limited space for diagnostic equipment. This study employs Bayesian inference and Gaussian process modeling to integrate data from pick-up coils, flux loops, and saddle coils, enabling a qualitative estimation of the plasma current density distribution relying on only external magnetic measurements. The methodology successfully infers total plasma current, plasma centroid position, and six plasma–wall gap positions, while adhering to DEMO’s stringent accuracy standards. Additionally, the interchangeability between normal pick-up coils and saddle coils was assessed, revealing a clear preference for saddle coils. Initial steps were taken to utilize Bayesian experimental design for optimizing the orientation (normal or tangential) of pick-up coils within DEMO’s design constraints to improve the diagnostic setup’s inference precision. Our approach indicates the feasibility of Bayesian integrated data analysis in achieving precise and accurate probability distributions of plasma parameter crucial for the successful operation of DEMO.
Two outstanding problems of particle physics and cosmology, namely the strong-CP problem and the nature of dark matter, can be solved with the discovery of a single new particle, the axion. The modular high magnetic field and flux hybrid magnet platform of LNCMI-Grenoble, which was recently put in operation up to 42 T, offers unique opportunities for axion/axion-like particle search using Sikivie-type haloscopes. In this paper, the focus will be on the 350–600 MHz frequency range corresponding to the 1–3 μeV axion mass range requiring a large-bore RF-cavity. It will be built by DMAG and integrated within the large-bore superconducting hybrid magnet outsert, providing a central magnetic field up to 9 T in 812 mm warm bore diameter. The progress achieved by Néel Institute in the design of the complex cryostat with its double dilution refrigerators to cooldown below 50 mK the ultra-light Cu RF-cavity of 650 mm inner diameter and the first stage of the RF measurement chain are presented. Perspectives for the targeted sensitivity, assuming less than 2-year integration time, are recalled.
In this article, we present an extension of the Full Bayesian Significance Test (FBST) for nonparametric settings, termed NP-FBST, which is constructed using the limit of finite dimension histograms. The test statistics for NP-FBST are based on a plug-in estimate of the cross-entropy between the null hypothesis and a histogram. This method shares similarities with Kullback–Leibler and entropy-based goodness-of-fit tests, but it can be applied to a broader range of hypotheses and is generally less computationally intensive. We demonstrate that when the number of histogram bins increases slowly with the sample size, the NP-FBST is consistent for Lipschitz continuous data-generating densities. Additionally, we propose an algorithm to optimize the NP-FBST. Through simulations, we compare the performance of the NP-FBST to traditional methods for testing uniformity. Our results indicate that the NP-FBST is competitive in terms of power, even surpassing the most powerful likelihood-ratio-based procedures for very small sample sizes.
This letter aims to investigate thermodynamic processes in small systems in the Onsager region by showing that fundamental quantities such as total entropy production can be discretized on the mesoscopic scale. Even thermodynamic variables can conjugate to thermodynamic forces, and thus, Glansdorff–Prigogine’s dissipative variable may be discretized. The canonical commutation rules (CCRs) valid at the mesoscopic scale are postulated, and the measurement process consists of determining the eigenvalues of the operators associated with the thermodynamic quantities. The nature of the quantized quantity β , entering the CCRs, is investigated by a heuristic model for nano-gas and analyzed through the tools of classical statistical physics. We conclude that according to our model, the constant β does not appear to be a new fundamental constant but corresponds to the minimum value.
The evolution of space technology in recent years, fueled by advancements in computing such as Artificial Intelligence (AI) and machine learning (ML), has profoundly transformed our capacity to explore the cosmos. Missions like the James Webb Space Telescope (JWST) have made information about distant objects more easily accessible, resulting in extensive amounts of valuable data. As part of this work-in-progress study, we are working to create an atmospheric absorption spectrum prediction model for exoplanets. The eventual model will be based on both collected observational spectra and synthetic spectral data generated by the ROCKE-3D general circulation model (GCM) developed by the climate modeling program at NASA's Goddard Institute for Space Studies (GISS). In this initial study, spline curves are used to describe the bin heights of simulated atmospheric absorption spectra as a function of one of the values of the planetary parameters. Bayesian Adaptive Exploration is then employed to identify areas of the planetary parameter space for which more data are needed to improve the model. The resulting system will be used as a forward model so that planetary parameters can be inferred given a planet's atmospheric absorption spectrum. This work is expected to contribute to a better understanding of exoplanetary properties and general exoplanet climates and habitability.
Axions and axion-like particles (ALPs) are well-motivated dark matter (DM) candidates that couple with photons in external magnetic fields. The parameter space around m a ∼ 50 meV remains largely unexplored by haloscope experiments. We present the first prototype of Weakly Interacting Sub-eV Particles (WISP) Searches on a Fiber Interferometer (WISPFI), a table-top, model-independent scheme based on resonant photon–axion conversion in a hollow-core photonic crystal fiber (HC-PCF) integrated into a Mach–Zehnder interferometer (MZI). Operating near a dark fringe with active phase-locking, combined with amplitude modulation, the interferometer converts axion-induced photon disappearance into a measurable signal. A 2 W, 1550 nm laser is coupled with a 1 m-long HC-PCF placed inside a ∼2 T permanent magnet array, probing a fixed axion mass of m a ≃ 49 meV with a projected sensitivity of g a γ γ ≳ 1.3× 10 − 9 GeV−1 for a measurement time of 30 days. Future upgrades, including pressure tuning of the effective refractive index and implementation of a Fabry–Pérot cavity, could extend the accessible mass range and improve sensitivity, establishing WISPFI as a scalable platform to explore previously inaccessible regions of the axion parameter space.
When the quantities of interest remain underdetermined a posteriori, we would like to draw inferences for at least one particular solution. Can we do so in a Bayesian way? What is a probability distribution over an underdetermined quantity? How do we get a posterior for one particular solution from a posterior for infinitely many underdetermined solutions? Guided by discrete invariant underdetermined ill-posed inverse problems, we find that a probability distribution over an underdetermined quantity is non-absolutely continuous, partially improper with respect to the initial reference measure but proper with respect to its restriction to its support. Thus, it is necessary and sufficient to choose the prior restricted reference measure to assign partially improper priors using e.g., the principle of maximum entropy and the posterior restricted reference measure to obtain the proper posterior for one particular solution. We can then work with underdetermined models like Hoeffding–Sobol expansions seamlessly, especially to effectively counter the curse of dimensionality within discrete nonparametric inverse problems. We show Singular Bayesian Inference (SBI) at work in an advanced Bayesian optimization application: dynamic pricing. Such a nice generalization of Bayesian–maxentropic inference could motivate many theoretical and practical developments.
This paper explores the application of the Value of Information, (VoI), based on the Claude Shannon/Ruslan Stratonovich framework within economic contexts. Unlike previous studies that examine circular settings and strategic interactions, we focus on a non-strategic linear setting. We employ standard economically motivated utility functions, including linear, quadratic, constant absolute risk aversion (CARA), and constant relative risk aversion (CRRA), across various priors of the stochastic environment, and analyse the resulting specific VoI forms. The curvature of these VoI functions play a decisive role in determining whether acquiring additional costly information enhances the efficiency of the decision making process. We also outline potential implications for broader decision-making frameworks.
Using Bayesian inference to calibrate constitutive model parameters has recently seen a rise in interest. The Markov chain Monte Carlo (MCMC) algorithm is one of the most commonly used methods to sample from the posterior. However, the choice of which MCMC algorithm to apply is typically pragmatic and based on considerations such as software availability and experience. We compare three commonly used MCMC algorithms: Metropolis-Hastings (MH), Affine Invariant Stretch Move (AISM) and No-U-Turn sampler (NUTS). For the comparison, we use the Kullback-Leibler (KL) divergence as a convergence criterion, which measures the statistical distance between the sampled and the ‘true’ posterior. We apply the Bayesian framework to a Newtonian squeeze flow problem, for which there exists an analytical model. Furthermore, we have collected experimental data using a tailored setup. The ground truth for the posterior is obtained by evaluating it on a uniform reference grid. We conclude that, for the same number of samples, the NUTS results in the lowest KL divergence, followed by the AISM sampler and last the MH sampler.
Bayesian analysis is particularly useful for inferring models and their parameters given data. This is a common task in metabolic modeling, where models of varying complexity are used to interpret data. Nested sampling is a class of probabilistic inference algorithms that are particularly effective for estimating evidence and sampling the parameter posterior probability distributions. However, the practicality of nested sampling for metabolic network inference has yet to be studied. In this technical report, we explore the amalgamation of nested sampling, specifically diffusive nested sampling, with reversible jump Markov chain Monte Carlo. We apply the algorithm to two synthetic problems from the field of metabolic flux analysis. We present run times and share insights into hyperparameter choices, providing a useful point of reference for future applications of nested sampling to metabolic flux problems.
Better understanding of the fitness and flexibility of microbial platform organisms is central to biotechnological process development. Live-cell experiments uncover the phenotypic heterogeneity of living cells, emerging even within isogenic cell populations. However, how this observed heterogeneity in growth relates to the variability of intracellular processes that drive cell growth and division is less understood. We here approach the question, how the observed phenotypic variability in single-cell growth rates links to metabolic processes, specifically intracellular reaction rates (fluxes). To approach this question, we employ the Maximum Entropy (MaxEnt) principle that allows us to bring together the phenotypic solution space, derived from metabolic network models, to single-cell growth rates observed in live-cell experiments. We apply the computational machinery to first-of-its-kind data of the microorganism Corynebacterium glutamicum, grown on different substrates under continuous medium supply. We compare the MaxEnt-based estimates of metabolic fluxes with estimates obtained by assuming that the average cell operates at its maximum growth rate, which is the current predominant practice in biotechnology.
