Firefly and PC GAMESS-related discussion club



Learn how to ask questions correctly


Re: Once again about XMCQDPT

Alex Granovsky
gran@classic.chem.msu.su


Dear Pasha,

a couple of additional comments in addition to what Sanya has already suggested.

1. Be careful to include all pairs of degenerate states into state-averaging procedure (and use equal weights for each of two states in pair). This is important to keep symmetry of MOs and to keep states
degenerate. Otherwise, results can be senseless. In some cases,
you may need to lower symmetry to a subgroup of the full symmetry group or even to C1 group. Below, I'm assuming you are working in C1 group.

2.  

>   STATE                       1ST ORDER                       2ND ORDER
>     1     E(MCSCF)=    -1727.1284797779     E(MP2)=    -1731.6455874427

>     2     E(MCSCF)=    -1726.9913511748     E(MP2)=    -1731.6292601411
>     3     E(MCSCF)=    -1726.9913511741     E(MP2)=    -1731.5827124526

This is clearly non-sense as degenerate states should preserve their degeneracy in PT treatment. Check if the symmetry of MOs is not  broken, and that you included each component of degenerate pair in PT treatment and in various averaging procedures with equal weights.

Check for correct syntax. Possibly, FF can do some strange things if you specify weights for states not included into PT:

>$xmcqdpt state(1)=1,1,1,1,0,0                                                
> avecoe(1)=1,1,1,1,1,-0
>wstate(1)=1,1,1,1,1,-0  $end  

Is this really the input you have used? I.e., four states in PT and five weights?


>Q1:
>If only 1-n states corresponds to single electron excitation, then we don't need to include another state for XMCQDPT. Or we should consider only state from CASSCF that corresponds to single electron excitation. Am i right?

Perhaps you are right, but could you be a bit more specific?

>Q2:
>What is a difference between
>kstate(1)=1,1,1,1,1,0
>avecoe(1)=1,1,1,0,0,0
>wstate(1)=1,1,1,0,0,0

>and

>kstate(1)=1,1,1,0,0,0
>avecoe(1)=1,1,1,0,0,0
>wstate(1)=1,1,1,0,0,0

>In both cases we have same number of states but weight "0"
>in second corresponding state in kstate array is "0".

This implies that the same averaging procedures are used to obtain semi-canonical orbitals and orbital energies, but different number of states enter PT treatment (i.e., the dimensions of the model space and the effective Hamiltonian differ).


>Q3:
>What should we look in an output file to check XMCQDPT results
>"EIGENVECTORS OF THE EFFECTIVE HAMILTONIAN"
>"EIGENVALUES OF THE NON-SYMMETRIC EFFECTIVE HAMILTONIAN"
>anything else?

Basically, this was already answered here:

http://classic.chem.msu.su/cgi-bin/ceilidh.exe/gran/gamess/forum/?C35b4d3df7Uee-7819-1031-00.htm

All the best,
Alex Granovsky





On Fri Aug 19 '11 6:43pm, Solntsev Pasha wrote
----------------------------------------------
>Dear FF users.

>I have been trying to investigate an excited state of a molecule with the XMCQDPT. From CIS,TDDFT is known that formation of first excited state rises from HOMO->LUMO transition and 2nd excited state from HOMO-1->LUMO. The LUMO orbitals are double degenerated.

>To start my work, i decided consider only HOMO-1, HOMO, 2xLUMO orbitals in my active space. So, i had done  SA-CASSCF(4,4)/6-311G(d,p) with nstate=6
>wstate(1)=1,1,1,1,1,0
>ntrack=5
>and found after MCSCF minimization (only the biggest coefficients are shown)
>1st state corresponds to my ground state
>STATE #    1  ENERGY =   -1727.128479778
>
>
>      CSF      COEF    OCCUPANCY (IGNORING CORE)

>      ---      ----    --------- --------- -----

>        1    0.964072  2200

>.....
>2nd state corresponds to single electron excitation (HOMO->LUMO)
>STATE #    2  ENERGY =   -1726.991351175
>
>
>
>      CSF      COEF    OCCUPANCY (IGNORING CORE)

>      ---      ----    --------- --------- -----

>        1    0.001215  2200

>        2    0.942511  2101

>..............................
>        8    0.261968  1210

>..............................
>       13    0.136427  1012

>.............................

>3rd state also corresponds to single electron excitation (HOMO->LUMO)
> STATE #    3  ENERGY =   -1726.991351174
>
>
>      CSF      COEF    OCCUPANCY (IGNORING CORE)

>      ---      ----    --------- --------- -----

>.................................
>        3    0.261719  1201

>.......................................
>        7   -0.942700  2110

>......................................
>       10    0.129259  2011

>....................................
>       18   -0.137203  1021

>...............................
>Rest of states don't have single electron excitation and corresponds to two electron excitation. The are not presented here.

>It is of interest to note, there is no state corresponding to excitation from HOMO-1->LUMO. Such states exist only as a minor component of another transition (HOMO->LUMO). Is it problem with active space?

>After that i did XMCQDPT and used
>$xmcqdpt state(1)=1,1,1,1,0,0                                                
> avecoe(1)=1,1,1,1,1,-0
>wstate(1)=1,1,1,1,1,-0  $end  

>I got following results:
> *** MC-XQDPT2 ENERGIES ***

> -----------------------------------------------------------------------

>   STATE                       1ST ORDER                       2ND ORDER

>     1     E(MCSCF)=    -1727.1284797779     E(MP2)=    -1731.6455874427

>     2     E(MCSCF)=    -1726.9913511748     E(MP2)=    -1731.6292601411

>     3     E(MCSCF)=    -1726.9913511741     E(MP2)=    -1731.5827124526

>     4     E(MCSCF)=    -1726.9163689912     E(MP2)=    -1730.2704096821

>EIGENVECTORS OF THE EFFECTIVE HAMILTONIAN
>
>
>
>                      1          2          3          4
>
>
>
>               ********************************************
>
>
>
>    1             0.990372   0.120637   0.025590   0.062892

>    2             0.107653  -0.972246   0.185035   0.094396

>    3             0.044631  -0.178233  -0.982213   0.038706

>    4            -0.074714   0.091749   0.019079   0.992792
>
>
>
> EIGENVALUES OF THE NON-SYMMETRIC EFFECTIVE HAMILTONIAN
>
>
>
>     1    -1731.63989383779100000000 + I *        0.00000000000000000000

>     2    -1731.62858390071600000000 + I *        0.00000000000000000000

>     3    -1731.58269032210300000000 + I *        0.00000000000000000000

>     4    -1730.27680165785000000000 + I *        0.00000000000000000000

>As we can see from "EIGENVECTORS OF THE EFFECTIVE HAMILTONIAN
>" part i steel have contribution from 2nd CASSCF state to the 1st XMCQDPT. And the same situation for 2nd and 3rd states for XMCQDPT.

>Experimental value for 1st excited state is 564 nm or 17730 cm^-1 or 0.08078 a.u.

>It looks i very close to finish, but maybe i missed something.

>Also i have couple questions to clarify my understanding.

>Q1:
>If only 1-n states corresponds to single electron excitation, then we don't need to include another state for XMCQDPT. Or we should consider only state from CASSCF that corresponds to single electron excitation. Am i right?

>Q2:
>What is a difference between
>kstate(1)=1,1,1,1,1,0
>avecoe(1)=1,1,1,0,0,0
>wstate(1)=1,1,1,0,0,0

>and

>kstate(1)=1,1,1,0,0,0
>avecoe(1)=1,1,1,0,0,0
>wstate(1)=1,1,1,0,0,0

>In both cases we have same number of states but weight "0"
>in second corresponding state in kstate array is "0".

>Q3:
>What should we look in an output file to check XMCQDPT results
>"EIGENVECTORS OF THE EFFECTIVE HAMILTONIAN"
>"EIGENVALUES OF THE NON-SYMMETRIC EFFECTIVE HAMILTONIAN"
>anything else?

>I will be very appreciate for you help.
>If you need extra information about my current problem just let me know.

>Best wishes.

>Pavel.
>
>
>
>


[ Previous ] [ Next ] [ Index ]           Sat Aug 20 '11 6:18pm
[ Reply ] [ Edit ] [ Delete ]           This message read 840 times