Electronic quantum motions in hydrogen molecule ion

The hydrogen molecule ion is a two‐center force system expressed under the prolate spheroidal coordinates, whose quantum motions and quantum trajectories have never been addressed in the literature before. The momentum operators in this coordinate system are derived for the first time from the Hamilton equations of motion and used to construct the Hamiltonian operator. The resulting Hamiltonian comprises a kinetic energy T and a total potential V_Total consisting of the Coulomb potential and a quantum potential. It is shown that the participation of the quantum potential and the accompanied quantum forces in the force interaction within H₂⁺ is essential to develop an electronic motion consistent with the prediction of the probability density function |Ψ|². The motion of the electron in H₂⁺ can be either described by the Hamilton equations derived from the Hamiltonian H = T_K + V_Total or by the Lagrange equations derived from the Lagrangian H = T_K ? V_Total. Solving the equations of motion with different initial positions, we show that the solutions yield an assembly of electronic quantum trajectories whose distribution and concentration reconstruct the σ and π molecular orbitals in H₂⁺. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Generalized electronic diabatic scheme: Diagonalizing the electronic Hamiltonian for artificial molecular systems. How do molecular meccanos move?

A quantum–classic model is presented and used to describe systems ranging from normal molecules up to electronic systems sensed in real space. The quantum system is a set of n‐electrons; a positive background in real space completes the model. A generalized electronic diabatic (GED) theory is introduced. The diabatic functions diagonalize the electronic Hamiltonian for any arrangement of the positive background. Physical quantum states are represented as linear superpositions in the diabatic basis; this latter is always fixed. For systems sensed in real space, the coefficients of the linear superposition are functions of the real space configuration coordinates. Physical changes are produced by interactions with external sources/sinks of energy. An interaction couples different diabatic states; diagonalizing the electronic Hamiltonian plus the couplings leads to new coefficients describing physical states. Among other things, these couplings can be used to simulate the effects produced by scanning tunneling microscopy, atomic force, and transmission electron microscopy on substrates located in real space. The important thing is that time–evolution in electronic Hilbert space can be related to actual motion in real space. The experiment of lateral hopping of a substrate on a metallic surface induced by vibration excitation and followed with scanning tunneling microscope is discussed. A result of the present work is that motion of molecular meccanos reflects then time–evolution in electronic Hilbert space. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Rovibrational molecular hamiltonian in mixed bond-angle and umbrella-like coordinates

A new exact quantum mechanical rovibrational Hamiltonian operator for molecules exhibiting large amplitude inversion and torsion motions is derived. The derivation is based on a division of a molecule into two parts: a frame and a top. The nuclei of the frame only are used to construct a molecular system of axes. The inversion motion of the frame is described in the umbrella-like coordinates, whereas the torsion motion of the top is described by the nonstandard torsion angle defined in terms of the nuclear vectors and one of the molecular axes. The internal coordinates chosen take into account the properties of the inversion and torsion motions. Vibrational s and rotational Omega vectors obtained for the introduced internal coordinates determine the rovibrational tensor G defined by simple scalar products of these vectors. The Jacobian of the transformation from the Cartesian to the internal coordinates considered and the G tensor specify the rovibrational Hamiltonian. As a result, the Hamiltonian for penta-atomic molecules like NH2OH with one inverter is presented and a complete set of the formulas necessary to write down the Hamiltonian of more complex molecules, like NH2NH2 with two inverters, is reported. The approach considered is essentially general and sufficiently simple, as demonstrated by derivation of a polyatomic molecule Hamiltonian in polyspherical coordinates, obtained by other methods with much greater efforts.     

Expressing kinetic‐energy for vibration–rotation motions in general rotating systems of axes and introducing quasirectilinear vibrational coordinates to simplify Hamiltonian forms

To understand how the internal and rotational motions of a polyatomic system depend on which rotating system of axes is selected, we derived the explicit form of the atomic velocities determined by an observer stationed on the general rotating system of axes. Using the derived velocities, we formulated the kinetic energy expression for vibration–rotation motions with respect to the rotating system of axes. From this expression, we clarified covariant metric tensors under zero angular momentum, which have been confused with an erroneous expression even in the professional literature, and the relationship between the kinetic energy expression and the rotating system of axes. Furthermore, to simplify the Hamiltonian form, we introduced quasirectilinear vibrational coordinates to describe the Hamiltonian. The resulting Hamiltonian form is superior to those of the previous studies in that the kinetic and potential energy expressions are simple and the vibrational frequencies are independent of the original internal coordinates used. In fact, we show that its application for three examples is useful. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 83: 22–29, 2001     

Generalization of classical mechanics for nuclear motions on nonadiabatically coupled potential energy surfaces in chemical reactions

Classical trajectory study of nuclear motion on the Born-Oppenheimer potential energy surfaces is now one of the standard methods of chemical dynamics. In particular, this approach is inevitable in the studies of large molecular systems. However, as soon as more than a single potential energy surface is involved due to nonadiabatic coupling, such a naive application of classical mechanics loses its theoretical foundation. This is a classic and fundamental issue in the foundation of chemistry. To cope with this problem, we propose a generalization of classical mechanics that provides a path even in cases where multiple potential energy surfaces are involved in a single event and the Born-Oppenheimer approximation breaks down. This generalization is made by diagonalization of the matrix representation of nuclear forces in nonadiabatic dynamics, which is derived from a mixed quantum-classical representation of the electron-nucleus entangled Hamiltonian [Takatsuka, K. J. Chem. Phys. 2006, 124, 064111]. A manifestation of quantum fluctuation on a classical subsystem that directly contacts with a quantum subsystem is discussed. We also show that the Hamiltonian thus represented gives a theoretical foundation to examine the validity of the so-called semiclassical Ehrenfest theory (or mean-field theory) for electron quantum wavepacket dynamics, and indeed, it is pointed out that the electronic Hamiltonian to be used in this theory should be slightly modified.     

A constant of the motion for the two-centre Kepler problem

The one-centre Coulombic potential gives rise, in non-relativistic mechanics, to additional constants of the motion which form the components of the Runge–Lenz vector. By a study of this vector, an extra constant of the motion is derived for the corresponding two-centre problem. The result holds quite generally in a space of arbitrary dimension and is applicable to both classical and quantum mechanics; but breaks down when relativistic corrections, which destroy the extra symmetry of the one-centre Coulombic potential, are taken into account. The effect of further Coulombic centres and of varying the form of the potential is briefly discussed. In particular a constant of the motion is derived for a two-centre potential which has both Coulombic and simple harmonic terms. The relationship between these constants of the motion and the separation of the Hamiltonian into spheroidal coordinates is noted (this had previously been known only for the two-centre Coulomb problem in three-dimensional space). Finally the application to the hydrogen molecule ion, treated in the adiabatic approximation, is pointed out. The extra constant of the motion is seen to account for an observed apparent breakdown in the noncrossing rule for the potential energy curves.     

A density functional representation of quantum chemistry. III. Rigorous realization of the program in lattice space

A mathematically rigorous reformulation of molecular quantum mechanics in terms of the particle density operator and a canonically conjugated phase field is given. Using a momentum cutoff, it is shown that the usual molecular Hamiltonian can be expressed in terms of the particle density operator and a rigorously defined phase operator. It is shown that this Hamiltonian converges strongly to the cutoff-free Hamiltonian. In spite of the fact that this Hamiltonian is of second order in the phase operators, all hitherto published expressions are not correct. Unfortunately, the correct formulation destroys the intuitive appeal of using the particle density operator as a coordinate for the many-body problems of quantum chemistry. Unless somebody provides an essential new and clever idea, we propose to resist the fascination of a local quantum field theory of molecular matter in terms of the particle density operator.     

A simple, direct derivation and proof of the validity of the SLLOD equations of motion for generalized homogeneous flows

We present a simple and direct derivation of the SLLOD equations of motion for molecular simulations of general homogeneous flows. We show that these equations of motion (1) generate the correct particle trajectories, (2) conserve the total thermal momentum without requiring the center of mass to be located at the origin, and (3) exactly generate the required energy dissipation. These equations of motion are compared with the g-SLLOD and p-SLLOD equations of motion, which are found to be deficient. Claims that the SLLOD equations of motion are incorrect for elongational flows are critically examined and found to be invalid. It is confirmed that the SLLOD equations are, in general, non-Hamiltonian. We derive a Hamiltonian from which they can be obtained in the special case of a symmetric velocity gradient tensor. In this case, it is possible to perform a canonical transformation that results in the well-known DOLLS tensor Hamiltonian.     

Density functional theory of born couplings: Consequences for electron flow in Jahn–Teller molecules and superconductors

The dynamics of Jahn–Teller systems has recently been discussed in terms of generalized electronic charge and current densities in nuclear-coordinate space. The introduction of the electronic phase as a function of both electronic and nuclear coordinates, in addition to the electronic density, was a crucial component of this formulation. Here, a densitybased treatment of Born couplings is derived from first-principles quantum mechanics beyond the Born–Oppenheimer approximation. Because of the degenerate electronic configuration of a Jahn–Teller molecule, there are an infinite number of ways in which the charge distribution can be oriented for the same energy, leading to a vanishing bond hardness for the molecule in the symmetric nuclear configuration. Further, the moving nuclear framework serves as the perturbation necessary to define the orientation of the charge density, leading to unhindered rotation of the charge cloud. This leads to the dynamical Jahn–Teller problem, namely, the coupling of electronic and nuclear motions through the Born coupling terms. Applications to superconductivity theory are discussed. © 1995 John Wiley & Sons, Inc.     

A Basis for Orbital Symmetry Rules

The influence of molecular symmetry on reaction rates is examined with an approach in which reactions are viewed as electronic transitions between states of reacants and products (described, in turn by quasiadiabatic potential surface). The moleculer Hamiltonian is used to derive selection rules for these transitions. The complete Hamilatonian has no useful symmetery. Neglect of non-Born-Oppenheimer and spin-orbit terms (and of other angular momentum coupling terms) leads to an apporixmate Hamiltonian and to selection rules which from the basis of the Woodward-Hoffmann rules. This apporch provides an alternative to the adiabatic potantial surfaces, reaction coordinates, and transition state theory used in more familiar discussions of the Woodward-Hoffmann rules. Further, it provides a particulary clear method for discussing violations of these symmetry rules, and for differentiating concerted and nonconcerted reactions.     

Electron trajectories in molecular orbitals

The time-dependent Schrödinger equation can be rewritten so that its interpretation is no longer probabilistic. Two well-known and related reformulations are Bohmian mechanics and quantum hydrodynamics. In these formulations, quantum particles follow real, deterministic trajectories influenced by a quantum force. Generally, trajectory methods are not applied to electronic structure calculations as they predict that the electrons in a ground-state, real, molecular wavefunction are motionless. However, a spin-dependent momentum can be recovered from the nonrelativistic limit of the Dirac equation. Therefore, we developed new, spin-dependent equations of motion for the quantum hydrodynamics of electrons in molecular orbitals. The equations are based on a Lagrange multiplier, which constrains each electron to an isosurface of its molecular orbital, as required by the spin-dependent momentum. Both the momentum and the Lagrange multiplier provide a unique perspective on the properties of electrons in molecules.     

A density functional representation of quantum chemistry. I. Motivation and general formalism

The molecular N-electron problem is formulated in terms of the particle density ??( r ) = ∑σ( r ? q _n) and its canonically conjugated field. Weyl's representation of the canonical commutation relations is used and a characteristic functional of the molecular system is introduced as the expectation value of Weyl's functional exp {i(α, Q) + i(β, P)}. The Hamiltonian equation of motion for the characteristic functional is derived and its cumulant approximation is discussed. The first paper is restricted to a general discussion, independent of a special choice of the Hamiltonian. The density functional representation may be considered as a step towards a formulation of a renormalized quantum chemistry. The lowest order cumulants of the characteristic functional are accessible to actual precision measurements; on the other hand, most experiments are correctly represented by a characteristic functional of the Gaussian type. The significance of the higher order correlations is neither theoretically nor experimentally disclosed yet, their possible importance for biological systems is mentioned, and the question is raised of how far quantum mechanics is empirically confirmed for extremely complex systems.     

Prediction of the mechanisms of enzyme-catalysed reactions using hybrid quantum mechanical/molecular mechanical methods

The use of hybrid methods, involving both quantum mechanics and molecular mechanics, to model the mechanism of enzyme-catalysed reactions, is discussed. Two alternative approaches to treating the electrostatic interactions between the quantum mechanical and molecular mechanical regions are studied, involving either the inclusion of this term in the electronic Hamiltonian (QM/MM), or evaluating it purely classically (MO + MM). In the latter scheme, possible problems of using force fields that are standard for macromolecular modelling are identified. The use of QM/MM schemes to investigate the mechanism of the enzymes thymidine phosphorylase (ThdPase) and protein tyrosine phosphatase (PTP) is described. For both systems, transition states have been identified using a PM3 Hamiltonian. For ThdPase, concerted motion of the enzyme during the course of the reaction is suggested and, for PTP, a two-step dephosphorylation reaction is indicated, both with quite low barriers.     

Complete error potential

The complete error potential is a function of electronic and nuclear coordinates which has to be added to the physical potential in order to obtain the eigenfunctions of the molecular Hamiltonian in an a priori specified approximation. This potential is essential for the error theory, in particular it furnishes a simple justification of the Born-Oppenheimer approximation. It is given as a sum of the electronic error potential and that of the intramolecular forces, whose importance is discussed.Dedicated to Professor J. Koutecký on the occasion of his 65th birthday     

Solution of the “Classical” Schrödinger equation for a driven symmetric triple well: A comparison with its classical counterpart

The behavior of a driven symmetric triple well potential has been studied by developing an algorithm where the well‐established Bohmian mechanics and time‐dependent Fourier Grid Hamiltonian method are incorporated and the quantum theory of motion (QTM) phase space structures of the particle are constructed, both in “nonclassical” and “classical” limits. Comparison of QTM phase space structures with their classical analogues shows both similarity as well as dissimilarities. The temporal nature and the spatial symmetry of applied perturbation play crucial roles in having similar phase space structures. © 2016 Wiley Periodicals, Inc.     

A Continuum Reaction Field Theory of Polarizable, Nondipolar, Quadrupolar Solvents: Ab Initio Study of Equilibrium Solvation in Benzene

A continuum theory to describe solvation in nondipolar quadrupolar solvents is developed by accounting for electronic polarizability. A general Hamiltonian for a solute–solvent system in an arbitrary nonequilibrium configuration is obtained in terms of two field variables—densities of the solvent quadrupole and induced dipole moments. Equilibrium solvation is studied by optimizing this Hamiltonian with account of cavity boundaries. As an application, electronic structures and free energies of small molecules in benzene are examined with ab initio methods. Solvation stabilization due to solvent quadrupole moments is found to be substantial; for the solutes considered here, it is comparable to and often in excess of that arising from solvent-induced dipole moments.     

Nonadiabatic electron wavepacket dynamics of molecules in an intense optical field: an ab initio electronic state study

A theory of quantum electron wavepacket dynamics that nonadiabatically couples with classical nuclear motions in intense optical fields is studied. The formalism is intended to track the laser-driven electron wavepackets in terms of the linear combination of configuration-state functions generated with ab initio molecular orbitals. Beginning with the total quantum Hamiltonian for electrons and nuclei in the vector potential of classical electromagnetic field, we reduce the Hamiltonian into a mixed quantum-classical representation by replacing the quantum nuclear momentum operators with the classical counterparts. This framework gives equations of motion for electron wavepackets in an intense laser field through the time dependent variational principle. On the other hand, a generalization of the Newtonian equations provides a matrix form of forces acting on the nuclei for nonadiabatic dynamics. A mean-field approximation to the force matrix reduces this higher order formalism to the semiclassical Ehrenfest theory in intense optical fields. To bring these theories into a practical quantum chemical package for general molecules, we have implemented the relevant ab initio algorithms in it. Some numerical results in the level of the semiclassical Ehrenfest-type theory with explicit use of the nuclear kinematic (derivative) coupling and the velocity form for the optical interaction are presented.