<div dir="ltr"><div>Dear Paolo and Lorrenzo,</div>Thanks so much for your analysis of this issue.  Â  Â The patch to setup.f90 in the 6.4.1 release does seem to fix this problem at least for the simple NaCl supercell.  Â  Â Thanks again.  Â  Natalie<div><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature">N. A. W. Holzwarth Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  email: <a href="mailto:natalie@wfu.edu" target="_blank">natalie@wfu.edu</a><div>Department of Physics Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â web: <a href="http://www.wfu.edu/~natalie" target="_blank">http://www.wfu.edu/~natalie</a></div><div>Wake Forest University Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  phone: 1-336-758-5510 </div><div>Winston-Salem, NC 27109 USA Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  office: Rm. 300 Olin Physical Lab</div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Apr 18, 2019 at 5:59 AM Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com">p.giannozzi@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">Oh well. I see: it's an unintended side effect of a bug fix to 6.3, present in 6.3-backports. <br></div><div dir="ltr">There is a simple solution, not guarantee to be unintended-side-effect-free, though: <br></div><div dir="ltr"><br></div><div dir="ltr">diff --git a/PW/src/setup.f90 b/PW/src/setup.f90<br>index 0918ac017..135e58789 100644<br>--- a/PW/src/setup.f90<br>+++ b/PW/src/setup.f90<br>@@ -519,7 +519,10 @@ SUBROUTINE setup()<br>   !<br>   ! ... nosym: do not use any point-group symmetry (s(:,:,1) is the identity)<br>   !<br>-  IF ( nosym ) nsym = 1<br>+  IF ( nosym ) THEN<br>+     nsym = 1<br>+     invsym = .FALSE.<br>+  END IF<br>   !<br>   IF ( nsym > 1 .AND. ibrav == 0 ) CALL infomsg('setup', &<br>        'DEPRECATED: symmetry with ibrav=0, use correct ibrav instead')<br></div><div><br></div><div>Paolo <br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Apr 18, 2019 at 11:11 AM Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com" target="_blank">paulatz@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">> symmetry detection?  Â  We did put nosym=.true. but the code seems to <br>
> have set invsym=.true. anyway. <br>
<br>
I see what you mean,<br>
with the exact same input version 6.3 set the variable invsym, from <br>
symm_base to false, while the current version sets it to true.<br>
<br>
This breaks start_therm because the only possible initial velocities <br>
that respect the constraint in this case are all zero.<br>
<br>
The bug is that version 6.3 sets invsym AFTER enforcing nosym=.true., <br>
while git version sets invsym BEFORE enforcing nosym=.true., the <br>
subroutine find_sym is identical<br>
<br>
The change seems to have occurred in commit 23fb73f9b7 inside <br>
PW/src/setup.f90 I'll let Paolo Giannozzi have a look before going <br>
deeper, as I do not immediately understand what has changed, I think it <br>
has something to do with nrot_ vs. nrot, as the former is 1 but the <br>
latter is 48 before the call to find_sym in setup.f90 line 513, but <br>
find_sym internally uses nrot.<br>
<br>
cheers<br>
<br>
<br>
Perhaps it is intended that users should<br>
> move atoms off of their symmetry positions before running MD <br>
> simulations?  Â  Â  Â Thanks for your advice about this.  Â  Natalie<br>
> <br>
> N. A. W. Holzwarth Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  email: <br>
> <a href="mailto:natalie@wfu.edu" target="_blank">natalie@wfu.edu</a> <mailto:<a href="mailto:natalie@wfu.edu" target="_blank">natalie@wfu.edu</a>><br>
> Department of Physics Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â web: <br>
> <a href="http://www.wfu.edu/~natalie" rel="noreferrer" target="_blank">http://www.wfu.edu/~natalie</a><br>
> Wake Forest University Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  phone: <br>
> 1-336-758-5510<br>
> Winston-Salem, NC 27109 USA Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  office: Rm. 300 Olin <br>
> Physical Lab<br>
> <br>
> <br>
> On Wed, Apr 10, 2019 at 6:29 PM Holzwarth, Natalie <<a href="mailto:natalie@wfu.edu" target="_blank">natalie@wfu.edu</a> <br>
> <mailto:<a href="mailto:natalie@wfu.edu" target="_blank">natalie@wfu.edu</a>>> wrote:<br>
> <br>
>  Â  Â Dear Quantum Espresso developers,<br>
>  Â  Â  Â  Â  Â  Â Hopefully this might be something that is easy to track down<br>
>  Â  Â but I cannot figure it out.  Â We noticed that we cannot set the<br>
>  Â  Â initial temperature in our md simulations in QE 6.4.1  (or 6.4) in<br>
>  Â  Â the same way that works fine in QE 6.3.  Â For example, in a simple<br>
>  Â  Â test , the comparison between the results for " grep temperature <br>
>  Â  Â *out"  Â is:<br>
>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  QE  Â 6.3  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  QE 6.4.1<br>
>  Â  Â  Â  Â  Â temperature  Â  =  1800.00000000 K  Â  Â  Â temperature  Â  =  Â <br>
>  Â  Â  Â 0.00000000 K<br>
>  Â  Â  Â  Â  Â temperature  Â  =  1798.79298998 K  Â  Â  Â temperature  Â  =  Â <br>
>  Â  Â  Â 0.00000000 K<br>
>  Â  Â  Â  Â  Â temperature  Â  =  1795.18924725 K  Â  Â  Â temperature  Â  =  Â <br>
>  Â  Â  Â 0.00000001 K<br>
>  Â  Â  Â  Â  Â temperature  Â  =  1789.20477087 K  Â  Â  Â temperature  Â  =  Â <br>
>  Â  Â  Â 0.00000002 K<br>
>  Â  Â  Â  Â  Â temperature  Â  =  1780.85173941 K  Â  Â  Â temperature  Â  =  Â <br>
>  Â  Â  Â 0.00000048 K<br>
>  Â  Â  Â  Â  Â temperature  Â  =  1770.15802096 K  Â  Â  Â temperature  Â  =  Â <br>
>  Â  Â  Â 0.00000014 K<br>
>  Â  Â  Â  Â  Â temperature  Â  =  1757.15885289 K  Â  Â  Â temperature  Â  =  Â <br>
>  Â  Â  Â 0.00000032 K<br>
>  Â  Â  Â  Â  Â temperature  Â  =  1741.89911601 K  Â  Â  Â temperature  Â  =  Â <br>
>  Â  Â  Â 0.00000575 K<br>
>  Â  Â  Â  Â  Â temperature  Â  =  1724.42918183 K  Â  Â  Â temperature  Â  =  Â <br>
>  Â  Â  Â 0.00001204 K<br>
>  Â  Â  Â  Â  Â temperature  Â  =  1704.80836203 K  Â  Â  Â temperature  Â  =  Â <br>
>  Â  Â  Â 0.00001108 K<br>
>  Â  Â  Â  Â  Â temperature  Â  =  1683.10589232 K  Â  Â  Â temperature  Â  =  Â <br>
>  Â  Â  Â 0.00000872 K<br>
>  Â  Â The sample input file is pasted below.  Â The only difference between<br>
>  Â  Â the runs is the version of QE.<br>
>  Â  Â  Â  Thanks in advance for your suggestions.  Â  Sincerely, Natalie<br>
>  Â  Â -------------------sample input file for md<br>
>  Â  Â run-----------------------------<br>
>  Â  Â #!/bin/tcsh<br>
>  Â  Â #<br>
>  Â  Â #SBATCH --nodes=1<br>
>  Â  Â #SBATCH --ntasks-per-node=8<br>
>  Â  Â #SBATCH --cpus-per-task=1<br>
>  Â  Â #SBATCH --account="natalieGrp"<br>
>  Â  Â #SBATCH --output="JOB-%j.o"<br>
>  Â  Â #SBATCH --error="JOB-%j.e"<br>
>  Â  Â #SBATCH --mail-user=<a href="mailto:natalie@wfu.edu" target="_blank">natalie@wfu.edu</a> <mailto:<a href="mailto:natalie@wfu.edu" target="_blank">natalie@wfu.edu</a>><br>
>  Â  Â #SBATCH --mail-type=BEGIN,END,FAIL<br>
>  Â  Â #SBATCH --job-name="NaCl"<br>
>  Â  Â #SBATCH --time=0-10:00:00<br>
>  Â  Â #SBATCH --mem=24gb<br>
>  Â  Â #SBATCH --partition=small<br>
>  Â  Â umask 002<br>
>  Â  Â # Note: SLURM has no batch input for cputime, excluding.<br>
>  Â  Â #<br>
>  Â  Â source /etc/profile.d/modules.csh<br>
>  Â  Â module purge<br>
>  Â  Â module load rhel7/openmpi/3.1.1-intel-2018<br>
>  Â  Â #<br>
>  Â  Â echo 'hostname' `/bin/hostname`<br>
>  Â  Â echo 'job directory' `pwd`<br>
>  Â  Â #<br>
>  Â  Â setenv TMPDIR /scratch/$SLURM_JOBID<br>
>  Â  Â echo 'Reset TMPDIR for this job to ' $TMPDIR<br>
> <br>
>  Â  Â set PW=/home/natalie/EL7/publiccode/QE/qe-6.4.1/bin/pw.x<br>
> <br>
>  Â  Â #NOTE:SLURM defaults to running jobs in the directory where submitted;<br>
>  Â  Â #NOTE:Consider --workdir directive instead; and check functionality!<br>
>  Â  Â cd ${SLURM_SUBMIT_DIR}<br>
> <br>
>  Â  Â cat > NaCl.in << EOF<br>
> <br>
>  Â  Â &CONTROL<br>
>  Â  Â  Â  calculation = "md",<br>
>  Â  Â  Â  pseudo_dir  =<br>
>  Â  Â '/deac/natalieGrp/wfurc9/natalie/EL6/rc9/PAWatoms/poscorenhat/',<br>
>  Â  Â  Â  verbosity  Â = "high",<br>
>  Â  Â  Â  outdir  Â  Â  = "$TMPDIR/",<br>
>  Â  Â  Â  prefix  Â  Â  = "conv",<br>
>  Â  Â  Â  restart_mode = 'from_scratch',<br>
>  Â  Â  Â  nstep = 30,<br>
>  Â  Â  Â  dt = 20,<br>
>  Â  Â  Â  forc_conv_thr = 1.0D-5,<br>
>  Â  Â  Â  etot_conv_thr = 1.0D-6,<br>
>  Â  Â  Â  tstress = .true.,<br>
>  Â  Â  Â  tprnfor = .true.,<br>
>  Â  Â /<br>
>  Â  Â &SYSTEM<br>
>  Â  Â  Â  ibrav  Â  Â  Â = 1,<br>
>  Â  Â  Â  celldm(1)=1.0331018E+01,<br>
>  Â  Â  Â  nat  Â  Â  Â  Â = 8,<br>
>  Â  Â  Â  ntyp  Â  Â  Â  = 2,<br>
>  Â  Â  Â  nosym  Â  Â  Â =.TRUE.,<br>
>  Â  Â  Â  ecutwfc  Â  Â = 64.D0,<br>
>  Â  Â  Â  use_all_frac = .TRUE.,<br>
>  Â  Â  Â  occupations = "smearing",<br>
>  Â  Â  Â  smearing  Â  = "gaussian",<br>
>  Â  Â  Â  degauss  Â  Â = 0.001D0,<br>
>  Â  Â /<br>
>  Â  Â &ELECTRONS<br>
>  Â  Â  Â  conv_thr  Â  = 1.D-8,<br>
>  Â  Â  Â  electron_maxstep = 200,<br>
>  Â  Â /<br>
>  Â  Â &IONS<br>
>  Â  Â  Â ion_dynamics = 'verlet',<br>
>  Â  Â  Â ion_temperature = 'initial',<br>
>  Â  Â  Â tempw = 1800.d0,<br>
>  Â  Â  Â pot_extrapolation = 'second-order',<br>
>  Â  Â  Â wfc_extrapolation = 'second-order',<br>
>  Â  Â /<br>
>  Â  Â &CELL<br>
>  Â  Â  Â  cell_dynamics='none',<br>
>  Â  Â  Â  wmass = 1.0,<br>
>  Â  Â  Â  press = 0.0,<br>
>  Â  Â /<br>
>  Â  Â ATOMIC_SPECIES<br>
>  Â  Â Na  22.989769  Na.LDA-PW-paw.UPF<br>
>  Â  Â Cl  35.45  Â Cl.LDA-PW-paw.UPF<br>
>  Â  Â ATOMIC_POSITIONS {crystal}<br>
>  Â  Â Na  Â  Â  Â 0.000000000  0.000000000  0.000000000<br>
>  Â  Â Na  Â  Â  Â 0.500000000  0.500000000  0.000000000<br>
>  Â  Â Na  Â  Â  Â 0.000000000  0.500000000  0.500000000<br>
>  Â  Â Na  Â  Â  Â 0.500000000  0.000000000  0.500000000<br>
>  Â  Â Cl  Â  Â  Â 0.500000000  0.500000000  0.500000000<br>
>  Â  Â Cl  Â  Â  Â 0.000000000  0.000000000  0.500000000<br>
>  Â  Â Cl  Â  Â  Â 0.500000000  0.000000000  0.000000000<br>
>  Â  Â Cl  Â  Â  Â 0.000000000  0.500000000  0.000000000<br>
>  Â  Â K_POINTS AUTOMATIC<br>
>  Â  Â 2 2 2 1 1 1<br>
>  Â  Â EOF<br>
> <br>
>  Â  Â mpirun $PW  -in  NaCl.in  >  NaCl.out<br>
> <br>
>  Â  Â -------------------------------------------------------------------------------------<br>
> <br>
>  Â  Â N. A. W. Holzwarth Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  email:<br>
>  Â  Â natalie@ <mailto:<a href="mailto:natalie@wfu.edu" target="_blank">natalie@wfu.edu</a>><br>
>  Â  Â Department of Physics Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â web:<br>
>  Â  Â <a href="http://www.wfu.edu/~natalie" rel="noreferrer" target="_blank">http://www.wfu.edu/~natalie</a><br>
>  Â  Â Wake Forest University Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  phone:<br>
>  Â  Â 1-336-758-5510<br>
>  Â  Â Winston-Salem, NC 27109 USA Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  office: Rm. 300 Olin<br>
>  Â  Â Physical Lab<br>
> <br>
> <br>
> _______________________________________________<br>
> Quantum Espresso is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
> <br>
<br>
-- <br>
Lorenzo Paulatto - Paris<br>
_______________________________________________<br>
Quantum Espresso is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail-m_3973581211940267769gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</blockquote></div>