See Visualization page for information on ways to plot the different families of orbitals and electron/spin densities discussed here. Mulliken/Löwdin/Hirshfeld charge and spin population analysis Mulliken, Löwdin and Mayer population analysis is automatically printed in a single-point energy calculation and the beginning and end of a geometry optimization. ORCA reduces printing in geometry optimization steps by default (prints full output in the first and last steps only). If more information is needed in each optimization step (for example to monitor the electronic structure through the population analysis during the optimization) use this: Density-based Hirshfeld population analysis is requested by the following output block option: ESP charges Coming... - CHELPG Other population analysis via other programs The Multiwfn program, a wavefunction analysis program, can be used to perform all sorts of wavefunction/density analysis. Program can do topological analysis of the electron density, AIM, ELF, Laplacian and population analysis like: Hirshfeld, VDD, Mulliken, Löwdin, Modified MPA (including three methods: SCPA, Stout & Politzer, Bickelhaupt), Becke, ADCH (Atomic dipole moment corrected Hirshfeld), CHELPG, Merz-Kollmann and AIM. Requires creation of a Molden file: orca_2mkl file -molden The Horton program can perform all sorts of real-space atoms-in-molecule analysis, ESP fitting and more. Requires creation of MKL file: orca_2mkl file Calculating localized orbitals The canonical orbitals from an SCF calculation are often difficult to interpret as they tend to be rather delocalized. Unitary transformations among the occupied orbitals can easily be performed without changing the total wavefunction, thus resulting in a different sets of orbitals that may be easier to interpret and are just as valid as the canonical ones. One can e.g. easily request a localized orbital analysis to be performed automatically after convergence of an SCF calculation. The disadvantage of localized orbitals is that they do not have well-defined orbital energies. Standard options: Pipek-Mezey, Foster-Boys Also: Intrinsic atomic orbitals and intrinsic bond orbitals (IAOIBO) by Gerald Knizia. See much more information about orca_loc in Chapter 9.32.5 in the ORCA 4.0 Manual. One can also perform the localization step afterwards using the orca_loc program that further allows one to only localize a certain range of orbitals. where loc-input.inp contains: To get a Löwdin orbital analysis of the localized orbitals you can read them in without iterations (Noiter) using a separate inputfile and print using Normalprint. Remember to specify the same basis set as in the previous calculation to prevent basis set projection: Analyzing corresponding orbitals Corresponding orbitals are automatically printed after a Broken-Symmetry-Flipspin calculation but can also be requested by: To get a Löwdin orbital analysis of the corresponding orbitals (stored in jobname.uco) you can read them in without iterations (Noiter) using a separate inputfile and print using Normalprint. Remember to specify the same basis set as in the previous calculation to prevent basis set projection: Using UNOs and QROs Both unrestricted natural orbitals (UNOs) and quasi-restricted orbitals (QROs) are calculated after SCF convergence if the UNO keyword is chosen: To get a Löwdin orbital analysis of the UNOs (stored in jobname.uno) or QROs (stored in jobname.qro) you can read them in without iterations (Noiter) in using a separate inputfile and print using Normalprint. Remember to specify the same basis set as in the previous calculation to prevent basis set projection: Natural Transition Orbitals (NTOs) Natural transition orbitals can be useful to understand the change in electronic structure between excited states. See the TDDFT page. Natural Bond Orbitals (NBO) and Natural Population analysis (NPA) NBO/NPA analysis in ORCA requires the NBO6 code that must be purchased from the NBO developers. For NBO analysis to work in ORCA, at runtime the NBOEXE variable must be set to the full path of the NBO executable. NBOEXE=/path/to/nbo.exe Note: The JANPA project allows the use of NBOs and NPA via a freeware program package. See also Use of JANPA with ORCA. If commercial NBO has been set up correctly, one can request NBO analysis like this: To get a Löwdin orbital analysis of the natural bond orbitals (NBOs) you can read them in (stored in jobname.nbo) without iterations (Noiter) using a separate inputfile and print using Normalprint. Remember to specify the same basis set as in the previous calculation to prevent basis set projection: Fractional occupation density (FOD) analysis For quantifying and visualizing static electron correlation or multireference effects. See FOD page. Atoms-In-Molecules (AIM) analysis It is possible to do AIM analysis on ORCA-calculated electron densities. Currently, only for closed-shell wavefunctions. The orca_2aim program reads an ORCA gbw file and converts it into a WFN file that can then be read by other programs for a topological analysis of the electron density (such as AIMAll). It is also possible to have ORCA create the WFN file automatically in an ORCA inputfile. Simply add the AIM keyword to the Simpleinput line in this case. If only AIM (Bader) charges are wanted, one could first create a .cube file of the electron density using orca_plot, See Visualization and printing, and then use the Bader analysis program by Henkelman and Jónsson. Note that a fine grid-density may be required for accurate charges. Electrostatic potential The orca_vpot program calculates the electrostatic potential for a set of points provided by the user. Example: ORCA filename.inp orca_vpot usage: Here filename.gbw is the GBW file from the previous energy calculation (that has to include the Keepdens keyword), filename.scfp contains the electron density matrix (exists because of Keepdens keyword; alternatively orca_plot can recalculate the scfp file), filename.vpot.xyz contains the set of user-defined points where the potential should be calculated and filename.vpot.out contains the electrostatic potential (4th column contains the value) filename.vpot.xyz should look like below where the first line contains the number of points, followed by XYZ coordinates (in Å) for the desired points. The outputfile, here filename.vpot.out, will contain the electrostatic potential as the fourth column in atomic units. To visualize the electostatic potential a Cube file can be created for visualization in various programs. There is a convenient Python program written by Marius Retegan to help with this. See here |
|
|