Hierarchies of approximations

Election nuclear dynamics theory is a direct nonadiababc dynamics approach to molecular processes and uses an electi onic basis of atomic orbitals attached to dynamical centers, whose positions and momenta are dynamical variables. Although computationally intensive, this approach is general and has a systematic hierarchy of approximations when applied in an ab initio fashion. It can also be applied with semiempirical treatment of electronic degrees of freedom [4]. It is important to recognize that the reactants in this approach are not forced to follow a certain reaction path but for a given set of initial conditions the entire system evolves in time in a completely dynamical manner dictated by the inteiparbcle interactions. 

Equations (169) and (171), together with Eqs. (170), fomi the basic equations that enable the calculation of the non-adiabatic coupling matrix. As is noticed, this set of equations creates a hierarchy of approximations starting with the assumption that the cross-products on the right-hand side of Eq. (171) have small values because at any point in configuration space at least one of the multipliers in the product is small [115]. 

The GIAO-MP2/TZP calculated 13C NMR chemical shifts of the cyclopropylidene substituted dienyl cation 27 show for almost all carbon positions larger deviations from the experimental shifts than the other cations 22-26. The GIAO-MP2/TZP method overestimates the influence of cr-delocalization of the positive charge into the cyclopropane subunit on the chemical shifts. Electron correlation corrections for cyclopropylidenemethyl cations such as 27 and 28 are too large to be adequately described by the GIAO-MP2 perturbation theory method and higher hierarchies of approximations such as coupled cluster models are required to rectify the problem. 

Starting from the normal mode approximation, one can introduce anharmonicity in different ways. Anharmonic perturbation theory [206] and local mode models [204] may be useful in some cases, where anharmonic effects are small or mostly diagonal. Vibrational self-consistent-field and configuration-interaction treatments [207, 208] can also be powerful and offer a hierarchy of approximation levels. Even more rigorous multidimensional treatments include variational calculations [209], diffusion quantum Monte Carlo, and time-dependent Hartree approaches [210]. 

In terms of these conditions, a fc-particle hierarchy of approximations can be defined, with Hartree-Fock as the one-particle approximation for closed-shell states. Unfortunately, the stationarity conditions do not determine the fully, and for their constmction additional information is required, which essentially guarantees -representability. Nevertheless, the fe-particle hierarchy based on the irreducible stationarity conditions opens a promising way for the solution of the -electron problem. 

Time-correlation functions can also be computed from their continued fraction representations (see section III.D) by exploiting a hierarchy of approximations of the following kind. Suppose that the nth order random force has a white spectrum.34 It then follows that... 

It is clear that the various density functional schemes for molecular applications rely on physical aiguments pertaining to specific systems, such as an electron gas, and fitting of parameters to produce eneigy functionals, which are certainly not universal. By focusing on the energy functional one has given up the connection to established quantum mechanics, which employs Hamiltonians and Hilbert spaces. One has then also abandoned the tradition of quantum chemistry of the development of hierarchies of approximations, which allows for step-wise systematic improvements of the description of electronic properties. 

This hierarchy of approximations constitutes a guide for judging the quality of the methods discussed below that aim to improve the efficiency of quasiparticle electron propagator calculations. 

The computational effort of a molecular calculation can be reduced significantly, if only a few electrons are taken into account explicitly and the interaction with the rest is approximated by means of effective Hamiltonians. The first step in the hierarchy of approximations is the so-called frozen-core approximation. 

Of paramount importance in this latter category is the Hartree-Fock approximation. The so-called Hartree-Fock limit represents a well-defined plateau, in terms of its methematical and physical properties, in the hierarchy of approximate solutions to Schrodinger s electronic equation. In addition, the Hartree-Fock solution serves as the starting point for many of the presently employed methods whose ultimate goal is to achieve solutions to equation (5) of chemical accuracy. A discussion of the Hartree-Fock method and its associated concept of a self-consistent field thus provides a natural starting point for the discussion of the calculation of potential surfaces. 

A hierarchy of approximations now exists for calculating interactions between a charged particle and a charged, planar interface in electrolyte solutions. At moderate surface potentials less than approximately 2(kT/e the linear Poisson-Boltzmann equation provides a good approximation in many circumstances, provided the solution is a 1 1 electrolyte at low to moderate ionic strength. The relative simplicity of the linear equation makes it particularly useful for examining problems that are complicated in other ways, such as interactions involving many particles, interactions with deformable interfaces, and interactions where the detailed structure and properties of the particle (or macromolecule) play an important role. 

Perdew has described a hierarchy of approximate treatments of the exchange-correlation term as Jacob s Ladder. The ladder is grounded in HF theory here on earth, and reaches to heaven, where the exact functional is found. Along the way are five rungs, each defining a set of assumptions made in creating an exchange-correlation expression. 

It is important to use the exact strain tensor definition, Eq. (6), to achieve rotational invariance with respect to lattice rotation the conventional linear strain tensor only provides differential rotational invariance of u in Eq. (7).hierarchy of approximations may be used for the elastic tensor 7. The most rigorous approach is to transform the bulk elastic tensor c according to... 

Interpreting bulk properties qualitatively on the basis of microscopic properties requires only consideration of the long-range attractive forces and short-range repulsive forces between molecules it is not necessary to take into account the details of molecular shapes. We have already shown one kind of potential that describes these intermolecular forces, the Lennard-Jones 6-12 potential used in Section 9.7 to obtain corrections to the ideal gas law. In Section 10.2, we discuss a variety of intermolecular forces, most of which are derived from electrostatic (Coulomb) interactions, but which are expressed as a hierarchy of approximations to exact electrostatic calculations for these complex systems. 

The application of this strategy to solvation problems means to introduce in the full quantum description of the problem a hierarchy of approximations which, step by step, leads to simpler algorithms. At each step the loss of information, and the errors introduced are controlled. The sequence of approximations we shall summarize here has been accompanied by a set of protocols to check the quality of the results. Our exposition will be rather concise. It is a summary of those given in our recent review on continuum solvation (Tomasi and Persico, 1994). All information about, the quality of the results may be found there, as well as in the source papers. 

Table 2. Hierarchie of approximate quantum-mechanical theories of unsaturated...

Using successively larger and larger operator manifolds we may define a hierarchy of approximate propagator methods based on Eq. (58). Also the choice of reference state influences Eq. (58) through the definition of the superoperator binary products in Eqs (52) and (60). It has been our experience that it is important to maintain a balance between the level of sophistication of the operator manifold and of the reference state. Better results are generally obtained this way, as we will see examples of in the subsequent sections. 

A solution of this equation (Heully et al. 1986) for X in closed form is not known and a hierarchy of approximations for X is used in practice instead (see, for example, the atomic balance in Section 2.4). 

However, this approach is too expensive because, even for very small molecules such as HF in moderately small orbital spaces, Ndet becomes very large—several billions or more. Thus, we must instead be content to work with approximate FCI wave functions. Fortunately, several useful hierarchies of approximations to the FCI solution have been developed, enabling us to approach the FCI wave function closely, even for rather large systems. 

The convergence of the quantum chemical calculations can be studied in terms of two types of hierarchies. First, the quality of the calculations depends on the flexibility of the MO space the AOs that are used to expand the MOs may be extended in a well-defined and systematic manner, thereby establishing a one-electron hierarchy. Second, we can increase the excitation level in coupled-cluster theory or the order of perturbation expansion, thus setting up an n-electron hierarchy of approximate electronic wave functions. In Fig. 5, the roles of the one-electron and the ra-electron hierarchies are illustrated. 

A systematic development of relativistic ECPs is possible within wavefunction-based approaches, which offer a hierarchy of approximations and the possiblity to converge, at least in principle, to the exact solution. Thus it is possible to... 

Clearly, the second of these is of little value unless the first is satisfied within some well-defined and well-understood hierarchy of approximation we do not wish to have explanations of the strength of a chemical bond using theories which are not capable of giving a good quantitative calculation of bond energies, for example. 

However, j-coupling is retained (with the attendant large number of coupled equations). We make a further approximation in which the rotational motion of the diatom is described in the strong interaction region by an uncoupled adiabatic description. The details of this formulation have been given elsewhere and so we proceed with only a brief outline of the approximation. The reduced dimensionality hierarchy of approximations is then presented. 

Methods for solving the electronic equation (1) have evolved into sophisticated codes that incorporate a hierarchy of approximations that can be used as black boxes to achieve accurate descriptions for the PES for ground states of molecular systems. Popular codes include Gaussian [12], GAMESS [13], and Jaguar [14] for finite molecules and VASP [15], CRYSTAL [16], CASTEP [17], and Sequest [18] for periodic systems. 

Historically, the first derivations of approximate relativistic operators of value in molecular science have become known as the Pauli approximation. Still, the best-known operators to capture relativistic corrections originate from those developments which provided well-known operators such as the spin-orbit or the mass-velocity or the Darwin operators. Not all of these operators are variationally stable, and therefore they can only be employed within the framework of perturbation theory. Nowadays, these difficulties have been overcome by, for instance, the Douglas-Kroll-Hess hierarchy of approximate Hamiltonians and the regular approximations to be introduced in a later section, so that operators such as the mass-velocity and Darwin terms are no... 

---