(fwd) [molpro-user] Natural Orbitals for UHF, MP2, and QCISD

H. -J. Werner werner at theochem.uni-stuttgart.de
Wed Mar 31 07:37:12 BST 2004


Dear Jun,
here is the answer:

1.) Using the current program (2002.6), you need to set
    option,cphf=1 for mp2 or qcisd:

    mp2;natorb,9000.2,5;option,cphf=1

    I just fixed this so that the cphf option is set automatically
    if a natorb directive is present. There will be a patch 
    2002.6.natorb.sh.gz available later today. For more details
    about the natorb directive, please see manual.

2.) There is a bug in generating the qcisd natural orbitals. This is
    also fixed in patch 2002.6.natorb.sh.gz

3.) To print uhf natural orbitals, use either

    gprint,orbitals=1
    uhf

or

    uhf
    print,1

Note that UHF natural orbitals are not printed in closed-shell cases (i.e., if
they are identical to the spin-orbitals). Also note that uhf cannot be
used as reference for mp2, qcisd, ccsd or any other correlation method.

H.-J. Werner
----- Forwarded message from "Li, Jun" <Jun.Li at pnl.gov> -----

From: "Li, Jun" <Jun.Li at pnl.gov>
To: molpro-user at molpro.net
Subject: [molpro-user] Natural Orbitals for UHF, MP2, and QCISD
Date: Tue, 30 Mar 2004 20:00:52 -0800
Message-id: <7931A5BD58A7014F93A09EFC25866C650154AE02 at pnlmse27.pnl.gov>

Dear MOLPRO community,

I am interested in printing the natural orbitals (NO) from calculations
like UHF, MP2, and QCISD.  I played with several options, including the
following two, but cannot figure it out.  If any of you have experience
with this, could you please advise?  

dm,9000.2;natorb,9000.2,print=5;
matrop;load,den,den,5000.2;natorb,orbs,den; 

Also I was trying to modify the UNO-CAS criteria (unomax=1.99;
unomin=0.01), but it seems there is no effect with these options.  Any
of your help would be greatly appreciated!

Best regards,

Jun Li

 

----- End forwarded message -----

-- 
Prof. Hans-Joachim Werner
Institute for Theoretical Chemistry
University of Stuttgart
Pfaffenwaldring 55
D-70569 Stuttgart, Germany
Tel.: (0049) 711 / 685 4400
Fax.: (0049) 711 / 685 4442
e-mail: werner at theochem.uni-stuttgart.de



More information about the Molpro-user mailing list