
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">

<HTML><HEAD>

<META http-equiv=Content-Type content="text/html; charset=gb2312">

<META content="MSHTML 6.00.2900.3314" name=GENERATOR>

<STYLE></STYLE>

</HEAD>

<BODY bgColor=#ffffff>

<DIV><FONT face=Arial size=2>Dear all:</FONT></DIV>

<DIV><FONT face=Arial size=2>&nbsp;&nbsp;&nbsp; I performed a calculation using 

cp.x. However I always get "ortho went bananas" errors either using calculation 

= 'scf' or 'cp'. Previous post on this issue says it's convergence problem. My 

system converges well using pw.x. I totally have no idea about this bananas 

error, can someone help me on this? Thanks.</FONT></DIV>

<DIV><FONT face=Arial size=2>Best,</FONT></DIV>

<DIV><FONT face=Arial size=2>ZP</FONT></DIV>

<DIV><FONT face=Arial size=2>-----cp.in--------</FONT></DIV>

<DIV><FONT face=Arial size=2>&amp;CONTROL<BR>&nbsp; calculation = 

'scf',<BR>&nbsp; restart_mode= 'from_scratch',<BR>&nbsp; pseudo_dir&nbsp; = 

'/home/xuzp/bin/lib/pwscf_pot/',<BR>&nbsp; outdir&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

= './',<BR>&nbsp; prefix&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 'cp',<BR>&nbsp; 

max_seconds = 20000.,<BR>&nbsp; 

dt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 50.,<BR>&nbsp; 

nstep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 20,<BR>&nbsp; 

iprint&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 20,<BR>&nbsp; 

isave&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 20,<BR>&nbsp; 

ndr&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 90, !input_restart 

unit<BR>&nbsp; ndw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 91, 

!output_restart unit<BR>/<BR>&amp;SYSTEM<BR>&nbsp; 

ibrav&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 14,<BR>&nbsp; 

a&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 

40.000,<BR>&nbsp; b&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

=&nbsp; 2.4519514633,<BR>&nbsp; 

c&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 

20.000,<BR>&nbsp; cosbc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.,<BR>&nbsp; 

cosac&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.,<BR>&nbsp; 

cosab&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.,<BR>&nbsp; 

nat&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp; = 20,<BR>&nbsp; 

ntyp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; = 4,<BR>&nbsp; 

ecutwfc&nbsp;&nbsp;&nbsp;&nbsp;= 30.,<BR>&nbsp; ecutrho&nbsp;&nbsp;&nbsp;&nbsp; 

= 240.,</FONT></DIV>

<DIV><FONT face=Arial size=2>&nbsp; nspin&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 

1,<BR>&nbsp; nr1b&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 16,<BR>&nbsp; 

nr2b&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 16,<BR>&nbsp; 

nr3b&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 16,<BR>/&nbsp; 

<BR>&amp;ELECTRONS<BR>&nbsp; electron_maxstep = 100,<BR>&nbsp; 

conv_thr&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.D-6,<BR>&nbsp; 

electron_dynamics = 'damp',<BR>&nbsp; electron_damping&nbsp; = 0.2,<BR>&nbsp; 

electron_velocities = 'default',<BR>&nbsp; 

emass&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 

400.,<BR>&nbsp; emass_cutoff&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 2.5,<BR>&nbsp; 

orthogonalization = 'ortho',<BR>&nbsp; 

ortho_eps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.D-8,<BR>&nbsp; 

ortho_max&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 

100,<BR>/<BR>&amp;IONS<BR>&nbsp; ion_dynamics = 'verlet',<BR>&nbsp; 

ion_nstepe&nbsp;&nbsp; = 1,<BR>&nbsp; ion_temperature = 'nose',<BR>&nbsp; 

tempw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 300.D0,<BR>&nbsp; 

fnosep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.D0,<BR>&nbsp; 

nhpcl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1, <BR>&nbsp; 

nhptyp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0,<BR>&nbsp; 

ndega&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0,<BR>&nbsp; 

tolp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 100.D0,<BR>&nbsp; 

tranp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = .false.,<BR>&nbsp; remove_rigid_rot 

= .true.,<BR>/</FONT></DIV>

<DIV><FONT face=Arial size=2></FONT>&nbsp;</DIV>

<DIV><FONT face=Arial size=2>ATOMIC_SPECIES<BR>C 12.0107&nbsp; 

C.pbe-rrkjus.UPF<BR>H 1.00794&nbsp; H.pbe-rrkjus.UPF<BR>N 14.00674 

N.pbe-rrkjus.UPF<BR>O 15.9994&nbsp; O.pbe-rrkjus.UPF</FONT></DIV>

<DIV><FONT face=Arial size=2></FONT>&nbsp;</DIV>

<DIV><FONT face=Arial size=2>K_POINTS automatic<BR>&nbsp;1 8 1 0 0 

0</FONT></DIV>

<DIV><FONT face=Arial size=2></FONT>&nbsp;</DIV>

<DIV><FONT face=Arial size=2>ATOMIC_POSITIONS 

angstrom<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5.8598630590&nbsp;&nbsp;&nbsp;&nbsp; 

1.2261622260&nbsp;&nbsp;&nbsp; 10.0000000000&nbsp; 

<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4.4498769380&nbsp;&nbsp;&nbsp;&nbsp; 

1.2259421420&nbsp;&nbsp;&nbsp; 10.0000000000&nbsp; <BR></FONT><FONT face=Arial 

size=2>--------------------------------------------------------------</FONT></DIV>

<DIV><FONT face=Arial size=2>-- <BR>XU Zhiping, Ph. D<BR>Civil and Environmental 

Engineering<BR>MIT, Cambridge, MA 02139</FONT></DIV></BODY></HTML>