Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy use behavior”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

BEAM CORE: A Flexible Ecosystem for Freight, Demographics, and Vehicle Analysis

The Behavior, Energy, Autonomy, and Mobility Comprehensive Regional Evaluator (BEAM CORE) is an open-source, modular ecosystem of highly refined, agent-based modeling tools developed by Lawrence Berkeley National Laboratory and the National Laboratory of the Rockies. Organizations such as metropolitan planning organizations, agencies, and companies can use BEAM CORE to analyze freight movement and optimize logistics solutions, assess the impacts of emerging freight technologies or e-commerce trends, model dynamic population growth and evolution, and understand the drivers and impacts of electric vehicle adoption across households in a given region. Users can choose from a flexible suite of modeling modules according to their needs and priorities. The modules integrate with travel demand models to support enhanced analysis of diverse scenarios involving freight movement, vehicle technologies, and other key factors relevant to regional planning.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Characterizing laser-heated polymer foams with simultaneous x-ray fluorescence spectroscopy and Thomson scattering at the Matter in Extreme Conditions Endstation at LCLS

Understanding the behavior of polymer foams at high energy density conditions is crucial to advance inertial fusion energy research. Here, we present a new experimental platform designed to measure the thermodynamic state of these materials at megabar pressures. At the Matter in Extreme Conditions Endstation of the Linac Coherent Light Source, we heat samples using an optical, high-intensity, femtosecond laser and dynamically probe them with ultra-short, coherent x-ray pulses of high peak brightness. We perform x-ray Thomson scattering measurements in forward and backward scattering geometries to capture both collective and non-collective electron behavior in the sample. Simultaneously, x-ray fluorescence spectroscopy is used to measure the emission from a mid-Z dopant, providing complementary information on the plasma conditions. By combining these techniques, we obtain temporally resolved temperature measurements of the transient warm dense matter states. Our initial experiments designed to benchmark the platform with carbon samples yielded data resolving the ultrafast response to laser heating with sub-picosecond resolution, measuring plasma temperatures exceeding 50 eV. These findings lay the foundation for precision studies of the dynamic evolution of laser-heated polymer foams.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2014

A Technique for Mapping Characteristic Lengths to Preserve Energy Dissipated via Strain Softening in a Multiscale Analysis

It is often advantageous to account for the microstructure of the material directly using multiscale modeling. For computational tractability, an idealized repeating unit cell (RUC) is used to capture all of the pertinent features of the microstructure. Typically, the RUC is dimensionless and depends only on the relative volume fractions of the different phases in the material. This works well for non-linear and inelastic behavior exhibiting a positive-definite constitutive response. Although, once the material exhibits strain softening, or localization, a mesh objective failure theories, such as smeared fracture theories, nodal and element enrichment theories (XFEM), cohesive elements or virtual crack closure technique (VCCT), can be utilized at the microscale, but the dimensions of the RUC must then be defined. One major challenge in multiscale progressive damage modeling is relating the characteristic lengths across the scales in order to preserve the energy that is dissipated via localization at the microscale. If there is no effort to relate the size of the macroscale element to the microscale RUC, then the energy that is dissipated will remain mesh dependent at the macroscale, even if it is regularized at the microscale. Here, a technique for mapping characteristic lengths across the scales is proposed. The RUC will be modeled using the generalized method of cells (GMC) micromechanics theory, and local failure in the matrix constituent subcells will be modeled using the crack band theory. The subcell characteristic lengths used in the crack band calculations will be mapped to the macroscale finite element in order to regularize the local energy in a manner consistent with the global length scale. Examples will be provided with and without the regularization, and they will be compared to a baseline case where the size and shape of the element and RUC are coincident (ensuring energy is preserved across the scales).

softening↗
DOE OSTI2021

Effect of temperature and aging duration on ethylene propylene diene monomer (EPDM) nonmetallic components used in caustic liquid waste transfer lines

Nonmetallic materials are used in waste transfer lines at the United States Department of Energy’s Hanford Site Tank Farm in Benton County, Washington, USA. During use, the ethylene propylene diene monomer (EPDM) inner hose of the hose-in-hose transfer lines (HIHTLs) are exposed to β and γ radiation, caustic solutions as well as high temperatures and high pressures. Aging behavior of specimens of EPDM HIHTLs and dog-bone shaped specimens were evaluated by exposing to a solution of 25% sodium hydroxide (NaOH) at 38 °C, 54 °C and 77 °C for 6 and 12-months as well as water only at 77 °C for 12-months. Tensile strength and burst pressure of the specimens were characterized and compared with the unexposed (baseline) specimens. Both the tensile strength of the EPDM material and the burst pressure of the HIHTLs significantly decreased with the higher temperatures and longer exposure times. Analyses using scanning electron microscopy with energy dispersive X-ray spectroscopy (SEM-EDS) of the HIHTLs and the dog-bone specimens were conducted. Our results show that the surface degradation increased with increasing exposure temperature. Penetration of NaOH into the EPDM material also increased with increasing temperature and exposure duration.

36 MATERIALS SCIENCE↗
DOE OSTI2022

HIVE™ [SWR-19-36]

The HIVE™ platform is a mobility services simulation platform developed to provide insight on the energy, infrastructure, service, and economic outcomes of various mobility as a service (MaaS) options. The HIVE platform takes a set of spatiotemporal travel origin-destination pairs and simulates the operation of a predefined mobility service fleet, incorporating request pooling, and various operational and charging behaviors. Hive specializes at modeling fleets of automated electric vehicles (AEVs) and can be used to site and size direct current fast charge (DCFC) stations and measure grid impacts of large-scale AEV fleets serving real-world MaaS trip demand (similar to taxis, Uber, Lyft, etc.). Potential outcomes from a Hive simulation include level of service, total vehicle miles traveled (VMT), deadheading (zero passenger) miles, simultaneous and total energy loads, average occupancy, and more. Hive is developed to generalize to new regions and can be customized to handle many scenarios and operating conditions.

Rames, Clement↗
NASA NTRS2008

Space-based Ornithology-Studying Bird Migration and Environmental Change in North America

Natural fluctuations in the availability of critical stopover sites coupled with anthropogenic destruction of wetlands, land-use change, and anticipated losses due to climate change present migratory birds with a formidable challenge. We have developed an individual-based, spatially explicit bird migration model that simulates the migration routes, timing and energy budgets of individual birds under dynamic weather and land surface conditions. Our model incorporates biophysical constraints, individual bird energy status, bird behavior, and flight aerodynamics. We model the speed, direction, and timing of individual birds moving through a user specified Lagrangian grid. The model incorporates environmental properties including wind speed and direction, topography, dynamic hydrologic properties of the landscape, and environmental suitability. The model is driven by important variables estimated from satellite observations of the land surface, by data assimilation products from weather and climate models, and biological field data. We illustrate the use of the model to study the impact of both short- and long-term environmental variatios, e.g. climate, drought, anthropogenic, on migration timing (phenology), spatial pattern, and fitness (survival and reproductive success). We present several theoretical simulations of the spring migration of Pectoral Sandpiper (Calidris melanotos) in North America with emphasis on the Central flyway from the Gulf of Mexico to Alaska.

Smith, James↗
DOE OSTI2022

Adaptive Urban Traffic Signal Control for Multiple Intersections: An LQR Approach

Traffic congestion leads to severe problems especially in urban traffic networks. It increases the chance of accidents, energy waste, and social costs. In order to address these problems, an adaptive linear quadratic regulator (LQR) approach is developed for traffic signal control at multiple intersections in an urban area. The proposed method controls the green time of the traffic signals to reduce traffic congestion and smooth traffic flow. Real-world data from vision-based traffic sensors are used to build the traffic network model, which mimics the real-world traffic behavior. In addition, the proposed control utilizes recursive least square parameter estimation, which is capable of tracking dynamic changes in traffic conditions. Simulation of Urban MObility (SUMO) is used to analyze the efficacy of the proposed method. Results of the simulation show that the proposed method outperforms pretimed control in various aspects.

adaptive LQR control↗
DOE OSTI2021

Sensitivity analysis of in-pile critical heat flux experiments in TREAT for characterization of RIA power-transient effects

A reactivity-initiated accident (RIA) is one type of postulated design basis accident (DBA) that can cause a departure from nucleate boiling (DNB) event in pressurized-water reactors (PWRs). A DNB occurrence and its consequences depend on the thermophysical properties of the fuel components and coolant, characteristics of the transient energy insertion into the fuel rod, and the onset of the critical heat flux (CHF) phenomenon. To leverage the restart of the Transient Reactor Test (TREAT) Facility, an effort is currently underway to better understand the cladding-to-coolant heat transfer mechanisms and the CHF phenomenon under fast-transient irradiation conditions. This paper characterizes the impact of power transients on the thermal-hydraulic behavior of a TREAT Facility reactor heater rodlet CHF experiment to provide the priority of parameters that need to be investigated for an improved CHF model. Sobol sensitivity analysis methods and the Reactor Excursion and Leak Analysis Program (RELAP5-3D) code were used to identify key input parameters on the uncertainty in the prediction of peak outer- and inner-surface temperatures of the heater tube, as well as the time of the DNB event. A series of sensitivity analyses revealed the total energy deposition on the tube and the transient effects of power pulse had large impacts on the maximum temperatures. The CHF multiplier had the largest impact on the time occurrence of CHF. The overall results show the energy deposition rate in the tube is the most influencing factor to the manifestation of CHF and the resulting thermal-hydraulic behaviors of the tube. The multiplier for the CHF, which is interpreted as the predicted CHF value, has the largest Sobol indices for the time of the CHF in all cases, since it directly determines the occurrence of CHF. It is inferred that the uncertainties in the thermal-hydraulic behaviors of fuels increase with respect to the key parameters as the power pulse becomes broader, and an accurate estimation of the energy deposition rate is required to reduce the uncertainty in the evaluation of the integrity of fuel if the CHF is expected to occur near the peak power. Therefore, the outputs are expected to provide rigorous interpretation of ongoing in-pile CHF experiments in the TREAT Facility reactor regarding thermal-hydraulic behavior of the fuel system aiming for a new transient in-pile CHF model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Atomic-Scale Behavior of Radiation-Resistant ZnO under High-Energy Electron Bombardment

Understanding the atomic structure and defect characteristics of ZnO thin films is crucial for optimizing their electronic properties and performance in advanced applications. Here, we investigate the atomic structure and defect characteristics of atomic layer deposition (ALD)-grown ZnO thin films by using aberration-corrected scanning transmission electron microscopy (STEM). Atomic-resolution imaging identifies prevalent stacking faults, dipole disorder, and various grain boundary types, which are believed to influence the electronic properties of ZnO. Additionally, real-time electron beam exposure experiments demonstrate structural transformations, including crystal growth and surface rearrangements. These findings provide insights into the growth mechanisms of ALD ZnO under high-energy electron irradiation conditions, an important finding for the use of polycrystalline ZnO wide bandgap semiconductors in space-like conditions. In conclusion, our results underscore the capability of STEM in directly visualizing and quantifying atomic-scale defects and beam-induced transformations in radiation-resistant ZnO.

Defects↗
DOE OSTI2025

EVs@Scale Next-Gen Profiles - Fleet Utilization 2024

As part of the U.S. Department of Energy’s EVs@Scale initiative, the Next-Gen Profiles (NGP) project provides a comprehensive, data-driven analysis of electric vehicle (EV) and electric vehicle supply equipment (EVSE) operations across real-world fleet deployments. This paper presents findings from the NGP’s Fleet Utilization study, which investigates operational behavior and asset usage across seventeen EV fleets and two EVSE fleets, encompassing a wide range of vehicle types and use cases. Data collected from diverse sources—varying in format and temporal resolution—are first reformatted into a unified structure. From this harmonized dataset, a suite of rigorously defined performance metrics is calculated at an hourly cadence, enabling consistent cross-comparison of charging, routing, and other key operational behaviors. Amid rapidly increasing EV adoption and growing demands for energy-efficient fleet operations, the analysis reveals clear utilization trends—including diurnal and weekly activity cycles, differences in short versus long charging session dependencies, and route-specific energy usage patterns. These findings highlight the need for tailored infrastructure strategies and the deployment of advanced energy management systems, such as Distributed Energy Resource Management Systems (DERMS) and Site Energy Management Systems (SEMS), which can optimize charging schedules and mitigate peak loads. By leveraging anonymized, harmonized datasets and standardized metrics, this study offers critical insights into fleet behavior and performance, providing a foundation to improve operational efficiency, reduce costs, and enable the scalable deployment of electrified transportation.

Wells, Landon↗
DOE OSTI2026

Building Energy Analysis of Manufactured and Multifamily Housing Types in Juneau, Alaska

This report details the results of building energy modeling analysis evaluating the potential energy savings, economic outcomes, and grid-level electricity reduction associated with cold climate air source heat pump (ccASHP) adoption across multifamily and manufactured housing (MMFH) building typologies in the City and Borough of Juneau (CBJ). Three building archetypes were evaluated: multifamily 4-plex apartments, multifamily 8-plex apartments, and manufactured housing units. Building energy models were developed using OpenStudio-HPXML and calibrated to actual utility consumption data and local meteorological data from the Juneau International Airport weather station using an automated calibration tool following the BPI-2400-S-2015 v.2 standard for model calibration. Occupant behaviors present the greatest variability in successful calibrations. Calibrated models were benchmarked against a baseline electric resistance heating condition, with the selected ccASHP modeled as the retrofit condition and typical meteorological year weather data for all results generation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Role of pH and Eh in geothermal systems: Thermodynamic examples and impacts on scaling and corrosion

The global demand for sustainable and clean renewable energy sources is rapidly increasing to meet growing energy demands while curtailing the CO2 emissions from fossil fuels. Harnessing the natural heat or geothermal energy of the Earth's interior can simultaneously address both issues by offering a clean and renewable energy source. Enhanced geothermal systems require the movement and cycling of fluid through a porous rock media. Further, this process requires understanding the chemistry occurring at the fluid-rock and fluid-metal interfaces both in the deep subsurface and surface, where fluid is circulated for energy transfer or storage. In this review paper, we evaluate the current state of research on how two key variables, pH and Eh, control the fluid-rock reactions controlling mineral scaling and the steel corrosion reactions associated with geothermal systems. These variables are frequently used in geochemical modeling to predict chemical reactions in geothermal systems. Prediction of the fluid-rock behavior can be utilized to formulate new strategies for extending the lifespan and efficiency of heat extraction and utilization and maximizing the safety of geothermal systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Ligand-Structure-Dependent Coherent Vibrational Wavepacket Dynamics in Pyrazolate-Bridged Pt(II) Dimers

Bimetallic transition metal complexes have gained increasing attention because of their versatile functions in solar energy conversion and photonics applications arising from inter-metal electronic coupling. In bimetallic platinum (Pt) complexes, electronic communication between the Pt-centered and ligand-centered moieties have been shown to be critical for defining their excited-state dynamic trajectories undergoing either localized ligand centered (LC)/metal-to-ligand-charge-transfer (MLCT) transitions or delocalized metal-metal-to-ligand-charge-transfer (MMLCT) transitions. The branching of the excited-state intersystem crossing (ISC) trajectories are modulated through structural factors that alter the relative energies of the different states. In this study, we investigated the correlation of the structural factors influencing the excited state trajectories. Using femtosecond broadband transient absorption (fs-BBTA) spectroscopy, ultrafast dynamics in the excited state of two select Pt(II) dimers have been mapped out using their coherent vibrational wavepacket signatures in corresponding transient absorption spectra. To examine how the ligand moieties of the Pt(II) dimers influence excited-state dynamics and the coherent vibrational wavepacket behavior, here we carried out comparative studies on two pyrazolate-bridged Pt(II) dimers of the general formula [Pt( t Bu 2 Pz)(N^C)] 2 ( t Bu 2 Pz = 3,5-di-tert-butylpyrazole); N^C = 7,8-benzoquinoline (bzq, 1) or 1-phenylisoquinoline (piq, 2)). We found that photoexcitation into the low energy absorption bands of 1 and 2 respectively induce the formation of 1 MMLCT states from which ultrafast ISC proceeds, resulting in stimulated emission quenching and decoherence of the vibrational wavepacket motions. The results obtained in this study suggest that both energetics and the structural rigidity of the aromatic cyclometalating ligands in 1 and 2 can significantly influence dynamics along the excited state trajectory characterized by dephasing of the coherent oscillations. The collective results provide direct evidence of how ligand structure alters electronic dynamics along excited state trajectories associated with ISC processes, providing insight into using ligand design to steer photochemical processes.

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

EMI Mitigation of a Ćuk-Based Power-Electronic System Using Switching-Sequence-Based Control

Switching-sequence-based control (SBC) laws when designed based on topological switching behavior can have positive effects on slow- and fast-scale dynamics of a power-electronic system (PES). The slow-scale control can encompass fast PES state regulation and tracking, based on predefined objective, while fastscale control can address differential-mode (DM) and commonmode (CM) spectral-peak energy associated with PES switching operation. Such control laws may offer enhanced programmability to conventional PES design where bulky electromagnetic interference (EMI) filters have been traditionally used to reduce EMI of switching power converters to meet EMI regulatory standards. An EMI filter is always a less programmable solution since it is usually designed for the worst-case EMI mitigation and usually overkill for a PES operating under reduced load condition. The control scheme outlined in this article offers EMI mitigation across wide operating regions without compromising PES regulation. Moreover, it does so by use of switching sequences that guarantee the reachability of the PES dynamics using an advanced Lyapunov-function-based approach. SBC is a powerful tool to generate control actions for a PES based on multivariate PES state constraints. Hence, contemporary EMI regulatory standards are used as constraints in the SBC formulation to operate the PES under wide operating regime while autonomously mitigating the EMI levels. The work may be of paramount importance for operating the ultra-fast-transition recent wide-bandgap semiconductor devices like GaN–FET and SiC MOSFET under higher power with increasing switching frequencies, which is usually desirable for increased power density and reduced switching losses. Here, a hardware Cuk–PES operated ´ with GaN–FETs is fabricated and is used for case illustration. It is shown by experimental results how SBC mitigate DM and CM EMI noise of the PES while maintaining regulation even for the higher order nonminimum phase PES, while reducing sensor requirements using state observer derived from the switching model of the PES.

42 ENGINEERING↗
DOE OSTI2026

The solubility and speciation of REE phosphate endmembers (CePO 4 and YPO 4 ) in Cl-rich aqueous fluids from 350 to 450 °C and implications for natural systems

The rare earth elements (REE) are important metals used increasingly in advanced technologies. Within the crust, the elements Ce and Y are commonly more abundant compared to other lanthanides and comprise important end-member constituents of REE-bearing minerals. Specifically, Ce is part of the light (L) REE which have larger ionic radii than the heavy (H) REE, which are grouped together with Y. These differences in ionic radius can lead to important physico-chemical trends within the lanthanide group. Despite a recent increase in experimental and thermodynamic data for the REE at high temperature and pressure, there is still a significant lack of these data at supercritical conditions. In this study we conducted batch-type experiments to measure the solubility of REE phosphates (CePO 4 and YPO 4 ) at varying starting pH (1.5–10), and salinity (0.01–1.4 mol/kg NaCl) at 350 and P sat , and from 400 to 450 °C at 700 bar. Results show that the solubility of Ce (33–0.14 ppb) is generally higher than Y (13–0.13 ppb) and that Ce complexes more strongly with both chloride and hydroxyl ligands compared to Y. The solubilities of both REE phosphates are highly pH-dependent and, to a lesser extent, depend on salinity at the studied conditions. The solubility data from this study were implemented into the GEMSFITS program to optimize the thermodynamic properties of Ce and Y hydroxyl and chloride species. The updated standard partial molal Gibbs energies of formation (Δ f G 0 T,P ) are used within the experimental temperature and pressure range to accurately predict the CePO 4 and YPO 4 solubility and Ce and Y speciation behavior. Based on the updated thermodynamic properties we also provide formation constants (log β n Cl,OH ) for Ce and Y hydroxyl and chloride species. Updated thermodynamic properties are applied to model REE-apatite dissolution and REE mobility based on the Pea Ridge iron oxide apatite deposit in Missouri, USA. The apatite dissolution model replicates natural observations including the replacement of monazite and xenotime after apatite and is an example of the utility of the new thermodynamic constants applied to supercritical crustal fluids. Furthermore, the findings of this study advance the predictive capabilities of geochemical models, our understanding of the behavior of individual REE, and permit modeling the overarching fractionation trends between LREE and HREE in supercritical crustal fluids.

58 GEOSCIENCES↗
DOE OSTI2023

Modeling of fission gas diffusion and release for Gd 2 O 3 doped UO 2

Uranium dioxide (UO 2 ) is the primary nuclear fuel in light water reactors, and its excess neutronic reactivity can be controlled by adding burnable absorbers, such as Gd 2 O 3 . This burnable absorber has a large neutron absorption cross-section, lowering the high reactivity of the reactor's initial fuel load. However, there needs to be more understanding of how added Gd 2 O 3 influences the properties of UO 2 under irradiation. To understand the behavior of defects and fission gas in the UO 2 /Gd 2 O 3 system under irradiation, we use cluster dynamics modeling supported by density functional theory calculations. First, we calculate the formation energies of Gd point and cluster defects, and evaluate the temperature-dependent defect concentrations using the defect formation energies and entropies. We show that Gd is soluble in UO 2 , introducing a negative charge in the system. Using this information, we adapted the cluster dynamics code Centipede to model the influence of Gd on U self-diffusion and Xe diffusion in UO 2 with 10 wt% Gd 2 O 3 . Also, we analyzed the Xe diffusion as a function of Gd 2 O 3 concentration, showing that the Xe diffusivity is decreased, which means that the athermal diffusivity due to electronic stopping persists at higher temperatures. In conclusion, the decrease in Xe diffusion means that more Xe stays in the matrix, decreasing the Xe release, and lowering its influence of fission gas release on the thermomechanical properties of UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A Comparative Tribo-Mechanical Behavior Analysis of Laser Cladded Nitronic 60 Coating Against Wrought Nitronic 60 Alloy

In this study the tribological and mechanical behavior of a laser cladded Nitronic 60 coating on maraging steel is compared to that of a wrought Nitronic 60 alloy. A multi-layer Nitronic 60 coating was deposited using laser directed energy deposition, and its microstructure, porosity, and microhardness were characterized. The microhardness of the coated region in the laser cladded Nitronic 60 varied between 270-300 HV whereas that for wrought form was 230 HV. Controlled tribological tests under motor oil lubricated conditions were conducted at both room and elevated temperatures. The laser cladded Nitronic coating exhibited poorer wear resistance at room temperature but higher wear resistance at elevated temperature as compared to the wrought Nitronic alloy. Severe plastic deformation and fracture was observed on the laser cladded coating after the wear tests at both temperatures. Wear mechanisms were investigated based on the morphological and compositional characterization of the wear tracks.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Large Scale Study of Ligand–Protein Relative Binding Free Energy Calculations: Actionable Predictions from Statistically Robust Protocols

The accurate and reliable prediction of protein–ligand binding affinities can play a central role in the drug discovery process as well as in personalized medicine. Of considerable importance during lead optimization are the alchemical free energy methods that furnish an estimation of relative binding free energies (RBFE) of similar molecules. Recent advances in these methods have increased their speed, accuracy, and precision. This is evident from the increasing number of retrospective as well as prospective studies employing them. However, such methods still have limited applicability in real-world scenarios due to a number of important yet unresolved issues. Here, we report the findings from a large data set comprising over 500 ligand transformations spanning over 300 ligands binding to a diverse set of 14 different protein targets which furnish statistically robust results on the accuracy, precision, and reproducibility of RBFE calculations. We use ensemble-based methods which are the only way to provide reliable uncertainty quantification given that the underlying molecular dynamics is chaotic. These are implemented using TIES (Thermodynamic Integration with Enhanced Sampling). Results achieve chemical accuracy in all cases. Ensemble simulations also furnish information on the statistical distributions of the free energy calculations which exhibit non-normal behavior. We find that the “enhanced sampling” method known as replica exchange with solute tempering degrades RBFE predictions. We also report definitively on numerous associated alchemical factors including the choice of ligand charge method, flexibility in ligand structure, and the size of the alchemical region including the number of atoms involved in transforming one ligand into another. Our findings provide a key set of recommendations that should be adopted for the reliable application of RBFE methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗