ChemicalBonding - Maple Help

Chemical Bonding

 Overview The nature of the chemical bond was elucidated through quantum mechanics in the 1920s and 1930s.  Here we explore the chemical bond by computing the following:     (1) the potential energy curve of a diatomic molecule (2) the polarity of the bond for a series of diatomic molecules

Chemical Bond

We load the QuantumChemistry package with Maple's with command

 > $\mathrm{with}\left(\mathrm{QuantumChemistry}\right);$
 $\left[{\mathrm{AOLabels}}{,}{\mathrm{ActiveSpaceCI}}{,}{\mathrm{ActiveSpaceSCF}}{,}{\mathrm{AtomicData}}{,}{\mathrm{BondAngles}}{,}{\mathrm{BondDistances}}{,}{\mathrm{Charges}}{,}{\mathrm{ChargesPlot}}{,}{\mathrm{CorrelationEnergy}}{,}{\mathrm{CoupledCluster}}{,}{\mathrm{DensityFunctional}}{,}{\mathrm{DensityPlot3D}}{,}{\mathrm{Dipole}}{,}{\mathrm{DipolePlot}}{,}{\mathrm{Energy}}{,}{\mathrm{ExcitationEnergies}}{,}{\mathrm{ExcitationSpectra}}{,}{\mathrm{ExcitationSpectraPlot}}{,}{\mathrm{ExcitedStateEnergies}}{,}{\mathrm{ExcitedStateSpins}}{,}{\mathrm{FullCI}}{,}{\mathrm{GeometryOptimization}}{,}{\mathrm{HartreeFock}}{,}{\mathrm{Interactive}}{,}{\mathrm{Isotopes}}{,}{\mathrm{MOCoefficients}}{,}{\mathrm{MODiagram}}{,}{\mathrm{MOEnergies}}{,}{\mathrm{MOIntegrals}}{,}{\mathrm{MOOccupations}}{,}{\mathrm{MOOccupationsPlot}}{,}{\mathrm{MOSymmetries}}{,}{\mathrm{MP2}}{,}{\mathrm{MolecularData}}{,}{\mathrm{MolecularGeometry}}{,}{\mathrm{NuclearEnergy}}{,}{\mathrm{NuclearGradient}}{,}{\mathrm{OscillatorStrengths}}{,}{\mathrm{Parametric2RDM}}{,}{\mathrm{PlotMolecule}}{,}{\mathrm{Populations}}{,}{\mathrm{RDM1}}{,}{\mathrm{RDM2}}{,}{\mathrm{RTM1}}{,}{\mathrm{ReadXYZ}}{,}{\mathrm{Restore}}{,}{\mathrm{Save}}{,}{\mathrm{SaveXYZ}}{,}{\mathrm{SearchBasisSets}}{,}{\mathrm{SearchFunctionals}}{,}{\mathrm{SkeletalStructure}}{,}{\mathrm{Thermodynamics}}{,}{\mathrm{TransitionDipolePlot}}{,}{\mathrm{TransitionDipoles}}{,}{\mathrm{TransitionOrbitalPlot}}{,}{\mathrm{TransitionOrbitals}}{,}{\mathrm{Variational2RDM}}{,}{\mathrm{VibrationalModeAnimation}}{,}{\mathrm{VibrationalModes}}{,}{\mathrm{Video}}\right]$ (2.1)

Energy of the chemical bond

Define the hydrogen fluoride (HF) molecule as a function of its internuclear distance R

 >
 ${\mathrm{HF}}{≔}\left[\left[{"H"}{,}{0}{,}{0}{,}{0}\right]{,}\left[{"F"}{,}{0}{,}{0}{,}{R}\right]\right]$ (2.1.1)

Setting R to 1.0 Angstrom (1 angstrom = 10-10 m) with the subs command, we have

 >
 ${\mathrm{HF1}}{≔}\left[\left[{"H"}{,}{0}{,}{0}{,}{0}\right]{,}\left[{"F"}{,}{0}{,}{0}{,}{1.00000000}\right]\right]$ (2.1.2)

With the Energy command we can compute the energy of HF at R = 1 angstrom

 >
 ${\mathrm{E1}}{≔}{-100.15637329}$ (2.1.3)

We can also compute the Energy as a list of R values (calculation may take a minute or so)

 >
 ${\mathrm{E_ls}}{≔}\left[\left[{0.70000000}{,}{-100.04258858}\right]{,}\left[{0.80000000}{,}{-100.12769696}\right]{,}\left[{0.90000000}{,}{-100.15593391}\right]{,}\left[{1.00000000}{,}{-100.15637329}\right]{,}\left[{1.10000000}{,}{-100.14341514}\right]{,}\left[{1.20000000}{,}{-100.12470651}\right]{,}\left[{1.60000000}{,}{-100.04862055}\right]{,}\left[{1.80000000}{,}{-100.02176025}\right]{,}\left[{2.20000000}{,}{-99.99225294}\right]{,}\left[{2.60000000}{,}{-99.98235671}\right]\right]$ (2.1.4)

We can make a plot of the energy as a function of R known as the potential energy curve

 >

(a) From the plot of the energy for HF, estimate the internuclear distance R in angstroms at which the energy reaches its minimum?

(b) Convert your result in (a) to meters.

(c)  From the plot of the energy for HF, estimate the bond dissociation energy, the energy to dissociate the HF molecule from its energy minimum?

(d) Convert the energy in (c) from hartrees (energy unit in atomic units) to kJ/mol (kilojoules/mole).

Hint: We can easily convert units in Maple.  For example, if the dissociation energy is 1.5 hartrees (or denoted as 1.5 E0), we can convert it to kJ/mol

 >
 ${{E}}_{{\mathrm{dissociation}}}{≔}{1.50000000}{}⟦{\mathrm{E0}}⟧$
 ${{N}}_{{A}}{≔}{6.02214086}{}{{10}}^{{23}}{}⟦\frac{{1}}{{\mathrm{mol}}}⟧$
 ${{E}}_{{\mathrm{dissociation}}}{≔}{3938.24945900}{}⟦{\mathrm{kJ}}⟧{}⟦\frac{{1}}{{\mathrm{mol}}}⟧$ (2.1.5)

 > 

Polarity of the chemical bond

A chemical bond is described as polar if there exists a separation of partial negative and partial positive charges.  For a diatomic molecule the polarity of the bond can be estimated in at least two ways: (1) the difference in the electronegativities of the two atoms connected by the bond and (2) the dipole moment of the molecule. We can compute and plot the dipole moment of a molecule with the Dipole and DipolePlot commands of the Quantum Chemistry Toolbox.  Note that the Toolbox uses the standard convention that the direction of the dipole moment is from the negative partial charge to the positive partial charge, which differs from some general chemistry textbooks.

The dipole moment of HF is

 >
 ${\mathrm{HF}}{≔}\left[\left[{"H"}{,}{0}{,}{0}{,}{0}\right]{,}\left[{"F"}{,}{0}{,}{0}{,}{0.91700000}\right]\right]$
 $\left[\begin{array}{cc}{"X"}& 0.000000000\\ {"Y"}& 0.000000000\\ {"Z"}& -2.02931607\end{array}\right]$
 ${\mathrm{total_dipole}}{≔}{2.02931607}$ (2.2.1)

The plot of the dipole moment of HF is (click on the plot to rotate)

 >

The dipole moment of HCl is

 >
 ${\mathrm{HCl}}{≔}\left[\left[{"H"}{,}{0}{,}{0}{,}{0}\right]{,}\left[{"Cl"}{,}{0}{,}{0}{,}{1.27400000}\right]\right]$
 $\left[\begin{array}{cc}{"X"}& 0.000000000\\ {"Y"}& 0.000000000\\ {"Z"}& -1.45323422\end{array}\right]$
 ${\mathrm{total_dipole}}{≔}{1.45323422}$ (2.2.2)

The plot of the dipole moment of HCl is (click on the plot to rotate)

 >

The dipole moment of HBr is

 >
 ${\mathrm{HBr}}{≔}\left[\left[{"H"}{,}{0}{,}{0}{,}{0}\right]{,}\left[{"Br"}{,}{0}{,}{0}{,}{1.41400000}\right]\right]$
 $\left[\begin{array}{cc}{"X"}& 0.000000000\\ {"Y"}& 0.000000000\\ {"Z"}& -1.40698952\end{array}\right]$
 ${\mathrm{total_dipole}}{≔}{1.40698952}$ (2.2.3)

The plot of the dipole moment of HBr is (click on the plot to rotate)

 >

The dipole moment of HI is

 >
 ${\mathrm{HI}}{≔}\left[\left[{"H"}{,}{0}{,}{0}{,}{0}\right]{,}\left[{"I"}{,}{0}{,}{0}{,}{1.60900000}\right]\right]$
 $\left[\begin{array}{cc}{"X"}& 0.000000000\\ {"Y"}& 0.000000000\\ {"Z"}& -0.994852279\end{array}\right]$
 ${\mathrm{total_dipole}}{≔}{0.99485228}$ (2.2.4)

The plot of the dipole moment of HI is (click on the plot to rotate)

 >

(e) Use the total (or net) dipole moments to arrange HF, HCl, HBr, and HI from lowest to highest bond polarity.

(f) Use the differences in electronegativity of A and B in AB to arrange HF, HCl, HBr, and HI from lowest to highest bond polarity.

(g) Are the results in (e) consistent with the results in (f)?

Hint: For part (f) the electronegtaivities of atoms can be directly obtained in the Quantum Chemistry Toolbox with the AtomicData command

 > $\mathrm{AtomicData}\left(F\right)\left[\mathrm{electronegativity}\right];$
 ${3.98000000}$ (2.2.5)

 >