Ask for temeparture command

Questions
#1

This is my input file for mixed calculation. I want to add temperature in this file. How to add temperature in qchem input file?

$molecule
0 1
P 0.40335 3.63431 3.68942
P 0.09807 1.34448 5.98624
O -0.04719 4.93666 4.69593
O -0.22359 3.86081 2.09024
H -1.17045 3.84348 2.12701
C -0.27358 1.99522 4.27180
C 0.30161 1.00306 3.27115
H 1.34898 0.88662 3.46676
H 0.15176 1.38231 2.29129
C -0.40282 -0.34837 3.37776
H -0.25167 -0.77636 4.35028
H -1.44686 -0.20323 3.21598
C 0.17751 -1.27760 2.29768
H 1.21598 -1.44522 2.49851
H 0.06763 -0.80556 1.34404
O -1.69760 2.06223 4.18062
H -1.93938 2.34756 3.30377
N -0.54301 -2.55472 2.29743
H -0.16229 -3.14257 1.58220
H -0.45332 -3.00224 3.18278
O 1.78931 1.34497 6.31106
H 1.93654 0.93105 7.17422
Zn 1.46050 5.24005 5.77394
Zn 0.44746 3.30313 8.37828
O -0.69822 2.19365 7.28774
C 3.93278 5.87932 9.31020
C 3.96562 6.54253 10.71217
H 4.55012 5.96347 11.41278
H 2.99405 6.65350 11.09504
C 4.63811 7.90351 10.44323
H 4.04651 8.42041 9.74671
H 5.61843 7.70424 10.05333
C 4.82849 8.81436 11.66034
H 3.88791 9.06222 12.12924
H 5.46007 8.35847 12.39095
N 5.47101 10.03883 11.12305
H 5.64693 10.69324 11.85933
H 6.35042 9.77002 10.67881
O 5.27995 5.91731 8.84605
H 5.57382 6.83421 8.86781
O 2.96037 5.54190 6.84632
O 1.51928 4.43460 9.52371
P 3.20456 4.24273 9.17935
P 2.81628 6.63789 8.15489
O 3.95321 3.04218 10.15512
H 4.88028 2.96009 9.88489
O 2.85139 5.86737 4.62504
H 2.54467 6.63808 4.14179
O 0.36945 2.17219 9.91548
H -0.50236 2.23145 10.31301
O 3.07272 8.07178 7.77726
O 3.34695 3.84461 7.73518
O -0.44297 -0.05969 5.99044
O 1.90357 3.59726 3.57815
O 1.25334 6.57525 8.84582
H 1.00943 7.44672 9.16620
$end

$rem
basis = gen
GUI = 2
JOBTYPE = opt
MEM_TOTAL = 42000
exchange=PBE
correlation=PBE
ECP=gen
ECP_FIT = TRUE
$end

$basis
C 0
6-31G(d,p)

H 0
6-31G(d,p)

N 0
6-31G(d,p)

O 0
6-31G(d,p)

P 0
6-31G(d,p)

Zn 0
LANL2DZ

$end

$ecp
Zn 0
LANL2DZ

$end

@@@@@@@
$molecule
read
$end

$rem
basis = gen
GUI = 2
JOBTYPE = freq
MEM_TOTAL = 42000
exchange=PBE
correlation=PBE
ECP=gen
ECP_FIT = TRUE
$end

$basis
C 0
6-31G(d,p)

H 0
6-31G(d,p)

N 0
6-31G(d,p)

O 0
6-31G(d,p)

P 0
6-31G(d,p)

Zn 0
LANL2DZ

$end

$ecp
Zn 0
LANL2DZ

$end

#2

Temperature only makes sense if you are doing (a) ab initio MD, or (b) vibrational frequencies, in which case it affects the harmonic partition function for the vibrational entropy and free energy G = H - TS.

For AIMD, see
https://manual.q-chem.com/latest/Ch9.S9.SS1.html
For vibrational frequencies, see
https://manual.q-chem.com/latest/Ch10.S7.SS1.html

In your vibrational frequency calculation, see this section to run it at finite temperature:
https://manual.q-chem.com/latest/ex_Isotope-Mass.html