# Molecular orbitals used with CCSD(T) geometries

Suppose I am given some geometry data (say, of a water dimer) from CCSD(T).
If I were to do a single-point energy calculation to generate orbitals for the system, wouldn't CCSD give something more representative of the correlation in the system than DFT?

I see "CCSD(T) orbitals don't exist - used DFT orbitals instead" kind of statements in papers (quotations below), and wasn't sure why they don't even mention CCSD. Perhaps it's just due to speed/availability, but that would lead to another question "why is it not worth it to use CCSD orbitals with CCSD(T) geometry".

The same question can be asked for DFT orbitals being used with MP2 geometries, but MP2 perturbs about HF ... so clearly DFT is needed to capture correlation.

[edit: more context]
What I'm trying to do is "export" these orbitals to a diffusion quantum Monte Carlo (DMC) calculation which doesn't do geometry optimization by itself but augments the imported orbitals with explicit correlations (eg. something that depends directly on electron pair distances). Usually the final DMC energy is not very sensitive to whether the orbitals are from CC or DFT (I have outputs from both), but I was checking here to make sure I wasn't missing any "deeper" reason.

Since MP2 orbitals are nonexistent, the charge-transfer interactions were computed using the B3LYP functional with the aug'-cc-pVTZ basis set at the MP2/aug'-cc-pVTZ geometries, so that at least some electron correlation effects could be included.

Due to the nonexistence of CCSD(T) orbitals, ΔE(del) was calculated with ωB97XD/aug-cc-pVTZ.

• I'm guessing that a lot of times, this comes down to the necessary methods being very costly and/or not implemented. The methods that then consume these orbitals are considered to be not even semi-quantitatively correct by many, so DFT orbitals are probably sufficient most of the time. Mar 11 '19 at 17:02
• I thought I had a straightforward answer to this, but I believe the difficulty is associated with the fact both these methods are perturbative in nature. In principle, one should be able to generate a new set of orbitals, but I believe in practice the orbitals from perturbative QM methods are actually somewhat arbitrary as the energy is only minimized with respect to excitation amplitudes out of the reference determinant which does not guarantee the orbitals are stationary. Thus, properties in these methods are computed via linear response rather than as expectation values over the orbitals. Mar 11 '19 at 17:45
• The real question is "what do you mean by orbitals?", which is partially related to "what do you want to look at them for?" Certainly all of these methods have natural orbitals which can be obtained by diagonalizing the corresponding 1-particle density matrix (which will have non-integer occupation numbers), and also certainly there are reference orbitals used to compute these methods. Those cited statements are, to some extent, meaningless without the answers to those questions. Mar 11 '19 at 23:13
• And to clarify, if you run a single-point, you have generated an energy, not orbitals. All orbitals are an implementational tool to obtain this energy and, in that sense, no orbitals "exist", in that they are not observable and can be replaced with an alternative method to obtain the same energy. So in order to obtain orbitals you need to clarify what you mean by that. Mar 11 '19 at 23:23
• I also don't know what you mean by using "X orbitals with Y geometries". All of these calculations are done with the Born-Oppenheimer approximation. You fix the nuclei, then solve for the electronic energy. How you obtained the geometry has nothing to do with what method you subsequently use to compute energies or properties. Mar 11 '19 at 23:23

The orbitals used for the coupled clusters variants that you mentioned are (almost always) those from the Hartree-Fock method, which is performed as a first step. Note that the latter does not account for correlation energy. The case of MP$$n$$ is exactly the same.