Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Covariance Matrix”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Measuring σ 8 using DESI Legacy Imaging Surveys Emission-Line galaxies and Planck CMB lensing, and the impact of dust on parameter inference

Measuring the growth of structure is a powerful probe for studying the dark sector, especially in light of the σ 8 tension between primary CMB anisotropy and low-redshift surveys. This paper provides a new measurement of the amplitude of the matter power spectrum, σ 8 , using galaxy-galaxy and galaxy-CMB lensing power spectra of Dark Energy Spectroscopic Instrument Legacy Imaging Surveys Emission-Line Galaxies and the Planck 2018 CMB lensing map. We create an ELG catalog composed of 24 million galaxies and with a purity of 85%, covering a redshift range 0 < z < 3, with z mean = 1.09. We implement several novel systematic corrections, such as jointly modeling the contribution of imaging systematics and photometric redshift uncertainties to the covariance matrix. We also study the impacts of various dust maps on cosmological parameter inference. We measure the cross-power spectra over f sky = 0.25 with a signal-to-background ratio of up to 30σ. We find that the choice of dust maps to account for imaging systematics in estimating the ELG overdensity field has a significant impact on the final estimated values of σ 8 and Ω M , with far-infrared emission-based dust maps preferring σ 8 to be as low as 0.702 ± 0.030, and stellar-reddening-based dust maps preferring as high as 0.719 ± 0.030. The highest preferred value is at ∼ 3 σ tension with the Planck primary anisotropy results. These findings indicate a need for tomographic analyses at high redshifts and joint modeling of systematics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

DESI DR1 Ly α 1D power spectrum: the Fast Fourier Transform estimator measurement

Here, we present the one-dimensional Lyman-α forest power spectrum measurement derived from the data release 1 (DR1) of the Dark Energy Spectroscopic Instrument (DESI). The measurement of the Lyman-α forest power spectrum along the line of sight from high-redshift quasar spectra provides information on the shape of the linear matter power spectrum, neutrino masses, and the properties of dark matter. In this work, we use a Fast Fourier Transform (FFT)-based estimator, which is validated on synthetic data in a companion paper. Compared to the FFT measurement performed on the DESI early data release, we improve the noise characterization with a cross-exposure estimator and test the robustness of our measurement using various data splits. We also refine the estimation of the uncertainties and now present an estimator for the covariance matrix of the measurement. Furthermore, we compare our results to previous high-resolution and eBOSS measurements. In another companion paper, we present the same DR1 measurement using the Quadratic Maximum Likelihood Estimator (QMLE). These two measurements are consistent with each other and constitute the most precise one-dimensional power spectrum measurement to date, while being in good agreement with results from the DESI early data release.

Lyman alpha forest↗
DOE OSTI2022

Solving quasi-free and quadratic Lindblad master equations for open fermionic and bosonic systems

The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), here we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k-point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Implementation of a practical Markov chain Monte Carlo sampling algorithm in PyBioNetFit

Abstract Summary Bayesian inference in biological modeling commonly relies on Markov chain Monte Carlo (MCMC) sampling of a multidimensional and non-Gaussian posterior distribution that is not analytically tractable. Here, we present the implementation of a practical MCMC method in the open-source software package PyBioNetFit (PyBNF), which is designed to support parameterization of mathematical models for biological systems. The new MCMC method, am, incorporates an adaptive move proposal distribution. For warm starts, sampling can be initiated at a specified location in parameter space and with a multivariate Gaussian proposal distribution defined initially by a specified covariance matrix. Multiple chains can be generated in parallel using a computer cluster. We demonstrate that am can be used to successfully solve real-world Bayesian inference problems, including forecasting of new Coronavirus Disease 2019 case detection with Bayesian quantification of forecast uncertainty. Availability and implementation PyBNF version 1.1.9, the first stable release with am, is available at PyPI and can be installed using the pip package-management system on platforms that have a working installation of Python 3. PyBNF relies on libRoadRunner and BioNetGen for simulations (e.g. numerical integration of ordinary differential equations defined in SBML or BNGL files) and Dask.Distributed for task scheduling on Linux computer clusters. The Python source code can be freely downloaded/cloned from GitHub and used and modified under terms of the BSD-3 license (https://github.com/lanl/pybnf). Online documentation covering installation/usage is available (https://pybnf.readthedocs.io/en/latest/). A tutorial video is available on YouTube (https://www.youtube.com/watch?v=2aRqpqFOiS4&t=63s). Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Pre-conditioned BFGS-based uncertainty quantification in elastic full-waveform inversion

SUMMARY Full-waveform inversion has become an essential technique for mapping geophysical subsurface structures. However, proper uncertainty quantification is often lacking in current applications. In theory, uncertainty quantification is related to the inverse Hessian (or the posterior covariance matrix). Even for common geophysical inverse problems its calculation is beyond the computational and storage capacities of the largest high-performance computing systems. In this study, we amend the Broyden–Fletcher–Goldfarb–Shanno (BFGS) algorithm to perform uncertainty quantification for large-scale applications. For seismic inverse problems, the limited-memory BFGS (L-BFGS) method prevails as the most efficient quasi-Newton method. We aim to augment it further to obtain an approximate inverse Hessian for uncertainty quantification in FWI. To facilitate retrieval of the inverse Hessian, we combine BFGS (essentially a full-history L-BFGS) with randomized singular value decomposition to determine a low-rank approximation of the inverse Hessian. Setting the rank number equal to the number of iterations makes this solution efficient and memory-affordable even for large-scale problems. Furthermore, based on the Gauss–Newton method, we formulate different initial, diagonal Hessian matrices as pre-conditioners for the inverse scheme and compare their performances in elastic FWI applications. We highlight our approach with the elastic Marmousi benchmark model, demonstrating the applicability of pre-conditioned BFGS for large-scale FWI and uncertainty quantification.

58 GEOSCIENCES↗
DOE OSTI2025

Evaluation of a generalized least squares algorithm for infrasound beamforming with coherent background noise

Infrasonic signals of interest can occur during periods with persistent, coherent, background noise, which may be natural or anthropogenic. For high signal-to-noise (SNR) ratio transient signals, an ‘overprinting’ of the coherent background may occur, and the signal may still be detected. However, this approach fails for low SNR signals of interest, which may be obscured by coherent noise. An infrasound beamforming method based on generalized least squares (GLS) is investigated for detecting transient signals of interest in the presence of coherent and incoherent background noise. This approach relies on an estimate of the noise covariance, captured in a covariance matrix, to effectively null contributions to the array response from noisy directions of arrival. Synthetic array data is used to investigate the performance of the GLS beamformer compared to the Bartlett beamformer when coherent and incoherent backgrounds are present. Additionally, the effects of array element number and relative strength of the interfering signal on the GLS estimates is investigated. GLS empirical area under the curve estimates suggest that the beamformer can recover coherent power for a signal of interest lower in amplitude than the coherent background, but this effectiveness degrades more quickly with SNR for a four element array compared to a six or eight element infrasound array. Finally, infrasound from the Forensic Surface Experiment, a bolide signal observed at IMS array I37NO, and a volcanic signal recorded at the Alaska Volcano Observatory array ADKI are used to evaluate GLS performance on recorded data. A ten minute window was used to capture the background noise, and the coherent background signal was nulled in all three examples.

58 GEOSCIENCES↗
DOE OSTI2020

The stellar mass Fundamental Plane: the virial relation and a very thin plane for slow rotators

ABSTRACT Early-type galaxies – slow and fast rotating ellipticals (E-SRs and E-FRs) and S0s/lenticulars – define a Fundamental Plane (FP) in the space of half-light radius Re, enclosed surface brightness Ie, and velocity dispersion σe. Since Ie and σe are distance-independent measurements, the thickness of the FP is often expressed in terms of the accuracy with which Ie and σe can be used to estimate sizes Re. We show that: (1) The thickness of the FP depends strongly on morphology. If the sample only includes E-SRs, then the observed scatter in Re is $\sim 16{{\ \rm per\ cent}}$, of which only $\sim 9{{\ \rm per\ cent}}$ is intrinsic. Removing galaxies with M* < 1011 M⊙ further reduces the observed scatter to $\sim 13{{\ \rm per\ cent}}$ ($\sim 4{{\ \rm per\ cent}}$ intrinsic). The observed scatter increases to $\sim 25{{\ \rm per\ cent}}$ usually quoted in the literature if E-FRs and S0s are added. If the FP is defined using the eigenvectors of the covariance matrix of the observables, then the E-SRs again define an exceptionally thin FP, with intrinsic scatter of only 5 per cent orthogonal to the plane. (2) The structure within the FP is most easily understood as arising from the fact that Ie and σe are nearly independent, whereas the Re−Ie and Re−σe correlations are nearly equal and opposite. (3) If the coefficients of the FP differ from those associated with the virial theorem the plane is said to be ‘tilted’. If we multiply Ie by the global stellar mass-to-light ratio M*/L and we account for non-homology across the population by using Sérsic photometry, then the resulting stellar mass FP is less tilted. Accounting self-consistently for M*/L gradients will change the tilt. The tilt we currently see suggests that the efficiency of turning baryons into stars increases and/or the dark matter fraction decreases as stellar surface brightness increases.

Bernardi, M.↗
DOE OSTI2020

Dark Energy Survey Year 3 results: cosmology with moments of weak lensing mass maps – validation on simulations

ABSTRACT We present a simulated cosmology analysis using the second and third moments of the weak lensing mass (convergence) maps. The second moment, or variances, of the convergence as a function of smoothing scale contains information similar to standard shear two-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. The analysis is geared towards the third year (Y3) data from the Dark Energy Survey (DES), but the methodology can be applied to other weak lensing data sets. We present the formalism for obtaining the convergence maps from the measured shear and for obtaining the second and third moments of these maps given partial sky coverage. We estimate the covariance matrix from a large suite of numerical simulations. We test our pipeline through a simulated likelihood analyses varying 5 cosmological parameters and 10 nuisance parameters and identify the scales where systematic or modelling uncertainties are not expected to affect the cosmological analysis. Our simulated likelihood analysis shows that the combination of second and third moments provides a 1.5 per cent constraint on S8 ≡ σ8(Ωm/0.3)0.5 for DES Year 3 data. This is 20 per cent better than an analysis using a simulated DES Y3 shear two-point statistics, owing to the non-Gaussian information captured by the inclusion of higher order statistics. This paper validates our methodology for constraining cosmology with DES Year 3 data, which will be presented in a subsequent paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Towards a self-consistent analysis of the anisotropic galaxy two- and three-point correlation functions on large scales: application to mock galaxy catalogues

We establish a practical method for the joint analysis of anisotropic galaxy two- and three-point correlation functions (2PCF and 3PCF, respectively) on the basis of the decomposition formalism of the 3PCF using tripolar spherical harmonics. We perform such an analysis with MultiDark-Patchy mock catalogues to demonstrate and understand the benefit of the anisotropic 3PCF. We focus on scales above $80\, h^{-1}\, {\rm Mpc}$, and use information from the shape and the baryon acoustic oscillation (BAO) signals of the 2PCF and 3PCF. We also apply density field reconstruction to increase the signal-to-noise ratio of BAO in the 2PCF measurement, but not in the 3PCF measurement. In particular, we study in detail the constraints on the angular diameter distance and the Hubble parameter. We build a model of the bispectrum or 3PCF that includes the non-linear damping of the BAO signal in redshift space. We carefully account for various uncertainties in our analysis including theoretical models of the 3PCF, window function corrections, biases in estimated parameters from the fiducial values, the number of mock realizations to estimate the covariance matrix, and bin size. The joint analysis of the 2PCF and 3PCF monopole and quadrupole components shows a $30{{\ \rm per\ cent}}$ and $20{{\ \rm per\ cent}}$ improvement in Hubble parameter constraints before and after reconstruction of the 2PCF measurements, respectively, compared to the 2PCF analysis alone. Furthermore, this study clearly shows that the anisotropic 3PCF increases cosmological information from galaxy surveys and encourages further development of the modelling of the 3PCF on smaller scales than we consider.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Cosmic shear cosmology beyond two-point statistics: a combined peak count and correlation function analysis of DES-Y1

ABSTRACT We constrain cosmological parameters from a joint cosmic shear analysis of peak-counts and the two-point shear correlation functions, as measured from the Dark Energy Survey (DES-Y1). We find the structure growth parameter $S_8\equiv \sigma _8\sqrt{\Omega _{\rm m}/0.3} = 0.766^{+0.033}_{-0.038}$ which, at 4.8 per cent precision, provides one of the tightest constraints on S8 from the DES-Y1 weak lensing data. In our simulation-based method we determine the expected DES-Y1 peak-count signal for a range of cosmologies sampled in four w cold dark matter parameters (Ωm, σ8, h, w0). We also determine the joint covariance matrix with over 1000 realizations at our fiducial cosmology. With mock DES-Y1 data we calibrate the impact of photometric redshift and shear calibration uncertainty on the peak-count, marginalizing over these uncertainties in our cosmological analysis. Using dedicated training samples we show that our measurements are unaffected by mass resolution limits in the simulation, and that our constraints are robust against uncertainty in the effect of baryon feedback. Accurate modelling for the impact of intrinsic alignments on the tomographic peak-count remains a challenge, currently limiting our exploitation of cross-correlated peak counts between high and low redshift bins. We demonstrate that once calibrated, a fully tomographic joint peak-count and correlation functions analysis has the potential to reach a 3 per cent precision on S8 for DES-Y1. Our methodology can be adopted to model any statistic that is sensitive to the non-Gaussian information encoded in the shear field. In order to accelerate the development of these beyond-two-point cosmic shear studies, our simulations are made available to the community upon request.

Harnois-Déraps, Joachim↗
DOE OSTI2021

Galaxy clustering in harmonic space from the dark energy survey year 1 data: compatibility with real-space results

ABSTRACT We perform an analysis in harmonic space of the Dark Energy Survey Year 1 Data (DES-Y1) galaxy clustering photometric data using products obtained for the real-space analysis. We test our pipeline with a suite of lognormal simulations, which are used to validate scale cuts in harmonic space as well as to provide a covariance matrix that takes into account the DES-Y1 mask. We then apply this pipeline to DES-Y1 data taking into account survey property maps derived for the real-space analysis. We compare with real-space DES-Y1 results obtained from a similar pipeline. We show that the harmonic space analysis we develop yields results that are compatible with the real-space analysis for the bias parameters. This verification paves the way to performing a harmonic space analysis for the upcoming DES-Y3 data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Advances in constraining intrinsic alignment models with hydrodynamic simulations

We use galaxies from the illustristng, massiveblack-ii, and illustris-1 hydrodynamic simulations to investigate the behaviour of large scale galaxy intrinsic alignments. Our analysis spans four redshift slices over the approximate range of contemporary lensing surveys z = 0–1. We construct comparable weighted samples from the three simulations, which we then analyse using an alignment model that includes both linear and quadratic alignment contributions. Our data vector includes galaxy–galaxy, galaxy–shape, and shape–shape projected correlations, with the joint covariance matrix estimated analytically. In all of the simulations, we report non-zero IAs at the level of several σ. For a fixed lower mass threshold, we find a relatively strong redshift dependence in all three simulations, with the linear IA amplitude increasing by a factor of ~2 between redshifts z = 0 and z = 1. We report no significant evidence for non-zero values of the tidal torquing amplitude, A 2 , in TNG, above statistical uncertainties, although MBII favours a moderately negative A 2 ~ –2. Examining the properties of the TATT model as a function of colour, luminosity and galaxy type (satellite or central), our findings are consistent with the most recent measurements on real data. We also outline a novel method for constraining the TATT model parameters directly from the pixelized tidal field, alongside a proof-of-concept exercise using TNG. This technique is shown to be promising, although comparison with previous results obtained via other methods is non-trivial.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The first Hubble diagram and cosmological constraints using superluminous supernovae

In this work, we present the first Hubble diagram of superluminous supernovae (SLSNe) out to a redshift of two, together with constraints on the matter density, Ω M , and the dark energy equation-of-state parameter, w(≡p/ρ). We build a sample of 20 cosmologically useful SLSNe I based on light curve and spectroscopy quality cuts. We confirm the robustness of the peak–decline SLSN I standardization relation with a larger data set and improved fitting techniques than previous works. We then solve the SLSN model based on the above standardization via minimization of the χ2 computed from a covariance matrix that includes statistical and systematic uncertainties. For a spatially flat Λ cold dark matter (ΛCDM) cosmological model, we find $\Omega _{\rm M}=0.38^{+0.24}_{-0.19}$, with an rms of 0.27 mag for the residuals of the distance moduli. For a w 0 w a CDM cosmological model, the addition of SLSNe I to a ‘baseline’ measurement consisting of Planck temperature together with Type Ia supernovae, results in a small improvement in the constraints of w 0 and wa of 4 per cent. We present simulations of future surveys with 868 and 492 SLSNe I (depending on the configuration used) and show that such a sample can deliver cosmological constraints in a flat ΛCDM model with the same precision (considering only statistical uncertainties) as current surveys that use Type Ia supernovae, while providing a factor of 2–3 improvement in the precision of the constraints on the time variation of dark energy, w 0 and w a . This paper represents the proof of concept for superluminous supernova cosmology, and demonstrates they can provide an independent test of cosmology in the high-redshift (z > 1) universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The impact and mitigation of broad-absorption-line quasars in Lyman α forest correlations

ABSTRACT Correlations in and with the flux transmission of the Lyman α (Ly α) forest in the spectra of high-redshift quasars are powerful cosmological tools, yet these measurements can be compromised if the intrinsic quasar continuum is significantly uncertain. One particularly problematic case is broad-absorption-line (BAL) quasars, which exhibit blueshifted absorption associated with many spectral features that are consistent with outflows of up to ∼0.1c. As these absorption features can both fall in the forest region and be difficult to distinguish from Ly α absorption, cosmological analyses eliminate the ∼12–16 per cent of quasars that exhibit BALs. In this paper, we explore an alternate approach that includes BALs in the Ly α autocorrelation function, with the exception of the expected locations of the BAL absorption troughs. This procedure returns over 95 per cent of the path-length that is lost by the exclusion of BALs, as well as increasing the density of sightlines. We show that including BAL quasars reduces the fractional uncertainty in the covariance matrix and correlation function by 12 per cent and does not significantly change the shape of the correlation function relative to analyses that exclude BAL quasars. We also evaluate different definitions of BALs, masking strategies, and potential differences in the quasar continuum in the forest region for BALs with different amounts of absorption.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Covariance matrices for variance-suppressed simulations

ABSTRACT Cosmological N-body simulations provide numerical predictions of the structure of the Universe against which to compare data from ongoing and future surveys, but the growing volume of the Universe mapped by surveys requires correspondingly lower statistical uncertainties in simulations, usually achieved by increasing simulation sizes at the expense of computational power. It was recently proposed to reduce simulation variance without incurring additional computational costs by adopting fixed-amplitude initial conditions. This method has been demonstrated not to introduce bias in various statistics, including the two-point statistics of galaxy samples typically used for extracting cosmological parameters from galaxy redshift survey data, but requires us to revisit current methods for estimating covariance matrices of clustering statistics for simulations. In this work, we find that it is not trivial to construct covariance matrices analytically for fixed-amplitude simulations, but we demonstrate that ezmock (Effective Zel’dovich approximation mock catalogue), the most efficient method for constructing mock catalogues with accurate two- and three-point statistics, provides reasonable covariance matrix estimates for such simulations. We further examine how the variance suppression obtained by amplitude-fixing depends on three-point clustering, small-scale clustering, and galaxy bias, and propose intuitive explanations for the effects we observe based on the ezmock bias model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Baryon acoustic oscillations in thin redshift shells from BOSS DR12 and eBOSS DR16 galaxies

ABSTRACT In an age of large astronomical data sets and severe cosmological tensions, the case for model independent analyses is compelling. We present a set of 14 baryon acoustic oscillations measurements in thin redshift shells with $3\,\mathrm{ per} \,\mathrm{ cent}$ precision that were obtained by analysing BOSS DR12 and eBOSS DR16 galaxies in the redshift range 0.32 < z < 0.66. Thanks to the use of thin shells, the analysis is carried out using just redshifts and angles so that the fiducial model is only introduced when considering the mock catalogues, necessary for the covariance matrix estimation and the pipeline validation. We compare our measurements, with and without supernova data, to the corresponding constraints from Planck 2018, finding good compatibility. A Monte Python module for this likelihood is available at github.com/ranier137/angularBAO.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

IGM damping wing constraints on the tail end of reionization from the enlarged XQR-30 sample

ABSTRACT The attenuation of Lyα photons by neutral hydrogen in the intergalactic medium (IGM) at z ≳ 5 continues to be a powerful probe for studying the epoch of reionization. Given a framework to estimate the intrinsic (true) Lyα emission of high-z sources, one can infer the ionization state of the IGM during reionization. In this work, we use the enlarged XQR-30 sample of 42 high-resolution and high signal-to-noise quasar spectra between $5.8\lesssim \, z\lesssim \, 6.6$ obtained with VLT/X-shooter to place constraints on the IGM neutral fraction. This is achieved using our existing Bayesian QSO reconstruction framework which accounts for uncertainties such as the: (i) posterior distribution of predicted intrinsic Lyα emission profiles (obtained via covariance matrix reconstruction of the Lyα and N v emission lines from unattenuated high-ionization emission line profiles; C iv, Si iv + O iv], and C iii]) and (ii) distribution of ionized regions within the IGM using synthetic damping wing profiles drawn from a 1.63 Gpc3 reionization simulation. Following careful quality control, we used 23 of the 42 available QSOs to obtain constraints/limits on the IGM neutral fraction during the tail-end of reionization. Our median and 68th percentile constraints on the IGM neutral fraction are: $0.20\substack{+0.14 -0.12}$ and $0.29\substack{+0.14 -0.13}$ at z = 6.15 and 6.35. Further, we also report 68th percentile upper limits of $\bar{x}_{\mathrm{H\, {\small I}}{}} \lt 0.21$, 0.20, 0.21, and 0.18 at z = 5.8, 5.95, 6.05, and 6.55. These results imply reionization is still ongoing at $5.8\lesssim \, z\lesssim \, 6.55$, consistent with previous results from XQR-30 (dark fraction and Lyα forest) along with other observational probes considered in the literature.

Greig, B. (ORCID:0000000240852094)↗
DOE OSTI2021

eQuilibrator 3.0: a database solution for thermodynamic constant estimation

Abstract eQuilibrator (equilibrator.weizmann.ac.il) is a database of biochemical equilibrium constants and Gibbs free energies, originally designed as a web-based interface. While the website now counts around 1,000 distinct monthly users, its design could not accommodate larger compound databases and it lacked a scalable Application Programming Interface (API) for integration into other tools developed by the systems biology community. Here, we report on the recent updates to the database as well as the addition of a new Python-based interface to eQuilibrator that adds many new features such as a 100-fold larger compound database, the ability to add novel compounds, improvements in speed and memory use, and correction for Mg2+ ion concentrations. Moreover, the new interface can compute the covariance matrix of the uncertainty between estimates, for which we show the advantages and describe the application in metabolic modelling. We foresee that these improvements will make thermodynamic modelling more accessible and facilitate the integration of eQuilibrator into other software platforms.

59 BASIC BIOLOGICAL SCIENCES↗