<div dir="ltr">The OPTION,ITRDM=1 option for the CI program is recognized by the RS2 program, but doesn't seem to do anything when run.<div><br></div><div>The RS2 program outputs the following</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">1PROGRAM * RS2 (Multireference RS Perturbation Theory)     Authors: H.-J. Werner (1993), P. Celani (1998)</blockquote>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"> </blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
            XMS-CASPT2 extension: T. Shiozaki (2011)</blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
 </blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">  Level shift= 0.20</blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
 </blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">    Transition moment calculation</blockquote>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">    ============================</blockquote></blockquote>
<div><br></div><div>And that's it. By this point, I'm fairly certain that you have to mess with the source code if you want to get MS-CASPT2 transition moments. </div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Thu, Jul 31, 2014 at 4:59 PM, Alan Chien <span dir="ltr"><<a href="mailto:alandc@umich.edu" target="_blank">alandc@umich.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div><div>Unfortunately, I have become a little stumped.</div><div><br></div><div>I've found that density matrices for various states can be saved by way of the DM card within RS2. </div></div><div>And I can access the particular density set I need using the DENSITY card from section 4.11</div>

<div><br></div><div>But I can't determine a card that will force Molpro to calculate the transition moments. The closest I've found is the TRANS card in the CI program, but that requires wavefunctions stored on records as the inputs, whereas I only have the density matrices. </div>

<div><br></div><div>I will keep trawling the manual and update any progress.</div><div>Any help in the meantime would be appreciated.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>~Alan Chien</div></font></span></div>
<div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><br><div class="gmail_quote">
On Thu, Jul 31, 2014 at 10:17 AM, Alan Chien <span dir="ltr"><<a href="mailto:alandc@umich.edu" target="_blank">alandc@umich.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div dir="ltr">Is it possible to calculate transition moments in a state-averaged CASPT2 calculation?<div><br></div><div>The manual states in section 21.1 "<span style="font-family:NimbusRomNo9L;font-size:11pt">the calculation of CASPT2/RASPT2 density matrices (and therefore molecu-
lar properties) is presently possible only with the RS2 command..."</span></div><div><span style="font-family:NimbusRomNo9L;font-size:11pt"><br></span></div><div>And so I would think that transition moments could be calculated (with rs2), but there aren't any explicit instructions in chapter 21. Careful reading seems to hint that a lot of the MRCI commands will work in CAPT2 calcs however. </div>


<div><br></div><div>Thus, I'm currently trying the commands that are explicitly mentioned for one-electron properties in the CI chapter, chapter 20. Specifically sections 20.3.13 and 20.3.14 (EXPEC and TRANS). Although now I wonder where CASPT2 saves its wavefunctions for use in TRANS. Or if I have to do that manually with a SAVE command. </div>


<div><br></div><div>Anyways, any help/confirmation would be nice so that I don't fumble around in the dark too much.</div><div><br></div><div>Thanks,</div><div>ALan CHien</div><div><span style="font-family:NimbusRomNo9L;font-size:11pt"><br>


</span></div><div><span style="font-family:NimbusRomNo9L;font-size:11pt">ps - All section numbers are from Molpro User Manual 2012.1 </span></div></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>