Search NASA⌕ Search

SEARCH · Search NASA

Results for “reactive force fields”

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 181 records · Page 10

Bond order predictions using deep neural networks

Machine learning is an extremely powerful tool for the modern theoretical chemist since it provides a method for bypassing costly algorithms for solving the Schrödinger equation. Already, it has proven able to infer molecular and atomic properties such as charges, enthalpies, dipoles, excited state energies, and others. Most of these machine learning algorithms proceed by inferring properties of individual atoms, even breaking down total molecular energy into individual atomic contributions. In this paper, we introduce a modified version of the Hierarchically Interacting Particle Neural Network (HIP-NN) capable of making predictions on the bonds between atoms rather than on the atoms themselves. We train the modified HIP-NN to infer bond orders for a large number of small organic molecules as computed via the Natural Bond Orbital package. We demonstrate that the trained model is extensible to molecules much larger than those in the training set by studying its performance on the COMP6 dataset. This method has applications in cheminformatics and force field parameterization and opens a promising future for machine learning models to predict other quantities that are defined between atoms such as density matrix elements, Hamiltonian parameters, and molecular reactivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Li + solvation and kinetics of Li + –BF 4 - /PF 6 - ion pairs in ethylene carbonate. A molecular dynamics study with classical rate theories

Using our polarizable force-field models and employing classical rate theories of chemical reactions, we examine the ethylene carbonate (EC) exchange process between the first and second solvation shells around Li+ and the dissociation kinetics of ion pairs Li+-[BF4] and Li+-[PF6] in this solvent. We calculate the exchange rates using transition state theory and correct them with transmission coefficients computed by the reactive flux; Impey, Madden, and McDonald approaches; and Grote-Hynes theory. We found the residence times of EC around Li+ ions varied from 70 to 450 ps, depending on the correction method used. We found the relaxation times changed significantly from Li+-[BF4] to Li+-[PF6] ion pairs in EC. Our results also show that, in addition to affecting the free energy of dissociation in EC, the anion type also significantly influence the dissociation kinetics of ion pairing. This work was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. The calculations were carried out using computer resources provided by the Office of Basic Energy Sciences.

Chang, Tsun-Mei↗

Plasma physics effects on thermonuclear burn rate in the presence of hydrodynamic mix

Hydrodynamic mix can significantly degrade thermonuclear burn rate in an inertial confinement fusion (ICF) target. Successful mitigation requires a detailed understanding of the physical mechanisms by which mix affects burn. Here we summarize the roles of three distinct plasma physics effects on burn rate. The first is the well-known effect of enhanced thermal energy loss from the hot spot and the mitigating role of self-generated or externally-applied magnetic field. The second is the fuel ion separation via inter-species ion diffusion driven by the powerful thermodynamic forces exacerbated by mix during the implosion process. The third is the fusion reactivity modification by fast ion transport in a mix-dominated ICF target, where hot plasma is intermingled with cold fuel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Influence of annealing on the photodeposition of silver on periodically poled lithium niobate

The preferential deposition of metal nanoparticles onto periodically poled lithium niobate surfaces, whereby photogenerated electrons accumulate in accordance with local electric fields and reduce metal ions from solution, is known to depend on the intensity and wavelength of the illumination and the concentration of the solution used. Here, it is shown that for identical deposition conditions (wavelength, intensity, concentration), post-poling annealing for 10 h at 200 °C modifies the surface reactivity through the reorientation of internal defect fields. Whereas silver nanoparticles deposit preferentially on the +z domains on unannealed crystals, the deposition occurs preferentially along 180° domain walls for annealed crystals. In neither case is the deposition selective; limited deposition occurs also on the unannealed –z domain surface and on both annealed domain surfaces. The observed behavior is attributed to a relaxation of the poling-induced defect frustration mediated by Li{sup +} ion mobility during annealing, which affects the accumulation of electrons, thereby changing the surface reactivity. The evolution of the defect field with temperature is corroborated using Raman spectroscopy.

36 MATERIALS SCIENCE↗

Enabling Technologies for High-Throughput Screening of Nano-Porous Materials: Collaboration with the Nanoporous Materials Genome Center

The overarching goal of this research was to develop new methodologies to enable the accurate and efficient modeling of complex materials using computer simulations. Using inter-molecular interaction energies calculated via an accurate but computationally expensive approach (symmetry-adapted perturbation theory), we parameterized efficient next-generation “force fields” to utilize in subsequent simulations. Since the resulting force fields incorporate much of the relevant physics of inter-molecular interactions, they consequently exhibit high transferability from one material to another. This transferability enables the modeling of a wide range of novel materials without additional computational cost. While this approach is quite general, a particular emphasis of this research involved applications to so-called “metal-organic framework”(MOF) materials relevant to energy-intensive gas separations. We focused specifically on CO 2 /N 2 selectivity, which is a key metric for post combustion CO 2 capture efforts at coal-fired power plants. The gas adsorption capacities and selectivity of the MOFs can be tailored via careful functionalization. We have demonstrated that our force fields exhibit predictive accuracy for a wide variety of functionalized MOFs, thus opening the door for the computational design of “tailored” materials for particular separations. Finally, we have also demonstrated the importance of accounting for the presence of reactive contaminant species when evaluating the performance of MOFs in practical applications.

36 MATERIALS SCIENCE↗

Full-Depth Coadds of the WISE and First-Year NEOWISE-Reactivation Images

The Near Earth Object Wide-field Infrared Survey Explorer (NEOWISE) Reactivation mission released data from its first full year of observations in 2015. This data set includes ~2.5 million exposures in each of W1 and W2, effectively doubling the amount of WISE imaging available at 3.4 μm and 4.6 μm relative to the AllWISE release. In this paper, we have created the first ever full-sky set of coadds combining all publicly available W1 and W2 exposures from both the AllWISE and NEOWISE-Reactivation (NEOWISER) mission phases. We employ an adaptation of the unWISE image coaddition framework, which preserves the native WISE angular resolution and is optimized for forced photometry. By incorporating two additional scans of the entire sky, we not only improve the W1/W2 depths, but also largely eliminate time-dependent artifacts such as off-axis scattered moonlight. We anticipate that our new coadds will have a broad range of applications, including target selection for upcoming spectroscopic cosmology surveys, identification of distant/massive galaxy clusters, and discovery of high-redshift quasars. In particular, our full-depth AllWISE+NEOWISER coadds will be an important input for the Dark Energy Spectroscopic Instrument selection of luminous red galaxy and quasar targets. Our full-depth W1/W2 coadds are already in use within the DECam Legacy Survey (DECaLS) and Mayall z-band Legacy Survey (MzLS) reduction pipelines. Finally, much more work still remains in order to fully leverage NEOWISER imaging for astrophysical applications beyond the solar system.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of Transport in the Lower Tropical Stratosphere in a Global Chemistry and Transport Model

A general circulation model (GCM) relies on various physical parameterizations and provides a solution to the atmospheric equations of motion. A data assimilation system (DAS) combines information from observations with a GCM forecast and produces analyzed meteorological fields that represent the observed atmospheric state. An off-line chemistry and transport model (CTM) can use winds and temperatures from a either a GCM or a DAS. The latter application is in common usage for interpretation of observations from various platforms under the assumption that the DAS transport represents the actual atmospheric transport. Here we compare the transport produced by a DAS with that produced by the particular GCM that is combined with observations to produce the analyzed fields. We focus on transport in the tropics and middle latitudes by comparing the age-of-air inferred from observations of SF6 and CO2 with the age-of-air calculated using GCM fields and DAS fields. We also compare observations of ozone, total reactive nitrogen, and methane with results from the two simulations. These comparisons show that DAS fields produce rapid upward tropical transport and excessive mixing between the tropics and middle latitudes. The unrealistic transport produced by the DAS fields may be due to implicit forcing that is required by the assimilation process when there is bias between the GCM forecast and observations that are combined to produce the analyzed fields. For example, the GCM does not produce a quasi-biennial oscillation (QBO). The QBO is present in the analyzed fields because it is present in the observations, and systematic implicit forcing is required by the DAS. Any systematic bias between observations and the GCM forecast used to produce the DAS analysis is likely to corrupt the transport produced by the analyzed fields. Evaluation of transport in the lower tropical stratosphere in a global chemistry and transport model.

Douglass, Anne R.↗

Interface modulated currents in periodically proton exchanged Mg doped lithium niobate

Conductivity in Mg doped lithium niobate (Mg:LN) plays a key role in the reduction of photorefraction and is therefore widely exploited in optical devices. However, charge transport through Mg:LN and across interfaces such as electrodes also yields potential electronic applications in devices with switchable conductivity states. Furthermore, the introduction of proton exchanged (PE) phases in Mg:LN enhances ionic conductivity, thus providing tailorability of conduction mechanisms and functionality dependent on sample composition. To facilitate the construction and design of such multifunctional electronic devices based on periodically PE Mg:LN or similar ferroelectric semiconductors, fundamental understanding of charge transport in these materials, as well as the impact of internal and external interfaces, is essential. In order to gain insight into polarization and interface dependent conductivity due to band bending, UV illumination, and chemical reactivity, wedge shaped samples consisting of polar oriented Mg:LN and PE phases were investigated using conductive atomic force microscopy. In Mg:LN, three conductivity states (on/off/transient) were observed under UV illumination, controllable by the polarity of the sample and the externally applied electric field. Measurements of currents originating from electrochemical reactions at the metal electrode–PE phase interfaces demonstrate a memresistive and rectifying capability of the PE phase. Furthermore, internal interfaces such as domain walls and Mg:LN–PE phase boundaries were found to play a major role in the accumulation of charge carriers due to polarization gradients, which can lead to increased currents. The insight gained from these findings yield the potential for multifunctional applications such as switchable UV sensitive micro- and nanoelectronic devices and bistable memristors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Crack propagation in silica from reactive classical molecular dynamics simulations

Mechanistic insight into the process of crack growth can be obtained through molecular dynamics (MD) simulations. In this investigation of fracture propagation, a slit crack was introduced into an atomistic amorphous silica model and mode I stress was applied through far–field loading until the crack propagates. Atomic displacements and forces and an Irving–Kirkwood method with a Lagrangian kernel estimator were used to calculate the J –integral of classical fracture mechanics around the crack tip. The resulting fracture toughness (K IC ), 0.76 ± 0.16 MPa√m, agrees with experimental values. In addition, the stress fields and dissipation energies around the slit crack indicate the development of an inelastic region ~30Å in diameter. This is one of the first reports of KIC values obtained from up–scaled atomic–level energies and stresses through the J –integral. In conclusion, the application of the ReaxFF classical MD force field in this study provides the basis for future research into crack growth in multicomponent oxides in a variety of environmental conditions.

36 MATERIALS SCIENCE↗

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effect of Domain Size, Boundary, and Loading Conditions on Mechanical Properties of Amorphous Silica: A Reactive Molecular Dynamics Study

Mechanical properties are very important when choosing a material for a specific application. They help to determine the range of usefulness of a material, establish the service life, and classify and identify materials. The size effect on mechanical properties has been well established numerically and experimentally. However, the role of the size effect combined with boundary and loading conditions on mechanical properties remains unknown. In this paper, by using molecular dynamics (MD) simulations with the state-of-the-art ReaxFF force field, we study mechanical properties of amorphous silica (e.g., Young’s modulus, Poisson’s ratio) as a function of domain size, full-/semi-periodic boundary condition, and tensile/compressive loading. We found that the domain-size effect on Young’s modulus and Poisson’s ratio is much more significant in semi-periodic domains compared to full-periodic domains. The results, for the first time, revealed the bimodular and anisotropic nature of amorphous silica at the atomic level. We also defined a “safe zone” regarding the domain size, where the bulk properties of amorphous silica can be reproducible, while the computational cost and accuracy are in balance.

bimodular materials↗

The influence of particle size on the fluid dynamics of a laser-induced plasma

The interaction of a laser-induced shock wave with nanoparticles and microparticles of aluminum oxide is investigated through experiments and modeling. The chemistry and physics of the interaction between the particles and plasma generated from laser ablation shows similarities and discrete differences for the two particle sizes. For both particle sizes, early stage (< 10 μs) ionization was dominant and evidenced by higher concentrations of Al II. While both sizes exhibit ionization over the same duration, the intensity of emission was greater for nanoparticles indicating greater concentrations of ionized species. Moreover, the dispersion of species was notably more elongated for microparticles while radial dispersion was more pronounced for nanoparticles with elevated drag forces. At later stages (i.e., > 10 μs), oxidation reactions were dominant for both particle sizes, but the same distinctions in flow field were observed and attributed to particle drag. In all stages of interaction, microparticles expand axially with less drag that suppresses their radial expansion. As a result, the dispersion of reactive species was mapped over an up to 80% larger area for nanoparticles relative to microparticles. Results shown here can be applied toward advancing experimental diagnostics and particle-shock wave modeling and simulation efforts for energetic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nanoscale wetting controls reactive Pd ensembles in synthesis of dilute PdAu alloy catalysts

The performance of bimetallic dilute alloy catalysts is largely determined by the size of minority metal ensembles on the nanoparticle surface. By analyzing the synthesis of catalysts comprising Pd 8 Au 92 nanoparticles supported on silica using surface-sensitive techniques, we report that whether Pd overgrowth occurs before or after Au nanoparticle deposition onto the support controls the surface Pd ensemble size and abundance. These differences in Pd ensembles influence catalytic reactivity in H 2 –D 2 isotope exchange and benzaldehyde hydrogenation, which, in correlation with theoretical calculations, is used to elucidate the active site(s) in each reaction. To clarify how the synthetic sequence controls the formation of Pd ensembles, we combine numerical wetting calculations and molecular dynamics simulations (with a machine-learned force field) to visualize Pd deposition and migration on the nanoparticle surface, respectively. Our results suggest that the nanoparticle–support interface restricts nanoparticle accessibility to Pd deposition, which consequently controls the Pd ensemble size, illustrating the critical role of nanoscale wetting phenomena during bimetallic catalyst preparation.

36 MATERIALS SCIENCE↗

Reactive optical matter: light-induced motility in electrodynamically asymmetric nanoscale scatterers

From Newton’s third law, which is known as the principle of actio et reactio, we expect the forces between interacting particles to be equal and opposite for closed systems. Otherwise, “nonreciprocal” forces can arise. This has been shown theoretically in the interaction between dissimilar optically trapped particles that are mediated by an external field. As a result, despite the incident external field not having a transverse component of momentum, the particle pair experiences a force in a direction that is transverse to the light propagation direction. In this letter, we directly measure the net nonreciprocal forces in electrodynamically interacting asymmetric nanoparticle dimers and nanoparticle structures that are illuminated by plane waves and confined to pseudo one-dimensional geometries. We show via electrodynamic theory and simulations that interparticle interactions cause asymmetric scattering from heterodimers. Therefore, the nonreciprocal forces are a consequence of momentum conservation. Our study demonstrates that asymmetric scatterers exhibit directed motion due to the breakdown of mirror symmetry in open systems that involve incident external fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Exploring the notion of space coupling propulsion

All existing methods of space propulsion are based on expelling a reaction mass (propellant) to induce motion. Alternatively, 'space coupling propulsion' refers to speculations about reacting with space-time itself to generate propulsive forces. Conceivably, the resulting increases in payload, range, and velocity would constitute a breakthrough in space propulsion. Such speculations are still considered science fiction for a number of reasons: (1) it appears to violate conservation of momentum; (2) no reactive media appear to exist in space; (3) no 'Grand Uniform Theories' exist to link gravity, an acceleration field, to other phenomena of nature such as electrodynamics. The rationale behind these objectives is the focus of interest. Various methods to either satisfy or explore these issues are presented along with secondary considerations. It is found that it may be useful to consider alternative conventions of science to further explore speculations of space coupling propulsion.

Millis, Marc G.↗

Modular Approach for the Synthesis of Bottlebrush Diblock Copolymers from Poly(Glycidyl Methacrylate)-block-Poly(Vinyldimethylazlactone) Backbones

A strategy that enables the facile synthesis of bottlebrush block copolymers with flexible backbones was developed. A demonstration of the strategy’s utility was carried out by grafting end-functionalized polymethylmethacrylate (PMMA) and polystyrene (PS) to the dually reactive block copolymer, poly(glycidyl methacrylate)-block-poly(vinyldimethylazlactone) (PGMA-b-PVDMA). Five different bottlebrush diblock copolymers were investigated by size-exclusion chromatography (SEC), 1H NMR, Fourier transform infrared (FT-IR), differential scanning calorimetry (DSC), X-ray scattering methods, atomic force microscopy (AFM), rheology and computational simulations using molecular dynamics (MD), and self-consistent field theory (SCFT). A relationship between the glass transition temperature and the fraction of chain ends was demonstrated by DSC and highlights the potential of this synthetic method to tailor polymer properties. All five samples were found to be in a disordered phase exhibiting multiscale structures revealed by two broad peaks in small-angle X-ray scattering (SAXS) that can be attributed to graft-to-graft and backbone-to-backbone density correlations using MD simulations. The SCFT-based simulations justify the observation of a disordered phase due to its stabilization by the grafts. Additionally, this modular approach can be easily extended to other grafts, including responsive, conducting, and charged polymers with the prerequisite end groups. The versatility and ease of assembling these functional bottlebrushes constitute a powerful “toolbox” method for the rapid and scalable synthesis of novel bottlebrush block copolymers with desired properties.

36 MATERIALS SCIENCE↗

Accurate simulation of surfaces and interfaces of ten FCC metals and steel using Lennard–Jones potentials

Abstract The earlier integration of validated Lennard–Jones (LJ) potentials for 8 fcc metals into materials and biomolecular force fields has advanced multiple research fields, for example, metal–electrolyte interfaces, recognition of biomolecules, colloidal assembly of metal nanostructures, alloys, and catalysis. Here we introduce 12-6 and 9-6 LJ parameters for classical all-atom simulations of 10 further fcc metals (Ac, Ca (α), Ce (γ), Es (β), Fe (γ), Ir, Rh, Sr (α), Th (α), Yb (β)) and stainless steel. The parameters reproduce lattice constants, surface energies, water interfacial energies, and interactions with (bio)organic molecules in 0.1 to 5% agreement with experiment, as well as qualitative mechanical properties under standard conditions. Deviations are reduced up to a factor of one hundred in comparison to earlier Lennard–Jones parameters, embedded atom models, and density functional theory. We also explain a quantitative correlation between atomization energies from experiments and surface energies that supports parameter development. The models are computationally very efficient and applicable to an exponential space of alloys. Compatibility with a wide range of force fields such as the Interface force field (IFF), AMBER, CHARMM, COMPASS, CVFF, DREIDING, OPLS-AA, and PCFF enables reliable simulations of nanostructures up to millions of atoms and microsecond time scales. User-friendly model building and input generation are available in the CHARMM-GUI Nanomaterial Modeler. As a limitation, deviations in mechanical properties vary and are comparable to DFT methods. We discuss the incorporation of reactivity and features of the electronic structure to expand the range of applications and further increase the accuracy.

Kanhaiya, Krishan (ORCID:0000000236222655)↗

Sodium Halide Adsorption and Water Structure at the α-Alumina(0001)/Water Interface

Alumina is one of the most abundant minerals and has a wide range of industrial applications, with catalysis as one of the most important. Of particular relevance for catalysis is the structure of the mineral/water interface. In this work, water structure and sodium halide adsorption at the neutral α-alumina(0001)/water interface are investigated using molecular dynamics simulations. This work demonstrates the accuracy of the chosen model of the alumina/water interface and shows that high charge density monovalent ions, such as Na+ and F–, have a strong affinity for the interface due to the specific pattern of alumina surface OH groups, such that the adsorbed ions displace waters that are hydrogen-bonded to the surface in their absence. A significant portion of the driving force for anion adsorption arises from surface bound Na+, which reverse the intrinsic surface dipole field and drive the accumulation of halides at the interface. The resulting electrolytic interfacial structure reorients water molecules as far as 1 nm from the surface. Although ion adsorption does not alter the global orientation of surface OH groups, it significantly affects their local geometry. This in turn may affect the reactivity of surface groups and thus play a role in chemical processes occurring at the interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗