<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 9pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'>
Dear Users<br><br>I'm a new user of pwscf, and I have some experience of castep, which seems to be similar with pwscf.<br>First I tested paralell vision of pwscf with anatase-TiO2 lattice cell optimization calculation. the input file as follows:<br> &CONTROL<br> title = 'Anatase lattice BFGS' ,<br> calculation = 'vc-relax' ,<br> restart_mode = 'from_scratch' ,<br> outdir = '/home/vega/tmp/' ,<br> pseudo_dir = '/home/vega/espresso-4.0/pseudo/' ,<br> prefix = 'Anatase lattice default' ,<br> etot_conv_thr = 0.000000735 ,<br> forc_conv_thr = 0.0011668141375 ,<br> nstep = 1000 ,<br> /<br> &SYSTEM<br> ibrav = 6,<br> celldm(1) = 7.135605333,<br> celldm(3) = 2.5121822033898305084745762711864,<br> nat = 12,<br> ntyp = 2,<br> ecutwfc = 25 ,<br> ecutrho = 200 ,<br> /<br> &ELECTRONS<br> conv_thr = 7.3D-8 ,<br> /<br> &IONS<br> /<br> &CELL<br> /<br>ATOMIC_SPECIES<br> Ti 47.86700 Ti.pw91-sp-van_ak.UPF <br> O 15.99940 O.pw91-van_ak.UPF <br>ATOMIC_POSITIONS angstrom <br> Ti 0.000000000 0.000000000 0.000000000 <br> Ti 1.888000000 1.888000000 4.743000000 <br> Ti 0.000000000 1.888000000 2.372000000 <br> Ti 1.888000000 0.000000000 7.115000000 <br> O 0.000000000 0.000000000 1.973000000 <br> O 1.888000000 1.888000000 6.716000000 <br> O 0.000000000 1.888000000 4.345000000 <br> O 1.888000000 0.000000000 9.088000000 <br> O 1.888000000 0.000000000 5.141000000 <br> O 0.000000000 1.888000000 0.398000000 <br> O 1.888000000 1.888000000 2.770000000 <br> O 0.000000000 0.000000000 7.513000000 <br>K_POINTS automatic <br> 7 7 3 1 1 1 <br><br>then I checked the initial structure in Xcrysden, and type paralell command like this: mpiexec -n 20 pw.x < inputfile > outputfile<br><br>the calculation ran immediately. After calculation finished, the last CELL_PARAMETERS section in outputfile is <br>CELL_PARAMETERS (alat)<br> 1.000000000 0.000000000 0.000000000<br> 0.000000000 1.000000000 0.000000000<br> 0.000000000 0.000000000 2.512182203<br>this indicate the 'vc-relax' didn't take effect, no cell parameters was not optimized during the vc-relax strategy. <br><br>And the atom positons also changed a lot, the ATOMIC_POSITIONS sections like this:<br><br>ATOMIC_POSITIONS (angstrom)<br>Ti 0.000000000 0.000000000 0.000108235<br>Ti 1.888000000 1.888000000 4.743111514<br>Ti 0.000000000 1.888000000 2.371065034<br>Ti 1.888000000 0.000000000 7.115503582<br>O 0.000000000 0.000000000 1.981862393<br>O 1.888000000 1.888000000 6.726003352<br>O 0.000000000 1.888000000 4.353977562<br>O 1.888000000 0.000000000 9.096778202<br>O 1.888000000 0.000000000 5.132185229<br>O 0.000000000 1.888000000 0.389375062<br>O 1.888000000 1.888000000 2.760157286<br>O 0.000000000 0.000000000 7.504294082<br><br>Both x,y remain as the initial value.<br><br>All of this indicate my above calculation is a 'relax' calculation not a 'vc-relax' strategy.<br>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'. <br><br>Does anyone encounter the same problem?<br><br>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:<br> <br>&CONTROL<br> title = 'Anatase101_2x' ,<br> calculation = 'relax' ,<br> restart_mode = 'from_scratch' ,<br> outdir = '/home/vega/tmp/' ,<br> pseudo_dir = '/home/vega/espresso-4.0/pseudo/' ,<br> prefix = 'Anatase101_2x' ,<br> etot_conv_thr = 0.000000735 ,<br> forc_conv_thr = 0.0011668141375 ,<br> nstep = 1000 ,<br> /<br> &SYSTEM<br> ibrav = 8,<br> celldm(1) = 19.293941717,<br> celldm(2) = 1.109509628630600744472407458828,<br> celldm(3) = 1.7484889932972017886534804648227,<br> nat = 72,<br> ntyp = 2,<br> ecutwfc = 30 ,<br> ecutrho = 240 ,<br> /<br> &ELECTRONS<br> conv_thr = 7.3D-8 ,<br> /<br> &IONS<br> ion_dynamics = 'bfgs' ,<br> /<br>ATOMIC_SPECIES<br> Ti 47.86700 Ti.pw91-sp-van_ak.UPF <br> O 15.99940 O.pw91-van_ak.UPF <br>ATOMIC_POSITIONS angstrom <br> Ti 2.649000000 -1.888000000 0.733000000 <br> Ti -1.060000000 0.000000000 4.242000000 <br> Ti 0.446000000 0.000000000 1.610000000 <br> Ti -3.264000000 -1.888000000 5.119000000 <br> O 0.807000000 -1.888000000 1.466000000 <br> O -2.902000000 0.000000000 4.975000000 <br> O -1.397000000 0.000000000 2.344000000 <br> O -5.105000000 -1.888000000 5.852000000 <br> O -1.421000000 -1.888000000 4.386000000 <br> O 2.287000000 0.000000000 0.877000000 <br> O 0.783000000 0.000000000 3.508000000 <br> O -4.323000000 -1.888000000 3.509000000 <br> Ti 2.649000000 5.664000000 0.733000000 <br> Ti -1.060000000 7.552000000 4.242000000 <br> Ti 0.446000000 7.552000000 1.610000000 <br> Ti -3.264000000 5.664000000 5.119000000 <br> O 0.807000000 5.664000000 1.466000000 <br> O -2.902000000 7.552000000 4.975000000 <br> O -1.397000000 7.552000000 2.344000000 <br> O -5.105000000 5.664000000 5.852000000 <br> O -1.421000000 5.664000000 4.386000000 <br> O 2.287000000 7.552000000 0.877000000 <br> O 0.783000000 7.552000000 3.508000000 <br> O -4.323000000 5.664000000 3.509000000 <br> Ti 2.649000000 1.888000000 0.733000000 <br> Ti -1.060000000 3.776000000 4.242000000 <br> Ti 0.446000000 3.776000000 1.610000000 <br> Ti -3.264000000 1.888000000 5.119000000 <br> O 0.807000000 1.888000000 1.466000000 <br> O -2.902000000 3.776000000 4.975000000 <br> O -1.397000000 3.776000000 2.344000000 <br> O -5.105000000 1.888000000 5.852000000 <br> O -1.421000000 1.888000000 4.386000000 <br> O 2.287000000 3.776000000 0.877000000 <br> O 0.783000000 3.776000000 3.508000000 <br> O -4.323000000 1.888000000 3.509000000 <br> Ti 4.045000000 -1.888000000 4.242000000 <br> Ti 1.842000000 0.000000000 5.119000000 <br> O 2.204000000 -1.888000000 4.975000000 <br> O 0.000000000 0.000000000 5.852000000 <br> O 3.684000000 0.000000000 4.386000000 <br> Ti 4.045000000 5.664000000 4.242000000 <br> Ti 1.842000000 7.552000000 5.119000000 <br> O 2.204000000 5.664000000 4.975000000 <br> O 0.000000000 7.552000000 5.852000000 <br> O 3.684000000 7.552000000 4.386000000 <br> Ti 4.045000000 1.888000000 4.242000000 <br> Ti 1.842000000 3.776000000 5.119000000 <br> O 2.204000000 1.888000000 4.975000000 <br> O 0.000000000 3.776000000 5.852000000 <br> O 3.684000000 3.776000000 4.386000000 <br> Ti 7.754000000 3.776000000 0.733000000 <br> Ti 5.549000000 1.888000000 1.611000000 <br> O 5.911000000 3.776000000 1.467000000 <br> O 3.708000000 1.888000000 2.344000000 <br> O 7.392000000 1.888000000 0.877000000 <br> O 9.596000000 3.776000000 -0.000000000 <br> O 4.490000000 1.888000000 0.000000000 <br> Ti 7.754000000 7.552000000 0.733000000 <br> Ti 5.549000000 5.664000000 1.611000000 <br> O 5.911000000 7.552000000 1.467000000 <br> O 3.708000000 5.664000000 2.344000000 <br> O 7.392000000 5.664000000 0.877000000 <br> O 9.596000000 7.552000000 -0.000000000 <br> O 4.490000000 5.664000000 0.000000000 <br> Ti 7.754000000 0.000000000 0.733000000 <br> Ti 5.549000000 -1.888000000 1.611000000 <br> O 5.911000000 0.000000000 1.467000000 <br> O 3.708000000 -1.888000000 2.344000000 <br> O 7.392000000 -1.888000000 0.877000000 <br> O 9.596000000 0.000000000 -0.000000000 <br> O 4.490000000 -1.888000000 0.000000000 <br>K_POINTS gamma <br><br>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.<br>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. <br>Is there any one can try my input file, and give me some suggestions. thank you so much for reading.<br><br>by the way, I also will tell you something about my hardware and softeware:<br>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 )<br>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.<br>And I the fftw in MKL can't be recognized by 'configure' of espresso, even the lib dictionary is included in LD_LIBRARY_PATH.<br><br>thank you all again for reading and I'm looking forward to your suggestions.<br><br><br>Vega Lew<br>PH.D Candidate in Chemical Engineering<br>College of Chemistry and Chemical Engineering<br>Nanjing University of Technology, 210009, Nanjing, Jiangsu, China<br /><hr />用 Windows Live Spaces 展示个性自我，与好友分享生活！ <a href='http://spaces.live.com/?page=HP' target='_new'>了解更多信息！</a></body>
</html>