Model core potential techniques

Model core potential (MCP) methods replace core orbitals by a potential just as in ECP. On the other hand, MCP valence orbitals preserve the nodal structure of valence orbitals, unlike ECP valence orbitals. The expectation values of (r ) for the valence orbitals show that the results of MCP are closer to those calculated with all-electron orbitals when comparing MCP, ECP, and the all electron case. Comparisons between MCP and an all electron basis utilizing the full Breit-Pauli spin-orbit Hamiltonian based on multiconfigura-tional quasidegenerate perturbation theory (MCQDPT) calculations show good agreement between the two methods for hydrides of P, As, and Sb. The MCP based spin-orbit calculation appears to be a promising technique, but systematic studies of many different molecular systems are still needed to assess its characteristics and accuracy. 

The final consequence of such a strategy would be to try to eliminate the degrees of freedom of the core electrons as well and to introduce a possibly nonlocal effective potential (pseudopotential), the parameters of which are adjusted either to experiments, which are relativistic from the very beginning, or suitable atomic properties derived from relativistic calculations. This method has developed to the real working horse of relativistic quantum chemistry, and several variants are known as relativistic pseudopotentials, effective core potentials (ECPs) or ab initio model potentials. See Relativistic Effective Core Potential Techniques for Molecules Containing Very Heavy Atoms. 

The drive toward reliable quantum mechanical predictions for large molecular systems is well represented in ECC articles by George Bacskay Solvation Modeling), Krishnan Balasub-ramanian Relativistic Effective Core Potential Techniques for Molecules Containing Very Heavy Atoms), Margareta Blomberg Configuration Interaction PCI-X and Applications), and Keiji Morokuma Hybrid Methods). 

All of the measurements employed the technique described above that involves the analysis of the isotope composition of 02 released from the carrier complexes in preequilibrated solutions. In addition, an established DFT method (mPWPW91)34 with the atomic orbital basis functions, Co, Fe, and Cl (the compact relativistic effective core potential basis CEP-31G),35 N and O (6-311G ), P (6-311G ), C(6-31G), and H (STO-3G),36 were used to calculate the 180 EIE in terms of actual and model structures. The latter approach has also been employed for hypothetical intermediates in enzymes as described below. 

Before any computational study on molecular properties can be carried out, a molecular model needs to be established. It can be based on an appropriate crystal structure or derived using any technique that can produce a valid model for a given compound, whether or not it has been prepared. Molecular mechanics is one such technique and, primarily for reasons of computational simplicity and efficiency, it is one of the most widely used technique. Quantum-mechanical modeling is far more computationally intensive and until recently has been used only rarely for metal complexes. However, the development of effective-core potentials (ECP) and density-functional-theory methods (DFT) has made the use of quantum mechanics a practical alternative. This is particularly so when the electronic structures of a small number of compounds or isomers are required or when transition states or excited states, which are not usually available in molecular mechanics, are to be investigated. However, molecular mechanics is still orders of magnitude faster than ab-initio quantum mechanics and therefore, when large numbers of... 

Two levels of theory are commonly used in the design of the nickel-based catalysts shown in Figure 11 Density Functional Theory (B3LYP functional used with effective core potentials for Ni and 6-3IG for everything else in the complex) and molecular mechanics (both the UFF (4) and reaction force field, RFF (85,86) are used) (87). All these methods are complementary, and the experiments are guided from the results of several calculations using different molecular modeling techniques. 

A further reduction of the computational effort in investigations of electronic structure can be achieved by the restriction of the actual quantum chemical calculations to the valence electron system and the implicit inclusion of the influence of the chemically inert atomic cores by means of suitable parametrized effective (core) potentials (ECPs) and, if necessary, effective core polarization potentials (CPPs). Initiated by the pioneering work of Hellmann and Gombas around 1935, the ECP approach developed into two successful branches, i.e. the model potential (MP) and the pseudopotential (PP) techniques. Whereas the former method attempts to maintain the correct radial nodal structure of the atomic valence orbitals, the latter is formally based on the so-called pseudo-orbital transformation and uses valence orbitals with a simplified radial nodal structure, i.e. pseudovalence orbitals. Besides the computational savings due to the elimination of the core electrons, the main interest in standard ECP techniques results from the fact that they offer an efficient and accurate, albeit approximate, way of including implicitly, i.e. via parametrization of the ECPs, the major relativistic effects in formally nonrelativistic valence-only calculations. A number of reviews on ECPs has been published and the reader is referred to them for details (Bala-subramanian 1998 Bardsley 1974 Chelikowsky and Cohen 1992 Christiansen et... 

To account for the multielectron interaction in alkali systems (alkali-like ions and alkali metal atoms) we successfully used model potential techniques. The basic idea of model potentials V(r), with r been the radial coordinate, is to simulate the multielectron core interaction with the single valence electron by an analytic modification of the Coulomb potential, with the following properties... 

The results of these calculations for this rather simple system indicate that our proposed technique produces the sought for effects, i.e., an effective interaction of simple form that yields a reasonable numerical result in a tractable model space. The need to treat co as a variable parameter, in order to obtain the same result as an exact calculation for He using the same potential, indicates the known deficiencies in our calculations, such as the use of the Reid-soft core potential, instead of one of the more attractive meson-exchange potentials, and the lack of a self-consistent basis space. 

Pseudopotentials describe the interaction of a valence electron with the core of the atoms. They are known in the literature under various names, such as model potentials, effective core potentials,. Model potentials are generally parametrized from atomic spectroscopic data whereas effective core potentials and pseudopotentials are most often derived from ab initio calculations. There is a huge literature on the subject and several review articles. " The recent paper by Krauss and Stevens is recommended for an overall survey of the subject with applications and comparisons with all-electron calculations. The recent review paper of Pelissier et al is devoted to transition elements. In the following we shall only review the main characteristics of the determination of atomic pseudopotentials by the ab initio simulation techniques of Section II.B. 

To date, the only applications of these methods to the solution/metal interface have been reported by Price and Halley, who presented a simplified treatment of the water/metal interface. Briefly, their model involves the calculation of the metal s valence electrons wave function, assuming that the water molecules electronic density and the metal core electrons are fixed. The calculation is based on a one-electron effective potential, which is determined from the electronic density in the metal and the atomic distribution of the liquid. After solving the Schrddinger equation for the wave function and the electronic density for one configuration of the liquid atoms, the force on each atom is ciculated and the new positions are determined using standard molecular dynamics techniques. For more details about the specific implementation of these general ideas, the reader is referred to the original article. ... 

The question of methanol protonation was revisited by Shah et al. (237, 238), who used first-principles calculations to study the adsorption of methanol in chabazite and sodalite. The computational demands of this technique are such that only the most symmetrical zeolite lattices are accessible at present, but this limitation is sure to change in the future. Pseudopotentials were used to model the core electrons, verified by reproduction of the lattice parameter of a-quartz and the gas-phase geometry of methanol. In chabazite, methanol was found to be adsorbed in the 8-ring channel of the structure. The optimized structure corresponds to the ion-paired complex, previously designated as a saddle point on the basis of cluster calculations. No stable minimum was found corresponding to the neutral complex. Shah et al. (237) concluded that any barrier to protonation is more than compensated for by the electrostatic potential within the 8-ring. 

In the present work we study binary soft-sphere mixtures with a core-size ratio mass ratio m2/m = 2.0 and an equimolar concentration (ij = 0.5) for both 3-d and 2-d systems. Using a constant-temperature MD technique and the peiodic boundary conditions, we have carried out MD simulations for the models. The pair potential, Eq. (1), was cut off over the distance rlnumber density wm kept constant, i.e. n" = 0.8 the temperature was varied to achieve a desired Teff. The microscopic time scale was chosen to... 

---