Category:Density mixing: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
Line 1: Line 1:
In each iteration of a DFT cycle, we start from a given charge density <math>\rho_{in}</math> and obtain the corresponding Kohn-Sham Hamiltonian and its eigenstates (wavefunctions). From the occupied states, we can compute a new charge density <math>\rho_{out}</math>, so that conceptionally we need to solve a multidimensional fixed-point problem. In the case of magnetism and MetaGGAs, the problem becomes even more complex, because in addition also the magnetization density and the kinetic-energy density are relevant.


== Charge Density ==
To solve this problem, we use nonlinear solvers that work with the input vector <math>\rho_{in}</math> and the residual <math>R = \rho_{out} - \rho_{in}</math>. In these methods a subspace is built from the input vectors and the optimal solution within this subspace is obtained. The most efficient solutions are the Broyden mixing and the Pulay mixing scheme. In the former method, an approximate of the Hessian matrix is iteratively improved to find the optimal solution. In the latter method, the input vectors are combined assuming linearity to minimize the residual.
 
 
== Magnetization Density ==
 
 
== Kinetic-energy density ==


The implementation in VASP is based on the





Revision as of 11:00, 10 May 2019

In each iteration of a DFT cycle, we start from a given charge density and obtain the corresponding Kohn-Sham Hamiltonian and its eigenstates (wavefunctions). From the occupied states, we can compute a new charge density , so that conceptionally we need to solve a multidimensional fixed-point problem. In the case of magnetism and MetaGGAs, the problem becomes even more complex, because in addition also the magnetization density and the kinetic-energy density are relevant.

To solve this problem, we use nonlinear solvers that work with the input vector and the residual . In these methods a subspace is built from the input vectors and the optimal solution within this subspace is obtained. The most efficient solutions are the Broyden mixing and the Pulay mixing scheme. In the former method, an approximate of the Hessian matrix is iteratively improved to find the optimal solution. In the latter method, the input vectors are combined assuming linearity to minimize the residual.

The implementation in VASP is based on the



Pages in category "Density mixing"

The following 12 pages are in this category, out of 12 total.