Hello, I would like to print the matrices A and B of the Principal Propagator in the RPA aproximation, defined as eq. 3.6 and 3.7 in the work of Oddershede et. al. Computer Physics Reports 2 (1984) 3392.
So far, I printed the whole principal propagator (REDUCED E(2)W*S(2) MATRIX) with the .dal (for example, for a HF molecule), but I need the matrix A and B separately.
**DALTON
.DIRECT
.RUN RESPONSE
**INTEGRALS
.FC
**WAVE FUNCTIONS
.HF
**RESPONSE
.TRPFLG
*LINEAR
.PROPRT
FC F1001
.PROPRT
FC H1002
.PRINT
31
*END OF
Impression of the A and B matrices of the Polarization Propagator in RPA

 Posts: 5
 Joined: 01 Oct 2020, 03:53
 First name(s): Daniel
 Middle name(s): Fernando Esteban
 Last name(s): Bajac
 Affiliation: Universidad Nacional del Nordeste. Corrientes.
 Country: Argentina
Impression of the A and B matrices of the Polarization Propagator in RPA
Last edited by DaniBajac on 01 Oct 2020, 19:47, edited 1 time in total.

 Posts: 395
 Joined: 27 Jun 2013, 18:44
 First name(s): Hans Jørgen
 Middle name(s): Aagaard
 Last name(s): Jensen
 Affiliation: Universith of Southern Denmark
 Country: Denmark
Re: Impression of the A and B values of the Polarization Propagator in RPA
There is a debug keyword, which will calculate the triplet A and B separately in the original orbital basis:
Code: Select all
**RESPONSE
.TRPFLG ! singlet A and B without this keyword
*LINEAR
.ABCHK

 Posts: 5
 Joined: 01 Oct 2020, 03:53
 First name(s): Daniel
 Middle name(s): Fernando Esteban
 Last name(s): Bajac
 Affiliation: Universidad Nacional del Nordeste. Corrientes.
 Country: Argentina
Re: Impression of the A and B values of the Polarization Propagator in RPA
Sorry but I can't print the matrices A and B. Specifically, I only need the twobody components of that matrices: A(1) y B(1).
I used the keyword .ABOCHK instead of .ABCHK, because the last is used with .SINGLE RESIDUE
I used the keyword .ABOCHK instead of .ABCHK, because the last is used with .SINGLE RESIDUE

 Posts: 5
 Joined: 01 Oct 2020, 03:53
 First name(s): Daniel
 Middle name(s): Fernando Esteban
 Last name(s): Bajac
 Affiliation: Universidad Nacional del Nordeste. Corrientes.
 Country: Argentina
Re: Impression of the A and B matrices of the Polarization Propagator in RPA
Hi Mr. Hans Jensen.
It took me a while to understand that the matrices FCOCO and FCOEX represents the Coulomb and Exchange contributions to the Fock matrix.
However, the FCOCO matrix always returns the zero matrix (for the HF and H2O2 molecules). Why could this be happening?
On the other hand, in the file fone.F, I saw that FCOEX = FCOCO + FCOEX. This may indicates that the matrix FCOEX of the output represents both contribution?
Finally, can I consider that the matrices FCOCO and FCOEX are the A and B matrices of the principal propagator?
Best regards
Daniel
It took me a while to understand that the matrices FCOCO and FCOEX represents the Coulomb and Exchange contributions to the Fock matrix.
However, the FCOCO matrix always returns the zero matrix (for the HF and H2O2 molecules). Why could this be happening?
On the other hand, in the file fone.F, I saw that FCOEX = FCOCO + FCOEX. This may indicates that the matrix FCOEX of the output represents both contribution?
Finally, can I consider that the matrices FCOCO and FCOEX are the A and B matrices of the principal propagator?
Best regards
Daniel

 Posts: 395
 Joined: 27 Jun 2013, 18:44
 First name(s): Hans Jørgen
 Middle name(s): Aagaard
 Last name(s): Jensen
 Affiliation: Universith of Southern Denmark
 Country: Denmark
Re: Impression of the A and B matrices of the Polarization Propagator in RPA
Note that the A and B matrices for .SINGLE RESIDUE are the same as for linear response (e.g. for FC).
The FCOCO and FCOEX contain Fock matrix type contributions to the A and B matrices. Note that the program will iteratively solve the linear response equations and not calculate the A and B matrices explicitly. The .ABCHK will calculate the A and B matrices for debugging a column at the time by using trial vectors with one and only one nonzero element (corresponding to the current column in A and B). FCOCO and FCOEX are intermediates for efficiency for trial vectors having many nonzero coefficients.
There is no easy way of getting only the twoelectron contributions to A and B.
The FCOCO and FCOEX contain Fock matrix type contributions to the A and B matrices. Note that the program will iteratively solve the linear response equations and not calculate the A and B matrices explicitly. The .ABCHK will calculate the A and B matrices for debugging a column at the time by using trial vectors with one and only one nonzero element (corresponding to the current column in A and B). FCOCO and FCOEX are intermediates for efficiency for trial vectors having many nonzero coefficients.
There is no easy way of getting only the twoelectron contributions to A and B.

 Posts: 5
 Joined: 01 Oct 2020, 03:53
 First name(s): Daniel
 Middle name(s): Fernando Esteban
 Last name(s): Bajac
 Affiliation: Universidad Nacional del Nordeste. Corrientes.
 Country: Argentina
Re: Impression of the A and B matrices of the Polarization Propagator in RPA
OK Mr. Jensen, I understand. Thank you for the clarification about this themes.
However I'm going to keep looking this matrices because I need the two electron contributions to A and B for my investigation.
Thanks a lot for your answers!
Daniel
However I'm going to keep looking this matrices because I need the two electron contributions to A and B for my investigation.
Thanks a lot for your answers!
Daniel
Who is online
Users browsing this forum: No registered users and 6 guests