Greeting, QE users!<div>I meet problems when I want to calculate IR spectrum and free energy of isolated CO molecule using fqha.x code.</div><div><br></div><div>At first I run opt with high precision:</div><div><div><div>    &CONTROL</div><div>    calculation = 'relax'</div><div>    forc_conv_thr = 3.8e-5</div><div>    outdir = './CO/'</div><div>    prefix = 'scf'</div><div>    pseudo_dir = './'</div><div>    restart_mode = 'restart'</div><div>    tefield = .FALSE.</div><div>    verbosity = 'high'</div><div>    wf_collect = .TRUE.</div><div>    /</div><div>    &SYSTEM</div><div>    degauss = 0.002</div><div>    ecutrho = 300</div><div>    ecutwfc = 50</div><div>    ibrav = 0</div><div>    nat = 2</div><div>    noncolin = .FALSE.</div><div>    nspin = 1</div><div>    ntyp = 2</div><div>    vdw_corr = 'grimme-d3'</div><div>    /</div><div>    &ELECTRONS</div><div>    conv_thr = 1e-10</div><div>    electron_maxstep = 100</div><div>    mixing_beta = 0.1</div><div>    mixing_mode = 'plain'</div><div>    scf_must_converge = .TRUE.</div><div>    startingwfc = 'file'</div><div>    /</div><div>    &IONS</div><div>    ion_dynamics = 'bfgs'</div><div>    /</div><div>    ATOMIC_SPECIES</div><div>    C 12.0107 C.pbe-n-kjpaw_psl.1.0.0.UPF</div><div>    O 15.9994 O.pbe-n-kjpaw_psl.0.1.UPF</div><div><br></div><div>    K_POINTS {automatic}</div><div>    1 1 1 0 0 0</div><div><br></div><div>    CELL_PARAMETERS {angstrom}</div><div>    20.00000 0.00000 0.00000</div><div>    0.00000 20.00000 0.00000</div><div>    0.00000 0.00000 20.00000</div><div><br></div><div>    ATOMIC_POSITIONS {angstrom}</div><div>    C 10.50000000 10.00000000 10.00000000</div><div>    O 9.50000000 10.00000000 10.00000000</div></div></div><div><br></div><div>then I run ph.x to get dynamic matrix:</div><div><div>    &INPUTPH</div><div>    outdir='./CO/'</div><div>    prefix='scf'</div><div>    tr2_ph=1.0d-12</div><div>    ldisp=.TRUE.</div><div>    nq1=1</div><div>    nq2=1</div><div>    nq3=1</div><div>    epsil=.TRUE.</div><div>    trans=.TRUE.</div><div>    asr=.TRUE.</div><div>    /</div></div><div><br></div><div>I want to get IR spectrum by dynmat.x:</div><div><div>    &input</div><div>    fildyn='./matdyn1'</div><div>    asr='zero-dim'</div><div>    /</div></div><div><br></div><div>,and q2r.x, matdyn.x to obtain phonon DOS:</div><div>input of q2r.x:</div><div><div>    &INPUT</div><div>    fildyn='./matdyn'</div><div>    zasr='zero-dim'</div><div>    flfrc='default.fc'</div><div>    /</div></div><div>input of matdyn.x:</div><div><div>    &INPUT</div><div>    asr='zero-dim'</div><div>    flfrc='default.fc'</div><div>    dos=.TRUE.</div><div>    nk1=1</div><div>    nk2=1</div><div>    nk3=1</div><div>    amass(1) = 12.0107</div><div>    amass(2) = 15.9994</div><div>    /</div></div><div><br></div><div>However, </div><div>1. I obtained no IR information in ouptut of dynmat.x, The output file ends with:"</div><div><div>     ...Force constants read</div><div>     ...epsilon and Z* read</div><div>     Acoustic Sum Rule: || Z*(ASR) - Z*(orig)|| =    1.057617E+00</div><div>     Acoustic Sum Rule: ||dyn(ASR) - dyn(orig)||=    8.730790E-01</div></div><div> ", and I have no idea about it. :(</div><div><br></div><div>2. The DOS value, second column of DOS file (obtained from matdyn.x) is full of 0.000000E+00, although I've got correct vibration frequencies in the output of ph.x. After consideration, I doubt whether it is meaningful to obtain "regular" DOS pattern of isolated molecule, cuz according to my understanding, DOS of isolated molecule is discrete and may contain several delta functions, rather than a continuum function of frequencies of bulk material. Is this the reason why I get full-of-zero phonon DOS? If so, should I smooth DOS manually? </div><div>There's also a relevant question: I've calculated another system in which CO is absorbed. Vibration information of CO is obtained as usual, however, DOS value of frequencies near CO stretching mode frequency are zero, which is similar with what I get in my isolated CO case. DO these two problems have the same reason? Or I made the same mistake when I calculate these two systems?</div><div><br></div><div>Answers will be appreciated. :)</div><div><br></div><div>Best regards!</div><div><br></div><div>Yike Huang, PhD candidate.</div><div>Dalian Institute of Chemical Physics.</div><div><br></div>