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 649 records · Page 36

Nature-based climate solutions can help mitigate the radiative forcing that follows deforestation

Widespread expansion of agriculture and forestry has altered the surface of the Earth, the composition of the atmosphere, and as a result, the climate. Here we quantify the radiative forcing caused by the historical deforestation of an ecoregion in the U.S. Upper Midwest and the adoption of eight nature-based climate solutions. We combined regional forest inventory data with over three decades of remote sensing and in situ data from a replicated land use change experiment. Deforestation of the region caused net global warming (1626 ± 44 µW m -2 ), mainly from the 76 % reduction of ecosystem carbon stocks, but also from the 84 % reduction of the soil methane sink and the 115 % increase in soil nitrous oxide emissions. The associated albedo increase offset 24 % of this greenhouse gas induced warming. For the adoption of nature-based climate solutions, we found that conservation agriculture can provide -39 to -76 ± 31 µW m -2 of climate mitigation over a 100-year time period while short/medium length forestry rotations can provide more at -296 to -881 ± 44 µW m -2 and natural forest regeneration can provide the most at -1555 ± 44 µW m -2 . As the impacts of climate change on nature and society intensify, consideration should be given to the climate mitigation, habitat, and ecosystem services that nature-based climate solutions can provide.

54 ENVIRONMENTAL SCIENCES↗

A coarse-grain force field based on quantum mechanics (CGq FF) for molecular dynamics simulation of poly(ethylene glycol)- block -poly(ε-caprolactone) (PEG- b -PCL) micelles

In order to provide the means to predict from molecular dynamics (MD) simulations the structures of copolymer-based micelles in solution, we developed coarse grain force field (CGq FF) parameters for poly(ethylene glycol) (PEG) and for poly(ε-caprolactone) (PCL). A key advance here is the use of quantum mechanics to train the parameters describing the non-bonded (NB) interactions between the CG beads. The functional forms are the same as the MARTINI CG FF so standard MD codes can be used. Our CGq FF describes well the experimentally observed properties for the polymer–air and polymer–water interfaces, indicating the accuracy of the NB interactions. Additionally, the structural properties (density, radius of gyration ( R g ), and end-to-end distance ( h )) match both experiment and all atom (AA) simulations. We illustrate the application of this CGq FF by following the formation of a spherical micelle from 250 chains of PEG 23 - b -PCL 9 diblock copolymer, each block with molecular weight of 1000 Daltons (10 500 beads, corresponding to 123 250 atoms), in a water box with 119 139 water beads (426 553 water molecules).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrogen bonds and dispersion forces serving as molecular locks for tailored Group 11 bis(amidine) complexes

A flexible polydentate bis(amidine) ligand LH 2 , LH 2 = {CH 2 NH( t Bu)C=N-2-(6-MePy)} 2 , operates as a molecular lock for various coinage metal fragments and forms the dinuclear complexes [LH 2 (MCl) 2 ], M = Cu (1), Au (2), the coordination polymer [{(LH 2 ) 2 (py) 2 (AgCl) 3 }(py) 3 ] n (3), and the dimesityl-digold complex [LH 2 (AuMes) 2 ] (4) by formal insertion of MR fragments (M = Cu, Ag, Au; R = Cl, Mes) into the N–H∙∙∙N hydrogen bonds of LH 2 in yields of 43–95%. Complexes 1, 2, and 4 adopt C 2 -symmetrical structures in the solid state featuring two interconnected 11-membered rings that are locked by two intramolecular N–H∙∙∙R–M hydrogen bonds. QTAIM analyses of the computational geometry-optimized structures 1a, 2a, and 4a reveal 13, 11, and 22 additional bond critical points, respectively, all of which are related to weak intramolecular attractive interactions, predominantly representing dispersion forces, contributing to the conformational stabilization of the C2-symmetrical stereoisomers in the solid state. Variable-temperature 1 H NMR spectroscopy in combination with DFT calculations indicate a dynamic conformational interconversion between two C 2 -symmetrical ground state structures in solution (ΔG ‡ c = 11.1–13.8 kcal mol –1 ), which is accompanied by the formation of an intermediate possessing C i symmetry that retains the hydrogen bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microscopic origin of tunable assembly forces in chiral active environments

Across a variety of spatial scales, from nanoscale biological systems to micron-scale colloidal systems, equilibrium self-assembly is entirely dictated by—and therefore limited by—the thermodynamic properties of the constituent materials. In contrast, nonequilibrium materials, such as self-propelled active matter, expand the possibilities for driving the assemblies that are inaccessible in equilibrium conditions. Recently, a number of works have suggested that active matter drives or accelerates self-organization, but the emergent interactions that arise between solutes immersed in actively driven environments are complex and poorly understood. Here, we analyze and resolve two crucial questions concerning actively driven self-assembly: (i) how, mechanistically, do active environments drive self-assembly of passive solutes? (ii) Under which conditions is this assembly robust? We employ the framework of odd hydrodynamics to theoretically explain numerical and experimental observations that chiral active matter, i.e., particles driven with a directional torque, produces robust and long-ranged assembly forces. Overall, these developments constitute an important step towards a comprehensive theoretical framework for controlling self-assembly in nonequilibrium environments.

36 MATERIALS SCIENCE↗

Closed-form theoretical model of the secondary drop size in partial coalescence—Capturing pertinent timescales and viscous forces

One of the most important outcomes of partial drop coalescence is the ratio of the secondary drop radius to the primary drop radius, known as the drop ratio, ri. This ratio is thought to be approximately constant and independent of physical parameters of the fluids involved. However, here we reveal that ambient fluid viscosity can alter the size of the secondary drop and the drop ratio consequently. Using scaling analysis, we derive a model that predicts the behavior of the drop ratio as a function of the Ohnesorge number, a dimensionless ratio of viscous to inertial forces. In addition, we present our experimental results of coalescing drops on a planar interface under the influence of surface tension gradients. A high-speed digital camera is used to observe the evolution of drops as they coalesce with a bulk liquid. We show that this process is influenced by the surface tension gradient between the drop and the bulk liquid. The ratio of the secondary drop to the primary drop in partial coalescence is smaller than the reported values for coalescence without a surface tension gradient. The analytical model derived through this study is based on a new modified Ohnesorge number that includes surface tension gradients. Our analytical model is compared against other models and the results illustrate good agreement with our experimental findings and experimental data in the literature.

42 ENGINEERING↗

A forced Korteweg–de Vries model for nonlinear mixing of oscillations in a dusty plasma

Nonlinear mixing of oscillations in a dusty plasma due to the harmonic time varying modulation of a nonlinear compressional oscillation is analyzed using a simple mathematical model consisting of a forced Korteweg-de Vries equation. An exact analytic solution of this equation is found to exhibit nonlinear mixing in the system. Here, the model solution can be usefully employed to predict the existence of nonlinear mixing of oscillations in a two-dimensional dusty plasma system of a particular experimental configuration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probing polarization dynamics at specific domain configurations: Computer-vision based automated experiment in piezoresponse force microscopy

Topological defects in ferroelectric materials have attracted much attention due to the emergence of conductive, ferroic, and magnetic functionalities. However, many topological configurations dynamically evolve during the switching processes, making them a challenge to characterize via traditional techniques. Here, we implement an automated experimentation approach for the exploration of functional properties in BiFeO 3 thin films. Specifically, we visualize the ferroelectric domain structures via single frequency piezoresponse force microscopy and implement a computer vision-based algorithm to discover features of interest at which spectroscopic measurements are taken. Subsequently, we employ dimensionality reduction techniques to reveal characteristic polarization behaviors at these features. This approach can be extended to other spectroscopies and modalities to probe only specific features of interest, ultimately enabling dynamical processes in ferroelectrics to be studied.

36 MATERIALS SCIENCE↗

A nearsighted force-training approach to systematically generate training data for the machine learning of large atomic structures

A challenge of atomistic machine-learning (ML) methods is ensuring that the training data are suitable for the system being simulated, which is particularly challenging for systems with large numbers of atoms. Most atomistic ML approaches rely on the nearsightedness principle (“all chemistry is local”), using information about the position of an atom’s neighbors to predict a per-atom energy. Here, in this work, we develop a framework that exploits the nearsighted nature of ML models to systematically produce an appropriate training set for large structures. We use a per-atom uncertainty estimate to identify the most uncertain atoms and extract chunks centered around these atoms. It is crucial that these small chunks are both large enough to satisfy the ML’s nearsighted principle (that is, filling the cutoff radius) and are large enough to be converged with respect to the electronic structure calculation. We present data indicating when the electronic structure calculations are converged with respect to the structure size, which fundamentally limits the accuracy of any nearsighted ML calculator. These new atomic chunks are calculated in electronic structures, and crucially, only a single force—that of the central atom—is added to the growing training set, preventing the noisy and irrelevant information from the piece’s boundary from interfering with ML training. The resulting ML potentials are robust, despite requiring single-point calculations on only small reference structures and never seeing large training structures. We demonstrated our approach via structure optimization of a 260-atom structure and extended the approach to clusters with up to 1415 atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

MBX: A many-body energy and force calculator for data-driven many-body simulations

Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the “many-body energy” (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. In this study, MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum mechanics based non-bonded force field functions for use in molecular dynamics simulations of materials and systems: The nitrogen and oxygen columns

Accurate Force Fields (FFs) are essential for Molecular Dynamics (MD) simulations of the dynamics of realistic materials in terms of atomic-level interactions. The FF parameters of short-range valence interactions can be derived through Quantum Mechanical (QM) calculations on model systems practical for QM (<300 atoms). Similarly, the dynamic electrostatic interactions can be described with methods such as QEq or PQEq that allow charges and polarization to adjust dynamically. However, accurately extracting long-range van der Waals (vdW) interactions from QM calculations poses challenges due to the absence of a definitive method to distinguish between the different energetic components of electrostatics, polarization, vdW, hydrogen bonding, and valence interactions. To do this we use the Perdew–Burke–Ernzerhof flavor of Density Functional Theory, including empirical D3 vdW corrections, to predict the Equation of State for each element (keeping any covalent bonds fixed), from which we obtain the two-body vdW nonbond potential. Here, we extend these calculations to include non-bonded parameters for the N and O columns of the periodic table so that we now describe columns 15 (N), 16 (O), 17 (F), and 18 (Ne) of the periodic table. For these 20 elements, we find that the two-body vdW potentials can all be mapped to a single universal two-body curve, with just three scaling parameters: Re, De, and L. We refer to this as the Universal NonBond (UNB) potential. We expect this to be useful for new MD simulations and a helpful starting point to obtain UNB parameters for the remainder of the periodic table.

Chemistry↗

Compton rocket effect due to the action of radiation reaction force in degenerate plasma

A closed set of fluid equations with radiation reaction force (RRF) are constructed from the moments of the appropriate single particle kinetic equation describing a relativistic degenerate (high density) electron plasma. The closure, in analogy with the Maxwellian closure for non-degenerate plasmas, is affected via a parametrized Fermi-Dirac distribution. It is shown that the degeneracy increases RRF just as will be predicted from the so-called “Compton Rocket” effect.

Physics↗

Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid–base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Potentials of mean force fail to describe chemical bond-breaking in solution

Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗