Peter E. Blöchl,Clemens J. Först,Johannes Schimpl
Abstract:A brief introduction to the projector augmented wave method is given and recent developments are reviewed. The projector augmented wave method is an all-electron method for efficient ab-initio molecular dynamics simulations with the full wave functions. It extends and combines the traditions of existing augmented wave methods and the pseudopotential approach. Without sacrificing efficiency, the PAW method avoids transferability problems of the pseudopotential approach and it has been valuable to predict properties that depend on the full wave functions.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop an efficient and accurate method for ab - initio molecular dynamics (AIMD) simulations. Specifically, the paper introduces a new method called Projector Augmented Wave (PAW), which combines the advantages of the augmented - wave method and the pseudopotential method and aims to solve the following problems:
1. **Balance between all - electron accuracy and computational efficiency**:
- **All - electron accuracy**: The PAW method can retain the information of the all - electron wave function, which is crucial for predicting properties that depend on the full wave function. For example, the properties of surfaces and interfaces, magnetic properties, etc.
- **Computational efficiency**: By using a smooth auxiliary wave function, the PAW method maintains computational efficiency while avoiding the transferability problems in traditional pseudopotential methods.
2. **Numerical stability and convergence**:
- **Numerical stability**: The PAW method makes numerical calculations more stable, especially when dealing with the high - momentum part, by introducing a smooth auxiliary wave function.
- **Convergence**: The convergence of the PAW method can be systematically controlled, thus ensuring the accuracy of the calculation results.
3. **Combining the advantages of multiple methods**:
- **Augmented - wave method**: The PAW method inherits the high accuracy of the augmented - wave method in the atomic region.
- **Pseudopotential method**: The PAW method borrows the efficiency of the pseudopotential method in the bonding region and at the same time avoids the problem of information loss in the core region.
4. **Ability of dynamic simulation**:
- **Car - Parrinello method**: The PAW method is combined with the Car - Parrinello molecular dynamics method, making ab - initio molecular dynamics simulations possible and enabling the study of dynamic phenomena and finite - temperature effects.
### Specific implementation
To achieve the above goals, the PAW method adopts the following key techniques and steps:
1. **Wave - function transformation**:
- Through a transformation operator \( T \), the real wave function \( |\Psi_n\rangle \) is mapped onto the smooth auxiliary wave function \( |\tilde{\Psi}_n\rangle \).
- The transformation operator \( T \) is defined as:
\[
T = 1+\sum_i(|\phi_i\rangle - |\tilde{\phi}_i\rangle)\langle\tilde{p}_i|
\]
- Here, \( |\phi_i\rangle \) and \( |\tilde{\phi}_i\rangle \) are the real and auxiliary partial wave functions respectively, and \( \langle\tilde{p}_i| \) is the projection operator.
2. **Expectation - value calculation**:
- The expectation value can be calculated using either the real wave function or the auxiliary wave function. For an operator \( A \), the expectation value can be expressed as:
\[
\langle A\rangle=\sum_n f_n\langle\tilde{\Psi}_n|T^\dagger A T|\tilde{\Psi}_n\rangle+\sum_{n = 1}^{N_c}\langle\phi_c^n|A|\phi_c^n\rangle
\]
- Here, \( f_n \) is the valence - state occupancy number, and \( N_c \) is the number of core states.
3. **Total - energy calculation**:
- The total energy \( E \) can be decomposed into a smooth part and an augmented part:
\[
E([\tilde{\Psi}_n], R_i)=\tilde{E}+\sum_R(E_1^R-\tilde{E}_1^R)
\]
- Smooth part