SEARCH · Search NASA
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.
Morse potential function determined from viscosity and second virial coefficient, used to approximate pair interaction potential for hydrocarbons
Morse potential function used to approximate pair interaction potential for wide variety of hydrocarbons determined from viscosity and second virial coefficient
Morse potential functions for nonpolar gases determined from experimental viscosity and second virial coefficient
Computer program for calculation of transport collision integrals with tables for Morse potential
Computer program for calculating classical and quantal JWKB Chapman-Enskog transport collision integrals with application to Morse potential
Collision integrals, effective cross sections, and angles of deviation for Morse potential
Relationship between Pitzer acentric factor and Morse intermolecular potential
Pitzer acentric factor relationship to Morse intermolecular potential function appears to hold only for normal fluids
The desorption of O/CO from graphitic carbon surfaces is investigated using a one-dimensional model describing the adsorbate interactions with the surface phonon bath. The kinetics of desorption are described through the solution of a master equation for the time-dependent population of the adsorbate in an oscillator state, which is modified through thermal fluctuations at the surface. The interaction of the adsorbate with the surface phonons is explicitly captured by using the computed phonon density Of states (PDOS) of the surface. The coupling of the adsorbate with the phonon bath results in the transition of the adsorbate up and down a vibrational ladder. The adsorbate-surface interaction is represented in the model using a Morse potential, which allows for the desorption process to be directly modeled as a transition from bound to free (continuum) state. The PDOS is a property of the material and the lattice; and is highly sensitive to the presence of defects. The effect of etch pits along with random surface defects on the PDOS is considered in the present work. The presence of defects causes a redshift and broadening of the PDOS, which in turn changes the phonon frequency modes available for adsorbate coupling at the surface. Using the realistic PDOS distributions, the phonon-induced desorption (PID) model was used to compute the transition and desorption rates for both pristine and defective systems. Mathissen’s rule is used to compute the phonon relaxation time for pristine and defective systems based on the phonon scattering times for each of the different scattering processes. First, the desorption rates of the pristine system is fitted against the experimental values to obtain the Morse potential parameters for each of the observed adatoms. These Morse potential parameters are used along with the defective PDOS and phonon relaxation time to compute the desorption rates for the defective system. The defective system rates (both transition and desorption) were consistently lower in comparison with the pristine system. The difference between the transition rates is more significant at lower initial states due to higher energy spacing between the levels. In the case of the desorption rates, the difference between the defective and pristine system is more significant at higher temperatures. The desorption rates for each of the system shows an order of magnitude decrease with the strongly bound systems exhibiting the greatest reduction in the desorption rates.
In the field of computational chemistry, predicting bond dissociation energies (BDEs) presents well-known challenges, particularly due to the multireference character of reactive systems. Many chemical reactions involve configurations where single-reference methods fall short, as the electronic structure can significantly change during bond breaking. As generating training data for partially broken bonds is a challenging task, even state-of-the-art reactive machine learning interatomic potentials (MLIPs) often fail to predict reliable BDEs and smooth dissociation curves. By contrast, simple and inexpensive physics-based models, such as the well-established Morse potential, do not suffer from any such limitations. This work leverages the Morse potential to improve reactive MLIPs by augmenting the training data set with inexpensive Morse data along the dissociation pathways. Further, this physics-constrained data augmentation (PCDA) approach results in MLIPs with smooth bond dissociation curves as well as near coupled-cluster level BDEs, all without requiring any expensive multireference quantum mechanical calculations. A case study for methane combustion demonstrates how the PCDA approach can improve an existing reactive MLIP, namely, ANI-1xnr. In conclusion, not only are the BDEs and bond dissociation curves for all radicals and molecules significantly improved compared to ANI-1xnr but the PCDA-trained MLIP retains the reliability of ANI-1xnr when performing reactive molecular dynamics simulations.
Ultrasonic‐velocity dispersion measurements have been performed in parahydrogen and its mixtures with helium, neon, and argon, all at 300°K. In each case the experimental dispersion curves can be matched successfully to those calculated under the assumption that the 0–2 rotational transition relaxes separately from the 2–4 and higher‐order terms. For pure pH 2 we find a relaxation time 𝜏 20 of 1.30 × 10 −8 sec for the 2 → 0 transition and a 𝜏 42 of 3.90 × 10 −8 sec for the 4 → 2 transition. Comparison with the quantum‐mechanical theories of Roberts and of Davison for H 2 –H 2 collisions using Morse potentials shows good agreement for 𝜏 20 over the temperature range of 75°–300°K. The Morse‐potential asymmetry parameter yielding the best fit is β = 0.113 for Roberts' calculation and 0.108 for Davison's. It is found that He–pH 2 collisions are more effective than pH 2 –pH 2 in producing the J = 0 to J = 2 transition, but less effective for higher‐order transitions. Collisions of neon with pH 2 are found to be more effective at room temperature for inducing the 0 → 2 and 2 → 4 transitions than either helium or argon.
Theoretical calculation of the ground state properties of spin-alined atomic hydrogen by the Monte Carlo method. The interatomic interaction, as described by the results of Kolos and Wolniewicz (1965), is made fit to a Morse potential form. An appropriate trail wavefunction is formed from the short-range part of the WKB solution for a pair of atoms interacting through a Morse potential.
A generic method is developed to investigate the fracture behavior of an epoxy resin using atomistic molecular (MD) dynamics simulations. The epoxy system consisted of the stoichiometric mixture of tetra-/tri-functional epoxies and di-functional hardener molecules and was cured with a step-growth crosslinking algorithm. A hybrid force field (hFF) in which Morse bond function were added in second generation-general Amber Force Field (GAFF2) was used to capture bond breaking of the crosslinked epoxy system under uniaxial deformation. The Morse bond parameters for covalent bonds in the crosslinked backbone were fit to dissociation curves computed by CASPT2/6-311+G**. After the systems in various crosslink densities were obtained and equilibrated, the full length of stress-strain (𝜎-𝜀) curve was produced including initial elastic regime, yielding, plastic flow, strain hardening, and progressive failure. We found that the hFF was effective in studying the dissociation of crosslinked polymers and understanding the failure mechanism. In the results, the hFF delivered the same elastic property obtained from the unmodified force field GAFF. The stress-strain curve, however, began to be distinct in the plastic regime due to broken covalent bonds capable of altering the tensile behavior. Once the stress reached maximum, it progressively decreased until the ultimate failure (𝜎=0). Since hFF employed the existing GAFF functions with including a minimal change of its bond potential, this approach is easier to perform and computationally more efficient for unveiling the fracture behavior of polymer materials at the molecular level than other approaches with reactive force fields. We expect this approach will be utilized to accelerate the material-by-design process for thermosets by incorporate data from molecular models.
Quantum spectral simulations of the yellow excimer emission band of CsXe are presented. Synthetic spectra as a function of wave number are calculated for the 2 Sigma 1/2 + (7s) - 2 Sigma 1/2 + (6s) transition by the use of the equation of Tellinghuisen et al. (1976) with a theoretical potential for the ground state and a Morse potential curve with an electron frequency of 32/cm for the excited state. Results based on emission studies at 450 K and 200 and 800 torr, are found to be consistent with absorption studies. The undulatory structure observed in the spectrum is attributed not to the vibrational spacing in the excited state, but rather to a characteristic reflection structure associated with nearly parallel upper and lower potential curves.
A time dependent wave packet method is presented for the rapid calculation of the properties of systems in thermal equilibrium and is applied, as an illustration, to electronic spectra. The thawed Gaussian approximation to quantum wave packet dynamics combined with evaluation of the density matrix operator by imaginary time propagation is shown to give exact electronic spectra for harmonic potentials and excellent results for both a Morse potential and for the band contours of the three transitions of the visible electronic absorption spectrum of the iodine molecule. The method, in principle, can be extended to many atoms (e.g., condensed phases) and to other properties (e.g., infrared and Raman spectra and thermodynamic variables).
First measurements of electron impact excitation on the O2 a delta-one sub g metastable state to a higher bound state, pi-one sub u, is reported. The data were taken at an incident electron energy of 500 eV and scattering angles of 0 deg to 15 deg. A superelastic electron scattering experiment was performed on the pure oxygen target gas in order to establish the composition of the discharged oxygen and to confirm the existence of delta-one sub g molecules in the interaction region. The electron-impact energy loss spectrum obtained is shown. The electronic assignment of the upper state which produces the vibrational structure was determined to be pi-one sub u by examining the potential energy curves of O2 and applying the electric dipole selection rules. The relative line intensities were obtained and used to determine the Morse potential for the pi-one sub u state.
The activation energies for diffusion were determined for gold, platinum and iridium adatoms on plane and plane PT surfaces and were found to be in good agreement with the measurements reported by Bassett and Webber. The Lennard-Jones pair potentials were used to model the interatomic forces, and relaxation of the substrate atoms in near proximity to the adatom was considered in detail. The present calculations clarify the mechanism of the observed two-dimensional diffusion of platinum and iridium atoms on a plane PT surface. The results are compared with those obtained using Morse potential functions and different relaxation techniques.
While the low thermal conductivities of silica aerogels have made them of interest to the aerospace community as lightweight thermal insulation, the application of conformal polymer coatings to these gels increases their strength significantly, making them potentially useful as structural materials as well. In this work we perform multiscale computer simulations to investigate the tensile and compressive strain behavior of silica and polymer-coated silica aerogels. Aerogels are made up of clusters of interconnected particles of amorphous silica of less than bulk density. We simulate gel nanostructure using a Diffusion Limited Cluster Aggregation (DLCA) procedure, which produces aggregates that exhibit fractal dimensions similar to those observed in real aerogels. We have previously found that model gels obtained via DLCA exhibited stress-strain curves characteristic of the experimentally observed brittle failure. However, the strain energetics near the expected point of failure were not consistent with such failure. This shortcoming may be due to the fact that the DLCA process produces model gels that are lacking in closed-loop substructures, compared with real gels. Our model gels therefore contain an excess of dangling strands, which tend to unravel under tensile strain, producing non-brittle failure. To address this problem, we have incorporated a modification to the DLCA algorithm that specifically produces closed loops in the model gels. We obtain the strain energetics of interparticle connections via atomistic molecular statics, and abstract the collective energy of the atomic bonds into a Morse potential scaled to describe gel particle interactions. Polymer coatings are similarly described. We apply repeated small uniaxial strains to DLCA clusters, and allow relaxation of the center eighty percent of the cluster between strains. The simulations produce energetics and stress-strain curves for looped and nonlooped clusters, for a variety of densities and interaction parameters.