7
$\begingroup$

When using most molecular electron structure codes that use GTOs for the density or wavefunction, it is possible to project the first-order density matrix onto 'natural orbitals' with MOs expressed using GTOs.

What open-source codes are currently available to fit (presumably in a least-squares sense) MOs expressed using GTOs to:

  1. An arbitrary 3-D gridded density (with either vacuum or 3-D periodic boundary conditions)
  2. Wavefunction/density representations (plane waves, e.g. Quantum Espresso; Augmented plane waves, e.g. EXCITEing, etc.). In the latter case, presumably the 'augmented' part involves GTOs that could be used as at least a starting point for fitting to a 'GTOs only' representation.

Additional background: I want to be able to accurately determine the topology of an arbitrary 3-D distribution of electron density rho(r), together with all the first and second partial derivatives with respect to the spatial coordinates. I also need an accurate quantum stress tensor, also at arbitrary points. Though gridded data introduces a discretization error, e.g. plane waves also introduce such a error as one uses a finite set of g-vectors from the FFT grid. Software I want to use expects GTOs, but if code exists to get the quantities above by other means with minimal spurious effects, then that's OK. And yes, everyone, no matter their chosen representation, is fighting against numerical noise and effects caused by the electronic structure algorithm producing the density. 'Best effort' is what I am going for.

Pointers to code repos gratefully received!

UPDATE: some progress - Frank Jensen below has kindly pointed out his recent work projecting onto GTOs, of which I was previously aware. The author of MultiWFN confirmed that this code can directly use wavefunctions from CP2K, and it seems that 'projwfc' from the Quantum-Espresso suite generates GTOs as an intermediate step to get Lowdin population analysis numbers from a plane-wave wavefunction, though it is currently not clear exactly how to get these GTOs out in Molden/WFX/other useable form. No significant progress on the 'fitting GTOs to an arbitrary 3-D density grid' yet, though.

$\endgroup$
9
  • 1
    $\begingroup$ Least square fitting sounds like a very bad idea… Also, it would be helpful if you can share what you would like to accomplish (not the technical details how you imagine it, but you goal) $\endgroup$
    – Greg
    Apr 3, 2022 at 15:47
  • $\begingroup$ @Greg I want to be able to accurately determine the topology of an arbitrary 3-D distribution of electron density rho(r), together with all the first and second partial derivatives with respect to the spatial coordinates. I also need an accurate quantum stress tensor, also at arbitrary points. Though gridded data introduces a discretization error, e.g. plane waves also introduce such a error as one uses a finite set of g-vectors from the FFT grid. Software I want to use expects GTOs, but if code exists to get the quantities above by other means with minimal spurious effects, then that's OK. $\endgroup$
    – srk
    Apr 4, 2022 at 3:19
  • $\begingroup$ @Greg And yes, everyone, no matter their chosen representation, is fighting against numerical noise and effects caused by the electronic structure algorithm producing the density. 'Best effort' is what I am going for. $\endgroup$
    – srk
    Apr 4, 2022 at 3:26
  • $\begingroup$ maybe you want to put these details in the original question (so people do not need to go through the comment section to see what is the actual question) $\endgroup$
    – Greg
    Apr 4, 2022 at 19:16
  • $\begingroup$ @Greg OK, done. $\endgroup$
    – srk
    Apr 5, 2022 at 0:09

1 Answer 1

4
$\begingroup$

Not clear to me exactly what you want, but we have shown that any valid electron density can be represented by a set of orbitals in a single Slater determinant: P. Jakobsen, F. Jensen "Representing exact electron densities by a single Slater determinant in finite basis sets." J. Chem. Theory Comp. 17 (2021) 269-276. One can think of these as "natural orbitals" but with occupation numbers of exactly 2 (and 0). The code for determining these orbitals is available, but is 'research grade'. We used this for determining a set of orbitals reproducing the full-CI electron density. This is using the reference density expressed in GTOs and the orbitals also expanded in GTOs, where the latter usually needs to be larger than the former for achieving a sufficient accuracy.

$\endgroup$
1
  • $\begingroup$ Hello Frank, yes, that's exactly the paper I linked to in the comments above. $\endgroup$
    – srk
    Apr 6, 2022 at 10:08

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .