Overtone spectroscopy and dynamics in monodeuteroacetylene (C2HD)

Complementary experimental, ab initio and dynamical investigations are reported on monodeuteroacetylene, C₂HD ( ¹Σ⁺). All experimental spectroscopic results previously reported in the literature on C₂HD, i.e. from 500 to 16000 cm⁻¹ are gathered. New results are included, which are obtained from the analysis of absorption data recorded with a Fourier transform interferometer at high resolution between 4600 and 9000 cm⁻¹. The presence of numerous weak bands along the whole spectral range is analysed in terms of systematic anharmonic couplings. The entire set of energy data is then used to produce thirty-five vibrational frequencies and anharmonicities from a fit of the vibrational energies to a Dunham-type expansion, and the vibrational level density is extrapolated, up to higher energy. One- and two-dimensional potential energy and dipole moment surfaces refined from new ab initio results are fitted to a selected set among those experimental data, associated to the stretch overtones. The iterative procedure involving an original package of computer programs is described. The evolution of the overtone intensities of the CH and CD stretches, up to η = 4, is interpreted on that basis in terms of electric and mechanical anharmonicity contributions. Eventually, dynamical aspects are studied thanks to the newly introduced vibrograms, which allow to obtain the time recurrences of the vibrational dynamics. Using the Gutzwiller and Berry-Tabor trace formulas, these vibrational recurrences are semiclassically assigned to periodic orbits of the classical Hamiltonian given by the Dunham expansion.     

Graphical approach to correlations in high‐spin systems

We present a graphical technique for generating and indexing spin monomials of high‐spin systems. The procedure consists of developing a graph with at most n line segments from each node in a given row to the one immediately lower, where n is the multiplicity of single‐particle spin function. The paths lead to monomials with definite M values at each node. This technique has been used to generate and diagonalize the model spin Hamiltonian. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 73: 389–393, 1999     

Vibrational spectrum of Li3 first‐excited electronic doublet state: Geometric‐phase effects and statistical analysis

We present J=0 calculations of all bound and pseudobound vibrational states of Li₃ in its first‐excited electronic doublet state by using a realistic double many‐body expansion potential‐energy surface and a minimum‐residual filter diagonalization technique. The action of the system Hamiltonian on the wave function was evaluated by the spectral transform method in hyperspherical coordinates. Calculations of the vibrational spectra were carried out both without consideration and with consideration of geometric‐phase effects. Dynamic Jahn–Teller and geometric‐phase effects are found to play a significant role, while the calculated fundamental symmetric stretching frequency is larger by 8.3% than its reported experimental value of 326 cm⁻¹. From the neighbor‐spacing distributions of the levels, it is observed that the title vibrational spectrum is quasiregular in the short range and quasi‐irregular in the long range. By the Δ₂ standard defined in this article, it is found that the spectra are more nonuniform than those of the “trough” states for the ground electronic state. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 75: 89–109, 1999     

Mapped grid methods for long-range molecules and cold collisions

The paper discusses ways of improving the accuracy of numerical calculations for vibrational levels of diatomic molecules close to the dissociation limit or for ultracold collisions, in the framework of a grid representation. In order to avoid the implementation of very large grids, Kokoouline et al. [J. Chem. Phys. 110, 9865 (1999)] have proposed a mapping procedure through introduction of an adaptive coordinate x subjected to the variation of the local de Broglie wavelength as a function of the internuclear distance R. Some unphysical levels ("ghosts") then appear in the vibrational series computed via a mapped Fourier grid representation. In the present work the choice of the basis set is reexamined, and two alternative expansions are discussed: Sine functions and Hardy functions. It is shown that use of a basis set with fixed nodes at both grid ends is efficient to eliminate "ghost" solutions. It is further shown that the Hamiltonian matrix in the sine basis can be calculated very accurately by using an auxiliary basis of cosine functions, overcoming the problems arising from numerical calculation of the Jacobian J(x) of the R-->x coordinate transformation.     

An SCF-CI method for determining the potential energy surface of a triatomic molecule

A self-consistent-field (SCF)-configuration interaction (CI) (SCF-CI) method for determining the potential energy surface of a triatomic molecule from the observed vibrational band origins has been suggested. By this method, the SCF-CI procedure in the internal coordinates is used to calculate the vibrational bond origins and their first derivatives with respect to parameters in the potential energy function using the exact vibrational Hamiltonian, and the optimizer LMF in the nonlinear-squares problem is employed to optimize parameters in the potential energy function. This approach is used to optimize the potential energy function of the water molecule. The standard deviation of this fitting to the 70 observed band origins is 1.154cm-1.     

Global analysis of composite particles

The theory of vibrations of a composite particle when vibrational amplitudes are not constrained to be small according to the Eckart conditions is developed using the methods of differential topology. A global classical Hamiltonian appropriate for this system is given, and for the case of the molecular vibration–rotation problem, it is transformed into a global quantum Hamiltonian operator. It is shown that the zeroth-order term in the global Hamiltonian operator is identical to the Wilson–Howard Hamiltonian; higher-order terms are shown to give successively better approximations to the large amplitude problem. Generalized Eckart conditions are derived for the global classical Hamiltonian; the quantum equivalent of these conditions along with the quantum equivalent of the Eckart conditions are given. The spectrum of the global Hamiltonian operator is discussed and it is shown that the calculation of the vibration–rotation energy states of the system reduces to the same straight-forward procedure, the solution of a secular determinant, as was carried out for the Wilson–Howard Hamiltonian at a later time by Nielsen.     

Internal dynamics of nonrigid molecules. I. Application to acetone

The group theory for nonrigid molecules is used for studying the internal dynamics of the two equivalent C_3v rotor “bent” molecules. Special emphasis is given to the deduction of the symmetry basis vectors which represent in box form the Hamiltonian operator. It is shown that these basis vectors may be advantageously employed in order to simplify the resolution of the two-rotor equation. The procedure is applied to the acetone molecule. It is found that the lowest solutions are clustered into groups of four. The four lowest levels are related to vibrational states, the upper 64 to vibro–rotational states, in which the rotors are rotating in a restricted manner. Only few states show some cogwheel effect. Internal rotation contributions to the principal thermodynamic parameters of acetone are also computed.     

The quasiparticle concept in vibrational–electronic problems in molecules. I. Partitioning of the vibrational–electronic Hamiltonian

The partitioning of the vibrational–electronic Hamiltonian is presented. This partitioning is based on a new quasiparticle transformation that is constructed in such a way that the adiabatic approximation is included into the unperturbed Hamiltonian; nonadiabacity, anharmonicity, and electron correlation are treated as perturbations. We also present the second quantization treatment for bosons. The many body perturbation theory expansion for the vibrational–electronic Hamiltonian is suggested. A comparison of this approach is made with gradient techniques.     

Discrete variable representation method applied to the determination of rotation-vibration bound states of NO2

The discrete variable representation method is applied to the determination of the rotation-vibration energy levels of the fundamental electronic state of NO₂. The Hamiltonian is expressed in Johnson hyperspherical coordinates and developed on a DVR basis for each internal coordinate, while parity-adapted linear combinations of Wigner functions are used to describe the rotational motion. The diagonalization of the Hamiltonian matrix is performed using the Lanczos algorithm for large symmetric and Hermitian matrices. Results for rovibrational states up to J = 11 for the first five vibrational energy levels are presented. © 1997 John Wiley & Sons, Inc.     

Separation of the vibrational,rotational, and translational motions of polyatomic molecules using curvilinear vibrational coordinates

A new method is suggested for separating the vibrational, rotational, and translational motions of polyatomic molecules using curvilinear vibrational coordinates that are linear with respect to the natural vibrational coordinates. It is shown that, in this case, Coriolis interactions between the vibrational and rotational motions are absent. The solutions of the anharmonic vibrational-rotational problems in the curvilinear and linear vibrational coordinates are compared. The absence of Coriolis interactions between the vibrational and rotational motions in the curvilinear vibrational coordinates is proved numerically. The same conclusion is additionally supported by calculations of the anharmonic vibrational energy levels for the H₂O, H₂S, NO₂, SO₂, and ClO₂ molecules in the linear and curvilinear vibrational coordinates using the Hamiltonian designed in the curvilinear vibrational coordinates with and without Coriolis vibrational-rotational interactions. Volgograd Pedagogical University. Translated fromZhurnal Strukturnoi Khimii, Vol. 36, No. 2, pp. 239–254, March–April, 1995. Translated by I. Izvekova     

From concepts to algorithms for the treatment of large amplitude internal motions and unimolecular reactions

The formalism of the large amplitude path Hamiltonian provides a theoretical framework for the study of dynamical problems ranging from anharmonic vibrations to unimolecular reactions. A hierarchy of models at different degrees of sophistication can be elaborated and some of them have been put into practical use through the development of the DiNa package. The simplest level requires just the characterization of all the stationary points encountered along a large amplitude path and their quadratic environments. Perturbation theory can next be used to analyze the role of anharmonicity in the vibrational modulation of physico-chemical observables for semirigid systems. An extension of the same approach to saddle points allows the computation of reliable reaction rates taking into account curvature and tunneling effects. Finally, full characterization of the harmonic valley surrounding the path allows study of vibrational modulation in flexible systems and energy flow between different degrees of freedom during chemical reactions.     

Probing the spin multiplicity of gas-phase polycyclic aromatic hydrocarbons through their infrared emission spectrum: A theoretical study

The anharmonic infrared emission spectrum following an optical excitation has been calculated for a variety of polycyclic aromatic hydrocarbon molecules in their ground singlet electronic state or in their triplet state. The computational protocol relies on second-order perturbation theory and involves a quartic vibrational Hamiltonian, the vibrational quantum numbers being sampled according to a Monte Carlo procedure. In the case of neutral naphthalene, the IR spectrum obtained in the (ground) singlet state differs significantly from the spectrum in the triplet state, especially for out-of-plane CH bending modes. Although not as prominent, spectral differences in larger molecules are still observable.     

On exact calculation of response properties of oscillators in static electric field: A Fourier grid Hamiltonian approach. I. One-dimensional systems

The Fourier grid Hamiltonian method is used to calculate the response properties of different types of 1-d (one-dimensional) quantum oscillators in a uniform static electric field. The calculations are potentially exact. Excepting the harmonic oscillator, the other model oscillators studied are seen to possess nonlinear polarizabilities. In general, the polarizabilities are not monotonic functions of appropriate vibrational quantum numbers. The exact nature of this vibrational-state dependence of polarizabilities is shown to depend on the type of mechanical anharmonicity in which the nuclei move and the nature of electrical anharmonicity characterizing the field–oscillator coupling. The large vibrational contribution to nonlinear polarizabilities often predicted for real diatomics could therefore originate from the mechanical and electrical anharmonicities of the potential in which the nuclei move when placed in a static electric field. © 1994 John Wiley & Sons, Inc.     

Vibrational Stark effect: Theoretical determination through the semiempirical AM1 method

The vibrational Stark effect of a series of small molecules has been calculated by means of the semiempirical AM1 method through addition of the electron–field interaction term in the one-electron Hamiltonian. Optimized geometrical parameters along with harmonic frequencies and line intensities are determined for different strengths of the applied uniform electric field. The perturbed spectra are compared with theoretical studies carried out at the ab initio level and with experimental results. The vibrational Stark effect of the retinal molecule is also computed, showing that this kind of study is feasible in systems of biochemical interest. © 1992 by John Wiley & Sons, Inc.     

Diagrammatic formulation of the second‐order many‐body multipartitioning perturbation theory

The second‐order multireference perturbation theory employing multiple partitioning of the many‐electron Hamiltonian into a zero‐order part and a perturbation is formulated in terms of many‐body diagrams. The essential difference from the standard diagrammatic technique of Hose and Kaldor concerns the rules of evaluation of energy denominators which take into account the dependence of the Hamiltonian partitioning on the bra and ket determinantal vectors of a given matrix element, as well as the presence of several two‐particle terms in zero‐order operators. The novel formulation naturally gives rise to a “sum‐over‐orbital” procedure of correlation calculations on molecular electronic states, particularly efficient in treating the problems with large number of correlated electrons and extensive one‐electron bases. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 73: 395–401, 1999     

Mechanism and control of molecular energy flow: a modeling perspective

 Vibrational energy flow in organic molecules occurs by a multiple-time-scale mechanism that can be modeled by a single exponential only in its initial stages. The mechanism is a consequence of the hierarchical structure of the vibrational Hamiltonian, which leads to diffusion of vibrational wavepackets on a manifold with far fewer than the 3N−6 dimensions of the full vibrational state space. The dynamics are controlled by a local density of states, which does not keep increasing with molecular size. In addition, the number of vibrational coordinates severely perturbed during chemical reaction is small, leading to preservation of the hierarchical structure at chemically interesting energies. This regularity opens up the possibility of controlling chemical reactions by controlling the vibrational energy flow. Computationally, laser control of intramolecular vibrational energy redistribution can be modeled by quantum-classical, or by purely quantum-mechanical models of the molecule and control field. Received: 26 July 2002 / Accepted: 30 September 2002 / Published online: 2 December 2002 Electronic Supplementary Material to this paper can be obtained by using the Springer Link server located at http://dx.doi.org/10.1007/s00214-002-0394-2. Acknowledgements. This work was supported by NSF grant CHE 9986670. Correspondence to: M. Gruebele e-mail: gruebele@scs.uiuc.edu     

A truncation/recoupling method for eigenvalues and eigenvectors ideal for parallel computation

Summary We review a new truncation/recoupling method to obtain eigenvalues and eigenvectors of anF-mode vibrational Hamiltonian. In particular we focus on the block diagonalization aspects of the method which are ideal for massive parallelization, and we demonstrate this using H₂O₂ as an example. We then present vibrational energies for non-rotating HO₂ and HCN, which illustrate several key advantages of this method.     

Valence : A Massively Parallel Implementation of the Variational Subspace Valence Bond Method

This work describes the software package, Valence , for the calculation of molecular energies using the variational subspace valence bond (VSVB) method. VSVB is an ab initio electronic structure method based on nonorthogonal orbitals. Important features of practical value include high parallel scalability, wave functions that can be constructed automatically by combining orbitals from previous calculations, and ground and excited states that can be modeled with a single configuration or determinant. The open-source software package includes tools to generate wave functions, a database of generic orbitals, example input files, and a library build intended for integration with other packages. We also describe the interface to an external software package, enabling the computation of optimized molecular geometries and vibrational frequencies. © 2019 Wiley Periodicals, Inc.