Level 4 — ApplicationPhysical Chemistry (Advanced)

Physical Chemistry (Advanced)

60 marksprintable — key stays hidden on paper

Time: 60 minutes Total marks: 60 Instructions: Answer all questions. No hints given. Use standard constants: h=6.626×1034Jsh = 6.626\times10^{-34}\,\mathrm{J\,s}, =1.055×1034Js\hbar = 1.055\times10^{-34}\,\mathrm{J\,s}, kB=1.381×1023JK1k_B = 1.381\times10^{-23}\,\mathrm{J\,K^{-1}}, c=2.998×1010cms1c = 2.998\times10^{10}\,\mathrm{cm\,s^{-1}}, NA=6.022×1023mol1N_A = 6.022\times10^{23}\,\mathrm{mol^{-1}}, F=96485Cmol1F = 96485\,\mathrm{C\,mol^{-1}}, R=8.314Jmol1K1R = 8.314\,\mathrm{J\,mol^{-1}K^{-1}}.


Q1. (Quantum chemistry + variational method) [14 marks]

A conjugated dye electron is modelled as a particle confined in a 1-D box of length LL.

(a) Using the trial function ϕ(x)=x(Lx)\phi(x) = x(L-x) on the interval [0,L][0,L] (which vanishes at both walls), apply the variational principle to estimate the ground-state energy. Show that Etrial=52mL2E_{\text{trial}} = \dfrac{5\hbar^2}{mL^2}. (7)

(b) Compare this to the exact ground-state energy E1=π222mL2E_1 = \dfrac{\pi^2\hbar^2}{2mL^2}. Compute the percentage error. (3)

(c) The dye absorbs at λ=620nm\lambda = 620\,\mathrm{nm} for a transition n=1n=2n=1\to n=2. Treating the electron mass as me=9.109×1031kgm_e = 9.109\times10^{-31}\,\mathrm{kg}, determine the box length LL in nm. (4)


Q2. (Spectroscopy + statistical thermodynamics) [14 marks]

For gaseous 12C16O^{12}\mathrm{C}^{16}\mathrm{O}: rotational constant B=1.931cm1B = 1.931\,\mathrm{cm^{-1}}, vibrational wavenumber ν~=2143cm1\tilde\nu = 2143\,\mathrm{cm^{-1}}.

(a) Compute the bond length (in pm). The reduced mass of CO is μ=1.139×1026kg\mu = 1.139\times10^{-26}\,\mathrm{kg}. (4)

(b) At T=300KT = 300\,\mathrm{K}, estimate the rotational partition function qrot=kBTσhcBq_{\text{rot}} = \dfrac{k_B T}{\sigma\, h c B} (symmetry number σ=1\sigma=1). (4)

(c) Compute the vibrational partition function qvibq_{\text{vib}} (measured from the ground vibrational level) at 300K300\,\mathrm{K}, and state whether vibration contributes significantly to CV,mC_{V,m} at this temperature. Justify. (6)


Q3. (Surface chemistry) [10 marks]

Adsorption of a gas on a catalyst follows the Langmuir isotherm θ=Kp1+Kp\theta = \dfrac{Kp}{1+Kp}.

(a) Two data points are measured: at p1=20kPap_1 = 20\,\mathrm{kPa}, θ1=0.40\theta_1 = 0.40; at p2=60kPap_2 = 60\,\mathrm{kPa}, θ2=0.60\theta_2 = 0.60. Determine whether these data are consistent with a single Langmuir constant KK. Show your working with the two implied KK values. (6)

(b) A student instead fits the low-coverage data to a Freundlich isotherm θ=cp1/n\theta = c\,p^{1/n}. Given θ=0.10\theta = 0.10 at p=5kPap = 5\,\mathrm{kPa} and θ=0.20\theta = 0.20 at p=20kPap = 20\,\mathrm{kPa}, find the exponent 1/n1/n. (4)


Q4. (Advanced electrochemistry) [12 marks]

For an electrode reaction the Tafel equation (anodic branch) is η=a+blog10j\eta = a + b\log_{10} j, with transfer coefficient αa=0.5\alpha_a = 0.5 at T=298KT = 298\,\mathrm{K} and n=1n=1.

(a) Derive the Tafel slope bb (in mV per decade) from the Butler–Volmer equation in its high-overpotential limit. (6)

(b) A current density of j=10mAcm2j = 10\,\mathrm{mA\,cm^{-2}} requires an overpotential η=120mV\eta = 120\,\mathrm{mV}. Find the exchange current density j0j_0. (4)

(c) State one physical consequence of a large j0j_0 for an electrocatalyst. (2)


Q5. (Photochemistry + solid-state) [10 marks]

(a) A photochemical reaction is irradiated at λ=400nm\lambda = 400\,\mathrm{nm} with an incident power of 10mW10\,\mathrm{mW} for 60s60\,\mathrm{s} (all absorbed). 3.5×106mol3.5\times10^{-6}\,\mathrm{mol} of product forms. Compute the quantum yield Φ\Phi. (6)

(b) An intrinsic semiconductor has band gap Eg=1.1eVE_g = 1.1\,\mathrm{eV}. Explain qualitatively how its conductivity changes with rising temperature, contrasting this with a metal. (4)

Answer keyMark scheme & solutions

Q1 [14]

(a) Variational energy E=ϕH^ϕϕϕE = \dfrac{\langle\phi|\hat H|\phi\rangle}{\langle\phi|\phi\rangle}, with H^=22md2dx2\hat H = -\dfrac{\hbar^2}{2m}\dfrac{d^2}{dx^2} inside the box.

  • ϕ=xLx2\phi = xL - x^2, ϕ=2\phi'' = -2. (1)
  • Numerator: 22m0Lϕϕdx=22m(2)0L(xLx2)dx=2mL36=2L36m-\dfrac{\hbar^2}{2m}\displaystyle\int_0^L \phi\phi''\,dx = -\dfrac{\hbar^2}{2m}(-2)\int_0^L (xL-x^2)dx = \dfrac{\hbar^2}{m}\cdot\dfrac{L^3}{6} = \dfrac{\hbar^2 L^3}{6m}. (2) (0L(xLx2)dx=L32L33=L36\int_0^L(xL-x^2)dx = \tfrac{L^3}{2}-\tfrac{L^3}{3}=\tfrac{L^3}{6}).
  • Denominator: 0Lx2(Lx)2dx=L530\displaystyle\int_0^L x^2(L-x)^2 dx = \dfrac{L^5}{30}. (2)
  • E=2L3/6mL5/30=3026mL2=52mL2E = \dfrac{\hbar^2 L^3/6m}{L^5/30} = \dfrac{30\hbar^2}{6mL^2} = \dfrac{5\hbar^2}{mL^2}. (2)

(b) Exact E1=π222mL2=4.93482mL2E_1 = \dfrac{\pi^2\hbar^2}{2mL^2} = \dfrac{4.9348\,\hbar^2}{mL^2}. Error =54.93484.9348×100%=1.32%= \dfrac{5-4.9348}{4.9348}\times100\% = 1.32\% (trial is an upper bound, as required). (3)

(c) ΔE=E2E1=(41)h28mL2=3h28mL2=hcλ\Delta E = E_2 - E_1 = \dfrac{(4-1)h^2}{8mL^2} = \dfrac{3h^2}{8mL^2} = \dfrac{hc}{\lambda}. (2) L2=3hλ8mcL^2 = \dfrac{3h\lambda}{8mc}. L2=3(6.626×1034)(620×109)8(9.109×1031)(2.998×108)L^2 = \dfrac{3(6.626\times10^{-34})(620\times10^{-9})}{8(9.109\times10^{-31})(2.998\times10^{8})} =1.2324×10392.1847×1021=5.641×1019m2= \dfrac{1.2324\times10^{-39}}{2.1847\times10^{-21}} = 5.641\times10^{-19}\,\mathrm{m^2}. L=7.51×1010m=0.751nmL = 7.51\times10^{-10}\,\mathrm{m} = 0.751\,\mathrm{nm}. (2)


Q2 [14]

(a) B=h8π2cIB = \dfrac{h}{8\pi^2 c\, I} (in cm⁻¹ with cc in cm/s), I=μr2I = \mu r^2. I=h8π2cB=6.626×10348π2(2.998×1010)(1.931)I = \dfrac{h}{8\pi^2 c B} = \dfrac{6.626\times10^{-34}}{8\pi^2(2.998\times10^{10})(1.931)} =6.626×10344.574×1012=1.449×1046kgm2= \dfrac{6.626\times10^{-34}}{4.574\times10^{12}} = 1.449\times10^{-46}\,\mathrm{kg\,m^2}. (2) r=I/μ=1.449×1046/1.139×1026=1.272×1020=1.128×1010m=112.8pmr = \sqrt{I/\mu} = \sqrt{1.449\times10^{-46}/1.139\times10^{-26}} = \sqrt{1.272\times10^{-20}} = 1.128\times10^{-10}\,\mathrm{m} = 112.8\,\mathrm{pm}. (2)

(b) qrot=kBThcB=(1.381×1023)(300)(6.626×1034)(2.998×1010)(1.931)q_{\text{rot}} = \dfrac{k_B T}{hcB} = \dfrac{(1.381\times10^{-23})(300)}{(6.626\times10^{-34})(2.998\times10^{10})(1.931)} =4.143×10213.836×1023=108= \dfrac{4.143\times10^{-21}}{3.836\times10^{-23}} = 108. (4)

(c) hcν~kBT=(6.626×1034)(2.998×1010)(2143)(1.381×1023)(300)\dfrac{hc\tilde\nu}{k_BT} = \dfrac{(6.626\times10^{-34})(2.998\times10^{10})(2143)}{(1.381\times10^{-23})(300)} =4.257×10204.143×1021=10.27= \dfrac{4.257\times10^{-20}}{4.143\times10^{-21}} = 10.27. (2) qvib=11e10.27=113.5×1051.000031.0q_{\text{vib}} = \dfrac{1}{1-e^{-10.27}} = \dfrac{1}{1-3.5\times10^{-5}} \approx 1.00003 \approx 1.0. (2) Since qvib1q_{\text{vib}}\approx1, essentially all molecules are in the vibrational ground state; the vibrational mode is "frozen" and contributes negligibly to CV,mC_{V,m} at 300 K (need Tθvib=3084T\gtrsim \theta_{\text{vib}}=3084 K). (2)


Q3 [10]

(a) From θ=Kp1+Kp\theta=\dfrac{Kp}{1+Kp}, solve K=θp(1θ)K = \dfrac{\theta}{p(1-\theta)}. (1)

  • Point 1: K1=0.4020(0.60)=0.4012=0.0333kPa1K_1 = \dfrac{0.40}{20(0.60)} = \dfrac{0.40}{12} = 0.0333\,\mathrm{kPa^{-1}}. (2)
  • Point 2: K2=0.6060(0.40)=0.6024=0.0250kPa1K_2 = \dfrac{0.60}{60(0.40)} = \dfrac{0.60}{24} = 0.0250\,\mathrm{kPa^{-1}}. (2) K1K2K_1 \neq K_2 (0.0333 vs 0.0250, ~25% difference), so the data are not consistent with a single Langmuir constant — the simple Langmuir model does not describe both points. (1)

(b) 1n=ln(θ2/θ1)ln(p2/p1)=ln(0.20/0.10)ln(20/5)=ln2ln4=0.69311.3863=0.500\dfrac{1}{n} = \dfrac{\ln(\theta_2/\theta_1)}{\ln(p_2/p_1)} = \dfrac{\ln(0.20/0.10)}{\ln(20/5)} = \dfrac{\ln 2}{\ln 4} = \dfrac{0.6931}{1.3863} = 0.500. (4)


Q4 [12]

(a) Butler–Volmer: j=j0[eαafηeαcfη]j = j_0\left[e^{\alpha_a f\eta} - e^{-\alpha_c f\eta}\right], f=FRTf=\dfrac{F}{RT}. High anodic overpotential: second term negligible, j=j0eαafηj = j_0 e^{\alpha_a f\eta}. (2) lnj=lnj0+αafηη=RTαaFlnjj0\ln j = \ln j_0 + \alpha_a f\eta \Rightarrow \eta = \dfrac{RT}{\alpha_a F}\ln\dfrac{j}{j_0}. (2) In base-10: η=2.303RTαaFlog10jj0\eta = \dfrac{2.303RT}{\alpha_a F}\log_{10}\dfrac{j}{j_0}, so b=2.303RTαaF=2.303(8.314)(298)0.5(96485)=5706.648242.5=0.1183V=118mV/decadeb = \dfrac{2.303RT}{\alpha_a F} = \dfrac{2.303(8.314)(298)}{0.5(96485)} = \dfrac{5706.6}{48242.5} = 0.1183\,\mathrm{V} = 118\,\mathrm{mV/decade}. (2)

(b) η=blog10(j/j0)log10(j/j0)=120118.3=1.0144\eta = b\log_{10}(j/j_0) \Rightarrow \log_{10}(j/j_0) = \dfrac{120}{118.3} = 1.0144. j/j0=101.0144=10.34j/j_0 = 10^{1.0144} = 10.34. (2) j0=10/10.34=0.967mAcm20.97mAcm2j_0 = 10/10.34 = 0.967\,\mathrm{mA\,cm^{-2}} \approx 0.97\,\mathrm{mA\,cm^{-2}}. (2)

(c) Large j0j_0 ⇒ facile electron transfer, so the same current is delivered at much lower overpotential (lower energy loss / higher efficiency — a better electrocatalyst). (2)


Q5 [10]

(a) Photon energy E=hcλ=(6.626×1034)(2.998×108)400×109=4.966×1019JE = \dfrac{hc}{\lambda} = \dfrac{(6.626\times10^{-34})(2.998\times10^{8})}{400\times10^{-9}} = 4.966\times10^{-19}\,\mathrm{J}. (1) Total energy absorbed =Pt=(0.010)(60)=0.60J= P\,t = (0.010)(60) = 0.60\,\mathrm{J}. (1) Photons =0.60/4.966×1019=1.208×1018= 0.60/4.966\times10^{-19} = 1.208\times10^{18}. Moles of photons =1.208×1018/6.022×1023=2.006×106mol= 1.208\times10^{18}/6.022\times10^{23} = 2.006\times10^{-6}\,\mathrm{mol}. (2) Φ=mol productmol photons=3.5×1062.006×106=1.74\Phi = \dfrac{\text{mol product}}{\text{mol photons}} = \dfrac{3.5\times10^{-6}}{2.006\times10^{-6}} = 1.74. (2) (Φ>1\Phi>1 ⇒ a chain/secondary process amplifies product beyond the Stark–Einstein 1:1 limit.)

(b) In an intrinsic semiconductor, raising TT thermally promotes more electrons across the gap (neEg/2kBTn\propto e^{-E_g/2k_BT}), increasing carrier concentration and hence conductivity — conductivity rises with TT. In a metal, carrier concentration is fixed; increased lattice vibrations enhance scattering, so conductivity falls with TT. Opposite temperature dependences. (4)

[
  {"claim":"Q1a variational energy is 5 hbar^2/(m L^2)","code":"x,L,m,hbar=symbols('x L m hbar',positive=True); phi=x*(L-x); num=integrate(-hbar**2/(2*m)*phi*diff(phi,x,2),(x,0,L)); den=integrate(phi**2,(x,0,L)); E=simplify(num/den); result = simplify(E - 5*hbar**2/(m*L**2))==0"},
  {"claim":"Q1c box length approx 0.751 nm","code":"h=6.626e-34; lam=620e-9; me=9.109e-31; c=2.998e8; L=(3*h*lam/(8*me*c))**0.5; result = abs(L-7.51e-10) < 1e-11"},
  {"claim":"Q2b q_rot approx 108","code":"kB=1.381e-23; T=300; h=6.626e-34; c=2.998e10; B=1.931; q=kB*T/(h*c*B); result = abs(q-108) < 2"},
  {"claim":"Q4a Tafel slope approx 0.1183 V","code":"R=8.314; T=298; F=96485; a=0.5; b=2.303*R*T/(a*F); result = abs(b-0.1183) < 0.001"},
  {"claim":"Q5a quantum yield approx 1.74","code":"h=6.626e-34; c=2.998e8; lam=400e-9; E=h*c/lam; Nphot=0.6/E; mol_phot=Nphot/6.022e23; Phi=3.5e-6/mol_phot; result = abs(Phi-1.74) < 0.05"}
]