I am interested in microiteration calculations with external programs (not with Gaussian). Could you please show us what data and format should be written in xxx_OUT4GEN.rrm. It seems that the manual has no descirptions about the case of TASK: MICROITERATION in xxx_INP4GEN.rrm.
I would appreciate if you could help us.
The file xxx_OUT4GEN.rrm should be written in the same format in every case (TASK: "MAKE GUESS", "MICROITERATION", "ENERGY", "ENERGY and GRADIENT", or "ENERGY, GRADIENT, and HESSIAN").
The required format, for example, for 3-atom system is:
atom1 x1 y1 z1
atom2 x2 y2 z2
atom3 x3 y3 z3
ENERGY = E1 E2 E3
= E4 E5 E6
S**2 = <S2>
DIPOLE = μ(x) μ(y) μ(z)
h(z1,x1) h(z1,y1) h(z1,z1)
h(x2,x1) h(x2,y1) h(x2,z1) h(x2,x2)
h(y2,x1) h(y2,y1) h(y2,z1) h(y2,x2) h(y2,y2)
h(z2,x1) h(z2,y1) h(z2,z1) h(z2,x2) h(z2,y2)
h(x3,x1) h(x3,y1) h(x3,z1) h(x3,x2) h(x3,y2)
h(y3,x1) h(y3,y1) h(y3,z1) h(y3,x2) h(y3,y2)
h(z3,x1) h(z3,y1) h(z3,z1) h(z3,x2) h(z3,y2)
h(y3,z2) h(y3,x3) h(y3,y3)
h(z3,z2) h(z3,x3) h(z3,y3) h(z3,z3)
dμ(x,x1) dμ(y,x1) dμ(z,x1)
dμ(x,y1) dμ(y,y1) dμ(z,y1)
dμ(x,z1) dμ(y,z1) dμ(z,z1)
dμ(x,x2) dμ(y,x2) dμ(z,x2)
dμ(x,y2) dμ(y,y2) dμ(z,y2)
dμ(x,z2) dμ(y,z2) dμ(z,z2)
dμ(x,x3) dμ(y,x3) dμ(z,x3)
dμ(x,y3) dμ(y,y3) dμ(z,y3)
dμ(x,z3) dμ(y,z3) dμ(z,z3)
α(x,z) α(y,z) α(z,z)
If the file xxx_OUT4GEN.rrm is written in different formats, GRRM17 cannot read xxx_OUT4GEN.rrm. Therefore you need to put zero (0.0) into matrix elements of GRADIENT, DIPOLE, HESSIAN, DIPOLE DERIVATIVES, and/or POLARIZABILITY that cannot be obtained from your "sub link = aaa" program. Note: ENERGY requires 6 elements (E1 – E6). In usual cases, the energy from your "sub link = aaa" program have to be put into the first element (E1), and put zero into the others (E2 – E6). The other E2 – E6 elements are reserved for the multistate calculations or some special usages.
In the case "TASK: MICROITERATION", GRRM17 requests you to perform microiterations by your "sub link = aaa" program and to give xxx_OUT4GEN.rrm which contains the coordinates after microiterations, energy, and gradients of the external atoms (although δV/δQn is almost zero for coordinates of external atoms Qn after microiterations). GRRM17 does not perform microiterations, and this is a specification (this might seem to be strange, but this is designed for general use).
For further details, try a microiteration calculation with Gaussian program, and look at the xxx_LinkJOB.rrm file. You will understand how GRRM17 does a microiteration calculation. Even if you use the general interface, GRRM17 requires the same information.
Dear Prof. Saita,
Thank you very much for the info and taking your time. We will give it a try.