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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

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.↗

NASA NTRS1993

Calculation of molecular excitation rates

State-to-state collisional excitation rates for interstellar molecules observed by radio astronomers continue to be required to interpret observed line intensities in terms of local temperatures and densities. A problem of particular interest is collisional excitation of water which is important for modeling the observed interstellar masers. In earlier work supported by a different NASA Grant, excitation of water in collisions with He atoms was studied; after many years of successively more refined calculations that problem now seems to be well understood, and discrepancies with earlier experimental data for related (pressure broadening) phenomena are believed to reflect experimental errors. Because of interstellar abundances, excitation by H2, the dominant interstellar species, is much more important than excitation by He, although it has been argued that rates for excitation by these are similar. Under the current grant theoretical study of this problem has begun which is greatly complicated by the additional degrees of freedom which must be included both in determining the interaction potential and also in the molecular scattering calculation. We have now computed the interaction forces for nearly a thousand molecular geometries and are close to having an acceptable global fit to these points which is necessary for the molecular dynamics calculations. Also, extensive modifications have been made to the molecular scattering code, MOLSCAT. These included coding the rotational basis sets and coupling matrix elements required for collisions of an asymmetric top with a linear rotor. A new method for numerical solution of the coupled equations has been incorporated. Because of the long-ranged nature of the water-hydrogen interaction it is necessary to integrate the equations to rather large intermolecular separations, and the integration methods previously available in MOLSCAT are not ideal for such cases. However, the method used by Alexander in his HIBRIDON code is particularly suited for such cases. We have obtained this code and incorporated that part which solves the coupled differential equations as an option in the MOLSCAT program.

Flynn, George↗

DOE OSTI2021

Effect of chemical substituents attached to the zwitterion cation on dielectric constant

Materials with high dielectric constant, ε s , are desirable in a wide range of applications including energy storage and actuators. Recently, zwitterionic liquids have been reported to have the largest ε s of any liquid and, thus, have the potential to replace inorganic fillers to modulate the material ε s . Although the large ε s for zwitterionic liquids is attributed to their large molecular dipole, the role of chemical substituents attached to the zwitterion cation on ε s is not fully understood, which is necessary to enhance the performance of soft energy materials. Here, we report the impact of zwitterionic liquid cation chemical substituents on ε s (50 < ε s < 300 at room temperature). Dielectric relaxation spectroscopy reveals that molecular reorientation is the main contributor to the high ε s . The low Kirkwood factor g calculated for zwitterionic liquids (e.g., 0.1–0.2) suggests the tendency for the antiparallel zwitterion dipole alignment expected from the strong electrostatic intermolecular interactions. With octyl cation substituents, the g is decreased due to the formation of hydrophobic-rich domains that restrict molecular reorientation under applied electric fields. In contrast, when zwitterion cations are functionalized with ethylene oxide (EO) segments, g increases due to the EO segments interacting with the cations, allowing more zwitterion rotation in response to the applied field. The reported results suggest that high ε s zwitterionic liquids require a large molecular dipole, compositionally homogeneous liquids (e.g., no aggregation), a maximized zwitterion number density, and a high g, which is achievable by incorporating polar chemical substituents onto the zwitterion cations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Quantifying Outer- and Inner-Coordination Sphere Effects Using Uranium Redox Chemistry in Molten Salt Solutions

Defining the relative influence of intramolecular and intermolecular forces is a fundamental problem in chemistry that is difficult to quantify. To address this challenge, we developed a method to evaluate the relative impact of direct chemical bonding in the inner-coordination sphere vs effects from cations in the outer-coordination sphere by comparative analysis of uranium redox reactivity in various molten salts. We observed that outer-coordination sphere cations (M 1+ ) and inner-coordination sphere anions (X 1– ) both affected uranium redox reactivity, with more polarizing M 1+ and larger X 1– favoring uranium in low oxidation states. Changing M 1+ (Li, Na, K) shifted the U IV + e 1– ⇌ U III (U IV/III ) and U III + 3e 1– → U 0 metal potentials by +330 and +240 mV, respectively. Changing X 1– (Cl, Br, I) caused larger shifts of +440 mV for the U IV/III redox potential and +1060 mV for the U 0 metal deposition potential. Using Coulomb’s Law, we correlated these potentials with electrostatic interactions between UIII and the molten salt. This model provided a facile way of predicting redox chemistry within molten salts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2016

Synthesis, Structure, Characterization, and Decomposition of Nickel Dithiocarbamates: Effect of Precursor Structure and Processing Conditions on Solid-State Products

Single-crystal X-ray structures of four nickel dithiocarbamate complexes, the homoleptic mixed-organic bis-dithiocarbamates Ni[S2CN(isopropyl)(benzyl)]2, Ni[S2CN(ethyl)(n-butyl)]2, and Ni[S2CN(phenyl)(benzyl)]2, as well as the heteroleptic mixed-ligand complex NiCl[P(phenyl)3][(S2CN(phenyl)(benzyl)], were determined. Synthetic, spectroscopic, structural, thermal, and sulfide materials studies are discussed in light of prior literature. The spectroscopic results are routine. A slightly distorted square-planar nickel coordination environment was observed for all four complexes. The organic residues adopt conformations to minimize steric interactions. Steric effects also may determine puckering, if any, about the nickel and nitrogen atoms, both of which are planar or nearly so. A trans-influence affects the Ni-S bond distances. Nitrogen atoms interact with the CS2 carbons with a bond order of about 1.5, and the other substituents on nitrogen display transoid conformations. There are no strong intermolecular interactions, consistent with prior observations of the volatility of nickel dithiocarbamate complexes. Thermogravimetric analysis of the homoleptic species under inert atmosphere is consistent with production of 1:1 nickel sulfide phases. Thermolysis of nickel dithiocarbamates under flowing nitrogen produced hexagonal or -NiS as the major phase; thermolysis under flowing forming gas produced millerite (-NiS) at 300 C, godlevskite (Ni9S8) at 325 and 350 C, and heazlewoodite (Ni3S2) at 400 and 450 C. Failure to exclude oxygen results in production of nickel oxide. Nickel sulfide phases produced seem to be primarily influenced by processing conditions, in agreement with prior literature. Nickel dithiocarbamate complexes demonstrate significant promise to serve as single-source precursors to nickel sulfides, a quite interesting family of materials with numerous potential applications.

microscopy↗

DOE OSTI2025

Probing Hydrogen-Bond Anharmonicity in Cocrystals with Colossal Thermal Expansion

Organic cocrystals with sizable anisotropic thermal expansion (TE) are attractive building blocks for thermally responsive soft actuators. Such materials often display intermolecular H-bonding, whose thermal strength fluctuations are frequently overlooked if displacements and/or conformational changes are visible during TE. In this work, we employed variable-temperature single-crystal synchrotron diffraction to investigate the linear expansion coefficients of two organic cocrystals, 1 and 2, which exhibit colossal thermal expansivity exceeding 140 MK –1 . Despite such pronounced TE, solid-state 2 H and 13 C NMR experiments show no evidence of fast internal reorientations or large-amplitude molecular displacement during expansion. Although nondirectional π–π stacking dominates the colossal TE, our results indicate that hydrogen bonds play a non-negligible role in the anisotropic response. Molecular dynamics simulations reveal a temperature-dependent broadening of H-bond length distributions, which arises from the increased thermal sampling of an anharmonic potential energy surface. Our data support an expansion mechanism based on cooperativity in π-π interactions, accompanied by thermal softening and fluctuation-driven stretching of directional hydrogen bonds in these molecular solids. This mechanism underscores the significance of directional interactions in organic cocrystals, which is a key factor in the design of organic thermoresponsive crystalline materials.

Hernández-Morales, Ernesto A. [Univ. Nacional Auto↗

DOE OSTI2023

Y-Type Non-Fullerene Acceptors with Outer Branched Side Chains and Inner Cyclohexane Side Chains for 19.36% Efficiency Polymer Solar Cells

Raising the lowest unoccupied molecular orbital (LUMO) energy level of Y-type non-fullerene acceptors can increase the open-circuit voltage (V oc ) and thus the photovoltaic performance of the current top performing polymer solar cells (PSCs). One of the viable routes is demonstrated by the successful Y6 derivative of L8-BO with the branched alkyl chains at the outer side. This will introduce steric hindrance and reduce intermolecular aggregation, thus open up the bandgap and raise the LUMO energy level. To take further advantages of the steric hindrance influence on optoelectronic properties of Y6 derivatives, two Y-type non-fullerene acceptors of BTP-Cy-4F and BTP-Cy-4Cl are designed and synthesized by adopting outer branched side chains and inner cyclohexane side chains. An outstanding V oc of 0.937 V is achieved in the D18:BTP-Cy-4F binary blend devices along with a power conversion efficiency (PCE) of 18.52%. With the addition of BTP-eC9 to extend the absorption spectral coverage, a remarkable PCE of 19.36% is realized finally in the related ternary blend devices, which is one of the highest values for single-junction PSCs at present. The results illustrate the great potential of cyclohexane side chains in constructing high-performance non-fullerene acceptors and their PSCs.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Enhanced Photo-Cross-Linking of Thymines in DNA Holliday Junction-Templated Squaraine Dimers

Programmable self-assembly of dyes using DNA templates to promote exciton delocalization in dye aggregates is gaining considerable interest. New methods to improve the rigidity of the DNA scaffold and thus the stability of the molecular dye aggregates to encourage exciton delocalization are desired. In these dye–DNA constructs, one potential way to increase the stability of the aggregates is to create an additional covalent bond via photo-cross-linking reactions between thymines in the DNA scaffold. Specifically, we report an approach to increase the yield of photo-cross-linking reaction between thymines in the core of a DNA Holliday junction while limiting the damage from UV irradiation to DNA. We investigated the effect of the distance between thymines on the photo-cross-linking reaction yields by using linkers with different lengths to tether the dyes to the DNA templates. By comprehensively evaluating the photo-cross-linking reaction yields of dye–DNA aggregates using linkers with different lengths, we conclude that interstrand thymines tend to photo-cross-link more efficiently with short linkers. A higher cross-linking yield was achieved due to the shorter intermolecular distance between thymines influenced by strong dye–dye interactions. Furthermore, our method establishes the possibility of improving the stability of DNA-scaffolded dye aggregates, thereby expanding their use in exciton-based applications such as light harvesting, nanoscale computing, quantum computing, and optoelectronics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

A theoretical investigation of the hydrolysis of uranium hexafluoride: the initiation mechanism and vibrational spectroscopy

Depleted uranium hexafluoride (UF 6 ), a stockpiled byproduct of the nuclear fuel cycle, reacts readily with atmospheric humidity, but the mechanism is poorly understood. Here we compare several potential initiation steps at a consistent level of theory, generating underlying structures and vibrational modes using hybrid density functional theory (DFT) and computing relative energies of stationary points with double-hybrid (DH) DFT. A benchmark comparison is performed to assess the quality of DH-DFT data using reference energy differences obtained using a complete-basis-limit coupled-cluster (CC) composite method. The associated large-basis CC computations were enabled by a new general-purpose pseudopotential capability implemented as part of this work. Dispersion-corrected parameter-free DH-DFT methods, namely PBE0-DH-D3(BJ) and PBE-QIDH-D3(BJ), provided mean unsigned errors within chemical accuracy (1 kcal mol -1 ) for a set of barrier heights corresponding to the most energetically favorable initiation steps. The hydrolysis mechanism is found to proceed via intermolecular hydrogen transfer within van der Waals complexes involving UF 6 , UF 5 OH, and UOF 4 , in agreement with previous studies, followed by the formation of a previously unappreciated dihydroxide intermediate, UF 4 (OH) 2 . The dihydroxide is predicted to form under both kinetic and thermodynamic control, and, unlike the alternate pathway leading to the UO 2 F 2 monomer, its reaction energy is exothermic, in agreement with observation. Finally, harmonic and anharmonic vibrational simulations are performed to reinterpret literature infrared spectroscopy in light of this newly identified species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Photocatalytic Hydrogen Evolution Using Mesoporous Honeycomb Iron Titanate

Mesoporous honeycomb iron titanate using a sol-gel, evaporation-induced self-assembly method is synthesized. A triblock copolymer, F127, serves as a structure-directing agents, with iron chloride and titanium (IV) isopropoxide as inorganic precursors. The strong intermolecular force of attraction among urea, metal precursors, and polymer led to the formation of the mesoporous honeycomb structure. The study of physicochemical properties using different techniques reveals the formation of microstructures with a remarkable degree of porosity. The amorphous iron titanate outperforms the photochemical generation of H 2 due to its disorderly structural arrangement and incomplete crystal formation. The randomness on the structure provides more area for catalytic reaction by providing more contact with the reactant and superior light absorption capability. The high amount of hydrogen gas, 40.66 mmolg -1 h -1 , is observed in the investigation over 3 h of activity for the iron titanate honeycomb sample. This yield is a more significant amount compared to the obtained for the commercially available TiO 2 (23.78 mmolg -1 h -1 ). The iron titanate materials synthesized with low-cost materials and methods are very effective and have the potential for hydrogen generation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2020

Contribution of the Hydroxyl Group on Interfacial Heat Conduction of Monohydric Alcohols: A Molecular Dynamics Study

Monohydric alcohols have been used as promising phase change materials (PCMs) for low-temperature latent heat storage. However, the heat storage/retrieval rates are limited due to the low thermal conductivity of such alcohols. In this work, nonequilibrium molecular dynamics (NEMD) simulations were performed to study the microscopic heat conduction in example monohydric alcohols, i.e., 1-dodecanol (C12H26O), 1-tetradecanol (C14H30O), and 1-hexadecanol (C16H34O). A simplified ideal crystal model was proposed to exploit the potential for improving the thermal conductivity of monohydric alcohols. The effect of ideal crystalline structures, especially the contribution of the hydroxyl group, on the microscopic heat conduction process was analyzed. The thermal conductivity of the ideal crystals of the various monohydric alcohols was predicted to be more than twice as compared to that of their respective solids. The major thermal resistance in the ideal crystals was found around the molecular interfaces, as a result of the excellent heat conduction performance along the linear molecular chains. The calculated vibrational density of states (VDOS) and interfacial heat transfer were then investigated. When the interfaces are surrounded by hydroxyl groups as walls, strong hydrogen bond (HB) interactions were observed. The interfacial heat transfer coefficient of the ideal crystalline structures of 1-tetradecanol was found to reach up to ~735.6 MW/m2 W. It was elucidated that the high interfacial heat transfer rate is clearly related to the stronger intermolecular interactions.

Engineering↗

DOE OSTI2022

Optical readout of singlet fission biexcitons in a heteroacene with photoluminescence detected magnetic resonance

Molecular spin systems based on photoexcited triplet pairs formed via singlet fission (SF) are attractive as carriers of quantum information because of their potentially pure and controllable spin polarization, but developing systems that offer optical routes to readout as well as initialization is challenging. Herein, we characterize the electron spin magnetic resonance change in the photoluminescence intensity for a tailored organic molecular crystal while sweeping a microwave drive up to 10 GHz in a broadband loop structure. We observe resonant transitions for both triplet and quintet spin sublevel populations showing their optical sensitivity and revealing the zero-field parameters for each. We map the evolution of these spectra in both microwave frequency and magnetic field, producing a pattern of optically detected magnetic resonance (ODMR) peaks. Fits to these data using a suitable model suggest significant spin polarization in this system with orientation selectivity. Unusual excitation intensity dependence is also observed, which inverts the sign of the ODMR signal for the triplet features, but not for the quintet. These observations demonstrate optical detection of the spin sublevel population dictated by SF and intermolecular geometry, and highlight anisotropic and multi-scale dynamics of triplet pairs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Gold-Standard Chemical Database 137 (GSCDB137): A Diverse Set of Accurate Energy Differences for Assessing and Developing Density Functionals

We present GSCDB137, a rigorously curated benchmark library of 137 data sets (8377 entries) covering main-group and transition-metal reaction energies and barrier heights, (intra- and intermolecular) noncovalent interactions, dipole moments, polarizabilities, electric-field response energies, and vibrational frequencies. Legacy data from GMTKN55 and MGCDB84 have been updated to today's best reference values; redundant or low-quality points were removed, and many new, property-focused sets were added. Testing 29 popular density functional approximations (DFAs) confirms the expected Jacob's-ladder hierarchy overall but also reveals notable exceptions: functional performance for frequencies and electric-field properties correlates poorly with that for other ground-state energetics. ωB97M-V and ωB97X-V are the most balanced hybrid meta-GGA and hybrid GGA, respectively; B97M-V and revPBE-D4 lead the meta-GGA and GGA classes. Double hybrids lower mean errors by about 30% versus their hybrid analogues but demand careful frozen-core, basis set, and spin contamination treatment. GSCDB137 offers a comprehensive, openly documented platform for rigorous validation of DFA and universal machine learning potentials, and training of the next generation of exchange-correlation functionals.

Liang, Jiashu [University of California, Berkeley,↗

DOE OSTI2024

Structural basis for inhibition of coagulation factor VIII reveals a shared antigenic hotspot on the C1 domain

Hemophilia A arises from dysfunctional or deficient coagulation factor (F)VIII and leads to inefficient fibrin clot formation and uncontrolled bleeding events. The development of antibody inhibitors is a clinical complication in hemophilia A patients receiving FVIII replacement therapy. LE2E9 is an anti-C1 domain inhibitor previously isolated from a mild/moderate hemophilia A patient and disrupts FVIII interactions with von Willebrand factor and FIXa, though the intermolecular contacts that underpin LE2E9-mediated FVIII neutralization are undefined. To determine the structure of the complex between FVIII and LE2E9 and characterize its mechanism of inhibition. FVIII was bound to the antigen binding fragment (Fab) of NB2E9, a recombinant construct of LE2E9, and its structure was determined by cryogenic electron microscopy. Here, this report communicates the 3.46 Å structure of FVIII bound to NB2E9, with its epitope comprising FVIII residues S2040 to Y2043, K2065 to W2070, and R2150 to H2155. Structural analysis reveals that the LE2E9 epitope overlaps with portions of the epitope for 2A9, a murine-derived inhibitor, suggesting that these residues represent a shared antigenic region on the C1 domain between FVIII –/– mice and hemophilia A patients. Furthermore, the FVIII:NB2E9 structure elucidates the orientation of the LE2E9 glycan, illustrating how the glycan sterically blocks interactions between the FVIII C1 domain and the von Willebrand factor D' domain. A putative model of the FVIIIa:FIXa complex suggests potential clashing between the NB2E9 glycan and FIXa light chain. These results describe an antigenic “hotspot” on the FVIII C1 domain and provide a structural basis for engineering FVIII replacement therapeutics with reduced antigenicity.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Rotational and translational diffusion of liquid n-hexane: EFP-based molecular dynamics analysis

Molecular Dynamics (MD) simulations based on the Effective Fragment Potential (EFP) method are utilized to provide a comprehensive assessment of diffusion in liquid n-hexane. We decompose translational diffusion into components along and orthogonal to the long axis of the molecule. Rotational diffusion is decomposed into tumbling and spinning motions about this axis. Our analysis yields four corresponding diffusion coefficients which are related to diagonal entries in the complete 6 × 6 diffusion tensor accounting for the three rotational and three translational degrees of freedom and for the potential coupling between them. However, coupling between different degrees of freedom is expected to be minimal for a natural choice of the molecular body-fixed axis, so then off-diagonal entries in the tensor are negligible. This expectation is supported by a hydrodynamic analysis of the diffusion tensor which treats the liquid surrounding the molecule being tracked as a viscous continuum. Thus, the EFP MD analysis provides a comprehensive characterization of diffusion and also reveals expected shortcomings of the hydrodynamic treatment, particularly for rotational diffusion, when applied to neat liquids.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Hydrogen Bonding in Amorphous Indomethacin

Amorphous Indomethacin has enhanced bioavailability over its crystalline forms, yet amorphous forms can still possess a wide variety of structures. Here, Empirical Potential Structure Refinement (EPSR) has been used to provide accurate molecular models on the structure of five different amorphous Indomethacin samples, that are consistent with their high-energy X-ray diffraction patterns. It is found that the majority of molecules in amorphous Indomethacin are non-bonded or bonded to one neighboring molecule via a single hydrogen bond, in contrast to the doubly bonded dimers found in the crystalline state. The EPSR models further indicate a substantial variation in hydrogen bonding between different amorphous forms, leading to a diversity of chain structures not found in any known crystal structures. The majority of hydrogen bonds are associated with the carboxylic acid group, although a significant number of amide hydrogen bonding interactions are also found in the models. Evidence of some dipole–dipole interactions are also observed in the more structurally ordered models. The results are consistent with a distribution of Z-isomer intramolecular type conformations in the more disordered structures, that distort when stronger intermolecular hydrogen bonding occurs. The findings are supported by 1H and 2H NMR studies of the hydrogen bond dynamics in amorphous Indomethacin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree–Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. Here, the accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, R cut . R cut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small R cut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Long-Range Fit: A Software Package for the Representation and Study of Long-Range Molecular Interactions

Describing intermolecular forces is fundamental to modeling and predicting the behavior of molecular systems. In particular, long-range molecular interactions─with electrostatic, induction, and dispersion as the main components─play a critical role, especially for low-temperature and low-density regimes. Long-range interactions are often described through perturbation theory, representing the electronic charge distribution via a multipolar series of the moments and polarizability tensors corresponding to each molecule. However, while the theory is well established, obtaining the resulting analytical expressions (and their practical implementation) constitutes a highly complex and system-dependent task. To address this challenge, we developed long-range-fit (LRF), an interactive and user-friendly software package designed to automate the generation and fitting of long-range interaction terms for arbitrary molecules in nondegenerate (ground or excited) electronic states. We have derived and implemented all terms up to 15th order, without approximations, via a spherical tensor representation, with symmetry adaptation to all molecular point-group symmetries. The resulting potential energy surface is compatible with most representations of the close interaction region.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗