Hi Francois,
Please see the output files attached.
I tested the same input file using Qbox 1.63.8 and 1.60.4.
It seems for the PBE0 calculation, the old version of Qbox converges much better than the new one. Why?
For the PBE calculation, they are the same.
Thank you!
Best wishes,
Ding
PBE0 converges differently between 2 Qbox versions
Forum rules
You must be a registered user to post in this forum. Registered users may also post new topics if they consider that their subject does not correspond to any topic already present on the forum.
You must be a registered user to post in this forum. Registered users may also post new topics if they consider that their subject does not correspond to any topic already present on the forum.
-
- Posts: 8
- Joined: Mon Mar 23, 2015 9:48 pm
PBE0 converges differently between 2 Qbox versions
- Attachments
-
- result.old.r
- qbox 1.60.4
- (187.73 KiB) Downloaded 1727 times
-
- result.r
- qbox 1.63.8
- (264.56 KiB) Downloaded 1695 times
-
- Site Admin
- Posts: 167
- Joined: Tue Jun 17, 2008 7:03 pm
Re: PBE0 converges differently between 2 Qbox versions
Hi Ding,
The default tolerance used to control the iterative diagonalization algorithm in the bisection calculation has been modified in release 1.60.9 (i.e. between these two versions). The value used in 1.60.4 was 1.e-6, and the value used in 1.63.8 is 1.0, which is likely too large. This can be modified using the following command:
The large default value (1.0) causes the bisection algorithm to stop too early in its iterations, which results in slight changes in the truncation domains of orbitals between successive iterations. This causes small fluctuations in the energy, and it prevents the calculation from reaching full convergence. Using the value 1.e-6 fixes the problem and full convergence can be reached. The default value will be decreased in the next release.
Francois
The default tolerance used to control the iterative diagonalization algorithm in the bisection calculation has been modified in release 1.60.9 (i.e. between these two versions). The value used in 1.60.4 was 1.e-6, and the value used in 1.63.8 is 1.0, which is likely too large. This can be modified using the following command:
Code: Select all
set debug BISECTION_TOL 1.e-6
Francois