<div dir="ltr">Dear Professor Knizia<br><br>Thank you very much for your reply. I have a question concerning the coefficients of <br>the doubly external configurations.<br><br>My aim is obtain the CI vector in terms of determinants. I know this can not be done <br>directly in the mrci module of Molpro, so I wrote a perl script.<br><br>I have no problems in the case of the reference function and simple external configurations. <br>I can go from CSF to determinants by the genealogical coupling scheme. However, this <br>is more difficult in the doubly external configurations case.<br><br>First, i transform each contracted configuration in terms of uncontracted CSF´s<br><br>        \Psi_{ijp}^{ab}=\frac{1}{2}(\hat{E}_{ai,bj}+p\hat{E}_{bi,aj})\Psi_{0}<br>                                 =\sum_{P\nu}<\Psi_{ijp}^{ab}|\Psi_{P\nu}^{ab}>\Psi_{P\nu}^{ab}<br>                       <br>where \Psi_{0} is the reference function (p=1 singlet coupling, p=-1 triplet coupling).<br>Then, I multiply each uncontracted CSF by the associated coefficient with \Psi_{ijp}^{ab} <br>(listed in the output).  Later I express each uncontracted CSF in terms of determinants. <br>As I understand it, the labels i, j, a, b, p in  \Psi_{ijp}^{ab} refer to <br>"I    J -> K    L     NP" in the output of Molpro.<br><br>However, performing this procedure does not achieve a square norm that matches the output<br><br>           CLASS          SQ.NORM        ECORR1        ECORR2<br>           +++++++++++++++++++++++++++++++++++++++++++++++++++<br>           Internals    0.00001717    0.00000000   -0.00000069<br>           Singles      0.00065315   -0.00263386   -0.00263924<br>           Pairs          <b> <u>0.00787253</u></b>   -0.03783725   -0.03783119<br><br>I do not understand that I'm doing wrong.<br>Perhaps, the coefficients associated with \Psi_{ijp}^{ab} are refer to a non-orthonormal basis?<br>I would appreciate any help that you could provide me.<br><br><br>regards    <br><br>Jose Jara<br>Universidad Nacional<br>Autónoma de México</div>