<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Dear Hyungjun,
<div class=""><br class="">
</div>
<div class="">since the last program run before mrcc was multi, where the occ and closed cards do not correspond to a single determinant, mrcc got a little confused. Before the call to mrcc, simply insert:</div>
<div class=""><br class="">
</div>
<div class="">
<div class="">{rhf;</div>
<div class=""> occ,4,1,1,0,1,0,0;</div>
<div class="">closed,2,1,1,0,1,0,0;</div>
<div class=""> wf,12,1,2;</div>
<div class=""> start,2140.2;maxit,1}</div>
<div class=""><br class="">
</div>
<div class="">Note that I've use the multi orbitals as initial guess and set the maximum number of iterations to 1. This will result in orbitals unchanged from your symmetry equivalenced ones from multi. You can then leave out all the options to mrcc except
method=ccsdt.</div>
<div class=""><br class="">
</div>
<div class="">regards,</div>
<div class=""><br class="">
</div>
<div class="">-Kirk</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div>
<blockquote type="cite" class="">
<div class="">On Nov 13, 2014, at 1:28 AM, Hyungjun Kim <<a href="mailto:hyungjun96@kaist.ac.kr" class="">hyungjun96@kaist.ac.kr</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class=""><style type="text/css" class="">
p{margin:0;padding:0;}
</style>
<div style="font-size: 11pt; font-family: 돋움; background-color: transparent;" class="">
<p class="">Dear Peterson, </p>
<p class=""><br class="">
</p>
<p class=""><span style="font-size: 11pt; line-height: 1.2;" class="">These are input and output files.</span></p>
<p class=""><br class="">
</p>
<p class="">Actually the input file for Zr atom is made by the simple modification on the input file you sent to me for Hf atom.</p>
<p class=""><br class="">
</p>
<p class="">Yours sincerely,</p>
<p class=""><br class="">
</p>
<p class="">-Hyungjun</p>
<br class="">
<table width="99%" align="center" class="">
<tbody class="">
<tr valign="top" class="">
<td class="">
<p class=""><span style="font-family: Arial;" class="">---------------------------------------------------</span><br class="">
<span style="font-family: Arial;" class=""><b class="">Hyungjun Kim</b>, Ph. D. Student,</span><br class="">
<span style="font-family: Arial;" class="">Quantum Chemistry Laboratory, Department of Chemistry,</span></p>
<p class=""><span style="font-family: Arial;" class="">Korea Advanced Institute of Science and Technology [KAIST],</span><br class="">
<span style="font-family: Arial;" class="">Daejeon 305-701, Korea</span><br class="">
<span style="font-family: Arial;" class="">e-mail: </span><a href="mailto:hyungjun96@kaist.ac.kr" class=""><span style="font-family: Arial;" class="">hyungjun96@kaist.ac.kr</span></a><span style="font-family: Arial;" class=""> or
<a href="mailto:hyungjun96@gmail.com" class="">hyungjun96@gmail.com</a></span></p>
<span style="font-family: Arial;" class=""></span>
<p class=""><span style="font-family: Arial;" class="">phone : +82-42-350-2861</span></p>
<span style="font-family: Arial;" class=""></span>
<p class=""><span style="font-family: Arial;" class="">mobile phone : +82-10-8537-5051</span></p>
<span style="font-family: Arial;" class=""></span>
<p class=""><span style="font-family: Arial;" class="">web : <a href="http://qclab.kaist.ac.kr" class="">
http://qclab.kaist.ac.kr</a></span><br class="">
</p>
</td>
</tr>
</tbody>
</table>
<p class=""><br class="">
</p>
<blockquote style="border-left: #000 2px solid;padding:0 0 0 5px;margin-left:5px;" class="">
<div style="line-height:18px;font-size:9pt" class="">-----------------------원본메일-----------------------<br class="">
보낸사람 : "Peterson, Kirk" <<a href="mailto:kipeters@wsu.edu" class="">kipeters@wsu.edu</a>><br class="">
받는사람: Hyungjun Kim <<a href="mailto:hyungjun96@kaist.ac.kr" class="">hyungjun96@kaist.ac.kr</a>><br class="">
참 조: "<a href="mailto:molpro-user@molpro.net" class="">molpro-user@molpro.net</a>" <<a href="mailto:molpro-user@molpro.net" class="">molpro-user@molpro.net</a>><br class="">
보낸날짜: 2014-11-12 21:43:09 GMT +0900 (Asia/Seoul)<br class="">
제목: Re: [molpro-user] Number of electrons in MRCC</div>
<br class="">
<br class="">
<div class="">Dear Hyungjiun,</div>
<div class=""><br class="">
</div>
<div class="">Can you post your actual input? There are some tricks to using orbitals from multi in mrcc. </div>
<div class=""><br class="">
</div>
<div class="">Regards,</div>
<div class=""><br class="">
</div>
<div class="">Kirk</div>
<div class=""><br class="">
Sent from my iPhone</div>
<div class=""><br class="">
On Nov 12, 2014, at 6:54 AM, Hyungjun Kim <<a href="mailto:hyungjun96@kaist.ac.kr" class="">hyungjun96@kaist.ac.kr</a>> wrote:<br class="">
<br class="">
</div>
<blockquote type="cite" class="">
<div class="">
<div style="font-size: 11pt; font-family: 돋움; background-color: transparent;" class="">
<p class=""><span style="font-family: Arial; font-size: 11pt;" class="">Dear Molpro developers,</span></p>
<p class=""><br class="">
</p>
<p class=""><span style="font-family: Arial; font-size: 11pt;" class="">I have reported the same problem to the MRCC forum and the developer of MRCC, Mihaly Kallay, suggests me to report this error to Molpro developer.</span></p>
<p class=""><br class="">
</p>
<p class=""><span style="font-family: Arial; font-size: 11pt;" class="">This is my post on the MRCC forum.</span></p>
<p class=""><br class="">
</p>
<p class=""><span style="font-family: Arial; font-size: 11pt;" class="">***</span></p>
<p class=""><span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">I have tried to calculate the electronic structure of Zr atom using CCSDT in MRCC program interfaced to Molpro. Stuttgart ECP is
employed during this calculation. (Therefore only 12 electrons are left.)</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">I have symmetrized HF orbitals and have done CCSDT calculation using these orbitals.</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">While 4s and 4p orbitals are left as core electrons, 2 electrons in 5s, and 2 electrons in 4d orbitals are distributed, respectively.</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">When MRCC reads CASSCF MO, there are 4 electrons.</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">Frozen orbitals: 4 ( 1 1 1 0 1 0 0 0)</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">Active orbitals: 6 ( 3 0 0 1 0 1 1 0)</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">Active electrons: 4</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">Spin quantum number: 1.0</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">However, the output created by MRCC gives there are only 2 electrons (1alpha and 1beta).</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">Number of alpha electrons: 1</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">Number of beta electrons: 1</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">I don't know why the number of electrons is reduced from 4 to 2. Could you give me a guide to include '4d' electrons in MRCC calculation?</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">Any comments are welcome.</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">Yours sincerely,</span></p>
<p class=""><span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class="">***</span></p>
<p class=""><span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 11pt; line-height: 17.3290004730225px;" class=""><br class="">
</span></p>
<p class=""><font color="#333333" face="Arial" class=""><span style="line-height: 17.3290004730225px;" class="">The reply of Kallay to this question is as follows:</span></font></p>
<p class=""><font color="#333333" face="Arial" class=""><span style="line-height: 17.3290004730225px;" class="">***</span></font></p>
<p class=""><span style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 11pt; line-height: 17.3290004730225px;" class="">Dear Hyungjun,</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 11pt; line-height: 17.3290004730225px;" class="">The occupation vector and the specification of active orbitals passed over to mrcc is not correct. You should contact
the molpro developers about this issue. The quick solution is to specify the directory where mrcc is executed, that is, the corresponding line in the molpro input should be extended, e.g., as</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 11pt; line-height: 17.3290004730225px;" class="">{MRCC,method=CCSDT,spatial=1,HF=0,nacto=0,dir='/scr/mrcc';</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 11pt; line-height: 17.3290004730225px;" class="">Then the scratch files required by mrcc will be saved to that directory. To correct the error edit file fort.56, and
replace the last two line by the following ones:</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 11pt; line-height: 17.3290004730225px;" class="">2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 11pt; line-height: 17.3290004730225px;" class="">1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 1 0 0 0</span><br style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 12.7272720336914px; line-height: 17.3290004730225px;" class="">
<span style="color: rgb(51, 51, 51); font-family: Arial, Helvetica, sans-serif; font-size: 11pt; line-height: 17.3290004730225px;" class="">Then execute dmrcc in that directory. Note that in this case an MRCCSDT calculation will be run. If you want to perform
a simple CCSDT calculation, the last line of fort.56 must be deleted.</span></p>
<p class=""><font color="#333333" face="Arial, Helvetica, sans-serif" size="2" class=""><span style="line-height: 17.3290004730225px; font-size: 11pt;" class="">***</span></font></p>
<p class=""><font color="#333333" face="Arial, Helvetica, sans-serif" size="2" class=""><span style="line-height: 17.3290004730225px; font-size: 11pt;" class=""><br class="">
</span></font></p>
<p class=""><font color="#333333" face="Arial, Helvetica, sans-serif" size="2" class=""><span style="line-height: 17.3290004730225px; font-size: 11pt;" class="">I have fixed the problem referring to Kallay's comment.</span></font></p>
<p class=""><font color="#333333" face="Arial, Helvetica, sans-serif" size="2" class=""><span style="line-height: 17.3290004730225px; font-size: 11pt;" class=""><br class="">
</span></font></p>
<p class=""><font color="#333333" face="Arial, Helvetica, sans-serif" class=""><span style="line-height: 17.3290004730225px;" class="">Could you please check this problem?</span></font></p>
<p class=""><font color="#333333" face="Arial, Helvetica, sans-serif" class=""><span style="line-height: 17.3290004730225px;" class=""><br class="">
</span></font></p>
<p class=""><font color="#333333" face="Arial, Helvetica, sans-serif" class=""><span style="line-height: 17.3290004730225px;" class="">Yours sincerely,</span></font></p>
<p class=""><font color="#333333" face="Arial, Helvetica, sans-serif" class=""><span style="line-height: 17.3290004730225px;" class=""><br class="">
</span></font></p>
<br class="">
<table width="99%" align="center" class="">
<tbody class="">
<tr valign="top" class="">
<td class="">
<p class=""><span style="font-family: Arial;" class="">---------------------------------------------------</span><br class="">
<span style="font-family: Arial;" class=""><b class="">Hyungjun Kim</b>, Ph. D. Student,</span><br class="">
<span style="font-family: Arial;" class="">Quantum Chemistry Laboratory, Department of Chemistry,</span></p>
<p class=""><span style="font-family: Arial;" class="">Korea Advanced Institute of Science and Technology [KAIST],</span><br class="">
<span style="font-family: Arial;" class="">Daejeon 305-701, Korea</span><br class="">
<span style="font-family: Arial;" class="">e-mail: </span><a href="mailto:hyungjun96@kaist.ac.kr" class=""><span style="font-family: Arial;" class="">hyungjun96@kaist.ac.kr</span></a><span style="font-family: Arial;" class=""> or
<a href="mailto:hyungjun96@gmail.com" class="">hyungjun96@gmail.com</a></span></p>
<p class=""><span style="font-family: Arial;" class="">phone : +82-42-350-2861</span></p>
<p class=""><span style="font-family: Arial;" class="">mobile phone : +82-10-8537-5051</span></p>
<p class=""><span style="font-family: Arial;" class="">web : <a href="http://qclab.kaist.ac.kr/" class="">
http://qclab.kaist.ac.kr</a></span><br class="">
</p>
</td>
</tr>
</tbody>
</table>
<p class=""><br class="">
</p>
</div>
<img id="mailexp" width="0" heigh="0" border="0" src="https://mail.kaist.ac.kr/Mail?act=RECEIPT_CHECK&ukey=5462ea753fd530fc6d23a290&userid=hyungjun96&mhost=kaist.ac.kr&ahost=d0001" class=""></div>
</blockquote>
<blockquote type="cite" class="">
<div class=""><span class="">_______________________________________________</span><br class="">
<span class="">Molpro-user mailing list</span><br class="">
<span class=""><a href="mailto:Molpro-user@molpro.net" class="">Molpro-user@molpro.net</a></span><br class="">
<span class=""><a href="http://www.molpro.net/mailman/listinfo/molpro-user" class="">http://www.molpro.net/mailman/listinfo/molpro-user</a></span></div>
</blockquote>
</blockquote>
<p class=""><br class="">
</p>
</div>
<img id="mailexp" width="0" heigh="0" border="0" src="https://mail.kaist.ac.kr/Mail?act=RECEIPT_CHECK&ukey=5463fcb73feaef2567042a44&userid=hyungjun96&mhost=kaist.ac.kr&ahost=d0001" class=""><span id="cid:07CBACA8-2307-4F54-AF4B-984B3B3F4312"><zr-atom.inp></span><span id="cid:BFA446D5-B007-4C02-9E8B-9378246C0F78"><zr-atom.out></span></div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>