Car-Parrinello approach

Remler D K and Madden P A 1990 Molecular dynamics without effective potentials via the Car-Parrinello approach Mol. Phys. 70 921-66... 

Regarding mechanical properties of polymers, the efficiency of the Car-Parrinello approach has enabled us to evaluate the ultimate Young s modulus of orthorhombic polyethylene, and demonstrate basis set convergence for that property. 

The Car-Parrinello approach combines an electronic structure method with a classical molecular dynamics scheme and thus unifies two major fields of computational chemistry, which have hitherto been essentially orthogonal. Through this unification a... 

The major drawback for employing the Car-Parrinello approach in dynamics simulations is that since a variational wavefunction is required, the electronic energy should in principle be minimized before the forces on the atoms are calculated. This greatly increases the amount of computer time required at each step of the simulation. Furthermore, the energies calculated with the electronic structure methods currently used in this approach are not exceptionally accurate. For example, it is well established that potential energy barriers, which are of importance to chemical reactivity, often require sophisticated methods to be accurately determined. Nonetheless, the Tirst-principles calculation of the forces during the dynamics is an appealing idea, and will continue to be developed as computer resources expand. 

Rentier DK, PA Madden (1990) Molecular-Dynamics without Effective Potentials Via the Car-Parrinello Approach. Mol. Phys. 70 (6) 921—966... 

Cucinotta et studied the equilibrium between the keto and enol tautomers of acetone, cf. Fig. 5, both in vacuum and in an aqueous solution. They used parameter-free, electronic-structure calculations for a periodically repeated supercell (see section IIB) with the Car-Parrinello approach (see, e.g., ref. 1). In static calculations for solely the acetone molecule they found a total-energy difference between the enol and keto forms of 11.8 kcal/mol. The barrier between the two forms was found to be around 58 kcal/mol (Fig. 5). 

Computational advances which have led to pseudopotential methods becoming efficient are Fast Fourier Transfoms and the approach of Car and Parrinello in which wavefunction coefficients are treated as d5mamical variables. Parts of the calculations on periodic systems are easier to calculate in real space and others in reciprocal space. Fourier transforms are needed to convert between the two spaces. The introduction of Fast Fourier Transforms substantially reduced the computational costs of this conversion. The key feature of the Car-Parrinello approach is the use of minimisation rather than diagonalisation of the Kohn-Sham equations to reach the ground state. 

We also expect that more studies will use MD and the Car-Parrinello approach, while the share of studies based solely on energy minimization and cluster approximations will diminish. It seems reasonable to expect that some of the difficulties encountered with embedding will be solved in the next couple of years. If this becomes true, one might expect that embedding, in fierce competition with the Car-Parrinello approach, will be used extensively in studies on zeolite catalysis. For embedding to become a success, it will be necessary to develop force fields that are compatible with the quantum chemistry model that is used in the description of the catalytic site. Such a compatibility seems possible only if these force fields allow for polarization. 

There are of course methods which attempt to tackle these problems in their complexity, and avoid some of the approximations we have described we quote, as an outstanding example, the Car-Parrinello approach (Car and Parrinello, 1985), that in recent extensions has aimed to give a coherent QM description of such models (Laasonen et al., 1993 Fois et al., 1994). However, these approaches are still in their infancy, and it is advisable to look at other models for liquids and chemical reactions in solutions. 

D. Sebastiani and U. Rothlisberger (2003) Advances in Density Functional Based Modelling Techniques Recent Extensions of the Car-Parrinello Approach. In Medicinal Quantum, Chemistry, eds. P. Carloni and F. Alber Wiley-VCH, Weinheim. Methods and Principles in Medicinal Chemistry, pp. 5-40... 

An elegant way of carrying out AIMD is the Car-Parrinello approach [20,21], in which the dynamics of the nuclei as well as the temporal evolution of the electronic wave function are described by Newtonian equations of motion. Using massively parallel computers, this approach allows the direct simulation of up to a few... 

A preliminary ab initio molecular dynamics simulation has been carried using the Car-Parrinello approach [9]. The volume of the simulation cell was kept constant at the experimental value of the moganite type PON. The cell contained 72 atoms ( P24O24N24 ) The weighted density of states of PON obtained from the ab intio molecular dynamics simulation is compared with the infrared spectrum in Fig. 4. It can be seen that the gross features of the experimental spectrum are reasonably reproduced in the simulation. 

This algorithm alternates between the electronic structure problem and the nuclear motion. It turns out that to generate an accurate nuclear trajectory using this decoupled algorithm the electrons must be fully relaxed to the ground state at each iteration, in contrast to the Car-Parrinello approach, where some error is tolerated. This need for very accurate basis set coefficients means that the minimum in the space of the coefficients must be located very accurately, which can be computationally very expensive. However, conjugate gradients minimisation is found to be an effective way to find this minimum, especially if information from previous steps is incorporated [Payne et al. 1992]. This reduces the number of minimisation steps required to locate accurately the best set of basis set coefficients. 

---