<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">Hi,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I’m trying to compute a dipole moment using the finite-field approach, for a variety of correlated methods. This is working, but during testing I found results that appear buggy for the correlated density. There are only two non-zero
components (Cs point group). Here are the values I get. The molecule is (SO2)(CO), the calc is CCSD/aug-cc-pV(D+t)Z. <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Frozen-core<o:p></o:p></p>
<p class="MsoNormal">Code #1: mu_x = 0.3663 D, mu_z = 1.8332 D (E = -660.852312)<o:p></o:p></p>
<p class="MsoNormal">Code #2: mu_x = 0.3653 D, mu_z = 1.8334 D (E = -660.852312)<o:p></o:p></p>
<p class="MsoNormal">MOLPRO: mu_x = 0.4450 D, mu_z = 1.7298 D (E = -660.852312)<o:p></o:p></p>
<p class="MsoNormal">MOLPRO, finite-field: mu_x = 0.3656(6) D, mu_z = 1.8325(3) D<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Active core<o:p></o:p></p>
<p class="MsoNormal">Code #1: mu_x = 0.3672 D, mu_z = 1.8395 D (E = -660.916466)<o:p></o:p></p>
<p class="MsoNormal">Code #2: mu_x = 0.3662 D, mu_z = 1.8397 D (E = -660.916466)<o:p></o:p></p>
<p class="MsoNormal">MOLPRO: mu_x = 0.4455 D, mu_z = 1.7370 D (E = -660.916466)<o:p></o:p></p>
<p class="MsoNormal">MOLPRO, finite-field: mu_x = 0.3665(6) D, mu_z = 1.8388(3) D<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Everything is nicely consistent at the HF level. <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Here is my input file:<o:p></o:p></p>
<p class="MsoNormal">~~~~~~~~~~~~~~~~~~~~~~~~~~~<o:p></o:p></p>
<p class="MsoNormal">***, SO2 CO Complex<o:p></o:p></p>
<p class="MsoNormal">memory,220,m<o:p></o:p></p>
<p class="MsoNormal">angstrom<o:p></o:p></p>
<p class="MsoNormal">geomtyp=xyz<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">geometry={<o:p></o:p></p>
<p class="MsoNormal"> 5<o:p></o:p></p>
<p class="MsoNormal">UCCSD(T)/AUG-CC-PV(D+D)Z ENERGY=-660.89251446<o:p></o:p></p>
<p class="MsoNormal">S -1.2013598699 0.0000000000 0.3645914371<o:p></o:p></p>
<p class="MsoNormal">O -1.2280384976 1.2649632874 -0.3826423804<o:p></o:p></p>
<p class="MsoNormal">O -1.2280384976 -1.2649632874 -0.3826423804<o:p></o:p></p>
<p class="MsoNormal">C 2.1633891600 0.0000000000 0.1896882719<o:p></o:p></p>
<p class="MsoNormal">O 3.2695357051 0.0000000000 -0.1104089492<o:p></o:p></p>
<p class="MsoNormal">}<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">basis=AUG-CC-PV(D+d)Z<o:p></o:p></p>
<p class="MsoNormal">hf<o:p></o:p></p>
<p class="MsoNormal">{ccsd,thrden=1.0d-11,thrvar=1.0d-11;<o:p></o:p></p>
<p class="MsoNormal">expec,dm<o:p></o:p></p>
<p class="MsoNormal">}<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">show,dm*<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">deb(1)=dmx*TODEBYE<o:p></o:p></p>
<p class="MsoNormal">deb(2)=dmy*TODEBYE<o:p></o:p></p>
<p class="MsoNormal">deb(3)=dmz*TODEBYE<o:p></o:p></p>
<p class="MsoNormal">show,deb<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">~~~~~~~~~~~~~~~~~~~~~~~~~~~<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Please let me know if I am doing something wrong here. <o:p>
</o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks,<o:p></o:p></p>
<p class="MsoNormal">Karl<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif"">-------------------------------------------------------------------</span><o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif"">Dr. Karl K. Irikura</span><o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif"">National Institute of Standards and Technology</span><o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif"">100 Bureau Drive, mail stop 8320</span><o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif"">Gaithersburg, MD 20899-8320</span><o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif"">voice: 301-975-2510 fax: 301-975-3670</span><o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif"">email:
<a href="mailto:karl.irikura@nist.gov"><span style="color:blue">karl.irikura@nist.gov</span></a></span><o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif""><a href="http://www.nist.gov/mml/csd/informatics_research/index.cfm"><span style="color:blue">http://www.nist.gov/mml/csd/informatics_research/index.cfm</span></a>
</span><o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif"">-------------------------------------------------------------------</span><o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>