p2ptrans¶
Command line¶
Usage:
p2ptrans [-h] [-I INITIAL_STRUCTURE] [-F FINAL_STRUCTURE] [-n N_CELLS] [-i] [-d] [-p FILENAME] [-c CRYSTFILE]
[-o OUTDIR] [-u USEDIR] [-m] [-s] [-r] [-a] [-v] [-t] [-f N_STEPS] [--version]
- Optional arguments:
-h, --help Show help message and exit -I, --initial INITIAL_STRUCTURE Initial crystal structure. It must be a POSCAR. -F, --final FINAL_STRUCTURE Final Structure structure. It must also be a POSCAR. -n, --ncell N_CELL Minimum number of cells to tile DEFAULT: 300 - -cn N_CELL
- If specified, number of cells used for the final mapping after tmat has been adjusted to be commensurate with both structures. If not specified, the mapping stays the same.
-i, --interactive Enable interactive display -d, --disp Save figures -p, --param FILENAME Name of the parameter file (see detail in Extra parameters). DEFAULT: p2p.in Caution: if a file named p2p.in exists in the current directory, p2ptrans will read it. -c, --crystal CRYSTFILE Name of the parameter file for crystallography analysis DEFAULT: cryst.in -o, --outdir OUTDIR Output directory DEFAULT: current directory -u, --use USEDIR Use previously calculated data in USEDIR -m, --minimize Force new optimization even if data is available -s, --switch Map the larger cell on the smaller cell instead of the opposite -r, --noprim Do not try to find the primitive cell of the input structures -a, --anim Produce an animation of the transition (slow) -v, --vol Make the two (stoichiometric) cells equal in volume (not tested) -t, --test Test the input file and prepare the run. You can continue this run with the -u USEDIR -r option -f N_STEPS, --frames N_STEPS Number of intermediate structures to create --version Show installed version of p2p
Extra Parameters¶
More parameters can be specified to p2ptrans via a parameter file. This file is in the Fortran namelist format. By default any file named p2p.in in the current directory will be read. The name of the parameter file can be changed using the -p option. The syntax of this file is the following:
&input
init_class = 1.0d0
tol = 1.0d-4
tol_std = tol*1.0d-3
tol_class = 1.0d-3
fracA = 0.09d0
fracB = 0.3d0
n_iter = 1000
n_ana = 300
n_conv = 5
n_class = 30
n_adjust = 10
max_vol = 4.0d0
free = .false.
check = .false.
savebest = trim(outdir)//"/best.dat"
usebest = .false.
remap = .true.
/
This is file contains all the default parameters, if an entry is not specified, it will take the value shown above.
- Init_class
- Initial separation tolerance for displacement classes. At the initial classification step, if the norm of the difference between two vectors is larger than init_class they will be classified in different groups.
- tol
- Convergence criterion for the gradient descent
- tol_std
- Convergence criterion for the std minimization
- tol_class
- Convergence criterion for the classification abs(std - previous std)
- fracA
- Fraction of the mapped structure that constitutes core atoms
- fracB
- Fraction of the mapping structure that constitutes mapping atoms
- n_iter
- Number of random starts
- n_ana
- Maximum number of iterations in the gradient descent
- n_conv
- Maximum number of remappings per minimization
- n_class
- Maximum number of classification iterations
- n_adjsut
- Maximum number of unconstrained post-processing minimization iterations
- max_vol
- Maximum volume of the random starting tmat when using the unrestricted minimization (free = .true.)
- free
- Use unrestricted minimization. Not limited to rigid rotations. The optimal result is the one for which the sum of the unstrained distance and the strained (unrestricted) distance is minimal so that overly stretched results, where skipping is likely to occur, are penalized.
- check
- Use the mapping given by the order of the input structures directly (do not map atoms).
- savebest
- Name of the file to save the optimal result to at the end of the minimization, before the post-processing steps.
- remap
- If true, allows remapping during the post-processing steps.