<div>Dear all,</div>
<div> </div>
<div>I am trying to do a structural relaxation on a relatively large system, which </div>
<div>has 240 atoms in one unit cell. </div>
<div> </div>
<div>I know it may sound crazy, but the crystal contains only C and H atoms,</div>
<div>so I thought that my cluster might be able to handle it.</div>
<div> </div>
<div>Indeed, I was able to run pwscf structural relaxation calculations. The 1st</div>
<div>iteration finished ok. However, at the end of the 2nd iteration (which also</div>
<div>finished ok), after the program printed out:</div>
<div>------------------------------------------</div>
<div>...</div>
<div>number of scf cycles = 2</div>
<div>number of bfgs steps = 1</div>
<div>energy old =  - 2084.xxx</div>
<div>energy new = - 2086.xxx</div>
<div>CASE: energy_new < energy_old<br clear="all">------------------------------------------</div>
<div>The run stopped without giving any information expect saying that "mpi</div>
<div>stopped without calling mpi_finalize" or sth like that.<br> </div>
<div>I repeated the calculation and it always stopped at the same place. </div>
<div>It is wierd since I had no problem before when running smaller systems</div>
<div>(say, 150 atoms).</div>
<div> </div>
<div>Could some one give me a hint on what could be the reason? Is it related</div>
<div>to the system size? If so, I have no idea why the pw.x program is able to</div>
<div>finish the first and the second iteration, but does not continue...</div>
<div> </div>
<div>Thanks in advance for any suggestions that you may have!</div>
<div><br>-Peter<br><br>********************************************<br>Peter Winey, Ph.D., DuPont<br><a href="mailto:peterwiney@gmail.com">peterwiney@gmail.com</a><br>******************************************** </div>