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 37 records · Page 2

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↗

Hydrogen Dynamics in Supercritical Water Probed by Neutron Scattering and Computer Simulations

In this work, an investigation of supercritical water is presented combining inelastic and deep inelastic neutron scattering experiments and molecular dynamics simulations based on a machine-learned potential of ab initio quality. The local hydrogen dynamics is investigated at 250 bar and in the temperature range of 553–823 K, covering the evolution from subcritical liquid to supercritical gas-like water. The evolution of libration, bending, and stretching motions in the vibrational density of states is studied, analyzing the spectral features by a mode decomposition. Moreover, the hydrogen nuclear momentum distribution is measured, and its anisotropy is probed experimentally. It is shown that hydrogen bonds survive up to the higher temperatures investigated, and we discuss our results in the framework of the coupling between intramolecular modes and intermolecular librations. Results show that the local potential affecting hydrogen becomes less anisotropic within the molecular plane in the supercritical phase, and we attribute this result to the presence of more distorted hydrogen bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solution–Doped Donor–Acceptor Copolymers Based on Diketopyrrolopyrrole and 3, 3'–Bis (2–(2–(2–Methoxyethoxy) Ethoxy) ethoxy)–2, 2'–Bithiophene Exhibiting Outstanding Thermoelectric Power Factors with p –Dopants

The design of polymeric semiconductors exhibiting high electrical conductivity (σ) and thermoelectric power factor (PF) will be vital for flexible large-area electronics. In this work, four polymers based on diketopyrrolopyrrole (DPP), 2,3-dihydrothieno[3,4-b][1,4]dioxine (EDOT), thieno[3,2-b]thiophene (TT), and 3, 3'-bis (2-(2-(2-methoxyethoxy) ethoxy) ethoxy)-2, 2'-bithiophene (MEET) are investigated as side-chains, with the MEET polymers newly synthesized for this study. These polymers are systematically doped with tetrafluorotetracyanoquinodimethane ( F 4 TCNQ), CF3SO3H, and the synthesized dopant Cp(CN) 3 -(COOMe) 3 , differing in geometry and electron affinity. The DPP-EDOT-based polymer containing MEET as side-chains exhibits the highest conductivity (σ) ≈700 S cm–1 in this series with the acidic dopant (CF 3 SO 3 H). This polymer also shows the lowest oxidation potential by cyclic voltammetry (CV), the strongest intermolecular interactions evidenced by differential scanning calorimetry (DSC), and has the most oxygen-based functionality for possible hydrogen bonding and ionic screening. Other polymers exhibit high σ ≈300–500 S cm–1 and power factor up to 300 µW m –1 K –2 . The mechanism of conductivity is predominantly electronic, as validated by time-dependent conductance studies and transient thermo voltage monitoring over time, including for those doped with the acid. Furthermore, these materials maintain significant thermal stability and air stability over ≈6 weeks. Density functional theory calculations reveal molecular geometries and inform about frontier energy levels. Raman spectroscopy, in conjunction with scanning electron microscopy (SEM-EDS) and x-ray diffraction, provides insight into the solid-state microstructure and degree of phase separation of the doped polymer films. Infrared spectroscopy enables this study to further quantify the degree of charge transfer from polymer to dopant.

36 MATERIALS SCIENCE↗

Gaussian approximation of dispersion potentials for efficient featurization and machine-learning predictions of metal–organic frameworks

Energy-related descriptors in machine learning are a promising strategy to predict adsorption properties of metal–organic frameworks (MOFs) in the low-pressure regime. Interactions between hosts and guests in these systems are typically expressed as a sum of dispersion and electrostatic potentials. The energy landscape of dispersion potentials plays a crucial role in defining Henry’s constants for simple probe molecules in MOFs. To incorporate more information about this energy landscape, we introduce the Gaussian-approximated Lennard-Jones (GALJ) potential, which fits pairwise Lennard-Jones potentials with multiple Gaussians by varying their heights and widths. The GALJ approach is capable of replicating information that can be obtained from the original LJ potentials and enables efficient development of Gaussian integral (GI) descriptors that account for spatial correlations in the dispersion energy environment. GI descriptors would be computationally inconvenient to compute using the usual direct evaluation of the dispersion potential energy surface. Here, we show that these new GI descriptors lead to improvement in ML predictions of Henry’s constants for a diverse set of adsorbates in MOFs compared to previous approaches to this task.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Concordant Mode Approach (CMA): Vibrational Analysis of New and Upgraded Intermolecular Benchmarks for Noncovalent Bonding

The Concordant Mode Approach (CMA) is a novel method that offers tremendous potential for increasing the system size and the level of theory attainable in quantum chemical computations of molecular vibrational frequencies. To investigate the extension of CMA to intermolecular vibrations, computations with coupled cluster singles and doubles with perturbative triples theory [CCSD(T)] using two augmented correlation-consistent polarized-valence triple-ζ basis sets (aug-cc-pVTZ or h-aug-ccpVTZ) were performed on 17 prototypical loosely bound complexes of hydrogen-bonded, dispersion, and mixed character. These Level A results provide new and upgraded benchmarks for noncovalent bonding and a severe test for CMA vibrational analyses. The Level A target frequencies were recovered remarkably well using second-order Møller−Plesset perturbation theory (MP2) with h-aug-cc-pVTZ for generating the underlying (Level B) normal modes of the CMA scheme. Employing this Level B within the lowest-rung CMA-0A method reproduces the 435 benchmark frequencies with a mean absolute error (MAE) of 0.23 cm −1 and a corresponding standard deviation (σ) of 0.84 cm −1 ; strikingly, the corresponding subset of 106 interfragment frequencies exhibits MAE = 0.34 cm −1 and σ = 0.90 cm −1 . Subsequent application of the higher-rung CMA-2A scheme eliminates all outliers and reduces the overall MAE to a minuscule 0.08 cm −1 with the inclusion of only 3.0% of the off-diagonal couplings not accounted for by CMA-0A. Accordingly, the highly efficient CMA methodology proves to be robust even for vibrations on flat potential energy surfaces.

Aromatic compounds↗

Spin‐State and Clustering Effects in Fe‐Complex Negolytes for Near‐Neutral Aqueous Redox Flow Batteries

Cost-effective redox-active materials are essential for advancing redox flow batteries (RFBs). Iron, with its abundance and suitability as a redox couple, is a promising candidate; however, achieving stable and fast redox reactions in aqueous RFBs remains a challenge. Here, this study presents an Fe-based negolyte stabilized by a hexadentate ligand, where Fe–ligand bonds are enhanced through intermolecular interactions. The sulfonate-substituted Fe complex exhibits a formal potential of −0.44 V versus Ag/AgCl and an exceptionally high rate constant of 0.69 cm s −1 . Near-neutral RFBs incorporating 0.5 M Fe complex show excellent cycling stability, with no discernible capacity fading over 300 cycles. This performance is attributed to intermolecular hydrogen bonds that reinforce Fe–ligand coordination and promote the formation of stable trimeric clusters. Operando electrochemical Raman spectroscopy and density functional theory reveal that π-backdonation from Fe(II) to the imino-phenolate moiety further stabilizes the complex after reduction. In contrast, the hydroxyl-substituted complex exhibits inferior stability due to weaker hydrogen bonding and less pronounced π-backdonation. These findings underscore the importance of ligand design and intermolecular interactions in developing cost-effective, high-performance redox-active materials for aqueous RFBs.

36 MATERIALS SCIENCE↗

Diffusion Monte Carlo approaches for studying nuclear quantum effects in fluxional molecules

Abstract Diffusion quantum Monte Carlo (DMC) provides a powerful approach for obtaining the ground state energy and wave function of molecules, ions, and molecular clusters. The approach is uniquely well suited for studies of fluxional molecules, which undergo large amplitude vibrational motions even in their ground state. In contrast to the electronic structure problem, where the wave function must be antisymmetric with respect to exchange of any pair of electrons, the wave function for the ground vibrational state is nodeless. This greatly simplifies the application of DMC for vibrational problems. Because there is not a single potential function that can be used to describe the intramolecular and intermolecular interactions in all molecular systems, most methods that are used to describe nuclear quantum effects rely on a carefully chosen zero‐order description of the molecular vibrations. In contrast, DMC calculations can be performed in Cartesian coordinates, making the DMC algorithm easily transferable between different chemical systems. In this contribution, the theory that underlies DMC will be discussed along with important considerations for performing DMC calculations. Extensions for evaluating vibrationally excited states and molecular properties are also discussed. Insights that can be obtained from DMC calculations are illustrated in the context of the protonated water clusters. This article is categorized under: Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods Theoretical and Physical Chemistry > Spectroscopy

Chemistry↗

Engineering Structural Transitions in a Multilevel Molecular Switch via Intermolecular Coupling

Controlling molecular conformations with atomic precision is essential for advancing molecular functional electronics, as well as our understanding of molecular dynamics. While switching between bistable molecular conformers is common in nature, creating systems with multiple, addressable states remains synthetically challenging. Here, we demonstrate a bottom-up strategy in which intermolecular interactions give rise to multilevel functionality within a simple two-molecule assembly. Using low-temperature scanning tunneling microscopy, we show that a pyrrolidine dimer on Cu(100) exhibits six distinct adsorption conformations, exceeding the four expected from two independent bistable units. This unusual complexity arises from the interplay between intermolecular van der Waals attraction and steric repulsion, which reshapes the potential energy landscape and changes a single high-energy transition into a sequential two-step pathway. Each step is driven by low-energy inelastic electron excitations, achieving a switching efficiency an order of magnitude higher than that of the monomer. Here, by tuning the bias voltage and tip–molecule distance, we achieve deterministic control over multiple stable states, establishing a general design principle for on-demand engineering of collective molecular behavior and energy-efficient multilevel molecular devices.

Molecular interactions↗

Extensins: Self-Assembly, Crosslinking, and the Role of Peroxidases

The extensin (EXT) network is elaborated by the covalent intermolecular crosslinking of EXT glycoprotein monomers, and its proper assembly is important for numerous aspects of basic wall architecture and cellular defense. In this review, we discuss new advances in the secretion of EXT monomers and the molecular drivers of EXT network self-assembly. Many of the functions of EXTs are conferred through covalent crosslinking into the wall, so we also discuss the different types of known intermolecular crosslinks, the enzymes that are involved, as well as the potential for additional crosslinks that are yet to be identified. EXTs also function in wall architecture independent of crosslinking status, and therefore, we explore the role of non-crosslinking EXTs. As EXT crosslinking is upregulated in response to wounding and pathogen infection, we discuss a potential regulatory mechanism to control covalent crosslinking and its relationship to the subcellular localization of the crosslinking enzymes.

54 ENVIRONMENTAL SCIENCES↗

Crystal Nucleation Kinetics and Mechanism: Influence of Interaction Potential

Modulating liquid-to-solid transitions and the resulting crystalline structure for tailored properties is much desired. Colloidal systems are exemplary to this end, but the fundamental knowledge gaps in relating the influence of intermolecular interactions to crystallization behavior continue to hinder progress. In this study, we address this knowledge gap by studying nucleation and growth in systems with modified Lennard-Jones potential. Specifically, we study the commonly used 12-6 potential and a softer 7-6 potential. The thermodynamic state point for the study is chosen such that both systems are investigated at the same level of supercooling and pressure. Under these conditions, we find that the nucleation rate for both systems is comparable. Interestingly, the nucleation pathways and resulting crystal structures are different. In the 12-6 system, nucleation and growth occur predominantly through the FCC structure. Softening the potential alters the critical nucleus composition and introduces two distinct nucleation pathways. One pathway predominantly leads to the nucleus with a body-centered cubic (BCC) structure, while the other favors the face-centered cubic (FCC) arrangement. Our study illustrates that polymorph selection can be achieved through modifications to intermolecular interactions without impacting nucleation kinetics. The results have significant implications in designing approaches for polymorph selection and modulating self-assembly mechanisms.

36 MATERIALS SCIENCE↗

Guest-Host Interactions in Clathrate Hydrates: Benchmark MP2 and CCSD(T)/CBS Binding Energies of CH4, CO2 and H2S in (H2O)20 Cages

We present benchmark binding energies of naturally occurring gas molecules CH4, CO2, and H2S in the small cage, namely the pentagonal dodecahedron (512) (H2O)20, which is one of the constituent cages of the 3 major lattices (structures I, II and H) of clathrate hydrates. These weak interactions require higher levels of electron correlation and converge slowly with increasing basis set to the Complete Basis Set (CBS) limit, necessitating the use of large basis sets up to the augcc- pV5Z and subsequent correction for Basis Set Superposition Error (BSSE). For the host hollow (H2O)20 cages we have identified a most stable isomer with binding energy of -200.8 ± 2.1 kcal/mol at the CCSD(T)/CBS limit (-199.2 ± 0.5 kcal/mol at the MP2/CBS limit). Additionally, we report converged second order Moller-Plesset (MP2) CBS binding energies for the encapsulation of guests in the (H2O)20 cage of -4.3 ± 0.1 for CH4@(H2O)20, -6.6 ± 0.1 for CO2@(H2O)20 and -8.5 ± 0.1 kcal/mol for H2S@(H2O)20, respectively. For CH4@(H2O)20, exhibiting the weakest encapsulation affinity among the three, we report CCSD(T)/aug-cc-pVTZ binding energies and, based on them, a CCSD(T)/CBS estimate of -4.75 ± 0.1 kcal/mol. To the best of our knowledge, the CCSD(T)/aug-cc-pVTZ calculation for CH4@(H2O)20 is the largest one reported to date (168 valence electrons, 1978 basis functions and the correlation of 84 doubly occupied and 1873 virtual orbitals) and required a scalable implementation of the (T) module on 6144 nodes (350208 cores) of the “Cori” supercomputer at the National Energy Research Supercomputing Center (NERSC) for a total execution time of 195 minutes (for the (T) part). These efficient scalable implementations of highly correlated methods offer the capability to obtain long-lasting benchmarks of intermolecular interactions in complex systems. They also provide a path towards parametrizing classical potentials needed to study the dynamical and transport properties in these complex systems as well as assess the accuracy of lower scaling electronic structure methods such as Density Functional Theory (DFT) and MP2 including its spin-biased variants.

Heindel, Joseph↗

R –8 Dispersion Interaction: Derivation and Application to the Effective Fragment Potential Method

The anisotropic and isotropic R –8 dispersion contributions (disp8) are derived and implemented within the framework of the effective fragment potential (EFP) method formulated with imaginary frequency-dependent Cartesian polarizability tensors distributed at the centroids of the localized molecular orbitals (LMOs). Two forms of damping functions, intermolecular overlap-based and Tang–Toennies, are extended for disp8. To obtain LMO polarizability tensors centered at LMO centroids, an origin-shifting transformation is derived and implemented for the dipole–octopole polarizability tensor and the quadrupole–quadrupole polarizability tensor. The analytic gradient is derived and implemented for the isotropic disp8 contribution. Relative to the previously implemented empirical EFP disp8 energy, the isotropic disp8 component of the interaction energy improves the overall agreement of the EFP dispersion energies with the symmetry-adapted perturbation theory (SAPT) benchmarks, reducing the mean absolute errors (MAEs) and mean absolute percentage errors for most of the databases examined in this work. Further, while the anisotropic disp8 can further enhance the accuracy of the EFP dispersion energy and yield smaller MAEs, significantly overbound dispersion energies are predicted by the anisotropic disp8 when the maximum element in the intermolecular overlap matrix is greater than 0.1, possibly due to the breakdown of the approximations made in the EFP dispersion derivation at a short range. For potential energy scan databases, the newly developed EFP dispersion model with isotropic disp8 yields the overall correct curvature and good agreement with SAPT benchmarks around equilibrium and longer but overestimates the dispersion interactions at a short range. While the overlap-based dispersion-damping functions produce better MAEs than Tang–Toennies damping functions, further improvement is needed to better screen the large attractive dispersion energies at a short range (overlap >0.1).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simultaneous electrocatalytic hydrogenation of aldehydes and phenol over carbon-supported metals

Electrocatalytic hydrogenation of furfural, benzaldehyde, and phenol was studied on carbon-supported Rh, Ru, Pd, and Cu catalysts. All metals were active for hydrogenation of the carbonyl group in furfural and benzaldehyde, but only Rh/C was active for phenol hydrogenation. The intrinsic activities for furfural and benzaldehyde conversion were Pd/C < Ru/C ? Rh/C < Cu/C and Ru/C < Rh/C < Pd/C < Cu/C, respectively. While the trend in furfural hydrogenation remained the same in the presence of phenol, the trend in benzaldehyde hydrogenation in the presence of phenol changed to Rh/C < Ru/C < Pd/C < Cu/C. The Faradaic efficiencies for hydrogenation of both furfural (in the 10–40% range) and benzaldehyde (in the 45–100% range) followed the trend: Cu/C < Rh/C < Ru/C < Pd/C. When phenol and an aldehyde were present in the reactant solution, phenol hydrogenation (only observed on Rh) was suppressed. In contrast, enhancements in the rates of carbonyl reduction were observed for both aldehydes in presence of phenol. Such enhancements depend on both the aldehyde and the metal catalyst. The enhancing factors for hydrogenation of furfural and benzaldehyde were in the 1.5–2 and 2–4 ranges, respectively. We conclude that intermolecular interactions between adsorbed phenol and the aldehydes is a general phenomenon that can potentially enhance the rates of carbonyl hydrogenation. This work highlights co-adsorption effects that pose both challenges and advantages for electrocatalytic reduction of organic compounds.

Electrocatalytic hydrogenation, metal catalysts, p↗

A deep potential model with long-range electrostatic interactions

Machine learning models for the potential energy of multi-atomic systems, such as the deep potential (DP) model, make molecular simulations with the accuracy of quantum mechanical density functional theory possible at a cost only moderately higher than that of empirical force fields. However, the majority of these models lack explicit long-range interactions and fail to describe properties that derive from the Coulombic tail of the forces. To overcome this limitation, we extend the DP model by approximating the long-range electrostatic interaction between ions (nuclei + core electrons) and valence electrons with that of distributions of spherical Gaussian charges located at ionic and electronic sites. The latter are rigorously defined in terms of the centers of the maximally localized Wannier distributions, whose dependence on the local atomic environment is modeled accurately by a deep neural network. In the DP long-range (DPLR) model, the electrostatic energy of the Gaussian charge system is added to short-range interactions that are represented as in the standard DP model. The resulting potential energy surface is smooth and possesses analytical forces and virial. Missing effects in the standard DP scheme are recovered, improving on accuracy and predictive power. By including long-range electrostatics, DPLR correctly extrapolates to large systems the potential energy surface learned from quantum mechanical calculations on smaller systems. We illustrate the approach with three examples: the potential energy profile of the water dimer, the free energy of interaction of a water molecule with a liquid water slab, and the phonon dispersion curves of the NaCl crystal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computational Discovery of Intermolecular Singlet Fission Materials Using Many-Body Perturbation Theory

Intermolecular singlet fission (SF) is the conversion of a photogenerated singlet exciton into two triplet excitons residing on different molecules. SF has the potential to enhance the conversion efficiency of solar cells by harvesting two charge carriers from one high-energy photon, whose surplus energy would otherwise be lost to heat. The development of commercial SF-augmented modules is hindered by the limited selection of molecular crystals that exhibit intermolecular SF in the solid state. Computational exploration may accelerate the discovery of new SF materials. The GW approximation and Bethe–Salpeter equation (GW+BSE) within the framework of many-body perturbation theory is the current state-of-the-art method for calculating the excited-state properties of molecular crystals with periodic boundary conditions. In this Review, we discuss the usage of GW+BSE to assess candidate SF materials as well as its combination with low-cost physical or machine learned models in materials discovery workflows. We demonstrate three successful strategies for the discovery of new SF materials: (i) functionalization of known materials to tune their properties, (ii) finding potential polymorphs with improved crystal packing, and (iii) exploring new classes of materials. In addition, three new candidate SF materials are proposed here, which have not been published previously.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Balancing molecular level influences of intermolecular frustrated Lewis pairs (FLP) for successful design of FLP catalysts for hydrogen storage applications

The rise in energy demands and the deleterious environmental issues related to fossil fuels has led to a surge of interest in hydrogen as a "green" alternative. Hydrogen's extraordinary energy density makes it a potential energy and economic "power"-house. Significant research has been dedicated to materials-based hydrogen storage. One area, liquid organic hydrogen carriers (LOHC) is of substantial interest for the reversible transportation of hydrogen from production to end-use facilities. There are challenges associated with this technology including the dependency on precious metal-based catalysts. Recent work in frustrated Lewis pair (FLP) catalysis demonstrates promise for addressing these challenges. Here this review is focused on assessing recent literature on the utilization of intermolecular FLP main group catalysts for improved hydrogenation/dehydrogenation of various substrates including potential LOHC complexes. This review will present an overview of FLPs, highlight potential hydrogen storage applications, and propose areas where knowledge gaps exist that require further investigations.

08 HYDROGEN↗

From biowaste to BioPave: Biological pathways for sequestration of anthropogenic CO 2 and enhancing durability of roadway infrastructures

Biomass-derived bio-oils are emerging as sustainable, low-carbon alternatives for construction materials, particularly in innovative pavement applications, 'BioPave'. Here, this study evaluates bio-oils from various sources, including algae and wastewater biomass, revealing variations in carbon, nitrogen, and sulfur contents, along with viscosity. We focused on their application in outdoor construction, assessing their resistance to thermal and solar radiation, and moisture. The study also examines their interaction with siliceous surfaces and their influence on asphalt adhesion in BioPave systems under moisture conditions. Laboratory experiments and computational modeling demonstrated molecular composition significantly influenced bio-oils' responses to thermal and UV exposure. Bio-oils rich in polar groups showed hardening upon thermal treatment, while those with higher concentrations of saturated aliphatics remained more stable. Furthermore, our findings highlight the significance of dosage control in maximizing beneficial effects of bio-oils in intermolecular interactions at bitumen-aggregate interface in BioPave applications. These findings offer insights into the potential of bio-oils in sustainable construction and the need for dosage control to optimize pavement performance.

59 BASIC BIOLOGICAL SCIENCES↗