The DensityPlot3D procedure generates a 3D plot of a molecule's electron density.

The procedure can generate the electron density of the entire molecule or a designated molecular orbital.

Two arguments, molecule and data, are required where molecule defines the atom or molecule while data is the Maple table returned by one of the electronic structure methods including 'HartreeFock', 'DensityFunctional', 'ActiveSpaceCI', 'ActiveSpaceSCF', 'CoupledCluster', 'FullCI', 'MP2', 'Parametric2RDM', 'ContractedSchrodinger', and 'Variational2RDM'.

The orbitalindex integer (optional) -- denotes the molecular orbital to plot; default is 0 which plots the molecule's density.

The densitycutoff float (optional)  -- sets the cutoff for the electron density contour.  Typically, less than 0.01; default is 0.005.

The phase boolean (optional)  -- displays the phase of the molecular orbital when set to true; default is true.

The vanderwaalsradii float (optional)  -- is the weight (>1) for Van der Waals radii for defining the density grid; default is 1.8.

The gridspacing float (optional)  -- is the spacing between grid points; default is 0.1.

The maximumpoints integer (optional) -- is the maximum total number of grid points in the density plot; default is 100000.

The bonds set of sets (optional) is a set of sets where each set contains two integers denoting the atoms to be connected by a chemical bond.  The integers correspond to the positions of the atoms in molecule.  If the optional keyword bonds is not set, the bonds to be shown are determined automatically.

The colors list of strings (optional) is a list of two strings containing the colors of the positive and negative lobes of the molecular orbitals. Any color strings from Maple's ColorTools package is supported.  If the optional keyword colors is not set, the plot uses the deault colors, "green" and "purple".

$\mathrm{with}\left(\mathrm{QuantumChemistry}\right)\:$

$\mathrm{molecule} ≔ \left[\left[''H''\,0\,0\,0\right]\,\left[''F''\,0\,0\,0.95\right]\right]\;$

$\mathrm{molecule}≔\left[\left[''H''\,0\,0\,0\right]\,\left[''F''\,0\,0\,0.95000000\right]\right]$

$\mathrm{output\_hf} ≔ \mathrm{HartreeFock}\left(\mathrm{molecule}\, \mathrm{basis}\=''dz''\right)\;$

$\mathrm{table}\left(\mathrm{\%id}\=36893488450759149660\right)$

$\mathrm{DensityPlot3D}\left(\mathrm{molecule}\,\mathrm{output\_hf}\, \mathrm{basis}\=''dz''\, \mathrm{orbitalindex}\=3\, \mathrm{densitycutoff}\=0.0005\right)\;$

$\mathrm{DensityPlot3D}\left(\mathrm{molecule}\,\mathrm{output\_hf}\, \mathrm{basis}\=''dz''\, \mathrm{orbitalindex}\=3\, \mathrm{densitycutoff}\=0.0005\,\mathrm{colors}\=\left[''Nautical\; Red''\,''Nautical\; Light\; Blue''\right]\right)\;$