Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “intermolecular forces”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

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

Self-, N2- and Ar-broadening and line mixing in HCN and C2H2

Self-, N2- and Ar-broadening coefficients were measured for the stretch-bend infrared combination bands nu-1 + nu-1/2 (4004/cm) of HCN and nu-1 + nu-1/5 (4091/cm) of C2H2, using a tunable difference-frequency laser. At atmospheric pressures, the Q branches of these bands exhibit significant rotational narrowing or line mixing. The broadening coefficients are fit with empirical rotationally inelastic collision rate laws, which are then used to model the line mixing in the overlapped Q-branch profiles. Simple energy gap fitting laws appear to be suitable for the shorter-range intermolecular quadrupole-quadrupole and induction forces, whereas an energy-corrected-sudden scaling law works better for the longer-range dipole-dipole and dipole-quadrupole collision partners. In all cases, the line-coupling coefficients are substantially reduced from the rotationally inelastic rates fit to the broadening coefficients, indicating that 35-70 percent of the broadening may be due to other collisional mechanisms such as cross-relaxation to the degenerate H state vibrational level.

Pine, A. S.↗
DOE OSTI2022

Influence of intermolecular interactions on the infrared complex indices of refraction for binary liquid mixtures

This paper investigates the accuracy of deriving the composite optical constants of binary mixtures from only the complex indices of refraction of the neat materials. These optical constants enable the reflectance spectra of the binary mixtures to be modeled for multiple scenarios (e.g., different substrates, thicknesses, volume ratios), which is important for contact and standoff chemical detection. Using volume fractions, each mixture's complex index of refraction was approximated via three different mixing rules. To explore the impact of intermolecular interactions, these predictions are tested by experimental measurements for two representative sets of binary mixtures: (1) tributyl phosphate combined with n-dodecane, a non-polar medium, to represent mixtures which primarily interact via dispersion forces and (2) tributyl phosphate and 1-butanol to represent mixtures with polar functional groups that can also interact via dipole–dipole interactions, including hydrogen bonding. The residuals and the root-mean-square error between the experimental and calculated index values are computed and demonstrate that for miscible liquids in which the average geometry of the cross-interactions can be considered isotropic (e.g., dispersion), the refractive indices of the mixtures can be modeled using composite n and k values derived from volume fractions of the neat liquids. Conversely, in spectral regions where the geometry of the cross-interactions is more restricted and anisotropic (e.g., hydrogen bonding), the calculated n and k values vary from the measured values. The impact of these interactions on the reflectance spectra are then compared by modeling a thin film of the binary mixtures on an aluminum substrate using both the measured and the mathematically computed indices of refraction.

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

Simple and Accurate One-Body Energy and Dipole Moment Surfaces for Water and Beyond

Water is often the testing ground for new, advanced force fields. While advanced functional forms for intermolecular interactions have been integral to the development of accurate water models, less attention has been paid to a transferable model for intramolecular valence terms. In this work, we present a one-body energy and dipole moment surface model, named 1B-UCB, that is simple yet accurate and can be feasibly adapted for both standard and advanced potentials. 1B-UCB for water is comparable in accuracy to those with much more complex functional forms, despite having drastically fewer parameters. The parametrization protocol has been implemented as part of the Q-Force automated workflow and requires only a quantum mechanical Hessian calculation as reference data, hence allowing it to be easily extended to a variety of molecular systems beyond water, which we demonstrate on a selection of small molecules with different symmetries.

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

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ~8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation–π and anion–π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion–π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E int ) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent—these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)- based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation–π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E int values, and maximum errors of ±0.2–0.3 kcal/mol (i.e., ~±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

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

Third-body stabilization of supercritical CO 2 in CO oxidation: development and application of a ReaxFF force field for the CO/O/CO 2 system

Supercritical CO 2 (scCO 2 ) plays a crucial role as a solvent in separation processes, advanced power cycles, and materials processing. Nonetheless, the atomistic comprehension of how the dense scCO 2 matrix influences the fundamental reaction of carbon monoxide (CO) is still insufficiently explored. Experimental studies and molecular dynamics (MD) simulations frequently fail to detect the highly reactive, transient intermediates, such as atomic oxygen (O), that drive these reactions. Here, to address this issue, we have developed a novel ReaxFF reactive force field for the CO 2 /CO/O system. The force field parameters were calibrated using density functional theory and second-order Møller-Plesset calculations to model CO 2 crystal properties, intermolecular interactions, bond dissociation curves, and reaction energy barriers. The force field reproduces the cohesive energy of the CO 2 crystal, the pressure characteristics of bulk scCO 2 , the equation-of-state behavior over a wide pressure–density range, the pressure dependence of the C–O bond length under compression, and the structural properties of liquid and scCO 2 , as documented by experiments, ab-initio MD, and prominent non-reactive models. The force field was subsequently applied to study the CO + O → CO 2 reaction. In a dilute environment, the reaction is inefficient as the newly formed CO 2 rapidly dissociates due to excess kinetic and potential energy acquired from the exothermic reaction. Conversely, in a dense scCO 2 environment, the surrounding matrix acts as an efficient third body, stabilizing the emerging CO2 product via molecular collisions. Statistical analysis confirms an average excess energy dissipation of 133.9 ± 3.6 kcal/mol over 112.4 ± 17.9 ps. Kinetic energy decomposition reveals that ∼ 92% of the excess kinetic energy is stored in internal (rotational and vibrational) degrees of freedom. This ReaxFF force field establishes a mechanistic foundation for third-body stabilization in dense reactive environments.

Chowdhury, Emdadul Haque [Pennsylvania State Univ.↗
DOE OSTI2021

From Intermolecular Interaction Energies and Observable Shifts to Component Contributions and Back Again: A Tale of Variational Energy Decomposition Analysis

Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N 2 , and BF bound to a [Ru(II)(NH 3 ) 5 ] 2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential.

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

Hydrogen Bonding Optimizes Singlet Fission in Carboxylic Acid Functionalized Anthradithiophene Films

The rate of singlet fission, the process of generating two triplet excitons with photoexcitation of one singlet exciton, depends on a combination of singlet/triplet energy balance and intermolecular coupling. In this work, we perform carboxylic acid functionalization of anthradithiophene (ADT) derivatives that result in hydrogen bonds that drive molecular orientation and strong electronic coupling of polycrystalline ADT thin films, leading to ultrafast singlet fission without significant enthalpic driving force. ADT with a single carboxylic acid group exhibits weak intermolecular coupling and slow and inefficient singlet fission, much like the parent ADT, and substitution of different alkylsilyl solubilizing groups has little effect. However, the addition of two carboxylic acid groups on either end of the long axis favors significant coupling and crystallinity in as-deposited thin films that increase the effective singlet fission rate by roughly three orders of magnitude. The properties of the triplet pair, particularly its propensity to form long-lived independent triplets, are also influenced by the degree of long-range intermolecular coupling. The enhancement of intermolecular coupling specific to singlet fission using the ubiquitous cyclic hydrogen bonding motif could impact triplet pair utilization schemes in a variety of contexts.

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

Tautomerism induces bending and twisting of biogenic crystals

Understanding and exploiting material flexibility through phenomena such as the bending and twisting of molecular crystals has been a subject of increased interest owing to the number of applications that benefit from these properties, such as optoelectronics, mechanophotonics, soft robotics, and smart sensors. Here, we report the growth of spontaneously bent and twisted ammonium urate crystals induced by the keto–enol tautomerism of the urate molecule. The major tautomer is native to biogenic crystals, whereas the minor tautomer functions as an effective crystal growth modifier to induce naturally bent and twisted ammonium urate crystals. We show that the degree of curvature can be tailored based on the judicious selection of growth conditions. A combination of state-of-the-art microscopy and spectroscopy techniques are used to characterize the origin of bending. Spatially resolved nano-electron diffraction and high-resolution electron microscopy of naturally bent crystals show nearly single crystallinity with local lattice deformations generated by a combination of screw and edge dislocations. These observations are consistent with photoinduced force microscopy and contact resonance atomic force microscopy, which confirmed spatially resolved changes in the intermolecular interactions and the mechanical properties throughout the cross-sectional and axial regions of bent crystals. A mechanism of bending involving the generation of regionally specific dislocations is proposed as an alternative to more commonly reported models. These findings highlight a unique characteristic of tautomeric crystals that may have broader implications for other biogenic materials.

Science & Technology - Other Topics↗
NASA NTRS2003

Crystallization Physics in Biomacromolecular Systems

The crystals are built of molecules of protein, nucleic acid and their complexes, like viruses, approx. 5x10(exp 3)+ 3x10(exp 6) Da in weight and 2 + 20 nm in effective diameter. This size strongly exceeds action range of molecular forces and makes a big difference with inorganic crystals. Intermolecular contacts form patches on the biomacromolecular surface. Each patch may occupy only a small percent of the whole surface and vary from polymorph to polymorph of the same protein. Thus, under different conditions (pH, solution chemistry, temperature, any area on the macromolecular surface may form a contact. The crystal Young moduli, E approx. equals 0.1 + 0.5 GPa are more than 10 times lower than that of inorganics and the biomolecules themselves. Water within biocrystals (30-70%) is unable to flow unless typical deformation time is longer than approx. 10(exp -5)s. This explains the discrepancy between light scattering and static measurements of E. Nucleation and Growth requires typically concentrations exceeding the equilibrium ones up to 100 times - because of the new size scale results in 10 - 10(exp 3) times lower kinetic coefficients than that needed for inorganic solution growth. All phenomena observed in the latter occur with protein crystallization and are even better studied by AFM. Crystals are typically facetted. Among unexpected findings of general significance are - net molecular exchange flux at kinks is much lower than that expected from supersaturation, steps with low (< approx. 10(exp -2)) kink density at steps follow Gibbs-Thomson law only at very low supersaturations, step segment growth rate may be independent of step energy. Crystal perfection is a must of biocrystallization to achieve the major goal to find 3-D atomic structure of biomacromolecules by x-ray diffraction. Poor diffraction resolution (> 3Angstrom) makes crystallization a bottleneck for structural biology. All defects typical of small molecule crystals are found in biocrystals, but the defects responsible for poor resolution are not identified. Conformational changes are one of them. Biocrystallization in microgravity reportedly results in 20% cases of better crystals. The mechanism of how lack of convection can do this is still not clear. Lower supersaturation, self-purification &om preferentially trapped homologous impurities and step bunching are viable hypotheses.

Chernov, A. A.↗
DOE OSTI2024

How Accurate Are Simulations and Experiments for the Lattice Energies of Molecular Crystals?

Molecular crystals play a central role in a wide range of scientific fields, including pharmaceuticals and organic semiconductor devices. However, they are challenging systems to model accurately with computational approaches because of a delicate interplay of intermolecular interactions such as hydrogen bonding and Van der Waals dispersion forces. Here, by exploiting recent algorithmic developments, we report the first set of diffusion Monte Carlo lattice energies for all 23 molecular crystals in the popular and widely used X23 dataset. Comparisons with previous state-of-the-art lattice energy predictions (on a subset of the dataset) and a careful analysis of experimental sublimation enthalpies reveals that high-accuracy computational methods are now at least as reliable as (computationally derived) experiments for the lattice energies of molecular crystals. Overall, this work demonstrates the feasibility of high-level explicitly correlated electronic structure methods for broad benchmarking studies in complex condensed phase systems, and signposts a route towards closer agreement between experiment and simulation. Published by the American Physical Society 2024

Physics↗
DOE OSTI2023

The j and k dependencies of the line coupling and line mixing effects: Theoretical studies of the relaxation matrices of N 2 -broadened CH 3 D

Line coupling and line mixing effects in parallel and perpendicular bands of CH 3 D perturbed by N 2 have been studied. The work focuses on exhibiting the j and k dependencies of these two processes. The calculations were based on a previously reported anisotropic intermolecular potential including both the long-range multipole, induction, and dispersion forces and a short-range atom-atom model. Throughout this work, it is shown that components with L 1 = 3 of the atom-atom model are dominant contributions to the diffusion operator. As a consequence, in comparison with other molecular systems such as the CH 3 Cl-N 2 and CH 3 I-N 2 , theoretically predicted line coupling and line mixing effects exhibit completely new j and k dependencies. In general, the theoretically calculated halfwidths and intra-doublets’ off-diagonal elements of the relaxation matrix are in reasonable agreement with measurements.

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

Flow resolution and domain of influence in rarefied hypersonic blunt-body flows

The study assesses the effects of the upstream domain size and grid resolution upon flow properties and body aerodynamics computed for rarefied flows over cold blunt bodies with a direct simulation Monte Carlo (DSMC) particle method. Empirical correlations are suggested for aerodynamic coefficients for two-dimensional flows past a perpendicular flat plate. Free-stream parameters which were varied in the study include the Mach number, Knudsen number, surface temperature, and intermolecular potential. Insufficient grid resolution leads to overprediction of aerodynamic heating and forces in the DSMC method. Solution accuracy correlates well with the Reynolds number defined at the wall temperature and the stagnation mean free path relative to the cell dimension. Insufficient upstream domain size in the DSMC method leads to overprediction of heating and drag. Errors in aerodynamic coefficients correlate well with the distance ahead of the body where flow temperature reaches half of its peak value. Simulation of a hard-sphere gas is more sensitive to grid resolution, while simulation of a Maxwell gas is more sensitive to upstream domain size.

Haas, Brian L.↗
NASA NTRS2023

The J and K Dependencies of the Line Coupling and Line Mixing Effects: Theoretical Studies of the Relaxation Matrices of N 2 -Broadened Ch 3 D

Line coupling and line mixing effects in parallel and perpendicular bands of CH 3 D perturbed by N 2 have been studied. The work focuses on exhibiting the j and k dependencies of these two processes. The calculations were based on a previously reported anisotropic intermolecular potential including both the long-range multipole, induction, and dispersion forces and a short-range atom-atom model. It is shown that components with L 1 = 3 of the atom-atom model are dominant contributions to the diffusion operator. As a consequence, in comparison with other molecular systems such as the CH 3 Cl-N 2 and CH 3 I-N 2 , theoretically predicted line coupling and line mixing effects exhibit completely new j and k dependencies. In general, the theoretically calculated halfwidths and intra-doublets’ off-diagonal elements of the relaxation matrix are in reasonable agreement with measurements.

Methane-d↗
NASA NTRS1995

Determination of Quantum Chemistry Based Force Fields for Molecular Dynamics Simulations of Aromatic Polymers

Ab initio quantum chemistry calculations for model molecules can be used to parameterize force fields for molecular dynamics simulations of polymers. Emphasis in our research group is on using quantum chemistry-based force fields for molecular dynamics simulations of organic polymers in the melt and glassy states, but the methodology is applicable to simulations of small molecules, multicomponent systems and solutions. Special attention is paid to deriving reliable descriptions of the non-bonded and electrostatic interactions. Several procedures have been developed for deriving and calibrating these parameters. Our force fields for aromatic polyimide simulations will be described. In this application, the intermolecular interactions are the critical factor in determining many properties of the polymer (including its color).

Jaffe, Richard↗
DOE OSTI2023

Direct Observation of Twin van der Waals Molecular Chains

The van der Waals (vdW) assemblies are the most common structures of materials. However, directly mapping intermolecular electron clouds of a vdW assembly has never been obtained, even though the intramolecular electron clouds were visualized by atomic-resolution techniques. In this report, we unprecedently mapped the intermolecular electron cloud of the assembles of ethanol molecules via ethyl groups with noncontact atomic force microscopy at 5 K, leading to the first visualization of vdW molecular chains, in which ethanol molecules assemble into twin vdW molecular chains in a reverse parallel configuration on Ag(111) plane. Furthermore, spontaneous order-disorder transitions in the chain were dynamically observed, suggesting its unusual properties different from 2D vdW materials. In conclusion, these findings provide an "eye" to see the atomic world of vdW materials.

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

Emergent Above-Gap Photoluminescence in Molecularly Engineered Hybrid Bilayer Crystals

Bilayer crystals, built by stacking two-dimensional (2D) covalent monolayers, give rise to coupled excitonic states whose properties are constrained by fixed lattice symmetry and orientation. Replacing one covalent monolayer with a 2D molecular crystal––held together by non-covalent forces––overcomes this limitation, as molecular functional groups afford tunable in-plane lattice geometry, intermolecular spacing, and interlayer coupling, providing a powerful knob for exciton engineering. Here, in this study, we report four-atom-thick hybrid bilayer crystals (HBCs) synthesized by directly growing single-crystalline PDI molecular crystal atop WS 2 monolayers, which exhibit a robust photoluminescence (PL) peak 120 meV above the WS 2 optical band gap alongside a belowgap emission. Both peaks display strong polarization anisotropy—nearing unity for the above-gap emission—and maintain a perfectly linear power-law dependence up to an excitation density of ~10 7 mW/cm 2 , indicative of coexisting localized and delocalized excitonic states. Substituting PDI with a PTCDA monolayer on WS 2 fully quenches PL, demonstrating molecular control over excitonic emission. Lattice scale ab initio GW and GW-BSE calculations reveal a significantly hybridized bilayer band structure in PDI/WS 2 that supports interlayer excitonic species both above and below the WS 2 gap with strong polarization anisotropy, in excellent agreement with experiment. Our work introduces a molecule-based bilayer platform for the bottom-up design and control of excitonic phenomena in atomically thin optoelectronic and quantum materials.

2D materials↗
DOE OSTI2023

Many-Body Effects in Aqueous Systems: Synergies Between Interaction Analysis Techniques and Force Field Development

Interaction analysis techniques, including the many-body expansion (MBE), symmetry-adapted perturbation theory, and energy decomposition analysis, allow for an intuitive understanding of complex molecular interactions. We review these methods by first providing a historical context for the study of many-body interactions and discussing how nonadditivities emerge from Hamiltonians containing strictly pairwise-additive interactions. We then elaborate on the synergy between these interaction analysis techniques and the development of advanced force fields aimed at accurately reproducing the Born–Oppenheimer potential energy surface. In particular, we focus on ab initio–based force fields that aim to explicitly reproduce many-body terms and are fitted to high-level electronic structure results. These force fields generally incorporate many-body effects through ( a) parameterization of distributed multipoles, ( b) explicit fitting of the MBE, ( c) inclusion of many-atom features in a neural network, and ( d) coarse-graining of many-body terms into an effective two-body term. We also discuss the emerging use of the MBE to improve the accuracy and speed of ab initio molecular dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗