Search NASA⌕ Search

Engineering topics

Babikov, Dmitri

Publications and source records attributed to Babikov, Dmitri.

11 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mixed quantum/classical theory for rotational energy exchange in symmetric-top-rotor + linear-rotor collisions and a case study of the ND 3 + D 2 system

The extension of mixed quantum/classical theory (MQCT) to describe collisional energy transfer is developed for a symmetric-top-rotor + linear-rotor system and is applied to ND 3 + D 2 . State-to-state transition cross sections are computed in a broad energy range for all possible processes: when both ND 3 and D 2 molecules are excited or both are quenched, when one is excited while the other is quenched and vice versa, when the ND 3 state changes its parity while D 2 is excited or quenched, and when ND 3 is excited or quenched while D 2 remains in the same state, ground or excited. In all these processes the results of MQCT are found to approximately satisfy the principle of microscopic reversibility. For a set of sixteen state-to-state transitions available from the literature for a collision energy of 800 cm –1 the values of cross sections predicted by MQCT are within 8% of accurate full-quantum results. A useful time-dependent insight is obtained by monitoring the evolution of state populations along MQCT trajectories. It is shown that, if before the collision, D 2 is in its ground state, the excitation of ND 3 rotational states proceeds through a two-step mechanism in which the kinetic energy of molecule–molecule collision is first used to excite D 2 and only then is transferred to the excited rotational states of ND 3 . Finally, it is found that both potential coupling and Coriolis coupling play important roles in ND 3 + D 2 collisions.

inelastic cross-sections↗
DOE OSTI2023

Description of quantum interference using mixed quantum/classical theory of inelastic scattering

Manifestation of the quantum interference effect in the oscillation of scattering cross section is explored using the N 2 + O system as a case study. Calculations are carried out for two electronic PESs of the system, for various initial rotational states of N 2 , in a broad range of N 2 + O collision energies and using three theoretical methods: two versions of the approximate mixed quantum/classical theory (MQCT and AT-MQCT) and the accurate full-quantum coupled-channel method (implemented in MOLSCAT). A good agreement between different methods is observed, especially at high energies. Elastic scattering cross-sections oscillate as a function of collision energy, which is the result of quantum interference. The effects of initial rotational excitation and of the PES properties are studied in detail. For the final (thermally averaged) cross sections, both MOLSCAT and MQCT calculations predict a rather regular pattern of quantum oscillations that persist through a broad range of collision energies and expand into the low-energy regime where quantum scattering resonances are common. The difference between cross sections predicted by MQCT and MOLSCAT decreases from ~8% at low energies to ~2% at high energies. Finally, experimental data available at high collision energies are well reproduced.

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

Rate coefficients for rotational state-to-state transitions in H 2 O + H 2 O collisions for cometary and planetary applications, as predicted by mixed quantum-classical theory

We present new calculations of collision cross sections for state-to-state transitions between the rotational states in an H 2 O + H 2 O system, which are used to generate a new database of collisional rate coefficients for cometary and planetary applications. Calculations were carried out using a mixed quantum-classical theory approach that is implemented in the code MQCT. The large basis set of rotational states used in these calculations permits us to predict thermally averaged cross sections for 441 transitions in para- and ortho-H 2 O in a broad range of temperatures. It is found that all state-to-state transitions in the H 2 O + H 2 O system split into two well-defined groups, one with higher cross-section values and lower energy transfer, which corresponds to the dipole-dipole driven processes. The other group has smaller cross sections and higher energy transfer, driven by higher-order interaction terms. We present a detailed analysis of the theoretical error bars, and we symmetrized the state-to-state transition matrixes to ensure that excitation and quenching processes for each transition satisfy the principle of microscopic reversibility. We also compare our results with other data available from the literature for H 2 O + H 2 O collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Molecular dynamics on quantum annealers

Abstract In this work we demonstrate a practical prospect of using quantum annealers for simulation of molecular dynamics. A methodology developed for this goal, dubbed Quantum Differential Equations (QDE), is applied to propagate classical trajectories for the vibration of the hydrogen molecule in several regimes: nearly harmonic, highly anharmonic, and dissociative motion. The results obtained using the D-Wave 2000Q quantum annealer are all consistent and quickly converge to the analytical reference solution. Several alternative strategies for such calculations are explored and it was found that the most accurate results and the best efficiency are obtained by combining the quantum annealer with classical post-processing (greedy algorithm). Importantly, the QDE framework developed here is entirely general and can be applied to solve any system of first-order ordinary nonlinear differential equations using a quantum annealer.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

SpectrumSDT: A program for parallel calculation of coupled rotational-vibrational energies and lifetimes of bound states and scattering resonances in triatomic systems

In this work, we present SpectrumSDT – a program for calculations of energies and lifetimes of bound rotational-vibrational states below and scattering resonances above the dissociation threshold on a global potential energy surface of a triatomic system, which may include stable molecules, weekly-bound van-der-Waals complexes, and unbound atom + diatom scattering systems. Large-amplitude vibrational motion is treated explicitly using hyper-spherical coordinates. Three options for the rotational-vibrational interaction are supported: uncoupled (symmetric top rotor), partially coupled (to include interaction between several nearest states only) and full-coupled (vibrating asymmetric-top rotor). In addition to energies and lifetimes, SpectrumSDT is able to integrate ro-vibrational wave functions over the user-defined regions of potential energy surface, which helps to classify these states. In this release of the code, SpectrumSDT is limited to ABA-type molecules with wave functions that do not extend into the regions near Eckart singularities.

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

Four Isotope-Labeled Recombination Pathways of Ozone Formation

A theoretical approach is developed for the description of all possible recombination pathways in the ozone forming reaction, without neglecting any process a priori, and without decoupling the individual pathways one from another. These pathways become physically distinct when a rare isotope of oxygen is introduced, such as 18 O, which represents a sensitive probe of the ozone forming reaction. Each isotopologue of O 3 contains two types of physically distinct entrance channels and two types of physically distinct product wells, creating four recombination pathways. Calculations are done for singly and doubly substituted isotopologues of ozone, eight rate coefficients total. Two pathways for the formation of asymmetric ozone isotopomer exhibit rather different rate coefficients, indicating large isotope effect driven by ΔZPE-difference. Rate coefficient for the formation of symmetric isotopomer of ozone (third pathway) is found to be in between of those two, while the rate of insertion pathway is smaller by two orders of magnitude. These trends are in good agreement with experiments, for both singly and doubly substituted ozone. The total formation rates for asymmetric isotopomers are found to be somewhat larger than those for symmetric isotopomers, but not as much as in the experiment. Altogether, the distribution of lifetimes is found to be very similar for the metastable states in symmetric and asymmetric ozone isotopomers.

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

Solving complex eigenvalue problems on a quantum annealer with applications to quantum scattering resonances

Quantum computing is a new and rapidly evolving paradigm for solving chemistry problems. In previous work, we developed the Quantum Annealer Eigensolver (QAE) and applied it to the calculation of the vibrational spectrum of a molecule on the D-Wave quantum annealer. However, the original QAE methodology was applicable to real symmetric matrices only. For many physics and chemistry problems, the diagonalization of complex matrices is required. For example, the calculation of quantum scattering resonances can be formulated as a complex eigenvalue problem where the real part of the eigenvalue is the resonance energy and the imaginary part is proportional to the resonance width. In the present work, we generalize the QAE to treat complex matrices: first complex Hermitian matrices and then complex symmetric matrices. These generalizations are then used to compute a quantum scattering resonance state in a 1D model potential for O + O collisions. These calculations are performed using both a software (classical) annealer and hardware annealer (the D-Wave 2000Q). The results of the complex QAE are also benchmarked against a standard linear algebra library (LAPACK). Furthermore, this work presents the first numerical solution of a complex eigenvalue problem of any kind on a quantum annealer, and it is the first treatment of a quantum scattering resonance on any quantum device.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗