Quantum chemistry methods approximations

When the temperature ranges of k and k- data are not approximately coincident, it may be necessary to correct the activation parameters to the same mean temperature by using heat capacity data. This correction can be estimated by statistical mechanics, after finding (e.g., by quantum chemistry methods) a structure for the activated complex. 

The theoretical tools of quantum chemistry briefly described in the previous chapter are numerously implemented, sometimes explicitly and sometimes implicitly, in ab initio, density functional (DFT), and semi-empirical theories of quantum chemistry and in the computer program suits based upon them. It is usually believed that the difference between the methods stems from different approximations used for the one- and two-electron matrix elements of the molecular Hamiltonian eq. (1.177) employed throughout the calculation. However, this type of classification is not particularly suitable in the context of hybrid methods where attention must be drawn to the way of separating the entire molecular system (eventually - the universe itself) into parts, of which some are treated explicitly on a quantum mechanical/chemical level, while others are considered classically and the rest is not addressed at all. That general formulation allows us to cover both the traditional quantum chemistry methods based on the wave functions and the DFT-based methods, which generally claim... 

A major goal of fundamental research aiming to rationalize the interplay of structure, dynamics, and chemical reactivity, is to determine multidimensional potentials for nuclei in various environments. On the one hand, potential surfaces can be calculated with quantum chemistry methods at various levels of approximation. On the other hand, from the experimentalist viewpoint, vibrational spectroscopy techniques can probe dynamics of atoms, molecules and ions, in various states of the matter. However, there are fundamental and technical limitations to the determination of potential hypersurfaces from vibrational spectra of complex systems, and the confrontation of experiments with theory is far from being free of ambiguities. Consequently, the interpretation of vibrational spectra remains largely based on experiments. Recent progress in neutron scattering techniques have revealed new dynamics, specially for... 

Relativistic and electron correlation effects play an important role in the electronic structure of molecules containing heavy elements (main group elements, transition metals, lanthanide and actinide complexes). It is therefore mandatory to account for them in quantum mechanical methods used in theoretical chemistry, when investigating for instance the properties of heavy atoms and molecules in their excited electronic states. In this chapter we introduce the present state-of-the-art ab initio spin-orbit configuration interaction methods for relativistic electronic structure calculations. These include the various types of relativistic effective core potentials in the scalar relativistic approximation, and several methods to treat electron correlation effects and spin-orbit coupling. We discuss a selection of recent applications on the spectroscopy of gas-phase molecules and on embedded molecules in a crystal enviromnent to outline the degree of maturity of quantum chemistry methods. This also illustrates the necessity for a strong interplay between theory and experiment. 

We review the Douglas-Kroll-Hess (DKH) approach to relativistic density functional calculations for molecular systems, also in comparison with other two-component approaches and four-component relativistic quantum chemistry methods. The scalar relativistic variant of the DKH method of solving the Dirac-Kohn-Sham problem is an efficient procedure for treating compounds of heavy elements including such complex systems as transition metal clusters, adsorption complexes, and solvated actinide compounds. This method allows routine ad-electron density functional calculations on heavy-element compounds and provides a reliable alternative to the popular approximate strategy based on relativistic effective core potentials. We discuss recent method development aimed at an efficient treatment of spin-orbit interaction in the DKH approach as well as calculations of g tensors. Comparison with results of four-component methods for small molecules reveals that, for many application problems, a two-component treatment of spin-orbit interaction can be competitive with these more precise procedures. 

For most chemical applications, one is not interested in negative energy solutions of a four-component Dirac-type Hamiltonian. In addition, the computational expense of treating four-component complex-valued wave functions often limited such calculations to benchmark studies of atoms and small molecules. Therefore, much effort was put into developing and implementing approximate quantum chemistry methods which explicitly treat only the electron degrees of freedom, namely two- and one-component relativistic formulations [2]. This analysis also holds for a relativistic DFT approach and the solutions of the corresponding DKS equation. 

In summary, the two problem areas of state-of-the-art mobility calculations are the neglect of inelasticity of molecular collisions, especially with respect to rotation, and poor quality or absence of force fields for ion-molecule interactions. However, the impossibility of rigorously solving the Schrodinger equation for polyatomic molecules has stimulated rather than precluded continuous improvement and application of approximate quantum chemistry methods. 

For molecules, Hartree-Fock approximation is the central starting point for most ab initio quantum chemistry methods. It was then shown by Fock that a Slater determinant, a determinant of one-particle orbitals first used by Heisenberg and Dirac in 1926, has the same antisymmetric property as the exact solution and hence is a suitable ansatz for applying the variational principle. 

In the limit of a complete basis set, this equation becomes equivalent to the Schrodinger equation. For a finite basis set, Equation (2) represents the best wave function (in the sense of the variation principle) that can be obtained. It is called the Full Cl (FCI) wave function. It serves as a calibration point for all approximate wave-function methods. It is obvious that many of the coefficients in Equation (3) are very small. We can consider most approximate MO models in quantum chemistry as approximations in one way or the other, where one attempts to include the most important of the configurations in Equation (2). We notice that the FCI wave function and... 

Semi-empirical quantum chemistry methods are based on the Hartree-Fock formalism, but make many approximations and obtain some parameters from empirical data. They are very important in computational chemistry for treating large molecules where the full Hartree-Fock method without the approximations is too expensive. The use of empirical parameters appears to allow some inclusion of electron correlation effects into the methods. Within the framework of Hartree-Fock calculations, some pieces of information (such as two-electron integrals) are sometimes approximated or completely omitted. 

Later, he was instrumental in the development and use of semiempirical quantum chemistry methods, like MINDO, MNDO, and AMI, for analysis of organic reactions [24]. Semiempirical methods are generally based on the Hartree-Fock formalism, but the computational effort is reduced by various approximations to the two-electron and overlap integrals that appear in Hartree-Fock... 

The rigorous Hartree-Fock method without approximations is too expensive to treat large systems such as large organic molecules. Thus semiempirical quantum chemistry methods, which are based on approximated Hartree-Fock formalism by inclusion of some parameters from empirical data, have been introduced to study systems that do not necessarily require the exact quantum solutions to understand the physicochemical properties and are, therefore, very important in simulating large molecular systems. 

Abstract Methods of computational chemistry seem to often be simply a melange of undecipherable acronyms. Frequently, the ability to characterize methods with respect to their quality and applied approximations or to ascribe the proper methodology to the physicochemical property of interest is sufficient to perform research. However, it is worth knowing the fundamental ideas underlying the computational techniques so that one may exploit the approximations intentionally and efficiently. This chapter is an introduction to quantum chemistry methods based on the wave function search in one-electron approximation. 

The conceptual framework for the - semiclassical simulation of ultrafast spectroscopic observables is provided by the Wigner representation of quantum mechanics [2, 3]. Specifically, for the ultrafast pump-probe spectroscopy using classical trajectories, methods based on the semiclassical limit of the Liouville-von Neumann equation for the time evolution of the vibronic density matrix have been developed [4-8]. Our approach [4,6-8] is related to the Liouville space theory of nonlinear spectroscopy developed by Mukamel et al. [9]. It is characterized by the ability to approximately describe quantum phenomena such as optical transitions by averaging over the ensemble of classical trajectories. Moreover, quantum corrections for the nuclear dynamics can be introduced in a systematic manner, e.g. in the framework of the entangled trajectory method [10,11]. Alternatively, these effects can be also accounted for in the framework of the multiple spawning method [12]. In general, trajectory-based methods require drastically less computational effort than full quantum mechanical calculations and provide physical insight in ultrafast processes. Additionally, they can be combined directly with quantum chemistry methods for the electronic structure calculations. 

As proposed by Munn and Hurst [74, 75], an elegant alternative is to employ the generalization of the electrostatic dipole interaction scheme proposed by SUber-stein for atoms [76, 77] to molecules, which consists in evaluating first the local field via the Lorentz-factor tensor [78] and then the macroscopic linear and NLO susceptibilities from the molecular responses calculated using quantum chemistry methods. To alleviate the limitations of the point dipole approximation, the molecule and its molecular responses are usually partitioned in submolecules [79]. Within this approach, the linear optical susceptibility tensor for a crystal with Z molecules labeled k (or Z submolecules labeled kj) per unit cell with a volume V reads ... 

Quantum chemical methods, exemplified by CASSCF and other MCSCF methods, have now evolved to an extent where it is possible to routinely treat accurately the excited electronic states of molecules containing a number of atoms. Mixed nuclear dynamics, such as swarm of trajectory based surface hopping or Ehrenfest dynamics, or the Gaussian wavepacket based multiple spawning method, use an approximate representation of the nuclear wavepacket based on classical trajectories. They are thus able to use the infoiination from quantum chemistry calculations required for the propagation of the nuclei in the form of forces. These methods seem able to reproduce, at least qualitatively, the dynamics of non-adiabatic systems. Test calculations have now been run using duect dynamics, and these show that even a small number of trajectories is able to produce useful mechanistic infomiation about the photochemistry of a system. In some cases it is even possible to extract some quantitative information. 

---