<div dir="ltr"><div><span style="font-size:12.8px">Dear </span><span class="" style="font-size:12.8px">Molpro</span><span style="font-size:12.8px"> Users and Developers</span><br></div><div><br></div><div>I have a doubt about the printing of the CI vector coefficients  (CSF and contracted bielectronic functions) using Multireference Rayleigh Schrödinger perturbation theory. In the 2010 version, the above can only be performed with the "RS2" module.</div><div><br></div><div>Reviewing the source code, I found that the CI vector is printed on the output file using the "ciprin" subroutine present in Molpro/src/mrci/ciexit.F, and the amount of information is controlled with "Kref, intern , ksing, kpair" variables.</div><div><br></div><div>In the RS2C module, "ciexit" it's not called, but it can be done by changing a few lines, for example by modifying</div><div><br></div><div>    1202|   c.....calculate properties and print results</div><div>    1203|   c</div><div>    1204|         if(iusecs.ne.0.or.iuseci.ne.0) then</div><div>    1205|           call cirs2exit</div><div>    1206|           goto 9000</div><div>    1207|         else</div><div>    1208|           call gptran(0)</div><div>    1209|           call cigpi(-1)</div><div>    1210|           call cigsi(-1)</div><div>    1211|           call ciexit(nstatt)</div><div>    1202|         end if</div><div><br></div><div>in the cipt2.F (cipro.F) routine, to</div><div><br></div><div>    1202|   c.....calculate properties and print results</div><div>    1203|   c</div><div>    1204|           call ciexit(nstatt)</div><div>    1205|</div><div><br></div><div>the CI vector is printed for systems where "core,0" and "closed,0" in the definition of the reference configuration space. If core or closed are not set to zero, an error occurs and the calculation ends with an error without print the CI vector.</div><div><br></div><div>This is fine for small systems, but it is impractical for medium or large systems.</div><div>My questions are:</div><div>         there is an error when the CI vector is printed in this way?<br></div><div>         Why the calculations end abruptly when core and closed are not equal to 0? </div><div>         Or in their respective case, Is there any other way to print the CI vector with the RS2C module?</div><div><br></div><div>I would greatly appreciate any information you could provide.</div><div><br></div><div>Best regards</div><div>José</div><div>Universidad Nacional</div><div>Autónoma de México</div></div>