Introduction To Quantum Chemistry Simulations With ORCA: HPRC Short Course Series
Introduction To Quantum Chemistry Simulations With ORCA: HPRC Short Course Series
%pal
nprocs 2
end
*gzmt 0 1
H
H 1 0.8
*
Details about the input file lines:
(ORCA Input files are pretty much free-format. Blank lines are allowed, input is
usually not case-sensitive. One can create comment lines by adding the "#" symbol
to a line.)
! start a keyword line. In the Simple input syntax, keywords are added in any order to
the line beginning with "!". Multiple "!" lines are also allowed.
B3LYP– use the B3LYP density functional for the calculation (method)
OPT – perform a geometry optimization
FREQ – perform an analytical frequency calculation (if opt and freq are both present,
an optimization is performed 1st and then a frequency calculation on the
optimized structure)
def2-SVP – use the def2-SVP basis set
# Comment lines can be added if the line begins with a #.
% Start a block input. Advanced settings are often specified using the block input for
different modules.
8. On your desktop, double click to open the orca/water directory. Right click to open the
terminal and load the module (moduleVloadVsc/orca¿). Open an Avogadro window
(avogadro&¿) and draw a water molecule. Make sure you only have one water
molecule in the view. Click Extensions ->Orca ->Generate Orca input. In the preview
window, make changes manually so that your input looks similar to this. (Your
coordinates may be different from mine.)
!BP86 def2-SVP opt freq
!largeprint
%pal
© 2020 Xin Yang LMS & HPRC 2
Introduction to Quantum Chemistry Simulations with ORCA
HPRC Short Course Series
nprocs 2
end
* xyz 0 1
O -3.56626 1.77639 0.00000
H -2.59626 1.77639 0.00000
H -3.88959 1.36040 -0.81444
*
9. Click Generate and find your home/username/Desktop/orca/water/ directory. Save the
file as water.inp. Close Avogadro and discard changes.
10. Submit the calculation using terminal
($ORCA_HOME/orcaVwater.inp>water.out¿).
11. Wait for the job to finish (<1 minutes) and open water.out in file browser. You may
need to reload to see the output. Check if your job terminate normally. Close the
output before move to the terminal.
12. To visualize results with Molden, type in prompt: moldenVwater.out¿
13. Change the draw mode by clicking Solid->Stick under
Draw Mode. To change a perspective, left-click, hold, and
drag to rotate the molecule.
14. Click on the Distance button and then select the
hydrogen atom and then the oxygen atom (You need to
click at the very center of the atom or Molden does not
respond). In the pop-up Calculate window, click the
monitor button. You should see the bond distance
appear in the window on the molecule.
15. To visualize the change in the geometry that occurred
during the geometry optimization, click on the Geom. conv.
button. The Geom. Convergence window has a chart of
energy vs. point, which shows how the energy of the
molecule change as it was optimized and how many
geometry optimization steps it took. Move the new window
over such that you can see it and the molecule at the same
time and then click on the Next button. You should see the
geometry change. NOTE: Molden always starts by
displaying the input structure (starting structure) and NOT the optimized structure.
After you click the Movie button, molden will display the last geometry in the out file,
which will be the optimized structure if the optimization finished normally. Close the
Geom. Convergence window.
16. Exit molden by clicking on the button with the skull and cross-bones.
17. (Only for demonstration/Optional) Convert result into Molden file to visualize
molecular orbitals. Type in prompt: $ORCA_HOME/orca_2mklVwaterV-molden¿ ;
then moldenVwater.molden.input¿; In molden control panel, click Dens. Mode,
Orbital, Space, change contour value to 0.2 and apply. You should be able to see plots
of MOs. Switch back to Mol. Mode and exit molden.
21. Extract the geometry optimization convergence criteria at each geometry step and after
the frequency calculation by typing at the prompt: grepV-A17V'Geometry convergence'
water.out¿
The last set of lines is the convergence for the last step of geometry optimization.
Did the geometry optimization reach convergence? You should be able to find
this pattern in the output.
***********************HURRAY********************
*** THE OPTIMIZATION HAS CONVERGED ***
*************************************************
22. Extract the vibrational frequencies by typing at the prompt: grepV-
A15VFREQUENCIESV water.out¿
The first six frequencies are 0 (ORCA counts from 0). You should have all real
frequencies for a minimum (or 1 imaginary frequency labeled with a minus sign
for a 1st order saddle point, i.e. transition state). If you have more than 1
imaginary mode then you have optimized to a higher order saddle point on the
potential energy surface or the geometry optimization was not converged.
23. Extract pure electronic energy E, Ho and Go by typing at the prompt (In ORCA, G
=H – T*S is referred to as Gibbs free enthalpy, not energy):
grepV'ElectronicVenergy'Vwater.out >> water.txt¿
grepV'TotalVEnthalpy' water.out >> water.txt¿
© 2020 Xin Yang LMS & HPRC 4
Introduction to Quantum Chemistry Simulations with ORCA
HPRC Short Course Series
grepV'FinalVGibbsVfreeVenthalpy' water.out >> water.txt¿
24. Reload the files , view the water.txt file and exit the terminal (type in prompt:
¿
exit ). (Please take the poll in zoom and take a short break:)
%pal
nprocs 2
end
%TDDFT
NROOTS 30
END
*xyzfile 0 1 propenal.xyz
27. Right click to open a terminal under the propenal/ directory and load the orca
module (moduleVloadVsc/orca¿). Submit the td.inp to ORCA.
$ORCA_HOME/orcaVtd.inp > td.out¿ )
28. When the calculation finshes, open
propenal_td.out with Avogadro (type in
prompt: avogadroVtd.out&¿)and check
the spectra from memu Extensions >
Spectra. Change to Absorption
spectrum and adjust the range and peak
width as you need. (These resuls
usually does not match the experiment
very well. It is much more prevalent that
the predicted spectra is "shifted" by ±0.1
to ±0.5 eV and it is costumary in the
literature to fix that, if there is any
reference, to correct systematic errors. )
29. To find the excitation information in the output file of the TD-DFT calculations, go
to the file browser and double click to open td.out, use search tool to find EXCITED.
30. Below the energies (given in au, eV, and cm-1), the contributions of each single
excitation are printed, first with its relative contribution, followed by the eigenvector
value.
STATE 1: E= 0.131767 au 3.586 eV 28919.5 cm**-1
14a -> 15a : 0.979239 (c= -0.98956485)
14a -> 16a : 0.018814 (c= -0.13716449)
© 2020 Xin Yang LMS & HPRC 5
Introduction to Quantum Chemistry Simulations with ORCA
HPRC Short Course Series
Excited State 1 composed of about 98% a HOMO to LUMO transition, from orbital 14
to orbital 15 (the "a" there means alpha orbital). Therefore, you could visualize MO 15
and 16 in Avogadro to determine their character and symmetry. (Note: ORCA counts
orbital from 0 while Avogadro counts from 1)
31. Scroll down a few pages, you will be able to find the states and their respective
vertical transition energies under the “ABSORPTION SPECTRUM VIA TRANSITION
ELECTRIC DIPOLE MOMENTS” section, together with the oscillator strength (fosc,
proportional to the extinction coefficient) and transition dipole moments. From the
calculation, we only get transition lines, not bands, as spectrum measured from
experiment. (Please take the poll in zoom:)
32. (Optional) Find excited states with the relatively large oscillator strength, identify
excitation(s) that has major contribution to the excited state, and visualize
corresponding orbitals with Avogadro.
33. When finished, close Avogadro and exit the terminal. Click System from the
menu listed on top of the screen to log out your account.