Skip to main content
Course and syllabus guidance

University Advanced Quantum Mechanics · unit UAQM-08

Negative Energies, Hole Theory and Charge Conjugation

Connect the ideas, make your working visible, and test the result. Use the hints when you need a starting point; open the full solution after your attempt.

Extra challenging questions

A selection from this pathway. Open a topic below for its complete practice and challenge set.

Extra challenging · 01

Negative-energy components and Zitterbewegung: infer an input with uncertainty

Zitterbewegung is the oscillation of a free Dirac electron's position at angular frequency 2E/ħ, with amplitude of order ħ/2mc, produced by interference between the positive- and negative-energy components of its spinor. Use the supplied topic model Y_out = (1-2*w)*q/sqrt(1+q^2), where Y_out is DRIFT VELOCITY (c). The output symbol Y_out is distinct from all input symbols. ^ denotes a power; π = 3.141592653589793; sqrt(x) means √x. q is Momentum p, measured in mc. Keep w (Negative-energy weight w₋) = 0.2 fixed. For this exercise use only 0.31875 mc ≤ q ≤ 2.68125 mc on this continuous branch. A grid of test settings is an input restriction for this exercise, not by itself a claim of physical quantisation. This calculation retains the topic model's stated assumptions; it does not extend them to a different physical setup. A reading gives Y_out = 0.52166489 c, with a bounded reading uncertainty ±0.0060763 c. Treat the model and all fixed inputs as exact for this exercise; the bound is not a standard deviation. First estimate the unknown q by straight-line interpolation between the two endpoint predictions. Then solve the actual model on the stated branch to an absolute input tolerance of 0.00023625 mc, using a numerical bracket if required. Find the full input interval consistent with the bounded reading uncertainty and compare it with a symmetric input-error assumption. Do not assume that the supplied calibration points prove monotonicity outside this branch.

Explore the question data

DRIFT VELOCITY against q

Predictions from the supplied Negative-energy components and Zitterbewegung model, with all other inputs fixed as stated in the question. 5 calibration points support a trend comparison; they are not measured data or a straight-line law.

Vertical axis (y): DRIFT VELOCITY / c

DRIFT VELOCITY against qPredictions from the supplied Negative-energy components and Zitterbewegung model, with all other inputs fixed as stated in the question. 5 calibration points support a trend comparison; they are not measured data or a straight-line law. Horizontal axis: q: Momentum p / mc. Vertical axis: DRIFT VELOCITY / c. Exact values are available in the data table below.0.1670.26960.37220.47480.57740.22420.86211.52.1382.776xy

Horizontal axis (x): q: Momentum p / mc

Axes show the supplied calibration interval; the origin may lie outside the plot.

Reading 1q: Momentum p / mc: 0.3187DRIFT VELOCITY / c: 0.1822
Read the exact data as a table
DRIFT VELOCITY against q — supplied values
q: Momentum p / mcDRIFT VELOCITY / c
0.318750.18221711748636246
0.90937499999999990.40367271724603726
1.50.4992301766027062
2.09062499999999970.541266840295488
2.681250.5621735795698961

These are the supplied values, not a worked answer. No fitted trend is drawn.

Show a hint

The model's endpoint predictions bracket the reading. Solve Y_out(q) = Y_measured and then solve again at Y_measured − u and Y_measured + u; a nonlinear inverse can give unequal input deviations.

Reveal the worked solution
  1. Endpoint predictions are Y_out(0.31875) = 0.182217 c and Y_out(2.68125) = 0.562174 c. The target lies between them. On the selected branch the response is strictly increasing.
  2. Endpoint interpolation gives q_line = 0.31875 + (0.521665 − 0.182217)(2.68125 − 0.31875)/(0.562174 − 0.182217) = 2.42937 mc. This estimates the inverse of the endpoint chord, rather than assuming that chord is the physical model.
  3. Use the actual equation, not the chord, and bracket its root. Evaluating the midpoint and retaining the sign-changing half-interval gives q = 1.75987 mc; substituting back gives Y_out = 0.52166489 c.
  4. The interpolation error in the input is |q_line − q| = 0.669499 mc. This is larger than the requested numerical tolerance.
  5. Invert both limits Y_out = 0.51558859 c and 0.52774119 c. Sort the roots even if the response decreases: 1.68016 mc ≤ q ≤ 1.84872 mc.
  6. The lower and upper input deviations are 0.0797148 mc and 0.0888439 mc. They differ because the inverse response is nonlinear; a symmetric error bar would only be a local approximation. This interval includes reading uncertainty only, not model error or uncertainty in fixed inputs.

Answer q ≈ 1.75987 mc; compatible interval [1.68016, 1.84872] mc. Endpoint interpolation gives 2.42937 mc.

Compare the assumptions, method and units with your attempt. Explain any correction before trying another question.

Extra challenging · 02

The Klein step: transmission through a step higher than 2mc²: infer an input with uncertainty

The Klein step is a sharp electrostatic step met by a Dirac electron; once the step exceeds E + mc², the electron's energy lies in the far side's negative-energy continuum and a finite fraction of its flux is transmitted, however high the step. Use the supplied topic model Y_out = 4*sqrt(max((k-v)*(k+2)*(k+2-v)*k,0))/(sqrt(abs((k-v)*(k+2)))+sqrt(abs((k+2-v)*k)))^2, where Y_out is T, q BY GROUP VELOCITY (dimensionless). The output symbol Y_out is distinct from all input symbols. ^ denotes a power; π = 3.141592653589793; sqrt(x) means √x; abs(x) means |x|; max(a,b) selects the larger value. v is Step height V₀, measured in mc². Keep k (Kinetic energy E − mc²) = 1 mc² fixed. For this exercise use only 3.16 mc² ≤ v ≤ 7.15 mc² on this continuous branch. A grid of test settings is an input restriction for this exercise, not by itself a claim of physical quantisation. This calculation retains the topic model's stated assumptions; it does not extend them to a different physical setup. A reading gives Y_out = 0.84408556, with a bounded reading uncertainty ±0.0043123. Treat the model and all fixed inputs as exact for this exercise; the bound is not a standard deviation. First estimate the unknown v by straight-line interpolation between the two endpoint predictions. Then solve the actual model on the stated branch to an absolute input tolerance of 0.000399 mc², using a numerical bracket if required. Find the full input interval consistent with the bounded reading uncertainty and compare it with a symmetric input-error assumption. Do not assume that the supplied calibration points prove monotonicity outside this branch.

Explore the question data

T, q BY GROUP VELOCITY against v

Predictions from the supplied The Klein step: transmission through a step higher than 2mc² model, with all other inputs fixed as stated in the question. 5 calibration points support a trend comparison; they are not measured data or a straight-line law.

Vertical axis (y): T, q BY GROUP VELOCITY

T, q BY GROUP VELOCITY against vPredictions from the supplied The Klein step: transmission through a step higher than 2mc² model, with all other inputs fixed as stated in the question. 5 calibration points support a trend comparison; they are not measured data or a straight-line law. Horizontal axis: v: Step height V₀ / mc². Vertical axis: T, q BY GROUP VELOCITY. Exact values are available in the data table below.0.45330.56220.67110.780.88934.0785.1556.2327.31xy

Horizontal axis (x): v: Step height V₀ / mc²

Axes show the supplied calibration interval; the origin may lie outside the plot.

Reading 1v: Step height V₀ / mc²: 3.16T, q BY GROUP VELOCITY: 0.4694
Read the exact data as a table
T, q BY GROUP VELOCITY against v — supplied values
v: Step height V₀ / mc²T, q BY GROUP VELOCITY
3.160.4694231892919099
4.15750.7677157339571586
5.1549999999999990.8297318790176073
6.15250.8572776702469411
7.14999999999999950.8728340170656331

These are the supplied values, not a worked answer. No fitted trend is drawn.

Show a hint

The model's endpoint predictions bracket the reading. Solve Y_out(v) = Y_measured and then solve again at Y_measured − u and Y_measured + u; a nonlinear inverse can give unequal input deviations.

Reveal the worked solution
  1. Endpoint predictions are Y_out(3.16) = 0.469423 and Y_out(7.15) = 0.872834. The target lies between them. On the selected branch the response is strictly increasing.
  2. Endpoint interpolation gives v_line = 3.16 + (0.844086 − 0.469423)(7.15 − 3.16)/(0.872834 − 0.469423) = 6.86566 mc². This estimates the inverse of the endpoint chord, rather than assuming that chord is the physical model.
  3. Use the actual equation, not the chord, and bracket its root. Evaluating the midpoint and retaining the sign-changing half-interval gives v = 5.5939 mc²; substituting back gives Y_out = 0.84408556.
  4. The interpolation error in the input is |v_line − v| = 1.27176 mc². This is larger than the requested numerical tolerance.
  5. Invert both limits Y_out = 0.83977326 and 0.84839786. Sort the roots even if the response decreases: 5.44719 mc² ≤ v ≤ 5.75634 mc².
  6. The lower and upper input deviations are 0.146714 mc² and 0.162438 mc². They differ because the inverse response is nonlinear; a symmetric error bar would only be a local approximation. This interval includes reading uncertainty only, not model error or uncertainty in fixed inputs.

Answer v ≈ 5.5939 mc²; compatible interval [5.44719, 5.75634] mc². Endpoint interpolation gives 6.86566 mc².

Compare the assumptions, method and units with your attempt. Explain any correction before trying another question.

Extra challenging · 03

Dirac's sea and the radiative catastrophe: infer an input with uncertainty

Dirac's sea is hole theory's vacuum: the state in which every negative-energy solution of the free Dirac equation, two spin states for each momentum, holds one electron, so that the Pauli principle forbids any further electron from radiating into the negative-energy continuum. Use the supplied topic model Y_out = wp+sqrt(1-1/137.036^2)-sqrt(wp^2-1), where Y_out is MOMENTUM TO NUCLEUS ≥ (mc). The output symbol Y_out is distinct from all input symbols. ^ denotes a power; π = 3.141592653589793; sqrt(x) means √x. wp is Landing energy |E′| of the electron, measured in mc². Keep we (Sea empty down to |E| (filled below)) = 2 mc² fixed. For this exercise use only 1.53125 mc² ≤ wp ≤ 5.46875 mc² on this continuous branch. A grid of test settings is an input restriction for this exercise, not by itself a claim of physical quantisation. This calculation retains the topic model's stated assumptions; it does not extend them to a different physical setup. A reading gives Y_out = 1.1292224 mc, with a bounded reading uncertainty ±0.0055565 mc. Treat the model and all fixed inputs as exact for this exercise; the bound is not a standard deviation. First estimate the unknown wp by straight-line interpolation between the two endpoint predictions. Then solve the actual model on the stated branch to an absolute input tolerance of 0.00039375 mc², using a numerical bracket if required. Find the full input interval consistent with the bounded reading uncertainty and compare it with a symmetric input-error assumption. Do not assume that the supplied calibration points prove monotonicity outside this branch.

Explore the question data

MOMENTUM TO NUCLEUS ≥ against wp

Predictions from the supplied Dirac's sea and the radiative catastrophe model, with all other inputs fixed as stated in the question. 5 calibration points support a trend comparison; they are not measured data or a straight-line law.

Vertical axis (y): MOMENTUM TO NUCLEUS ≥ / mc

MOMENTUM TO NUCLEUS ≥ against wpPredictions from the supplied Dirac's sea and the radiative catastrophe model, with all other inputs fixed as stated in the question. 5 calibration points support a trend comparison; they are not measured data or a straight-line law. Horizontal axis: wp: Landing energy |E′| of the electron / mc². Vertical axis: MOMENTUM TO NUCLEUS ≥ / mc. Exact values are available in the data table below.1.0811.1561.2321.3071.3831.3742.4373.54.5635.626xy

Horizontal axis (x): wp: Landing energy |E′| of the electron / mc²

Axes show the supplied calibration interval; the origin may lie outside the plot.

Reading 1wp: Landing energy |E′| of the electron / mc²: 1.531MOMENTUM TO NUCLEUS ≥ / mc: 1.372
Read the exact data as a table
MOMENTUM TO NUCLEUS ≥ against wp — supplied values
wp: Landing energy |E′| of the electron / mc²MOMENTUM TO NUCLEUS ≥ / mc
1.531251.371599916763656
2.5156251.207272326672661
3.51.1458714077189383
4.4843751.112893336955775
5.468751.0921792643662478

These are the supplied values, not a worked answer. No fitted trend is drawn.

Show a hint

The model's endpoint predictions bracket the reading. Solve Y_out(wp) = Y_measured and then solve again at Y_measured − u and Y_measured + u; a nonlinear inverse can give unequal input deviations.

Reveal the worked solution
  1. Endpoint predictions are Y_out(1.53125) = 1.3716 mc and Y_out(5.46875) = 1.09218 mc. The target lies between them. On the selected branch the response is strictly decreasing.
  2. Endpoint interpolation gives wp_line = 1.53125 + (1.12922 − 1.3716)(5.46875 − 1.53125)/(1.09218 − 1.3716) = 4.94675 mc². This estimates the inverse of the endpoint chord, rather than assuming that chord is the physical model.
  3. Use the actual equation, not the chord, and bracket its root. Evaluating the midpoint and retaining the sign-changing half-interval gives wp = 3.93313 mc²; substituting back gives Y_out = 1.1292224 mc.
  4. The interpolation error in the input is |wp_line − wp| = 1.01363 mc². This is larger than the requested numerical tolerance.
  5. Invert both limits Y_out = 1.1236659 mc and 1.1347789 mc. Sort the roots even if the response decreases: 3.77645 mc² ≤ wp ≤ 4.10413 mc².
  6. The lower and upper input deviations are 0.156676 mc² and 0.171002 mc². They differ because the inverse response is nonlinear; a symmetric error bar would only be a local approximation. This interval includes reading uncertainty only, not model error or uncertainty in fixed inputs.

Answer wp ≈ 3.93313 mc²; compatible interval [3.77645, 4.10413] mc². Endpoint interpolation gives 4.94675 mc².

Compare the assumptions, method and units with your attempt. Explain any correction before trying another question.

Extra challenging · 04

Holes as positrons: charge, energy, momentum and spin: infer an input with uncertainty

A hole is an empty negative-energy state in the filled Dirac sea. Measured from the vacuum it carries minus the charge, energy, momentum and angular momentum of the missing electron, which makes it a particle of charge +e, positive energy, mass m and spin ½: the positron. Use the supplied topic model Y_out = 0.511*(sqrt(1+p^2)+1), where Y_out is MIN PAIR COST |E| + mc² (MeV). The output symbol Y_out is distinct from all input symbols. ^ denotes a power; π = 3.141592653589793; sqrt(x) means √x. p is Emptied state's momentum p_z, measured in mc. Keep sz (Emptied state's spin J_z) = 1 ħ/2 fixed. For this exercise use only -1.575 mc ≤ p ≤ -0.105 mc on this continuous branch. A grid of test settings is an input restriction for this exercise, not by itself a claim of physical quantisation. This calculation retains the topic model's stated assumptions; it does not extend them to a different physical setup. A reading gives Y_out = 1.1284629 MeV, with a bounded reading uncertainty ±0.015548 MeV. Treat the model and all fixed inputs as exact for this exercise; the bound is not a standard deviation. First estimate the unknown p by straight-line interpolation between the two endpoint predictions. Then solve the actual model on the stated branch to an absolute input tolerance of 0.000147 mc, using a numerical bracket if required. Find the full input interval consistent with the bounded reading uncertainty and compare it with a symmetric input-error assumption. Do not assume that the supplied calibration points prove monotonicity outside this branch.

Explore the question data

MIN PAIR COST |E| + mc² against p

Predictions from the supplied Holes as positrons: charge, energy, momentum and spin model, with all other inputs fixed as stated in the question. 5 calibration points support a trend comparison; they are not measured data or a straight-line law.

Vertical axis (y): MIN PAIR COST |E| + mc² / MeV

MIN PAIR COST |E| + mc² against pPredictions from the supplied Holes as positrons: charge, energy, momentum and spin model, with all other inputs fixed as stated in the question. 5 calibration points support a trend comparison; they are not measured data or a straight-line law. Horizontal axis: p: Emptied state's momentum p_z / mc. Vertical axis: MIN PAIR COST |E| + mc² / MeV. Exact values are available in the data table below.1.0071.1261.2451.3631.482-1.634-1.237-0.84-0.4431-0.0462xy

Horizontal axis (x): p: Emptied state's momentum p_z / mc

Axes show the supplied calibration interval; the origin may lie outside the plot.

Reading 1p: Emptied state's momentum p_z / mc: -1.575MIN PAIR COST |E| + mc² / MeV: 1.464
Read the exact data as a table
MIN PAIR COST |E| + mc² against p — supplied values
p: Emptied state's momentum p_z / mcMIN PAIR COST |E| + mc² / MeV
-1.5751.4643437368677679
-1.20751.312155481823753
-0.841.1783589570838169
-0.47249999999999991.0761706779869686
-0.104999999999999981.024809165960476

These are the supplied values, not a worked answer. No fitted trend is drawn.

Show a hint

The model's endpoint predictions bracket the reading. Solve Y_out(p) = Y_measured and then solve again at Y_measured − u and Y_measured + u; a nonlinear inverse can give unequal input deviations.

Reveal the worked solution
  1. Endpoint predictions are Y_out(-1.575) = 1.46434 MeV and Y_out(-0.105) = 1.02481 MeV. The target lies between them. On the selected branch the response is strictly decreasing.
  2. Endpoint interpolation gives p_line = -1.575 + (1.12846 − 1.46434)(-0.105 − -1.575)/(1.02481 − 1.46434) = -0.451664 mc. This estimates the inverse of the endpoint chord, rather than assuming that chord is the physical model.
  3. Use the actual equation, not the chord, and bracket its root. Evaluating the midpoint and retaining the sign-changing half-interval gives p = -0.6783 mc; substituting back gives Y_out = 1.1284629 MeV.
  4. The interpolation error in the input is |p_line − p| = 0.226636 mc. This is larger than the requested numerical tolerance.
  5. Invert both limits Y_out = 1.1129149 MeV and 1.1440109 MeV. Sort the roots even if the response decreases: -0.731128 mc ≤ p ≤ -0.622483 mc.
  6. The lower and upper input deviations are 0.052828 mc and 0.0558169 mc. They differ because the inverse response is nonlinear; a symmetric error bar would only be a local approximation. This interval includes reading uncertainty only, not model error or uncertainty in fixed inputs.

Answer p ≈ -0.6783 mc; compatible interval [-0.731128, -0.622483] mc. Endpoint interpolation gives -0.451664 mc.

Compare the assumptions, method and units with your attempt. Explain any correction before trying another question.

Choose a topic