Requests for technical support from the VASP group should be posted in the VASP-forum.

Category:Interface pinning

From Vaspwiki
Jump to navigation Jump to search

Interface pinning[1] is used to determine the melting point from a molecular-dynamics simulation of the interface between a liquid and a solid phase. The typical behavior of such a simulation is to freeze or melt, while the interface is pinned with a bias potential. This potential applies an energy penalty for deviations from the desired two-phase system. It is preferred simulating above the melting point because the bias potential prevents melting better than freezing.

The Steinhardt-Nelson[2] order parameter discriminates between the solid and the liquid phase. With the bias potential

penalizes differences between the order parameter for the current configuration and the one for the desired interface . is an adjustable parameter determining the strength of the pinning.

Under the action of the bias potential, the system equilibrates to the desired two-phase configuration. An important observable is the difference between the average order parameter in equilibrium and the desired order parameter . This difference relates to the the chemical potentials of the solid and the liquid phase

where is the number of atoms in the simulation.

Computing the forces requires a differentiable . In the VASP implementation a smooth fading function is used to weight each pair of atoms at distance for the calculation of the order parameter. This fading function is given as

Here and are the near- and far-fading distances, respectively. The radial distribution function of the crystal phase yields a good choice for the fading range. To prevent spurious stress, should be small where the derivative of is large. Set the near fading distance to the distance where goes below 1 after the first peak. Set the far fading distance to the distance where goes above 1 again before the second peak.

How to

Interface pinning uses the ensemble where the barostat only acts along the direction. This ensemble uses a Langevin thermostat and a Parrinello-Rahman barostat with lattice constraints in the remaining two dimensions. The solid-liquid interface must be in the - plane perpendicular to the action of the barostat.

Set the following tags for the interface pinning method:

Defines the near-fading distance .
Defines the far-fading distance .
Defines the coupling strength of the bias potential.
Defines the desired value of the order parameter .

The following example INCAR file calculates the interface pinning in sodium[1]:

TEBEG = 400                   # temperature in K
POTIM = 4                     # timestep in fs
IBRION = 0                    # run molecular dynamics
ISIF = 3                      # use Parrinello-Rahman barostat for the lattice
MDALGO = 3                    # use Langevin thermostat
LANGEVIN_GAMMA_L = 3.0        # friction coefficient for the lattice degree of freedoms (DoF)
LANGEVIN_GAMMA = 1.0          # friction coefficient for atomic DoFs for each species
PMASS = 100                   # mass for lattice DoFs
LATTICE_CONSTRAINTS = F F T   # fix x-y plane, release z lattice dynamics
OFIELD_Q6_NEAR = 3.22         # near fading distance for function w(r) in Angstrom
OFIELD_Q6_FAR = 4.384         # far fading distance for function w(r) in Angstrom
OFIELD_KAPPA = 500            # strength of bias potential in eV/(unit of Q)^2
OFIELD_A = 0.15               # desired value of the Q6 order parameter



Pages in category "Interface pinning"

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