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.
Learning intermolecular forces at liquid–vapor interfaces
By adopting a perspective informed by contemporary liquid-state theory, we consider how to train an artificial neural network potential to describe inhomogeneous, disordered systems. Here, we find that neural network potentials based on local representations of atomic environments are capable of describing some properties of liquid-vapor interfaces but typically fail for properties that depend on unbalanced long-ranged interactions that build up in the presence of broken translation symmetry. These same interactions cancel in the translationally invariant bulk, allowing local neural network potentials to describe bulk properties correctly. By incorporating explicit models of the slowly varying long-ranged interactions and training neural networks only on the short-ranged components, we can arrive at potentials that robustly recover interfacial properties. We find that local neural network models can sometimes approximate a local molecular field potential to correct for the truncated interactions, but this behavior is variable and hard to learn. Generally, we find that models with explicit electrostatics are easier to train and have higher accuracy. We demonstrate this perspective in a simple model of an asymmetric dipolar fluid, where the exact long-ranged interaction is known, and in an ab initio water model, where it is approximated.
Intermolecular vibrational energy transfer enabled by microcavity strong light–matter coupling
Selective vibrational energy transfer between molecules in the liquid phase, a difficult process hampered by weak intermolecular forces, is achieved through polaritons formed by strong coupling between cavity photon modes and donor and acceptor molecules. Using pump-probe and two-dimensional infrared spectroscopy, we found that the excitation of the upper polariton, which is composed mostly of donors, can efficiently relax to the acceptors within ~5 picoseconds. The energy-transfer efficiency can be further enhanced by increasing the cavity lifetime, suggesting that the energy transfer is a polaritonic process. This vibrational energy-transfer pathway opens doors for applications in remote chemistry, sensing mechanisms, and vibrational polariton condensation.
Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces
Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems
Scale translation yields insights into gas adsorption under nanoconfinement
This work describes a scale-translating simulation framework to investigate gas adsorption behavior in nanoconfined pores. The framework combines molecular simulations (MSs), equation of state (EoS), and lattice Boltzmann (LB) simulations. MSs reveal the physics of methane adsorption in nano-sized pores, where input values of fugacity coefficients are optimized based on EoS predictions. Then, an LB free-energy model, which incorporates a viral EoS, upscales intermolecular forces and estimates adsorption behavior via a proposed fluid–wall interaction model. Armed with the values of the LB interaction parameter as a function of pressure, the LB model is used to predict fluid behavior in irregular nanopores, and the results are validated against reference MS data. The LB model is then used to study adsorption behavior at a continuum scale in representative organic shale nanopores based on finely characterized Vaca Muerta shale samples. Furthermore, the results show that methane adsorption could significantly increase contained fluids by 10%–25% in pores smaller than 20 nm. However, in larger pores (40 nm to 90 nm), adsorption's impact diminishes to 2%–3%, suggesting sorption's negligible role beyond a 40 nm pore size.
Long timescale solvation dynamics and confinement: The case of non-ionic deep eutectic solvents of lauric acid and N-methylacetamide
Microscopic segregation and molecular heterogeneities in complex liquids are the result of the interplay between different intermolecular forces, all of which contribute to the energy landscape of the system. A consequence of the intricate energy landscape is the nontrivial effect on the solvation dynamics. Here, the impact of molecular heterogeneities on the solvation dynamics is studied using infrared spectroscopies and molecular dynamics simulations. In particular, this study focuses on the dynamical effect of nanoscopic heterogeneities present in deep eutectic solvents (DESs) composed of lauric acid (LA) and N-methylacetamide (NMA). To this end, a molecular probe containing a carbon triple bond is used as an infrared reporter. The results show that the vibrational probe is likely to be located in the NMA polar domains. Furthermore, the probe solvation dynamics derived from the 2DIR spectra presents a slowdown of its timescale with increasing LA concentration in the DES. Kubo modeling of the probe solvation dynamics shows a correlation between the amplitude of its long time component and the presence of molecular heterogeneities in the sample. Semiclassical modeling of the vibrational line shape of the triple bond stretch demonstrates that the heterogeneities affect the whole solvation dynamics of the system through the amplitudes of the frequency fluctuations. Molecular dynamics simulations confirm the experimental results and their interpretation by showing a slowdown of the solvation dynamics when the LA heterogeneities are present. Overall, the study presents a molecular framework to explain the effect of confinement created by nanoscopic LA heterogeneities on the solvation dynamics of the system.
First-principles study of an S = 1 quasi one-dimensional quantum molecular magnetic material
We use density functional theory to study the structural, magnetic, and electronic structures of the organometallic quantum magnet NiCl 2-4 SC (NH 2 ) 2 (DTN). Recent work has demonstrated the quasi one-dimensional nature of the molecular crystal and studied its quantum phase transitions at low temperatures. The system includes a magnetoelectric (ME) coupling and, when doped with Br, the presence of an exotic Bose-glass state. Using the generalized gradient approximation with inclusion of a van der Waals term to account for weak intermolecular forces and by introducing a Hubbard U term to the total energy, we systematically show that our calculations reproduce the magnetic anisotropy, the intermolecular exchange coupling strength, and the magnetoelectric effect in DTN, which have been observed in previous experiments. Further analysis of the electronic structure gives insight into the underlying magnetic interactions, including what mechanisms may be causing the ME effect. Using this computationally efficient model, we predict what effect applying an electric field might have on the magnetic properties of this quantum magnet.
Reinterpreting π-stacking
The nature of π–π interactions has long been debated. The term “π-stacking” is considered by some to be a misnomer, in part because overlapping π-electron densities are thought to incur steric repulsion, and the physical origins of the widely-encountered “slip-stacked” motif have variously been attributed to either sterics or electrostatics, in competition with dispersion. Here, we use quantum-mechanical energy decomposition analysis to investigate π–π interactions in supramolecular complexes of polycyclic aromatic hydrocarbons, ranging in size up to realistic models of graphene, and for comparison we perform the same analysis on stacked complexes of polycyclic saturated hydrocarbons, which are cyclohexane-based analogues of graphane. Our results help to explain the short-range structure of liquid hydrocarbons that is inferred from neutron scattering, trends in melting-point data, the interlayer separation of graphene sheets, and finally band gaps and observation of molecular plasmons in graphene nanoribbons. Analysis of intermolecular forces demonstrates that aromatic π–π interactions constitute a unique and fundamentally quantum-mechanical form of non-bonded interaction. Not only do stacked π–π architectures enhance dispersion, but quadrupolar electrostatic interactions that may be repulsive at long range are rendered attractive at the intermolecular distances that characterize π-stacking, as a result of charge penetration effects. The planar geometries of aromatic sp 2 carbon networks lead to attractive interactions that are “served up on a molecular pizza peel”, and adoption of slip-stacked geometries minimizes steric (rather than electrostatic) repulsion. The slip-stacked motif therefore emerges not as a defect induced by electrostatic repulsion but rather as a natural outcome of a conformational landscape that is dominated by van der Waals interactions (dispersion plus Pauli repulsion), and is therefore fundamentally quantum-mechanical in its origins. This reinterpretation of the forces responsible for π-stacking has important implications for the manner in which non-bonded interactions are modeled using classical force fields, and for rationalizing the prevalence of the slip-stacked π–π motif in protein crystal structures.
A Solid State Zwitterionic Plastic Crystal With High Static Dielectric Constant
Solid materials with a high dielectric constant have a wide range of applications in the energy storage field. In this research, an imidazolium-based zwitterion is designed, synthesized, and confirmed to have a plastic crystal phase based on the following experimental and computational evidence: (i) the presence of long-range order with weak intermolecular forces and competing attractive-repulsive interactions along different crystallographic directions; (ii) the observation of more than one endotherm on heating including a solid-solid phase transition at T s-s = −26 °C and melting of the plastic crystal at T m = 72°C; (iii) a low entropy of fusion at melting (2.1 JK −1 mol −1 ); (iv) a strongly anisotropic morphology; (v) relatively fast dynamics originating from short-range degrees of freedom. Furthermore, it exhibits a very high dielectric constant in the plastic crystal solid state (147 at −10°C and 103 at 70°C) due to the rotational degrees of freedom of plastic crystals that arise from weak net intermolecular interactions of zwitterions due to only having two carbons between the anion and cation. This material conveniently remains in the plastic crystal phase within 50 K of ambient temperature. This discovery opens new opportunities in the search for solid-state high dielectric constant materials.
Examination of Solvent- Room Temperature Ionic Liquid Interactions using 2D NMR and Computational Chemistry
Rare earth metals (REMs) are the key functioning material found in microelectronics, advanced batteries, and high-efficiency electrical motors. Current methods for producing REMs require high temperatures (i.e. energy intensive) and produce an abundance of hazardous fluoride waste. Room temperature ionic liquids (RTILs) are an ideal alternative electrolyte because they allow for a low-temperature, energy-efficient electrochemical process. RTIL properties such as viscosity and conductivity determine how effective they can be in this process. Since viscosity is directly related to both intra and intermolecular forces, this study sets out to investigate two RTILs-Trihexyltetradecylphosphonium dicyanamide (Cyphos 105) and Trihexyltetradecylphosphonium bistrifimide (Cyphos 109) and their interactions. Results from 2-dimensional nuclear magnetic resonance (2D NMR) spectroscopy show that there is an interaction between the tetradecyl tail and the phosphorus atom of the Cyphos RTIL. In order to establish whether this is an intra or an inter-molecular interaction, DFT-based computational calculations were performed to determine the optimized geometry. Calculations were done both in the gas phase and with propylene carbonate as an augmenting solvent. Preliminary results show that there is minimal intramolecular interaction between the tetradecyl tail and the phosphorus atom of the RTILs, and that the interaction seen in the 2D NMR experiment is likely of an intermolecular nature instead.
The Crystalline Behavior and Device Function of Nonfullerene Acceptors in Organic Solar Cells
The current research investigates the structure features and intermolecular interactions of nonfullerene acceptors (NFAs) in single crystal and thin films, as well as their solar cell applications. Guiding parameters and key intermolecular forces that lead to 2D brickwork or 3D web packing are identified. The atomic modification is shown as the key to induce hydrogen bonding or π–π stacking column, which results in different crystalline packing. The molecular assembly in thin film is initiated by hydrogen bonding and completed by π–π stacking reorganization. The packing energy is seen as a guiding parameter that dictates the NFA crystalline morphology in blended thin films. The crystalline packing motif is not directly related with device efficiency. However, the crystalline morphology is the key parameter to influence exciton/carrier dynamics and device performance. In conclusion, a broader picture on the scaling behavior of organic semiconductor crystals ranging from oligoacenes to NFAs is established.
Thin‐Film Organic Heteroepitaxy
Abstract Incorporating crystalline organic semiconductors into electronic devices requires understanding of heteroepitaxy given the ubiquity of heterojunctions in these devices. However, while rules for commensurate epitaxy of covalent or ionic inorganic material systems are known to be dictated by lattice matching constraints, rules for heteroepitaxy of molecular systems are still being written. Here, it is found that lattice matching alone is insufficient to achieve heteroepitaxy in molecular systems, owing to weak intermolecular forces that describe molecular crystals. It is found that, in addition, the lattice matched plane also must be the lowest energy surface of the adcrystal to achieve one‐to‐one commensurate molecular heteroepitaxy over a large area. Ultraviolet photoelectron spectroscopy demonstrates the lattice matched interface to be of higher electronic quality than a disordered interface of the same materials.
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.
Bridging adsorption behavior of confined CH 4 -CO 2 binary mixtures across scales
An accurate understanding of the competitive adsorption of CH 4 -CO 2 binary mixtures in nano-confined systems is critical for engineering CO 2 storage in shale gas reservoirs. Due to difficulties in making reliable experimental observations in nano-scale, atomistic simulations (ASs), such as the Grand Canonical Monte Carlo (GCMC) method, provide a viable approach to studying the adsorption behavior of confined fluids. ASs are, however, limited in the size of the compositional domain due to the high computational cost. This work proposes a framework that combines AS and the lattice Boltzmann (LB) method to bridge the physics of confined fluids across scales. The Peng–Robinson equation of state (PR-EoS) produces fugacity coefficients, which serve as input for conducting multi-component GCMC simulations. These GCMC simulations explore the competitive adsorption behavior of CH 4 -CO 2 in nano-slits at various composition, pressure, and channel-width conditions. Both components generate adsorption layers with high densities near the walls with CO 2 preferentially adsorbing compared to CH 4 on the organic walls of carbon sheets. At the mesoscale, a pseudopotential model represents the intermolecular forces in multi-component, multiple-relaxation-time LB simulations. The LB simulations are in good agreement with the GCMC results, allowing us to obtain values for tunable LB parameters. We then extend the use of LB to simulate adsorption behavior in complex networks with nano-sized channels. The phase behavior and fluid properties in the complex geometries of nano-channels differ from nano-slits and bulk systems. Furthermore, the bridging of physics from GCMC (microscale) to LB (mesoscale) via the macroscale PR-EoS connects the adsorption behavior of binary systems across scales.
Bridging confined phase behavior of CH 4 -CO 2 binary systems across scales
Phase behavior of confined fluids may deviate significantly from that of the bulk fluid due to the fluid-wall interactions being a significant portion of all intermolecular interactions under confinement. Despite recent advancements in understanding confined phase behavior of pure fluids, confined phase behavior of mixtures remains an understudied topic. In this work, we examine the confined phase behavior of a CH 4 -CO 2 binary system by combining Monte Carlo (MC) simulations, a cubic equation of state (EoS), and the lattice Boltzmann method (LBM). First, the effects of confinement on density and phase distribution in nano-size pores are established using Gibbs Ensemble MC calculations, which produce precise results of liquid and vapor confined pressures and account for the modification of the phase change location. By comparing the phase envelopes of bulk and confined mixtures at a fixed temperature, here it is observed that the phase envelopes shrink with reductions in pore size. Based on this observation, we extend a modified Peng-Robinson EoS, which was originally developed for pure fluids under confinement, to mixtures via van-der-Waals-type mixing rules and by accounting for shifts in the critical properties of confined CH 4 -CO 2 . The resulting phase envelopes are in good agreement with the MC data. In addition, a local density model is used in combination with the confined EoS to calculate adsorption isotherms of CH 4 -CO 2 mixtures and to characterize the behavior of confined matter in nanopores. Finally, we incorporate this EoS in a multicomponent multiphase LBM that uses a pseudopotential model to represent intermolecular forces. This workflow utilizes multiscale simulation techniques to bridge the behavior of multicomponent systems across scales and to shed light on the confined phase behavior of CH 4 -CO 2 binary systems.
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.
Microscopic Theory of Long-Time Center-of-Mass Self-Diffusion and Anomalous Transport in Ring Polymer Liquids
We construct a microscopic theory at the level of segment-scale correlated space–time intermolecular forces for the long-time center-of-mass (CM) diffusion constant and intermediate-time non-Fickian transport in dense solutions and melts of ring polymers. The approach combines ideas of polymer, colloid, and liquid-state statistical mechanics to quantify how the multifractal intra-ring conformational structure and inter-ring packing correlations determine dynamic caging constraints and time-dependent friction. Breakdown of Rouse theory is predicted to occur due to length scale-dependent temporal correlation of forces exerted on pairs of tagged ring segments from surrounding polymers. At large enough degrees of polymerization (N), a stronger scaling of the CM diffusion constant (D ∝ N –2 ) is predicted, with a crossover N D proportional to the product of the system-specific macromolecular volume fraction and dimensionless compressibility. In analogy with caging effects in glass-forming fluids, the theory does appear to begin to fail at sufficiently high N/N D for the center-of-mass diffusivity, likely due to another crossover to an even slower activated transport regime. However, use of N/N D with the theoretically predicted N D , in conjunction with dynamic blob scaling ideas for local Rouse friction, collapses semidilute and concentrated solution simulation data onto a master curve. Based on the same physical ideas employed to predict the diffusion constant, a generalized Langevin equation description is formulated for intermediate-time CM transport. It predicts two subdiffusive regimes that emerge due to the self-similar nature of internal ring structure. Analytic and numerical predictions for the apparent non-Fickian exponent as a function of time, N/N D and dimensionless compressibility, the evolution of the maximum degree of subdiffusive motion with N/N D , and the time scale for recovering Fickian diffusion are made, all of which are in good accordance with melt simulations. Furthermore, the present work sets the stage to address activated dynamics and glass formation on the macromolecular scale.
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.