LocalOrbitals - Maple Help
For the best experience, we recommend viewing online help using Google Chrome or Microsoft Edge.

Online Help

All Products    Maple    MapleSim


QuantumChemistry

  

LocalOrbitals

  

creates a set of spatially localized orbitals from a set of molecular orbitals

  

 

Calling Sequence

Parameters

Description

References

Examples

Calling Sequence

LocalOrbitals(molecule, mos, options)

Parameters

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

Description

• 

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

References

1. 

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"

2. 

C. Edmiston and K. Ruedenberg, Reviews of Modern Physics. 35 (3): 457–465 (1963). "Localized Atomic and Molecular Orbitals"

3. 

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"

Examples

After loading the QuantumChemistry package

withQuantumChemistry:

and defining the molecule,

 mol  MolecularGeometryhydroxyurea;

molO,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

(1)

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

QuantumChemistry

 

 

 


Download Help Document