CC2 Energy: wrong MP2 values leading to failure

If you think you found a bug this is where you can report it
Post Reply
Jean Q.
Posts: 7
Joined: 18 Aug 2016, 13:55
First name(s): Jean
Last name(s): Quertinmont
Affiliation: UNamur
Country: Belgium

CC2 Energy: wrong MP2 values leading to failure

Post by Jean Q. » 02 May 2017, 10:42

Hi,

I'm trying to run 'simple' CC2 energy calculations on 'smallish' systems (~100 electrons), starting with a small basis set (cc-pVDZ) but the direct coupled cluster energy program's initial iteration gives a incorrect MP2 value: -63693.1 a.u. vs -630.1 a.u. obtained with a separate MP2 calculation with the same basis set. The following CC2 iterations energies are out of format and the calculation stops after the following message: "ZERO PIVOT - MATRIX SINGULAR".

I'm using Dalton2016.1 compiled with intel 17.0.1, 64-bit integers, and MKL lib. I have the same issue with intel 14.0.2.

Thanks in advance,

Jean.
Attachments
Emp2_anil-dz.txt
(98.12 KiB) Downloaded 7 times
Ecc2_anil-dz.txt
(91.69 KiB) Downloaded 8 times

Jean Q.
Posts: 7
Joined: 18 Aug 2016, 13:55
First name(s): Jean
Last name(s): Quertinmont
Affiliation: UNamur
Country: Belgium

Re: CC2 Energy: wrong MP2 values leading to failure

Post by Jean Q. » 05 May 2017, 11:49

More information:
if I use the Choleski decomposition, the MP2 energy is correct and the calculation finishes normally.
Attachments
Ecc2_chol_anil-dz.txt
(101 KiB) Downloaded 7 times

sonco
Posts: 86
Joined: 06 Sep 2013, 13:49
First name(s): Sonia
Last name(s): Coriani
Affiliation: DTU Chemistry
Country: Denmark

Re: CC2 Energy: wrong MP2 values leading to failure

Post by sonco » 05 May 2017, 15:16

Hello,
I can try and have a look by the end of next week.
Sonia

Jean Q.
Posts: 7
Joined: 18 Aug 2016, 13:55
First name(s): Jean
Last name(s): Quertinmont
Affiliation: UNamur
Country: Belgium

Re: CC2 Energy: wrong MP2 values leading to failure

Post by Jean Q. » 08 May 2017, 10:36

Hello,
Out of curiosity I tried to compute excitation energies with "RUN PROPERTIES" and surprisingly (to me at least), the CC2 energy converges properly and the calculation finishes normally.
Jean.
Attachments
EXCcc2_prop_anil-dz.txt
(106.09 KiB) Downloaded 8 times

kennethruud
Posts: 228
Joined: 27 Aug 2013, 16:42
First name(s): Kenneth
Last name(s): Ruud
Affiliation: UiT The Arctic University of Norway
Country: Norway

Re: CC2 Energy: wrong MP2 values leading to failure

Post by kennethruud » 08 May 2017, 12:21

Hmmm, somewhat surprised myself. Could you try a Hartree-Fock calculation and see if the end result is the same. I am somewhat concerned it continues without an error/warning, using a Hartree-Fock wave function and density. If this is the case, then we need to insert a stop so that this does not happen.

Thanks and best regards,

Kenneth

Jean Q.
Posts: 7
Joined: 18 Aug 2016, 13:55
First name(s): Jean
Last name(s): Quertinmont
Affiliation: UNamur
Country: Belgium

Re: CC2 Energy: wrong MP2 values leading to failure

Post by Jean Q. » 08 May 2017, 13:44

Looks like you're right, the excitation energies (and oscillator strengths) are the same for the CC2 calculation and HF calculation:
exc HF (eV) CC2 (eV)
1 4.709882 4.709882
2 5.334553 5.334552
3 5.701641 5.701640
4 5.799008 5.799008
5 6.130453 6.130453

(sorry for the bad formatting, I couldn't figure out how to make a table or tabulations)
Attachments
EXChf_prop_anil-dz.txt
(97.56 KiB) Downloaded 7 times

Jean Q.
Posts: 7
Joined: 18 Aug 2016, 13:55
First name(s): Jean
Last name(s): Quertinmont
Affiliation: UNamur
Country: Belgium

Re: CC2 Energy: wrong MP2 values leading to failure

Post by Jean Q. » 09 May 2017, 09:00

Since there seems to be an issue with the properties module, I also tried the response one. As expected, the calculation failed at the CC2 energy step, just like for the energy calculation (wrong MP2 then failure to converge).

I also have an excitation energies calculation with CC2 and Choleski running (input taken from the 2016 manual, p. 178). I'm not sure to fully understand that input, but it looks like it uses a specific module for CC linear response.
Attachments
EXCcc2_resp_anil-dz.txt
(92.03 KiB) Downloaded 6 times

Jean Q.
Posts: 7
Joined: 18 Aug 2016, 13:55
First name(s): Jean
Last name(s): Quertinmont
Affiliation: UNamur
Country: Belgium

Re: CC2 Energy: wrong MP2 values leading to failure

Post by Jean Q. » 10 May 2017, 09:31

CC2 Choleski excitation energies are different from the HF ones:
exc HF (eV) CC2/Chol (eV)
1 4.709882 4.05682
2 5.334553 4.55965
3 5.701641 5.05287
4 5.799008 5.17112
5 6.130453 5.89983
Attachments
EXCcc2_chol_anil-dz.txt
(440.2 KiB) Downloaded 6 times

Jean Q.
Posts: 7
Joined: 18 Aug 2016, 13:55
First name(s): Jean
Last name(s): Quertinmont
Affiliation: UNamur
Country: Belgium

Re: CC2 Energy: wrong MP2 values leading to failure

Post by Jean Q. » 21 Jun 2017, 16:09

Any news?

Post Reply

Who is online

Users browsing this forum: No registered users and 2 guests