Search NASA⌕ Search

SEARCH · Search NASA

Results for “Constraining”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Coupled Cluster Theory for Nonadiabatic Dynamics: Nuclear Gradients and Nonadiabatic Couplings in Similarity Constrained Coupled Cluster Theory

Coupled cluster theory is one of the most accurate electronic structure methods for predicting ground and excited state chemistry. However, the presence of numerical artifacts at electronic degeneracies, such as complex energies, has made it difficult to apply the method in nonadiabatic dynamics simulations. While it has already been shown that such numerical artifacts can be fully removed by using similarity constrained coupled cluster (SCC) theory [J. Phys. Chem. Lett. 2017, 8(19), 4801–4807], simulating dynamics requires efficient implementations of gradients and nonadiabatic couplings. Here, we present an implementation of nuclear gradients and nonadiabatic derivative couplings at the similarity constrained coupled cluster singles and doubles (SCCSD) level of theory, thereby making possible nonadiabatic dynamics simulations using a coupled cluster theory that provides a correct description of conical intersections between excited states. We present a few numerical examples that show good agreement with literature values and discuss some limitations of the method.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Observationally Constrained Modeling of the Reactive Uptake of Isoprene-Derived Epoxydiols under Elevated Relative Humidity and Varying Acidity of Seed Aerosol Conditions

Isoprene is the non-methane volatile organic compound (VOC) emitted in largest amounts to the atmosphere, and it is a significant source of secondary organic aerosol (SOA) mass. The uptake of isoprene oxidation products followed by multiphase chemistry in fine particles is the key pathway to form isoprene epoxydiol-derived SOA (IEPOX-SOA). However, many parameters that relate to diffusion and reaction of IEPOX in the particle phase remain uncertain, since reaction kinetics previously measured in bulk aqueous phase solutions might be different from atmospheric aerosols. Here, we use simultaneous environmental chamber measurements of multiple parameters governing IEPOX-SOA formation at timescales of ~hours: particle size distribution, composition, and volatility of IEPOX-SOA to constrain the key parameters governing IE-POX-SOA formation under humid (i.e., 50% relative humidity, RH) and varying seed aerosol acidity conditions. Reducing the 2-methyltetrol (tetrol) reaction rate constants by a factor of 4 brings the model predictions in agreement with the IEPOX-SOA measurements with acidified ammonium bisulfate seed aerosols. For less acidic ammonium sulfate aerosols both the organo-sulfate (OS) and tetrol reaction rate constants need to be reduced to bring model predictions closer to chamber observations. Using the measured non-volatile content of IEPOX-SOA we constrain the oligomerization timescale of 2-methyltetrols. We find that the oligomerization timescale is 4 hours with acidified seed aerosols, but a much longer time scale of 24 hours is needed for non-acidified seed aerosols, indicating that the aerosol acidity greatly affects the oligomerization rate of tetrols. Here we show that the actual kinetics of IEPOX-SOA formation rate on aerosol seeds consisting of both ammonium bisulfate and ammonium sulfate are a factor of 4~5 slower under 50-60% RH conditions compared to their application in previous models, which were based on bulk aqueous solution measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

A Strong Anthropogenic Black Carbon Forcing Constrained by Pollution Trends Over China

Estimates of the effective radiative forcing from aerosol-radiation interaction (ERF ari ) of anthropogenic Black Carbon (BC) have been disputable and require better constraints. Here we find a substantial decline in atmospheric absorption of −5.79Wm −2 decade −1 over eastern central China (ECC) responding to recent anthropogenic BC emission reductions. By combining the observational finding with advances from Coupled Model Intercomparison Project phase6 (CMIP6), we identify an emergent constraint on the ERF ari of anthropogenic BC. We show that across CMIP6 models the simulated trends correlate well with simulated annual mean shortwave atmospheric absorption by anthropogenic BC over China. Making use of this emergent relationship allows us to constrain the aerosol absorption optical depth of anthropogenic BC and further provide a constrained range of 2.4–3.0 Wm −2 for its top-of-atmosphere ERF ari over China, higher than existing estimates. Furthermore, our work supports a strong warming effect of BC over China, and highlights the need to improve BC simulations over source regions.

CMIP6 models↗
DOE OSTI2023

Emergent Constrained Projections of Mean and Extreme Warming in China

Reliable regional temperature projections including heat extremes are essential for climate change adaptation and mitigation. Taking China as an example, simple averages from Coupled Model Intercomparison Project Phase 6 (CMIP6) models project high warming due to sampling many high climate sensitivities in the ensemble. Here, we develop an emergent constraint (EC) framework to obtain constrained mean and daily maximum temperature (TXx) warming over China by using observed global warming and local residual warming. The constrained annual mean and TXx warming over China (2.33°C [1.61–3.05°C] and 2.31°C [1.21–2.99°C]) are 0.65°C [0.04–1.76°C] and 0.63°C [–0.50–2.39°C], respectively, lower than raw projections (2.98°C [1.85–4.22°C] and 2.94°C [2.04–4.39°C]) for 2080–2099 under the intermediate-emission scenario. Approximately half model uncertainty is reduced after constraint. The land area (population) experiencing temperature extremes in our metric is 78% (85%) of the raw projections. Our results imply a lower impact of extreme heat than implied by current raw CMIP6 projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Land Use Change Alters Soil Organic Carbon: Constrained Global Patterns and Predictors

Abstract Land use change (LUC) alters the global carbon (C) stock, but our estimation of the alteration remains uncertain and is a major impediment to predicting the global C cycle. The uncertainty is partly due to the limited number and geographical bias of observations, and limited exploration of its predictors. Here we generated a comprehensive global database of 5,980 observations from 790 articles. The number of sites evaluated is at least seven times larger than in previous meta‐analyses. Our constrained estimates of different LUC's effects on soil organic C (SOC) and their variations across global climates reveal underestimation/overestimation in previous estimates. Converting forests and grasslands to croplands reduced SOC by 24.5% ± 1.53% (−11.03 ± 1.06 Mg ha −1 ) and 22.7% ± 1.22% (−8.09 ± 0.67 Mg ha −1 ), while 28.0% ± 1.56% (4.46 ± 0.42 Mg ha −1 ) and 33.5% ± 1.68% (5.8 ± 0.38 Mg ha −1 ) increases, respectively, were obtained in the reverse processes. Converting forests to grasslands decreased SOC by 2.1% ± 1.22% (−1.13 ± 0.44 Mg ha −1 ), while the reverse process increased SOC by 18.6% ± 1.73% (3.31 ± 0.51 Mg ha −1 ). Modeled relative importance of 10 drivers of LUC's impact on SOC revealed that higher initial SOC (iSOC) does not solely determine SOC loss in SOC‐negative LUC scenarios as previously proposed. Across four decades, reconverting croplands to forests and grasslands recovered only 49.5% (6.1 ± 0.51 Mg ha −1 ) and 75.3% (7.0 ± 0.38 Mg ha −1 ) of the iSOC, respectively, indicating the need for protecting C‐rich ecosystems. Our global data set advances information on LUC's effect on SOC and can be valuable to constrain Earth system models to reliably estimate global SOC stocks and plan climate change mitigation strategies.

Environmental Sciences & Ecology↗
DOE OSTI2024

Learning Constitutive Relations From Soil Moisture Data via Physically Constrained Neural Networks

Abstract The constitutive relations of the Richardson‐Richards equation encode the macroscopic properties of soil water retention and conductivity. These soil hydraulic functions are commonly represented by models with a handful of parameters. The limited degrees of freedom of such soil hydraulic models constrain our ability to extract soil hydraulic properties from soil moisture data via inverse modeling. We present a new free‐form approach to learning the constitutive relations using physically constrained neural networks. We implemented the inverse modeling framework in a differentiable modeling framework, JAX, to ensure scalability and extensibility. For efficient gradient computations, we implemented implicit differentiation through a nonlinear solver for the Richardson‐Richards equation. We tested the framework against synthetic noisy data and demonstrated its robustness against varying magnitudes of noise and degrees of freedom of the neural networks. We applied the framework to soil moisture data from an upward infiltration experiment and demonstrated that the neural network‐based approach was better fitted to the experimental data than a parametric model and that the framework can learn the constitutive relations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Constraining maximum event magnitude during injection-triggered seismicity

Abstract Understanding mechanisms controlling fluid injection-triggered seismicity is key in defining strategies to ameliorate it. Recent triggered events (e.g. Pohang, Mw 5.5) have exceeded predictions of average energy release by a factor of >1000x, necessitating robust methodologies to both define critical antecedent conditions and to thereby constrain anticipated event size. We define maximum event magnitudes resulting from triggering as a function of pre-existing critical stresses and fluid injection volume. Fluid injection experiments on prestressed laboratory faults confirm these estimates of triggered moment magnitudes for varied boundary conditions and injection rates. In addition, observed ratios of shear slip to dilation rates on individual faults signal triggering and may serve as a measurable proxy for impending rupture. This new framework provides a robust method of constraining maximum event size for preloaded faults and unifies prior laboratory and field observations that span sixteen decades in injection volume and four decades in length scale.

58 GEOSCIENCES↗
DOE OSTI2022

Machine learning-based observation-constrained projections reveal elevated global socioeconomic risks from wildfire

Reliable projections of wildfire and associated socioeconomic risks are crucial for the development of efficient and effective adaptation and mitigation strategies. The lack of or limited observational constraints for modeling outputs impairs the credibility of wildfire projections. Here, we present a machine learning framework to constrain the future fire carbon emissions simulated by 13 Earth system models from the Coupled Model Intercomparison Project phase 6 (CMIP6), using historical, observed joint states of fire-relevant variables. During the twenty-first century, the observation-constrained ensemble indicates a weaker increase in global fire carbon emissions but higher increase in global wildfire exposure in population, gross domestic production, and agricultural area, compared with the default ensemble. Such elevated socioeconomic risks are primarily caused by the compound regional enhancement of future wildfire activity and socioeconomic development in the western and central African countries, necessitating an emergent strategic preparedness to wildfires in these countries.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Unleashed from constrained optimization: quantum computing for quantum chemistry employing generator coordinate inspired method

Hybrid quantum-classical approaches offer potential solutions to quantum chemistry problems, yet they often manifest as constrained optimization problems. Here, we explore the interconnection between constrained optimization and generalized eigenvalue problems through the Unitary Coupled Cluster (UCC) excitation generators. Inspired by the generator coordinate method, we employ these UCC excitation generators to construct non-orthogonal, overcomplete many-body bases, projecting the system Hamiltonian into an effective Hamiltonian, which bypasses issues such as barren plateaus that heuristic numerical minimizers often encountered in standard variational quantum eigensolver (VQE). Diverging from conventional quantum subspace expansion methods, we introduce an adaptive scheme that robustly constructs the many-body basis sets from a pool of the UCC excitation generators. This scheme supports the development of a hierarchical ADAPT quantum-classical strategy, enabling a balanced interplay between subspace expansion and ansatz optimization to address complex, strongly correlated quantum chemical systems cost-effectively, setting the stage for more advanced quantum simulations in chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Constrained quantum optimization for extractive summarization on a trapped-ion quantum computer

Abstract Realizing the potential of near-term quantum computers to solve industry-relevant constrained-optimization problems is a promising path to quantum advantage. In this work, we consider the extractive summarization constrained-optimization problem and demonstrate the largest-to-date execution of a quantum optimization algorithm that natively preserves constraints on quantum hardware. We report results with the Quantum Alternating Operator Ansatz algorithm with a Hamming-weight-preserving XY mixer (XY-QAOA) on trapped-ion quantum computer. We successfully execute XY-QAOA circuits that restrict the quantum evolution to the in-constraint subspace, using up to 20 qubits and a two-qubit gate depth of up to 159. We demonstrate the necessity of directly encoding the constraints into the quantum circuit by showing the trade-off between the in-constraint probability and the quality of the solution that is implicit if unconstrained quantum optimization methods are used. We show that this trade-off makes choosing good parameters difficult in general. We compare XY-QAOA to the Layer Variational Quantum Eigensolver algorithm, which has a highly expressive constant-depth circuit, and the Quantum Approximate Optimization Algorithm. We discuss the respective trade-offs of the algorithms and implications for their execution on near-term quantum hardware.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Beyond the two-point correlation: Constraining primordial non-Gaussianity with density-perturbation moments

Constraining primordial non-Gaussianity (PNG) on the large-scale cosmic structure (LSS) is an important step in understanding properties of the early Universe, specifically in distinguishing between different inflationary models. Measuring PNG relies on evaluating the scale-dependent correlations in the density field. New summary statistics beyond the two- and three-point correlation functions in configuration space and their Fourier-space counterparts, the power- and bispectrum may provide increased sensitivity. We introduce a new method for extracting the PNG signal imprinted on the LSS by using the first three Gaussian moments of the normalized correlation in density perturbations, evaluated on varying distance scales. We aim to assess this method’s sensitivity to local PNG, parameterized by f NL . We performed spherical convolutions on a range of scales on dark-matter-halo simulations to measure the scale-dependent correlations in the density field. From these, we computed the first three moments and compared them to a model expectation vector, parameterized to the second power in f NL . Our method provides about 21% improvement in sensitivity to f NL with respect to using the two-point correlation function alone. Notably, we find that the second moment alone carries nearly as much constraining power as the mean, highlighting the potential of higher order statistics. Given its simplicity and efficiency, this framework is well suited for application to current and upcoming large-scale surveys such as the Dark Energy Spectroscopic Instrument (DESI).

early universe↗
DOE OSTI2023

Constraining Nucleosythesis in Neutrino-Driven Winds Using the Impact of ({alpha}, xn) Reaction Rates

The lighter heavy elements of the first r-process peak, between strontium and silver, can be synthesized in the moderately neutron-rich neutrino-driven ejecta of either core-collapse supernovae or neutron star mergers via the weak r-process. This nucleosynthesis scenario exhibits uncertainties from the absence of experimental data from (alpha, xn) reactions on neutron-rich nuclei, which are currently based on statistical model estimates. We have performed a new impact study to identify the most important (alpha, xn) reactions that can affect the production of the lighter heavy elements under different astrophysical conditions using new, constrained (alpha, xn) reaction rates based on the Atomki-V2 alpha OMP. Our results show how when reducing the nuclear physics uncertainties, we can use abundance ratios to constrain the astrophysical conditions/environment. This can be achieved in the near future, when the key (alpha, xn) reaction rates will be measured experimentally in radioactive beam facilities.

Psaltis, Athanasios↗
DOE OSTI2025

Fast methods for multisite charge transfer processes. I. Constrained, state averaged CASSCF(1,n) and CASSCF(2n − 1,n) simulations

We design a dynamically weighted state-averaged constrained complete active space self-consistent field (DW-SA-cCASSCF) algorithm to treat electrons or holes moving between n molecular fragments (where n can be larger than 2). Within such a so-called eDSCn/hDSCn approach, we consider configurations that are mutually single excitations of each other, and we apply a generalized set of constraints to tailor the method for studying charge transfer problems. The constrained optimization problem is efficiently solved using a DIIS-SQP algorithm, thus maintaining computational efficiency. We demonstrate the method for a finite Su–Schrieffer–Heeger chain, successfully reproducing the expected exponential decay of diabatic couplings with distance. When combined with a gradient, the current extension immediately enables efficient nonadiabatic dynamics simulations of complex multi-state charge transfer processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Constraining the sources of ultra-high-energy cosmic rays across and above the ankle with the spectrum and composition data measured at the Pierre Auger Observatory

In this work we present the interpretation of the energy spectrum and mass composition data as measured by the Pierre Auger Collaboration above 6 × 10$^{17}$ eV. We use an astrophysical model with two extragalactic source populations to model the hardening of the cosmic-ray flux at around 5 × 10$^{18}$ eV (the so-called “ankle” feature) as a transition between these two components. We find our data to be well reproduced if sources above the ankle emit a mixed composition with a hard spectrum and a low rigidity cutoff. The component below the ankle is required to have a very soft spectrum and a mix of protons and intermediate-mass nuclei. The origin of this intermediate-mass component is not well constrained and it could originate from either Galactic or extragalactic sources.To the aim of evaluating our capability to constrain astrophysical models, we discuss the impact on the fit results of the main experimental systematic uncertainties and of the assumptions about quantities affecting the air shower development as well as the propagation and redshift distribution of injected ultra-high-energy cosmic rays (UHECRs).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Constraining multi-field inflation using the SPHEREx all-sky survey power spectra

Abstract We investigate how well the SPHEREx all-sky survey can constrain local primordial non-Gaussianity beyond the parameterf NL using galaxy power spectra. We forecast joint constraints on the parametersf NL ,g NL andτ NL obtained assuming a simple two-field curvaton model of inflation. The parametersf NL andg NL characterise the squeezed limits of the primordial bispectrum and trispectrum respectively, and lead to a characteristic scale-dependence of the galaxy bias that increases out to arbitrarily large scales. Values of the parameterτ NL > (6/5f NL ) 2 cause the galaxy power spectrum to have a stochastic component which also increases out to arbitrarily large scales. Our MCMC forecasts indicate that SPHEREx can provide joint constraints on any two of the three parametersf NL ,g NL andτ NL . Due to strong degeneracies among these parameters, measurements of the galaxy power spectra alone may not be sufficient to jointly constrain all three. Constraints onf NL ,g NL andτ NL obtained from galaxy power spectrum observations depend on the modelling of underlying nuisance parameters. We study the robustness of our forecast constraints to modelling choices and note that even with relatively conservative modelling assumptions, SPHEREx galaxy power spectra can provide strong evidence of local non-Gaussianity, even if the particular values off NL andg NL cannot be measured precisely.

Astronomy & Astrophysics↗
DOE OSTI2021

Asynchronous quadratic control for constrained hidden markov jump linear systems with incomplete MTPM and MOCPM

Abstract This paper investigates the quadratic optimal control problem for constrained Markov jump linear systems with incomplete mode transition probability matrix (MTPM). Considering original system mode is not accessible, observed mode is utilized for asynchronous controller design where mode observation conditional probability matrix (MOCPM), which characterizes the emission between original modes and observed modes is assumed to be partially known. An LMI optimization problem is formulated for such constrained hidden Markov jump linear systems with incomplete MTPM and MOCPM. Based on this, a feasible state-feedback controller can be designed with the application of free-connection weighting matrix method. The desired controller, dependent on observed mode, is an asynchronous one which can minimize the upper bound of quadratic cost and satisfy restrictions on system states and control variables. Furthermore, clustering observation where observed modes recast into several clusters, is explored for simplifying the computational complexity. Numerical examples are provided to illustrate the validity.

Zhu, Jin↗
DOE OSTI2021

Blast from the past: constraining progenitor models of SN 1972E

ABSTRACT We present a novel technique to study Type Ia supernovae (SNe Ia) by constraining surviving companions of historical extragalactic SN by combining archival photographic plates and Hubble Space Telescope(HST) imaging. We demonstrate this technique for Supernova 1972E, the nearest known SN Ia in 125 yr. Some models of SNe Ia describe a white dwarf with a non-degenerate companion that donates enough mass to trigger thermonuclear detonation. Hydrodynamic simulations and stellar evolution models show that these donor stars will survive the explosion, and show increased luminosity for at least a 1000 yr. Thus, late-time observations of the exact location of a supernova can constrain the presence of a surviving donor star and progenitor models. We find the explosion site of SN 1972E by analysing 17 digitized photographic plates taken with the European Southern Observatory 1-m Schmidt and 1 plate taken with the Cerro Tololo Inter-American Observatory 1.5-m telescope. Using the Gaia eDR3 catalogue to determine Supernova 1972E’s location yields: α = 13h39m52${_{.}^{\rm s}}$708 ± 0${_{.}^{\rm s}}$004 and δ = −31°40’9${_{.}^{\prime\prime}}$00 ± 0${_{.}^{\prime\prime}}$04 (ICRS). In 2005, HST/ACS imaged the host galaxy of SN 1972E with the F435W, F555W, and F814W filters covering the explosion site. The nearest detected source is offset by 3.0 times our positional precision, and is inconsistent with the colours expected of a surviving donor star. Thus, the limiting magnitude of the HST observation (F555W > 28 mag) rules out all He star companion models and the most luminous main-sequence companion model currently in the literature. The remaining main-sequence companion models could be tested with a 10 orbit HST exposure in the F606W filter.

Do, Aaron↗
DOE OSTI2021

Dark energy survey year 1 results: Constraining baryonic physics in the Universe

ABSTRACT Measurements of large-scale structure are interpreted using theoretical predictions for the matter distribution, including potential impacts of baryonic physics. We constrain the feedback strength of baryons jointly with cosmology using weak lensing and galaxy clustering observables (3 × 2pt) of Dark Energy Survey (DES) Year 1 data in combination with external information from baryon acoustic oscillations (BAO) and Planck cosmic microwave background polarization. Our baryon modelling is informed by a set of hydrodynamical simulations that span a variety of baryon scenarios; we span this space via a Principal Component (PC) analysis of the summary statistics extracted from these simulations. We show that at the level of DES Y1 constraining power, one PC is sufficient to describe the variation of baryonic effects in the observables, and the first PC amplitude (Q1) generally reflects the strength of baryon feedback. With the upper limit of Q1 prior being bound by the Illustris feedback scenarios, we reach $\sim 20{{\ \rm per\ cent}}$ improvement in the constraint of $S_8=\sigma _8(\Omega _{\rm m}/0.3)^{0.5}=0.788^{+0.018}_{-0.021}$ compared to the original DES 3 × 2pt analysis. This gain is driven by the inclusion of small-scale cosmic shear information down to 2.5 arcmin, which was excluded in previous DES analyses that did not model baryonic physics. We obtain $S_8=0.781^{+0.014}_{-0.015}$ for the combined DES Y1+Planck EE+BAO analysis with a non-informative Q1 prior. In terms of the baryon constraints, we measure $Q_1=1.14^{+2.20}_{-2.80}$ for DES Y1 only and $Q_1=1.42^{+1.63}_{-1.48}$ for DESY1+Planck EE+BAO, allowing us to exclude one of the most extreme AGN feedback hydrodynamical scenario at more than 2σ.

79 ASTRONOMY AND ASTROPHYSICS↗