Engineering Papers⌕ Search

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.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1971

Shock wave structure using nonlinear model Boltzmann equations

The structure of a strong plane shock wave in a monatomic rarefied perfect gas is one of the simplest problems able to be posed in kinetic theory, and one of the hardest to solve. Its simplicity lies in the absence of solid boundaries, geometrical complications, or internal molecular energy. Its difficulty arises from the great departure of the gas from equilibrium within the shock, which invalidates many of the techniques used successfully elsewhere in kinetic theory. In addition to this theoretical challenge, the modern development of ballistics and hypersonic flight has helped to stimulate extensive theoretical and experimental interest in the shock problem. The experimenters in turn have encountered great difficulties on account of the very small physical dimensions of shocks. In fact, until very recently indeed, any close comparisons of theoretical and experimental shock structure results have been rather unprofitable due to the inadequacies of both theory and experiment. During the last few years this situation has been appreciably improved by development of the Monte Carlo method. This allows idealized 'experiments' to be performed on large computers instead of in wind tunnels, using a known intermolecular force law. The most developed of these methods has been shown to be equivalent theoretically to the Boltzmann equation and to give results which agree extremely closely with measurements of high accuracy. Thus Monte Carlo results not only form the soundest basis for our present theoretical knowledge of shock wave structure, but, for purposes of developing other theories, can also be considered a very valuable experimental resource. However, such results remain very expensive to obtain. In this thesis we develop more economical kinetic theory methods for the approximate prediction of shock structure, and compare our results with those of the Monte Carlo method.

Segal, Ben Maurice↗
NASA NTRS1996

Cross-Effects in Microgravity Flows

Film growth by chemical/physical vapor deposition is a process of considerable interest in microgravity experiments. The absence of natural convection should allow better control of film growth processes but, in highly non-isothermal ampoules, thermal slip (creep) can become a matter of significant concern. The reported research is a theoretical and experimental investigation of the flow of gas/vapor mixtures under non-continuum conditions. The Boltzmann equation has been solved for a monatomic gas under non-condensing conditions and the various phenomenological coefficients have been computed. Computations for realistic potentials as well as for velocity and creep slip have been completed and the creep slip has been found to be dependent on the type of gas confirming the accuracy of previous variational results. The variational technique has been extended and planar flows calculated via the Burnett solutions. Velocity, diffusion and creep slips have been computed for gas mixtures and previously unknown dependencies of the creep slip on the mixture properties have been observed. Also for gas mixtures, an integral representation of the linearized Boltzmann operator has been developed for use in numerical and variational calculations for all intermolecular force laws. Two, two-bulb capillary systems have been designed, built and tested for the measurements of cross-flows; one of glass for isothermal measurements and one of stainless steel for non-isothermal measurements. Extensive data have been collected for Ar-He and N2-He mixtures at a variety of pressures and mole ratios. Viscosity, velocity slip coefficients and tangential momentum accommodation coefficients have been obtained from measurements with a spinning rotor gauge via a new theory that has been formulated for the spinning rotor gauge in the slip regime. The FIDAP fluid dynamics code has been applied to condensing flows in ampoules in the continuum regime and agreement obtained with the earlier work of Duval.

Loyalka, Sudarshan K.↗
NASA NTRS2001

Viscosity Coefficient Curve Fits for Ionized Gas Species Grant Palmer

Viscosity coefficient curve fits for neutral gas species are available from many sources. Many do a good job of reproducing experimental and computational chemistry data. The curve fits are usually expressed as a function of temperature only. This is consistent with the governing equations used to derive an expression for the neutral species viscosity coefficient. Ionized species pose a more complicated problem. They are subject to electrostatic as well as intermolecular forces. The electrostatic forces are affected by a shielding phenomenon where electrons shield the electrostatic forces of positively charged ions beyond a certain distance. The viscosity coefficient for an ionized gas species is a function of both temperature and local electron number density. Currently available curve fits for ionized gas species, such as those presented by Gupta/Yos, are a function of temperature only. What they did was to assume an electron number density. The problem is that the electron number density they assumed was unrealistically high. The purpose of this paper is two-fold. First, the proper expression for determining the viscosity coefficient of an ionized species as a function of both temperature and electron number density will be presented. Then curve fit coefficients will be developed using the more realistic assumption of an equilibrium electron number density. The results will be compared against previous curve fits and against highly accurate computational chemistry data.

Palmer, Grant↗
NASA NTRS2004

Theoretical and Experimental Investigation of the Stability of an Evaporating Constrained Vapor Bubble

The major accomplishments of the experimental portion of the research were documented in Ling Zheng's doctoral dissertation. Using Pentane, he obtained a considerable amount of data on the stability and heat transfer characteristics of an evaporating meniscus. The important points are that experimental equipment to obtain data on the stability and heat transfer characteristics of an evaporating meniscus were built and successfully operated. The data and subsequent analyses were accepted by the Journal of Heat Transfer for publication in 2004 [PU4]. The work was continued by a new graduate student using HFE-7000 [PU3] and then Pentane at lower heat fluxes. The Pentane results are being analyzed for publication. The experimental techniques are currently being used in our other NASA Grant. The oscillation of the contact line observed in the experiments involves evaporation (retraction part) and spreading. Since both processes occur with finite contact angles, it is important to derive a precise equation of the intermolecular forces (disjoining pressure) valid for non-zero contact angles. This theoretical derivation was accepted for publication by Journal of Fluid Mechanics [PU5]. The evaporation process near the contact line is complicated, and an idealized micro heat pipe has been proposed to help in elucidating the detailed evaporation process [manuscripts in preparation].

Wayner, P. C., Jr.↗
NASA NTRS2004

Microscale Investigation of Thermo-Fluid Transport in the Transition FIL, Region of an Evaporating Capillary Meniscus Using a Microgravity Environment

In order to enhance the fundamental understanding of thin film evaporation and thereby improve the critical design concept for two-phase heat transfer devices, microscale heat and mass transport is to be investigated for the transition film region using state-of-the-art optical diagnostic techniques. By utilizing a microgravity environment, the length scales of the transition film region can be extended sufficiently, from submicron to micron, to probe and measure the microscale transport fields which are affected by intermolecular forces. Extension of the thin film dimensions under microgravity will be achieved by using a conical evaporator made of a thin silicon substrate under which concentric and individually controlled micro-heaters are vapor-deposited to maintain either a constant surface temperature or a controlled temperature variation. Local heat transfer rates, required to maintain the desired wall temperature boundary condition, will be measured and recorded by the concentric thermoresistance heaters controlled by a Wheatstone bridge circuit, The proposed experiment employs a novel technique to maintain a constant liquid volume and liquid pressure in the capillary region of the evaporating meniscus so as to maintain quasi-stationary conditions during measurements on the transition film region. Alternating use of Fizeau interferometry via white and monochromatic light sources will measure the thin film slope and thickness variation, respectively. Molecular Fluorescence Tracking Velocimetry (MFTV), utilizing caged fluorophores of approximately 10-nm in size as seeding particles, will be used to measure the velocity profiles in the thin film region. An optical sectioning technique using confocal microscopy will allow submicron depthwise resolution for the velocity measurements within the film for thicknesses on the order of a few microns. Digital analysis of the fluorescence image-displacement PDFs, as described in the main proposal, can further enhance the depthwise resolution.

Kihm, K. D.↗
DOE OSTI2020

On the origin of ground-state vacuum-field catalysis: Equilibrium consideration

Recent experiments suggest that vibrational strong coupling (VSC) may significantly modify ground-state chemical reactions and their rates even without external pumping. The intrinsic mechanism of this “vacuum-field catalysis” remains largely unclear. Generally, modifications of thermal reactions in the ground electronic states can be caused by equilibrium or non-equilibrium effects. The former are associated with modifications of the reactant equilibrium distribution as expressed by the transition state theory of chemical reaction rates, while the latter stem from the dynamics of reaching and leaving transition state configurations. In this work, we examine how VSC can affect chemical reactions rates in a cavity environment according to transition state theory. Our approach is to examine the effect of coupling to cavity mode(s) on the potential of mean force (PMF) associated with the reaction coordinate. Within the context of classical nuclei and classical photons and also assuming no charge overlap between molecules, we find that while the PMF can be affected by the cavity environment, this effect is negligible for the usual micron-length cavities used to examine VSC situations.

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

A deep potential model with long-range electrostatic interactions

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

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Proton radiation effects in indium oxide using cascade molecular dynamics simulations

Metal oxide (MO) semiconductors, characterized by their wide band gaps and notable charge transport properties, are promising candidates for electronic applications in extreme environments, including near-Earth space. However, atomistic simulations of radiation–matter interactions in MOs remain challenging due to the limitations of existing interatomic potentials, which often fail to capture both the short-range repulsive forces essential for radiation damage modeling and the long-range electrostatic effects governing defect evolution. In this work, we develop a customized interatomic potential tailored for radiation damage simulations in indium oxide (In 2 O 3 ) as a model system, a representative MO material. Our potential integrates the Ziegler-Biersack-Littmark potential to accurately describe short-range interactions with Buckingham and Coulombic potentials to account for long-range forces. We perform molecular dynamics simulations of low-energy proton irradiation using this custom potential. We employ the primary knock-on atom (PKA) cascade method to study atomic displacements and primary defect formation. Simulations were conducted for 1 keV proton irradiation in four randomly chosen directions, and PKA-driven defect analyses at 5, 10, and 15 keV to examine the effects of direction and energy level on damage generation. Our results provide insight into the impact of irradiation direction and energy level on the cascade evolution and defect formation mechanisms.

Atomistic simulations↗
DOE OSTI2020

Examining the order-of-limits problem and lattice constant performance of the Tao–Mo functional

In their recent communication a semi-local density functional derived from the density matrix expansion of the exchange hole localised by a general coordinate transformation. We show that the order-of-limits problem present in the functional, dismissed as harmless in the original publication, causes severe errors in predicted phase transition pressures. We also show that the claim that lattice volume prediction accuracy exceeds that of existing similar functionals was based on comparison to reference data that misses anharmonic zero-point expansion and consequently overestimates accuracy. Here, by highlighting these omissions, we give a more accurate assessment of the Tao-Mo functional and show a possible direction for resolving the order-of-limits problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Machine learning Frenkel Hamiltonian parameters to accelerate simulations of exciton dynamics

In this manuscript, we develop multiple machine learning (ML) models to accelerate a scheme for parameterizing site-based models of exciton dynamics from all-atom configurations of condensed phase sexithiophene systems. This scheme encodes the details of a system’s specific molecular morphology in the correlated distributions of model parameters through the analysis of many single-molecule excited-state electronic-structure calculations. These calculations yield excitation energies for each molecule in the system and the network of pair-wise intermolecular electronic couplings. In this work, we demonstrate that the excitation energies can be accurately predicted using a kernel ridge regression (KRR) model with Coulomb matrix featurization. We present two ML models for predicting intermolecular couplings. The first one utilizes a deep neural network and bi-molecular featurization to predict the coupling directly, which we find to perform poorly. The second one utilizes a KRR model to predict unimolecular transition densities, which can subsequently be analyzed to compute the coupling. We find that the latter approach performs excellently, indicating that an effective, generalizable strategy for predicting simple bimolecular properties is through the indirect application of ML to predict higher-order unimolecular properties. Such an approach necessitates a much smaller feature space and can incorporate the insight of well-established molecular physics.

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

Ligand structure and adsorption free energy of nanocrystals on solid substrates

In this paper, we present an investigation on the absorption of alkylthiolated nanocrystals on a solid substrate. We calculate adsorption free energies and report a number of effects induced by the substrate. Nearest neighbor distances and bonding free energies are significantly different than for a free floating case, there is a weakening of bonding free energies among nanocrystals, and the adsorption is manifestly anisotropic, i.e., stronger along certain directions of the nanocrystal core. We contend that this last result accounts for the Bain transition (fcc → bcc) observed in experimental results. We report the presence of vortices induced by the substrate, which explain the increased nearest neighbor distance among nanocrystals, which is in excellent quantitative agreement with experimental results and with the predictions of the Orbifold Topological Model. Implications for the assembly of nanostructures and future experiments are also discussed.

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

Phonon gas model for thermal conductivity of dense, strongly interacting liquids

Developing predictive thermal property models for liquids based on microscopic principles has been elusive. The difficulty is that liquids have gas-like and solid-like attributes that are at odds when considering the frameworks of microscopic models: Models for gases are simple due to randomness and low density, whereas models for crystalline solids rely on symmetry and long-range order for easier calculation. The short-range order in liquids does, however, provide structure to neighboring molecules similar to amorphous solids, and there have been recent advances indicating that collective vibrational modes store heat in liquids. Models combining Debye approximations from solid-state physics and Frenkel’s theory of liquids can accurately predict the heat capacity of liquids. Phonon-like dispersions in liquids have also been widely observed in neutron scattering experiments. These developments motivate us to propose a model where high-frequency vibrational modes, which travel at the speed of sound and have a mean free path on the order of the average intermolecular distance, conduct heat in liquids. We use this liquid phonon gas model to calculate the thermal conductivity of liquids with varying intermolecular interaction energies from strongest to weakest—Coulomb, hydrogen-bonding, Keesom, and London dispersion energy. Generally, the model is more accurate as the intermolecular interaction energy and density of liquids increase. The calculated thermal conductivity of Coulombic-bound molten sodium nitrate and hydrogen-bonded water is within 1.46% and 2.98% of the experimentally measured values, respectively, across their entire temperature ranges. Further modal analysis of the velocity and the mean free path of collective vibrations could establish the liquid phonon gas model as an accurate model for weakly interacting liquids as well.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quadrupole-mediated dielectric response and the charge-asymmetric solvation of ions in water

Treating water as a linearly responding dielectric continuum on molecular length scales allows very simple estimates of the solvation structure and thermodynamics for charged and polar solutes. While this approach can successfully account for basic length and energy scales of ion solvation, computer simulations indicate not only its quantitative inaccuracies but also its inability to capture some basic and important aspects of microscopic polarization response. Here, we consider one such shortcoming, a failure to distinguish the solvation thermodynamics of cations from that of otherwise-identical anions, and we pursue a simple, physically inspired modification of the dielectric continuum model to address it. The adaptation is motivated by analyzing the orientational response of an isolated water molecule whose dipole is rigidly constrained. Its free energy suggests a Hamiltonian for dipole fluctuations that accounts implicitly for the influence of higher-order multipole moments while respecting constraints of molecular geometry. Finally, we propose a field theory with the suggested form, whose nonlinear response breaks the charge symmetry of ion solvation. An approximate variational solution of this theory, with a single adjustable parameter, yields solvation free energies that agree closely with simulation results over a considerable range of solute size and charge.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 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 OSTI2022

Polariton relaxation under vibrational strong coupling: Comparing cavity molecular dynamics simulations against Fermi’s golden rule rate

Under vibrational strong coupling (VSC), the formation of molecular polaritons may significantly modify the photo-induced or thermal properties of molecules. In an effort to understand these intriguing modifications, both experimental and theoretical studies have focused on the ultrafast dynamics of vibrational polaritons. Here, following our recent work [Li et al., J. Chem. Phys. 154, 094124 (2021)], we systematically study the mechanism of polariton relaxation for liquid CO 2 under a weak external pumping. Classical cavity molecular dynamics (CavMD) simulations confirm that polariton relaxation results from the combined effects of (i) cavity loss through the photonic component and (ii) dephasing of the bright-mode component to vibrational dark modes as mediated by intermolecular interactions. The latter polaritonic dephasing rate is proportional to the product of the weight of the bright mode in the polariton wave function and the spectral overlap between the polariton and dark modes. Both these factors are sensitive to parameters such as the Rabi splitting and cavity mode detuning. Compared to a Fermi’s golden rule calculation based on a tight-binding harmonic model, CavMD yields a similar parameter dependence for the upper polariton relaxation lifetime but sometimes a modest disagreement for the lower polariton. We suggest that this disagreement results from polariton-enhanced molecular nonlinear absorption due to molecular anharmonicity, which is not included in our analytical model. Finally, we also summarize recent progress on probing nonreactive VSC dynamics with CavMD.

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

Density functional theory of water with the machine-learned DM21 functional

The delicate interplay between functional-driven and density-driven errors in density functional theory (DFT) has hindered traditional density functional approximations (DFAs) from providing an accurate description of water for over 30 years. Recently, the deep-learned DeepMind 21 (DM21) functional has been shown to overcome the limitations of traditional DFAs as it is free of delocalization error. To determine if DM21 can enable a molecular-level description of the physical properties of aqueous systems within Kohn–Sham DFT, we assess the accuracy of the DM21 functional for neutral, protonated, and deprotonated water clusters. Furthermore, we find that the ability of DM21 to accurately predict the energetics of aqueous clusters varies significantly with cluster size. Additionally, we introduce the many-body MB-DM21 potential derived from DM21 data within the many-body expansion of the energy and use it in simulations of liquid water as a function of temperature at ambient pressure. We find that size-dependent functional-driven errors identified in the analysis of the energetics of small clusters calculated with the DM21 functional result in the MB-DM21 potential systematically overestimating the hydrogen-bond strength and, consequently, predicting a more ice-like local structure of water at room temperature.

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

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

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

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗