<div dir="ltr"><div><div><div>Dear PWSCF users, <br><br></div>I am trying to realize a cp 
dynamics of a Ca2+ + 2Cl- + 64 H2O solution, and I successfully realized
 NVE dynamics, and NVT dynamics with a single thermostat ( nhpcl  = 1,  
 nhptyp = 0,). However, when I try with more thermostats (nhpcl>1 
and/or nhptyp>0), the energy systematically  diverges. More 
precisely, the quantity called "vnhp" is the first to diverge and the 
one showing the more obvious divergency. Would you have any hints on 
what to do? I would be very grateful if you could guide me on which 
parameters may have the most effect on this divergency? <br><br></div>Below is a typical input of this divergent computation<br><br> &control<br>       calculation = 'cp',<br>      restart_mode = 'reset_counters',<br>            prefix = 'CaCl2-62wrl',<br>           disk_io = 'default' ,<br>    pseudo_dir     = '$WORKDIR',<br>    outdir         = '${WORKDIR}', <br>    tprnfor        = .true.,<br>    tstress        = .true.,<br>    dt             =  5.0,<br>    ndr=53,<br>    ndw = 54,<br>    nstep          =  30000,<br>    iprint         =  30,<br>    isave          =  30,<br>    etot_conv_thr = 1.d-4,<br>    ekin_conv_thr = 1.d-4,<br>    forc_conv_thr = 1.d-3,<br>/&end<br> &system<br>    ibrav = 0, celldm(1)=23.4752,<br>    nat = 189, ntyp = 4, ecutwfc = 80.0, ecutrho=320.0,<br>    nr1b=20, nr2b=20,nr3b=20, <br>/&end<br> &electrons<br>    emass = 400.d0,<br>    emass_cutoff = 2.5d0,<br>    orthogonalization = 'ortho',<br>    ortho_eps = 5.d-8,<br>    ortho_max = 25,<br>    electron_dynamics = 'verlet', ! 'damp', !'verlet',<br>    electron_velocities='default',<br>    electron_temperature='not_cont<wbr>rolled',<br> /<br> &ions<br>    ion_dynamics = 'verlet',<br>    ion_radius(1) = 1.0d0,<br>    ion_radius(2) = 1.0d0,<br>    ion_radius(3) = 0.8d0,<br>    ion_radius(4) = 0.5d0,<br>    ion_velocities = 'default',<br>    ion_temperature = 'nose',<br>    tempw=300.0,<br>    fnosep =  60.0,<br>    nhpcl  = 2,<br>    nhptyp = 1, <br> /<br> &cell<br>    cell_dynamics = 'none',<br>    cell_velocities = 'zero',<br>    press = 0.0d0,<br>    wmass = 70000.0d0<br> /<br>ATOMIC_SPECIES<br>  Ca   39.9625   Ca.blyp-spn-rrkjus_psl.1.0.0.U<wbr>PF<br>  Cl   34.9689   Cl.blyp-nl-rrkjus_psl.1.0.0.UP<wbr>F<br>  O    15.9949   O.blyp.UPF            <br>  H     1.0079   H.blyp2.UPF <br><br>CELL_PARAMETERS (alat= 23.47520000)<br>   1.000000000   0.000000000   0.000000000<br>   0.000000000   1.000000000   0.000000000<br>   0.000000000   0.000000000   1.000000000<br><br>ATOMIC_POSITIONS (bohr)<br>(...)<br><br></div>Thanks in advance, <br clear="all"><div><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Merlin Méheut<br>adresse labo:<br>GET - OMP   - Université Paul Sabatier <br>14 avenue Edouard Belin<br>31400 Toulouse<br>FRANCE <br>tel: (+33) 5 61 33 26 17<br><br></div></div>
</div></div>