Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Force on Force”

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 OSTI2021

Significant impact of forcing uncertainty in a large ensemble of climate model simulations

Forcing due to solar and volcanic variability, on the natural side, and greenhouse gas and aerosol emissions, on the anthropogenic side, are the main inputs to climate models. Reliable climate model simulations of past and future climate change depend crucially upon them. Here we analyze large ensembles of simulations using a comprehensive Earth System Model to quantify uncertainties in global climate change attributable to differences in prescribed forcings. The different forcings considered here are those used in the two most recent phases of the Coupled Model Intercomparison Project (CMIP), namely CMIP5 and CMIP6. We show significant differences in simulated global surface air temperature due to volcanic aerosol forcing in the second half of the 19th century and in the early 21st century. The latter arise from small-to-moderate eruptions incorporated in CMIP6 simulations but not in CMIP5 simulations. We also find significant differences in global surface air temperature and Arctic sea ice area due to anthropogenic aerosol forcing in the second half of the 20th century and early 21st century. These differences are as large as those obtained in different versions of an Earth System Model employing identical forcings. In simulations from 2015 to 2100, we find significant differences in the rates of projected global warming arising from CMIP5 and CMIP6 concentration pathways that differ slightly but are equivalent in terms of their nominal radiative forcing levels in 2100. Our results highlight the influence of assumptions about natural and anthropogenic aerosol loadings on carbon budgets, the likelihood of meeting Paris targets, and the equivalence of future forcing scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Active liquid crystals powered by force-sensing DNA-motor clusters

Significance Single-molecule techniques have elucidated how isolated molecular motors generate piconewton forces with unprecedented detail. However, in diverse biological and synthetic settings, force-generating proteins collectively power nonequilibrium dynamics, including continuous large-scale rearrangements and persistent fluid flows. Characterizing motor-generated forces in these dense and dynamical environments remains a challenge. We assembled a reversible DNA-based force-sensing probe that, by an optical readout, reveals the molecular arrangements and the force loads experienced by kinesin motors. These probes provide insight into motor-generated forces that collectively power the unique dynamics of microtubule-based active nematics, a noteworthy example of an internally driven active matter system. DNA-based force probes can be extended to study forces and stresses in various synthetic systems as well as diverse cellular environments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A Feasibility Study on Use of an Accelerometer to Measure the Dynamic Forces in Turning

Metal cutting is a highly dynamic process that generates continuously varying forces. Measurement of these forces is essential to characterizing a cutting process and to fully utilize the machine tool. A force dynamometer is often used to measure these varying forces; however, it is limited by the natural frequency of the sensor and can sometimes be hard to set up on a machine. Therefore, this study aims to estimate the dynamic component of the cutting force using an accelerometer placed directly below the cutting edge. The placement of the sensor helps achieving better signal to noise ratio (S/N). The cutting tool is modelled as a single degree of freedom (SDOF) system relative to the workpiece and dynamic loads are experimentally simulated on it. The measured acceleration is numerically integrated to obtain velocity and deflection, which along with the natural frequency (ω n ) and damping ratio (ζ), is used to predict the dynamic load using the system equation. The dynamic forces tested over a wide range of frequencies, show good agreement with the data from simultaneous dynamometer measurement and the force estimated using the proposed method. The study shows the feasibility of this method in a real cutting scenario and the capability to apply it to a multi DOF system. One major limitation is the inability to capture the static or quasi-static component of the cutting force.

Cutting vibrations↗
DOE OSTI2023

Optimal white-noise stochastic forcing for linear models of turbulent channel flow

In the present study an optimisation problem is formulated to determine the forcing of an eddy-viscosity-based linearised Navier–Stokes model in channel flow at $Re_\tau \approx 5200$ ( $Re_\tau$ is the friction Reynolds number), where the forcing is white-in-time and spatially decorrelated. The objective functional is prescribed such that the forcing drives a response to best match a set of velocity spectra from direct numerical simulation (DNS), as well as remaining sufficiently smooth. Strong quantitative agreement is obtained between the velocity spectra from the linear model with optimal forcing and from DNS, but only qualitative agreement between the Reynolds shear stress co-spectra from the model and DNS. The forcing spectra exhibit a level of self-similarity, associated with the primary peak in the velocity spectra, but they also reveal a non-negligible amount of energy spent in phenomenologically mimicking the non-self-similar part of the velocity spectra associated with energy cascade. By exploiting linearity, the effect of the individual forcing components is assessed and the contributions from the Orr mechanism and the lift-up effect are also identified. Finally, the effect of the strength of the eddy viscosity on the optimisation performance is investigated. The inclusion of the eddy viscosity diffusion operator is shown to be essential in modelling of the near-wall features, while still allowing the forcing of the self-similar primary peak. In particular, reducing the strength of the eddy viscosity results in a considerable increase in the near-wall forcing of wall-parallel components.

Mechanics↗
DOE OSTI2021

Atomic Force Microscopy and Infrared Nanospectroscopy of COVID-19 Spike Protein for the Quantification of Adhesion to Common Surfaces

The COVID-19 pandemic has claimed millions of lives worldwide, sickened many more, and has resulted in severe socio-economic consequences. As society returns to normal, understanding the spread and persistence of SARS CoV-2 on commonplace surfaces can help to mitigate future outbreaks of coronaviruses and other pathogens. We hypothesize that such an understanding can be aided from studying the binding and interaction of viral proteins with non-biological surfaces. Here, we propose a methodology for investigating the adhesion of the SARS CoV-2 spike glycoprotein on common inorganic surfaces such as aluminum-, copper-, iron-, silica-, and ceria-oxides as well as metallic gold. Quantitative adhesion was obtained from the analysis of measured forces at the nanoscale using an atomic force microscope operated under ambient conditions. Without imposing further constraints on the measurement conditions, our preliminary findings suggest that spike glycoproteins interact with similar adhesion forces across the majority of the metal oxides tested with the exception to gold, for which attraction forces ~10 times stronger than all other materials studied were observed. Ferritin, which was used as a reference protein, was found to exhibit similar adhesion forces as SARS CoV-2 spike protein. Further, this study results show that glycoprotein adhesion forces for similar ambient humidity, tip shape, and contact surface are non-specific to the properties of metal oxide surfaces, which are expected to be covered by a thin water film. The findings suggest that at ambient conditions glycoprotein adhesion to metal oxides is primarily controlled by the water capillary forces and they depend on the surface tension of the liquid (water). We discuss further strategies warranted to decipher the intricate nanoscale forces for improved quantification of the adhesion.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Interplay of the forces governing steroid hormone micropollutant adsorption in vertically-aligned carbon nanotube membrane nanopores

Vertically-aligned carbon nanotube (VaCNT) membranes allow water to conduct rapidly at low pressures and open up the possibility for water purification and desalination, although the ultralow viscous stress in hydrophobic and low-tortuosity nanopores prevents surface interactions with contaminants. In this experimental investigation, steroid hormone micropollutant adsorption by VaCNT membranes is quantified and explained via the interplay of the hydrodynamic drag and friction forces acting on the hormone, and the adhesive and repulsive forces between the hormone and the inner carbon nanotube wall. It is concluded that a drag force above 2.2 × 10 —3 pN overcomes the friction force resulting in insignificant adsorption, whereas lowering the drag force from 2.2 × 10 —3 to 4.3 × 10 —4 pN increases the adsorbed mass of hormones from zero to 0.4 ng cm —2 . At a low drag force of 1.6 × 10 —3 pN, the adsorbed mass of four hormones is correlated with the hormone–wall adhesive (van der Waals) force. These findings explain micropollutant adsorption in nanopores via the forces acting on the micropollutant along and perpendicular to the flow, which can be exploited for selectivity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Detectable ship tracks account for just 5% of aerosol indirect forcing from ship emissions

Ship emissions are a major source of aerosols over oceans, affecting both air quality and energy balance of the climate. However, estimates of their climate forcing diverge between studies relying on visible ship-tracks and those based on models. Here we show that forcing due to visible ship-tracks accounts for just 5% of the total forcing over the southeast Atlantic shipping-lane. Most forcing from ship emissions comes from aerosols that do not form detectable ship-tracks. They are only tips of the iceberg. We make three forcing calculations, one bottom-up based on visible ship-tracks, one top-down based on spatial relationships, and a hybrid approach that combines top-down or model estimated cloud droplet number concentration changes and cloud adjustments. Although the forcing based on machine learning detected ship tracks is an order of magnitude greater than prior results using manually detected ship-tracks, it remains only 5% of that inferred by top-down or cloud adjustment based methods for pre-2020 shipping. The top-down and the combined cloud adjustments methods show similar forcing for the post-2020 reduction in ships’ sulfur emission, although the methods have important regional differences in cloud adjustments that need further investigation. Our results reconcile a long-standing discrepancy in the literature and have important implications for aerosol indirect forcing and marine cloud brightening.

Yuan, Tianle [NASA Goddard Space Flight Center (GS↗
DOE OSTI2023

Accurate Hellmann–Feynman forces from density functional calculations with augmented Gaussian basis sets

The Hellmann–Feynman (HF) theorem provides a way to compute forces directly from the electron density, enabling efficient force calculations for large systems through machine learning (ML) models for the electron density. The main issue holding back the general acceptance of the HF approach for atom-centered basis sets is the well-known Pulay force which, if naively discarded, typically constitutes an error upward of 10 eV/Å in forces. In this work, we demonstrate that if a suitably augmented Gaussian basis set is used for density functional calculations, the Pulay force can be suppressed, and HF forces can be computed as accurately as analytical forces with state-of-the-art basis sets, allowing geometry optimization and molecular dynamics to be reliably performed with HF forces. Our results pave a clear path forward for the accurate and efficient simulation of large systems using ML densities and the HF theorem.

Pathak, Shivesh (ORCID:0000000172798195)↗
DOE OSTI2025

Multi-fidelity learning for interatomic potentials: low-level forces and high-level energies are all you need

The promise of machine learning interatomic potentials (MLIPs) has led to an abundance of public quantum mechanical (QM) training datasets. The quality of an MLIP is directly limited by the accuracy of the energies and atomic forces in the training dataset. Unfortunately, most of these datasets are computed with relatively low-accuracy QM methods, e.g. density functional theory with a moderate basis set. Due to the increased computational cost of more accurate QM methods, e.g. coupled-cluster theory with a complete basis set (CBS) extrapolation, most high-accuracy datasets are much smaller and often do not contain atomic forces. The lack of high-accuracy atomic forces is quite troubling, as training with force data greatly improves the stability and quality of the MLIP compared to training to energy alone. Because most datasets are computed with a unique level of theory, traditional single-fidelity (SF) learning is not capable of leveraging the vast amounts of published QM data. In this study, we apply multi-fidelity learning (MFL) to train an MLIP to multiple QM datasets of different levels of accuracy, i.e. levels of fidelity. Specifically, we perform three test cases to demonstrate that MFL with both low-level forces and high-level energies yields an extremely accurate MLIP—far more accurate than a SF MLIP trained solely to high-level energies and almost as accurate as a SF MLIP trained directly to high-level energies and forces. Therefore, MFL greatly alleviates the need for generating large and expensive datasets containing high-accuracy atomic forces and allows for more effective training to existing high-accuracy energy-only datasets. Indeed, low-accuracy atomic forces and high-accuracy energies are all that are needed to achieve a high-accuracy MLIP with MFL.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Wake forces in a background of quadratically coupled mediators

Two particles can exert forces on each other when embedded in a sea of weakly coupled particles. These “wake forces” occur whenever the source and target particles have quadratic interactions with the mediating particles; they are proportional to the ambient energy density and typically have a range of order the characteristic de Broglie wavelength of the background. The effect can be understood as source particles causing a disturbance in the background waves—a wake—which subsequently interacts with the target particles. Wake forces can be mediated by bosons or fermions, can have spin dependence, may be attractive or repulsive, and have a generally anisotropic spatial profile and range that depends on the phase-space distribution of the ambient particles. In this work, I investigate the application of wake forces to dark matter searches, recast existing limits on short-range forces into leading constraints on dark matter with quadratic couplings, and sketch out potential experimental modifications to optimize sensitivity. Wake forces occur in the Standard Model: the presence of the cosmic neutrino background induces a millimeter-range force about 22 orders of magnitude weaker than gravity. Wake forces may also be relevant in condensed-matter and atomic physics. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2021

Influence of Tool Runout on Force Measurement During Internal Void Monitoring for Friction Stir Welding of 6061-T6 Aluminum

Abstract The goal of this research was to examine how altering the amount of friction stir tool eccentricity while controlling the amount of slant in the tool shoulder (drivers of oscillatory process forces) effects the generation of process force transients during sub-surface void interaction. The knowledge gained will help improve the accuracy of force-based void monitoring methods that have the potential to reduce the need for post-weld inspection. Process force transients during sub-surface void formation were examined for multiple tools with varying magnitudes of kinematic runout. The eccentric motion of the tool produced oscillations in the process forces at the tools rotational frequency that became distorted when features (flats) on the tool probe interacted with voided volumes, generating an amplitude in the force signals at three times the tool rotational frequency (for three-flat tools). A larger tool eccentricity generates a larger amplitude in the force signals at the tool’s rotational frequency that holds a larger potential to create a distortion during void interaction. It was determined that once void becomes large enough to produce an interaction that generates an amplitude at the third harmonic larger than 30% of the amplitude at the rotational frequency in a weld with no interaction (amplitude solely at rotational frequency), the trailing edge of the tool shoulder cannot fully consolidate the void, i.e., it will remain in the final weld. Additionally, once the void exceeds a certain size, the amplitudes of the third harmonics saturate at 70% of the amplitude at the rotational frequency during full consolidation. The interaction between the eccentric probe and sub-surface void was isolated by ensuring any geometric imperfection in the shoulder (slant) with respect to the rotational axis was removed. The results suggest that geometric imperfections (eccentricity and slant) with respect to the tool’s rotational axis must be known when developing a void monitoring method from force transients of this nature.

Engineering↗
DOE OSTI2023

SKS Polarization Anomalies Due to the Coriolis Force

The Earth’s Coriolis force has been well–known to impact surface waves and normal modes, which is essential to accurately interpret these waves. However, the Coriolis force on body waves has been assumed to be negligible and mostly ignored. It has been previously shown that the Coriolis force impacts polarizations of shear waves, whereas the wavefronts remain unaffected. We expand on the potential influences of Earth’s Coriolis force on shear–wave polarization measurements by conducting 3D numerical simulations for elastic waves generated by earthquake and explosive sources in a radially symmetric, and 3D mantle and crustal models. The Coriolis force can produce polarization anomalies of mantle shear waves up to 7° and core phases, such as SKS and SKKS, up to 4°. Uncorrected shear–wave polarizations due to the Coriolis force can cause an additional source of error (5°–10° in fast direction, and 0.2–0.3 s delay time depending on the method and seismic phase), inaccurate interpretation of station misalignments, and imprecise estimates of the core–mantle boundary topography. We show how to correct for the Coriolis force on teleseismic shear waves using 1D ray tracing for well–isolated phases. Here, we recommend the use of full waveform simulations to accurately account for earthquake sources parameters, poorly isolated phases that could include interfering phase arrivals within the measurement time window, and the effect of the Coriolis force on the polarizations of shear waves.

58 GEOSCIENCES↗
DOE OSTI2021

ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution using the Community Ice Sheet Model

Abstract. The future retreat rate for marine-based regions of the Antarctic Ice Sheet is one of the largest uncertainties in sea-level projections. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) aims to improve projections and quantify uncertainties by running an ensemble of ice sheet models with atmosphere and ocean forcing derived from global climate models. Here, the Community Ice Sheet Model (CISM) is used to run ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution. Using multiple combinations of sub-ice-shelf melt parameterizations and calibrations, CISM is spun up to steady state over many millennia. During the spin-up, basal friction parameters and basin-scale thermal forcing corrections are adjusted to optimize agreement with the observed ice thickness. The model is then run forward for 550 years, from 1950–2500, applying ocean thermal forcing anomalies from six climate models. In all simulations, the ocean forcing triggers long-term retreat of the West Antarctic Ice Sheet, especially in the Filchner–Ronne and Ross sectors. Mass loss accelerates late in the 21st century and then rises steadily for several centuries without leveling off. The resulting ocean-forced sea-level rise at year 2500 varies from about 150 to 1300 mm, depending on the melt scheme and ocean forcing. Further experiments show relatively high sensitivity to the basal friction law, moderate sensitivity to grid resolution and the prescribed collapse of small ice shelves, and low sensitivity to the stress-balance approximation. The Amundsen sector exhibits threshold behavior, with modest retreat under many parameter settings but complete collapse under some combinations of low basal friction and high thermal forcing anomalies. Large uncertainties remain, as a result of parameterized sub-shelf melt rates, simplified treatments of calving and basal friction, and the lack of ice–ocean coupling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Acoustic force spectroscopy reveals subtle differences in cellulose unbinding behavior of carbohydrate-binding modules

Protein adsorption to solid carbohydrate interfaces is critical to many biological processes, particularly in biomass deconstruction. To engineer more-efficient enzymes for biomass deconstruction into sugars, it is necessary to characterize the complex protein–carbohydrate interfacial interactions. A carbohydrate-binding module (CBM) is often associated with microbial surface-tethered cellulosomes or secreted cellulase enzymes to enhance substrate accessibility. However, it is not well known how CBMs recognize, bind, and dissociate from polysaccharides to facilitate efficient cellulolytic activity, due to the lack of mechanistic understanding and a suitable toolkit to study CBM–substrate interactions. Our work outlines a general approach to study the unbinding behavior of CBMs from polysaccharide surfaces using a highly multiplexed single-molecule force spectroscopy assay. Here, we apply acoustic force spectroscopy (AFS) to probe a Clostridium thermocellum cellulosomal scaffoldin protein (CBM3a) and measure its dissociation from nanocellulose surfaces at physiologically relevant, low force loading rates. An automated microfluidic setup and method for uniform deposition of insoluble polysaccharides on the AFS chip surfaces are demonstrated. The rupture forces of wild-type CBM3a, and its Y67A mutant, unbinding from nanocellulose surfaces suggests distinct multimodal CBM binding conformations, with structural mechanisms further explored using molecular dynamics simulations. Applying classical dynamic force spectroscopy theory, the single-molecule unbinding rate at zero force is extrapolated and found to agree with bulk equilibrium unbinding rates estimated independently using quartz crystal microbalance with dissipation monitoring. However, our results also highlight critical limitations of applying classical theory to explain the highly multivalent binding interactions for cellulose–CBM bond rupture forces exceeding 15 pN.

09 BIOMASS FUELS↗
DOE OSTI2022

Data-model files associated with the manuscript "The Effects of Spatial and Temporal Resolution of Gridded Meteorological Forcing on Watershed Hydrological Responses" (Shuai et al., 2022 HESS)

This data package contains the model inputs and outputs used in "The Effects of Spatial and Temporal Resolution of Gridded Meteorological Forcing on Watershed Hydrological Responses" (Shuai et al., 2022 HESS). The data.zip file contains the data used to drive the model simulations. The model.zip file contains the XML input file for ATS. The notebook.zip file contains the Jupyter notebooks for pre- and post- processing model results. The figures.zip file contains the raw figures associated with the manuscript.Meteorological forcing plays a critical role in accurately simulating the watershed hydrological cycle. With the advancement of high-performance computing and the development of integrated watershed models, simulating the watershed hydrological cycle at high temporal (hourly to daily) and spatial resolution (10s of meters) has become efficient and computationally affordable. These hyperresolution watershed models require high resolution of meteorological forcing as model input to ensure the fidelity and accuracy of simulated responses. In this study, we utilized the Advanced Terrestrial Simulator (ATS), an integrated watershed model, to simulate surface and subsurface flow and land surface processes using unstructured meshes at the Coal Creek Watershed near Crested Butte (Colorado). We compared simulated watershed hydrologic responses including streamflow, and distributed variables such as evapotranspiration, snow water equivalent (SWE), and groundwater table driven by three publicly available, gridded meteorological forcings (GMFs) -- Daily Surface Weather and Climatological Summaries (Daymet), Parameter-elevation Regressions on Independent Slopes Model (PRISM), and North American Land Data Assimilation System (NLDAS). By comparing various spatial resolutions (ranging from 400 m to 4 km) of PRISM, the simulated streamflow only becomes marginally worse when spatial resolution of meteorological forcing is coarsened to 4 km (or 30% of the watershed area). However, the 4 km resolution has much worse performance than finer resolution in spatially distributed variables such as SWE. Using temporally disaggregated PRISM, we compared models forced by different temporal resolutions (hourly to daily), sub-daily resolution preserves the dynamic watershed responses (e.g., diurnal fluctuation of streamflow) that are absent in results forced by daily resolution. Conversely, the simulated streamflow shows better performance using daily resolution compared to that using sub-daily resolution. Our findings suggest that the choice of GMF and its spatiotemporal resolution depends on the quantity of interest and its spatial and temporal scale, which may have important implications on model calibration and watershed management decisions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Seasonal delay of Sahelian rainfall driven by an east–west contrast in radiative forcing in idealized CESM experiments

Recent studies suggest the observed seasonal delay of rainfall over the Sahel is mainly driven by anthropogenic aerosol forcing, which features a robust east–west contrasting changes in recent decades, with negative and positive top-of-the-atmosphere shortwave radiative forcing in South and East Asia (SA&EA) and in North America and Europe (NA&EU), respectively. Their individual effects on the Sahel rainfall annual cycle remain unclear. Here, by designing idealized sensitivity experiments based on Community Earth System Model (CESM1.2), we show that both negative radiative forcing over SA&EA and positive radiative forcing over NA&EU contribute to the delayed phase and enhanced amplitude of Sahel rainfall annual cycle. To understand the underlying physical processes, both the convective quasi-equilibrium (CQE) and atmospheric energetic frameworks are utilized. Both frameworks can well explain the seasonal delay in the rainfall annual cycle over the Sahel under the regional radiative forcings, with the CQE framework showing stronger explanatory power. Based on the insights from both frameworks, we conducted a moisture budget analysis and found that the negative radiative forcing over SA&EA causes an anomalous low-level anticyclonic circulation, which transports more moisture to the Sahel. As the anomalous anticyclone is stronger in summer and fall compared to spring, rainfall is delayed in the Sahel. Further, under the influence of the positive radiative forcing over NA&EU, more warming in the northern hemisphere causes northward transport of moisture to the northern tropics associated with the cross-equatorial transport of energy. The moisture transport increases relativehumidity over the Sahel and enhances the effective atmospheric heat capacity that delays the seasonal rainfall.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Benchmark calculations of radiative forcing by greenhouse gases

Changes in concentrations of greenhouse gases lead to changes in radiative fluxes throughout the atmosphere. The value of this change, the instantaneous radiative forcing, varies across climate models, due partly to differences in the distribution of clouds, humidity, and temperature across models, and partly due to errors introduced by approximate treatments of radiative transfer. This paper describes an experiment within the Radiative Forcing Model Intercomparision Project that uses benchmark calculations made with line–by–line models to identify parameterization error in the representation of absorption and emission by greenhouse gases. Clear–sky instantaneous forcing by greenhouse gases is computed using a set of 100 profiles, selected from a re–analysis of present–day conditions, that represent the global annual mean forcing from pre–industrial times to the present–day with sampling errors of less than 0.01 W m –2 . Six contributing line–by–line models agree in their estimate of this forcing to within 0.025 W m –2 while even recently–developed parameterizations have typical errors four or more times larger, suggesting both that the samples reveal true differences among line–by–line models and that parameterization error will be readily identifiable. Agreement among line–by–line models is better in the longwave than in the shortwave where differing treatments of the water vapor vapor continuum affect estimates of forcing by carbon dioxide and methane. The impacts of clouds on instantaneous radiative forcing are estimated from climate model simulations. Here, the adjustment due to stratospheric temperature change by assuming fixed dynamical heating. Adjustments are large only for ozone and for carbon dioxide, for which stratospheric cooling introduces modest non–linearity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]

Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗