QuantumChemistry
LocalOrbitals
creates a set of spatially localized orbitals from a set of molecular orbitals
Calling Sequence
Parameters
Description
References
Examples
LocalOrbitals(molecule, mos, options)
molecule
-
list of lists; each list has 4 elements, the string of an atom's symbol and atom's x, y, and z coordinates
mos
Matrix; canonical molecular orbitals coefficients returned as an entry of the output table with key mo_coeff from one of the electronic structure methods including 'HartreeFock', 'DensityFunctional', 'RDMFunctional', 'ActiveSpaceCI', 'ActiveSpaceSCF', 'CoupledCluster', 'FullCI', 'MP2', 'Variational2RDM', 'Parametric2RDM', and 'ContractedSchrodinger'
options
(optional) equation(s) of the form options = value where options is localization, mo_indices, or an option of the molecule such as basis or charge
The LocalOrbitals procedure generates a set of spatially localized orbitals from a set of canonical molecular orbitals (MOs).
Local orbitals offer an alternative, spatially localized rotation of the canonical MOs, which can be particularly useful for visualizing chemical bonding, improving the convergence and cost of correlated methods, and selecting active spaces in active-space methods like CI-CASSCF (ActiveSpaceCI or ActiveSpaceSCF) and V2RDM-CASSCF (Variational2RDM).
Two arguments, molecule and mos, are required where molecule defines the atom or molecule while mos is a Matrix of canonical molecular orbitals coefficients returned as the entry of the output table with the key mo_coeff from one of the electronic structure methods such as 'HartreeFock', 'DensityFunctional', 'RDMFunctional', 'ActiveSpaceCI', 'ActiveSpaceSCF', 'CoupledCluster', 'FullCI', 'MP2', 'Parametric2RDM', 'ContractedSchrodinger', or 'Variational2RDM'.
The localization string (optional) -- selects the localization technique: "Boys", "PM", and "ER". Boys Localization ("Boys") creates orbitals by keeping electron density as close together as possible in space; Pipek-Mezey Localization ("PM") generates orbitals that resemble atomic orbitals as much as possible, and Edmiston-Ruedenberg Localization ("ER") produces orbitals that minimize electron repulsion.
The mo_indices list of integers (optional) -- denotes the MOs to be localized. If not provided, the command localizes all MOs.
The command also supports options to describe the molecular calculation that produced the mos such as basis, charge, spin, symmetry, ghost, and unit. Note that it is critical to provide these keywords if you are not using their default values; for example, if mos was not produced with the default "sto-3g" basis set, you must include the basis keyword with the correct basis set.
On the Windows operating system the LocalOrbitals command requires the installation of Microsoft's Windows Subsystem for Linux (WSL). For Windows 10 (version 2004 and higher) and Windows 11 you can install the WSL by opening the Command Prompt in administrator mode and entering the command: wsl --install -d Ubuntu For additional details, please refer to: https://learn.microsoft.com/en-us/windows/wsl/install
S. F. Boys, Reviews of Modern Physics 32 (2): 296–299 (1960). "Construction of molecular orbitals to be minimally variant for changes from one molecule to another"
C. Edmiston and K. Ruedenberg, Reviews of Modern Physics. 35 (3): 457–465 (1963). "Localized Atomic and Molecular Orbitals"
J. Pipek and P. G. Mezey, The Journal of Chemical Physics. 90 (9): 4916 (1989). "A fast intrinsic localization procedure applicable for ab initio and semiempirical linear combination of atomic orbital wave functions"
After loading the QuantumChemistry package
withQuantumChemistry:
and defining the molecule,
mol ≔ MolecularGeometryhydroxyurea;
mol≔O,1.88790000,−0.00950000,−0.00270000,O,−0.48550000,−1.30440000,−0.00030000,N,0.67610000,0.69460000,0.00260000,N,−1.61350000,0.69390000,−0.00290000,C,−0.46490000,−0.07450000,0.00330000,H,0.71730000,1.68970000,0.18780000,H,−1.57330000,1.70750000,−0.00260000,H,−2.53340000,0.26600000,−0.00670000,H,2.08350000,−0.00760000,−0.95500000
let's try a DFT calculation of hydroxyurea with the "pc0" basis set
data_dft ≔ DensityFunctionalmol, basis=pc0;
Next, we use the result from the DFT calculation to compute the localized orbitals with the "ER" option
local_mo_coeff ≔ LocalOrbitalsmol, data_dftmo_coeff, basis=pc0,localization=ER;
To visualize the 20th local orbital, we replace the canonical MO coefficients in data_dft with the local MOs
data_dftmo_coeff ≔ local_mo_coeff;
And then we use the DensityPlot3D command
DensityPlot3Dmol,data_dft,basis=pc0,orbitalindex=20;
Note that the localized orbital is a p orbital on the oxygen atom connected to the central carbon atom.
See Also
Download Help Document