<div dir="ltr"><div><div><div><div><div><div>Dear PWSCF users, <br><br></div>I am trying to realize a molecular dynamics using cp.x, and I find some difficulties. In particular, the kinetic energy is slightly deriving, and I would like to cool down the electrons, using the cg electron_dynamics option, as proposed in the espresso-5.1/CPV/Doc/INPUT_CP.def: <br><br>             'cg'      : conjugate gradient is used to converge the<br>                         wavefunction at each ionic step. 'cg' can be used<br>                         interchangeably with 'verlet' for a couple of ionic<br>                         steps in order to "cool down" the electrons and<br>                         return them back to the Born-Oppenheimer surface.<br>                         Then 'verlet' can be restarted again. This procedure<br>                         is useful when electronic adiabaticity in CP is lost<br>                         yet the ionic velocities need to be preserved.<br><br><br></div><div>So I do molecular dynamics (with  electron_dynamics = 'verlet' ) then I change to 'cg'. <br></div>Unfortunately, the calculation is systematically stuck .  The last lines of the output files being: <br><br>-------------------------------------------------------------<br>...  <br>   reading restart file: /tmpdir/mmeheut/Ca++/CaCl2-62wrl_56.save<br>   restart file read in    2.062 sec.<br><br><br>   formf: eself=  1290.97722<br>   formf:     vps(g=0)=  -0.0010930     rhops(g=0)=  -0.0007593<br>   formf: sum_g vps(g)=  -0.7938444 sum_g rhops(g)=  -0.8983221<br>   formf:     vps(g=0)=  -0.0011291     rhops(g=0)=  -0.0005315<br>   formf: sum_g vps(g)=  -1.8515153 sum_g rhops(g)=  -0.6288255<br>   formf:     vps(g=0)=  -0.0005757     rhops(g=0)=  -0.0004585<br>   formf: sum_g vps(g)=  -1.5145512 sum_g rhops(g)=  -1.0525002<br>   formf:     vps(g=0)=  -0.0000296     rhops(g=0)=  -0.0000770<br>   formf: sum_g vps(g)=  -0.0789760 sum_g rhops(g)=  -0.7182652<br>   Delta V(G=0):   0.071172Ry,    1.936684eV<br> PERFORMING CONJUGATE GRADIENT MINIMIZATION OF EL. STATES<br>-------------------------------------------------------------<br></div><div><br></div><div>In the log file of the job I have the following message: <br></div><div> forrtl: severe (174): SIGSEGV, segmentation fault occurred</div><div><br></div>Here is my input: <br><br><br> &control<br>       calculation = 'cp',<br>      restart_mode = 'restart', ! 'reset_counters',<br>            prefix = 'CaCl2-62wrl',<br>           disk_io = 'default' ,<br>    pseudo_dir     = '/tmpdir/mmeheut/Ca++',<br>    outdir         = '/tmpdir/mmeheut/Ca++', <br>    tprnfor        = .true.,<br>    tstress        = .true.,<br>    dt             =  5.0, ! 0.2 fs <br>    ndr = 56,   ! 51: relaxe ini ; 51-52: heat50K; 52-53: heat150K ; 56: simmer10 <br>    ndw = 57,<br>    nstep          =  800,  <br>    iprint         =  10,<br>    isave          =  100,<br>    etot_conv_thr = 3.d-5,<br>    ekin_conv_thr = 1.d-6,<br>    forc_conv_thr = 3.d-4,<br>/&end<br> &system<br>    ibrav = 0, celldm(1)=23.4752,<br>    nat = 189, ntyp = 4, ecutwfc = 40.0, ecutrho=200.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 = 'cg',<br>    electron_velocities='default',<br>    electron_temperature='not_controlled',<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 = 'not_controlled',<br>    tempw=300.0,<br>    fnosep = 100.0,<br>    nhpcl  = 1,<br>    nhptyp = 0, <br> /<br> &cell<br>    cell_dynamics = 'none',<br>    cell_velocities = 'zero',<br>    press = 0.0d0,<br>    wmass = 70000.0d0<br> /<br><br>ATOMIC_SPECIES<br>  Ca   39.9625   ca_pbe_v1.uspp.F.UPF<br>  Cl   34.9689   cl_pbe_v1.4.uspp.F.UPF<br>  O    15.9949   o_pbe_v1.2.uspp.F.UPF <br>  H     1.0079   h_pbe_v1.4.uspp.F.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 (angstrom)<br>Ca       2.246613456   8.592421018   8.991482765<br>...<br><br></div>I tried with versions 5.1. and 5.3. and acheved the same result. Would you have any hints to help me move on? <br><br></div>Thank you in advance for your help, <br><br></div>Best regards, <br clear="all"><div><div><div><div><div><div><div><br>-- <br><div class="gmail_signature"><div dir="ltr">Merlin Méheut (associate professor)<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></div></div></div></div></div></div>