SEARCH · Search NASA

Results for “Structural error”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,693 records · Page 94

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Non-empirical, low-cost recovery of exact conditions with model-Hamiltonian inspired expressions in jmDFT

Density functional theory (DFT) is widely applied to both molecules and materials, but well known energetic delocalization and static correlation errors in practical exchange-correlation approximations limit quantitative accuracy. Common methods that correct energetic delocalization errors, such as the Hubbard U correction in DFT+U or Hartree-Fock exchange in global hybrids, do so at the cost of worsening static correlation errors. We recently introduced an alternate approach [Bajaj et al., J. Chem. Phys. 147, 191101 (2017)] known as judiciously modified DFT (jmDFT), wherein the deviation from exact behavior of semilocal functionals over both fractional spin and charge, i.e., the so-called flat plane, was used to motivate functional forms of second order analytic corrections. In this work, we introduce fully nonempirical expressions for all four coefficients in a DFT+U+J-inspired form of jmDFT, where all coefficients are obtained only from energies and eigenvalues of the integer-electron systems. We show good agreement for U and J coefficients obtained nonempirically as compared with the results of numerical fitting in a jmDFT U+J/J' correction. Incorporating the fully nonempirical jmDFT correction reduces and even eliminates the fractional spin error at the same time as eliminating the energetic delocalization error. We show that this approach extends beyond s-electron systems to higher angular momentum cases including p- and d-electrons. Finally, we diagnose some shortcomings of the current jmDFT approach that limit its ability to improve upon DFT results for cases such as weakly bound anions due to poor underlying semilocal functional behavior.

Delocalization↗

DOE OSTI2020

Reaction Profiles and Kinetics for Radical–Radical Hydrogen Abstraction via Multireference Coupled Cluster Theory

Radical-radical abstractions in hydrocarbon oxidation chemistry are disproportionation reactions that are generally exothermic with little or no barrier yet are underappreciated and poorly studied. Such challenging multireference electronic structure problems are tackled here using the recently developed state-specific multireference coupled cluster methods Mk-MRCCSD and Mk-MRCCSD(T), as well as the companion perturbation theory Mk-MRPT2 and the established MRCISD, MRCISD+Q, and CASPT2 approaches. Reaction paths are investigated for five prototypes involving radical-radical hydrogen abstraction: H + BeH → H 2 + Be, H + NH 2 → H 2 + NH, CH 3 + C 2 H 5 → CH 4 + C 2 H 4 , H + C2H 5 → H 2 + C 2 H 4 , and H + HCO → H 2 + CO. Full configuration interaction (FCI) benchmark computations for the H + BeH, H + NH 2 , and H + HCO reactions prove that Mk-MRCCSD(T) provides superior accuracy for the interaction energies in the entrance channel, with mean absolute errors less than 0.3 kcal mol(-1) and percentage deviations less than 10% over the fragment separations of relevance to kinetics. To facilitate combustion studies, energetics for the CH 3 + C 2 H 5 , H + C 2 H 5 , and H + HCO reactions were computed at each level of theory with correlation-consistent basis sets (cc-pVXZ, X = T, Q 5) and extrapolated to the complete basis set (CBS) limit. These CBS energies were coupled with CASPT2 projected vibrational frequencies along a minimum energy path to obtain rate constants for these three reactions. The rigorous Mk-MRCCSD(T)/CBS results demonstrate unequivocally that these three reactions proceed with no barrier in the entrance channel, contrary to some earlier predictions. Mk-MRCCSD(T) also reveals that the economical CASPT2 method performs well for large interfragment separations but may deteriorate substantially at shorter distances.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Flow-based likelihoods for non-Gaussian inference

We investigate the use of data-driven likelihoods to bypass a key assumption made in many scientific analyses, which is that the true likelihood of the data is Gaussian. In particular, we suggest using the optimization targets of flow-based generative models, a class of models that can capture complex distributions by transforming a simple base distribution through layers of nonlinearities. We call these flow-based likelihoods (FBL). We analyze the accuracy and precision of the reconstructed likelihoods on mock Gaussian data, and show that simply gauging the quality of samples drawn from the trained model is not a sufficient indicator that the true likelihood has been learned. We nevertheless demonstrate that the likelihood can be reconstructed to a precision equal to that of sampling error due to a finite sample size. We then apply FBLs to mock weak lensing convergence power spectra, a cosmological observable that is significantly non-Gaussian (NG). We find that the FBL captures the NG signatures in the data extremely well, while other commonly used data-driven likelihoods, such as Gaussian mixture models and independent component analysis, fail to do so. This suggests that works that have found small posterior shifts in NG data with data-driven likelihoods such as these could be underestimating the impact of non-Gaussianity in parameter constraints. By introducing a suite of tests that can capture different levels of NG in the data, we show that the success or failure of traditional data-driven likelihoods can be tied back to the structure of the NG in the data. Here, unlike other methods, the flexibility of the FBL makes it successful at tackling different types of NG simultaneously. Because of this, and consequently their likely applicability across datasets and domains, we encourage their use for inference when sufficient mock data are available for training.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Structure of the Wall Sheath in a Gas-Discharge Plasma for an Arbitrary Orientation of a Flat Probe Relative to the Electric Field in the Plasma

We investigate the structure of the wall sheath of a gas discharge near a flat surface at a negative potential for high mean electron energy. It is shown that in the conditions where the mean energy of ions in the plasma is much lower than the mean electron energy, the parameters of the wall sheath weakly depend on the mutual orientation of the normal to the surface and the electric field in the plasma for an arbitrary ratio of the Debye radius to the ion mean free path relative to the resonant charge exchange process. It is found that for inert gases (He, Ar) for ratio E/P of the electric field to pressure exceeding 10 V/(cm Torr) in the plasma, the disregard of ionization in the perturbed wall sheath can lead to substantial errors in the calculation of its parameters. It is shown that the ionization leads to an increase in the electric field in the wall sheath and, as a consequence, to an increase in the mean velocity of ions at the boundary between the quasi-neutral presheath and the part of the perturbed wall sheath in which quasi-neutrality is substantially violated. The parameters of the wall sheath where quasi-neutrality is significantly violated depend on the ionization rate much less strongly than the corresponding parameters of the quasi-neutral presheath. We have determined the relation for concentration of charged particles in the unperturbed plasma from the ion saturation current considering the actual ion energy distribution function in the plasma as well as ionization in the presheath and the part of the perturbed wall sheath in which quasi-neutrality is violated significantly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Reaction mechanisms of a cyclic ether intermediate: Ethyloxirane

Oxiranes are a class of cyclic ethers formed in abundance during low-temperature combustion of hydrocarbons and biofuels, either via chain-propagating steps that occur from unimolecular decomposition of β -hydroperoxyalkyl radicals ( β -˙QOOH) or from reactions of HOÓ with alkenes. Ethyloxirane is one of four alkyl-substituted cyclic ether isomers produced as an intermediate from n -butane oxidation. While rate coefficients for β -˙QOOH → ethyloxirane + ÓH are reported extensively, subsequent reaction mechanisms of the cyclic ether are not. As a result, chemical kinetics mechanisms commonly adopt simplified chemistry to describe ethyloxirane consumption by convoluting several elementary reactions into a single step, which may introduce mechanism truncation error—uncertainty derived from missing or incomplete chemistry. Here we provide fundamental insight on reaction mechanisms of ethyloxirane in support of ongoing efforts to minimize mechanism truncation error. Reaction mechanisms are inferred from the detection of products during chlorine atom-initiated oxidation experiments using multiplexed photoionization mass spectrometry conducted at 10 Torr and temperatures of 650 K and 800 K. To complement the experiments, calculations of stationary point energies were conducted using the ccCA-PS3 composite method on ˙R + O 2 potential energy surfaces for the four ethyloxiranyl radical isomers, which produced barrier heights for 24 reaction pathways. In addition to products from ˙QOOH → cyclic ether + ÓH and ˙R + O 2 → conjugate alkene + HOÓ, both of which were significant pathways and are prototypical to alkane oxidation, other species were identified from ring-opening of both ethyloxiranyl and ˙QOOH radicals. The latter occurs when the unpaired electron is localized on the ether group, causing the initial ˙QOOH structure to ring-open and form a resonance-stabilized ketohydroperoxide-type radical. The present work provides the first analysis of ethyloxirane oxidation chemistry, which reveals that consumption pathways are complex and may require an expansion of submechanisms to increase the fidelity of chemical kinetics mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Shape restricted additive hazards models: Monotone, unimodal, and U‐shape hazard functions

We consider estimation of the semiparametric additive hazards model with an unspecified baseline hazard function where the effect of a continuous covariate has a specific shape but otherwise unspecified. Such estimation is particularly useful for a unimodal hazard function, where the hazard is monotone increasing and monotone decreasing with an unknown mode. A popular approach of the proportional hazards model is limited in such setting due to the complicated structure of the partial likelihood. Our model defines a quadratic loss function, and its simple structure allows a global Hessian matrix that does not involve parameters. Thus, once the global Hessian matrix is computed, a standard quadratic programming method can be applicable by profiling all possible locations of the mode. However, the quadratic programming method may be inefficient to handle a large global Hessian matrix in the profiling algorithm due to a large dimensionality, where the dimension of the global Hessian matrix and number of hypothetical modes are the same order as the sample size. We propose the quadratic pool adjacent violators algorithm to reduce computational costs. The proposed algorithm is extended to the model with a time‐dependent covariate with monotone or U‐shape hazard function. In simulation studies, our proposed method improves computational speed compared to the quadratic programming method, with bias and mean square error reductions. We analyze data from a recent cardiovascular study.

Mathematical & Computational Biology↗

DOE OSTI2020

Remaining Useful Strength (RUS) Prediction of SiCf-SiCm Composite Materials Using Deep Learning and Acoustic Emission

Prognosis techniques for prediction of remaining useful life (RUL) are of crucial importance to the management of complex systems for they can lead to appropriate maintenance interventions and improvements in reliability. While various data-driven methods have been introduced to predict the remaining useful life (RUL) of machinery systems or batteries, no research has been reported on the remaining useful strength (RUS) prediction of silicon carbide fiber reinforced silicon carbide matrix (SiCf-SiCm) materials with pivotal role in its potential usage as a structural material in nuclear reactors and turbine engines. Knowledge of its degradation process is of the utmost importance to the manufacturers. For this purpose, two approaches based on the machine-learning techniques of random-forest (RF) and convolutional neural network (CNN) are proposed to predict the RUS of SiCf-SiCm using only acoustic emission (AE) signals generated during the material’s stress applying process. Experimental results show that the CNN models achieved better predictive performance than the RF models but the latter with expert-engineered features achieves better prediction for AE signals in the early stage of degradation. Additionally, our results demonstrate that both models can correctly predict the SiCf-SiCm RUS as evaluated by our robust testing method from which the best average root mean square error (RMSE) and Pearson correlation coefficient of 3.55 ksi units and 0.85 were obtained.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Local gravitational redshifts can bias cosmological measurements

Measurements of cosmological parameters via the distance-redshift relation usually rely on models that assume a homogenous universe. It is commonly presumed that the large-scale structure evident in our Universe has a negligible impact on the measurement if distances probed in observations are sufficiently large (compared to the scale of inhomogeneities) and are averaged over different directions on the sky. This presumption does not hold when considering the effect of the gravitational redshift caused by our local gravitational potential, which alters light coming from all distances and directions in the same way. Despite its small magnitude, this local gravitational redshift gives rise to noticeable effects in cosmological inference using SN Ia data. Assuming conservative prior knowledge of the local potential given by sampling a range of gravitational potentials at locations of Milky-Way-like galaxies identified in cosmological simulations, we show that ignoring the gravitational redshift effect in a standard data analysis leads to an additional systematic error of ∼1% in the determination of density parameters and the dark energy equation of state. We conclude that our local gravitational field affects our cosmological inference at a level that is important in future observations aiming to achieve percent-level accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Delayed Gamma Measurements in Different Nuclear Research Reactors Bringing Out the Importance of the Delayed Contribution in Gamma Flux Calculations

Neutron and gamma flux levels are key parameters in nuclear research reactors. In Material Testing Reactors, such as the future Jules Horowitz Reactor, under construction at the French Alternative Energies and Atomic Energy Commission (CEA Cadarache, France), the expected gamma flux levels are very high (nuclear heating is of the order of 20 W/g at 100 MWth). As gamma rays deposit their energy in the reactor structures and structural materials it is important to take them into account when designing irradiation devices. There are only a few sensors which allow measurements of the nuclear heating ; a recent development at the CEA Cadarache allows measurements of the gamma flux using a miniature ionization chamber (MIC). The measured MIC response is often compared with calculation using modern Monte Carlo (MC) neutron and photon transport codes, such as TRIPOLI-4 and MCNP6. In these calculations only the production of prompt gamma rays in the reactor is usually modelled thus neglecting the delayed gamma rays. Hence calculations and measurements are usually in better accordance for the neutron flux than for the gamma flux. In this paper we study the contribution of delayed gamma rays to the total MIC signal in order to estimate the systematic error in gamma flux MC calculations. In order to experimentally determine the delayed gamma flux contributions to the MIC response, we performed gamma flux measurements with CEA developed MIC at three different research reactors: the OSIRIS reactor (MTR - 70 MWth at CEA Saclay, France), the TRIGA MARK II reactor (TRIGA - 250 kWth at the Jozef Stefan Institute, Slovenia) and the MARIA reactor (MTR - 30 MWth at the National Center for Nuclear Research, Poland). In order to experimentally assess the delayed gamma flux contribution to the total gamma flux, several reactor shut down (scram) experiments were performed specifically for the purpose of the measurements. Results show that on average about 30 % of the MIC signal is due to the delayed gamma rays. In this paper we describe experiments in each of the three reactors and how we estimate delayed gamma rays with MIC measurements. The results and perspectives are discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2016

Marker-free registration of forest terrestrial laser scanner data pairs with embedded confidence metrics

Terrestrial laser scanning (TLS) has emerged as an effective tool for rapid comprehensive measurement of object structure. Registration of TLS data is an important prerequisite to overcome the limitations of occlusion. However, due to the high dissimilarity of point cloud data collected from disparate viewpoints in the forest environment, adequate marker-free registration approaches have not been developed. The majority of studies instead rely on the utilization of artificial tie points (e.g., reflective tooling balls) placed within a scene to aid in coordinate transformation. We present a technique for generating view-invariant feature descriptors that are intrinsic to the point cloud data and, thus, enable blind marker-free registration in forest environments. To overcome the limitation of initial pose estimation, we employ a voting method to blindly determine the optimal pairwise transformation parameters, without an a priori estimate of the initial sensor pose. To provide embedded error metrics, we developed a set theory framework in which a circular transformation is traversed between disjoint tie point subsets. This provides an upper estimate of the Root Mean Square Error (RMSE) confidence associated with each pairwise transformation. Output RMSE errors are commensurate with the RMSE of input tie points locations. Thus, while the mean output RMSE=16.3cm, improved results could be achieved with a more precise laser scanning system. This study 1) quantifies the RMSE of the proposed marker-free registration approach, 2) assesses the validity of embedded confidence metrics using receiver operator characteristic (ROC) curves, and 3) informs optimal sample spacing considerations for TLS data collection in New England forests. Furthermore, while the implications for rapid, accurate, and precise forest inventory are obvious, the conceptual framework outlined here could potentially be extended to built environments.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Collagen Turnover in Relation to Risk Factors and Hemodynamics in Human Intracranial Aneurysms

Determinants for molecular and structural instability, that is, impending growth or rupture, of intracranial aneurysms (IAs) remain uncertain. To elucidate this, we endeavored to estimate the actual turnover rates of the main molecular constituent in human IA (collagen) on the basis of radiocarbon ( 14 C) birth dating in relation to IA hemodynamics. Collagen turnover rates in excised human IA samples were calculated using mathematical modeling of 14 C birth dating data of collagen in relation to risk factors and histological markers for collagen maturity/turnover in selected IA. Hemodynamics were simulated using image-based computational fluid dynamics. Correlation, logistic regression, and receiver operating characteristic analyses were performed. Collagen turnover rates were estimated in 46 IA (43 patients); computational fluid dynamics could be performed in 20 IA (20 patients). The mean collagen turnover rate (γ) constituted 126% (±1% error) per year. For patients with arterial hypertension, γ was greater than 2600% annually, whereas γ was distinctly lower with 32% (±1% error) per year for patients without risk factors, such as smoking and hypertension. There was a distinct association between histological presence of rather immature collagen in human IA and the presence of modifiable risk factors. Spatial-temporal averaged wall shear stress predicted rapid collagen turnover (odds ratio, 1.6 [95% CI, 1.0–2.7]). Receiver operating characteristic analysis demonstrated a good test accuracy (area under the curve, 0.798 [95% CI, 0.598–0.998]) for average wall shear stress with a threshold ≥4.9 Pa for rapid collagen turnover. Our data indicate that turnover rates and stability of collagen in human IA are strongly associated with the presence of modifiable risk factors and aneurysmal hemodynamics. Finally, these findings underline the importance of strict risk factor modification in patients with unruptured IA. Future should include more detailed risk factor data to establish a more causal understanding of hemodynamics and the rupture risk of individual IA.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

Correction of quantification errors in pelvic and spinal lesions caused by ignoring higher photon attenuation of bone in [{sup 18}F]NaF PET/MR

Purpose: MR-based attenuation correction (MRAC) in routine clinical whole-body positron emission tomography and magnetic resonance imaging (PET/MRI) is based on tissue type segmentation. Due to lack of MR signal in cortical bone and the varying signal of spongeous bone, standard whole-body segmentation-based MRAC ignores the higher attenuation of bone compared to the one of soft tissue (MRAC{sub nobone}). The authors aim to quantify and reduce the bias introduced by MRAC{sub nobone} in the standard uptake value (SUV) of spinal and pelvic lesions in 20 PET/MRI examinations with [{sup 18}F]NaF. Methods: The authors reconstructed 20 PET/MR [{sup 18}F]NaF patient data sets acquired with a Philips Ingenuity TF PET/MRI. The PET raw data were reconstructed with two different attenuation images. First, the authors used the vendor-provided MRAC algorithm that ignores the higher attenuation of bone to reconstruct PET{sub nobone}. Second, the authors used a threshold-based algorithm developed in their group to automatically segment bone structures in the [{sup 18}F]NaF PET images. Subsequently, an attenuation coefficient of 0.11 cm{sup −1} was assigned to the segmented bone regions in the MRI-based attenuation image (MRAC{sub bone}) which was used to reconstruct PET{sub bone}. The automatic bone segmentation algorithm was validated in six PET/CT [{sup 18}F]NaF examinations. Relative SUV{sub mean} and SUV{sub max} differences between PET{sub bone} and PET{sub nobone} of 8 pelvic and 41 spinal lesions, and of other regions such as lung, liver, and bladder, were calculated. By varying the assigned bone attenuation coefficient from 0.11 to 0.13 cm{sup −1}, the authors investigated its influence on the reconstructed SUVs of the lesions. Results: The comparison of [{sup 18}F]NaF-based and CT-based bone segmentation in the six PET/CT patients showed a Dice similarity of 0.7 with a true positive rate of 0.72 and a false discovery rate of 0.33. The [{sup 18}F]NaF-based bone segmentation worked well in the pelvis and spine. However, it showed artifacts in the skull and in the extremities. The analysis of the 20 [{sup 18}F]NaF PET/MRI examinations revealed relative SUV{sub max} differences between PET{sub nobone} and PET{sub bone} of (−8.8% ± 2.7%, p = 0.01) and (−8.1% ± 1.9%, p = 2.4 × 10{sup −8}) in pelvic and spinal lesions, respectively. A maximum SUV{sub max} underestimation of −13.7% was found in lesion in the third cervical spine. The averaged SUV{sub mean} differences in volumes of interests in lung, liver, and bladder were below 3%. The average SUV{sub max} differences in pelvic and spinal lesions increased from −9% to −18% and −8% to −17%, respectively, when increasing the assigned bone attenuation coefficient from 0.11 to 0.13 cm{sup −1}. Conclusions: The developed automatic [{sup 18}F]NaF PET-based bone segmentation allows to include higher bone attenuation in whole-body MRAC and thus improves quantification accuracy for pelvic and spinal lesions in [{sup 18}F]NaF PET/MRI examinations. In nonbone structures (e.g., lung, liver, and bladder), MRAC{sub nobone} yields clinically acceptable accuracy.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Electronic structure factors and the importance of adsorbate effects in chemisorption on surface alloys

The chemisorption energy is an integral aspect of surface chemistry, central to numerous fields such as catalysis, corrosion, and nanotechnology. Electronic-structure-based methods such as the Newns-Anderson model are therefore of great importance in guiding the engineering of material surfaces with optimal properties. However, existing methods are inadequate for interpreting complex, multi-metallic systems. Herein, we introduce a physics-based chemisorption model for alloyed transition metal surfaces employing primarily metal d-band properties that accounts for perturbations in both the substrate and adsorbate electronic states upon interaction. Importantly, we show that adsorbate-induced changes in the adsorption site interact with its chemical environment leading to a second-order response in chemisorption energy with the d-filling of the neighboring atoms. We demonstrate the robustness of the model on a wide range of transition metal alloys with O, N, CH, and Li adsorbates yielding a mean absolute error of 0.13 eV versus density functional theory reference chemisorption energies.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Multitiered computational methodology for extracting three-dimensional rotational diffusion coefficients from x-ray photon correlation spectroscopy data without structural information

X-ray photon correlation spectroscopy (XPCS) is a powerful technique for analyzing particle systems by investigating their dynamics in suspensions across a broad range of temporal and spatial scales. This is done by illuminating samples with coherent x-ray beams and calculating the correlation function of the obtained x-ray scattering images. XPCS is uniquely suited for studying Brownian dynamics, consisting of translational and rotational diffusion. While traditional XPCS image analysis techniques can extract translational diffusion components, they are unable to estimate rotational diffusion coefficients. Here, we introduce a methodology that combines the angular-temporal cross-correlation analysis and a algorithmic framework called Multi-Tiered Estimation for Correlation Spectroscopy in 3D for estimating three-dimensional rotational diffusion coefficients from XPCS images of three-dimensional particle systems. We demonstrate our methodology for extracting rotational diffusion coefficients from XPCS data by applying it to simulated noisy x-ray images of systems of crossing nanotubes and proteins that evolve under translational and rotational Brownian motion for different diffusion rates. Furthermore, our results show that our approach determines rotational diffusion coefficients within a few percent error.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Predicting oxidation damage in ultra high-temperature borides: A machine learning approach

Ultra-high temperature (UHT) borides are ceramics materials with melting points above 3000 °C for structural applications in extreme environments. However, at temperatures exceeding 1600 °C and under oxidizing conditions, the material suffers from detrimental degradation. Optimized design and performance of diboride materials under such extreme conditions requires filling the missing composition-microstructure-oxidation gap. This study proposes a computational data-driven framework to connect the processing and microstructure of Ultra-high temperature borides with the oxidation damage. Random Forest Regressor (RFR) model is adopted to forecast the oxide scale thickness developed after oxidation testing based on processing variables and microstructural features. The model trained on a dataset consisting of 107 samples of experimental data extracted from the literature aims to predict oxidation damage. With proper data manipulation and fine model tuning, the predictor could forecast the oxide scale thickness of UHT diborides with a Mean Absolute Error of 37.45 μm and an R-square of 0.83. This model could be used as a high-throughput scheme to design and test new UHT diborides materials computationally. Furthermore, a model with larger composition capabilities could also be developed in the future as more experimental data become available.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics

Ligand binding free energy simulations (LB-FES) that involve sampling of protein functional conformations have been longstanding challenges in research on molecular recognition. Particularly, modeling of the conformational transition pathway and design of the heuristic biasing mechanism are severe bottlenecks for the existing enhanced configurational sampling (ECS) methods. Inspired by the key role of hydration in regulating conformational dynamics of macromolecules, this report proposes a novel ECS approach that facilitates binding-associated structural dynamics by accelerated hydration transitions in combination with the lambda-exchange of free energy perturbation (FEP). Further, two challenging protein-ligand binding processes involving large configurational transitions of the receptor are studied, with hydration transitions at binding sites accelerated by Hamiltonian-simulated annealing of the hydration layer. Without the need for pathway analysis or ad hoc barrier flattening potential, LB-FES were performed with FEP/lambda-exchange molecular dynamics simulation at a minor overhead for annealing of the hydration layer. The LB-FES studies showed that the accelerated rehydration significantly enhances the collective conformational transitions of the receptor, and convergence of binding affinity calculations is obtained at a sweet-spot simulation time scale. Alchemical LB-FES with the proposed ECS strategy is free from the effort of trial and error for the setup and realizes efficient on-the-fly sampling for the collective functional response of the receptor and bound water and therefore presents a practical approach to high-throughput screening in drug discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Improving gas adsorption modeling for MOFs by local calibration of Hubbard U parameters

While computational screening with density functional theory (DFT) is frequently employed for the screening of metal–organic frameworks (MOFs) for gas separation and storage, commonly applied generalized gradient approximations (GGAs) exhibit self-interaction errors, which hinder the predictions of adsorption energies. We investigate the Hubbard U parameter to augment DFT calculations for full periodic MOFs, targeting a more precise modeling of gas molecule–MOF interactions, specifically for N2, CO2, and O2. We introduce a calibration scheme for the U parameter, which is tailored for each MOF, by leveraging higher-level calculations on the secondary building unit (SBU) of the MOF. When applied to the full periodic MOF, the U parameter calibrated against hybrid HSE06 calculations of SBUs successfully reproduces hybrid-quality calculations of the adsorption energy of the periodic MOF. The mean absolute deviation of adsorption energies reduces from 0.13 eV for a standard GGA treatment to 0.06 eV with the calibrated U, demonstrating the utility of the calibration procedure when applied to the full MOF structure. Furthermore, attempting to use coupled cluster singles and doubles with perturbative triples calculations of isolated SBUs for this calibration procedure shows varying degrees of success in predicting the experimental heat of adsorption. It improves accuracy for N2 adsorption for cases of overbinding, whereas its impact on CO2 is minimal, and ambiguities in spin state assignment hinder consistent improvements of O2 adsorption. Our findings emphasize the limitations of cluster models and advocate the use of full periodic MOF systems with a calibrated U parameter, providing a more comprehensive understanding of gas adsorption in MOFs.

Chemistry↗

DOE OSTI2020

Magnetic Navigation for GPS-Denied Airborne Applications

Most current flight systems are dependent on GPS for navigation. Recently, however, navigation in GPS-denied environments has become an area of intensive research. Additional navigation sensor data can be obtained from visual observations (stars or terrain), inertial measurement units, radar, measurements of the local magnetic field, or perhaps even gravity. Absolute and relative positioning via magnetic field measurements have been shown to be viable in many applications including ground navigation, low altitude aircraft flight, and spaceflight. There is greater variability in the magnetic field over shorter distances when flying at low altitude and in ground applications, leading to more accurate positioning. However, ground-based magnetic navigation is often heavily influenced by man-made structures, especially in urban environments. This is not the case for airborne magnetic navigation since the influence of buildings, roads, etc. is negligible for typical aircraft altitudes. For absolute magnetic navigation, the positioning accuracy decreases as altitude increases for a given vehicle velocity, but the observed time variability in the field can be reclaimed by traveling faster through the field. Thus, navigation accuracy becomes a balance of speed and altitude since the higher altitude can be counterbalanced by higher velocity. To understand these effects quantitatively, we explored various techniques to aid a simulated inertial measurement unit with magnetic information. Using a technique known as two-dimensional magnetic map matching, we simulated the performance of airborne magnetic navigation at fixed speed while varying the altitude, flight direction, magnetometer data collection time, reference magnetic map bias error, and type of trajectory (over land or over ocean).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗