CI without Hartree-Fock SCF

Usually textbooks on CI takes the SCF Hartree-Fock orbitals as a starting point. Can one use the configuration interaction method starting directly from the single-particle Slater determinants avoiding the iterative procedure. This question concerns mainly small molecules like the hydrogen molecule. Do I get a high inaccuracy in this case? If yes, why? Could this inaccuracy be decreased taking into account more configurations?

• ummmm… why would you want to do that? :-) I guess I'm wondering if you are doing singles and doubles or toward full CI? I think you would get toward the Full CI limit faster (with fewer excitations) if you start with more reasonable orbitals. Commented May 25, 2013 at 22:08
• Thank you for the reply. That is actually what I am asking. How faster do I get convergence in the case of SCF orbitals. How many additional excitations I have to take into account in the case of avoiding preliminary HF computations. Commented May 25, 2013 at 22:19
• @freude: I don't think it's really possible to make an accurate but general statement about it: it'll vary a lot by system, and will be less important in simple systems where the difference between cheap approximations and accurate models is smaller. Commented May 26, 2013 at 9:16
• Szabo and Ostlund describes CI calculations using natural orbitals where the convergence is faster than with Hartree-Fock orbitals. You should read that book for more information. Commented Oct 20, 2013 at 21:45
• "How faster do I get convergence in the case of SCF orbitals?" Start with a better initial guess. Commented Oct 23, 2013 at 20:57

If you were to perform Full CI (within a finite basis) then I don't think that it matters whether you start with SCF orbitals or some other choice. This solution will have the lowest conceivable energy no matter what.

But if you are doing truncated (less than full) CI, then you are very much at the mercy of your starting orbitals. In general, one would want to start with orbitals that resemble those of the actual state that you are trying to describe via CI.

Otherwise, your precious excitations (singles and doubles) may be consumed correcting the "slop" that you provided, rather than moving toward the full CI limit.

• (As an example, I know some people using CI-related approaches (CASSCF and CASPT2, specifically) on heavier elements who prefer to start with DFT orbitals because they produce a better result than starting with HF.) Commented May 25, 2013 at 22:46
• Right. I had a little trouble answering this question, because practically speaking, there is always optimization of the inital vectors. It's sort of a "textbook" question. Very interesting about DFT orbitals for metals. I've never done MCSCF on metals--they probably need all the help they can get :-) Commented May 25, 2013 at 23:01
• Isn't there a theory that the exact natural orbitals (and approximations to them) should give you the quickest convergence, though? I believe that's the idea behind Iterative Natural Orbital approaches, for example. Commented May 26, 2013 at 0:02
• @EricBrown: Natural orbitals are the eigenfunctions of the 1st order reduced density matrix of a wavefunction. If you take n strongly orthogonal wavefunctions and mix them, your leading NOs will not make up the Slater determinant with the largest prefactor. Of course that's contrived, since the wavefunctions also are orthogonal in the energy product. Commented May 30, 2013 at 16:05
• If by "some other choice" you mean a basis that can be obtained from the other one by unitary rotation then of course you are right. Otherwise it is not true that the Full CI will give you the same result. Energy minimization by Full CI only depends on your basis. Commented Oct 21, 2013 at 15:16