Dplot display of HOMO and LUMO
Clicked A Few Times
4:46:06 PM PDT - Mon, Oct 7th 2013
I have similar problems trying to visualize the MO. I'm only interested in the HOMO and LUMO using Gaussian cube file and view with ECCE. I'm not sure which orbital should I list in the "orbitals density". Furthermore, the cube file that was generated only show the charge density, not the HOMO and LUMO view. Below is my input file:

Title "Calculation"

Start Calculation

echo

charge 0

geometry autosym units angstrom
C     -2.05329     1.37998     0.0348500

C     -3.36029     0.863980     0.0358500

C     -3.59129     -0.518020     0.0168500

C     -2.52329     -1.43302     -0.00515000

C     -1.23029     -0.924020     -0.00715000

C     -1.00229     0.463980     0.0118500

C     0.456710     0.631980     0.00385000

C     0.993710     -0.668020     -0.0201500

N     -0.0282900     -1.57602     -0.0271500

C     1.27371     1.76298     0.0148500

C     2.66271     1.55798     0.00185000

C     3.20171     0.262980     -0.0231500

C     2.36871     -0.870020     -0.0341500

H     0.0807100     -2.52602     -0.0421500

H     -1.86229     2.46498     0.0498500

H     -4.21629     1.55998     0.0538500

H     -4.62829     -0.896020     0.0188500

H     -2.70029     -2.52102     -0.0211500

H     0.841710     2.77598     0.0338500

H     3.33871     2.42898     0.00985000

H     4.29671     0.129980     -0.0331500

H     2.78871     -1.88902     -0.0531500
end

ecce_print ecce.out

basis "ao basis" cartesian print
 H library "6-31G"

 C library "6-31G"

 N library "6-31G"
END

dft
 mult 1

 XC b3lyp

 mulliken
end

scf
vectors output Calculation.movecs
end

dplot
 TITLE Calculation

 vectors Calculation.movecs

  LimitXYZ

-6.0 6.0 20  

-6.0 6.0 20 

-6.0  6.0  20

 spin total

 orbitals density

 10

 30 31 32 33 34 35 36 37 38 39 

 gaussian

 output lumo.cube
end
driver
 default
end

task dft energy
task dplot

Forum Vet
5:21:06 PM PDT - Mon, Oct 7th 2013
This is what your input should look like to plot HOMO and LUMO orbitals
More details at
http://nwchemgit.github.io/index.php/Release62:DPLOT
Title "Calculation"

Start Calculation

echo

charge 0

geometry autosym units angstrom
C     -2.05329     1.37998     0.0348500
C     -3.36029     0.863980     0.0358500
C     -3.59129     -0.518020     0.0168500
C     -2.52329     -1.43302     -0.00515000
C     -1.23029     -0.924020     -0.00715000
C     -1.00229     0.463980     0.0118500
C     0.456710     0.631980     0.00385000
C     0.993710     -0.668020     -0.0201500
N     -0.0282900     -1.57602     -0.0271500
C     1.27371     1.76298     0.0148500
C     2.66271     1.55798     0.00185000
C     3.20171     0.262980     -0.0231500
C     2.36871     -0.870020     -0.0341500
H     0.0807100     -2.52602     -0.0421500
H     -1.86229     2.46498     0.0498500
H     -4.21629     1.55998     0.0538500
H     -4.62829     -0.896020     0.0188500
H     -2.70029     -2.52102     -0.0211500
H     0.841710     2.77598     0.0338500
H     3.33871     2.42898     0.00985000
H     4.29671     0.129980     -0.0331500
H     2.78871     -1.88902     -0.0531500
end

ecce_print ecce.out

basis "ao basis" cartesian print
 H library "6-31G"
 C library "6-31G"
 N library "6-31G"
END

dft
direct
 mult 1
 XC b3lyp
 vectors output Calculation.movecs
 mulliken
end

task dft

dplot
 TITLE Calculation
 vectors Calculation.movecs
  LimitXYZ
-6.0 6.0 20  
-6.0 6.0 20 
-6.0  6.0  20
 spin total
 orbitals view
 1
 45
 gaussian
 output lumo.cube
end

task dplot 

dplot
 TITLE Calculation
 vectors Calculation.movecs
  LimitXYZ
-6.0 6.0 20  
-6.0 6.0 20 
-6.0  6.0  20
 spin total
 orbitals view
 1
 44
 gaussian
 output humo.cube
end

task dplot

Clicked A Few Times
6:12:38 PM PDT - Fri, Oct 25th 2013
Thanks! I have another problem with the ecce_print. When I open the ecce output, it shows the unoptimized structure, instead of the optimized one. Do you know what happened?

Forum Vet
8:43:56 AM PDT - Mon, Oct 28th 2013
Quote:Bk.ong Oct 25th 5:12 pm
Thanks! I have another problem with the ecce_print. When I open the ecce output, it shows the unoptimized structure, instead of the optimized one. Do you know what happened?

It's hard to tell without having seen the full NWChem output

Clicked A Few Times
8:55:13 AM PST - Wed, Jan 28th 2015
I would like to make a cube file with HOMO / LUMO density that takes into consideration CCSD(T) electron correlation (not just SCF density), is there a way to do it?

I have tried to use a

vectors output filename.movecs

command inside the TCE block, but it doesn't work, is there another way?


Thanks,
Andreas

Forum Vet
12:51:23 PM PST - Wed, Jan 28th 2015
CCSD(T) is a perturbative method that generates just an energy result.

Quote:Extremis Jan 28th 7:55 am
I would like to make a cube file with HOMO / LUMO density that takes into consideration CCSD(T) electron correlation (not just SCF density), is there a way to do it?

I have tried to use a

vectors output filename.movecs

command inside the TCE block, but it doesn't work, is there another way?


Thanks,
Andreas

Just Got Here
3:03:24 PM PDT - Wed, Jun 12th 2019
Hi,

How did you determine which orbital is the HOMO and which orbital is the LUMO?

Thank you very much for your time.

Quote:Edoapra Oct 7th 5:21 pm
This is what your input should look like to plot HOMO and LUMO orbitals
More details at
http://nwchemgit.github.io/index.php/Release62:DPLOT
Title "Calculation"

Start Calculation

echo

charge 0

geometry autosym units angstrom
C     -2.05329     1.37998     0.0348500
C     -3.36029     0.863980     0.0358500
C     -3.59129     -0.518020     0.0168500
C     -2.52329     -1.43302     -0.00515000
C     -1.23029     -0.924020     -0.00715000
C     -1.00229     0.463980     0.0118500
C     0.456710     0.631980     0.00385000
C     0.993710     -0.668020     -0.0201500
N     -0.0282900     -1.57602     -0.0271500
C     1.27371     1.76298     0.0148500
C     2.66271     1.55798     0.00185000
C     3.20171     0.262980     -0.0231500
C     2.36871     -0.870020     -0.0341500
H     0.0807100     -2.52602     -0.0421500
H     -1.86229     2.46498     0.0498500
H     -4.21629     1.55998     0.0538500
H     -4.62829     -0.896020     0.0188500
H     -2.70029     -2.52102     -0.0211500
H     0.841710     2.77598     0.0338500
H     3.33871     2.42898     0.00985000
H     4.29671     0.129980     -0.0331500
H     2.78871     -1.88902     -0.0531500
end

ecce_print ecce.out

basis "ao basis" cartesian print
 H library "6-31G"
 C library "6-31G"
 N library "6-31G"
END

dft
direct
 mult 1
 XC b3lyp
 vectors output Calculation.movecs
 mulliken
end

task dft

dplot
 TITLE Calculation
 vectors Calculation.movecs
  LimitXYZ
-6.0 6.0 20  
-6.0 6.0 20 
-6.0  6.0  20
 spin total
 orbitals view
 1
 45
 gaussian
 output lumo.cube
end

task dplot 

dplot
 TITLE Calculation
 vectors Calculation.movecs
  LimitXYZ
-6.0 6.0 20  
-6.0 6.0 20 
-6.0  6.0  20
 spin total
 orbitals view
 1
 44
 gaussian
 output humo.cube
end

task dplot


Forum >> NWChem's corner >> General Topics