SEARCH · Search NASA
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.
A coarse-graining approach for molecular simulation that retains the dynamics of the all-atom reference system by implementing hydrodynamic interactions
We report on a new approach for deriving coarse-grained intermolecular forces that retains the frictional contribution that is often discarded by conventional coarse-graining methods. The approach is tested for water and an aqueous glucose solution, and the results from the new implementation for coarse-grained molecular dynamics simulation show remarkable agreement with the dynamics obtained from reference all-atom simulations. The agreement between the structural properties observed in the coarse-grained and all-atom simulations is also preserved. We discuss how this approach may be applied broadly to any existing coarse-graining method where the coarse-grained models are rigorously derived from all-atom reference systems.
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
A quantitative analysis of weak intermolecular interactions & quantum chemical calculations (DFT) of novel chalcone derivatives
The novel chalcone derivatives have widespread applications in material science and medicinal industries. The density functional theory (DFT) is used to optimized the molecular structure of the three chalcone derivatives (M-I, II, III). The observed discrepancies between the theoretical and experimental (X-ray data) results attributed to different environments of the molecules, the experimental values are of the molecule in solid state there by subjected to the intermolecular forces, like non-bonded hydrogen bond interactions, where as isolated state in gas phase for theoretical studies. The lattice energy of all the molecules have been calculated using PIXELC module in Coulomb –London –Pauli (CLP) package and is partitioned into corresponding coulombic, polarization, dispersion and repulsion contributions. Lattice energy data confirm and strengthen the finding of the X-ray results that the weak but significant intermolecular interactions like C-H…O, Π- Π and C-H… Π plays an important role in the stabilization of crystal packing.
Static electric dipole polarizabilities of An{sup 5+/6+} and AnO{sub 2}{sup +/2+} (An = U, Np, and Pu) ions
The parallel components of static electric dipole polarizabilities have been calculated for the lowest lying spin-orbit states of the penta- and hexavalent oxidation states of the actinides (An) U, Np, and Pu, in both their atomic and molecular diyl ion forms (An{sup 5+/6+} and AnO{sub 2}{sup +/2+}) using the numerical finite-field technique within a four-component relativistic framework. The four-component Dirac-Hartree-Fock method formed the reference for MP2 and CCSD(T) calculations, while multireference Fock space coupled-cluster (FSCC), intermediate Hamiltonian Fock space coupled-cluster (IH-FSCC) and Kramers restricted configuration interaction (KRCI) methods were used to incorporate additional electron correlation. It is observed that electron correlation has significant (∼5 a.u.{sup 3}) impact upon the parallel component of the polarizabilities of the diyls. To the best of our knowledge, these quantities have not been previously reported and they can serve as reference values in the determination of various electronic and response properties (for example intermolecular forces, optical properties, etc.) relevant to the nuclear fuel cycle and material science applications. The highest quality numbers for the parallel components (α{sub zz}) of the polarizability for the lowest Ω levels corresponding to the ground electronic states are (in a.u.{sup 3}) 44.15 and 41.17 for UO{sub 2}{sup +} and UO{sub 2}{sup 2+}, respectively, 45.64 and 41.42 for NpO{sub 2}{sup +} and NpO{sub 2}{sup 2+}, respectively, and 47.15 for the PuO{sub 2}{sup +} ion.
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.
Energy transfer upon collision of selectively excited CO{sub 2} molecules: State-to-state cross sections and probabilities for modeling of atmospheres and gaseous flows
Carbon dioxide molecules can store and release tens of kcal/mol upon collisions, and such an energy transfer strongly influences the energy disposal and the chemical processes in gases under the extreme conditions typical of plasmas and hypersonic flows. Moreover, the energy transfer involving CO{sub 2} characterizes the global dynamics of the Earth-atmosphere system and the energy balance of other planetary atmospheres. Contemporary developments in kinetic modeling of gaseous mixtures are connected to progress in the description of the energy transfer, and, in particular, the attempts to include non-equilibrium effects require to consider state-specific energy exchanges. A systematic study of the state-to-state vibrational energy transfer in CO{sub 2} + CO{sub 2} collisions is the focus of the present work, aided by a theoretical and computational tool based on quasiclassical trajectory simulations and an accurate full-dimension model of the intermolecular interactions. In this model, the accuracy of the description of the intermolecular forces (that determine the probability of energy transfer in molecular collisions) is enhanced by explicit account of the specific effects of the distortion of the CO{sub 2} structure due to vibrations. Results show that these effects are important for the energy transfer probabilities. Moreover, the role of rotational and vibrational degrees of freedom is found to be dominant in the energy exchange, while the average contribution of translations, under the temperature and energy conditions considered, is negligible. Remarkable is the fact that the intramolecular energy transfer only involves stretching and bending, unless one of the colliding molecules has an initial symmetric stretching quantum number greater than a threshold value estimated to be equal to 7.
The effect of surfaces on molecular ordering in thin liquid-crystal systems
A theoretical method for analyzing the interplay between pair long-range intermolecular forces and nonlocal, relatively short-range, surface interactions in liquid crystals, confined between plates of thin planar cells, is developed. It is shown that this method, as involving the concept of local orientational and translational order parameters, enables detailed investigations of the emergence of smectic A, nematic, and isotopic phases, as well as yields an insight into phase transitions between them, in cases of systems strongly affected by surfaces. The evidence of various surface effects, including the coexistence of different phases and the inward propagation of surface melting under the increase of temperature, is also given. The underlying numerical procedure, based on the algorithm of self-consistent calculations of local order parameters, is found to be very effective, allowing one to consider model systems of rather large thicknesses, corresponding to thicknesses of real sample cells.
Generalized and efficient algorithm for computing multipole energies and gradients based on Cartesian tensors
Accurate representation of intermolecular forces has been the central task of classical atomic simulations, known as molecular mechanics. Recent advancements in molecular mechanics models have put forward the explicit representation of permanent and/or induced electric multipole (EMP) moments. The formulas developed so far to calculate EMP interactions tend to have complicated expressions, especially in Cartesian coordinates, which can only be applied to a specific kernel potential function. For example, one needs to develop a new formula each time a new kernel function is encountered. The complication of these formalisms arises from an intriguing and yet obscured mathematical relation between the kernel functions and the gradient operators. Here, I uncover this relation via rigorous derivation and find that the formula to calculate EMP interactions is basically invariant to the potential kernel functions as long as they are of the form f(r), i.e., any Green’s function that depends on inter-particle distance. I provide an algorithm for efficient evaluation of EMP interaction energies, forces, and torques for any kernel f(r) up to any arbitrary rank of EMP moments in Cartesian coordinates. The working equations of this algorithm are essentially the same for any kernel f(r). Recently, a few recursive algorithms were proposed to calculate EMP interactions. Depending on the kernel functions, the algorithm here is about 4–16 times faster than these algorithms in terms of the required number of floating point operations and is much more memory efficient. I show that it is even faster than a theoretically ideal recursion scheme, i.e., one that requires 1 floating point multiplication and 1 addition per recursion step. This algorithm has a compact vector-based expression that is optimal for computer programming. The Cartesian nature of this algorithm makes it fit easily into modern molecular simulation packages as compared with spherical coordinate-based algorithms. A software library based on this algorithm has been implemented in C++11 and has been released.
Flexible all-carbon photovoltaics with improved thermal stability
The structurally robust nature of nanocarbon allotropes, e.g., semiconducting single-walled carbon nanotubes (SWCNTs) and C{sub 60}s, makes them tantalizing candidates for thermally stable and mechanically flexible photovoltaic applications. However, C{sub 60}s rapidly dissociate away from the basal of SWCNTs under thermal stimuli as a result of weak intermolecular forces that “lock up” the binary assemblies. Here, we explore use of graphene nanoribbons (GNRs) as geometrically tailored protecting layers to suppress the unwanted dissociation of C{sub 60}s. The underlying mechanisms are explained using a combination of molecular dynamics simulations and transition state theory, revealing the temperature dependent disassociation of C{sub 60}s from the SWCNT basal plane. Our strategy provides fundamental guidelines for integrating all-carbon based nano-p/n junctions with optimized structural and thermal stability. External quantum efficiency and output current–voltage characteristics are used to experimentally quantify the effectiveness of GNR membranes under high temperature annealing. Further, the resulting C{sub 60}:SWCNT:GNR ternary composites display excellent mechanical stability, even after iterative bending tests. - Graphical abstract: The incorporation of solvent resistant, mechanically flexible and electrically addressable 2-D soft graphene nanoribbons facilitates the assembly of photoconductive carbon nano-p/n junctions for thermally stable and flexible photovoltaic cells.
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.
The Odd Power of Dispersion
In ancient China, water has been regarded as one of the five vital components of life. It has been observed that water has many fascinating properties: water is ‘soft’ yet it can penetrate a hard rock; water is ‘pure’ yet it can tolerate other beings. Because of its unique properties, water is often associated with good quality and has been given the highest praise by Laozi in his book Tao Te Ching saying: the highest/best quality that one can have is being like water. However, little did people understand why and how water possesses such fascinating properties. Modern scientific developments made people realize that the macroscopic liquid water is made of a large number of water molecules held together via a network of hydrogen bonds. And those wonderful properties of water are merely the macroscopic manifestations of the interactions between water molecules and other molecules. For example, the dissolving ability of water is due to the fact that the interaction between a water molecule and the other molecular species is stronger than the interactions among their own molecular species. In fact the interactions between any two molecules are governed by the same physics and are termed intermolecular interaction (or intermolecular forces in some literature, although technically ‘force’ is incorrect usage here). Although the very existence of the intermolecular interactions is easily proved, e.g. the mere presence of the solid phase of matter, and scientists today have recognized that the seemingly weak intermolecular interactions essentially hold the world together through a delicate and cooperative process, the theoretical understanding of various intermolecular interactions is still far from satisfactory. On the practical side, theoreticians need to balance computational cost and accuracy. Because of the relatively small magnitudes of the intermolecular interactions, errors that appear tiny compared to the usual chemical (covalent) bonding may change conclusions qualitatively. High-level ab initio methods including explicit description of electron correlation can achieve the desired accuracy at very high computational cost. (Chapter 5 and 6) However the cooperative network of hundreds of thousands of molecules that reflects the true power of intermolecular interactions cannot be modeled easily by ab initio methods. Deeper understanding of intermolecular interactions yields better theoretical models; better theoretical models facilitate and even deepen the understanding of intermolecular interactions. With the aforementioned motivation in mind, a significant portion of this dissertation is dedicated to developing a method to describe the intermolecular interactions accurately with affordable computational resources.
Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab initio water at near-zero added cost
It is now established that nuclear quantum motion plays an important role in determining water's hydrogen bonding, structure, and dynamics. Such effects are important to include in density functional theory (DFT) based molecular dynamics simulation of water. The standard way of treating nuclear quantum effects, path integral molecular dynamics (PIMD), multiplies the number of energy/force calculations by the number of beads required. In this work we introduce a method whereby PIMD can be incorporated into a DFT simulation with little extra cost and little loss in accuracy. The method is based on the many body expansion of the energy and has the benefit of including a monomer level correction to the DFT energy. Our method calculates intramolecular forces using the highly accurate monomer potential energy surface developed by Partridge–Schwenke, which is cheap to evaluate. Intermolecular forces and energies are calculated with DFT only once per timestep using the centroid positions. We show how our method may be used in conjunction with a multiple time step algorithm for an additional speedup and how it relates to ring polymer contraction and other schemes that have been introduced recently to speed up PIMD simulations. We show that our method, which we call “monomer PIMD”, correctly captures changes in the structure of water found in a full PIMD simulation but at much lower computational cost.
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.
Concerted Mitigation of O···H and C(π)···H Interactions Prospects Sixfold Gain in Optical Nonlinearity of Ionic Stilbazolium Derivatives
DAST (4-dimethylamino-N-methyl-4-stilbazolium tosylate) is the most commercially successful organic nonlinear optical (NLO) material for frequency-doubling, integrated optics, and THz wave applications. Its success is predicated on its high optical nonlinearity with concurrent sufficient thermal stability. Many chemical derivatives of DAST have therefore been developed to optimize their properties; yet, to date, none have surpassed the overall superiority of DAST for NLO photonic applications. This is perhaps because DAST is an ionic salt wherein its NLO-active cation is influenced by multiple types of subtle intermolecular forces that are hard to quantify, thus, making difficult the molecular engineering of better functioning DAST derivatives. Here, we establish a model parameter, ηinter, that isolates the influence of intermolecular interactions on second-order optical nonlinearity in DAST and its derivatives, using second-harmonic generation (SHG) as a qualifier; by systematically mapping intercorrelations of all possible pairs of intermolecular interactions to ηinter, we uncover a relationship between concerted intermolecular interactions and SHG output. This correlation reveals that a sixfold gain in the intrinsic second-order NLO performance of DAST is possible, by eliminating the identified interactions. This prediction offers the first opportunity to systematically design next-generation DAST-based photonic device nanotechnology to realize such a prospect.
Atomization and dense-fluid breakup regimes in liquid rocket engines
Until recently, modern theory has lacked a fundamentally based model to predict the operating pressures where classical sprays transition to dense-fluid mixing with diminished surface tension. In this paper, such a model is presented to quantify this transition for liquid-oxygen–hydrogen and n-decane–gaseous-oxygen injection processes. The analysis reveals that respective molecular interfaces break down not necessarily because of vanishing surface tension forces but instead because of the combination of broadened interfaces and a reduction in mean free molecular path. When this occurs, the interfacial structure itself enters the continuum regime, where transport processes rather than intermolecular forces dominate. Using this model, regime diagrams for the respective systems are constructed that show the range of operating pressures and temperatures where this transition occurs. The analysis also reveals the conditions where classical spray dynamics persists even at high supercritical pressures. As a result, it demonstrates that, depending on the composition and temperature of the injected fluids, the injection process can exhibit either classical spray atomization, dense-fluid diffusion-dominated mixing, or supercritical mixing phenomena at chamber pressures encountered in state-of-the-art liquid rocket engines.