<div dir="ltr"><div>Dear users,</div><div>I want to calculate phonon dos for different lattice constants at different temperatures using thermo_pw module of QE. Following is my input:</div><div><br></div><div>&input_thermo<br>what='mur_lc_t',<br>lmurn=TRUE<br>nq1_d=128, nq2_d=128, nq3_d=128<br>ngeo=7<br>step_ngeo=0.5<br>tmin=1<br>tmax=800<br>deltat=3.<br>/<br> &control<br>    calculation = 'scf',<br>    prefix = '${a}'<br>    verbosity= 'high'<br>    tstress= .true.<br>    tprnfor= .true.<br>    outdir = '/home/user/thermo/'<br>    pseudo_dir = '/home/user/thermo/pseudo/'<br> /<br> &system<br>    ibrav =  2,<br>    celldm(1) = $a,<br>    nat =  2,<br>    ntyp = 2,<br>    ecutwfc = 100,<br>/<br>&electrons<br>    mixing_beta = 0.7<br> /<br><br>ATOMIC_SPECIES<br><br>Ca 40.078  Ca.pbe-nsp-van.UPF<br>O 15.999 O.pbe-van_ak.UPF<br><br>ATOMIC_POSITIONS<br>Ca 0.0 0.0 0.0<br>O 0.50 0.50 0.50<br><br>K_POINTS (automatic)<br>  $NK $NK $NK 1 1 1<br>---<br>ciao<br> &inputph<br>  tr2_ph=1.0d-12,<br>  prefix='${a}',<br>  fildyn='${a}.dyn',<br>  ldisp=.TRUE.<br>  fildyn='${a}.dyn',<br>  nq1=4, nq2=4, nq3=4,<br>/<br>&input<br>fildyn='${a}.dyn',<br>zasr='simple',<br>flfrc='${a}.fc'<br>/</div><div><br></div><div>All this blocks are placed in a single input file. But when I try to run, I face the following error:</div><div>./thermo_control: line 1: syntax error near unexpected token `&'<br>./thermo_control: line 1: `&input_thermo'</div><div><br></div><div>Can anyone explain, what is the issue and how do I solve it.<br></div><div>Thanks and regards.<br></div></div>