This chapter presents a brief account on both approaches, analytic and synthetic, within chemistry. It reviews the different levels at which molecular systems are traditionally studied. The chapter explores the Born–Oppenheimer approximation and then revisiting the methodologies available to deal with both electronic configuration and dynamics, distinguishing in this latter case between molecular dynamics for a few degree-of-freedom systems and statistical mechanics. It focuses on some fundamental aspects of Bohmian mechanics and describes some specific illustrations of the application of Bohmian mechanics to problems typical of the levels mentioned. Molecular systems typically include a rather large number of degrees of freedom, which makes unaffordable their direct study. In the wave model each atom is assumed to be a point-like, positively charged nucleus surrounded by an electronic cloud, which is associated with a certain orbital and represents the electron probability distribution in configuration space.
This chapter presents a brief account on both approaches, analytic and synthetic, within chemistry. It reviews the different levels at which molecular systems are traditionally studied. The chapter explores the Born–Oppenheimer approximation and then revisiting the methodologies available to deal with both electronic configuration and dynamics, distinguishing in this latter case between molecular dynamics for a few degree-of-freedom systems and statistical mechanics. It focuses on some fundamental aspects of Bohmian mechanics and describes some specific illustrations of the application of Bohmian mechanics to problems typical of the levels mentioned. Molecular systems typically include a rather large number of degrees of freedom, which makes unaffordable their direct study. In the wave model each atom is assumed to be a point-like, positively charged nucleus surrounded by an electronic cloud, which is associated with a certain orbital and represents the electron probability distribution in configuration space.
Perhaps because of the popularity that trajectory-based methodologies have always had in Chemistry and the important role they have played, Bohmian mechanics has been increasingly accepted within this community, particularly in those areas of the theoretical chemistry based on quantum mechanics, e.g., quantum chemistry, chemical physics, or physical chemistry. From a historical perspective, this evolution is remarkably interesting, particularly when the scarce applications of Madelung's former hydrodynamical formulation, dating back to the late 1960s and the 1970s, are compared with the many different applications available at present. As also happens with classical methodologies, Bohmian trajectories are essentially used to described and analyze the evolution of chemical systems, to design and implement new computational propagation techniques, or a combination of both. In the first case, Bohmian trajectories have the advantage that they avoid invoking typical quantum-classical correspondence to interpret the corresponding phenomenon or process, while in the second case quantum-mechanical effects appear by themselves, without the necessity to include artificially quantization conditions. Rather than providing an exhaustive revision and analysis of all these applications (excellent monographs on the issue are available in the literature for the interested reader, which can be consulted in the bibliography here supplied), this Chapter has been prepared in a way that it may serve the reader to acquire a general view (or impression) on how Bohmian mechanics has permeated the different traditional levels or pathways to approach molecular systems in Chemistry: electronic structure, molecular dynamics and statistical mechanics. This is done with the aid of some illustrative examples -- theoretical developments in some cases and numerical simulations in other cases.
In this chapter, a glimpse of the way how the concept of quantum trajectories has been used in different areas of quantum chemistry and chemical physics is provided. Putting aside ontological issues, Bohmian mechanics has been exploited within these fields with practical purposes, both analytical and computational. Analytically, although quantum trajectories provide the same information at a predictive level as standard quantum mechanics, they have the advantage that one can follow and describe the evolution of quantum systems by means of streamlines, whose results are less vague and more intuitive than standard wave packet simulations. On the other hand, computationally, Bohmian mechanics has been the source of a number of numerical algorithms devised as an alternative to other methods based on the wave function or the probability density. This will be seen throughout the analysis and discussion of different applications, which go from electronic structure to chemical dynamics, thus covering essentially the research spectrum of quantum chemistry and chemical physics. Furthermore, in the form of a final summary, we present a unifying scheme based on minimum principles, which allows us to connect electronic structurewith classical and quantum (Bohmian) chemical dynamics.
Chemical dynamics, in principle, should be understood by solving the time-dependent Schrödinger equation for a molecular system, describing motion of the nuclei and electrons. However, the computational efforts to solve this partial second-order differential equation scales exponentially with the system size, which prevents us from getting exact numerical solutions for systems larger than 4-5 atoms. Thus, approximations simplifying the picture are necessary. The so-called Born-Oppenheimer approximation, separating motion of the electrons and nuclei is the central one: solution to the electronic Schrödinger equation defines the potential energy surface on which the nuclear motion unfolds, and there are standard quantum chemistry software packages for solving the electronic Schrödinger equation. For the nuclear Schrödinger equation, however, there are no widely applicable quantum-mechanical approaches, and most simulations are performed using classical Newtonian mechanics which is often adequate due to large nuclear masses. However, the nuclear quantum effects are significant for chemical processes involving light nuclei at low energies, and including these effects into simulation, even approximately, is highly desirable. In this dissertation, an approximate methodology of including quantum-mechanical effects within the quantum trajectory or the de Broglie-Bohm formulation of the Schrödinger equations is developed. Use of the trajectory framework makes the approach scalable to hundreds of degrees of freedom. The methodology is applied to study high-dimensional systems (solid He4 and others) relevant to chemistry.
Since the early days of quantum mechanics, the wave function theory has proven to be very successful in describing many different quantum processes and phenomena. However, in many problems of quantum chemistry and solid-state physics, where the dimensionality of the systems studied is relatively high, ab initio calculations of the structure of atoms, molecules, clusters, and crystals, and their interactions are very often prohibitive. Hence, alternative formulations based on the direct use of the probability density, gathered under what is generally known as the density matrix theory [1], were also developed since the very beginning of the new mechanics. The independent electron approximation or Thomas-Fermi model, and the Hartree and Hartree-Fock approaches are former statistical models developed in that direction [2]. These models can be considered direct predecessors of the more recent density functional theory (DFT) [3], whose principles were established by Hohenberg, Kohn, and Sham [4,5] in the mid-1960s. According to this theory, the fundamental physical information about a many-body system is provided by single-particle densities in a three-dimensional space, which are obtained variationally within a time-independent framework. When compared with other previous formalisms, DFT presents two clear advantages: (1) it is able to treat many-body problems in a sufficiently accurate way and (2) it is computationally simple. This explains why it is one of the most widely used theories to deal with electronic structure-the electronic ground-state energy as a function of the position of the atomic nuclei determines molecular structures and solids, providing at the same time the forces acting on the atomic nuclei when they are not at their equilibrium positions. At present, DFT is used routinely to solve many problems in gas phase and condensed matter. Furthermore, it has made possible the development of accurate molecular dynamics schemes in which the forces are evaluated quantum mechanically ‘‘on the fly.’’ Nonetheless, DFT is a fundamental tool provided the systems studied are relatively large; for small systems, standard methods based on the use of the wave function render quite accurate results [6]. Moreover, it is also worth stressing that all practical applications of DFT rely on essentially uncontrolled approximations [7] (e.g., the local density approximation [4,5], the local spin-density approximation, or generalized gradient approximations [8]), and therefore the validity of DFT is conditioned to its capability to provide fairly good values of the experimental data.
Bohm theory is a formulation of Quantum Mechanics that characterises the state of a quantum system according to both the wave function, as in the conventional formulation, and the coordinates (positions) of all the particles that evolve in time drawing quantum continuous trajectories. Furthermore, a statistical ensemble of all the possible trajectories, raising from the impossibility to know the initial position of all the particles, establishes the exact correspondence with the traditional Quantum Mechanics. From a computational point of view, Bohm theory has found many applications in Chemical Physics especially to develop new methodologies for solving the Schrödinger equation and to address semi-classical approximations of Quantum Mechanics. From a theoretical point of view, the most appealing feature of Bohm theory is its capability to supply a conceptual map between the quantum formalism and our representation of what a chemical system is. Chemical systems are composed of molecules, but the same idea of molecule requires a specific arrangement in the space of particles, i.e., the nuclei of the atoms. The statistical description of conventional Quantum Mechanics on the basis of wave function alone is insufficient to establish a clear correspondence with such a picture of molecules. Indeed, chemists employ usually Classical Mechanics in order to overcome this drawback of the standard quantum theory. On the other hand, if the particles position is included in the quantum formalism, as Bohm theory does, the map can be defined in a self-consistent way. In other words, Bohm theory appears to be the suitable quantum framework to represent molecules and their motion. The chemical representation of molecular systems finds a natural correspondence with a single Bohm trajectory, since it is always implicitly assumed that molecular components have specific spatial position independently of our knowledge about it. Consequently, we develop a quantum method whose fundamental assumption is that a single Bohm trajectory, i.e., a quantum molecular trajectory, describes the molecular systems and the molecular motion correctly. First of all, we examine the correspondence between a single Bohm trajectory and the conventional Quantum Mechanics, without using the ensemble of trajectories. We verify that such a correspondence exists through numerical simulations and we prove formally that the statistical properties of a single Bohm trajectory explain the probabilistic description of Quantum Mechanics. Once the consistency of this original approach has been established, we investigate the predicted properties. For instance, we take into account the constants of motion (such as the energy) corresponding to the time evolution of the coordinates and the behaviour of simple chemical systems, e.g., the vibrational motion of single molecules interacting with a resonant field. In this way, unexpected features of the molecular motion are found. Secondly, we tackle the challenge of describing many components systems (like the chemical systems in ordinary conditions). As a matter of fact, the computation of the Bohm trajectory and of the wave function is extremely demanding. However, the statistical properties of the Bohm trajectory allow the derivation of stochastic theories for examining the dynamics of open quantum systems, i.e., few molecules (or few degrees of freedom) interacting with their environment (the other molecules). One of the developed stochastic methods correlates the dynamics of the reduced density matrix, for the degrees of freedom of interest, to the evolution of the corresponding Bohm coordinates. In other words, the Bohm equation, determining the set of all the particles velocities according to the full wave function, is replaced with a stochastic one that approximates the velocity of a subset of coordinates according to the reduced density matrix. In such a way, the quantum fluctuations induced by the environment are taken into account. The advantage of this method concerns its capability of describing quantum systems, including open quantum systems, in terms of a quantum trajectory. This could allow the understanding of the molecular motion during a spectroscopical experiment. The possibility of investigating reactive systems, such as conformational changes, is particularly interesting. As a matter of fact, chemical reactions can be completely characterised only through the particles motion and we define the suit- able quantum methodology providing a self-consistent description of the molecular motion.