<![CDATA[Multiwfn forum / Multiwfn output differs from ORCA output]]> //www.umsyar.com/wfnbbs/viewtopic.php?id=795 The most recent posts in Multiwfn output differs from ORCA output. Tue, 21 Mar 2023 13:05:47 +0000 FluxBB <![CDATA[Re: Multiwfn output differs from ORCA output]]> //www.umsyar.com/wfnbbs/viewtopic.php?pid=3058#p3058

Thank you.
Now everything is in the correct magnitude.

dummy@example.com (tdog123) Tue, 21 Mar 2023 13:05:47 +0000 //www.umsyar.com/wfnbbs/viewtopic.php?pid=3058#p3058
<![CDATA[Re: Multiwfn output differs from ORCA output]]> //www.umsyar.com/wfnbbs/viewtopic.php?pid=3057#p3057

The unit is Bohr

dummy@example.com (sobereva) Tue, 21 Mar 2023 12:13:53 +0000 //www.umsyar.com/wfnbbs/viewtopic.php?pid=3057#p3057
<![CDATA[Re: Multiwfn output differs from ORCA output]]> //www.umsyar.com/wfnbbs/viewtopic.php?pid=3056#p3056

Thank you for that.

I took the coordinates of the atoms of the molden.input file but the value for the core dipole moment is way to large.
Do you know what the unit of the coordinates in the molden.input file is?

dummy@example.com (tdog123) Tue, 21 Mar 2023 10:53:00 +0000 //www.umsyar.com/wfnbbs/viewtopic.php?pid=3056#p3056
<![CDATA[Re: Multiwfn output differs from ORCA output]]> //www.umsyar.com/wfnbbs/viewtopic.php?pid=3041#p3041

Result:

 Note: In below output the case of i=j corresponds to contribution of electron t o dipole moment of excited state i Transition electric dipole moment between excited states (a.u.): i j X Y Z Diff.(eV) Oscil.str 1 1 8.1128566 27.5325457 -11.1031678 0.00000 0.00000 1 2 -0.3564785 -0.9319312 -0.2005617 0.49400 0.01254 1 3 0.1879258 0.6038274 0.1079457 1.57900 0.01592 2 2 9.0731804 27.3089606 -11.2001798 0.00000 0.00000 2 3 0.0290171 -0.4554603 -0.0817156 1.08500 0.00571 3 3 6.8035680 24.7108321 -11.8004531 0.00000 0.00000

As mentioned in the prompt, the data like 1->1 only corresponds to electron contribution to the transition dipole moment. Therefore, if the center of nuclear charges is not at origin, you need to manually add contribution of nuclear charges to it to obtain electric dipole moment of excited state.

dummy@example.com (sobereva) Tue, 14 Mar 2023 12:47:50 +0000 //www.umsyar.com/wfnbbs/viewtopic.php?pid=3041#p3041
<![CDATA[Re: Multiwfn output differs from ORCA output]]> //www.umsyar.com/wfnbbs/viewtopic.php?pid=3040#p3040

I got the same result.
Thank you.

And what do you get for the permanent dipole moments (0 -> 0, 1 -> 1, 2 -> 2) and for the transition dipole moment from S1 to S2?
The permanent dipole moments are rather big.

dummy@example.com (tdog123) Tue, 14 Mar 2023 11:46:31 +0000 //www.umsyar.com/wfnbbs/viewtopic.php?pid=3040#p3040
<![CDATA[Re: Multiwfn output differs from ORCA output]]> //www.umsyar.com/wfnbbs/viewtopic.php?pid=3036#p3036

I have checked your file, the result printed by subfunction 5 of main function 18 Multiwfn is in good agreement with ORCA

Those printed by Multiwfn:

 Transition electric dipole moment between ground state (0) and excited states ( a.u.) i j X Y Z Diff.(eV) Oscil.str 0 1 1.2200139 1.7805647 0.3778114 2.26800 0.26680 0 2 -0.6526576 -0.1042663 -0.0162139 2.76200 0.02958 0 3 -0.9657421 0.9433320 0.3632050 3.84700 0.18421

Those outputted by ORCA

----------------------------------------------------------------------------- ABSORPTION SPECTRUM VIA TRANSITION ELECTRIC DIPOLE MOMENTS ----------------------------------------------------------------------------- State Energy Wavelength fosc T2 TX TY TZ (cm-1) (nm) (au**2) (au) (au) (au) ----------------------------------------------------------------------------- 1 18292.8 546.7 0.265860026 4.78462 -1.21955 -1.77634 -0.37675 2 22279.3 448.8 0.029599100 0.43737 0.65295 0.10384 0.01572 3 31027.8 322.3 0.183715970 1.94927 0.96454 -0.94235 -0.36181

For example, the Y component of transition of S0-S1 given by Multiwfn is 1.7805647, the counterpart given by ORCA is -1.77634, the difference is quite small. The sign is different, because Multiwfn and ORCA employ different convention of electric dipole operator, this doesn't matter. If you hope to make the agreement even better, you can use TPRINT 1E-10, but the output file of ORCA will be larger.

dummy@example.com (sobereva) Fri, 10 Mar 2023 10:32:51 +0000 //www.umsyar.com/wfnbbs/viewtopic.php?pid=3036#p3036
<![CDATA[Re: Multiwfn output differs from ORCA output]]> //www.umsyar.com/wfnbbs/viewtopic.php?pid=3035#p3035

Ok.

I send you a dropbox link. Is this ok?

https://www.dropbox.com/s/6gxb1ad56w85q … input?dl=0

https://www.dropbox.com/s/vs9pqxmem11wl … 2.out?dl=0

dummy@example.com (tdog123) Fri, 10 Mar 2023 09:56:10 +0000 //www.umsyar.com/wfnbbs/viewtopic.php?pid=3035#p3035
<![CDATA[Re: Multiwfn output differs from ORCA output]]> //www.umsyar.com/wfnbbs/viewtopic.php?pid=3034#p3034

There is not evident problem from your input file.

To reproduce your observation, I need your .molden.input file and ORCA output file, you may compress and share them via file hosting service and paste downloading link here, I will check. Alternatively, send them to my E-mail.

dummy@example.com (sobereva) Fri, 10 Mar 2023 08:52:00 +0000 //www.umsyar.com/wfnbbs/viewtopic.php?pid=3034#p3034
<![CDATA[Re: Multiwfn output differs from ORCA output]]> //www.umsyar.com/wfnbbs/viewtopic.php?pid=3033#p3033

Hi Tian Lu,

thank you for your quick response.
I tried your advices but the results are still not as aspected.

This is the input-file I used:

! wB97X-D3BJ def2-SVP def2/J RIJCOSX TightSCF


%pal nprocs 32
end

%TDDFT 
    NROOTS    3
    TDA    true
TPRINT 1E-8
END

%output
PrintLevel Normal
#Print[ P_MOs ] 1
#Print[ P_Overlap ] 1
end

#%cpcm
#epsilon 4.0
#end


* xyzfile 0 1 cla2_opt.xyz

dummy@example.com (tdog123) Fri, 10 Mar 2023 08:42:39 +0000 //www.umsyar.com/wfnbbs/viewtopic.php?pid=3033#p3033
<![CDATA[Re: Multiwfn output differs from ORCA output]]> //www.umsyar.com/wfnbbs/viewtopic.php?pid=3022#p3022

Two possible reasons:
(1) You didn't specify "tprint 1E-8" in %tddft
(2) You performed TDDFT rather than TDA-DFT, in the former case the result should be inaccurate if de-excitation configurations contribute notably

For details, please check "ORCA users" in Part 3 of Section "3.21.A Basic information about electron excitation analysis module" in Multiwfn manual.

dummy@example.com (sobereva) Wed, 08 Mar 2023 12:40:03 +0000 //www.umsyar.com/wfnbbs/viewtopic.php?pid=3022#p3022
<![CDATA[Re: Multiwfn output differs from ORCA output]]> //www.umsyar.com/wfnbbs/viewtopic.php?pid=3020#p3020

I meant that I calculated with ORCA the electric transition dipole moments and afterwards all dipole moments with Multiwfn,

dummy@example.com (tdog123) Wed, 08 Mar 2023 11:55:48 +0000 //www.umsyar.com/wfnbbs/viewtopic.php?pid=3020#p3020
<![CDATA[Multiwfn output differs from ORCA output]]> //www.umsyar.com/wfnbbs/viewtopic.php?pid=3019#p3019

Hello everybody,

I calculated transition and permanent, magnetic and electric dipole moments with ORCA and afterwards with Multiwfn.
My problem is that the results differ from the ORCA output file, i.e. the electric transition dipole moments (0 -> 1, 0 -> 2, ...).
Can someone tell me what Multiwfn is exactly doing and why those results are different?
The other results are not in the expected range (comparison with the experiment) as well.

Thanks in advance.

best regards,
Thomas

dummy@example.com (tdog123) Wed, 08 Mar 2023 11:54:10 +0000 //www.umsyar.com/wfnbbs/viewtopic.php?pid=3019#p3019
Baidu
map