Theoretical methods Hamiltonian approximation

The ASEP/MD method, acronym for Averaged Solvent Electrostatic Potential from Molecular Dynamics, is a theoretical method addressed at the study of solvent effects that is half-way between continuum and quantum mechanics/molecular mechanics (QM/MM) methods. As in continuum or Langevin dipole methods, the solvent perturbation is introduced into the molecular Hamiltonian through a continuous distribution function, i.e. the method uses the mean field approximation (MFA). However, this distribution function is obtained from simulations, i.e., as in QM/MM methods, ASEP/MD combines quantum mechanics (QM) in the description of the solute with molecular dynamics (MD) calculations in the description of the solvent. 

Numerical solutions of the Schrodinger equation can be obtained within several degrees of approximation, for almost any system, using its exact Hamiltonian. Density functional theory has proven to be one of the most effective techniques, because it provides significantly greater accuracy than Hartree-Fock theory with just a modest increase in computational cost.io> 3-45 The accuracy of DFT method is comparable, and even greater than other much more expensive theoretical methods that also include electron correlation such as second and higher order perturbation theory. 

The forces binding the atoms A and B together in AB are chemical in nature and must be introduced, at least approximately, in the Hamiltonian. Then it should be possible to apply the same theoretical methods (e.g., HNC and MS approximations) used to study strong electrolytes to investigate incomplete dissociation in weak electrolytes as well. The binding between A and B is quite distinct from the ion pair formation observed for higher valence electrolytes (Fig. 9). In these cases no alterations in the Hamiltonian models already discussed were required to account qualitatively for the experimental observations. 

The examples included in this chapter use some different theoretical models for the interpretation of, primarily. UPS valence band data, for pristine and doped systems as well as for the initial stages of interface formation for metals on conjugated molecules. Among the theoretical methods used in the examples are sem-iempirical Hartree-Fock methods such as the modified neglect of diatomic overlap (MNDO) [48,49] and Austin Model 1 (AMI) [50], the nonempirical valence effective Hamiltonian (VEH) pseudopotential method [51,52], ab initio Hartree-Fock techniques, and the local spin density (LSD) approximation [53,54],... 

In this review we shall first establish the theoretical foundations of the semi-classical theory that eventually lead to the formulation of the Breit-Pauli Hamiltonian. The latter is an approximation suited to make the connection to phenomenological model Hamiltonians like the Heisenberg Hamiltonian for the description of electronic spin-spin interactions. The complete derivations have been given in detail in Ref. (21), but turn out to be very involved and are thus scattered over many pages in Ref. (21). For this reason, we aim here at a summary that is as brief and concise as possible so that all relevant connections between different levels of approximation are evident. This allows us to connect present-day quantum chemical methods to phenomenological Hamiltonians and hence to establish and review the current status of these first-principles methods applied to transition-metal clusters. 

The parameters of Hamiltonians (1) and (2) are determined in our approach by pure theoretical way using different quantum chemical models and calculations unlike the traditional fitting the experimental thermodynamic and dielectric data. Our method of the many-pseudospin clusters [ 1,4] seems to be the most reliable way of determination. The latter are obtained in this case within the static approximation from the system of equations for a typical crystal fragment (cluster) for all possible proton distributions on H-bonds. The left-hand side of any equation expresses the cluster total energy in terms of Jy, while the right-hand side is determined by means of the quantum chemical calculation of this energy. 

The MFA [1] introduces the perturbation due to the solvent effect in an averaged way. Specifically, the quantity that is introduced into the solute molecular Hamiltonian is the averaged value of the potential generated by the solvent in the volume occupied by the solute. In the past, this approximation has mainly been used with very simplified descriptions of the solvent, such as those provided by the dielectric continuum [2] or Langevin dipole models [3], A more detailed description of the solvent has been used by Ten-no et al. [4], who describe the solvent through atom-atom radial distribution functions obtained via an extended version of the interaction site method. Less attention has been paid, however, to the use of the MFA in conjunction with simulation calculations of liquids, although its theoretical bases are well known [5]. In this respect, we would refer to the papers of Sese and co-workers [6], where the solvent radial distribution functions obtained from MD [7] calculations and its perturbation are introduced a posteriori into the molecular Hamiltonian. 

The approximation techniques described in the earlier sections apply to any (non-relativistic) quantum system and can be universally used. On the other hand, the specific methods necessary for modeling molecular PES that refer explicitly to electronic wave function (or other possible tools mentioned above adjusted to describe electronic structure) are united under the name of quantum chemistry (QC).15 Quantum chemistry is different from other branches of theoretical physics in that it deals with systems of intermediate numbers of fermions - electrons, which preclude on the one hand the use of the infinite number limit - the number of electrons in a system is a sensitive parameter. This brings one to the position where it is necessary to consider wave functions dependent on spatial r and spin s variables of all N electrons entering the system. In other words, the wave functions sought by either version of the variational method or meant in the frame of either perturbational technique - the eigenfunctions of the electronic Hamiltonian in eq. (1.27) are the functions D(xi,..., xN) where. r, stands for the pair of the spatial radius vector of i-th electron and its spin projection s to a fixed axis. These latter, along with the... 

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... 

Ab initio calculations of electronic wave functions are well established as useful and powerful theoretical tools to investigate physical and chemical processes at the molecular level. Many computational packages are available to perform such calculations, and a variety of mathematical methods exist to approximate the solutions of the electronic hamiltonian. Each method is based (or should be) on a well defined physical model, specified by a certain partition of the electronic hamiltonian, in such a way as to include a subset of all the interactions present in the exact one. It is expected that this subset contains the most important effects to describe consistently the situation of interest. The identification of which physical interactions to include is a major step in developing and applying quantum chemical theory to the study of real problems. 

From the conceptual point of view, there are two general approaches to the molecular structure problem the molecular orbital (MO) and the valence bond (VB) theories. Technical difficulties in the computational implementation of the VB approach have favoured the development and the popularization of MO theory in opposition to VB. In a recent review [3], some related issues are raised and clarified. However, there still persist some conceptual pitfalls and misinterpretations in specialized literature of MO and VB theories. In this paper, we attempt to contribute to a more profound understanding of the VB and MO methods and concepts. We briefly present the physico-chemical basis of MO and VB approaches and their intimate relationship. The VB concept of resonance is reformulated in a physically meaningful way and its point group symmetry foundations are laid. Finally it is shown that the Generalized Multistructural (GMS) wave function encompasses all variational wave functions, VB or MO based, in the same framework, providing an unified view for the theoretical quantum molecular structure problem. Throughout this paper, unless otherwise stated, we utilize the non-relativistic (spin independent) hamiltonian under the Bom-Oppenheimer adiabatic approximation. We will see that even when some of these restrictions are removed, the GMS wave function is still applicable. 

Section III is devoted to illustrating the first theoretical tool under discussion in this review, the GME derived from the Liouville equation, classical or quantum, through the contraction over the irrelevant degrees of freedom. In Section III.A we illustrate Zwanzig s projection method. Then, in Section III.B, we show how to use this method to derive a GME from Anderson s tight binding Hamiltonian The second-order approximation yields the Pauli master equation. This proves that the adoption of GME derived from a Hamiltonian picture requires, in principle, an infinite-order treatment. The case of a vanishing diffusion coefficient must be considered as a case of anomalous diffusion, and the second-order treatment is compatible only with the condition of ordinary... 

The DIM method is most commonly employed as a semi-empirical technique. The fragment Hamiltonian matrices are usually related to atomic and diatomic energies by making various approximations for the overlap matrices. Both the form of the DIM equation and the chosen set of PBF must be sufficient to account for all the qualitative features of the system being studied. Under such circumstances the approach may offer acceptable accuracy for modest computational effort. Given the input of experimental and accurate theoretical data for the fragments, it is not unreasonable to suppose that the method can yield results comparable to those from larger... 

---