Slow SCF convergence?
See original GitHub issueIt looks like SCF convergence is slower in the current version of PySCF than in previous releases.
For instance, the example calculation on water
from pyscf import gto, scf
mol = gto.Mole()
mol.verbose = 4
mol.atom = '''
O 0 0 0
H 0 -2.757 2.587
H 0 2.757 2.587'''
mol.basis = 'ccpvdz'
mol.symmetry = 1
mol.build()
for guess in ['1e', 'atom', 'minao', 'huckel']:
mf = scf.RHF(mol)
mf.init_guess = guess
mf.kernel()
takes
- 29 iterations to converge from the
core
guess - 27 iterations to converge from the
atom
guess - 17 iterations to converge from the
minao
guess - 26 iterations to converge from the
huckel
guess
Water in cc-pVDZ should converge in 6-8 iterations. What’s happened?!
Issue Analytics
- State:
- Created 2 years ago
- Reactions:1
- Comments:6 (3 by maintainers)
Top Results From Across the Web
SCF Convergence Issues - ORCA Input Library - Google Sites
If the SCF is oscillating wildly in the first iterations or just converging slowly, the grid (DFT or COSX) can sometimes be the...
Read more >Convergence in SCF Calculation - Mopac
(b) The equations are only very slowly convergent. This can be due to a long-lived oscillation or to a slow transfer of charge....
Read more >4.5 Converging SCF Calculations - Q-Chem Manual
SCF calculations for systems with zero or small HOMO-LUMO gap (such as metals) can exhibit very slow convergence or may even fail to...
Read more >Convergence in SCF calculation
The equations are only very slowly convergent. This can be due to a long-lived oscillation or to a slow transfer of charge. Oscillation...
Read more >SCF Convergence Guidelines for ADF
Finally, many SCF convergence problems are rooted in non-physical ... be used as a starting point for a slow but steady SCF iteration...
Read more >
Top Related Medium Post
No results found
Top Related StackOverflow Question
No results found
Troubleshoot Live Code
Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start Free
Top Related Reddit Thread
No results found
Top Related Hackernoon Post
No results found
Top Related Tweet
No results found
Top Related Dev.to Post
No results found
Top Related Hashnode Post
No results found
The same geometry is used in many other examples, as well.
dft/02-gks.py gto/10-atom_info.py gto/11-basis_info.py gto/12-serialization.py scf/01-h2o.py scf/02-ghf.py scf/02-rohf_uhf.py scf/03-level_shift.py
I’m not sure I’d call it H2O when the OH distance is 3.78 Å(!)
The odd thing is that the calculations in #1000 don’t work anymore, even though they should work…