Category:Forces: Difference between revisions

From VASP Wiki
Line 111: Line 111:


So, the pressure is the average of the normal components of the stress tensor in the three spatial directions. In electronic structure calculations, finite basis sets are used to express
So, the pressure is the average of the normal components of the stress tensor in the three spatial directions. In electronic structure calculations, finite basis sets are used to express
the electron density. Due to this finiteness of the basis set errors on the stress tensor and the pressure are introduced. The error in the pressure is refered to as [[https://www.vasp.at/wiki/index.php/Energy_vs_volume_Volume_relaxations_and_Pulay_stress|Pulay stress]].
the electron density. Due to this finiteness of the basis set errors on the stress tensor and the pressure are introduced. The error in the pressure is refered to as [[https://www.vasp.at/wiki/index.php/Energy_vs_volume_Volume_relaxations_and_Pulay_stress| Pulay stress]].


'''Force constants and phonons'''
'''Force constants and phonons'''

Revision as of 13:54, 19 October 2023

On this page the computation of forces acting on ions will be presented.

Introduction

Forces on particles are a fundamental concept in condensed matter physics and chemistry. These forces describe the interactions that cause particles, such as atoms and molecules, to move and behave in specific ways. In VASP forces result from electromagnetic interactions which can be computed from DFT by the use of the Hellmann-Feynman theorem, the random-phase approximation or by the use of machine learning force fields. Understanding these forces is crucial in many aspects of science, as for example:

  • predicting the atomic structure of solids and molecules
  • to engineer and design new materials
  • predicting and optimizing chemical reactions
  • improving and understanding catalysis
  • predicting and understanding thermodynamic proerties

Formally the force can be defined as follows. Let be the position of the particle, then the velocity is defined as the change of position with time

and the momentum of the particle is the velocity times the particle mass m

Newton's second law of motion states that the change of motion of an object is proportional to the force acting on the object and oriented in the same direction as the force vector. Therefore the force is defined as the change of particle momentum with time

where is the acceleration of the particle. With this equation of motion, the knowledge of some starting conditions and and an algorithm to compute the forces the trajectory of a particle can be predicted for all times.

Theory

There is an important relation between forces and the negative gradient of the potential energy which can be computed from the Lagrangian of the particle system of interest. The Lagrangian for a N particle system is

where is the potential energy of the system. With Lagrange's equation of the second kind

the relation

Therefore to predict forces and particle trajectories a way to compute the negative gradient of the potential energy has to be established.


DFT Forces

One way to compute the potential energy's negative gradient is through DFT. In DFT there is no classical potential energy function but a Hamiltonian depending on the ionic positions and the electronic positions . The exact form of the Hamiltonian is given by

where denotes the electronic ground state density and are the Kohn-Sham orbitals. is the exchange correlation energy. To obtain the force acting on ion A the Hellmann-Feynman theorem has to be used.

where denotes the gradient with respect to ionic position . The DFT forces will depend on the chosen exchange correlation functional via the electronic ground state density . Therefore the choice of the proper exchange correlation functional for the system of interest is crucial for obtaining proper forces and hence the correct material properties.

RPA-forces

To obtain more accurate forces the RPA approximation can be used to get better estimates for the exchange correlation energy. The RPA adds the following correlation energy to the DFT ground state energy

.

With this, the Helmann-Feynman theorem can be rewritten as

Note that the RPA is a correction to the underlying DFT approximation of the forces. Therefore the choice of the proper exchange correlation functional is still crucial in the RPA approach for obtaining forces.

Machine-learning forces

A speedy but less accurate approach for obtaining forces is the machine learning approach. In this approach, a machine learning model is first fitted to either the DFT or RPA approach. During training the MLFF the forces, energies, and stresses are fitted with respect to a DFT or RPA approach. In the case of the RPA the stress tensor is not computed. The machine learning approach will be an approximation to the underlying method against which it was fitted. The machine learning force field decomposes the total DFT energy into local atomic contributions depending on all atomic positions in the system. Therefore the force acting on ion A is computed by

where is the kernel matrix which can be found on the machine learning theory page. The kernel matrix as the local energies depends on the positions of all atoms in the actual atomic configuration.


Related concepts

Stress and pressure

The stress tensor provides valuable information about how forces are distributed throughout a material, both in magnitude and direction. It includes normal stresses, which act perpendicular to a given plane, and shear stresses, which act parallel to the plane. Together, these components allow engineers and scientists to predict how materials will behave under various conditions, such as tension, compression, or shear. The stress tensor can be computed from a viral theorem including pair forces or with a finite difference approach deforming the simulation box. Pressure, often denoted as P, is a scalar component of the stress tensor. It represents the normal force per unit area acting on a surface within the material. In the stress tensor, pressure is related to the diagonal components σxx, σyy, and σzz:

So, the pressure is the average of the normal components of the stress tensor in the three spatial directions. In electronic structure calculations, finite basis sets are used to express the electron density. Due to this finiteness of the basis set errors on the stress tensor and the pressure are introduced. The error in the pressure is refered to as [Pulay stress].

Force constants and phonons

The forces are defined by the negative gradient of the potential energy. The force constant matrix is defined by

and is therefore the gradient of the force. The force constant matrix is a fundamental concept in solid-state physics and materials science, especially in the context of understanding the vibrational properties of crystals. It is a mathematical representation of the interatomic forces and their interactions within a crystal lattice. This matrix is used to describe the relationships between atomic displacements and the resulting forces that occur in a crystal. By Fourier transforming the Force constant matrix the dynamical matrix is obtained. By computing the eigenvalues of the dynamic matrix on various reciprocal lattice points the phonon dispersion relation can be obtained. Understanding phonons is essential as they influence materials properties such as the electrical conductivity, thermal conductivity, and mechanical properties of materials.

How To


Contents

Subcategories

This category has only the following subcategory.

Pages in category "Forces"

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