Search NASA⌕ Search

SEARCH · Search NASA

Results for “HARMONIC OSCILLATOR MODELS”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Improved Force-Field Parameters for QM/MM Simulations of the Energies of Adsorption for Molecules in Zeolites and a Free Rotor Correction to the Rigid Rotor Harmonic Oscillator Model for Adsorption Enthalpies

Quantum mechanics/molecular mechanics (QM/MM) simulations provide an efficient avenue for studying reactions catalyzed in zeolite systems; however, the accuracy of such calculations is highly dependent on the zeolite MM parameters used. Previously reported parameters (P1), which were chosen to minimize the root mean square (RMS) deviations of adsorption energies compared with full QM ωB97X-D/6-31+G** adsorption energies, are shown to overestimate binding energies compared with experimental values, particularly for larger substrates. To address this issue, a new parameter set (P2) is derived by rescaling the previously reported characteristic energies of the Lennard-Jones potential in P1. The accuracy of the thermal correction for adsorption enthalpies determined by the rigid rotor-harmonic oscillator approximation (RRHO) is examined and shown to be improved by treating low-lying vibrational modes as free translational and rotational modes via a quasi-RRHO model. Here, with P2 and quasi-RRHO, adsorption energies calculated with QM/MM agree with experimental values with an RMS error of 1.8 kcal/mol for both nonpolar and polar molecules adsorbed in MFI, H-MFI, and H-BEA. By contrast, the RMS error for the same test sets obtained using parameter set P1 is 8.3 kcal/mol. Glucose-fructose isomerization catalyzed by Sn-BEA is taken as an example to demonstrate that improved values for apparent activation energies can be obtained using the methodology reported here. With parameter set P2, the apparent activation energy calculated with QM/MM reproduces the experimental value to within 1 kcal/mol. By contrast, using parameter set P1, the error is -12.9 kcal/mol. (Graph Presented).

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

Generalized factorization and isospectral potentials

It is shown that the method of using excited-state wave functions in constructing exactly solvable potentials as suggested by Berger and Ussembayev [Phys. Rev. A 82, 022121 (2010)] can be generalized further. Using this method, exactly solvable potentials related to the harmonic oscillator and some other potentials are obtained. The nonlinear symmetry algebras underlying the harmonic-oscillator models are also obtained explicitly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Broadband terahertz time-domain spectroscopy of ferroelectric LiTaO{sub 3}: Phonon-polariton dispersion

The lowest frequency IR active A{sub 1}(z) and E(x) modes of a ferroelectric congruent lithium tantalate crystal were studied by the broadband polarized Terahertz Time-Domain Spectroscopy using the coherent THz radiation from 0.2 to 6.5 THz. The dispersion relations of the real and imaginary parts of a polariton wave vector were determined from the complex extraordinary and ordinary refractive indices for the A{sub 1}(z) and E(x) symmetry phonon-polaritons, respectively. The observed complex dispersion relations of two symmetries are in agreement within the experimental uncertainty with the calculated dispersion curves by the damped harmonic oscillator model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Independent bases on the spatial wavefunction of four-identical-particle systems

We construct the independent bases on the spatial wavefunction of four-identical-particle systems classified under the rotational group SO(3) and the permutation group S{sub 4} with the usage of transformation coefficients that relate wavefunctions described in one set of internal coordinates with those in another. The basis functions for N⩽ 2 are presented in the explicit expressions based on the harmonic oscillator model. Such independent bases are supposed to play a key role in the construction of the wavefunctions of the five-quark states and the variation calculation of four-body systems. Our prescription avoids the spurious states and can be programmed for arbitrary N.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Using the electron localization function to correct for confinement physics in semi-local density functional theory

We have previously proposed that further improved functionals for density functional theory can be constructed based on the Armiento-Mattsson subsystem functional scheme if, in addition to the uniform electron gas and surface models used in the Armiento-Mattsson 2005 functional, a model for the strongly confined electron gas is also added. However, of central importance for this scheme is an index that identifies regions in space where the correction provided by the confined electron gas should be applied. The electron localization function (ELF) is a well-known indicator of strongly localized electrons. We use a model of a confined electron gas based on the harmonic oscillator to show that regions with high ELF directly coincide with regions where common exchange energy functionals have large errors. This suggests that the harmonic oscillator model together with an index based on the ELF provides the crucial ingredients for future improved semi-local functionals. For a practical illustration of how the proposed scheme is intended to work for a physical system we discuss monoclinic cupric oxide, CuO. A thorough discussion of this system leads us to promote the cell geometry of CuO as a useful benchmark for future semi-local functionals. Very high ELF values are found in a shell around the O ions, and take its maximum value along the Cu–O directions. An estimate of the exchange functional error from the effect of electron confinement in these regions suggests a magnitude and sign that could account for the error in cell geometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Non-resonant dynamic stark control of vibrational motion with optimized laser pulses

The term dynamic Stark control (DSC) has been used to describe methods of quantum control related to the dynamic Stark effect, i.e., a time-dependent distortion of energy levels. Here, we employ analytical models that present clear and concise interpretations of the principles behind DSC. Within a linearly forced harmonic oscillator model of vibrational excitation, we show how the vibrational amplitude is related to the pulse envelope, and independent of the carrier frequency of the laser pulse, in the DSC regime. Furthermore, we shed light on the DSC regarding the construction of optimal pulse envelopes — from a time-domain as well as a frequency-domain perspective. Finally, in a numerical study beyond the linearly forced harmonic oscillator model, we show that a pulse envelope can be constructed such that a vibrational excitation into a specific excited vibrational eigenstate is accomplished. The pulse envelope is constructed such that high intensities are avoided in order to eliminate the process of ionization.

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

Including screening in van der Waals corrected density functional theory calculations: The case of atoms and small molecules physisorbed on graphene

The Density Functional Theory (DFT)/van der Waals-Quantum Harmonic Oscillator-Wannier function (vdW-QHO-WF) method, recently developed to include the vdW interactions in approximated DFT by combining the quantum harmonic oscillator model with the maximally localized Wannier function technique, is applied to the cases of atoms and small molecules (X=Ar, CO, H{sub 2}, H{sub 2}O) weakly interacting with benzene and with the ideal planar graphene surface. Comparison is also presented with the results obtained by other DFT vdW-corrected schemes, including PBE+D, vdW-DF, vdW-DF2, rVV10, and by the simpler Local Density Approximation (LDA) and semilocal generalized gradient approximation approaches. While for the X-benzene systems all the considered vdW-corrected schemes perform reasonably well, it turns out that an accurate description of the X-graphene interaction requires a proper treatment of many-body contributions and of short-range screening effects, as demonstrated by adopting an improved version of the DFT/vdW-QHO-WF method. We also comment on the widespread attitude of relying on LDA to get a rough description of weakly interacting systems.

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

Vibronic coupling simulations for linear and nonlinear optical processes: Theory

A comprehensive vibronic coupling model based on the time-dependent wavepacket approach is derived to simulate linear optical processes, such as one-photon absorbance and resonance Raman scattering, and nonlinear optical processes, such as two-photon absorbance and resonance hyper-Raman scattering. This approach is particularly well suited for combination with first-principles calculations. Expressions for the Franck-Condon terms, and non-Condon effects via the Herzberg-Teller coupling approach in the independent-mode displaced harmonic oscillator model are presented. The significance of each contribution to the different spectral types is discussed briefly.

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

iSPECTRON: a simulation interface for linear and nonlinear spectra with ab-initio quantum chemistry software

We introduce iSPECTRON, an open source (under the Educational Community License version 2.0) program that parses data from common quantum chemistry software (NWChem, OpenMolcas, Gaussian, Cobramm, etc.), produces the input files for the simulation of linear and nonlinear spectroscopy of molecules with the Spectron code, and analyzes the spectra with a broad range of tools. Vibronic spectra are expressed in term of the electronic eigenstates, obtained through quantum chemistry computations, and vibrational/bath effects are incorporated in the framework of the displaced harmonic oscillator model, where all required quantities are computed at the Franck-Condon point. The code capabilities are illustrated by simulating linear absorption, transient absorption and two dimensional electronic spectra of the pyrene molecule. Two levels of electronic structure theory, TDDFT (with NWChem) and RASSCF/RASPT2 (with OpenMolcas), are compared where possible. Acknowledgements: F.S., A.N., D.R.N., N.G., S.M, M.G. acknowledge support from the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division under Award Nos. DE-SC0019484, KC-030103172684. The Spectron code was developed with support from the National Science Foundation (Grant CHE- 1953045). This research benefited from computational resources provided by EMSL, a DOE Office of Science User Facility sponsored by the Office of Biological and Environmental Research and located at PNNL. PNNL is operated by Battelle Memorial Institute for the United States Department of Energy under DOE Contract No. DE-AC05-76RL1830.

Segatta, Francesco↗
DOE OSTI2022

Configuration Space Integration for Adsorbate Partition Functions: The Effect of Anharmonicity on the Thermophysical Properties of CO–Pt(111) and CH 3 OH–Cu(111)

A method for computing anharmonic thermophysical properties for adsorbates on metal surfaces has been extended to include libration, or frustrated rotation. Classical phase space integration is used with Monte Carlo sampling of the configuration space to obtain the partition function of CO on Pt(111) and CH 3 OH on Cu(111). A minima-preserving neural network potential energy surrogate is used within the integration routines. Direct state counting using discrete variable representation is used to benchmark the results. We find that the phase space integration approach is in excellent agreement with the direct state counting results. Comparison with standard models such as the harmonic oscillator indicates that anharmonicity contributes significantly to the thermodynamic properties of CH 3 OH on Cu(111). We find that there is also a considerable difference between the harmonic oscillator and phase space integration for CO on Pt(111), although the discrepancy can largely be attributed to the presence of multiple binding sites within the unit cell. We demonstrate that a multisite harmonic oscillator model might be sufficient for CO-Pt(111). A more thorough description of the potential energy surface, which can be achieved with phase space integration, is necessary for weakly bound adsorbates such as CH 3 OH. In conclusion, the thermophysical properties were used to calculate free energies of adsorption on the respective metals, and subsequently the equilibrium constants and Langmuir isotherms in relevant temperature ranges. The results show that the choice of model to obtain partition functions greatly affects the resulting surface coverages in kinetic models.

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

Communication: Rigidity of the molecular ion H{sub 5}{sup +}

The fourth-age quantum chemical code GENIUSH is used for the variational determination of rotational-vibrational energy levels corresponding to reduced- and full-dimensional models of H{sub 5}{sup +}, a molecular ion exhibiting several strongly coupled large-amplitude motions. The computations are supplemented with one- and two-dimensional analytic results which help to understand the peculiar rovibrational energy-level structure computed correctly for the first time. An unusual aspect of the results is that the canonical Eckart-embedding of molecule-fixed axes, a cornerstone of the computational spectroscopy of semirigid molecules, is found to be inadequate. Furthermore, it is shown that while the 1D “active torsion” model provides proper results when compared to the full, 9D treatment, models excluding the torsion have limited physical significance. The structure of the rovibrational energy levels of H{sub 5}{sup +} proves that this is a prototypical astructural molecule: the rotational and vibrational level spacings are of the same order of magnitude and the level structure drastically deviates from that computed via perturbed rigid-rotor and harmonic-oscillator models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

High fidelity modeling of thermal relaxation and dissociation of oxygen

A master equation study of vibrational relaxation and dissociation of oxygen is conducted using state-specific O{sub 2}–O transition rates, generated by extensive trajectory simulations. Both O{sub 2}–O and O{sub 2}–O{sub 2} collisions are concurrently simulated in the evolving nonequilibrium gas system under constant heat bath conditions. The forced harmonic oscillator model is incorporated to simulate the state-to-state relaxation of oxygen in O{sub 2}–O{sub 2} collisions. The system of master equations is solved to simulate heating and cooling flows. The present study demonstrates the importance of atom-diatom collisions due to the extremely efficient energy randomization in the intermediate O{sub 3} complex. It is shown that the presence of atomic oxygen has a significant impact on vibrational relaxation time at temperatures observed in hypersonic flow. The population of highly-excited O{sub 2} vibrational states is affected by the amount of atomic oxygen when modeling the relaxation under constant heat bath conditions. A model of coupled state-to-state vibrational relaxation and dissociation of oxygen is also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Observation of discrete oscillations in a model-independent plot of cosmological scale factor versus lookback time and scalar field model

We have observed damped longitudinal cosmological-scale oscillations in a unique model-independent plot of scale factor against lookback time for Type Ia supernovae data. We found several first-derivative relative maxima/minima spanning the range of reported transition redshifts. These extrema comprise two full cycles with a period of approximately 0.15 Hubble times (H{sub 0}=68 km s{sup −1} Mpc{sup −1}). This period corresponds to a fundamental frequency of approximately seven cycles over the Hubble time. Transition-z values quoted in the literature generally fall near these minima and may explain the reported wide spread up to the predicted ΛCDM value of approximately z = 0.77. We also observe second and third harmonics of the fundamental. The scale factor data is analyzed several different ways, including smoothing, Fourier transform, and autocorrelation. We propose a cosmological scalar field harmonic oscillator model for the observation. On this timescale, for a quantum scalar field, the scalar field mass is extraordinarily small at 3 × 10{sup −32} eV. Our scalar field density parameter precisely replaces the ΛCDM dark matter density parameter in the Friedmann equations, resulting in essentially identical data fits, and its present value matches the Planck value. Thus the wave is fundamentally a dark matter wave. We therefore posit that this scalar field manifests itself as the dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

What are the intensities and line-shapes of the twenty four polarization terms in coherent anti-Stokes Raman spectroscopy?

Coherent anti-Stokes Raman spectroscopy (CARS) is conventionally described by just one diagram/term where the three electric field interactions act on the ket side in a Feynman dual time-line diagram in a specific time order of pump, Stokes and probe pulses. In theory, however, any third-order nonlinear spectroscopy with three different electric fields interacting with a molecule can be described by forty eight diagrams/terms. They reduce to just 24 diagrams/terms if we treat the time ordering of the electric field interactions on the ket independently of those on the bra, i.e. the ket and bra wave packets evolve independently. The twenty four polarization terms can be calculated in the multidimensional, separable harmonic oscillator model to obtain the intensities and line-shapes. It is shown that in fs/ps CARS, for the two cases of off-resonance CARS in toluene and resonance CARS in rhodamine 6G, where we use a fs pump pulse, a fs Stokes pulse and a ps probe pulse, we obtain sharp vibrational lines in four of the polarization terms where the pump and Stokes pulses can create a vibrational coherence on the ground electronic state, while the spectral line-shapes of the other twenty terms are broad and featureless. The conventional CARS term with sharp vibrational lines is the dominant term, with intensity at least one order of magnitude larger than the other terms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Molecular reorganization of selected quinoline derivatives in the ground and excited states—Investigations via static DFT

Quinoline derivatives are interesting objects to study internal reorganizations due to the observed excited-state-induced intramolecular proton transfer (ESIPT). Here, we report on computations for selected 12 quinoline derivatives possessing three kinds of intramolecular hydrogen bonds. Density functional theory was employed for the current investigations. The metric and electronic structure simulations were performed for the ground state and first excited singlet and triplet states. The computed potential energy profiles do not show a spontaneous proton transfer in the ground state, whereas excited states exhibit this phenomenon. Atoms in Molecules (AIM) theory was applied to study the nature of hydrogen bonding, whereas Harmonic Oscillator Model of aromaticity index (HOMA) provided data of aromaticity evolution as a derivative of the bridge proton position. The AIM-based topological analysis confirmed the presence of the intramolecular hydrogen bonding. In addition, using the theory, we were able to provide a quantitative illustration of bonding transformation: from covalent to the hydrogen. On the basis of HOMA analysis, we showed that the aromaticity of both rings is dependent on the location of the bridge proton. Further, the computed results were compared with experimental data available. Finally, ESIPT occurrence was compared for the three investigated kinds of hydrogen bridges, and competition between two bridges in one molecule was studied.

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

Franck–Condon Models for Simulating the Band Shape of Electronic Absorption Spectra

Band shape is an essential ingredient in the simulation of electronic absorption spectra. The excitation of multiple series of vibrational levels during an electronic excitation is a main contributor to band shapes. Here we present two simple models based on the Franck–Condon displaced-harmonic-oscillator model. The models are both derived from the time-dependent formulation of electronic spectroscopy. They assume that the transition dipoles do not depend on geometry and that the potential energy surfaces are locally quadratic; one model is second order in time and is called LQ2, and the other is third order in time and is called LQ3. These models are suitable for simulating the unresolved vibronic band shapes of electronic spectra that involve many vibrational modes. The models are straightforward and can be easily applied to simulate absorption spectra that are composed of many electronic transitions. As compared to carrying out molecular dynamics simulations, they require relatively few electronic structure calculations, and the additional cost for constructing the spectra is negligible. Therefore, the models are suitable for simulating the spectra of complex systems such as transition-metal complexes.

UV/Vis spectroscopy↗
DOE OSTI2016

Efficacy of the SU(3) scheme for ab initio large-scale calculations beyond the lightest nuclei

We report on the computational characteristics of ab initio nuclear structure calculations in a symmetry-adapted no-core shell model (SA-NCSM) framework. We examine the computational complexity of the current implementation of the SA-NCSM approach, dubbed LSU3shell, by analyzing ab initio results for 6 Li and 12 C in large harmonic oscillator model spaces and SU3-selected subspaces. We demonstrate LSU3shell’s strong-scaling properties achieved with highly-parallel methods for computing the many-body matrix elements. Results compare favorably with complete model space calculations and significant memory savings are achieved in physically important applications. In particular, a well-chosen symmetry-adapted basis affords memory savings in calculations of states with a fixed total angular momentum in large model spaces while exactly preserving translational invariance.

74 ATOMIC AND MOLECULAR PHYSICS↗