Someone else has asked a related question and it hasn't gotten a response. I want to do a frequency calculation in the presence of spacially invariant electric field. So, I am using the multipole field add in for qchem. Here is my input file

$molecule
0 1
C 0.0000000 0.0000000 -0.5602591
N 0.0000000 0.0000000 0.5969240
H 0.0000000 0.0000000 -1.6188328
$end

$multipole_field
X 0.0
Y 0.0
Z 0.0005
$end

$rem
BASIS = 6-31+G*
GUI = 2
JOB_TYPE = freq
METHOD = B3LYP
SCF_CONVERGENCE = 8
$end

When I attempt to run the file, I get the following error:

Q-Chem fatal error occurred in module /home/scratch/svnadmin/11556_zgan/qchem/drvman/scf_hess.C, line 224:

Hessian with external multipole field is not yet implemented

p0_20284: p4_error: interrupt SIGx: 6


The user manual offers no help with this issue and is incredibly terse regarding how to use the multipole implementation. Any help with this would be great! Perhaps someone out there has done this calculation using qchem? Batista at yale has done this exact computation in Gaussian...