Multiwfn forum / Calculate the one-particle and two-particle density matrix

Re: Calculate the one-particle and two-particle density matrix

dummy@example.com (prasanta13) — Mon, 09 Oct 2023 11:44:08 +0000

Sorry for the late reply, I am incredibly sorry.
I have asked the PSI4 community to extract 2PDM at CCSD(T) level. However, I got nothing yet. Could you tell me how to obtain the 2PDM from the CCSD(T) calculation from PSI4 (preferred since I can access it) or Molpro.

Thanking You
Prasanta

Re: Calculate the one-particle and two-particle density matrix

dummy@example.com (jxzou) — Mon, 11 Sep 2023 04:02:48 +0000

Please note that (un)relaxed density of CCSD(T) in PySCF are problematic currently. The problem occurs at the (T) step, while the CCSD density is correct. See the PySCF github issue https://github.com/pyscf/pyscf/issues/1272.

My personal recommendation is to use Molpro or PSI4. Also note that CCSD(T) density is not supported in ORCA (so far for all 5.x).

Re: Calculate the one-particle and two-particle density matrix

dummy@example.com (prasanta13) — Fri, 11 Aug 2023 12:54:46 +0000

Many many thanks Prof Tian. I am working with PSI4, Lets see if I can find something.

Re: Calculate the one-particle and two-particle density matrix

dummy@example.com (sobereva) — Wed, 09 Aug 2023 13:25:00 +0000

Dear Prasanta,

I don't have much experience in extracting 2RDM. If you only need 1RDM for CCSD(T), I suggest using PSI4, the steps are described in Section 4.A.8 of Multiwfn manual.

Tian

Calculate the one-particle and two-particle density matrix

dummy@example.com (prasanta13) — Wed, 09 Aug 2023 12:07:22 +0000

Hi Prof Tian,
I am trying to obtain the HF and CCSD(T) one and two-particle density matrix from the respective wavefunction (like HF or Correlated). I can do this using PySCF, but it takes a lot of memory (1TB for ~200 basis functions for CCSD(T)). This is why I tried to do it with other programs like Psi4 or ORCA which are less memory hungry. However, I do not know how to obtain the DM from wfn file. Note: I am using a custom basis functions.
Thanks
Prasanta