Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “intermolecular potentials”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Optimization of the Reax force field for the lithium–oxygen system using a high fidelity charge model

Studies using molecular dynamics (MD) have long struggled to simulate the failure modes of materials, predicting unrealistically high ductility and failing to capture brittle fracture. The primary cause of this shortcoming is an inadequate description of bond breaking. While reactive force fields such as ReaxFF show improvements compared to traditional force fields, the charge models used yield unphysical partial charges, especially during dissociation of ionic bonds. This flaw may be remedied by using the atom-condensed Kohn–Sham density functional theory (DFT) approximated to a second order (ACKS2) charge model for determining partial charges. In this work, we present a new ACKS2-enabled Reax force field for fracture simulations of lithium oxide systems, which was obtained by training against an extensive set of DFT, multireference configuration interaction (MRCI), and MRCI+Q reference data using genetic optimization techniques. This new force field significantly improves the bond breaking behavior, but still cannot fully capture the brittle fracture in MD simulations, suggesting more research is needed to improve simulation of brittle fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A classical model for three-body interactions in aqueous ionic systems

We present a classical induction model to evaluate the three-body ion–water–water (I–W–W) and water–water–water (W–W–W) interactions in aqueous ionic systems. The classical description of the induction energy is based on electrostatic distributed multipoles up to hexadecapole and distributed polarizabilities up to quadrupole–quadrupole on the O and H atoms of water. The monatomic ions were described by a point charge and a dipole–dipole polarizability, while for the polyatomic ions, distributed multipoles up to hexadecapole and distributed polarizabilities up to quadrupole–quadrupole were used. The accuracy of the classical model is benchmarked against an accurate dataset of 936 (I–W–W) and 2184 (W–W–W) three-body terms for 13 different monatomic and polyatomic cation and anion systems. The classical model shows excellent agreement with the reference second order Moller–Plesset and coupled-cluster single double and perturbative triple [CCSD(T)] three-body energies. The Root-Mean-Square-Errors (RMSEs) for monatomic cations, monatomic anions, and polyatomic ions were 0.29, 0.25, and 0.12 kcal/mol, respectively. The corresponding RMSE for 1744 CCSD(T)/aVTZ three-body (W–W–W) energies, used to train MB-pol, was 0.12 kcal/mol. The accuracy of the proposed classical model demonstrates that the three-body term for aqueous ionic systems can be accurately modeled classically. This approach provides a fast, efficient, and as-accurate path toward modeling the three-body term in aqueous ionic systems that is fully transferable across systems with different ions without the need to fit to tens of thousands of ab initio calculations for each ion to extend existing many-body force fields to interactions between water and ions.

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

Proton radiation effects in indium oxide using cascade molecular dynamics simulations

Metal oxide (MO) semiconductors, characterized by their wide band gaps and notable charge transport properties, are promising candidates for electronic applications in extreme environments, including near-Earth space. However, atomistic simulations of radiation–matter interactions in MOs remain challenging due to the limitations of existing interatomic potentials, which often fail to capture both the short-range repulsive forces essential for radiation damage modeling and the long-range electrostatic effects governing defect evolution. In this work, we develop a customized interatomic potential tailored for radiation damage simulations in indium oxide (In 2 O 3 ) as a model system, a representative MO material. Our potential integrates the Ziegler-Biersack-Littmark potential to accurately describe short-range interactions with Buckingham and Coulombic potentials to account for long-range forces. We perform molecular dynamics simulations of low-energy proton irradiation using this custom potential. We employ the primary knock-on atom (PKA) cascade method to study atomic displacements and primary defect formation. Simulations were conducted for 1 keV proton irradiation in four randomly chosen directions, and PKA-driven defect analyses at 5, 10, and 15 keV to examine the effects of direction and energy level on damage generation. Our results provide insight into the impact of irradiation direction and energy level on the cascade evolution and defect formation mechanisms.

Atomistic simulations↗
NASA NTRS1973

Curve crossing for low-energy elastic scattering of He/+/ by Ne.

The perturbation seen in the experimental differential elastic-scattering cross section for the 40-eV He/+/ + Ne system has been attributed to a single crossing of two intermolecular potential-energy curves. A new theoretical treatment of the curve-crossing problem, namely, that of Delos and Thorson, is employed to obtain the crossing probabilities and phases associated with the crossing. These are determined by utilizing ab initio potentials involved in the crossing and are further used in a partial-wave calculation of the cross section, which is compared with our experiment. The origin of the oscillatory structure observed in the differential cross section is discussed in semiclassical terms by defining the problem in terms of two pseudo-deflection-functions. A rainbow effect is shown to be related to a particular feature (a maximum rather than a minimum) of these deflection functions.

Bobbio, S. M.↗
NASA NTRS1974

Short-range interaction between hydrogen molecules

Recent calculations of the ground state energy of a system of four hydrogen atoms are reviewed with the aim of discerning the short-range repulsive part of the interaction potential between two hydrogen molecules. Extended-basis CI calculations which include a diffuse 2p orbital appear to be capable of determining the total interaction energy between to hydrogen molecules for any separation. Consistent results of such calculations suggest that the potential for intermolecular separations in the region from 1 to 2.5 A can now be specified to better than 10% with considerable confidence. Analytic fits to spherical averages of these results are presented.

Mcmahan, A. K.↗
NASA NTRS1988

Thermodynamics of dense molecular hydrogen-helium mixtures at high pressure

Effective intermolecular pair potentials derived from liquid hydrogen and helium shock wave experiments are presently used in Monte Carlo simulations of mixtures of hydrogen and helium, at densities of up to 1.2 g/cu cm. The model interaction Helmholtz free energy derived accurately reproduces both the Monte Carlo calculation results and the experimental data obtained for densities of up to about 0.6 g/cu cm. An equation of state is derived from the free energy expression that could be useful in interior models of such Jovian planets as Saturn.

Marley, Mark S.↗
NASA NTRS2020

Line shape parameters of PH3 transitions: Theoretical studies of self-broadened widths and line mixing effects

Line mixing effects have been calculated in various parallel and perpendicular bands of self-broadened PH3 lines and compared with recent experimental data. The theoretical approach is an extension to symmetric tops with high inversion barrier of the formalism previously developed for NH3 [Q. Ma and C. Boulet, J. Chem. Phys. 144, 224303 (2016)]. The model takes into account the non-diagonality of the scattering operator within the line space as well as, in a correct way, the double degeneracy of the j, k levels when k ≠ 0. Transitions between such levels should be considered as doublets whose components may be coupled by the line mixing process. It has been shown that, at low pressure, the inversion of the experimental data will strongly depend on the splitting between the two components of a doublet. When it is significant, one can measure independently both the width of one component and the intra-doublet coupling matrix element. Otherwise, one can only measure the sum of these two elements. Comparisons with measurements show that the present formalism leads to accurate predictions of the experimental line shapes.

Remote sensing↗
DOE OSTI2024

Multiscale investigation of the mechanism of biomass deconstruction in the dimethyl isosorbide/water Co-solvent pretreatment system

In the context of promoting a circular bioeconomy, the development of green and efficient lignocellulosic biomass pretreatment technologies so as to realize high value-added biomass utilization is of intense interest. We demonstrated the potential of the bio-based green solvent dimethyl isosorbide (DMI) for the fractionation of Eucalyptus biomass with excellent performance. Here, to investigate the mechanisms involved in biomass fractionation, microimaging and microspectroscopic techniques were employed together with molecular dynamics (MD) simulation and COSMO-RS quantum chemical calculations to derive multiscale information. Both the microstructure and regional chemistry of the cell wall vary significantly with the volume ratio of DMI/H 2 O. The strongest effects were found at DMI/H 2 O = 9:1 and showed visible cell wall tearing cracks and cell wall deformation and collapse as well as the lowest values of cell wall thickness and circularity. From the MD simulations, lignin exhibits collapsed-like structure in pure H 2 O with low solvent accessibility surface area (SASA) and radius of gyration (R g ). In contrast, lignin in DMI/H 2 O shows extended structure with high SASA and solvent interactions dominated by van der Waals forces, with maximal contact in the 9:1 (v/v) system. Further, the COSMO-RS calculated sigma (σ-) potential suggests the intermolecular interactions in DMI and DMI/H 2 O co-solvent are weak, leading to stronger interaction with lignin and correspondingly higher lignin dissolution. The radial distribution functions and σ-potential all show that again DMI/H 2 O at 9:1 is an optimal volume ratio for high lignin dissolution. Finally, this study provides a solvent-ratio dependent mechanism for the action of polar aprotic solvents in the deconstruction of biomass.

09 BIOMASS FUELS↗
NASA NTRS1974

Correlation of the rate of ultrasound propagation with other characteristic magnitudes of the liquid state

The application of ultrasonic radiation to the procedure for determining the physical properties of liquids is discussed. The basis for linking the kinetic properties of liquids to the equilibrium properties through the molecular structure and the potential of the intermolecular forces is defined. The test procedure using ultrasonic radiation is described. Results of tests with specified fluids are tabulated.

Stetiu, C.↗
DOE OSTI2025

Collision integrals within the Chapman–Enskog theory for a generalized Lennard-Jones potential

Here we report the values of the collision integrals, needed for the calculation of the macroscopic transport properties such as viscosity (η) and diffusion coefficient (D) of gases within the Chapman–Enskog kinetic gas theory, for a generalized Lennard-Jones potential (gLJ), a more general potential with an adjustable long range 1/r dependence that can describe a wide range of intermolecular interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1980

Energy transfer in NH3-He collisions

The paper presents a new electron gas intermolecular NH3-He potential determined by the Hartree-Fock results of Davis, Boggs, and Mehrotra (1979) with semiempirical long-range induction and dispersion terms. The resulting potential agrees with beam scattering and pressure broadening measurements; it was also found that the quadrupole induced dipole interaction is not as effective as the short range anisotropy. It was concluded that discrepancies remain in the quantitative description of energy transfer in the NH3-He system; discrepancies result from improper consideration of the dependence of cross-sections on degenerate magnetic quantum levels in the experimental investigation.

Green, S.↗
DOE OSTI2022

Molecular interactions in diffusion-controlled aldol condensation with mesoporous silica nanoparticles

The aldol reaction of p-nitrobenzaldehyde in amino-catalyzed mesoporous silica nanoparticles (MSN) has revealed varying catalytic activity with the size of the pores of MSN. The pore size dependence related to the reactivity indicates that the diffusion process is important. A detailed molecular-level analysis for understanding diffusion requires assessment of the noncovalent interactions of the molecular species involved in the aldol reaction with each other, with the solvent, and with key functional groups on the pore surface. Such an analysis is presented here based upon the effective fragment potential (EFP). The EFP method can calculate the intermolecular interactions, decomposed into Coulomb, polarization, dispersion, exchange-repulsion, and charge-transfer interactions. In this study, the potential energy surfaces corresponding to each intermolecular interaction are analyzed for homo- and hetero-dimers with various configurations. The monomers that compose dimers are five molecules such as p-nitrobenzaldehyde, acetone, n-hexane, propylamine, and silanol. The results illustrate that the dispersion interaction is crucial in most dimers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1998

Modeling Intermolecular Interactions in Nanotubes, Fullerenes and Graphite using a New Long-Range Potential

The cohesive energy and compressibility of strands of a single-wall nanotube rope has been computed using a new long-range potential energy function derived from accurate ab initio quantum chemistry calculations of the benzene dimer and calibrated for energetic and mechanical properties of graphite (at pressures up to 12 GPa). We also use this potential to calculate a variety of properties of carbon nanotubes (both single- and multi-wall) and fullerenes. Extensive comparisons are made with previously published potentials.

Jaffe, Richard↗
DOE OSTI2022

Thermal neutron scattering cross sections for amorphous carbon

Carbon materials are commonly found in both nuclear reactors and experimental systems. Various carbon structures occur in nuclear applications ranging from crystalline and nuclear graphite to the amorphous carbon seen in next-generation advanced reactor designs. Amorphous carbon is based on a randomized graphite-like structure and offers the unique ability to disperse impurities throughout the bulk composition. A graphite-like amorphous carbon system was modeled using the classical molecular dynamics (MD) code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). An improved version of the temperature-dependent Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) potential was used to model the carbon-carbon atomic interactions for the temperature at 300 K along with densities 1.60, 1.70, 1.85, and 2.23 g/cm{sup 3}. From the normalized velocity autocorrelation function (VACF), the phonon density of state (DOS) was then calculated as the Fourier transform of the normalized VACF. This DOS was then used as the primary input for the evaluation of the thermal scattering law (TSL, i.e. S(α,β)) and associated neutron thermal scattering cross sections. The TSL was analyzed using the Full Law Analysis Scattering System Hub (FLASSH). The amorphous structure results in shifts of the phonon DOS to lower energy modes than typically displayed for ideal crystalline graphite. This impact on the DOS is directly reflected in the TSL. Furthermore, the typical optical peak at 0.25 eV for the ideal graphite disappears for amorphous carbon, in good agreement with the expected structure. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Thermal Neutron Scattering Cross Sections for Graphitic Amorphous Carbon

Carbon materials are commonly found in both nuclear reactors and experimental systems. Various carbon structures occur in nuclear applications ranging from crystalline and nuclear graphite to the amorphous carbon seen in next-generation advanced reactor designs. Amorphous carbon is based on a randomized graphite-like structure and offers the unique ability to disperse impurities throughout the bulk composition. A graphite-like amorphous carbon system was modeled using the classical molecular dynamics (MD) code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). An improved version of the temperature-dependent Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) potential was used to model the carbon-carbon atomic interactions for the temperature at 300 K along with densities 1.60, 1.70, 1.85, and 2.23 g/cm 3 . From the normalized velocity autocorrelation function (VACF), the phonon density of state (DOS) was then calculated as the Fourier transform of the normalized VACF. This DOS was then used as the primary input for the evaluation of the thermal scattering law (TSL, i.e. S(α,β)) and associated neutron thermal scattering cross sections. The TSL was analyzed using the Full Law Analysis Scattering System Hub (FLASSH). The amorphous structure results in shifts of the phonon DOS to lower energy modes than typically displayed for ideal crystalline graphite. This impact on the DOS is directly reflected in the TSL. Furthermore, the typical features and the optical graphitic peak at 0.25 eV for the ideal graphite DOS disappear for graphite-like amorphous carbon, which shows good agreement with the expected structure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Quantitative intramolecular fission in oligoacenes, materials, and methods of use thereof

The present invention provides soluble, stable singlet fission (SF) compounds, compositions, materials, methods of their use, and methods for their preparation that provide efficient intramolecular singlet fission (iSF) and multiple excitons. The SF compound may be a dimer, an oligomer, or a polymer of polyoligoacenes, where for example, the compound achieves a triplet yield reaching about 200% per absorbed photon. In this system, SF does not depend on intermolecular inter-actions. Instead, SF is an intrinsic property of the molecule and therefore occurs independent of intermolecular interactions. Singlet fission has the potential to significantly improve the photocurrent in single junction solar cells and thus raise the Shockley-Queisser power conversion efficiency limit from about 33% to about 46% or greater. Quantitative SF yield at room temperature has only been observed in crystalline solids or aggregates of higher acenes.

36 MATERIALS SCIENCE↗