Skip to main content
Chemistry LibreTexts

#7B Homework Solution: ab initio Calculations

  • Page ID
    120327
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \) \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)\(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\) \(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\)\(\newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    Name: ______________________________

    Section: _____________________________

    Student ID#:__________________________

    This homework is a little different. We will use the cool web base ab initio site by Perri at Sonoma State U. Download the following paper and review the concept (http://pubs.acs.org/doi/pdf/10.1021/ed5004228). Follow the directions on this tutorial (you will not need to install Avogadro for this HW). We will use this package multiple times in this class so make sure you can complete the online tutorial if you are confused.

    To start the calculations

    • click here: https://chemcompute.sonoma.edu/index.html
    • sign in: Username: "Chem110A" and password: "QuantumRocks1!"
    • click "student mode"
    • click on red/blue molecular orbital (under GAMESS)
    • click on "submit" on top right

    We are going to calculate the binding energy of the first two row of atoms with the ChemCompute software. Before beginning, we need to discuss an advanced aspect of electronic wavefuctions: Spin. You know from your general chemistry class that electrons have a "spin" that is related to the \(m_s\) quantum number that can be either "up" with \(m_s = +1/2\) or "down" with \(m_s = -1/2\). Review your electron configuration rules (especially Hund's rules) if you are not familiar with this.

    The multiplicity of an energy state is defined as \(2S+1\), where \(S\) is the total spin angular momentum and is determined by adding the \(m_s\) quantum numbers of all electrons in the system (atom or molecule):

    \[S =\sum m_s\]

    The spin multiplicity of the system will be discussed in detail later on in course. States with multiplicity 1, 2, 3, 4, 5 are respectively called singlets, doublets, triplets, quartets and quintets. The triplet state indicates that the multiplicity \(2S+1 = 3\), so that the total spin \(S = 1\). This spin is due to two unpaired electrons, as a result of Hund's rule which favors the single filling of degenerate orbitals.

    For example, the electronic configuration of the ground state of the nitrogen atom is \(1s^22S^22P^3\) and has one electron in each of the 2p orbitals (each with the same spin so same \(m_s\) value, e.g., +1/2. So for this system

    \[S = + 1/2 -1/2 +1/2 - 1/2 +1/2 + 1/2 +1/2 = 3/2\]

    and the multiplicity is \(2 \times 3/2 +1 = 4\); the ground state nitrogen atom is a quartet state.

    You will want to prepare the table below with the appropriate details before you do the quantum calculation.

    Building the first Two rows of the Periodic Table

    Run the single-point calculations on the first ten atoms of the periodic table (Hydrogen to Neon) and fill in the below table. You only need to extract the total energy (binding energy) of each system for the table, but you need to address the spin of the system explicitly (those are the rules of quantum calculations).

    • If the spin multiplicity of the atom is not 1, then use the UHR (Unrestricted Hartree-Fock) method,
    • If the spin multiplicity of the atom is 1, you can use either the RHF (Restricted Hartree-Fock) or UHR method (should give the same results).

    Use the 6-311G** basis set for both. If you select the RHF and have a non-unity spin multiplicity, the program will give you an error.

    Pull in the 03166 input file (job# 3166) to get the helium atom. The hydrogen atom is already done for you

    =3162.png

    When you recalculate a new atom you will have to change the multiplicity for the molecule (from the table you are making below) and then you "create" the input file and swap to the next atom in the input (see arrow above) by changing these two things

    1. the name of the atom to the standard chemical abbreviation and
    2. the number after that is the atomic number \(Z\).
    Atom Nuclear Charge (\(Z\)) Number of Electrons \(n\) Electronic Configuration Total electron spin \(S = \sum m_s\) Spin Multiplicity \(2S+1\) Binding Energy \(E_Z\) \(n Z^2 E_{1} \)
    Hydrogen 1 1 \(1s^1\) 1/2 2 -313
    Helium 2 2 \(1s^2\) 0 1 -1794
    Lithium 3 3 \(1s^22s^1\) 1/2 2 -4663
    Beryllium 4 4 \(1s^22s^2\) 0 1 -9143
    Boron 5 5 \(1s^22s^22p^1\) 1/2 2 -15392
    Carbon 6 6 \(1s^22s^22p^2\) 1 3 -23650
    Nitrogen 7 7 \(1s^22s^22p^3\) 3/2 4 -34135
    Oxygen 8 8 \(1s^22s^22p^4\) 1 3 -46940
    Fluorine 9 9 \(1s^22s^22p^5\) 1/2 2 -62371
    Neon 10 10 \(1s^22s^22p^6\) 0 1 -80648

    To be Submitted for Credit

    It is VERY useful to review these two pages to interpret your results for credit.

    1. http://chem.libretexts.org/Under_Con...rons_and_Light
    2. http://chem.libretexts.org/LibreText...Slater's_Rules
    • Submit the above table filled out for the total binding electronic energy of each atom (do not copy from others as if we will compare tables)
    • Plot up this binding energy \(E_Z\) vs. \(n\) on your favorite software (e.g., excel or matlab)
    • Plot up the \(nE_1\) vs. \(n\) on same plot
    • From these data, answer these questions:
      • What is the definition of binding energy (i.e., total electronic energy) of these atoms?
      • Why are your binding energies negative?
      • What is the meaning of \(nZ^2E_{1}\) expression?
        • What is the origin of the \(n\) in this expression?
        • What is the origin of the \(Z^2\) in this expression?
        • What is the origin of the \(E_{1}\) in this expression?
      • What is the origin of difference between \(nZ^2E_1\)and \(E_n\)?
    • We can make an effective nuclear charge \(Z_{eff}\) that can be added to this expression such that this equation holds \[E_z = n Z_{eff}^2 E_{Z=1} \].Calculate the \(Z_{eff}\) for each of the ten atoms and plot up \(Z_{eff}\) vs. \(Z\). Explain the origin of this difference (use may want to reference Slater's Rules).

    See Excel sheet


    #7B Homework Solution: ab initio Calculations is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by LibreTexts.

    • Was this article helpful?