How Much Computing Power is Required for Full CI Calculation of a Molecule?

AI Thread Summary
To determine the computing power required for a Full CI calculation of a molecule with n nuclei and x electrons, one must consider the basis set size, k, which defines the number of spin orbitals as 2k. The number of Slater determinants formed is approximately x^(2k), leading to a Hamiltonian matrix of size m x m that needs diagonalization. This matrix grows factorially with the system size, making exact calculations impractical for more than one electron. For example, using a finite basis set for a simple molecule like H2O (n=3, x=10) illustrates the significant computational demands, as processing time increases dramatically with the number of electrons. Thus, while Full CI offers exact solutions, the required resources often render it infeasible for larger systems.
Big-Daddy
Messages
333
Reaction score
1
For a molecule of n nuclei and x total electrons, how can I work out how much computing power or processing time is required for an exact Full CI calculation for that molecule?
 
Chemistry news on Phys.org
Big-Daddy said:
For a molecule of n nuclei and x total electrons, how can I work out how much computing power or processing time is required for an exact Full CI calculation for that molecule?

For an exact full CI you need a complete MO basis which contains an infinity of molecular orbitals and you can form an infinity of Slater determinants from it. Hence you will need infinite computing power for more than one electron in general.
So you have to be less ambitious. E.g. you can chose a certain finite basis set, e.g. VQZ. Let's say you have 100 basis functions per atom (quite a big basis already) (this will also depend on the type of atom). If you have 10 atoms with 10 electrons each, you can form 100 out of 1000 Slater determinants (roughly 10^300) so this will still take almost forever.
 
Hang on, please clarify. Ok so as you've pointed out an exact solution will take infinite time, i.e. impossible, except for one-electron systems. If you have a certain basis set with k basis functions per nucleus, and n nuclei with x total electrons ... I still don't get how you are proposing to work out the amount of computing power needed? Thanks for the guidance. Can you do an example, maybe, let there be k basis functions, but now we're talking about an H2O molecule (n=3, x=10), how much processing power/time is needed?
 
Big-Daddy said:
Hang on, please clarify. Ok so as you've pointed out an exact solution will take infinite time, i.e. impossible, except for one-electron systems. If you have a certain basis set with k basis functions per nucleus, and n nuclei with x total electrons ... I still don't get how you are proposing to work out the amount of computing power needed? Thanks for the guidance. Can you do an example, maybe, let there be k basis functions, but now we're talking about an H2O molecule (n=3, x=10), how much processing power/time is needed?

So if you have k basis functions and x electrons, you have 2k spin orbitals. As each Slater determinant has to be made up from x different spin orbitals, you can form $$m=\begin{pmatrix} 2k\\x \end{pmatrix}\approx x^{2k}$$ different Slater determinants (you could classify them by spin and symmetry to reduce the number, but we are interested in the gross picture). So basically your Hamiltonian becomes an $$m \times m$$ matrix which you have to diagonalize. The matrix will be sparse as the two electron matrix elements vanish between most of the Slater determinants but still the size increases factorially with the size of the system.
 
Ok thanks.
So what value of k would we typically use for something very simple like a He molecule or H- ion?

And with the same k, let's say processing time on a certain computer is 10 seconds for a H- ion (2 electrons); will it then be 10!/2!*10 seconds for a (10-electron) water molecule?
 
It seems like a simple enough question: what is the solubility of epsom salt in water at 20°C? A graph or table showing how it varies with temperature would be a bonus. But upon searching the internet I have been unable to determine this with confidence. Wikipedia gives the value of 113g/100ml. But other sources disagree and I can't find a definitive source for the information. I even asked chatgpt but it couldn't be sure either. I thought, naively, that this would be easy to look up without...
I was introduced to the Octet Rule recently and make me wonder, why does 8 valence electrons or a full p orbital always make an element inert? What is so special with a full p orbital? Like take Calcium for an example, its outer orbital is filled but its only the s orbital thats filled so its still reactive not so much as the Alkaline metals but still pretty reactive. Can someone explain it to me? Thanks!!
Back
Top