Dear all,

I want to calculate the derivative of the dipole moment using Q-Chem. Here is my simple input file:$molecule
$molecule
-1 2
O 0.0082665743052699 0.0264660201009203 -0.0169816216498090
H 0.9627580947076764 0.0231938084260272 0.0104758031541364
H -0.2418348863845040 0.9649523667325521 0.0440012990273091
$end

$external_charges
20.8970398275127636 20.6456862861180745 3.9279618930789715 -0.8340000000000000
21.8105842737005311 20.8356287532084465 4.2531574499543972 0.4170000000000000
20.4010288046450441 21.0858762747596700 4.6549778097773498 0.4170000000000000
$end

$rem
JOBTYPE freq
EXCHANGE wb97x
BASIS 6-31g*
$end
I can only get the derivative of the dipole moment for QM atoms in Q-Chem output file. I want to know if there is a way to get the derivative of dipole moment for the MM atoms. Looking forward to your relpy. Thanks !