[Pw_forum] Relaxation of low symmetry lattices
Dal Corso Andrea
dalcorso at sissa.it
Thu Oct 8 10:52:52 CEST 2009
On Thu, 2009-10-08 at 10:37 +0200, Stefano Baroni wrote:
> Hi Manoj. the forces you are getting seem really small to me.
> If you enforce cubic symmetry, forces vanish exactly. If you do not,
> the code itself has to "figure out" that the atoms are indeed at their
> proper equilibrium positions, and this can only be achieved within the
> given overall numerical accuracy. One thing that may explicitly
> (although very slightly) break cubic symmetry in your supercell is tha
> FFT grid. Not that it makes any physical sense, but if you really want
> to overkill zero fores, you may want to play a little bit with the
> kinetic-energy cutoff and/or FFT grid. SB
I would also decrease the threshold of the self-consistency conv_thr.
HTH
Andrea
>
> On Oct 7, 2009, at 11:23 PM, Manoj Srivastava wrote:
>
> > Dear PWSCF users,
> > I have a question regarding relaxation of low symmetry lattices. I
> > did
> > relaxation of 5 atoms FCC unit cell of copper. This is pure copper,
> > no
> > interface, no defect or anything. The unit cell is just not
> > conventional,
> > but a different way of choosing bulk FCC. For the default value of
> > convergence the calculation finishes at the end of 1st scf cycle,
> > but I
> > see different forces when I describe nosym=.TRUE. or cell_parameters
> > {cubic}. Following is the input file.
> >
> > &control
> > calculation='relax'
> > pseudo_dir = '/home/manoj/espresso-4.0.4/pseudo',
> > outdir='./',
> > prefix='lcu',
> > wf_collect=.TRUE.
> > /
> > &system
> > ibrav =0,
> > celldm(1)=6.81650937063832
> > nat= 5,
> > ntyp= 1,
> > ecutwfc =50.0,
> > occupations='smearing',
> > smearing='gaussian',
> > degauss=0.02,
> > ecutrho=500
> > /
> > &electrons
> > conv_thr = 1.0e-8
> > mixing_beta=0.7
> > /
> > &ions
> > ion_dynamics='bfgs'
> > pot_extrapolation='second_order'
> > wfc_extrapolation='second_order'
> > /
> > ATOMIC_SPECIES
> > Cu 63.55 Cu.pz-d-rrkjus.UPF
> > ATOMIC_POSITIONS
> > Cu 0.0 0.0 0.0
> > Cu 0.5 0.5 0.0
> > Cu 0.0 1.0 0.0
> > Cu 1.0 1.0 0.0
> > Cu 0.5 1.5 0.0
> >
> > K_POINTS (automatic)
> > 8 8 8 0 0 0
> > CELL_PARAMETERS {cubic}
> > 1.5 0.5 0.0
> > -0.5 1.5 0.0
> > 0.0 0.5 0.5
> >
> > Result of relaxation : The covergence is achieved with follwing
> > forces on
> > the atoms.
> > Forces acting on atoms (Ry/au):
> >
> > Forces acting on atoms (Ry/au):
> >
> > atom 1 type 1 force = 0.00000000 0.00000000
> > 0.00000000
> > atom 2 type 1 force = 0.00000567 0.00000185
> > -0.00000172
> > atom 3 type 1 force = 0.00001155 0.00000377
> > -0.00000372
> > atom 4 type 1 force = -0.00001155 -0.00000377
> > 0.00000372
> > atom 5 type 1 force = -0.00000567 -0.00000185
> > 0.00000172
> >
> >
> > I realized that I should not have chosen cubic symmetry group, as my
> > unit
> > cell is not cubic. So, as described in the input file parameters,
> > for the
> > low symmetry lattices, I specifed nosym=.TRUE., and again the
> > relaxation
> > stopped after 1 scf cycle, acheiving convergence. But the forces on
> > the
> > atoms are -
> >
> > Forces acting on atoms (Ry/au):
> >
> > atom 1 type 1 force = 0.00001355 0.00000524
> > -0.00000516
> > atom 2 type 1 force = -0.00001493 -0.00000395
> > 0.00000395
> > atom 3 type 1 force = 0.00000253 -0.00000043
> > 0.00000045
> > atom 4 type 1 force = 0.00000820 0.00000230
> > -0.00000239
> > atom 5 type 1 force = -0.00000935 -0.00000316
> > 0.00000315
> >
> > So, still forces are .001 ev/a.u, which i think is higher for a
> > perfect
> > FCC lattice. I can specify stricter convergence of force, but that
> > would
> > move atoms around and break symmetry. So, as long as the lattice
> > parameter
> > is correct, I think for perfect FCC lattice, forces on atoms should
> > be
> > smaller than what they are presently. I have checked the lattice
> > constant
> > for bulk FCC with a different calculation, so lattice constant is
> > not an
> > issue. Another possibility is symmetry group, as you can see with
> > the same
> > input variables, changing symmetry group changes results.I was
> > wondering
> > if there is a way i can choose correct symmetry group for this low
> > symmetry lattice which might give me smaller forces? Are there any
> > other
> > parameters I need to specify? On a broad view I am curious to know
> > how to
> > do relaxation of low symmetry lattices.
> >
> > Regards,
> > Manoj Srivastava
> > University of Florida
> > Gainesville, USA
> >
> > _______________________________________________
> > Pw_forum mailing list
> > Pw_forum at pwscf.org
> > http://www.democritos.it/mailman/listinfo/pw_forum
> >
>
> ---
> Stefano Baroni - SISSA & DEMOCRITOS National Simulation Center -
> Trieste
> http://stefano.baroni.me [+39] 040 3787 406 (tel) -528 (fax) /
> stefanobaroni (skype)
>
>
> La morale est une logique de l'action comme la logique est une morale
> de la pensée - Jean Piaget
>
>
> Please, if possible, don't send me MS Word or PowerPoint attachments
> Why? See: http://www.gnu.org/philosophy/no-word-attachments.html
>
>
>
>
>
>
>
>
>
>
>
>
>
>
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
--
Andrea Dal Corso Tel. 0039-040-3787428
SISSA, Via Beirut 2/4 Fax. 0039-040-3787528
34151 Trieste (Italy) e-mail: dalcorso at sissa.it
More information about the Pw_forum
mailing list