<div dir="ltr">Dear Aurelie,<div><br></div><div style>I am not sure about the error you get but take a look on a recent paper by Tatiana Korona (<a href="http://www.tandfonline.com/doi/suppl/10.1080/00268976.2011.638330#tabModule">http://www.tandfonline.com/doi/suppl/10.1080/00268976.2011.638330#tabModule</a>). In the supplementary material she explains how a vector can be added to a set of HF orbitals. While her usage is different than yours it should be easy for you to modify it to suit your needs.</div>
<div style><br></div><div style>I hope this helps.</div><div style><br></div><div style>Best,</div><div style>Shachar </div><div style><br></div><div style><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Tue, Aug 13, 2013 at 10:03 PM, aurelie perveaux <span dir="ltr"><<a href="mailto:aurelie.perveaux@u-psud.fr" target="_blank">aurelie.perveaux@u-psud.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
Dear All,</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
I want to add an orbital in a set of orbital already defined. I want to
add this atomic orbital to be consider as a molecular orbital, i don t
want it to mix with the other orbital that is why i m not using a user
defined basis set to put it in.</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
So i thought about using the move/add option but i am encouting some
trouble.</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
here my input:</p>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
memory,100,m</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
print,basis,orbitals,civector</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
symmetry,x</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
geometry={</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
O1</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
O2 , O1 , R1</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
O3 , O1 , R2 , O2 , A</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
}</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
R1=2.409</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
R2=2.409</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
A=116.9 Degree</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
spherical</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
basis=aug-cc-pvtz</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
hf</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
orbital,2140.2</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
{mcscf</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
start,2140.2</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
closed,3,0</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
occ,12,3</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
wf,24,2,0;state,2</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
wf,24,1,0;state,4</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
orbital,2141.2</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
}</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
symmetry,x</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
geometry={</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
O</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
}</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
spherical</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
basis={</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
spdf,O,vtz</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
s,1,0.0000001</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
}</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
hf</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
orbital,2142.2</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
symmetry,x</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
geometry={</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
O1</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
O2 , O1 , R1</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
O3 , O1 , R2 , O2 , A</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
}</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
R1=2.409</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
R2=2.409</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
A=116.9 Degree</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
{merge</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
orbital,2141.2</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
move,1.1,1.1,1.1</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
move,2.1,2.1,2.1</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
move,3.1,3.1,3.1</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
move,4.1,4.1,4.1</div>
<div style="font-size:13px;font-family:helvetica,arial,sans-serif">
<div>
move,5.1,5.1,5.1</div>
<div>
move,6.1,6.1,6.1</div>
<div>
move,7.1,7.1,7.1</div>
<div>
move,8.1,8.1,8.1</div>
<div>
move,9.1,9.1,9.1</div>
<div>
move,10.1,10.1,10.1</div>
<div>
move,11.1,11.1,11.1</div>
<div>
move,12.1,12.1,12.1</div>
<div>
move,1.2,1.2,1.2</div>
<div>
move,2.2,2.2,2.2</div>
<div>
move,3.2,3.2,3.2</div>
<div>
orbital,2142.2</div>
<div>
add,5.1,5.1,13.1</div>
<div>
save,2143.2}</div>
<div>
</div>
<div>
</div>
<div>
I don't want to touch to my first molecular orbitals that is why i m
redefining them., and after with the add card i m adding my extra
orbital save in the record 2142.2</div>
<div>
unfortunalty that generat this error message:</div>
<div>
<div>
</div>
<div>
PROGRAM * MERGE Author: H.-J. Werner</div>
<div>
</div>
<div>
</div>
<div>
Molecular orbitals read from record 2141.2
Type=MCSCF/NATURAL (state averaged)</div>
<div>
</div>
<div>
</div>
<div>
Error in MERGE/ADD/MOVE: Basis dimension or offset not correct:
1 90 61</div>
<div>
</div>
</div>
</div>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
==> i do understand that the problem must be the size of the basis
which is diferent in the output and the input (one orbital more in the
output than in the input). But i don t understand who to use the offset
(if the offset is usefull to fix this problem) or to define properly the
basis size.</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
Maybe i m not even doing the right strategy for what i want to achieve.
Feel free to let me know!</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
Thank you for your answers.</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
Best Regards</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
Aurelie Perveaux</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
------------------------------</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
Laboratoire de Chimie Physique d'Orsay, Université de
Paris-Sud</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
Bâtiment349 - Campus d'Orsay</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
15, avenue Jean PERRIN</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
91405 Orsay Cedex</p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
Tél : <a href="tel:%2B33%20%280%29%201%2069%2015%2075%2075" value="+33169157575" target="_blank">+33 (0) 1 69 15 75 75</a></p>
<p style="font-size:13px;font-family:helvetica,arial,sans-serif">
Fax : <a href="tel:%2B33%20%280%29%201%2069%2015%2061%2088" value="+33169156188" target="_blank">+33 (0) 1 69 15 61 88</a></p>
<br>_______________________________________________<br>
Molpro-user mailing list<br>
<a href="mailto:Molpro-user@molpro.net">Molpro-user@molpro.net</a><br>
<a href="http://www.molpro.net/mailman/listinfo/molpro-user" target="_blank">http://www.molpro.net/mailman/listinfo/molpro-user</a><br></blockquote></div><br></div>