Question about Geometry Optimization |
- Date: 2015/08/22 18:31
- Name: Eike F. Schwier
<schwier@hiroshima-u.ac.jp>
- Dear "MD.type RF" experts,
I have come across a rather strange behaviour during one of my bulk unit-cell optimisations and would like to ask, if anybody can explain to my why this is occurring and if there are any parameters I may change / set / optimize to prevent this from happening.
In my relaxation of the FRAC coordinates (RF) of TlBiSe2, I found that there was a lateral force on the Se atoms (not on Tl or Bi). If Se positions are allowed to relax laterally (MD.Fixed.XYZ 0 0 0) all Se move ~0.1 A away from their original positions until the convergence criterion is reached. This is due to a force Fx,Fy~0.008 Ha/bohr. If I the restart the optimisation from the "relaxed" coordinates, the Se atoms move back to (almost) their original positions before the first lateral relaxation. This seemed puzzling as I do not expect any more relaxation, if I restart the RF from its already relaxed coordinates. Btw: If I choose EF instead of RF I get the same behaviour.
Se positions in Ang (original, first relax, second relax) 1,x 1,y 2,x 2,y 3,x 3,y 0 0 0 0 2.12 1.22398 0.10537 0.04579 -0.10537 -0.04578 2.22532 1.26981 0.00031 0.00096 -0.00031 -0.00096 2.12024 1.22495
4,x 4,y 5,x 5,y 6,x 6,y 2.12 1.22398 0 2.44797 0 2.44797 2.01463 1.17812 0.10537 2.49384 -0.10533 2.40213 2.11978 1.22314 0.00022 2.4488 -0.00024 2.447
I tried to change the basisset, the cutoff and the k-mesh, but the initial force on the Se atoms remains nonzero regardless of my choice of parameters. However, if I displace the first Se atom (dX >= 0.+0.00001001 FRAC) the lateral forces Fx and Fy become zero (Fx,Fy~0.00005 Ha/bohr) If I choose a dX <= 0.00001 forces are non-zero again. My suspicion is, that I break some internal symmetry in the calculation with the Se displacement larger than 0.00001 FRAC. But I do not understand why a broken symmetry should change the forces to such an extend. Furthermore, if the groundstate is indeed one with the Se atoms displaced and the symmetry broken, then why do the Se atoms move back to the symmetry positions if I restart the calculation? Is the initial symmetry conserved in the RF (EF) run, even if the fractional coordinates are changed to break it?
best regards, Eike
# # Definition of Atomic Species #
Species.Number 3 <Definition.of.Atomic.Species Se Se7.0-s3p2d1 Se_PBE13 Bi Bi8.0-s3p2d2f1 Bi_PBE13 Tl Tl8.0-s3p3d2f1 Tl_PBE13 Definition.of.Atomic.Species>
Atoms.Number 12 Atoms.SpeciesAndCoordinates.Unit FRAC # Ang|AU <Atoms.SpeciesAndCoordinates 1 Se +0.00001500 -0.00000000 +0.24202732 3 3 2 Se -0.00000000 -0.00000000 +0.75797044 3 3 3 Se +0.66666666 +0.33333333 +0.57535558 3 3 4 Se +0.66666666 +0.33333333 +0.09131124 3 3 5 Se +0.33333333 +0.66666666 +0.90868607 3 3 6 Se +0.33333333 +0.66666666 +0.42464308 3 3 7 Bi -0.00000000 -0.00000000 +0.50000000 7.5 7.5 8 Bi +0.66666666 +0.33333333 +0.83333333 7.5 7.5 9 Bi +0.33333333 +0.66666666 +0.16666666 7.5 7.5 10 Tl -0.00000000 -0.00000000 -0.00000000 9.5 9.5 11 Tl +0.66666666 +0.33333333 +0.33333333 9.5 9.5 12 Tl +0.33333333 +0.66666666 +0.66666666 9.5 9.5 Atoms.SpeciesAndCoordinates>
Atoms.UnitVectors.Unit Ang # Ang|AU <Atoms.UnitVectors 4.24 0 0 -2.12 3.67194771204602 0 0 0 22.33000 Atoms.UnitVectors>
<MD.Fixed.XYZ 1 1 1 0 2 1 1 0 3 1 1 0 4 1 1 0 5 1 1 0 6 1 1 0 7 1 1 1 8 1 1 1 9 1 1 1 10 1 1 1 11 1 1 1 12 1 1 1 MD.Fixed.XYZ>
# # SCF or Electronic System #
scf.XcType GGA-PBE # LDA|LSDA-CA|LSDA-PW|GGA-PBE scf.ElectronicTemperature 300.0 # default=300 (K) scf.energycutoff 300.0 # default=150 (Ry) scf.maxIter 160 # default=40 scf.EigenvalueSolver band # Recursion|Cluster|Band scf.Kgrid 5 5 3 # means n1 x n2 x n3 scf.Generation.Kpoint MP # regular|MP scf.Mixing.Type rmm-diisk # Simple|Rmm-Diis|Gr-Pulay scf.Init.Mixing.Weight 0.3 # default=0.30 scf.Min.Mixing.Weight 0.001 # default=0.001 scf.Max.Mixing.Weight 0.600 # default=0.40 scf.Mixing.History 20 # default=5 scf.Mixing.StartPulay 10 # default=6 scf.criterion 1e-9 # default=1.0e-6 (Hartree) scf.restart off
# # MD or Geometry Optimization #
MD.Type RF # MD.maxIter 1 # default=1 MD.Opt.criterion 1.0e-5
scf.fixed.grid 0.00000000000000 0.00000000000000 0.00000000000000
![メンテ](./img/mente.gif)
| |