compute the spin-free 2-electron reduced density matrix (2-RDM)
RDM2(molecule, method, options)
list of lists; each list has 4 elements, the string of an atom's symbol and atom's x, y, and z coordinates
(optional) method = name/procedure where name is one of 'MP2' default, 'ActiveSpaceCI', 'ActiveSpaceSCF', 'CoupledCluster', 'FullCI','MP2', 'Variational2RDM', 'Parametric2RDM', 'ContractedSchrodinger'
(optional) equation(s) of the form option = value where option is any valid option of the chosen method
The RDM2(molecule) command computes the spin-free 2-electron reduced density matrix (2-RDM).
The output is a 4-index Array with dimensions r x r x r x r, where r is the number of the orbitals.
The element A[i,k,j,l] of the output Array correponds to the element of the 2-RDM D2[i,j,k,l].
The default method is the 'MP2' method.
The result depends upon the chosen molecule, method, and basis set among other options such as charge, spin, and symmetry.
Because the methods employ Maple cache tables, the procedure only computes the 2-RDM if it has not been previously computed by calling the method directly or indirectly through another property.
Consider the molecule hydrogen fluoride
molecule ≔ H, 0, 0, 0, F, 0, 0, 0.95 ;
Compute the 2-RDM from the parametric 2-RDM method in the double-zeta ("dz") basis set
d2 ≔ RDM2molecule, method=Parametric2RDM,basis=dz;
Download Help Document