3
#include <libdpd/dpd.h>
7
/* cc3(): calls libdpd functions to compute the triples contributions
9
* <S| H <T| (Uhat T2)c |0> |T> -> t1
10
* <D| Hhat <T| (Uhat T2)c |0> |T> -> t2
11
* Denominators are subsequently applied in denom()
16
dpdfile2 TIA_new, Tia_new, FME, Fme;
17
dpdbuf4 TIJAB_new, Tijab_new, TIjAb_new;
18
dpdbuf4 TIJAB, Tijab, TIjAb, TiJaB;
19
dpdbuf4 WABEI, WAbEi, Wabei, WaBeI;
20
dpdbuf4 WMBIJ, WMbIj, Wmbij, WmBiJ;
21
dpdbuf4 Dints, DIJAB_anti, Dijab_anti, DIjAb, DiJaB;
22
dpdbuf4 WAMEF, WAmEf, Wamef, WaMeF;
23
dpdbuf4 WMNIE, WMnIe, Wmnie, WmNiE;
25
if(params.ref == 0) { /* RHF */
26
dpd_file2_init(&TIA_new, CC_OEI, 0, 0, 1, "New tIA");
27
dpd_buf4_init(&TIjAb_new, CC_TAMPS, 0, 0, 5, 0, 5, 0, "New tIjAb");
29
dpd_buf4_init(&TIjAb, CC_TAMPS, 0, 0, 5, 0, 5, 0, "tIjAb");
30
dpd_buf4_init(&WAbEi, CC3_HET1, 0, 10, 5, 10, 5, 0, "CC3 WAbEi (iE,bA)");
31
dpd_buf4_init(&WMbIj, CC3_HET1, 0, 0, 10, 0, 10, 0, "CC3 WMbIj (Ij,Mb)");
32
dpd_buf4_init(&Dints, CC_DINTS, 0, 0, 5, 0, 5, 0, "D 2<ij|ab> - <ij|ba>");
33
dpd_file2_init(&FME, CC_OEI, 0, 0, 1, "FME");
34
dpd_buf4_init(&WAmEf, CC3_HET1, 0, 10, 5, 10, 5, 0, "CC3 WAmEf (mA,Ef)");
35
dpd_buf4_init(&WMnIe, CC3_HET1, 0, 0, 10, 0, 10, 0, "CC3 WMnIe (Mn,Ie)");
37
if (params.t3_Ws_incore)
38
cc3_sigma_RHF_ic(&TIjAb, &WAbEi, &WMbIj, 1, &Dints, &TIA_new, 1, &FME, &WAmEf,
39
&WMnIe, &TIjAb_new, moinfo.occpi, moinfo.occ_off, moinfo.virtpi,
40
moinfo.vir_off, 0.0, outfile);
42
cc3_sigma_RHF(&TIjAb, &WAbEi, &WMbIj, 1, &Dints, &TIA_new, 1, &FME, &WAmEf,
43
&WMnIe, &TIjAb_new, moinfo.occpi, moinfo.occ_off, moinfo.virtpi,
44
moinfo.vir_off, 0.0, outfile);
46
dpd_buf4_close(&TIjAb);
47
dpd_buf4_close(&WAbEi);
48
dpd_buf4_close(&WMbIj);
49
dpd_buf4_close(&Dints);
50
dpd_file2_close(&FME);
51
dpd_buf4_close(&WAmEf);
52
dpd_buf4_close(&WMnIe);
54
dpd_file2_close(&TIA_new);
55
dpd_buf4_close(&TIjAb_new);
57
else if(params.ref == 2) { /* UHF */
58
dpd_file2_init(&TIA_new, CC_OEI, 0, 0, 1, "New tIA");
59
dpd_file2_init(&Tia_new, CC_OEI, 0, 2, 3, "New tia");
60
dpd_buf4_init(&TIJAB_new, CC_TAMPS, 0, 0, 5, 2, 7, 0, "New tIJAB");
61
dpd_buf4_init(&Tijab_new, CC_TAMPS, 0, 10, 15, 12, 17, 0, "New tijab");
62
dpd_buf4_init(&TIjAb_new, CC_TAMPS, 0, 22, 28, 22, 28, 0, "New tIjAb");
64
/*** alpha-alpha-alpha */
66
dpd_buf4_init(&TIJAB, CC_TAMPS, 0, 0, 5, 2, 7, 0, "tIJAB");
67
dpd_buf4_init(&WABEI, CC3_HET1, 0, 20, 5, 20, 7, 0, "CC3 WABEI (IE,B>A)");
68
dpd_buf4_init(&WMBIJ, CC3_HET1, 0, 0, 20, 2, 20, 0, "CC3 WMBIJ (I>J,MB)");
69
dpd_buf4_init(&DIJAB_anti, CC_DINTS, 0, 0, 5, 0, 5, 0, "D <IJ||AB>");
70
dpd_file2_init(&FME, CC_OEI, 0, 0, 1, "FME");
71
dpd_buf4_init(&WAMEF, CC3_HET1, 0, 20, 5, 20, 7, 0, "CC3 WAMEF (MA,F>E)");
72
dpd_buf4_init(&WMNIE, CC3_HET1, 0, 0, 20, 2, 20, 0, "CC3 WMNIE (M>N,IE)");
74
cc3_sigma_UHF_AAA(&TIJAB, &WABEI, &WMBIJ, 1, &DIJAB_anti, &TIA_new,
75
1, &FME, &WAMEF, &WMNIE, &TIJAB_new, moinfo.aoccpi, moinfo.aocc_off,
76
moinfo.avirtpi, moinfo.avir_off, 0.0, outfile);
78
dpd_buf4_close(&TIJAB);
79
dpd_buf4_close(&WABEI);
80
dpd_buf4_close(&WMBIJ);
81
dpd_buf4_close(&DIJAB_anti);
82
dpd_file2_close(&FME);
83
dpd_buf4_close(&WAMEF);
84
dpd_buf4_close(&WMNIE);
86
/*** beta-beta-beta */
88
dpd_buf4_init(&Tijab, CC_TAMPS, 0, 10, 15, 12, 17, 0, "tijab");
89
dpd_buf4_init(&Wabei, CC3_HET1, 0, 30, 15, 30, 17, 0, "CC3 Wabei (ie,b>a)");
90
dpd_buf4_init(&Wmbij, CC3_HET1, 0, 10, 30, 12, 30, 0, "CC3 Wmbij (i>j,mb)");
91
dpd_buf4_init(&Dijab_anti, CC_DINTS, 0, 10, 15, 10, 15, 0, "D <ij||ab>");
92
dpd_file2_init(&Fme, CC_OEI, 0, 2, 3, "Fme");
93
dpd_buf4_init(&Wamef, CC3_HET1, 0, 30, 15, 30, 17, 0, "CC3 Wamef (ma,f>e)");
94
dpd_buf4_init(&Wmnie, CC3_HET1, 0, 10, 30, 12, 30, 0, "CC3 Wmnie (m>n,ie)");
96
cc3_sigma_UHF_BBB(&Tijab, &Wabei, &Wmbij, 1, &Dijab_anti, &Tia_new,
97
1, &Fme, &Wamef, &Wmnie, &Tijab_new, moinfo.boccpi, moinfo.bocc_off,
98
moinfo.bvirtpi, moinfo.bvir_off, 0.0, outfile);
100
dpd_buf4_close(&Tijab);
101
dpd_buf4_close(&Wabei);
102
dpd_buf4_close(&Wmbij);
103
dpd_buf4_close(&Dijab_anti);
104
dpd_file2_close(&Fme);
105
dpd_buf4_close(&Wamef);
106
dpd_buf4_close(&Wmnie);
108
/*** alpha-alpha-beta */
110
dpd_buf4_init(&TIJAB, CC_TAMPS, 0, 0, 5, 2, 7, 0, "tIJAB");
111
dpd_buf4_init(&TIjAb, CC_TAMPS, 0, 22, 28, 22, 28, 0, "tIjAb");
112
dpd_buf4_init(&TiJaB, CC_TAMPS, 0, 23, 29, 23, 29, 0, "tiJaB");
114
dpd_buf4_init(&WABEI, CC3_HET1, 0, 20, 5, 20, 7, 0, "CC3 WABEI (IE,B>A)");
115
dpd_buf4_init(&WaBeI, CC3_HET1, 0, 24, 28, 24, 28, 0, "CC3 WaBeI (Ie,Ba)");
116
dpd_buf4_init(&WAbEi, CC3_HET1, 0, 27, 29, 27, 29, 0, "CC3 WAbEi (iE,bA)");
117
dpd_buf4_init(&WMBIJ, CC3_HET1, 0, 0, 20, 2, 20, 0, "CC3 WMBIJ (I>J,MB)");
118
dpd_buf4_init(&WMbIj, CC3_HET1, 0, 22, 24, 22, 24, 0, "CC3 WMbIj (Ij,Mb)");
119
dpd_buf4_init(&WmBiJ, CC3_HET1, 0, 23, 27, 23, 27, 0, "CC3 WmBiJ (iJ,mB)");
120
dpd_buf4_init(&DIJAB_anti, CC_DINTS, 0, 0, 5, 0, 5, 0, "D <IJ||AB>");
121
dpd_buf4_init(&DIjAb, CC_DINTS, 0, 22, 28, 22, 28, 0, "D <Ij|Ab>");
123
dpd_file2_init(&FME, CC_OEI, 0, 0, 1, "FME");
124
dpd_file2_init(&Fme, CC_OEI, 0, 2, 3, "Fme");
125
dpd_buf4_init(&WAMEF, CC3_HET1, 0, 20, 5, 20, 7, 0, "CC3 WAMEF (MA,F>E)");
126
dpd_buf4_init(&WaMeF, CC3_HET1, 0, 24, 28, 24, 28, 0, "CC3 WaMeF (Ma,Fe)");
127
dpd_buf4_init(&WAmEf, CC3_HET1, 0, 27, 29, 27, 29, 0, "CC3 WAmEf (mA,fE)");
128
dpd_buf4_init(&WMNIE, CC3_HET1, 0, 0, 20, 2, 20, 0, "CC3 WMNIE (M>N,IE)");
129
dpd_buf4_init(&WMnIe, CC3_HET1, 0, 22, 24, 22, 24, 0, "CC3 WMnIe (Mn,Ie)");
130
dpd_buf4_init(&WmNiE, CC3_HET1, 0, 23, 27, 23, 27, 0, "CC3 WmNiE (mN,iE)");
132
cc3_sigma_UHF_AAB(&TIJAB, &TIjAb, &TiJaB, &WABEI, &WaBeI, &WAbEi,
133
&WMBIJ, &WMbIj, &WmBiJ, 1, &DIJAB_anti, &DIjAb, &TIA_new, &Tia_new,
134
1, &FME, &Fme, &WAMEF, &WaMeF, &WAmEf, &WMNIE, &WMnIe, &WmNiE,
135
&TIJAB_new, &TIjAb_new, moinfo.aoccpi, moinfo.aocc_off, moinfo.boccpi,
136
moinfo.bocc_off, moinfo.avirtpi, moinfo.avir_off, moinfo.bvirtpi,
137
moinfo.bvir_off, 0.0, outfile);
139
dpd_buf4_close(&TIJAB); dpd_buf4_close(&TIjAb); dpd_buf4_close(&TiJaB);
140
dpd_buf4_close(&WABEI); dpd_buf4_close(&WaBeI); dpd_buf4_close(&WAbEi);
141
dpd_buf4_close(&WMBIJ); dpd_buf4_close(&WMbIj); dpd_buf4_close(&WmBiJ);
142
dpd_buf4_close(&DIJAB_anti); dpd_buf4_close(&DIjAb);
143
dpd_file2_close(&FME); dpd_file2_close(&Fme);
144
dpd_buf4_close(&WAMEF); dpd_buf4_close(&WaMeF); dpd_buf4_close(&WAmEf);
145
dpd_buf4_close(&WMNIE); dpd_buf4_close(&WMnIe); dpd_buf4_close(&WmNiE);
147
/*** beta-beta-alpha term 1 */
149
dpd_buf4_init(&Tijab, CC_TAMPS, 0, 10, 15, 12, 17, 0, "tijab");
150
dpd_buf4_init(&TIjAb, CC_TAMPS, 0, 22, 28, 22, 28, 0, "tIjAb");
151
dpd_buf4_init(&TiJaB, CC_TAMPS, 0, 23, 29, 23, 29, 0, "tiJaB");
153
dpd_buf4_init(&Wabei, CC3_HET1, 0, 30, 15, 30, 17, 0, "CC3 Wabei (ie,b>a)");
154
dpd_buf4_init(&WaBeI, CC3_HET1, 0, 24, 28, 24, 28, 0, "CC3 WaBeI (Ie,Ba)");
155
dpd_buf4_init(&WAbEi, CC3_HET1, 0, 27, 29, 27, 29, 0, "CC3 WAbEi (iE,bA)");
156
dpd_buf4_init(&Wmbij, CC3_HET1, 0, 10, 30, 12, 30, 0, "CC3 Wmbij (i>j,mb)");
157
dpd_buf4_init(&WMbIj, CC3_HET1, 0, 22, 24, 22, 24, 0, "CC3 WMbIj (Ij,Mb)");
158
dpd_buf4_init(&WmBiJ, CC3_HET1, 0, 23, 27, 23, 27, 0, "CC3 WmBiJ (iJ,mB)");
159
dpd_buf4_init(&Dijab_anti, CC_DINTS, 0, 10, 15, 10, 15, 0, "D <ij||ab>");
160
dpd_buf4_init(&DiJaB, CC_DINTS, 0, 23, 29, 23, 29, 0, "D <iJ|aB>");
162
dpd_file2_init(&FME, CC_OEI, 0, 0, 1, "FME");
163
dpd_file2_init(&Fme, CC_OEI, 0, 2, 3, "Fme");
164
dpd_buf4_init(&Wamef, CC3_HET1, 0, 30, 15, 30, 17, 0, "CC3 Wamef (ma,f>e)");
165
dpd_buf4_init(&WaMeF, CC3_HET1, 0, 24, 28, 24, 28, 0, "CC3 WaMeF (Ma,Fe)");
166
dpd_buf4_init(&WAmEf, CC3_HET1, 0, 27, 29, 27, 29, 0, "CC3 WAmEf (mA,fE)");
167
dpd_buf4_init(&Wmnie, CC3_HET1, 0, 10, 30, 12, 30, 0, "CC3 Wmnie (m>n,ie)");
168
dpd_buf4_init(&WMnIe, CC3_HET1, 0, 22, 24, 22, 24, 0, "CC3 WMnIe (Mn,Ie)");
169
dpd_buf4_init(&WmNiE, CC3_HET1, 0, 23, 27, 23, 27, 0, "CC3 WmNiE (mN,iE)");
171
cc3_sigma_UHF_BBA(&Tijab, &TIjAb, &TiJaB, &Wabei, &WaBeI, &WAbEi,
172
&Wmbij, &WMbIj, &WmBiJ, 1, &Dijab_anti, &DiJaB, &TIA_new, &Tia_new,
173
1, &FME, &Fme, &Wamef, &WaMeF, &WAmEf, &Wmnie, &WMnIe, &WmNiE,
174
&Tijab_new, &TIjAb_new, moinfo.aoccpi, moinfo.aocc_off, moinfo.boccpi,
175
moinfo.bocc_off, moinfo.avirtpi, moinfo.avir_off, moinfo.bvirtpi,
176
moinfo.bvir_off, 0.0, outfile);
178
dpd_buf4_close(&Tijab); dpd_buf4_close(&TIjAb); dpd_buf4_close(&TiJaB);
179
dpd_buf4_close(&Wabei); dpd_buf4_close(&WaBeI); dpd_buf4_close(&WAbEi);
180
dpd_buf4_close(&Wmbij); dpd_buf4_close(&WMbIj); dpd_buf4_close(&WmBiJ);
181
dpd_buf4_close(&Dijab_anti); dpd_buf4_close(&DiJaB);
182
dpd_file2_close(&FME); dpd_file2_close(&Fme);
183
dpd_buf4_close(&Wamef); dpd_buf4_close(&WaMeF); dpd_buf4_close(&WAmEf);
184
dpd_buf4_close(&Wmnie); dpd_buf4_close(&WMnIe); dpd_buf4_close(&WmNiE);
187
dpd_file2_init(&t1, CC_OEI, 0, 0, 1, "New tIA");
188
dpd_file2_print(&t1, outfile);
189
dpd_file2_close(&t1);
190
dpd_file2_init(&t1, CC_OEI, 0, 2, 3, "New tia");
191
dpd_file2_print(&t1, outfile);
192
dpd_file2_close(&t1);
193
dpd_buf4_init(&T2, CC_TAMPS, 0, 2, 7, 2, 7, 0, "New tIJAB");
194
dpd_buf4_print(&T2, outfile, 1);
196
dpd_buf4_init(&T2, CC_TAMPS, 0, 12, 17, 12, 17, 0, "New tijab");
197
dpd_buf4_print(&T2, outfile, 1);
199
dpd_buf4_init(&T2, CC_TAMPS, 0, 22, 28, 22, 28, 0, "New tIjAb");
200
dpd_buf4_print(&T2, outfile, 1);