Difference between revisions of "Mod:Hunt Research Group:multiwfn mp2 density"

From wiki
Jump to navigation Jump to search
Line 1: Line 1:
 
==Using MP2 Densities in Surface Analysis with Multiwfn==
 
==Using MP2 Densities in Surface Analysis with Multiwfn==
  
* Suppose that a given molecular system has been optimized at the MP2/6-311G(d,p) level of theory and you want to calculate the value of (V_s,max) on a specific atom using MP2 density. For this either you need to generate a .wfn/.wfx file with density=current or density=mp2 keyword in the input file route line or you can do it from the .fchk file as shown below:
+
*for a system at the MP2 (or other non-HF, non-DFT) level of theory  
* The route card must have density=current or density=mp2 keyword, then create the .fchk file from chk file.
+
*either generate the wfn/wfx file as part of the G16 job or using the *.fchk file
* Load the fchk file into Multiwfn then follow these steps for Multiwfn to use the mp2 density instead of the SCF density.
+
*BOTH cases you must use density=current or density=mp2 keyword  
 +
*you can see what density is available by running
 +
<pre>
 +
grep -i '*.fck'
 +
<*pre>
 +
*for example
 +
<pre>
 +
grep -i 'dens' HNO3_Dimer_c6_MP2_augccpVT+dZ_SP_wfx.fchk
 +
# output=wfx density=mp2 mp2=full/gen int=ultrafine scf=conv
 +
RMS Density                                R      6.409540209043613E-10
 +
Total SCF Density                          R  N=      85905
 +
Total MP2 Density                          R  N=      85905
 +
</pre>
 +
*from the *.fchk, load the *.fchk file into Multiwfn  
 +
*then follow these steps
  
<pre>
+
::select <pre>200 Other functions (Part 2)</pre>
200 Other functions (Part 2)
+
::16 <pre>Generate natural orbitals based on the density matrix in .fch/.fchk file</pre>
16 Generate natural orbitals based on the density matrix in .fch/.fchk file
+
::Input the type of density matrix, e.g. SCF, CI, CC, MP2, MP3, MP4...
Input the type of density matrix, e.g. SCF, CI, CC, MP2, MP3, MP4...
 
 
  e.g. If the .fch was produced by Gaussian at MP2 level and "density" keyword has been used, you may input "SCF" or "MP2"
 
  e.g. If the .fch was produced by Gaussian at MP2 level and "density" keyword has been used, you may input "SCF" or "MP2"
 
  Hint: You can manually open .fch file by text editor and search "density" to see which density matrices are available in this file
 
  Hint: You can manually open .fch file by text editor and search "density" to see which density matrices are available in this file

Revision as of 23:31, 13 July 2025

Using MP2 Densities in Surface Analysis with Multiwfn

  • for a system at the MP2 (or other non-HF, non-DFT) level of theory
  • either generate the wfn/wfx file as part of the G16 job or using the *.fchk file
  • BOTH cases you must use density=current or density=mp2 keyword
  • you can see what density is available by running
grep -i '*.fck'
<*pre>
*for example
<pre>
grep -i 'dens' HNO3_Dimer_c6_MP2_augccpVT+dZ_SP_wfx.fchk
# output=wfx density=mp2 mp2=full/gen int=ultrafine scf=conv
RMS Density                                R      6.409540209043613E-10
Total SCF Density                          R   N=       85905
Total MP2 Density                          R   N=       85905
  • from the *.fchk, load the *.fchk file into Multiwfn
  • then follow these steps
select
200 Other functions (Part 2)
16
Generate natural orbitals based on the density matrix in .fch/.fchk file
Input the type of density matrix, e.g. SCF, CI, CC, MP2, MP3, MP4...
e.g. If the .fch was produced by Gaussian at MP2 level and "density" keyword has been used, you may input "SCF" or "MP2"
Hint: You can manually open .fch file by text editor and search "density" to see which density matrices are available in this file

MP2

Done! Basis function information now correspond to natural orbitals

If next you intend to analyze real space functions based on the NOs, you should export new.mwfn in current folder and then reload it, so that GTF information will also correspond to NOs
Would you like to do this immediately? (y/n)

n

0 Return

  • Then you can carry out the analysis you want. They density used now is the MP2 density