NPT production
 &cntrl
  imin = 1,	ntx = 1,	maxcyc=1,
  ntpr = 1,	ntwr = 1,	ntwx = 1,
  ntf = 1,	ntc = 1,
  ntb = 2,	cut = 12.0,
  temp0 = 293.0,	ntt = 3,	gamma_ln = 2,
  ntp = 1,	pres0 = 1.0,	taup = 1.0,
  ifqnt=1,
 /
&qmmm
  qmmask=':1',
  qmcharge= 0 ,
  qm_theory='EXTERN',
  qmshake=0,
  qm_ewald=0,
  writepdb=1,
/
&gau
method='blyp',
basis='6-31G(d)',
num_threads=8,
/

