[Pw_forum] something wrong with vc-relax & relax on my clusters
lewvega
vegalew at hotmail.com
Tue Jun 24 11:12:12 CEST 2008
Dear Users
I'm a new user of pwscf, and I have some experience of castep, which seems to be similar with pwscf.
First I tested paralell vision of pwscf with anatase-TiO2 lattice cell optimization calculation. the input file as follows:
&CONTROL
title = 'Anatase lattice BFGS' ,
calculation = 'vc-relax' ,
restart_mode = 'from_scratch' ,
outdir = '/home/vega/tmp/' ,
pseudo_dir = '/home/vega/espresso-4.0/pseudo/' ,
prefix = 'Anatase lattice default' ,
etot_conv_thr = 0.000000735 ,
forc_conv_thr = 0.0011668141375 ,
nstep = 1000 ,
/
&SYSTEM
ibrav = 6,
celldm(1) = 7.135605333,
celldm(3) = 2.5121822033898305084745762711864,
nat = 12,
ntyp = 2,
ecutwfc = 25 ,
ecutrho = 200 ,
/
&ELECTRONS
conv_thr = 7.3D-8 ,
/
&IONS
/
&CELL
/
ATOMIC_SPECIES
Ti 47.86700 Ti.pw91-sp-van_ak.UPF
O 15.99940 O.pw91-van_ak.UPF
ATOMIC_POSITIONS angstrom
Ti 0.000000000 0.000000000 0.000000000
Ti 1.888000000 1.888000000 4.743000000
Ti 0.000000000 1.888000000 2.372000000
Ti 1.888000000 0.000000000 7.115000000
O 0.000000000 0.000000000 1.973000000
O 1.888000000 1.888000000 6.716000000
O 0.000000000 1.888000000 4.345000000
O 1.888000000 0.000000000 9.088000000
O 1.888000000 0.000000000 5.141000000
O 0.000000000 1.888000000 0.398000000
O 1.888000000 1.888000000 2.770000000
O 0.000000000 0.000000000 7.513000000
K_POINTS automatic
7 7 3 1 1 1
then I checked the initial structure in Xcrysden, and type paralell command like this: mpiexec -n 20 pw.x < inputfile > outputfile
the calculation ran immediately. After calculation finished, the last CELL_PARAMETERS section in outputfile is
CELL_PARAMETERS (alat)
1.000000000 0.000000000 0.000000000
0.000000000 1.000000000 0.000000000
0.000000000 0.000000000 2.512182203
this indicate the 'vc-relax' didn't take effect, no cell parameters was not optimized during the vc-relax strategy.
And the atom positons also changed a lot, the ATOMIC_POSITIONS sections like this:
ATOMIC_POSITIONS (angstrom)
Ti 0.000000000 0.000000000 0.000108235
Ti 1.888000000 1.888000000 4.743111514
Ti 0.000000000 1.888000000 2.371065034
Ti 1.888000000 0.000000000 7.115503582
O 0.000000000 0.000000000 1.981862393
O 1.888000000 1.888000000 6.726003352
O 0.000000000 1.888000000 4.353977562
O 1.888000000 0.000000000 9.096778202
O 1.888000000 0.000000000 5.132185229
O 0.000000000 1.888000000 0.389375062
O 1.888000000 1.888000000 2.760157286
O 0.000000000 0.000000000 7.504294082
Both x,y remain as the initial value.
All of this indicate my above calculation is a 'relax' calculation not a 'vc-relax' strategy.
Also if I set the ionic and cell optimization both using BFGS schedule as the castep did, the cell parameters changed so much that an error occured said 'non orthogonal operation'.
Does anyone encounter the same problem?
Second, I uesed the optimized anatase cell via pwscf's default vc-relax damped MD schedule to cleave the 101 surface and make the further calculation. But another serious problem happened. the input file as follows:
&CONTROL
title = 'Anatase101_2x' ,
calculation = 'relax' ,
restart_mode = 'from_scratch' ,
outdir = '/home/vega/tmp/' ,
pseudo_dir = '/home/vega/espresso-4.0/pseudo/' ,
prefix = 'Anatase101_2x' ,
etot_conv_thr = 0.000000735 ,
forc_conv_thr = 0.0011668141375 ,
nstep = 1000 ,
/
&SYSTEM
ibrav = 8,
celldm(1) = 19.293941717,
celldm(2) = 1.109509628630600744472407458828,
celldm(3) = 1.7484889932972017886534804648227,
nat = 72,
ntyp = 2,
ecutwfc = 30 ,
ecutrho = 240 ,
/
&ELECTRONS
conv_thr = 7.3D-8 ,
/
&IONS
ion_dynamics = 'bfgs' ,
/
ATOMIC_SPECIES
Ti 47.86700 Ti.pw91-sp-van_ak.UPF
O 15.99940 O.pw91-van_ak.UPF
ATOMIC_POSITIONS angstrom
Ti 2.649000000 -1.888000000 0.733000000
Ti -1.060000000 0.000000000 4.242000000
Ti 0.446000000 0.000000000 1.610000000
Ti -3.264000000 -1.888000000 5.119000000
O 0.807000000 -1.888000000 1.466000000
O -2.902000000 0.000000000 4.975000000
O -1.397000000 0.000000000 2.344000000
O -5.105000000 -1.888000000 5.852000000
O -1.421000000 -1.888000000 4.386000000
O 2.287000000 0.000000000 0.877000000
O 0.783000000 0.000000000 3.508000000
O -4.323000000 -1.888000000 3.509000000
Ti 2.649000000 5.664000000 0.733000000
Ti -1.060000000 7.552000000 4.242000000
Ti 0.446000000 7.552000000 1.610000000
Ti -3.264000000 5.664000000 5.119000000
O 0.807000000 5.664000000 1.466000000
O -2.902000000 7.552000000 4.975000000
O -1.397000000 7.552000000 2.344000000
O -5.105000000 5.664000000 5.852000000
O -1.421000000 5.664000000 4.386000000
O 2.287000000 7.552000000 0.877000000
O 0.783000000 7.552000000 3.508000000
O -4.323000000 5.664000000 3.509000000
Ti 2.649000000 1.888000000 0.733000000
Ti -1.060000000 3.776000000 4.242000000
Ti 0.446000000 3.776000000 1.610000000
Ti -3.264000000 1.888000000 5.119000000
O 0.807000000 1.888000000 1.466000000
O -2.902000000 3.776000000 4.975000000
O -1.397000000 3.776000000 2.344000000
O -5.105000000 1.888000000 5.852000000
O -1.421000000 1.888000000 4.386000000
O 2.287000000 3.776000000 0.877000000
O 0.783000000 3.776000000 3.508000000
O -4.323000000 1.888000000 3.509000000
Ti 4.045000000 -1.888000000 4.242000000
Ti 1.842000000 0.000000000 5.119000000
O 2.204000000 -1.888000000 4.975000000
O 0.000000000 0.000000000 5.852000000
O 3.684000000 0.000000000 4.386000000
Ti 4.045000000 5.664000000 4.242000000
Ti 1.842000000 7.552000000 5.119000000
O 2.204000000 5.664000000 4.975000000
O 0.000000000 7.552000000 5.852000000
O 3.684000000 7.552000000 4.386000000
Ti 4.045000000 1.888000000 4.242000000
Ti 1.842000000 3.776000000 5.119000000
O 2.204000000 1.888000000 4.975000000
O 0.000000000 3.776000000 5.852000000
O 3.684000000 3.776000000 4.386000000
Ti 7.754000000 3.776000000 0.733000000
Ti 5.549000000 1.888000000 1.611000000
O 5.911000000 3.776000000 1.467000000
O 3.708000000 1.888000000 2.344000000
O 7.392000000 1.888000000 0.877000000
O 9.596000000 3.776000000 -0.000000000
O 4.490000000 1.888000000 0.000000000
Ti 7.754000000 7.552000000 0.733000000
Ti 5.549000000 5.664000000 1.611000000
O 5.911000000 7.552000000 1.467000000
O 3.708000000 5.664000000 2.344000000
O 7.392000000 5.664000000 0.877000000
O 9.596000000 7.552000000 -0.000000000
O 4.490000000 5.664000000 0.000000000
Ti 7.754000000 0.000000000 0.733000000
Ti 5.549000000 -1.888000000 1.611000000
O 5.911000000 0.000000000 1.467000000
O 3.708000000 -1.888000000 2.344000000
O 7.392000000 -1.888000000 0.877000000
O 9.596000000 0.000000000 -0.000000000
O 4.490000000 -1.888000000 0.000000000
K_POINTS gamma
after checking the initial structure via Xcrysden, I typed 'mpiexec -n 20 pw.x < inputfile > outputfile'. But unfortunately nothing happened, even the first few lines of pwscf was not printed, and all the processes of pw.x triggered by mpich2 are sleeping. No calculation is running at all. Then I run this task in serial version of pwscf by typing 'pw.x < inputfile > outputfile', the calculation ran immediately. then I tried again and again, the paralell calculation of the input file always failed with no message out. And to make sure whether the environment is suitable for paralell calculation, I ran the first input file again in 20 paralell process. Everything seems to be ok with the paralell environment.
I had tested the same task(same number of atoms, same atom, same convergence threshold et al.) in CASTEP successfully, so the error is not because of hardware.
Is there any one can try my input file, and give me some suggestions. thank you so much for reading.
by the way, I also will tell you something about my hardware and softeware:
hardware: five computers with intel coreII Q6600 CPU (quadcore) and 40G menmory (8G memory per computer) on intel S3000AH system board. Computers linked with 1Gbit/s network by only one Swich( TP-Link )
software: Readhat linux enterprise 4 Advanced Server Update 4, mpich2, intel fortran compiler 10.1.008 (64 bit-edition), intel c++ 10.1.008 (64 bit-edition), intel MKL 10.0.011 (64 bit-edition), fftw2.1.
And I the fftw in MKL can't be recognized by 'configure' of espresso, even the lib dictionary is included in LD_LIBRARY_PATH.
thank you all again for reading and I'm looking forward to your suggestions.
Vega Lew
PH.D Candidate in Chemical Engineering
College of Chemistry and Chemical Engineering
Nanjing University of Technology, 210009, Nanjing, Jiangsu, China
_________________________________________________________________
多个邮箱同步管理，live mail客户端万人抢用中
http://get.live.cn/product/mail.html
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20080624/51ea2d1e/attachment.html>
More information about the users
mailing list