{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Solutions Q17 - 30"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"# import all python add-ons etc that will be needed later on\n",
"%matplotlib inline\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"from sympy import *\n",
"init_printing() # allows printing of SymPy results in typeset maths format\n",
"plt.rcParams.update({'font.size': 14}) # set font size for plots"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Q17 answer\n",
"The nitrogen in ammonia is sp$^3$ hybridized with a lone pair of electrons taking up the fourth position; if there were a proton on this position making NH$_4^+$ the point group would be tetrahedral. Two views of ammonia are shown in Fig. 7.78. In the left figure the third mirror plane is in line with the axis and not visible behind the other two.\n",
"\n",
"From the right-hand view it is clear that the molecule has a threefold axis, the molecule can be rotated by $120^\\text{o}$, a $C_3^+$ operation, or by $-120^\\text{o}$ a $C_3^-$ operation, to become indistinguishable. There is no horizontal mirror plane or $C_2$ axes, but there are three vertical mirror planes.\n",
"\n",
"To form the product table, products $C_3^+ \\otimes C_3^-$ and $C_3 \\otimes \\sigma_V$ etc. have to be evaluated this way as well as the other way round. There are two ways at least of doing this; one is geometric, the other algebraic and is to define a basis set and to work out the matrix for each operation and then to multiply them together. To use the geometrical method the H atoms have to be labelled so that they can be distinguished as shown in figure 78.\n",
"\n",
"![Drawing](matrices-fig78.png)\n",
"\n",
"Figure 78 $C_{3V}$ geometry. The $C_3$ axis is also the principal axis.\n",
"______\n",
"\n",
"The effect of $C_3^+ \\otimes C_3^-$ is first to move atom H$_1$ by one place ($120^\\text{o}$) to H$_2$ then to move it back, therefore making the molecule identical $E$; as an equation $C_3^+ \\otimes C_3^- = E$. The reversed multiplication order produces the same result. In two consecutive $C_3$ operations, the first moves H$_1$ by one place and the second does the same producing a total rotation of $240^\\text{o}$ or $C_3^+ \\otimes C_3^+ =C_3^-$. Now consider a rotation and reflection; $C_3^+ \\otimes \\sigma_V$. The rotation moves H$_1$ to H$_2$ and H$_3$ to to H$_1$; reflection now exchanges the H$_3$ and H$_1$ atoms and the effect is the same as a $\\sigma_V'$ reflection: $C_3^+ \\otimes \\sigma_V = \\sigma_V'$. The reverse product is $\\sigma_V \\otimes C_3^+ = \\sigma_V''$ showing that these two operations do not commute; Figure 79.\n",
"\n",
"![Drawing](matrices-fig79.png)\n",
"\n",
"Figure 79 Top: $C_3^+ \\otimes \\sigma_V = \\sigma_V'$ and bottom $\\sigma_V \\otimes C_3^+ = \\sigma_V''$ showing that these operations do not commute.\n",
"______\n",
"\n",
"The complete table is shown below.\n",
"\n",
"$$\\displaystyle \\begin{array}{l|llllll}\n",
"\\hline\n",
" & E & C_3^+ & C_3^- & \\sigma_V & \\sigma_V' & \\sigma_V'' \\\\[2ex]\n",
"\\hline \n",
" \\text{(first) } E & E & C_3^+ & C_3^- & \\sigma_V & \\sigma_V' & \\sigma_V'' \\\\\n",
" C_3^+ & C_3^+ & C_3^- & E & \\sigma_V' & \\sigma_V'' &\\sigma_V \\\\\n",
" C_3^- & C_3^- & E & C_3^+ & \\sigma_V'' & \\sigma_V & \\sigma_V' \\\\\n",
" \\sigma_V & \\sigma_V & \\sigma_V'' & \\sigma_V' & E & C_3^- & C_3^+\\\\\n",
" \\sigma_V' & \\sigma_V' & \\sigma_V & \\sigma_V'' & C_3^+ & E & C_3^-\\\\\n",
" \\sigma_V'' & \\sigma_V'' & \\sigma_V' & \\sigma_V & C_3^- & C_3^+ &E\\\\\n",
" \\hline\n",
" \\end{array}$$\n",
"\n",
"\n",
"## Q18 answer\n",
"Making a sketch of ammonia, the bonds can be labelled $a, b, c$ to represent vectors along them. A $C_3^+$ rotation has the effect $abc \\to cab$ which as a matrix equation is\n",
"\n",
"$$\\displaystyle \\begin{bmatrix} c \\\\ a \\\\b \\end{bmatrix}= C_3^+ \\begin{bmatrix} a \\\\ b \\\\c \\end{bmatrix}= \\begin{bmatrix} 0 & 0 & 1 \\\\ 1 & 0 & 0 \\\\0 & 1 & 0 \\end{bmatrix}\\begin{bmatrix} a \\\\ b \\\\c \\end{bmatrix}$$\n",
"\n",
"![Drawing](matrices-figQ18.png)\n",
"\n",
"Other matrices are\n",
"\n",
"$$\\displaystyle \n",
" C_3^- \\equiv \\begin{bmatrix} 0 & 1 & 0 \\\\ 0 & 0 & 1 \\\\1 & 0 & 0 \\end{bmatrix} ,\\; \n",
" \\sigma_V \\equiv \\begin{bmatrix} 0 & 1 & 0 \\\\ 1 & 0 & 0 \\\\0 & 0 & 1 \\end{bmatrix} ,\\;\n",
" \\sigma_V' \\equiv \\begin{bmatrix} 1 & 0 & 0 \\\\ 0 & 0 & 1 \\\\0 & 1 & 0 \\end{bmatrix} ,\\;\n",
"\\sigma_V'' \\equiv \\begin{bmatrix} 0 & 0 & 1 \\\\ 0 & 1 & 0 \\\\1 & 0 & 0 \\end{bmatrix} ,\\;\n",
" E \\equiv \\begin{bmatrix} 1 & 0 & 0 \\\\ 0 & 1 & 0 \\\\0 & 0 & 1 \\end{bmatrix} $$\n",
"\n",
"A simple way to check the table is to use Sympy to multiply out the matrices and then check the product. The product always has to be member of the group so it is only necessary to check which one this is via a look up table of names."
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" E C3p C3m sigV sig_V' sigV'' \n",
"----------------------------------------------------------\n",
"\n",
"E | E C3p C3m sigV sig_V' sigV'' \n",
"\n",
"C3p | C3p C3m E sigV'' sigV sig_V'\n",
"\n",
"C3m | C3m E C3p sig_V' sigV'' sigV \n",
"\n",
"sigV | sigV sig_V' sigV'' E C3p C3m \n",
"\n",
"sig_V' | sig_V' sigV'' sigV C3m E C3p \n",
"\n",
"sigV'' | sigV'' sigV sig_V' C3p C3m E \n",
"\n"
]
}
],
"source": [
"# Algorithm Matrix Products\n",
"E, C3p, C3m, sigV, sigVd, sigVdd = symbols('E, C3p, C3m, sigV, sigVd, sigVdd' )\n",
"\n",
"E = Matrix([[1,0,0],[0,1,0],[0,0,1]])\n",
"C3p = Matrix([[0,0,1],[1,0,0],[0,1,0]]) # C3 plus \n",
"C3m = Matrix([[0,1,0],[0,0,1],[1,0,0]]) # C3 minus\n",
"sigV = Matrix([[0,1,0],[1,0,0],[0,0,1]])\n",
"sigVd = Matrix([[1,0,0],[0,0,1],[0,1,0]]) # sigma V'\n",
"sigVdd= Matrix([[0,0,1],[0,1,0],[1,0,0]])\n",
"\n",
"alist = [ E, C3p, C3m, sigV, sigVd, sigVdd ]\n",
"n = len(alist)\n",
"nlist=['E', 'C3p', 'C3m','sigV',\"sig_V'\",\"sigV'' \"]\n",
"print('{:12s}'.format(' '),end='')\n",
"print(*('{:6s}'.format(k) for i,k in enumerate(nlist)) )\n",
"print('----------------------------------------------------------\\n')\n",
"for i in range(n):\n",
" print('{:8s} |'.format(nlist[i]),end='')\n",
" for j in range(n):\n",
" prod = alist[i]*alist[j]\n",
" print(' {:6s}'.format( nlist[alist.index(prod) ] ) ,end='' ) # find index of match\n",
" print('\\n')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Q19 answer\n",
"Multiplying gives\n",
"\n",
"$$\\displaystyle \\begin{array}{l|lll}\n",
" & E & B & C \\\\\n",
"\\hline\n",
"E & E & B & C\\\\\n",
"B & B & BB & BC\\\\\n",
"C & C & CB & CC \\end{array}$$\n",
"\n",
"and $BB$ can only be $E$ or $C$. Suppose that $BB = E$ then the table is\n",
"\n",
"$$\\displaystyle \\begin{array}{l|lll}\n",
" & E & B & C \\\\\n",
"\\hline\n",
"E & E & B & C\\\\\n",
"B & B & E & C\\\\\n",
"C & C & C & CC \\end{array}$$\n",
"\n",
"with $BC = C$, which is not correct no matter what $CC$ is because $C$ appears twice in the bottom row and last column. Trying $BB = C$ produces the correct table.\n",
"\n",
"$$\\displaystyle \\begin{array}{l|lll}\n",
" & E & B & C \\\\\n",
"\\hline\n",
"E & E & B & C\\\\\n",
"B & B & C & E\\\\\n",
"C & C & E & B \\end{array}$$\n",
"\n",
"## Q20 answer\n",
"Using the matrices from the previous question the products are"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$\\displaystyle \\left[\\begin{matrix}1 & 0 & 0\\\\0 & 0 & 1\\\\0 & 1 & 0\\end{matrix}\\right]$"
],
"text/plain": [
"⎡1 0 0⎤\n",
"⎢ ⎥\n",
"⎢0 0 1⎥\n",
"⎢ ⎥\n",
"⎣0 1 0⎦"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"A1 = C3m.inv()* sigV*C3m\n",
"A1"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$\\displaystyle \\left[\\begin{matrix}0 & 0 & 1\\\\0 & 1 & 0\\\\1 & 0 & 0\\end{matrix}\\right]$"
],
"text/plain": [
"⎡0 0 1⎤\n",
"⎢ ⎥\n",
"⎢0 1 0⎥\n",
"⎢ ⎥\n",
"⎣1 0 0⎦"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"A2 = C3p.inv()*sigV*C3p\n",
"A2"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Looking at the matrices it follows that $A_1 = \\sigma_V'$ and $A_2 = \\sigma_V''$ or $\\sigma_V' = \\sigma_V^{-1}C_3^-\\sigma_V$ and $\\sigma_V'' = \\sigma_V^{-1}C_3^+\\sigma_V$.\n",
"\n",
"## Q21 answer\n",
"The highest symmetry axis, the principal axis,is a $90^\\text{o}$ rotation axis about the Xe atom, labelled $C_4$, and is perpendicular to the plane of the molecule. There is also a $C_2$ axis, a $180^\\text{o}$ rotation about the same axis. The molecule has a centre of inversion.\n",
"\n",
"There is a $C_2$ axis through each diagonal along the F-Xe-F bonds of the square and similarly bisecting adjacent F atoms. There is a horizontal mirror plane in the plane of the molecule; horizontal because it is perpendicular to the principal axis. This immediately suggests a $D$ point group. There are dihedral mirror planes bisecting adjacent F atoms as well as vertical mirror planes along the F-Xe-F atoms.\n",
"\n",
"The Xe atom is at the centre of inversion. (An inversion centre is present if any position on one side of an object can be passed through a single point in the object to an equivalent position on the other side.) Finally, there are two $90^\\text{o}$ rotation reflection operations labelled $S_4$. The fourfold rotation axis, the horizontal mirror plane and the inversion centre are sufficient to identify $D_{4h}$ as the point group.\n",
"\n",
"![Drawing](matrices-fig80.png)\n",
"\n",
"Figure 80. XeF$_4$ and some of its rotation and reflection symmetry elements.\n",
"______\n",
"\n",
"## Q22 answer\n",
"Looking at figure 11 the reflection in the mirror plane labelled $\\sigma_V'$ leaves the atoms unchanged. The symbolic equation is therefore,\n",
"\n",
"$$\\displaystyle [\\sigma']=\\begin{bmatrix} Cl \\\\O_a \\\\ O_b \\end{bmatrix}=\\begin{bmatrix} Cl \\\\O_a \\\\ O_b \\end{bmatrix}$$\n",
"\n",
"As the atoms are unchanged, the matrix should be a unit matrix,\n",
"\n",
"$$\\displaystyle [\\sigma']=\\begin{bmatrix} 1 & 0 & 0 \\\\0 & 1 & 0 \\\\ 0 & 0 & 1 \\end{bmatrix}$$\n",
"\n",
"The reflection in the other plane $\\sigma_V$ swaps atoms $a$ and $b$ therefore the generic equation is\n",
"\n",
"$$\\displaystyle [\\sigma]=\\begin{bmatrix} Cl \\\\O_a \\\\ O_b \\end{bmatrix}=\\begin{bmatrix} Cl \\\\O_b \\\\ O_a \\end{bmatrix}$$\n",
"\n",
"which has the same effect as a $C_2$ rotation so that the matrix is\n",
"\n",
"$$\\displaystyle [\\sigma]=\\begin{bmatrix} 1 & 0 & 0 \\\\0 & 0 & 1 \\\\ 0 & 1 & 0 \\end{bmatrix}$$\n",
"\n",
"## Q23 answer\n",
"(a) The rotation equation is \n",
"\n",
"$$\\displaystyle [C_2]\\begin{bmatrix} Cl \\\\O_a \\\\ O_b \\end{bmatrix}=\\begin{bmatrix} 1 & 0 & 0 \\\\0 & 0 & 1 \\\\ 0 & 1 & 0 \\end{bmatrix}\\begin{bmatrix} Cl \\\\O_a \\\\ O_b \\end{bmatrix}=\\begin{bmatrix} Cl \\\\O_b \\\\ O_a \\end{bmatrix}$$\n",
"\n",
"following this by a reflection $\\sigma_V'$ the equation is\n",
"\n",
"$$\\displaystyle [\\sigma_V'][C_2]\\begin{bmatrix} Cl \\\\O_a \\\\ O_b \\end{bmatrix}=\\begin{bmatrix} 1 & 0 & 0 \\\\0 & 1 & 0 \\\\ 0 & 0 & 1 \\end{bmatrix}\\begin{bmatrix} 1 & 0 & 0 \\\\0 & 0 & 1 \\\\ 0 & 1 & 0 \\end{bmatrix}\\begin{bmatrix} Cl \\\\O_a \\\\ O_b \\end{bmatrix}=\\begin{bmatrix} Cl \\\\O_b \\\\ O_a \\end{bmatrix}$$\n",
"\n",
"which shows that these two operations have had the same effect as reflection $\\sigma_V$ because the unit diagonal matrix, $\\sigma_V'$ is the same as the identity matrix and changes nothing:\n",
"\n",
"$$\\displaystyle \\begin{bmatrix} 1 & 0 & 0 \\\\0 & 1 & 0 \\\\ 0 & 0 & 1 \\end{bmatrix}\\begin{bmatrix} 1 & 0 & 0 \\\\0 & 0 & 1 \\\\ 0 & 1 & 0 \\end{bmatrix}=\\begin{bmatrix} 1 & 0 & 0 \\\\0 & 0 & 1 \\\\ 0 & 1 & 0 \\end{bmatrix}\\equiv [\\sigma]$$\n",
"\n",
"so the two operations are equivalent to a reflection $[\\sigma_V'][C_2]=[\\sigma_V]$\n",
"\n",
"Reversing the order of operations changes nothing in this case because $[\\sigma_V]$ is a unit diagonal matrix so that $[C_2][\\sigma_V']=[\\sigma_V]$\n",
"\n",
"(b) The fact that $[\\sigma'][C_2]=[\\sigma]$ and $[C_2][\\sigma']=[\\sigma]$ means that $[C_2,\\sigma']=0$ and therefore these two commute.\n",
"\n",
"## Q24 answer\n",
"The carbon atoms are on the principal axis so remain unchanged by any operation. The $C_2$ operation swaps $H_a$ with $H)b$ and $F_a$ with $F_b$ and leaves the two carbon atoms unchanged. The symbolic matrix is\n",
"\n",
"$$\\displaystyle [C_2]\\begin{bmatrix} H_a \\\\H_b \\\\ C \\\\C \\\\F_a\\\\F_b \\end{bmatrix}= \\begin{bmatrix} H_b \\\\H_a \\\\ C \\\\C \\\\F_b\\\\F_a \\end{bmatrix}$$\n",
"and following the previous examples we can try\n",
"\n",
"$$\\displaystyle [C_2]=\\begin{bmatrix} 0& 1 & 0 & 0 &0 & 0 \\\\ 1& 0 & 0 & 0 & 0 & 0 \\\\ 0& 0 & 1 & 0 &0&0 \\\\0& 0 & 0 & 1 &0&0 \\\\ 0& 0 & 0 & 0 &0&1 \\\\ 0& 0 & 0 & 0 &1&0 \\end{bmatrix}$$\n",
"\n",
"To check the calculation with Sympy, the matrix is formed containing only zeros and then the six non-zero elements are then added."
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"C2, M, Ha, Hb, C, Fa, Fb = symbols('C2, M, Ha, Hb, C, Fa, Fb')\n",
"C2 = zeros(6,6)\n",
"C2[0,1] = 1\n",
"C2[1,0] = 1\n",
"C2[2,2] = 1\n",
"C2[3,3] = 1\n",
"C2[4,5] = 1\n",
"C2[5,4] = 1"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$\\displaystyle \\left[\\begin{matrix}0 & 1 & 0 & 0 & 0 & 0\\\\1 & 0 & 0 & 0 & 0 & 0\\\\0 & 0 & 1 & 0 & 0 & 0\\\\0 & 0 & 0 & 1 & 0 & 0\\\\0 & 0 & 0 & 0 & 0 & 1\\\\0 & 0 & 0 & 0 & 1 & 0\\end{matrix}\\right]$"
],
"text/plain": [
"⎡0 1 0 0 0 0⎤\n",
"⎢ ⎥\n",
"⎢1 0 0 0 0 0⎥\n",
"⎢ ⎥\n",
"⎢0 0 1 0 0 0⎥\n",
"⎢ ⎥\n",
"⎢0 0 0 1 0 0⎥\n",
"⎢ ⎥\n",
"⎢0 0 0 0 0 1⎥\n",
"⎢ ⎥\n",
"⎣0 0 0 0 1 0⎦"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"C2"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$\\displaystyle \\left[\\begin{matrix}Hb\\\\Ha\\\\C\\\\C\\\\Fb\\\\Fa\\end{matrix}\\right]$"
],
"text/plain": [
"⎡Hb⎤\n",
"⎢ ⎥\n",
"⎢Ha⎥\n",
"⎢ ⎥\n",
"⎢C ⎥\n",
"⎢ ⎥\n",
"⎢C ⎥\n",
"⎢ ⎥\n",
"⎢Fb⎥\n",
"⎢ ⎥\n",
"⎣Fa⎦"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"M = Matrix([Ha, Hb, C, C, Fa, Fb] )\n",
"C2*M"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"If the ordering in the basis set is changed, say to $[Fb,C,Ha,C,Fa,Hb]$ and _provided the matrix is similarly changed_ then the same result is obtained.\n",
"\n",
"## Q25 answer\n",
"The structure is shown below; the principal axis passes through the N atoms and equally distant from each of the H atoms.\n",
"\n",
"![Drawing](matrices-fig81.png)\n",
"\n",
"Figure 81. Ammonia with $+120^\\text{o}$ rotation about the principal axis; H atoms $1\\to 2,2\\to 3,3\\to 1$ and one of three vertical mirror planes $\\sigma$. In this plan H atom $3$ is unchanged, but $1\\to 2$ and $2\\to 1$.\n",
"______\n",
"Rotation by $120^\\text{o}$ makes the molecule indistinguishable from its original shape. There are two rotation operations; clockwise rotation by $120^\\text{o}$ ($C_3$), and by $240^\\text{o}$, which is the same as two rotations of $120^\\text{o}$ ($C_3^2$), or one rotation anticlockwise. Each leaves the molecule indistinguishable. Rotation by $+120^\\text{o}$ about the principal axis leaves the N atom alone and swaps the H atoms about in the order $a \\to b,\\;b \\to c$ and $c \\to a$. The matrix equation is\n",
"\n",
"$$\\displaystyle [C_3]\\begin{bmatrix} N \\\\H_1 \\\\ H_2 \\\\H_3 \\end{bmatrix}= \\begin{bmatrix} N \\\\H_2 \\\\ H_3 \\\\H_1 \\end{bmatrix}$$\n",
"\n",
"The matrix to do this is $\\displaystyle C_3=\\begin{bmatrix} 1 & 0 & 0 &0 \\\\ 0 & 0 & 1 & 0 \\\\ 0 & 0 & 0 & 1 \\\\ 0& 1 & 0 & 0\\\\ \\end{bmatrix}$ which can be confirmed by multiplying it out.\n",
"\n",
"The matrix for the reflection shown above is simpler as only the H atoms change place $\\displaystyle \\sigma_V=\\begin{bmatrix} 1 & 0 & 0 &0 \\\\ 0 & 0 & 1 & 0 \\\\ 0 & 1 & 0 & 0 \\\\ 0& 0 & 0 & 1\\\\ \\end{bmatrix}$.\n",
"\n",
"## Q26 answer\n",
"(a) Rotation about the principal axis is by $2\\pi/n$ for a $C_n$ axis then the matrix is (see text)\n",
"\n",
"$$\\displaystyle \\sigma_V=\\begin{bmatrix} \\cos(2\\pi/n) & \\sin(2\\pi/n) & 0 \\\\ -\\sin(2\\pi/n) & \\cos(2\\pi/n) & 0 \\\\ 0 & 0 & 1 \\end{bmatrix}$$\n",
"\n",
"Reflection in an axis perpendicular to the principal or z-axis inverts the $z$ coordinates and leaves $x$ and $y$ unchanged. The matrix is therefore $\\displaystyle \\begin{bmatrix} 1 & 0 & 0 \\\\ 0 & 1 & 0 \\\\ 0 & 0 & -1 \\end{bmatrix}$ and multiplying these forms the combined operation giving \n",
"\n",
"$$\\displaystyle \\sigma_V=\\begin{bmatrix} \\cos(2\\pi/n) & \\sin(2\\pi/n) & 0 \\\\ -\\sin(2\\pi/n) & \\cos(2\\pi/n) & 0 \\\\ 0 & 0 & -1 \\end{bmatrix}$$\n",
"\n",
"(b) The trace is $T_{S_n}=2\\cos(2\\pi/n)-1$\n",
"\n",
"## Q27 answer\n",
"(a) The principal axis is taken as the axis of highest symmetry, which is the threefold axis through the centre and perpendicular to the plane of the molecule. The plane of the molecule is therefore a horizontal mirror plane and there are three, twofold axes perpendicular to the principal axis as well as three vertical mirror planes. These elements are shown in the figure. Only one of the three vertical planes is shown for clarity.\n",
"\n",
"![Drawing](matrices-fig82.png)\n",
"\n",
"Figure 82. Symmetry elements in a planar molecule belonging to the $D_{3h}$ point group.\n",
"_______\n",
"\n",
"(b) Considering just the three p orbitals produces $C_3$ symmetry. The horizontal mirror plane and the $C_2$ axes perpendicular to the principal axis are lost from $D_{3h}$ because of the phase of the orbitals.\n",
"\n",
"Operating on the orbitals with the C operators moves the location of each orbital so according to the rules these count zero. The $C_3^2$ operation moves the orbital twice by $120^\\text{o}$ or once by $-120^\\text{o}$ from its initial position. The identity counts three therefore the reducible representation is\n",
"\n",
"$$\\displaystyle \\begin{array}{c|ccc} \n",
" & E & C_3 & C_3^2 \\\\\n",
"\\hline\n",
"\\Gamma_R & 3 & 0 & 0 \\\\\n",
"\\hline \n",
"\\end{array}$$\n",
"\n",
"Using the tabular method this is reduced as\n",
"\n",
"$$\\displaystyle \\begin{array}{c|ccc} \n",
"C_3 & E & C_3 & C_3^2 & (\\sum )/h \\\\\n",
"\\hline\n",
"\\Gamma_R & 3 & 0 & 0 \\\\\n",
"\\hline\n",
"A & 3 & 0 & 0 & 3/3=1 \\\\\n",
"E & 3 & 0 & 0 & 3/3=1 \\\\\n",
" & 3 & 0 & 0 & 3/3=1 \\\\\n",
"\\hline \n",
"\\end{array}$$\n",
"\n",
"which means that the orbitals contain the $A$ and $E$ symmetry species. The normalized molecular orbitals are therefore constructed using the point group character table\n",
"\n",
"$$\\displaystyle \\varphi_A= \\frac{1}{\\sqrt{3}}(\\psi_{p_1} +\\psi_{p_2}+\\psi_{p_3} ) \\\\\n",
"\\varphi_{E_1}= \\frac{1}{\\sqrt{3}}(\\psi_{p_1} +e^{-2\\pi i/3}\\psi_{p_2}+e^{2\\pi i/3}\\psi_{p_3} ) \\\\\n",
"\\varphi_{E_2}= \\frac{1}{\\sqrt{3}}(\\psi_{p_1} +e^{2\\pi i/3}\\psi_{p_2}+e^{-2\\pi i/3}\\psi_{p_3} ) $$\n",
"\n",
"which contain complex numbers this is, however, not a problem. Note that any measurement involves an expectation value of the form $\\langle A \\rangle = \\int \\varphi^* A \\varphi d\\tau$ if the operator is $A$ and because of the complex conjugate the result is always a real number. However, to avoid having a complex number in the wavefunction a linear combination can be made such as $\\varphi_{E_1} +\\varphi_{E_2}$ and $i(\\varphi_{E_1} -\\varphi_{E_2})$ both of which will be real. The identity $e^{i\\theta}+e^{-i\\theta}=2$ and $e^{i\\theta}-e^{-i\\theta}=2i\\sin(\\theta)$ can be used. As the angle is $\\theta = 120$ the orbitals become\n",
"\n",
"$$\\displaystyle \\varphi_+ =\\frac{1}{\\sqrt{2}}\\left(\\varphi_{E_1}+\\varphi_{E_2}\\right)=\\frac{1}{\\sqrt{6}}\\left(2\\psi_{p_1}-\\psi_{p_2}-\\psi_{p_3}\\right) \\\\ \\varphi_- =i\\frac{1}{\\sqrt{2}}\\left(\\varphi_{E_1}-\\varphi_{E_2}\\right) =\\frac{1}{\\sqrt{2}}\\left(\\psi_{p_2}-\\psi_{p_3}\\right)$$\n",
"\n",
"which can easily be visualised.\n",
"\n",
"## Q28 answer\n",
"(a) The point group is $C_{2V}$. The mirror plane of the carbon atoms is lost because the orbital are $\\pm$ or black/white above and below this plane. The principal and z-axis is perpendicular to the plane of the carbons.\n",
"\n",
"(b) There are ten orbitals in the basis set;therefore,the character of $E = 10$. On rotation about the principal axis all the atoms are moved so the character of $C_2 = 0$. On reflection about the $x$ axis (carbons 9/10) eight orbitals are changed each of which counts zero and two are unchanged making the character of $\\sigma(xz) = 2$. The other reflection plane swaps all atoms so counts zero. The reducible representation is therefore\n",
"\n",
"$$\\displaystyle \\begin{array}{c|rrrr} \n",
" & E & C_2 & \\sigma(x,z) & \\sigma(y,z) &(\\sum )/4 \\\\\n",
"\\hline\n",
"\\Gamma_R & 10 & 0 & 2 & 2 \\\\\n",
"\\hline\n",
"A_1 & 10 & 0 & 2 & 0 & 12/4=3 \\\\\n",
"A_2 & 10 & 0 & -2 & 0 & 8/4=2 \\\\\n",
"B_1 & 10 & 0 & 2 & 0 & 3 \\\\\n",
"B_2 & 10 & 0 & -2 & 0 & 2 \\\\\n",
"\\hline \n",
"\\end{array}$$\n",
"\n",
"The reduced representation contains still contains ten orbitals but organized as $3A_1 + 2A_2 + 3B_1 + 2B_2$symmetry species.\n",
"\n",
"\n",
"## Q29 answer\n",
"Imagine vectors representing p orbitals with different phase above and below the plane of the molecule. Operate on these to form the reduced representation. (Recall that atoms moved count $0$, atoms unchanged $1$ and inverted $-1$.)\n",
"\n",
"$$\\displaystyle \\begin{array}{c|cccccccc}\n",
"D_{2h} & E & C_2(z) & C_2(y) & C_2(x) & i & \\sigma(xy) & \\sigma(xz) & \\sigma(yz) \\\\\n",
"\\hline \n",
"\\Gamma_R & 10 & 0 & -2 & 0 & 0 & 2 & 0 & -10 \\\\ \n",
"\\hline \n",
"\\end{array}$$\n",
"\n",
"The number of operations is $h = 8$ and the reduction table is shown below with columns that are zero deleted.\n",
"\n",
"$$\\displaystyle \\begin{array}{c|cccccccc}\n",
"D_{2h} & E & C_2(y) & \\sigma(xy) & \\sigma(yz) & (\\sum)/8\\\\\n",
"\\hline \n",
"A_g & 10 & -2 & 2 & -10 & 0\\\\ \n",
"B_{1g} & 10 & 2 & 2 & 10 & 24/8=3 \\\\\n",
"B_{2g} & 10 & -2 & -2 & 10 & 2 \\\\\n",
"B_{3g} & 10 & 2 & -2 & -10 & 0 \\\\\n",
"A_{u} & 10 & -2 & -2 & 10 & 2 \\\\\n",
"B_{1u} & 10 & 2 & -2 & -10 & 0 \\\\\n",
"B_{2u} & 10 & -2 & 2 & -10 & 0 \\\\\n",
"B_{3u} & 10 & 2 & 2 & 10 & 3 \\\\\n",
"\\hline\n",
"\\end{array}$$\n",
"\n",
"The reduced representation therefore contains $3B_{1g} + 2B_{2g} + 2A_u + 3B_{3u}$ which are $10$ in total, the same as the number of orbitals.\n",
"\n",
"Choosing atom 1 as victim, the $L_M$ formula Section 6.16 produces the following table. The row for the totally symmetric species $A_g$ shows that the sum is zero, no $A_g$ species is present in the reduced representation. Three rows are included for $B_{1g}$, starting with orbitals on atoms 1, 2, and 9. The remaining rows for $B_{2g}$ and $B_{3u}$ are missing and are left for you to work out. Note that if the $C_{2V}$ point group was used instead then the same results would have been obtained, and with less effort, showing that the $\\pi$ orbital symmetry really is important.\n",
"\n",
"$$\\displaystyle \\begin{array}{c|cccccccc}\n",
"D_{2h} & E & C_2(z) & C_2(y) & C_2(x) & i & \\sigma(xy) & \\sigma(xz) & \\sigma(yz) & L_M \\\\\n",
"\\hline \n",
"L_{Ag} & v_1 & -v_4 & -v_8 & v_5 & -v_5 & v_8 & v_4 & -v_1 & 0 \\\\ \n",
"\\hline \n",
"L_{B1g} & v_1 & -v_4 & v_8 & -v_5 & -v_5 & v_8 & -v_4 & v_1 & 2(v_1-v_4-v_5+v_8) \\\\\n",
"L_{B1g} & v_2 & -v_3 & v_7 & -v_6 & -v_6 & v_7 & -v_3 & v_2 & 2(v_2-v_3-v_6+v_7) \\\\\n",
"L_{B1g} & v_9 & -v_{10} & v_9 & -v_{10} & -v_{10} & v_9 & -v_{10} & v_9 & 4(v_9-v_{10}) \\\\\n",
"\\hline \n",
"L_{B2g} & v_1 & v_4 & -v_8 & -v_5 & v_5 & -v_8 & v_4 & v_1 & 2(v_1+v_4-v_8) \\\\\n",
"\\hline \n",
"L_{Au} & v_1 & -v_4 & -v_8 & v_5 & v_5 & -v_8 & -v_4 & v_1 & 2(v_1-v_4+v_5-v_8) \\\\\n",
"L_{Au} & v_2 & -v_3 & -v_7 & v_6 & v_6 & -v_7 & -v_3 & v_2 & 2(v_2-v_3+v_6-v_7) \\\\\n",
"\\hline\n",
"\\end{array}$$\n",
"\n",
"The degenerate symmetry species can form linear combinations with one another if they are normalized and orthogonal. If not orthogonal, they can be made so using the Gram - Schmidt method. The energy can be found by calculating the expectation values $\\langle \\psi_a|H|\\psi_a\\rangle,\\; \\langle \\psi_a|H|\\psi_b\\rangle;\\;\\langle \\psi_a|H|\\psi_c\\rangle;$ , etc. where $\\psi_{a,b,c}$ are different symmetry adapted orbitals, for example the three normalized and orthogonal $B_{1g}$ orbitals are $\\psi_a =(v_2 -v_3 -v_6 +v_7 )/2,\\;\n",
"\\psi_b =(v_1 -v_4 -v_5 +v_8)/2 $ and $\\psi_c =(v_9 -v_{10} )/\\sqrt{2}$.\n",
"\n",
"The calculation is shown below. Note that a vector such as $v_2-v_3+v_6-v_7$ when used in a programme has to be written with each term with a separate index in a 1D matrix or list, i.e. as $[v_2,-v_3,+v_6,-v_7]$. As each $v$ vector is real and has unit length then ${v^*\\cdot v}=v^2=1$ and this simplifies the normalisation terms."
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [],
"source": [
"v1,v2,v3,v4,v5,v6,v7,v8,v9,v10 = symbols('v1,v2,v3,v4,v5,v6,v7,v8,v9,v10',positive=True)\n",
"\n",
"psi1 = Matrix([2*v1,0,0,2*v4,-2*v5,0,0,+2*v8,0,0]) # each length 10\n",
"psi2 = Matrix([0,2*v2,-2*v3,0,0,-2*v6,+2*v7,0,0,0])\n",
"psi3 = Matrix([0,0,0,0,0,0,0,0,4*v9,-4*v10])"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$\\displaystyle \\left[ \\left[\\begin{matrix}\\frac{2 v_{1}}{\\sqrt{4 v_{1}^{2} + 4 v_{4}^{2} + 4 v_{5}^{2} + 4 v_{8}^{2}}}\\\\0\\\\0\\\\\\frac{2 v_{4}}{\\sqrt{4 v_{1}^{2} + 4 v_{4}^{2} + 4 v_{5}^{2} + 4 v_{8}^{2}}}\\\\- \\frac{2 v_{5}}{\\sqrt{4 v_{1}^{2} + 4 v_{4}^{2} + 4 v_{5}^{2} + 4 v_{8}^{2}}}\\\\0\\\\0\\\\\\frac{2 v_{8}}{\\sqrt{4 v_{1}^{2} + 4 v_{4}^{2} + 4 v_{5}^{2} + 4 v_{8}^{2}}}\\\\0\\\\0\\end{matrix}\\right], \\ \\left[\\begin{matrix}0\\\\\\frac{2 v_{2}}{\\sqrt{4 v_{2}^{2} + 4 v_{3}^{2} + 4 v_{6}^{2} + 4 v_{7}^{2}}}\\\\- \\frac{2 v_{3}}{\\sqrt{4 v_{2}^{2} + 4 v_{3}^{2} + 4 v_{6}^{2} + 4 v_{7}^{2}}}\\\\0\\\\0\\\\- \\frac{2 v_{6}}{\\sqrt{4 v_{2}^{2} + 4 v_{3}^{2} + 4 v_{6}^{2} + 4 v_{7}^{2}}}\\\\\\frac{2 v_{7}}{\\sqrt{4 v_{2}^{2} + 4 v_{3}^{2} + 4 v_{6}^{2} + 4 v_{7}^{2}}}\\\\0\\\\0\\\\0\\end{matrix}\\right], \\ \\left[\\begin{matrix}0\\\\0\\\\0\\\\0\\\\0\\\\0\\\\0\\\\0\\\\\\frac{4 v_{9}}{\\sqrt{16 v_{10}^{2} + 16 v_{9}^{2}}}\\\\- \\frac{4 v_{10}}{\\sqrt{16 v_{10}^{2} + 16 v_{9}^{2}}}\\end{matrix}\\right]\\right]$"
],
"text/plain": [
"⎡⎡ 2⋅v₁ ⎤ ⎡ 0 ⎤ \n",
"⎢⎢──────────────────────────────────⎥ ⎢ ⎥ \n",
"⎢⎢ _______________________________⎥ ⎢ 2⋅v₂ ⎥ \n",
"⎢⎢ ╱ 2 2 2 2 ⎥ ⎢──────────────────────────────────⎥ \n",
"⎢⎢╲╱ 4⋅v₁ + 4⋅v₄ + 4⋅v₅ + 4⋅v₈ ⎥ ⎢ _______________________________⎥ ⎡\n",
"⎢⎢ ⎥ ⎢ ╱ 2 2 2 2 ⎥ ⎢\n",
"⎢⎢ 0 ⎥ ⎢╲╱ 4⋅v₂ + 4⋅v₃ + 4⋅v₆ + 4⋅v₇ ⎥ ⎢\n",
"⎢⎢ ⎥ ⎢ ⎥ ⎢\n",
"⎢⎢ 0 ⎥ ⎢ -2⋅v₃ ⎥ ⎢\n",
"⎢⎢ ⎥ ⎢──────────────────────────────────⎥ ⎢\n",
"⎢⎢ 2⋅v₄ ⎥ ⎢ _______________________________⎥ ⎢\n",
"⎢⎢──────────────────────────────────⎥ ⎢ ╱ 2 2 2 2 ⎥ ⎢\n",
"⎢⎢ _______________________________⎥ ⎢╲╱ 4⋅v₂ + 4⋅v₃ + 4⋅v₆ + 4⋅v₇ ⎥ ⎢\n",
"⎢⎢ ╱ 2 2 2 2 ⎥ ⎢ ⎥ ⎢\n",
"⎢⎢╲╱ 4⋅v₁ + 4⋅v₄ + 4⋅v₅ + 4⋅v₈ ⎥ ⎢ 0 ⎥ ⎢\n",
"⎢⎢ ⎥ ⎢ ⎥ ⎢\n",
"⎢⎢ -2⋅v₅ ⎥ ⎢ 0 ⎥ ⎢\n",
"⎢⎢──────────────────────────────────⎥, ⎢ ⎥, ⎢\n",
"⎢⎢ _______________________________⎥ ⎢ -2⋅v₆ ⎥ ⎢\n",
"⎢⎢ ╱ 2 2 2 2 ⎥ ⎢──────────────────────────────────⎥ ⎢\n",
"⎢⎢╲╱ 4⋅v₁ + 4⋅v₄ + 4⋅v₅ + 4⋅v₈ ⎥ ⎢ _______________________________⎥ ⎢\n",
"⎢⎢ ⎥ ⎢ ╱ 2 2 2 2 ⎥ ⎢\n",
"⎢⎢ 0 ⎥ ⎢╲╱ 4⋅v₂ + 4⋅v₃ + 4⋅v₆ + 4⋅v₇ ⎥ ⎢\n",
"⎢⎢ ⎥ ⎢ ⎥ ⎢\n",
"⎢⎢ 0 ⎥ ⎢ 2⋅v₇ ⎥ ⎢\n",
"⎢⎢ ⎥ ⎢──────────────────────────────────⎥ ⎢\n",
"⎢⎢ 2⋅v₈ ⎥ ⎢ _______________________________⎥ ⎢\n",
"⎢⎢──────────────────────────────────⎥ ⎢ ╱ 2 2 2 2 ⎥ ⎢\n",
"⎢⎢ _______________________________⎥ ⎢╲╱ 4⋅v₂ + 4⋅v₃ + 4⋅v₆ + 4⋅v₇ ⎥ ⎢\n",
"⎢⎢ ╱ 2 2 2 2 ⎥ ⎢ ⎥ ⎢\n",
"⎢⎢╲╱ 4⋅v₁ + 4⋅v₄ + 4⋅v₅ + 4⋅v₈ ⎥ ⎢ 0 ⎥ ⎣\n",
"⎢⎢ ⎥ ⎢ ⎥ \n",
"⎢⎢ 0 ⎥ ⎢ 0 ⎥ \n",
"⎢⎢ ⎥ ⎢ ⎥ \n",
"⎣⎣ 0 ⎦ ⎣ 0 ⎦ \n",
"\n",
" ⎤\n",
" ⎥\n",
" ⎥\n",
" ⎥\n",
" 0 ⎤⎥\n",
" ⎥⎥\n",
" 0 ⎥⎥\n",
" ⎥⎥\n",
" 0 ⎥⎥\n",
" ⎥⎥\n",
" 0 ⎥⎥\n",
" ⎥⎥\n",
" 0 ⎥⎥\n",
" ⎥⎥\n",
" 0 ⎥⎥\n",
" ⎥⎥\n",
" 0 ⎥⎥\n",
" ⎥⎥\n",
" 0 ⎥⎥\n",
" ⎥⎥\n",
" 4⋅v₉ ⎥⎥\n",
"─────────────────────⎥⎥\n",
" __________________⎥⎥\n",
" ╱ 2 2 ⎥⎥\n",
"╲╱ 16⋅v₁₀ + 16⋅v₉ ⎥⎥\n",
" ⎥⎥\n",
" -4⋅v₁₀ ⎥⎥\n",
"─────────────────────⎥⎥\n",
" __________________⎥⎥\n",
" ╱ 2 2 ⎥⎥\n",
"╲╱ 16⋅v₁₀ + 16⋅v₉ ⎦⎥\n",
" ⎥\n",
" ⎥\n",
" ⎥\n",
" ⎦"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"GramSchmidt([psi1,psi2,psi3],True) # orthogonalise and normalise"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The expectation values are the product of the terms including the wavefunctions. Remember that these are really ten elements long with the coefficients for the missing ts set to zero. If a wavefunction is written as \n",
"\n",
"$$\\displaystyle \\psi = \\sum_{k=1}^n a_kv_k$$\n",
"\n",
"where in this example $n=10$. The for $\\psi_a$ coefficients $a$ are $0, 1/2, -1/2, 0, 0, -1/2, 1/2, 0, 0, 0$. The expectation is the double summation\n",
"\n",
"$$\\displaystyle\\langle \\psi_a|H|\\psi_b\\rangle= \\sum_{k=1}^n\\sum_{j=1}^n a_kb_k \\langle v_k|H|v_j\\rangle$$\n",
"\n",
"where $b$ are the coefficients of the other wavefunction. If the Hamiltonian is of the Huckel type, where interaction is only with adjacent orbitals and zero otherwise, the the bra-ket $\\langle v_k|H|v_j\\rangle=0$ or $1$ depending upon whether the atoms are adjacent or not. \n",
"\n",
"The nine energies (for three degenerate wavefunctions) are placed into a $3 \\times 3$ matrix and the eigenvalues found. The same procedure is followed for the other symmetry orbitals. This makes the calculation far smaller by than calculating a $10 \\times 10$ matrix. The eigenvalues for the $A_u$ species are $(\\pm 5 - 1)/2$.\n",
"\n",
"## Q30 answer\n",
"The reducible representation is A1′+E′ and is worked out using the vectors in the diagram. This accounts for three of the normal modes. The other three are the in-plane and out-of-plane bending modes.\n",
"\n",
"![Drawing](matrices-fig83.png)\n",
"\n",
"Figure 83. Vectors for in-plane stretches in BF$_3$.\n",
"_______\n",
"\n",
"By the method of Section 6.13 the reducible representation is found to be\n",
"\n",
"$$\\displaystyle \\begin{array}{l|rrrrrrr}\n",
"D_{3h} & E & 2C_3 & 3C_2 & \\sigma_h & 2S_3 & 3\\sigma_V \\\\\n",
"\\hline\n",
"\\Gamma_R & 3 & 0 & 1 & 3 & 0 & 1 \\\\\n",
"\\hline \\end{array}$$\n",
"\n",
"Using the tabular method, this works out as\n",
"\n",
"$$\\displaystyle \\begin{array}{l|rrrrrrr r }\n",
"D_{3h} & E & 2C_3 & 3C_2 & \\sigma_h & 2S_3 & 3\\sigma_V &(\\sum)/12\\\\\n",
"\\hline\n",
"\\Gamma_R & 3 & 0 & 1 & 3 & 0 & 1 \\\\\n",
"\\hline\n",
"A_1' & 3 & 0& 3 & 3 & 0 & 3 & 1 \\\\\n",
"A_2' & 3 & 0 & -3 & 3 & 0 & -3 &0 \\\\\n",
"E' & 6 & 0 & 0 & 6 & 0 & 0 & 1 \\\\\n",
"A_1'' & 3 & 0 & 3 & -3 & 0 & -3 & 0 \\\\\n",
"A_2'' & 3 & 0 &-3 & -3 & 0 & 3 & 0 \\\\\n",
"E'' & 6 & 0 & 0 & -6 & 0 & 0 & 0 \\\\\n",
"\\hline \\end{array}$$\n",
"\n",
"Choosing vector $v_1$ as victim, the projection operator method, Section 6.15, produces the following results but only when the full table is used. The classes of operators $2C_3$ etc. have to be expanded so that each is included separately. The (1), (2), and (3) refer to the axes along which the vectors are drawn; the superscripts $\\pm$ to rotation clockwise or anticlockwise by $120^\\text{o}$.\n",
"\n",
"$$\\displaystyle \\begin{array}{l|rrrrrrr r }\n",
"D_{3h} & E & C_3^+ & C_3^- & C_2(1) & C_2(2) & C_2(3) & \\sigma_h & S_3^+ & S_3^- & \\sigma_V(1) & \\sigma_V(2) & \\sigma_V(3) \\\\\n",
"\\hline \n",
"L_{A_1'} & v_1 & v_2 & v_3 & v_1 & v_3 & v_2 &v_1 & v_2 & v_3 & v_1 & v_3 & v_2 \\\\\n",
"L_{E'} & 2v_1 & -v_1 & -v_3 & 0 & 0 & 0 & 2v_1 & -v_2 & -v_3 & 0 & 0 & 0 \\\\\n",
"L_{E'} & 2v_1 & -v_3 & -v_1 & 0 & 0 & 0 & 2v_2 & -v_3 & -v_1 & 0 & 0 & 0 \\\\\n",
"\\hline\n",
"\\end{array}$$\n",
"\n",
"The top row produces the mode pattern for the totally symmetric vibration which is $v_1 + v_2 + v_3$ and shows that in this normal mode each bond stretches together, much as the vectors in the sketch. The normalized vector is $(v_1 + v_2 + v_3)/\\sqrt{3}$. The vectors for the doubly degenerate mode,$2v_1 -v_2 -v_3$ and $-v_1 +2v_2 -v_3$ are calculated by taking vector 1 then 2 but the results are not orthogonal. To check this form, the dot product which is $[2 -1 -1] [-1 2 -1]^T = -3$. Using Sympy these can be made orthogonal, they do not need to be normalized to do this, but the 'True' in the GramSchmidt command normalises the vectors."
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$\\displaystyle \\left[ \\left[\\begin{matrix}\\frac{\\sqrt{6}}{3}\\\\- \\frac{\\sqrt{6}}{6}\\\\- \\frac{\\sqrt{6}}{6}\\end{matrix}\\right], \\ \\left[\\begin{matrix}0\\\\\\frac{\\sqrt{2}}{2}\\\\- \\frac{\\sqrt{2}}{2}\\end{matrix}\\right]\\right]$"
],
"text/plain": [
"⎡⎡ √6 ⎤ ⎤\n",
"⎢⎢ ── ⎥ ⎡ 0 ⎤⎥\n",
"⎢⎢ 3 ⎥ ⎢ ⎥⎥\n",
"⎢⎢ ⎥ ⎢ √2 ⎥⎥\n",
"⎢⎢-√6 ⎥ ⎢ ── ⎥⎥\n",
"⎢⎢────⎥, ⎢ 2 ⎥⎥\n",
"⎢⎢ 6 ⎥ ⎢ ⎥⎥\n",
"⎢⎢ ⎥ ⎢-√2 ⎥⎥\n",
"⎢⎢-√6 ⎥ ⎢────⎥⎥\n",
"⎢⎢────⎥ ⎣ 2 ⎦⎥\n",
"⎣⎣ 6 ⎦ ⎦"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"e1, e2 = symbols('e1, e2')\n",
"\n",
"e1 = Matrix([2,-1,-1])\n",
"e2 = Matrix([-1,2,-1])\n",
"GramSchmidt([e1,e2],True)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The resulting orthogonal and normalised vectors are $(2v_1-v_2-v-3)\\sqrt{6}$ and $(v_2-v_3)\\sqrt{2}$, see figure 84.\n",
"\n",
"The symmetry method used with just three vectors does not allow for any translational motion of the B atom and this is essential to ensure that the centre of gravity does not move. The B atom moves in such a way as to preserve the symmetry. This motion can be imagined as shown by the dashed line in the sketch when the components of the other vectors are scaled for their differences in mass and then added together.\n",
"\n",
"![Drawing](matrices-fig84.png)\n",
"\n",
"Figure 84 Doubly degenerate, vibrational symmetry species. The solid lines show the motion of the F atoms the dashed line that estimated for the B atom to keep the centre of gravity fixed. Neither is to scale. The typical bond displacement of the $v = 0$ level of a vibration is $1$ to $5$% of the bond length, hardly large enough to show on the diagram."
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.9.6"
}
},
"nbformat": 4,
"nbformat_minor": 2
}