Engineering topics
Nyman, Gunnar
Publications and source records attributed to Nyman, Gunnar.
Reaction rate constant for radiative association of CF{sup +}
Reaction rate constants and cross sections are computed for the radiative association of carbon cations (C{sup +}) and fluorine atoms (F) in their ground states. We consider reactions through the electronic transition 1{sup 1}Π → X{sup 1}Σ{sup +} and rovibrational transitions on the X{sup 1}Σ{sup +} and a{sup 3}Π potentials. Semiclassical and classical methods are used for the direct contribution and Breit–Wigner theory for the resonance contribution. Quantum mechanical perturbation theory is used for comparison. A modified formulation of the classical method applicable to permanent dipoles of unequally charged reactants is implemented. The total rate constant is fitted to the Arrhenius–Kooij formula in five temperature intervals with a relative difference of <3%. The fit parameters will be added to the online database KIDA. For a temperature of 10–250 K, the rate constant is about 10{sup −21} cm{sup 3} s{sup −1}, rising toward 10{sup −16} cm{sup 3} s{sup −1} for a temperature of 30 000 K.
Application of a New Ensemble Conserving Quantum Dynamics Simulation Algorithm to Liquid para-Hydrogen and ortho-Deuterium
Here, we apply the Feynman-Kleinert Quasi-Classical Wigner (FK-QCW) method developed in our previous work [Smith et al., J. Chem. Phys. 142, 244112 (2015)] for the determination of the dynamic structure factor of liquid para-hydrogen and ortho-deuterium at state points of (T = 20.0 K, n = 21.24 nm -3 ) and (T = 23.0 K, n = 24.61 nm -3 ), respectively. When applied to this challenging system, it is shown that this new FK-QCW method consistently reproduces the experimental dynamic structure factor reported by Smith et al. [J. Chem. Phys. 140, 034501 (2014)] for all momentum transfers considered. Moreover, this shows that FK-QCW provides a substantial improvement over the Feynman-Kleinert linearized path-integral method, in which purely classical dynamics are used. Furthermore, for small momentum transfers, it is shown that FK-QCW provides nearly the same results as ring-polymer molecular dynamics (RPMD), thus suggesting that FK-QCW provides a potentially more appealing algorithm than RPMD since it is not formally limited to correlation functions involving linear operators.
A new class of ensemble conserving algorithms for approximate quantum dynamics: Theoretical formulation and model problems
We develop two classes of quasi-classical dynamics that are shown to conserve the initial quantum ensemble when used in combination with the Feynman-Kleinert approximation of the density operator. These dynamics are used to improve the Feynman-Kleinert implementation of the classical Wigner approximation for the evaluation of quantum time correlation functions known as Feynman-Kleinert linearized path-integral. As shown, both classes of dynamics are able to recover the exact classical and high temperature limits of the quantum time correlation function, while a subset is able to recover the exact harmonic limit. A comparison of the approximate quantum time correlation functions obtained from both classes of dynamics is made with the exact results for the challenging model problems of the quartic and double-well potentials. It is found that these dynamics provide a great improvement over the classical Wigner approximation, in which purely classical dynamics are used. In a special case, our first method becomes identical to centroid molecular dynamics.
The rate constant for radiative association of HF: Comparing quantum and classical dynamics
Radiative association for the formation of hydrogen fluoride through the A{sup 1}Π → X{sup 1}Σ{sup +} and X{sup 1}Σ{sup +} → X{sup 1}Σ{sup +} transitions is studied using quantum and classical dynamics. The total thermal rate constant is obtained for temperatures from 10 K to 20 000 K. Agreement between semiclassical and quantum approaches is observed for the A{sup 1}Π → X{sup 1}Σ{sup +} rate constant above 2000 K. The agreement is explained by the fact that the corresponding cross section is free of resonances for this system. At temperatures below 2000 K we improve the agreement by implementing a simplified semiclassical expression for the rate constant, which includes a quantum corrected pair distribution. The rate coefficient for the X{sup 1}Σ{sup +} → X{sup 1}Σ{sup +} transition is calculated using Breit–Wigner theory and a classical formula for the resonance and direct contributions, respectively. In comparison with quantum calculations the classical formula appears to overestimate the direct contribution to the rate constant by about 12% for this transition. Below about 450 K the resonance contribution is larger than the direct, and above that temperature the opposite holds. The biggest contribution from resonances is at the lowest temperature in the study, 10 K, where it is more than four times larger than the direct. Below 1800 K the radiative association rate constant due to X{sup 1}Σ{sup +} → X{sup 1}Σ{sup +} transitions dominates over A{sup 1}Π → X{sup 1}Σ{sup +}, while above that temperature the situation is the opposite.