Physical Chemistry 1 Quiz: Checking Limiting Cases
20 questions · exam conditions
0:00
Checking Limiting CasesQuestion 1 of 20

The Maxwell-Boltzmann distribution for molecular speeds is f(v)=4π(m2πkT)3/2v2emv2/2kTf(v) = 4\pi\left(\frac{m}{2\pi kT}\right)^{3/2} v^2 e^{-mv^2/2kT}. When checking limiting cases, a student finds that as T0T \to 0, f(v)0f(v) \to 0 for all v>0v > 0 (expected), but as TT \to \infty, f(v)f(v) does not approach a uniform distribution as anticipated. What is the flaw in the student's reasoning about the high-temperature limit?

The Maxwell-Boltzmann distribution becomes invalid at high temperatures due to relativistic effects that require the Maxwell-Jüttner distribution
The normalization condition 0f(v)dv=1\int_0^\infty f(v) dv = 1 prevents the distribution from becoming uniform, regardless of temperature
A uniform distribution over [0,)[0, \infty) cannot be normalized, so the student's expectation is mathematically impossible
At high temperatures, quantum effects become important, requiring Fermi-Dirac or Bose-Einstein statistics instead of classical Maxwell-Boltzmann
The distribution approaches uniform behavior only over a finite range; the exponential tail always dominates at sufficiently high speeds
← Back to quizzes

Physical Chemistry 1 Quiz

Physical Chemistry 1 Quiz: Checking Limiting Cases

Practice Checking Limiting Cases in Physical Chemistry 1 with focused quiz questions that help you check what you know, review explanations, and build confidence with test-style prompts.

What this quiz covers

This quiz focuses on Checking Limiting Cases, giving you a quick way to practice the rules, question types, and explanations that matter most for Physical Chemistry 1.

How to use this quiz

Try each quiz question before looking at the correct answer. Use the explanations to review missed ideas, then come back to similar questions until the pattern feels familiar.

All questions

Question 1

The Maxwell-Boltzmann distribution for molecular speeds is f(v)=4π(m2πkT)3/2v2emv2/2kTf(v) = 4\pi\left(\frac{m}{2\pi kT}\right)^{3/2} v^2 e^{-mv^2/2kT}. When checking limiting cases, a student finds that as T0T \to 0, f(v)0f(v) \to 0 for all v>0v > 0 (expected), but as TT \to \infty, f(v)f(v) does not approach a uniform distribution as anticipated. What is the flaw in the student's reasoning about the high-temperature limit?

  1. The Maxwell-Boltzmann distribution becomes invalid at high temperatures due to relativistic effects that require the Maxwell-Jüttner distribution
  2. The normalization condition 0f(v)dv=1\int_0^\infty f(v) dv = 1 prevents the distribution from becoming uniform, regardless of temperature
  3. A uniform distribution over [0,)[0, \infty) cannot be normalized, so the student's expectation is mathematically impossible (correct answer)
  4. At high temperatures, quantum effects become important, requiring Fermi-Dirac or Bose-Einstein statistics instead of classical Maxwell-Boltzmann
  5. The distribution approaches uniform behavior only over a finite range; the exponential tail always dominates at sufficiently high speeds
Explanation: When analyzing the behavior of probability distributions in limiting cases, you need to consider what types of distributions can actually exist mathematically. The key insight here is understanding what makes a probability distribution valid. The student's expectation that f(v)f(v) should approach a uniform distribution as TT \to \infty reveals a fundamental misunderstanding. A uniform distribution over an infinite interval [0,)[0, \infty) is mathematically impossible because it cannot satisfy the normalization condition 0f(v)dv=1\int_0^\infty f(v) dv = 1. If f(v)=cf(v) = c (constant) for all v0v \geq 0, then 0cdv=\int_0^\infty c \, dv = \infty unless c=0c = 0. Since probability distributions must integrate to 1, a uniform distribution over [0,)[0, \infty) simply cannot exist. This makes C correct - the student's expectation is mathematically impossible from the start. A is wrong because while relativistic effects do matter at extremely high temperatures, this doesn't address the fundamental mathematical flaw in expecting a uniform distribution over an infinite domain. B is incorrect because it suggests the normalization condition somehow "prevents" uniformity, when actually the normalization condition reveals that uniform distributions over infinite intervals are impossible by definition. D is wrong because quantum effects become important at low temperatures (when kTkT approaches the spacing between energy levels), not high temperatures where classical behavior dominates. Study tip: When analyzing limiting behavior of probability distributions, always check whether your expected limiting distribution can actually satisfy the basic requirements of a probability distribution, especially normalization over the given domain.

Question 2

For the van der Waals equation of state, (P+aV2)(Vb)=RT(P + \frac{a}{V^2})(V - b) = RT, a student derives that as VV \to \infty, the equation should reduce to the ideal gas law. However, when checking this limiting case by setting V=106V = 10^6 L/mol in their calculation, they obtain PV/RT=1.15PV/RT = 1.15 instead of 1.00. Which of the following best explains this discrepancy?

  1. The van der Waals equation is fundamentally incorrect at high volumes and cannot reproduce ideal gas behavior
  2. The student used an inappropriate finite value; the limiting behavior requires VV \to \infty mathematically, not just a very large number (correct answer)
  3. The temperature was too low for the ideal gas approximation to be valid, regardless of volume
  4. The van der Waals constants aa and bb are too large for this particular gas, preventing convergence to ideal behavior
  5. The student made an algebraic error; PV/RT=1.15PV/RT = 1.15 is actually consistent with ideal gas behavior within experimental uncertainty
Explanation: When analyzing limiting behavior in thermodynamics, you need to understand the difference between mathematical limits and finite approximations. The van der Waals equation contains correction terms that account for intermolecular forces and molecular volume, and these terms should become negligible as volume approaches infinity. Let's examine what happens mathematically. As VV \to \infty, the term aV20\frac{a}{V^2} \to 0 and bb becomes negligible compared to VV, so (Vb)V(V-b) \to V. This gives us PV=RTP \cdot V = RT, perfectly reproducing the ideal gas law. However, the key word here is "approaches" - it's a limit process, not a substitution. Answer B correctly identifies that using V=106V = 10^6 L/mol, while large, still leaves finite contributions from both correction terms. The aV2\frac{a}{V^2} term equals a1012\frac{a}{10^{12}}, and bb is still subtracted from VV. These seemingly small corrections can easily account for the 15% deviation observed. Answer A is wrong because the van der Waals equation does mathematically reduce to ideal behavior in the proper limit. Answer C incorrectly focuses on temperature when the issue is clearly about the volume approximation. Answer D suggests the constants prevent convergence, but any finite values of aa and bb will eventually become negligible as VV truly approaches infinity. Remember: when checking limiting behavior, don't substitute large finite numbers. Instead, analyze what happens to each term algebraically as the variable approaches its limit.

Question 3

For the Debye model of heat capacity, CV=9R(TΘD)30ΘD/Tx4ex(ex1)2dxC_V = 9R\left(\frac{T}{\Theta_D}\right)^3 \int_0^{\Theta_D/T} \frac{x^4 e^x}{(e^x - 1)^2} dx, a student checks two limiting cases: (1) TΘDT \ll \Theta_D should give CVT3C_V \propto T^3, and (2) TΘDT \gg \Theta_D should give CV=3RC_V = 3R. Their numerical integration confirms case (2) but shows CVT2.8C_V \propto T^{2.8} for case (1). Which error most likely caused this discrepancy?

  1. The integration limits are incorrect; the upper limit should be \infty for the low-temperature case rather than ΘD/T\Theta_D/T
  2. The numerical integration algorithm fails to converge properly when the upper limit ΘD/T\Theta_D/T becomes very large (correct answer)
  3. The student used an inappropriate temperature range; T3T^3 behavior only appears in the quantum limit where T0T \to 0 exactly
  4. The prefactor 9R(T/ΘD)39R(T/\Theta_D)^3 was incorrectly implemented, causing the overall temperature dependence to deviate from cubic
  5. The Debye model breaks down at very low temperatures due to anharmonic effects not captured in the harmonic approximation
Explanation: When analyzing limiting behavior in physical chemistry models, numerical integration challenges often arise when integration limits become extreme. The Debye model presents a classic case where mathematical theory meets computational reality. In the low-temperature limit (TΘDT \ll \Theta_D), the upper integration limit ΘD/T\Theta_D/T becomes very large. Theoretically, as this limit approaches infinity, the integral evaluates to 4π4/154\pi^4/15, giving the expected CVT3C_V \propto T^3 behavior. However, numerical integration algorithms struggle with large upper limits because the integrand x4ex(ex1)2\frac{x^4 e^x}{(e^x - 1)^2} involves exponential terms that can cause overflow errors or convergence issues when xx becomes large. This computational difficulty explains why the student observes T2.8T^{2.8} instead of T3.0T^{3.0} - the numerical integration isn't accurately capturing the tail behavior of the integral, leading to a systematic error that affects the temperature exponent. Option A is incorrect because the finite upper limit ΘD/T\Theta_D/T is mathematically correct for the Debye model. Option C misunderstands the physics - the T3T^3 law appears over a finite temperature range, not just at absolute zero. Option D is wrong because if the prefactor were implemented incorrectly, it would affect the magnitude but not necessarily change the temperature exponent from 3.0 to 2.8 in this systematic way. Study tip: When checking limiting behavior numerically, always verify that your integration algorithm handles extreme limits properly. Consider using adaptive quadrature or variable substitutions to improve convergence for large integration bounds.

Question 4

The Langmuir adsorption isotherm θ=bP1+bP\theta = \frac{bP}{1 + bP} predicts that as P0P \to 0, coverage should be θP\theta \propto P (Henry's law), and as PP \to \infty, θ1\theta \to 1 (saturation). Experimental data shows perfect agreement with the high-pressure limit but deviates from linearity at low pressures, showing θP0.7\theta \propto P^{0.7}. What is the most physically reasonable explanation?

  1. The Langmuir constant bb is pressure-dependent due to lateral interactions between adsorbed molecules that become important at low coverage
  2. Multiple types of adsorption sites exist with different binding energies, violating the Langmuir assumption of equivalent sites (correct answer)
  3. Desorption kinetics become rate-limiting at low pressures, preventing true equilibrium from being established in the experimental timescale
  4. Gas-phase non-ideality becomes significant at low pressures, requiring activity coefficients in place of pressures in the isotherm
  5. Multilayer adsorption occurs even at low pressures, requiring the more complex BET model rather than Langmuir
Explanation: When analyzing deviations from the Langmuir isotherm, you need to consider which fundamental assumptions might be violated. The Langmuir model assumes all adsorption sites are equivalent and independent, leading to the characteristic shape you see here. The key insight is recognizing what θP0.7\theta \propto P^{0.7} represents at low pressures. This fractional exponent is a hallmark of the Freundlich isotherm, which empirically describes adsorption on heterogeneous surfaces. When multiple site types exist with different binding energies, the strongest sites fill first, creating an effective distribution of adsorption energies. This produces the observed power-law behavior at low coverage, while high-pressure saturation remains unchanged since all sites eventually fill. Option A is incorrect because lateral interactions typically become important at high coverage, not low coverage where molecules are far apart. Option C misses the mark because the high-pressure behavior shows perfect equilibrium is achieved—kinetic limitations would affect both pressure regimes. Option D is backwards: gas non-ideality is significant at high pressures where intermolecular forces and finite molecular size matter, not at low pressures where gases behave most ideally. The perfect agreement at high pressures combined with low-pressure deviation strongly points to site heterogeneity (B). Different binding sites create a spectrum of adsorption energies, leading to the gradual filling pattern described by θP0.7\theta \propto P^{0.7}. Study tip: When you see fractional exponents in adsorption isotherms, immediately think "site heterogeneity." The Freundlich isotherm's power-law form is the signature of multiple site types with different binding strengths.

Question 5

For a reversible electrochemical cell, the Nernst equation gives E=ERTnFlnQE = E^\circ - \frac{RT}{nF} \ln Q. A student checks the limiting case where Q0Q \to 0 (reactants in vast excess) and expects E+E \to +\infty, but their calculation gives EE+2.3E \to E^\circ + 2.3 V at 298 K. Which of the following best explains this result?

  1. The student incorrectly used lnQ\ln Q instead of logQ\log Q in the Nernst equation, introducing a factor of 2.303
  2. The limiting case Q0Q \to 0 is physically unrealistic; the student used Q=0.1Q = 0.1 as an approximation, giving a finite correction term (correct answer)
  3. The standard potential EE^\circ includes temperature-dependent terms that become significant only in extreme concentration limits
  4. Activity coefficients become non-ideal at extreme concentrations, requiring the use of activities rather than concentrations in the reaction quotient
  5. The student used n=10n = 10 electrons instead of the correct value, making the logarithmic correction term much smaller than expected
Explanation: When analyzing limiting cases in electrochemistry, you need to consider both mathematical behavior and physical reality. The Nernst equation predicts cell potential based on concentration ratios, but extreme theoretical limits often clash with practical constraints. Let's examine what happens when Q0Q \to 0. Mathematically, ln(0)=\ln(0) = -\infty, so the term RTnFlnQ-\frac{RT}{nF}\ln Q would become ++\infty, making E+E \to +\infty. However, the student's result of EE+2.3E \to E^\circ + 2.3 V suggests they used a finite value for QQ, likely Q=0.1Q = 0.1. This gives RTnFln(0.1)=(8.314)(298)(1)(96485)ln(0.1)0.059-\frac{RT}{nF}\ln(0.1) = -\frac{(8.314)(298)}{(1)(96485)}\ln(0.1) \approx 0.059 V for n=1n=1, but the 2.3 V suggests a different scenario or multiple electrons. Answer B correctly identifies that true Q=0Q = 0 is physically impossible—you cannot have zero concentration of products or infinite excess of reactants in real systems. The student wisely used a reasonable approximation. Answer A is wrong because the Nernst equation correctly uses natural logarithm (ln\ln), not common logarithm. Answer C incorrectly suggests temperature effects in EE^\circ, but standard potentials are defined at standard conditions. Answer D mentions activity coefficients, which is a real consideration but doesn't explain the specific 2.3 V result. Study tip: When working with limiting cases in physical chemistry, always ask yourself: "Is this physically realistic?" Mathematical limits don't always correspond to achievable experimental conditions.

Question 6

For the Michaelis-Menten equation v=Vmax[S]Km+[S]v = \frac{V_{\text{max}}[S]}{K_m + [S]}, the limiting cases are: (1) [S]Km[S] \ll K_m gives v=Vmax[S]Kmv = \frac{V_{\text{max}}[S]}{K_m} (first-order), and (2) [S]Km[S] \gg K_m gives v=Vmaxv = V_{\text{max}} (zero-order). A student's experimental data shows perfect zero-order behavior at high [S][S] but exhibits v[S]1.3v \propto [S]^{1.3} at low [S][S]. Which explanation is most consistent with this deviation?

  1. Substrate inhibition occurs at low concentrations due to non-productive binding modes that reduce effective enzyme concentration
  2. The enzyme exists as multiple oligomerization states with different catalytic activities, creating cooperative binding effects
  3. Product inhibition becomes significant at low substrate concentrations where product accumulation is proportionally more important
  4. Multiple substrate binding sites exist with positive cooperativity, requiring the Hill equation rather than Michaelis-Menten (correct answer)
  5. Diffusion limitations become rate-determining at low substrate concentrations, altering the apparent kinetic order
Explanation: When you encounter deviations from simple Michaelis-Menten kinetics, focus on what the experimental behavior tells you about the underlying mechanism. The key clue here is the fractional power dependence: v[S]1.3v \propto [S]^{1.3} at low substrate concentrations. The correct answer is D because this fractional exponent (1.3) is the hallmark of cooperative binding. When an enzyme has multiple substrate binding sites that exhibit positive cooperativity, the binding of one substrate molecule increases the affinity for subsequent binding events. This creates a sigmoidal relationship between velocity and substrate concentration, described by the Hill equation: v=Vmax[S]nKd+[S]nv = \frac{V_{\max}[S]^n}{K_d + [S]^n}, where n > 1 indicates positive cooperativity. At low [S][S], this reduces to v[S]nv \propto [S]^n, matching your observed [S]1.3[S]^{1.3} dependence. Option A is incorrect because substrate inhibition typically occurs at high concentrations, not low ones, and wouldn't produce fractional power behavior. Option B mentions cooperativity but attributes it to oligomerization states rather than multiple binding sites—the mechanism matters for understanding catalysis. Option C is wrong because product inhibition would decrease the apparent rate but wouldn't create the specific [S]1.3[S]^{1.3} dependence observed. Study tip: Whenever you see fractional or integer exponents greater than 1 in kinetic data, immediately think cooperativity and the Hill equation. The exponent value directly relates to the degree of cooperative interaction between binding sites.

Question 7

For the Einstein model of lattice heat capacity, CV=3R(ΘET)2eΘE/T(eΘE/T1)2C_V = 3R\left(\frac{\Theta_E}{T}\right)^2 \frac{e^{\Theta_E/T}}{(e^{\Theta_E/T} - 1)^2}, the expected limiting behaviors are: (1) TΘET \gg \Theta_E gives CV=3RC_V = 3R, and (2) TΘET \ll \Theta_E gives CVeΘE/TC_V \propto e^{-\Theta_E/T}. A student's numerical model correctly reproduces limit (1) but shows CVT3eΘE/TC_V \propto T^{-3} e^{-\Theta_E/T} for limit (2). What is the most likely explanation for the extra T3T^{-3} factor?

  1. The student forgot to include the prefactor 3R(ΘE/T)23R(\Theta_E/T)^2 and only calculated the exponential portion of the heat capacity expression (correct answer)
  2. Numerical underflow causes eΘE/Te^{\Theta_E/T} in the denominator to be set to zero, leaving only the numerator terms
  3. The Einstein model becomes invalid at very low temperatures where quantum zero-point motion dominates classical harmonic oscillation
  4. The student incorrectly expanded (eΘE/T1)2(e^{\Theta_E/T} - 1)^{-2} using a Taylor series that is only valid for ΘE/T1\Theta_E/T \ll 1
  5. Temperature-dependent phonon-phonon interactions cause deviations from the simple Einstein model at low temperatures
Explanation: When analyzing limiting behaviors of physical models, you need to carefully examine what happens to each term in the expression as the limiting condition is applied. For the low-temperature limit where TΘET \ll \Theta_E, let's work through the correct limiting behavior. When ΘE/T\Theta_E/T becomes very large, eΘE/Te^{\Theta_E/T} dominates over 1 in the denominator, so (eΘE/T1)2(eΘE/T)2(e^{\Theta_E/T} - 1)^2 \approx (e^{\Theta_E/T})^2. This simplifies the expression to: CV3R(ΘET)2eΘE/T(eΘE/T)2=3R(ΘET)2eΘE/TC_V \approx 3R\left(\frac{\Theta_E}{T}\right)^2 \frac{e^{\Theta_E/T}}{(e^{\Theta_E/T})^2} = 3R\left(\frac{\Theta_E}{T}\right)^2 e^{-\Theta_E/T} The student's result shows CVT3eΘE/TC_V \propto T^{-3} e^{-\Theta_E/T} instead of the expected CVeΘE/TC_V \propto e^{-\Theta_E/T}. The extra T3T^{-3} factor comes from (ΘE/T)2=ΘE2T2(\Theta_E/T)^2 = \Theta_E^2 T^{-2}, which along with the missing 3R3R constant suggests the student calculated only the exponential terms and forgot the crucial prefactor 3R(ΘE/T)23R(\Theta_E/T)^2. This makes (A) correct. (B) is wrong because numerical underflow would affect the high-temperature limit too, but that works correctly. (C) is incorrect because the Einstein model remains mathematically valid at low temperatures, even if physically less accurate than the Debye model. (D) is wrong because the Taylor expansion issue would apply to the high-temperature limit (ΘE/T1\Theta_E/T \ll 1), not the low-temperature case. Study tip: When debugging limiting behaviors, check each term in your expression systematically—missing prefactors often reveal themselves through unexpected power-law dependencies.

Question 8

For the Lindemann mechanism of unimolecular reactions, the rate expression is kuni=k1k2[M]k1[M]+k2k_{\text{uni}} = \frac{k_1 k_2 [M]}{k_{-1}[M] + k_2}. In checking limiting cases: (1) high pressure (k1[M]k2k_{-1}[M] \gg k_2) should give kuni=k1k2k1=K1k2k_{\text{uni}} = \frac{k_1 k_2}{k_{-1}} = K_1 k_2, and (2) low pressure (k1[M]k2k_{-1}[M] \ll k_2) should give kuni=k1[M]k_{\text{uni}} = k_1[M]. A student's experimental data shows the expected high-pressure limit but finds kuni[M]1.5k_{\text{uni}} \propto [M]^{1.5} at low pressure instead of the predicted linear dependence. What is the most likely physical explanation?

  1. Falloff behavior occurs due to pressure-dependent energy transfer efficiency, requiring the Troe or other falloff expressions (correct answer)
  2. Multiple conformers of the activated complex exist, each with different pressure dependencies for stabilization
  3. Termolecular collisions become important at low pressures, adding a [M]2[M]^2 term to the activation step
  4. The collision partner [M][M] exhibits non-ideal gas behavior at the pressures studied, requiring activity coefficients
  5. Bath gas heating occurs due to energy release from the reaction, creating a temperature dependence that appears as anomalous pressure dependence
Explanation: When analyzing unimolecular reaction kinetics, deviations from the simple Lindemann mechanism often reveal important physical phenomena about energy transfer processes. The Lindemann model assumes that collisional activation and deactivation are single-step processes with constant efficiency. The correct answer is A because falloff behavior is a well-documented phenomenon where the efficiency of energy transfer during molecular collisions becomes pressure-dependent. At low pressures, not all collisions are equally effective at transferring the required vibrational energy, leading to more complex pressure dependencies than the simple linear relationship predicted by Lindemann theory. The Troe formulation and other falloff expressions specifically account for this by introducing additional parameters that describe how collision efficiency varies with pressure, often resulting in fractional power dependencies like the observed [M]1.5[M]^{1.5}. Option B is incorrect because multiple conformers would typically affect the pre-exponential factor or activation energy, not create a fractional pressure dependence. Option C misunderstands the pressure regime - termolecular collisions become less important at low pressures, not more, and would add complexity to deactivation, not activation. Option D incorrectly invokes non-ideal gas behavior, which typically becomes significant at high pressures where intermolecular forces matter, not at the low pressures described. When you encounter deviations from simple kinetic models, especially fractional power dependencies in pressure, think about the physical assumptions being violated. Energy transfer efficiency is rarely constant across all pressure ranges, making falloff expressions essential for accurate kinetic modeling.

Question 9

For the temperature dependence of equilibrium constants, dlnKdT=ΔHRT2\frac{d \ln K}{dT} = \frac{\Delta H^\circ}{RT^2}, integration gives lnK=ΔHRT+C\ln K = -\frac{\Delta H^\circ}{RT} + C. A student checks the limiting case where ΔH0\Delta H^\circ \to 0 (thermoneutral reaction) and expects KK to become temperature-independent. However, their experimental data for a nearly thermoneutral reaction (ΔH=0.5\Delta H^\circ = -0.5 kJ/mol) still shows significant temperature dependence. Which factor most likely explains this discrepancy?

  1. The student neglected the temperature dependence of ΔH\Delta H^\circ itself, which becomes significant even for small enthalpy changes
  2. Entropy changes (ΔS\Delta S^\circ) contribute to the temperature dependence through ΔG=ΔHTΔS\Delta G^\circ = \Delta H^\circ - T\Delta S^\circ, independent of enthalpy (correct answer)
  3. The reaction is close to equilibrium, making small experimental errors in KK appear as large relative changes with temperature
  4. Non-ideal solution behavior causes activity coefficients to vary with temperature, affecting the apparent equilibrium constant
  5. The van 't Hoff equation is only valid for gas-phase reactions and breaks down for condensed-phase systems near thermoneutrality
Explanation: When analyzing equilibrium constant temperature dependence, you need to consider the complete thermodynamic relationship, not just the enthalpy term. The van't Hoff equation dlnKdT=ΔHRT2\frac{d \ln K}{dT} = \frac{\Delta H^\circ}{RT^2} comes from the fundamental relationship ΔG=RTlnK=ΔHTΔS\Delta G^\circ = -RT \ln K = \Delta H^\circ - T\Delta S^\circ. The correct answer is B because even when ΔH0\Delta H^\circ \approx 0, the entropy term TΔST\Delta S^\circ still creates temperature dependence. Taking the derivative of ΔG=ΔHTΔS\Delta G^\circ = \Delta H^\circ - T\Delta S^\circ with respect to temperature gives contributions from both enthalpy and entropy changes. For a truly thermoneutral reaction, KK would only be temperature-independent if both ΔH=0\Delta H^\circ = 0 AND ΔS=0\Delta S^\circ = 0, which is extremely rare. Option A is incorrect because temperature dependence of ΔH\Delta H^\circ (heat capacity effects) is typically small over moderate temperature ranges and wouldn't explain "significant" temperature dependence. Option C misses the point—the issue isn't experimental error but fundamental thermodynamics. The relative changes observed are real, not artifacts. Option D introduces unnecessary complexity; while activity coefficients do vary with temperature, this doesn't explain why a nearly thermoneutral reaction shows temperature dependence when the student expected none. Remember: equilibrium constants depend on ΔG\Delta G^\circ, which has both enthalpy AND entropy components. Even if ΔH0\Delta H^\circ \approx 0, significant ΔS\Delta S^\circ values will cause KK to vary with temperature through the TΔST\Delta S^\circ term.

Question 10

The Stern-Volmer equation for fluorescence quenching is F0F=1+KSV[Q]\frac{F_0}{F} = 1 + K_{SV}[Q], where F0F_0 and FF are fluorescence intensities without and with quencher at concentration [Q][Q]. The limiting cases should be: (1) [Q]=0[Q] = 0 gives F=F0F = F_0, and (2) [Q][Q] \to \infty gives F0F \to 0. Experimental data confirms case (1) but shows that FF approaches a non-zero constant value at high [Q][Q]. Which explanation is most consistent with this observation?

  1. Static quenching occurs in addition to dynamic quenching, requiring a modified Stern-Volmer equation with exponential terms
  2. A fraction of the fluorophore population is protected from quenching by association with other molecules or interfaces (correct answer)
  3. Quencher aggregation at high concentrations reduces the effective quenching efficiency per molecule
  4. Inner filter effects cause significant reabsorption of emitted light, reducing the apparent quenching efficiency
  5. The excited state lifetime becomes comparable to the diffusion timescale, invalidating the steady-state approximation
Explanation: When analyzing deviations from ideal Stern-Volmer behavior, you need to consider what physical factors might protect some fluorophores from quenching. The classic Stern-Volmer equation assumes all fluorophores are equally accessible to quenchers, but real systems often have heterogeneous populations. The key insight is that if fluorescence approaches a non-zero plateau at high quencher concentrations, some fraction of fluorophores must be completely inaccessible to quenching. This creates a two-population system: one that can be quenched normally and another that remains protected regardless of quencher concentration. Answer B correctly identifies this scenario. When fluorophores associate with proteins, membranes, or other protective environments, they become shielded from quenchers. The modified equation becomes F0F=fa1+KSV[Q]+fp\frac{F_0}{F} = \frac{f_a}{1 + K_{SV}[Q]} + f_p, where faf_a is the accessible fraction and fpf_p is the protected fraction. As [Q][Q] \to \infty, fluorescence approaches fpF0f_p \cdot F_0, not zero. Answer A incorrectly invokes static quenching, which would actually enhance quenching efficiency, not reduce it. Answer C suggests quencher aggregation reduces efficiency, but this would typically show a downward curvature in Stern-Volmer plots, not a plateau. Answer D mentions inner filter effects, which cause apparent intensity reductions due to optical artifacts but don't create the specific plateau behavior described. Remember: when Stern-Volmer plots show residual fluorescence at high quencher concentrations, think about fluorophore accessibility and microenvironmental protection, not just quenching mechanisms.

Question 11

For the partition function of a two-level system, q=1+eΔE/kTq = 1 + e^{-\Delta E/kT}, the limiting cases are: (1) T0T \to 0 gives q1q \to 1 (ground state only), and (2) TT \to \infty gives q2q \to 2 (equal populations). A student calculates the heat capacity C=k(ΔEkT)2eΔE/kT(1+eΔE/kT)2C = k\left(\frac{\Delta E}{kT}\right)^2 \frac{e^{-\Delta E/kT}}{(1 + e^{-\Delta E/kT})^2} and finds that as TT \to \infty, C0C \to 0 instead of remaining finite. Is this result physically reasonable?

  1. No, the heat capacity of any finite system should approach Nk/2Nk/2 per degree of freedom at high temperature according to equipartition
  2. Yes, once thermal equilibrium between levels is achieved, no additional heat is required to maintain the temperature, so C0C \to 0
  3. No, the student made an algebraic error; the correct high-temperature limit should give CkC \to k for a two-level system
  4. Yes, the heat capacity measures the temperature dependence of energy, and for a two-level system, energy approaches a constant at high TT (correct answer)
  5. No, quantum effects become negligible at high temperature, requiring classical statistical mechanics where heat capacity diverges
Explanation: When analyzing heat capacity in statistical mechanics, focus on what happens to the energy distribution as temperature changes. Heat capacity measures how much energy is required to increase temperature, which depends on how the population distribution between energy levels responds to temperature changes. For this two-level system, at high temperatures both levels become equally populated (q2q \to 2). Once this equal distribution is reached, the average energy approaches a constant value of ΔE/2\Delta E/2. Since heat capacity measures dEdT\frac{dE}{dT}, and energy stops changing significantly with further temperature increases, C0C \to 0 is indeed physically correct. Let's examine why the other answers miss the mark: A incorrectly applies the equipartition theorem, which applies to continuous degrees of freedom (like translational motion), not discrete two-level systems. The Nk/2Nk/2 per degree of freedom rule doesn't apply here. B misunderstands what thermal equilibrium means. Thermal equilibrium refers to temperature being uniform throughout the system, not to achieving equal populations between energy levels. The reasoning about "no additional heat required" is also flawed. C suggests an algebraic error, but you can verify the math is correct. At high TT, eΔE/kT1e^{-\Delta E/kT} \to 1, making the numerator approach 01=00 \cdot 1 = 0 while the denominator approaches (1+1)2=4(1+1)^2 = 4, giving C0C \to 0. Key insight: For discrete energy systems, heat capacity can approach zero at high temperatures because energy becomes saturated. This differs fundamentally from continuous systems where equipartition applies. Always distinguish between discrete and continuous degrees of freedom in statistical mechanics problems.

Question 12

The Marcus equation for electron transfer rates is kET=Ae(ΔG+λ)2/(4λRT)k_{ET} = A e^{-(\Delta G^\circ + \lambda)^2/(4\lambda RT)}, where λ\lambda is the reorganization energy. When checking limiting cases, a student finds that as ΔG0\Delta G^\circ \to 0 (thermoneutral), kET=Aeλ/(4RT)k_{ET} = A e^{-\lambda/(4RT)}, and as ΔGλ\Delta G^\circ \to -\lambda (optimal driving force), kET=Ak_{ET} = A. However, for ΔGλ\Delta G^\circ \ll -\lambda (very exergonic), they calculate kET0k_{ET} \to 0 and question whether this "inverted region" is physically meaningful.

  1. The inverted region is an artifact of the Marcus model; real electron transfer rates continue to increase with more negative ΔG\Delta G^\circ
  2. The inverted region reflects the quantum mechanical nature of electron tunneling, which decreases when energy gaps become too large
  3. The Marcus equation becomes invalid in the highly exergonic regime due to breakdown of the harmonic approximation for nuclear motion
  4. The inverted region occurs because highly exergonic reactions require more nuclear reorganization to match the product state energy (correct answer)
  5. The apparent rate decrease reflects competing reaction pathways that become dominant when the primary electron transfer is too fast to measure
Explanation: When analyzing electron transfer kinetics, the Marcus equation reveals a fascinating phenomenon where reaction rates don't always increase with driving force. The key insight lies in understanding how nuclear reorganization couples to electron transfer. The correct answer is D because highly exergonic reactions create an energy mismatch problem. When ΔGλ\Delta G^\circ \ll -\lambda, the reaction releases so much energy that the nuclear coordinates must reorganize extensively to reach a configuration where reactant and product electronic states have the same energy (the transition state). This large reorganization requirement increases the activation barrier, paradoxically slowing the reaction despite its high thermodynamic favorability. Option A is wrong because the inverted region has been experimentally observed in many systems, particularly in photosynthetic reaction centers and artificial donor-acceptor complexes. It's a real physical phenomenon, not a theoretical artifact. Option B incorrectly attributes the effect to quantum tunneling. While electron transfer does involve quantum mechanical tunneling, the inverted region specifically arises from classical nuclear reorganization dynamics, not tunneling distance effects. Option C misidentifies the limitation. The harmonic approximation for nuclear motion generally remains valid even in highly exergonic regimes. The inverted region emerges naturally from the harmonic model itself. Remember this key principle: in Marcus theory, the optimal reaction rate occurs when the driving force exactly matches the reorganization energy (ΔG=λ\Delta G^\circ = -\lambda). Deviations in either direction—too little or too much driving force—increase the activation barrier and slow the reaction.

Question 13

For the Fermi-Dirac distribution f(E)=1e(Eμ)/kT+1f(E) = \frac{1}{e^{(E-\mu)/kT} + 1}, the limiting cases are: (1) T0T \to 0 gives a step function (f=1f = 1 for E<μE < \mu, f=0f = 0 for E>μE > \mu), and (2) TT \to \infty should approach the Maxwell-Boltzmann distribution. A student checks case (2) by expanding for kTEμkT \gg |E - \mu| and obtains f(E)12Eμ4kTf(E) \approx \frac{1}{2} - \frac{E - \mu}{4kT}. They conclude this linear form contradicts the expected exponential Maxwell-Boltzmann behavior. What is the flaw in their analysis?

  1. The expansion is only first-order; higher-order terms in (Eμ)/(kT)(E-\mu)/(kT) are needed to recover the exponential form
  2. The condition kTEμkT \gg |E - \mu| is too restrictive; the classical limit requires EμkTE - \mu \gg kT instead (correct answer)
  3. The chemical potential μ\mu becomes temperature-dependent at high TT, invalidating the expansion around a fixed energy
  4. The Fermi-Dirac distribution only approaches Maxwell-Boltzmann for EμE \gg \mu, not in the regime where EμE \approx \mu
  5. The student used an incorrect Taylor expansion; the correct first-order approximation should be f(E)e(Eμ)/kTf(E) \approx e^{-(E-\mu)/kT}
Explanation: This question tests your understanding of when and how the Fermi-Dirac distribution transitions to classical Maxwell-Boltzmann statistics. The key insight is recognizing the proper limiting conditions for this transition. The student's error lies in choosing the wrong regime for the classical limit. When kTEμkT \gg |E - \mu|, you're examining energies very close to the chemical potential, where quantum effects (Pauli exclusion) remain important. The correct classical limit requires EμkTE - \mu \gg kT, meaning you need energies well above the chemical potential where orbital occupancy becomes sparse. In the proper classical regime (EμkTE - \mu \gg kT), the exponential term e(Eμ)/kTe^{(E-\mu)/kT} becomes much larger than 1, so the "+1" in the denominator becomes negligible. This gives f(E)1e(Eμ)/kT=e(Eμ)/kTf(E) \approx \frac{1}{e^{(E-\mu)/kT}} = e^{-(E-\mu)/kT}, which is indeed the Maxwell-Boltzmann exponential form. Option A is wrong because higher-order terms won't convert the linear expansion into an exponential—the fundamental issue is the regime choice. Option C incorrectly suggests the chemical potential's temperature dependence invalidates the analysis, but this doesn't address the core issue of which energy regime to examine. Option D is partially correct about needing EμE \gg \mu, but it doesn't clearly identify that the student used the wrong limiting condition. Study tip: For distribution function limits, always identify the physical regime first. Classical behavior emerges when quantum restrictions become negligible, which happens at high energies relative to the Fermi level, not near it.

Question 14

A researcher models the temperature dependence of a reaction rate using the Arrhenius equation k=AeEa/RTk = A e^{-E_a/RT}. When checking the limiting case as T0T \to 0, they expect k0k \to 0, but their computational model gives undefined results. Simultaneously, as TT \to \infty, they expect kAk \to A, but their model predicts exponential divergence. What is the most likely explanation for both discrepancies?

  1. The Arrhenius equation is only valid within a narrow temperature range and breaks down at extreme temperatures due to quantum effects
  2. The activation energy EaE_a has been incorrectly assigned a negative value, causing the exponential behavior to invert (correct answer)
  3. The pre-exponential factor AA is temperature-dependent, violating the assumptions of the simple Arrhenius model at extreme temperatures
  4. Computational overflow occurs because the exponential function exceeds machine precision at very high and very low temperature limits
  5. The gas constant RR requires temperature-dependent corrections that become significant only at extreme temperatures
Explanation: When analyzing the Arrhenius equation's behavior at temperature extremes, you need to carefully examine how the exponential term eEa/RTe^{-E_a/RT} behaves based on the sign of the activation energy. If EaE_a is incorrectly assigned a negative value, the equation becomes k=Ae+Ea/RTk = A e^{+|E_a|/RT}. This completely inverts the expected temperature dependence. As T0T \to 0, the positive exponent +Ea/RT+|E_a|/RT approaches ++\infty, making kk approach infinity rather than zero—causing computational overflow and undefined results. As TT \to \infty, the exponent approaches zero from the positive side, so kk approaches AA, but the model shows exponential growth during the approach rather than the expected exponential decay. Option A is incorrect because while quantum effects do modify reaction kinetics at extreme temperatures, they wouldn't cause the specific computational issues described—undefined results at low T and exponential divergence at high T. Option C is wrong because temperature-dependent pre-exponential factors would cause gradual deviations from expected behavior, not the dramatic computational failures observed. Option D misidentifies the problem as purely computational. While overflow can occur with large exponentials, the issue here is fundamental—the wrong sign creates physically unrealistic behavior that leads to computational problems. Study tip: Always check that your activation energy is positive in Arrhenius calculations. A negative EaE_a is a red flag that immediately tells you something is wrong with your model or data entry, as it predicts faster reactions at lower temperatures.

Question 15

For the Arrhenius temperature dependence of diffusion, D=D0eEa/RTD = D_0 e^{-E_a/RT}, a student analyzes self-diffusion data for a crystalline solid and checks two limiting cases: (1) T0T \to 0 should give D0D \to 0 (no thermal motion), and (2) TTmeltT \to T_{\text{melt}} should show deviation from Arrhenius behavior. Their data confirms case (1) but shows perfect Arrhenius behavior right up to the melting point. Which explanation is most physically reasonable?

  1. The activation energy decreases with temperature due to thermal expansion, maintaining Arrhenius form but with temperature-dependent parameters
  2. Multiple diffusion mechanisms exist with different activation energies, but only one mechanism dominates over the entire temperature range studied
  3. The crystal structure remains stable without defect formation or phase transitions up to the melting point, maintaining a single diffusion pathway
  4. Cooperative effects between diffusing atoms become important near melting, but these effects exactly compensate for other deviations
  5. The temperature range studied is insufficient to observe pre-melting effects, which only occur within a few degrees of TmeltT_{\text{melt}} (correct answer)
Explanation: When analyzing diffusion data that follows Arrhenius behavior across an entire temperature range, you need to consider what physical processes could maintain this mathematical form even as conditions change dramatically. The most physically reasonable explanation is that multiple diffusion mechanisms with different activation energies are present, but only one dominates throughout the studied range (Answer B). In crystalline solids, diffusion can occur through vacancy mechanisms, interstitial mechanisms, or grain boundary diffusion, each with distinct activation energies. However, one mechanism often has significantly lower energy barriers, making it dominant across the entire temperature range. This single dominant pathway maintains the simple Arrhenius form D=D0eEa/RTD = D_0 e^{-E_a/RT} with constant parameters. Answer A is incorrect because while thermal expansion does affect activation energy, this would cause observable curvature in an Arrhenius plot (ln D vs 1/T), not perfect linearity. Answer C is physically unrealistic—all crystals develop increasing defect concentrations as temperature rises, and these defects are essential for diffusion to occur. Answer D represents an extremely unlikely coincidence where multiple competing effects would need to cancel exactly across a wide temperature range. Study tip: When you see perfect Arrhenius behavior over a large temperature range, think about mechanism dominance rather than the absence of other effects. Real materials always have multiple competing processes, but often one mechanism is so favorable that it masks the others, giving the appearance of simple behavior.

Question 16

The Planck distribution for blackbody radiation is u(ν)=8πhν3c31ehν/kT1u(\nu) = \frac{8\pi h\nu^3}{c^3} \frac{1}{e^{h\nu/kT} - 1}. In checking limiting cases, a student verifies that as TT \to \infty (high temperature), the distribution approaches the Rayleigh-Jeans law u(ν)=8πν2kTc3u(\nu) = \frac{8\pi\nu^2 kT}{c^3}. However, when checking the low-temperature limit (T0T \to 0), they expect u(ν)0u(\nu) \to 0 but find that their numerical calculation gives a finite, non-zero result. What is the most likely source of this error?

  1. Machine precision limitations cause ehν/kTe^{h\nu/kT} to be evaluated as exactly zero when the exponent becomes very large
  2. The student used a finite temperature value instead of taking the true mathematical limit T0T \to 0
  3. Floating-point overflow occurs when ehν/kTe^{h\nu/kT} exceeds the maximum representable number, causing undefined behavior (correct answer)
  4. The low-temperature limit requires quantum corrections to the Planck distribution that are not captured in the classical formula
  5. The student incorrectly applied Wien's displacement law, which shifts the peak frequency but doesn't affect the limiting behavior
Explanation: When analyzing limiting behavior of physical equations computationally, you need to consider how computers handle extreme numerical values, not just the mathematical theory. The Planck distribution contains the term ehν/kTe^{h\nu/kT} in the denominator. As T0T \to 0, this exponent becomes extremely large, making ehν/kTe^{h\nu/kT} astronomically big. When this exponential exceeds the maximum number your computer can represent (typically around 1030810^{308} for double precision), floating-point overflow occurs. Instead of gracefully approaching infinity, the calculation breaks down and returns special values like "infinity" or "NaN" (Not a Number). When these undefined values propagate through the rest of the calculation, you get nonsensical finite results instead of the expected zero. Answer A is backwards—the issue isn't the exponential being too small (underflow), but too large (overflow). Answer B misses the point; even using small finite temperatures should show the trend toward zero, so this isn't about mathematical vs. numerical limits. Answer D is incorrect because the Planck distribution already incorporates full quantum mechanics—no additional corrections are needed for low temperatures. The key insight is that ehν/kT1ehν/kTe^{h\nu/kT} - 1 \approx e^{h\nu/kT} when the exponent is large, so u(ν)8πhν3c3ehν/kTu(\nu) \approx \frac{8\pi h\nu^3}{c^3} e^{-h\nu/kT}, which should decay exponentially to zero as T0T \to 0. Study tip: When computational results contradict expected limiting behavior, always check for numerical overflow or underflow issues before questioning the underlying physics.

Question 17

The Schrödinger equation for a particle in a 1D infinite potential well gives energy levels En=n2π222mL2E_n = \frac{n^2\pi^2\hbar^2}{2mL^2}. When checking the classical limit, a student expects that as 0\hbar \to 0, the energy levels should become continuous (En0E_n \to 0). However, they realize this would violate energy conservation since the particle must have kinetic energy. Which resolution of this apparent paradox is most physically sound?

  1. The classical limit requires LL \to \infty simultaneously with 0\hbar \to 0 such that 2/L2\hbar^2/L^2 remains finite
  2. In the classical limit, the particle escapes the infinite potential well, so the quantum energy formula no longer applies
  3. The classical limit is achieved when nn \to \infty while keeping EnE_n finite, requiring \hbar to scale appropriately with nn
  4. The energy levels approach zero spacing (ΔE0\Delta E \to 0) rather than individual levels approaching zero, creating effectively continuous energy (correct answer)
  5. The infinite potential well is inherently quantum mechanical and has no valid classical analogue, making the classical limit undefined
Explanation: When you encounter questions about classical limits in quantum mechanics, the key insight is understanding how quantum discreteness transitions to classical continuity. The student's confusion here stems from conflating individual energy levels with the energy spectrum's density. The correct resolution (D) recognizes that the classical limit emerges when energy level spacing becomes infinitesimally small, not when individual energies approach zero. As quantum numbers increase, adjacent energy levels get closer: ΔE=En+1En=(2n+1)π222mL2\Delta E = E_{n+1} - E_n = \frac{(2n+1)\pi^2\hbar^2}{2mL^2}. For large nn, this spacing approaches zero relative to the total energy, creating a quasi-continuous spectrum that matches classical expectations where a particle can have any kinetic energy value. Option A incorrectly suggests manipulating both \hbar and LL simultaneously, which doesn't address the fundamental issue of energy quantization versus continuity. Option B misunderstands the infinite potential well model—the particle remains confined in both quantum and classical treatments; the walls don't disappear in the classical limit. Option C puts the focus on scaling \hbar with nn, but this artificially manipulates the fundamental constant rather than recognizing the natural emergence of continuity through dense energy spacing. The key study insight: classical limits in quantum mechanics typically involve making discrete quantities so densely packed that they appear continuous, not making individual quantum values approach zero. Remember that \hbar is a fundamental constant—the classical limit emerges through system parameters (like high quantum numbers), not by artificially changing physical constants.

Question 18

The Debye-Hückel equation for activity coefficients is lnγ±=Az+zI1+BaI\ln \gamma_\pm = -\frac{A z_+ z_- \sqrt{I}}{1 + B a \sqrt{I}}. When checking limiting cases, a student finds that as I0I \to 0, γ±1\gamma_\pm \to 1 (ideal solution, expected), but as II \to \infty, they calculate γ±eAz+z/Ba\gamma_\pm \to e^{-A z_+ z_-/Ba} instead of the expected approach to a constant value. What error did the student likely make?

  1. They incorrectly took the limit of I\sqrt{I} in the numerator while keeping the denominator finite, violating L'Hôpital's rule requirements
  2. They assumed the denominator approaches BaIBa\sqrt{I} as II \to \infty, leading to lnγ±Az+z/(Ba)\ln \gamma_\pm \to -A z_+ z_-/(Ba) (correct answer)
  3. They failed to recognize that the Debye-Hückel equation is only valid for low ionic strengths and breaks down as II \to \infty
  4. They confused the extended Debye-Hückel equation with the limiting law, using the wrong expression for high ionic strength
  5. They neglected higher-order terms in the virial expansion that become important at high ionic strength, invalidating the limiting analysis
Explanation: When analyzing limiting behavior of equations, you need to carefully track how each term behaves as the variable approaches its limit. The Debye-Hückel equation contains both numerator and denominator terms that depend on I\sqrt{I}. The correct analysis for II \to \infty requires examining the full expression. As ionic strength becomes very large, I\sqrt{I} grows without bound. In the numerator, you have Az+zIA z_+ z_- \sqrt{I}, which increases linearly with I\sqrt{I}. In the denominator, you have 1+BaI1 + B a \sqrt{I}. For large II, the BaIBa\sqrt{I} term dominates over the constant 1, so the denominator effectively becomes BaIBa\sqrt{I}. The student's error was assuming this approximation too hastily. They treated the denominator as simply BaIBa\sqrt{I}, making the expression lnγ±Az+zIBaI=Az+zBa\ln \gamma_\pm \approx -\frac{A z_+ z_- \sqrt{I}}{Ba\sqrt{I}} = -\frac{A z_+ z_-}{Ba}, which gives γ±eAz+z/Ba\gamma_\pm \to e^{-A z_+ z_-/Ba} as stated in choice B. Choice A incorrectly suggests L'Hôpital's rule violations, but this isn't an indeterminate form requiring that rule. Choice C mentions the equation's validity range, which is true but doesn't explain the specific mathematical error. Choice D confuses different forms of the equation, but the student used the correct extended form. Study tip: When taking limits of rational functions, always consider whether constant terms become negligible compared to growing terms, but don't drop them prematurely in your algebra—complete the limit process rigorously first.

Question 19

The Poisson-Boltzmann equation for electrostatic potential around a charged sphere is 2ψ=κ2sinh(eψkT)\nabla^2 \psi = \kappa^2 \sinh\left(\frac{e\psi}{kT}\right), where κ\kappa is the Debye-Hückel parameter. For weak potentials, this linearizes to 2ψ=κ2ψ\nabla^2 \psi = \kappa^2 \psi. A student checks the limiting case κ0\kappa \to 0 (very low ionic strength) and expects the solution to approach the vacuum Coulomb potential ψ1/r\psi \propto 1/r. However, their numerical solution shows exponential decay even at κ=106\kappa = 10^{-6} m1^{-1}. What is the most likely explanation?

  1. The boundary conditions were incorrectly implemented, preventing the solution from matching the Coulomb form at large distances
  2. Numerical precision limitations cause small non-zero values of κ2\kappa^2 to be treated as finite, maintaining exponential screening
  3. The linearization is invalid even for weak potentials when κ0\kappa \to 0, requiring the full nonlinear Poisson-Boltzmann equation
  4. The computational domain is finite, and exponential boundary conditions at the edges artificially impose screening (correct answer)
  5. Machine precision errors accumulate during the iterative solution, creating spurious exponential behavior that dominates at large distances
Explanation: When solving differential equations numerically, boundary conditions and computational domain size critically affect the solution behavior, especially when examining limiting cases. The linearized Poisson-Boltzmann equation 2ψ=κ2ψ\nabla^2 \psi = \kappa^2 \psi has the general solution ψ=Areκr+Breκr\psi = \frac{A}{r}e^{-\kappa r} + \frac{B}{r}e^{\kappa r} in spherical coordinates. As κ0\kappa \to 0, this should approach the Coulomb potential ψ=Ar\psi = \frac{A}{r}. However, numerical simulations require finite computational domains with boundary conditions at some outer radius rmaxr_{max}. The correct answer is D because when you impose boundary conditions that force ψ0\psi \to 0 at rmaxr_{max} (typical for computational convenience), you're artificially requiring exponential decay. Even tiny values of κ\kappa will dominate the solution form over large distances, preventing the 1/r1/r Coulomb behavior from emerging within your finite domain. A is incorrect because standard boundary conditions (ψ=0\psi = 0 at infinity) are physically reasonable—the issue is implementing "infinity" as a finite boundary. B is wrong because numerical precision of 10610^{-6} is easily handled by modern computers without significant roundoff errors affecting the qualitative solution form. C is incorrect because the linearization becomes increasingly valid as κ0\kappa \to 0 since the potential weakens and eψ/kT1+ψ/kTe\psi/kT \approx 1 + \psi/kT. Study tip: When numerically solving equations with known analytical limits, always check that your boundary conditions and domain size don't artificially constrain the solution away from the expected limiting behavior.

Question 20

A researcher derives the following expression for the heat capacity of a diatomic gas: CV=72nR+3nRθ2eθ/T2T2(eθ/T1)2C_V = \frac{7}{2}nR + \frac{3nR\theta^2 e^{\theta/T}}{2T^2(e^{\theta/T} - 1)^2}, where θ\theta is a vibrational characteristic temperature. When checking the high-temperature limit (TθT \gg \theta), which result would indicate an error in the derivation?

  1. CV52nRC_V \to \frac{5}{2}nR, indicating only translational and rotational contributions survive while vibrational terms vanish unexpectedly
  2. CV92nRC_V \to \frac{9}{2}nR, which correctly accounts for full activation of translational, rotational, and vibrational degrees of freedom
  3. CV72nRC_V \to \frac{7}{2}nR, suggesting vibrational modes become inactive at high temperatures contrary to classical equipartition predictions (correct answer)
  4. CV72nR+3nR2C_V \to \frac{7}{2}nR + \frac{3nR}{2}, which properly includes the classical limit contribution from vibrational modes
Explanation: When analyzing heat capacity expressions for diatomic gases, you need to consider how molecular motion behaves at different temperatures according to classical equipartition theory and quantum mechanics. To find the high-temperature limit, examine what happens to the vibrational term when TθT \gg \theta. As temperature increases, θ/T\theta/T becomes very small, so eθ/T1+θ/Te^{\theta/T} \approx 1 + \theta/T. Substituting this approximation into the vibrational term: 3nRθ2eθ/T2T2(eθ/T1)23nRθ2(1+θ/T)2T2(θ/T)2=3nR2\frac{3nR\theta^2 e^{\theta/T}}{2T^2(e^{\theta/T} - 1)^2} \approx \frac{3nR\theta^2(1 + \theta/T)}{2T^2(\theta/T)^2} = \frac{3nR}{2} Therefore, the correct high-temperature limit should be CV72nR+3nR2=102nR=5nRC_V \to \frac{7}{2}nR + \frac{3nR}{2} = \frac{10}{2}nR = 5nR. Option A incorrectly suggests vibrational contributions vanish at high temperatures, but classical equipartition predicts they should be fully active. Option B gives 92nR\frac{9}{2}nR, which doesn't match the correct mathematical limit. Option D correctly shows the proper classical limit structure with both the original 72nR\frac{7}{2}nR term and the 3nR2\frac{3nR}{2} vibrational contribution. Option C indicates an error because it suggests the vibrational term contributes nothing at high temperatures (CV72nRC_V \to \frac{7}{2}nR), contradicting classical equipartition theory, which predicts vibrational modes should contribute nRnR per mode at high temperatures. Remember: when checking quantum mechanical expressions in classical limits, vibrational contributions should approach their equipartition values, not disappear.