potfit wiki

open source force-matching

User Tools

Site Tools


models:analytic

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
Last revisionBoth sides next revision
models:analytic [2018/09/16 10:49] – ↷ Links adapted because of a move operation danielmodels:analytic [2019/02/23 09:48] daniel
Line 3: Line 3:
 ---- ----
  
-Analytic potentials can be enabled in //potfit// by compiling it with the ''apot'' switch.  +Analytic potentials use the parameters $\alpha_i$ as arguments to analytic function templates
-There are some additional options for analytic potentials like the chemical potentials or the smooth cutoff.+A well-known example for this would be the Lennard-Jones potential, which uses two parameters $\epsilon$ and $\sigma$ to define a potential:
  
-====  Global parameters  ====+$$V(r)=4\varepsilon\left[\left(\frac{\sigma}{r}\right)^{12}-\left(\frac{\sigma}{r}\right)^{6}\right]$$ 
  
-//potfit// supports global parameters, which can occur in different potentials. With this feature it is possible for any potential to "share" parameters with each other.  +{{ :analytic.png?300 }}
-To use them, you need a section called ''global n'' after the header in your potential file, where ''n'' is the number of global parameters.  +
-The format for the parameters is the same as for normal parameters''name value min max''+
-If you want to use a global parameter in a potential, use the name of the global parameter and add an exclamation mark.+
  
-<code potfit> +Using an analytic form for $V(r)$ allows to define the atomic interaction for all possible distances between atoms with very few parameters compared to tabulated potentialsThe problem of local parameters (see [[#Tabulated potentials]]) is also not present, parameters are usually global parametersIn certain cases it is even possible to attribute some physical meaning to individual parameters.
-#F 0 3 +
-#C Mg Zn +
-#I 0 0 0 +
-#E +
- +
-global 2 +
-glob1 1 0 2 +
-glob2 2 0 5 +
- +
-type lj +
-cutoff 8 +
-glob1! +
-glob2! +
- +
-type lj +
-cutoff 8 +
-glob2! +
-glob1! +
-</code> +
- +
-====  Smooth Cutoff  ==== +
- +
-The smooth cutoff function is very important for analytic potentials since it ensures that the  +
-potential and its gradient vanish at the desired cutoff radius. +
-Therefore the potential is multiplied with a cutoff function $\Psi$. +
- +
-$$V_{SC}(r)=\Psi\left(\frac{r-r_c}{h}\right)V(r)\qquad\text{where}\qquad\Psi(x)=\frac{x^4}{1+x^4}$+
- +
-To enable the smooth cutoff for a potential, ''_sc'' has to be added to the potential identifier. +
- +
-<code potfit> +
-type lj_sc +
-cutoff 7 +
-epsilon 0.1 0 1 +
-sigma 2.5 1 4 +
-h 1 0 2 +
-</code> +
- +
-Potentials fitted with the ''_sc'' option will have an additional parameter ''h'' at the end.  +
-If you omit it, it will be kept fixed at a value of 1. +
- +
-====  Implemented Potential Functions  ==== +
- +
-The currently implemented analytic potential functions can be found [[:analytic_functions|here]]. +
- +
-====  Adding New Potential Functions  ==== +
- +
-If you want to add other analytic potentials see [[:add_new_functions|this guide]].+
  
 +The analytic potential implementation in //potfit// has additional features which are only useful with this kind of interactions. They are described on the [[potfiles::format0|analytic format]] page.
models/analytic.txt · Last modified: 2019/02/23 09:52 by daniel