[Wien] constrain some atoms to minimize the forces

Laurence Marks laurence.marks at gmail.com
Thu Jul 14 14:43:31 CEST 2016


Your question is based upon a common misunderstanding of minimizing forces.
You are thinking that the time to optimize depends upon the number of
positions (3*Number of Pt atoms that move) which is wrong. With PORT it
scales as the number of clusters of phonon frequencies and weakly with the
width of the clusters; with MSR1a it scales as the number of clusters of
eigenvalues of the joint electron-phonon problem and their width.

You will gain far more by generating a proper symmetric model with P-3m1
(maybe P-31m) symmetry rather than what you appear to have, a P1 structure.
This will reduce the number of uniques phonon etc modes by something like
6-12, and will be much faster and converge better.

---
Professor Laurence Marks
"Research is to see what everybody else has seen, and to think what nobody
else has thought", Albert Szent-Gyorgi
http://www.numis.northwestern.edu
Corrosion in 4D http://MURI4D.numis.northwestern.edu
Partner of the CFW 100% gender equity project, www.cfw.org/100-percent
Co-Editor, Acta Cryst A

On Jul 14, 2016 06:25, "Saleem Ayaz" <saleemayaz_hu at yahoo.com> wrote:

Dear users
I want to minimize the forces for Pt111 surface. I want to restrict the
structure by minimizing the forces for first two layers and keeping the
other layers bulk-like. Please some one can help me how to change case.inM
file. Should I make some numbers equal to zero?
Looking forward

regards
Saleem

PORT 2.00 0.35      # PORT/NEWT;  tolf, Initial Trust Radius
1.0 1.0 1.0 1.0   #Atom    1 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom    2 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom    3 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom    4 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom    5 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom    6 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom    7 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom    8 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom    9 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   10 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   11 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   12 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   13 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   14 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   15 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   16 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   17 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   18 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   19 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   20 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   21 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   22 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   23 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   24 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   25 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   26 Generated by pairhess
1.0 1.0 1.0 1.0   #Atom   27 Generated by pairhess
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20160714/e5630849/attachment.html>


More information about the Wien mailing list