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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
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 OSTI2021

Optimizing the self-assembly of conjugated polymers and small molecules through structurally programmed non-covalent control

Organic conjugated polymers and oligomers are key electronic materials for applications such as transistors, photovoltaics, and light emitting devices due to their potential for solution processability, mechanical flexibility, and precise structure-based tuning compared to inorganic materials. In dilute environments, the optoelectronic properties of conjugated polymers are largely governed by their constitutional structure and, to a lesser degree, their solution-state intramolecular configuration. In the solid state, intramolecular conformation and intermolecular electronic coupling impact these properties substantially, especially in relation to device performance. Therefore, an increasingly important area of research concerning conjugated materials is developing design strategies aimed at optimizing the solid-state packing for electronic applications. Programming solid-state packing arrangements through discrete non-covalent interactions is an emerging strategy within the context of conjugated polymers. This review focuses on the use of the two most prevalent discrete and directional interactions used to dictate the self-assembly of conjugated polymers and oligomers—hydrogen bonds and chalcogen bonds. Further, we also discuss how these design motifs can imbue conjugated materials with appealing physical properties while simultaneously retaining or improving electronic capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Pressure–Modulated Luminescence Enhancement and Quenching in a Hydrogen–Bonded Organic Framework

Light emission in the solid state is central for illumination, sensing, and imaging applications. Unlike luminescence in dilute solutions, where the excited states are unimolecular in nature, intermolecular interaction plays a significant role in the quantum yield of solid-state luminophores, manifested as competing aggregation-caused quenching (ACQ) and aggregation-induced enhancement (AIE). Both effects are extensively studied in various systems; however, it remains unclear how their competition depends on molecular conformation and intermolecular stacking. Here the direct observation of pressure-modulated AIE-ACQ competition in a crystalline hydrogen-bonded organic framework (HOF) is reported. Using in situ spectroscopies and computational modeling, the intramolecular vibration and intermolecular π–π stacking directly responsible for the non-radiative decay of the excited state are identified. The extent of these two contributions is modulated by hydrostatic pressure and guest molecules in the HOF pores. Furthermore, this work demonstrates a physically neat model system to understand and control solid-state luminescence, and a potential material platform for piezoluminescent sensing.

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

Collisional excitation of interstellar PN by H 2 : New interaction potential and scattering calculations

Here, rotational excitation of interstellar PN molecules induced by collisions with H 2 is investigated. We present the first ab initio four-dimensional potential energy surface (PES) for the PN–H 2 van der Waals system. The PES was obtained using an explicitly correlated coupled cluster approach with single, double, and perturbative triple excitations [CCSD(T)-F12b]. The method of interpolating moving least squares was used to construct an analytical PES from these data. The equilibrium structure of the complex was found to be linear, with H 2 aligned at the N end of the PN molecule, at an intermolecular separation of 4.2 Å. The corresponding well-depth is 224.3 cm -1 . The dissociation energies were found to be 40.19 cm -1 and 75.05 cm -1 for complexes of PN with ortho-H 2 and para-H 2 , respectively. Integral cross sections for rotational excitation in PN–H 2 collisions were calculated using the new PES and were found to be strongly dependent on the rotational level of the H 2 molecule. These new collisional data will be crucial to improve the estimation of PN abundance in the interstellar medium from observational spectra.

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

Investigating intermolecular interactions among CO 2 , water and PEEK-ionene membrane using cryo ToF-SIMS and isotopic labeling

Cryogenic time-of-flight secondary ion mass spectrometry (cryo ToF-SIMS) has emerged as a powerful tool for investigating molecular interactions, speciation, and dynamics in materials for CO 2 capture. In this study, we apply cryo ToF-SIMS to probe interactions between CO 2 , water, and PEEK-ionene membranes—a promising material for direct CO 2 capture due to its selectivity, durability, and efficiency. Despite this potential, the mechanisms governing CO 2 diffusion and the influence of water vapor on CO 2 behavior remain unclear. To address this, we loaded PEEK-ionene membranes with 13 CO 2 and D 2 O and employed cryo ToF-SIMS to visualize the 3D distribution of CO 2 and water within the membrane. While prior studies suggest that 13 CO 2 is absorbed under ambient conditions, our cryo ToF-SIMS analysis revealed no enhancement of the 13 C/ 12 C ratio, suggesting weak CO 2 -membrane interactions. As a result, CO 2 vaporizes even at low temperatures (−140°C) under vacuum conditions. In contrast, D 2 O displayed a relatively homogeneous distribution in the membrane, suggesting stronger water-membrane interactions via hydrogen bonding (18–20 kJ/mol). Interestingly, CO 2 was not detected in D 2 O-loaded membranes, indicating minimal interference from water vapor on CO 2 diffusion. As a comparison, the cryo ToF-SIMS data show that 13 CO 2 can readily react with a basic Na 2 CO 3 aqueous solution to form NaH 13 CO 3 . These findings demonstrate cryo ToF-SIMS as a critical technique for understanding gas-water-membrane interactions, offering insights for membrane functionalization to improve CO 2 capture efficiency.

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

Structure, Hydrodynamics, and Phase Transition of Freely Suspended Liquid Crystals

Smectic liquid crystals are phases of rod shaped molecules organized into one dimensionally (1D) periodic arrays of layers, each layer being between one and two molecular lengths thick. In the least ordered smectic phases, the smectics A and C, each layer is a two dimensional (2D) liquid. Additionally there are a variety of more ordered smectic phases having hexatic short range translational order or 2D crystalline quasi long range translational order within the layers. The inherent fluid-layer structure and low vapor pressure of smectic liquid crystals enable the long term stabilization of freely suspended, single component, layered fluid films as thin as 30A, a single molecular layer. The layering forces the films to be an integral number of smectic layers thick, quantizing their thickness in layer units and forcing a film of a particular number of layers to be physically homogeneous with respect to its layer structure over its entire area. Optical reflectivity enables the precise determination of the number of layers. These ultrathin freely suspended liquid crystal films are structures of fundamental interest in condensed matter and fluid physics. They are the thinnest known stable condensed phase fluid structures and have the largest surface-to-volume ratio of any stable fluid preparation, making them ideal for the study of the effects of reduced dimensionality on phase behavior and on fluctuation and interface phenomena. Their low vapor pressure and quantized thickness enable the effective use of microgravity to extend the study of basic capillary phenomena to ultrathin fluid films. Freely suspended films have been a wellspring of new liquid crystal physics. They have been used to provide unique experimental conditions for the study of condensed phase transitions in two dimensions. They are the only system in which the hexatic has been unambiguously identified as a phase of matter, and the only physical system in which fluctuations of a 2D XY system and Kosterlitz Thouless phase transition has been observed and 2D XY quasi long range order verified. Smectic films have enabled the precise determination of smectic layer electron density and positional fluctuation profile and have been used to show that the interlayer interactions in anti-ferroelectric tilted smectics do not extend significantly beyond nearest neighbors. The interactions which are operative in liquid crystals are generally weak in comparison to those in crystalline phases, leading to the facile manipulation of the order in liquid crystals by external agents such as applied fields and surfaces. Effects arising from weak ordering are significantly enhanced in ultrathin free films and filaments wherein the intermolecular coupling is effectively reduced by loss of neighbors. Over the past four years this research, which we now detail, has produced a host of exciting new discoveries and unexpected results, maintaining the position of the study of freely suspended liquid crystal structures as one of most exciting and fruitful areas of complex fluid physics. In addition, several potentially interesting microgravity free film experiments have been identified.

Clark, Noel A.↗
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

Local structure analysis of low-temperature neutron pair distribution function coupled with molecular dynamics simulations of CH 4 and CO 2 hydrates from 2 to 210 K

CH 4 hydrates occur naturally and are an abundant potential fuel source with a corresponding risk of potent greenhouse gas release, due to their stability conditions at low temperature and high pressure. Byproduct CO 2 can be exchanged with CH 4 in these natural deposits and can potentially stabilize the hydrates at higher temperatures. CH 4 , CO 2 , and mixed CH 4 -CO 2 hydrates are studied with in situ neutron pair distribution function experiments from 2 to 210 K to investigate the impact of varying the CH 4 -CO 2 guest composition in the gas hydrate structure. These experiments combined with Reverse Monte Carlo analysis allow for the characterization of intermolecular CH 4 and CO 2 interactions with the water molecules which form the hydrate lattice and how they impact the local structure of the lattice itself. Results indicate that when CH 4 and CO 2 co-occupy the hydrate, the host is more strongly distorted than in the pure CH 4 and CO 2 hydrates, but this becomes less defined with increasing temperature. The presence of CO 2 in mixed hydrate increases the stability range and creates a barrier for CH 4 to completely leave the structure.

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

Highly Emissive Semi-Ladder-Type Copolymers, Aggregation State, and Solution-Processed Organic Light-Emitting Transistor

An organic light-emitting transistor (OLET) integrates the logic function of a field-effect transistor (FET) with electroluminescence of a light-emitting diode (LED) into a single device. This has the potential to play an important role in smart displays, electrically pumped organic laser, and other flexible electronic devices. For achieving high-performance solution-processed OLET, three semiladder copolymers—TPTI-CC, TPTI-C, and TPTI-F—were developed in this study. It was found that variation in the chemical structure lead to significant differences in energy level, emission color, and intermolecular aggregation. High photoluminescence quantum yields (PLQY) was obtained in these three polymers. Solution-processed trilayered OLET devices were fabricated, which exhibited strong electroluminescence, balanced charge mobility, and external quantum efficiency (EQE) of 2.8%.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Structural Characterisation and DFT-Aided Interpretation of Vibrational Spectra for Cyclo(l-Cys-d-Cys) Cyclic Dipeptide in a Solid State

Cyclic dipeptides with two intramolecular peptide bonds forming a six-membered 2,5-diketopiperazine ring are gaining significant attention due to their biological and chemical properties. Small changes in the local geometry of such molecules (from cis to trans) can lead to significant structural differences. This work presents the results of a study of cyclo(l-Cys-d-Cys), a dipeptide comprising two cysteine molecules in opposite chiral configurations, with the functional groups situated at both sides of the diketopiperazine ring. X-ray diffraction (XRD) experiment revealed that the molecule crystallises in the $P$-1 space group, which includes the centre of inversion. The IR and Raman vibrational spectra of the molecule were acquired and interpreted in terms of the potential energy distribution (PED) according to the results of density functional theory (DFT) calculations. The DFT-assisted analysis of energy frameworks for the hydrogen bond network within molecular crystals was performed to support the interpretation of X-ray structural data. The optimisation of the computational model based on three-molecule geometry sections from the crystallographic structure, selected to appropriately reflect the intermolecular interactions responsible for the formation of 1D molecular tapes in cyclo(l-Cys-d-Cys) crystal, allowed for better correspondence between theoretical and experimental vibrational spectra. This work can be considered the first complete structural characterisation of cyclo(l-Cys-d-Cys), complemented via vibrational spectroscopy results with full band assignment aided with the use of the DFT method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1994

Dielectric Properties of Poly(ethylene oxide) from Molecular Dynamics Simulations

The order, conformations and dynamics of poly(oxyethylene) (POE) melts have been investigated through molecular dynamics simulations. The potential energy functions were determined from detailed ab initio electronic structure calculations of the conformational energies of the model molecules 1,2-dimethoxyethane (DME) and diethylether. The x-ray structure factor for POE from simulation will be compared to experiment. In terms of conformation, simulations reveal that chains are extended in the melt relative to isolated chains due to the presence of strong intermolecular O...H interactions, which occur at the expense of intramolecular O...H interactions. Conformational dynamics about the C-C bond were found to be significantly faster than in polymethylene, while conformational dynamics about the C-O bond even faster than the C-C dynamics. The faster local dynamics in POE relative to polymethylene is consistent with C-13 NMR spin-lattice relaxation experiments. Conformational transitions showed significant second-neighbor correlation, as was found for polymethylene. This correlation of transitions with C-C neighbors was found to be reduced relative to C-O neighbors. Dielectric relaxation from simulation will also be compared with experiment.

Smith, Grant D.↗
DOE OSTI2023

Wetting of a Hydrophobic Surface: Far-IR Action Spectroscopy and Dynamics of Microhydrated Naphthalene

The interaction of water and polycyclic aromatic hydrocarbons is of fundamental importance in areas as diverse as materials science and atmospheric and interstellar chemistry. The interplay between hydrogen bonding and dipole–π interactions results in subtle dynamics that are challenging to describe from first principles. Here, we employ far-IR action vibrational spectroscopy with the infrared free-electron laser FELIX to investigate naphthalene with one to three water molecules. We observe diffuse bands associated with intermolecular vibrational modes that serve as direct probes of the loose binding of water to the naphthalene surface. These signatures are poorly reproduced by static DFT or Møller–Plesset computations. Instead, a rationalization is achieved through Born–Oppenheimer Molecular Dynamics simulations, revealing the active mobility of water over the surface, even at low temperatures. Therefore, our work provides direct insights into the wetting interactions associated with shallow potential energy surfaces while simultaneously demonstrating a solid experimental–computational framework for their investigation.

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

Raman spectrum and polarizability of liquid water from deep neural networks

We introduce a scheme based on machine learning and deep neural networks to model the environmental dependence of the electronic polarizability in insulating materials. Application to liquid water shows that training the network with a relatively small number of molecular configurations is sufficient to predict the polarizability of arbitrary liquid configurations in close agreement with ab initio density functional theory calculations. In combination with a neural network representation of the interatomic potential energy surface, the scheme allows us to calculate the Raman spectra along 2-nanosecond classical trajectories at different temperatures for H 2 O and D 2 O. The vast gains in efficiency provided by the machine learning approach enable longer trajectories and larger system sizes relative to ab initio methods, reducing the statistical error and improving the resolution of the low-frequency Raman spectra. Decomposing the spectra into intramolecular and intermolecular contributions elucidates the mechanisms behind the temperature dependence of the low-frequency and stretch modes.

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

Observation of a super-tetrahedral cluster of acetonitrile-solvated dodecaborate dianion via dihydrogen bonding

We launched a combined negative ion photoelectron spectroscopy and multiscale theoretical investigation on the geometric and electronic structures of a series of acetonitrile-solvated dodecaborate clusters, i.e., B 12 H 12 2- ·nCH 3 CN (n = 1–4). The electron binding energies of B 12 H 12 2- ·nCH 3 CN are observed to increase with cluster size, suggesting their enhanced electronic stability. B3LYP-D3(BJ)/ma-def2-TZVP geometry optimizations indicate each acetonitrile molecule binds to B 12 H 12 2- via a threefold dihydrogen bond (DHB) B3–H3 ⋮⋮⋮ H3C–CN unit, in which three adjacent nucleophilic H atoms in B 12 H 12 2- interact with the three methyl hydrogens of acetonitrile. The structural evolution from n = 1 to 4 can be rationalized by the surface charge redistributions through the restrained electrostatic potential analysis. Notably, a super-tetrahedral cluster of B 12 H 12 2- solvated by four acetonitrile molecules with 12 DHBs is observed. The post-Hartree–Fock domain-based local pair natural orbital- coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)] calculated vertical detachment energies agree well with the experimental measurements, confirming the identified isomers as the most stable ones. Furthermore, the nature and strength of the intermolecular interactions between B 12 H 12 2- and CH 3 CN are revealed by the quantum theory of atoms-in-molecules and the energy decomposition analysis. Ab initio molecular dynamics simulations are conducted at various temperatures to reveal the great kinetic and thermodynamic stabilities of the selected B 12 H 12 2- ·CH 3 CN cluster. The binding motif in B 12 H 12 2- ·CH 3 CN is largely retained for the whole halogenated series B 12 X 12 2- ·CH 3 CN (X = F–I). This study provides a molecular-level understanding of structural evolution for acetonitrile-solvated dodecaborate clusters and a fresh view by examining acetonitrile as a real hydrogen bond (HB) donor to form strong HB interactions.

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

Partition function corrections

The vibration-rotation interaction coupling in diatomic molecules is analyzed and first-order expressions for the energy levels are deduced. Corrections to the partition function to account for this coupling are derived. The effects of long-range intermolecular forces due to induced dipole-induced dipole interactions, dipole-induced dipole interactions, charge-dipole interactions, and charge-charge interactions are considered. London's quantized oscillator model of the induced dipole-induced dipole interaction is developed. The partition function including effects of such dispersion forces is derived and is shown to lead to the virial equation of state. For a model combining the hard sphere and long-range attractive type potentials, this is found to lead to the van der Waals equation of state. The equation of state near the critical point is presented in reduced form (thermodynamic variables expressed in units of critical point values). Finally, approximate corrections to the electronic partition functions are derived to account for perturbations of excited electronic states induced by neighboring neutral particles, heavy ions, and fast electron collisions.

Source record↗
NASA NTRS2004

Structure and Dynamics of Freely Suspended Liquid Crystals

Smectic liquid crystals are phases of rod shaped molecules organized into one dimensionally (1 D) periodic arrays of layers, each layer being between one and two molecular lengths thick. In the least ordered smectic phases, the smectics A and C, each layer is a two dimensional (2D) liquid. Additionally there are a variety of more ordered smectic phases having hexatic short range translational order or 2D crystalline or quasi long range translational order within the layers. The inherent fluid-layer structure and low vapor pressure of smectic liquid crystals enables the long term stabilization of freely suspended, single component, layered fluid films as thin as 30A, a single molecular layer. The layering forces the films to be an integral number of smectic layers thick, quantizing their thickness in layer units and forcing a film of a particular number of layers to be physically homogeneous with respect to its layer structure over its entire area. Optical reflectivity enables the precise determination of the number of layers. These ultrathin freely suspended liquid crystal films are structures of fundamental interest in condensed matter and fluid physics. They are the thinnest known stable fluid structures and have the largest surface-to-volume ratio of any stable fluid preparation, making them ideal for the study of the effects of reduced dimensionality on phase behavior and on fluctuation and interface phenomena. Their low vapor pressure and quantized thickness enable the effective use of microgravity to extend the study of basic capillary phenomena to ultrathin fluid films. Freely suspended films have been a wellspring of new LC physics. They have been used to provide unique experimental conditions for the study of condensed phase transitions in two dimensions. They are the only system in which the hexatic has been unambiguously identified as a phase of matter, and the only physical system in which fluctuations of a 2D XY system and Kosterlitz Thouless phase transition has been observed and 2D XY quasi long range order verified. Smectic films have enabled the precise determination of smectic layer electron density and positional fluctuation profiles and have been used to show that the interlayer interactions in antiferroelectric tilted smectics do not extend significantly beyond nearest neighbors. Freely suspended films played a pivotal role in the recent discovery of macroscopic chiral-polar ordering in fluids of achiral molecules. The interactions which are operative in liquid crystals are generally weak in comparison to those in crystalline phases, leading to the facile manipulation of the order in liquid crystals by external agents such as applied fields and surfaces. Effects arising from weak ordering are significantly enhanced in ultrathin free films and filaments, in which the intermolecular coupling is effectively further reduced by loss of neighbors. Over the past four years this research, which we now detail, has produced a host of exciting new discoveries and unexpected results, maintaining the study of freely suspended liquid crystal structures as one of most exciting and fruitful areas of complex fluid physics. In addition, a class of experiments on the behavior of 1D interfaces in 2D films have been pursued with results that point to potentially quite interesting effects in microgravity.

Clark, Noel A.↗
DOE OSTI2026

Predicting compatibility between ferredoxins and the Fe protein of nitrogenase using in silico protein modeling

Biological nitrogen fixation is the process by which certain bacteria and archaea use the enzyme nitrogenase to reduce atmospheric nitrogen into bioavailable ammonium. Engineering non‐nitrogen‐fixing organisms, like plants, to use nitrogenase could reduce dependency on synthetic fertilizer and mitigate the environmental impacts of industrial fertilizer production. However, nitrogenase activity requires delivery of reducing power by small electron carrying proteins known as ferredoxins and flavodoxins, and successfully engineering nitrogenase into new systems will require a mechanistic understanding of electron delivery by these proteins. Most organisms often have multiple ferredoxins, raising the question of which ferredoxin can support nitrogenase activity. The purpose of this study is to gain insight into how we can predict which ferredoxin is compatible with the Fe protein, the component of nitrogenase that interacts with ferredoxin or flavodoxin. Our in silico protein–protein docking simulations reveal that most ferredoxins and flavodoxins involved in nitrogen fixation have the shortest distance (≤10 Å) between their redox cofactor and the [4Fe‐4S] cluster of the Fe protein. We found shorter cofactor distance contributes to faster intermolecular electron tunneling rates. Bacterial ferredoxins that play a role in nitrogen fixation also exhibit more complementary interactions with the Fe protein than bacterial and plant ferredoxins not involved in this process. Heterologous expression of a set of ferredoxins from both nitrogen‐fixing and non‐nitrogen‐fixing bacteria in the diazotroph Rhodopseudomonas palustris supports our model‐derived prediction that shorter distances between the electron‐carrying cofactors favor nitrogenase compatibility. These findings offer a framework to predict and potentially enhance ferredoxin–nitrogenase compatibility, which will help to improve our ability to engineer nitrogen fixation into non‐nitrogen‐fixing organisms like plants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Effects of Hydrogen Bonding on Nuclear Data Development of Liquid Anhydrous HF

Anhydrous Hydrogen Fluoride (HF) at high temperatures and pressures is used to process and manufacture nuclear fuel. As HF is often used directly with uranium, correct neutron thermal scattering cross sections are crucial to criticality safety applications. Classical molecular dynamics (CMD) simulation of the flexible HF system was used to create the thermal scattering law (TSL) and cross sections. The initial 2-site model is used in LAMMPS, and it can not capture the H-bond. To correctly represent the H-bond effects, a second, 3-site model was constructed in GROMACS. The 3-site model handled H-bonds by connecting a massless charge to the molecule. Key model parameters were compared to experimental data to verify the approach and models. To get the normalized VACF, the model was compared using hydrogen and fluorine bond length, density, potential energy, and diffusion coefficient. The phonon DOSs for both models were derived from the normalized VACF. DOSs were used to estimate the TSL ( S ( α, β )) and neutron thermal scattering cross sections for hydrogen in HF. The TSLs were evaluated using the FLASSH code with the Schofield diffusion model. It was observed that the representation of the hydrogen bonding changes the TSL's diffusional contributions. This is represented in the low energy scattering cross section, where intermolecular binding effects shift the cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗