Multiwfn forum

Multiwfn official website: http://sobereva.com/multiwfn. Multiwfn forum in Chinese: http://bbs.keinsci.com/wfn. E-mail of admin: sobereva[at]sina.com

You are not logged in.

#1 2022-12-09 09:44:23

MhamadHantro
Member
Registered: 2022-12-09
Posts: 2

Difference in the transition density between Gaussian and Multiwfn.

Hello,

I wanted to calculate transition dipole moments using the NTOs, so I used Gaussian to output me a formchk containing all the NTOs, then I can construct the transition density and do a simple integral. But, this gave me a wrong answer. After checking my code, and not finding anything, I decided to use Mutltiwfn and feed it the original formchk and the log file from Gaussian, so that it gives me back the NTOs, and of course I took the IOp into account.
Surprisingly, the NTOs constructed by the formchk outputted by Multiwfn gave the right answer that I was waiting, and yet Gaussian didn't. Now, this is already enough for me, as I can continue my work.

Though, I'm just curious on how does Gaussian outputs its NTOs in the formchk? Because to my knowledge, it should just replace the MO coeffcients by the NTO coefficients. Yet, apparently, this is not the case.

Thank you.

Last edited by MhamadHantro (2022-12-09 09:54:25)

Offline

#2 2022-12-10 05:18:05

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 1,611
Website

Re: Difference in the transition density between Gaussian and Multiwfn.

By default, even if you added pop=NTO in Gaussian input file, the chk/fch file produced by CIS/TDDFT task still records MOs of reference state (usually ground state), which are not NTOs.

Offline

#3 2022-12-12 08:37:08

MhamadHantro
Member
Registered: 2022-12-09
Posts: 2

Re: Difference in the transition density between Gaussian and Multiwfn.

Note that if you mean by default, not adding the SaveNTO keyword, I did add it. That's why I'm finding it weird...

Offline

#4 2022-12-14 02:06:13

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 1,611
Website

Re: Difference in the transition density between Gaussian and Multiwfn.

If you indeed have correctly added pop=saveNTO, the fch file should contain NTOs instead of MOs, you can easily confirm by visualizing the orbital isosurfaces.

Offline

Board footer

Powered by FluxBB