<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">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_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>