Skip to main content
Course and syllabus guidance

University Quantum Field Theory I · unit UQFTI-04

The S-Matrix from Correlation Functions: Kallen-Lehmann and LSZ

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

In and out states, S = 1 + iT and the invariant amplitude M: infer an input with uncertainty

The invariant amplitude M is what an S-matrix element leaves once the no-scattering term and the momentum-conserving factor are removed: ⟨β|S − 1|α⟩ = (2π)⁴ δ⁴(p_α − p_β) iM(α → β). Use the supplied topic model Y_out = (exp(yc+eta)+exp(-yc-eta))/2, where Y_out is E₂ / m IN THIS FRAME (dimensionless). The output symbol Y_out is distinct from all input symbols. ^ denotes a power; π = 3.141592653589793; exp(x) means eˣ. eta is Observer's boost rapidity η, dimensionless. Keep yc (CM rapidities ±Y of the pair) = 1.4 fixed. For this exercise use only -1.26 ≤ eta ≤ 2.3625 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 = 5.289041, with a bounded reading uncertainty ±0.64188. Treat the model and all fixed inputs as exact for this exercise; the bound is not a standard deviation. First estimate the unknown eta 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.00036225, 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

E₂ / m IN THIS FRAME against eta

Predictions from the supplied In and out states, S = 1 + iT and the invariant amplitude M 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): E₂ / m IN THIS FRAME

E₂ / m IN THIS FRAME against etaPredictions from the supplied In and out states, S = 1 + iT and the invariant amplitude M 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: eta: Observer's boost rapidity η. Vertical axis: E₂ / m IN THIS FRAME. Exact values are available in the data table below.0.18865.73211.2716.8222.36-1.405-0.42680.55121.5292.507xy

Horizontal axis (x): eta: Observer's boost rapidity η

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

Reading 1eta: Observer's boost rapidity η: -1.26E₂ / m IN THIS FRAME: 1.01
Read the exact data as a table
E₂ / m IN THIS FRAME against eta — supplied values
eta: Observer's boost rapidity ηE₂ / m IN THIS FRAME
-1.26000000000000021.0098160171280166
-0.35437500000000021.5983241419312486
0.55124999999999983.589787636809125
1.45687499999999978.73224642129301
2.362521.53957848659048

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(eta) = 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.26) = 1.00982 and Y_out(2.3625) = 21.5396. The target lies between them. On the selected branch the response is strictly increasing.
  2. Endpoint interpolation gives eta_line = -1.26 + (5.28904 − 1.00982)(2.3625 − -1.26)/(21.5396 − 1.00982) = -0.504926. 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 eta = 0.949725; substituting back gives Y_out = 5.289041.
  4. The interpolation error in the input is |eta_line − eta| = 1.45465. This is larger than the requested numerical tolerance.
  5. Invert both limits Y_out = 4.647161 and 5.930921. Sort the roots even if the response decreases: 0.817621 ≤ eta ≤ 1.06614.
  6. The lower and upper input deviations are 0.132104 and 0.116418. 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 eta ≈ 0.949725; compatible interval [0.817621, 1.06614]. Endpoint interpolation gives -0.504926.

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

Extra challenging · 02

The Källén-Lehmann representation: infer an input with uncertainty

The Källén-Lehmann representation writes the exact two-point function of a scalar field as an integral over mass μ of free Feynman propagators, weighted by a spectral density ρ(μ²) ≥ 0 that measures how strongly φ(0) creates states of invariant mass μ from the vacuum. Use the supplied topic model Y_out = gm^2*r^2/(4*pi^2)*(atan(1/sqrt(4*r^2-1))/sqrt(4*r^2-1) - 1/(4*r^2)), where Y_out is ψψ̄ WEIGHT ∫ρ ds/2π (dimensionless). The output symbol Y_out is distinct from all input symbols. ^ denotes a power; π = 3.141592653589793; sqrt(x) means √x; unmarked trigonometric functions use radians. r is Mass ratio M/m, dimensionless. Keep gm (Coupling g/m) = 2 fixed. For this exercise use only 0.650938 ≤ r ≤ 1.39906 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.0041248765, with a bounded reading uncertainty ±0.00025807. Treat the model and all fixed inputs as exact for this exercise; the bound is not a standard deviation. First estimate the unknown r 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.0000748125, 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

ψψ̄ WEIGHT ∫ρ ds/2π against r

Predictions from the supplied The Källén-Lehmann representation 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): ψψ̄ WEIGHT ∫ρ ds/2π

ψψ̄ WEIGHT ∫ρ ds/2π against rPredictions from the supplied The Källén-Lehmann representation 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: r: Mass ratio M/m. Vertical axis: ψψ̄ WEIGHT ∫ρ ds/2π. Exact values are available in the data table below.0.0017090.0064010.011090.015780.020480.6210.8231.0251.2271.429xy

Horizontal axis (x): r: Mass ratio M/m

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

Reading 1r: Mass ratio M/m: 0.6509ψψ̄ WEIGHT ∫ρ ds/2π: 0.01978
Read the exact data as a table
ψψ̄ WEIGHT ∫ρ ds/2π against r — supplied values
r: Mass ratio M/mψψ̄ WEIGHT ∫ρ ds/2π
0.65093750000000010.019780563289828426
0.837968750.008492480110232003
1.0250.004981394343232566
1.212031250.0033323820519804742
1.39906249999999990.00240442944847799

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(r) = 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.650938) = 0.0197806 and Y_out(1.39906) = 0.00240443. The target lies between them. On the selected branch the response is strictly decreasing.
  2. Endpoint interpolation gives r_line = 0.650938 + (0.00412488 − 0.0197806)(1.39906 − 0.650938)/(0.00240443 − 0.0197806) = 1.32499. 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 r = 1.10729; substituting back gives Y_out = 0.0041248765.
  4. The interpolation error in the input is |r_line − r| = 0.217695. This is larger than the requested numerical tolerance.
  5. Invert both limits Y_out = 0.0038668065 and 0.0043829465. Sort the roots even if the response decreases: 1.07981 ≤ r ≤ 1.13765.
  6. The lower and upper input deviations are 0.0274834 and 0.0303579. 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 r ≈ 1.10729; compatible interval [1.07981, 1.13765]. Endpoint interpolation gives 1.32499.

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

Extra challenging · 03

Physical mass and field-strength renormalisation Z: infer an input with uncertainty

The physical mass m is the position of the exact propagator's one-particle pole, and the field-strength renormalisation Z is its residue, Z = |⟨Ω|φ(0)|p⟩|², the share of the field's total spectral weight that goes into making a single particle. Use the supplied topic model Y_out = gm^2/(16*pi^2)*(4*asin(r/2)/(r*sqrt(4-r^2))-1)/r^2, where Y_out is −dΣ/dp² AT p² = m² (dimensionless). The output symbol Y_out is distinct from all input symbols. ^ denotes a power; π = 3.141592653589793; sqrt(x) means √x; unmarked trigonometric functions use radians. r is Mass ratio m/M, dimensionless. Keep gm (Coupling g/M) = 6 fixed. For this exercise use only 0.564688 ≤ r ≤ 1.78531 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.058809071, with a bounded reading uncertainty ±0.0027318. Treat the model and all fixed inputs as exact for this exercise; the bound is not a standard deviation. First estimate the unknown r 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.000122062, 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

−dΣ/dp² AT p² = m² against r

Predictions from the supplied Physical mass and field-strength renormalisation Z 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): −dΣ/dp² AT p² = m²

−dΣ/dp² AT p² = m² against rPredictions from the supplied Physical mass and field-strength renormalisation Z 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: r: Mass ratio m/M. Vertical axis: −dΣ/dp² AT p² = m². Exact values are available in the data table below.0.037240.059910.082590.10530.12790.51590.84541.1751.5051.834xy

Horizontal axis (x): r: Mass ratio m/M

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

Reading 1r: Mass ratio m/M: 0.5647−dΣ/dp² AT p² = m²: 0.0406
Read the exact data as a table
−dΣ/dp² AT p² = m² against r — supplied values
r: Mass ratio m/M−dΣ/dp² AT p² = m²
0.56468750.04059680288487651
0.869843750.0448667659309379
1.17499999999999980.05298261520393529
1.48015624999999980.07013520088792928
1.78531249999999990.12458649026796594

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(r) = 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.564688) = 0.0405968 and Y_out(1.78531) = 0.124586. The target lies between them. On the selected branch the response is strictly increasing.
  2. Endpoint interpolation gives r_line = 0.564688 + (0.0588091 − 0.0405968)(1.78531 − 0.564688)/(0.124586 − 0.0405968) = 0.829367. 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 r = 1.30927; substituting back gives Y_out = 0.058809071.
  4. The interpolation error in the input is |r_line − r| = 0.479902. This is larger than the requested numerical tolerance.
  5. Invert both limits Y_out = 0.056077271 and 0.061540871. Sort the roots even if the response decreases: 1.25189 ≤ r ≤ 1.35901.
  6. The lower and upper input deviations are 0.0573796 and 0.0497367. 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 r ≈ 1.30927; compatible interval [1.25189, 1.35901]. Endpoint interpolation gives 0.829367.

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

Extra challenging · 04

The LSZ reduction formula derived: infer an input with uncertainty

The LSZ reduction formula writes every S-matrix element as the vacuum correlation function ⟨Ω|Tφ⋯φ|Ω⟩ with each external leg acted on by (□ + m²) at the physical mass, multiplied by i/√Z and Fourier-transformed with that particle's momentum. Use the supplied topic model Y_out = (1-z)*(1+t^2)^(-0.75), where Y_out is ITS ENVELOPE (dimensionless). The output symbol Y_out is distinct from all input symbols. ^ denotes a power; π = 3.141592653589793. t is Time mt after the probe state, dimensionless. Keep z (Residue Z of the one-particle pole) = 0.7; gap (Gap above the pole, (μₜₕ − m)/m) = 1 fixed. For this exercise use only 1.59375 ≤ t ≤ 13.4062 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.011383184, with a bounded reading uncertainty ±0.00079453. Treat the model and all fixed inputs as exact for this exercise; the bound is not a standard deviation. First estimate the unknown t 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.00118125, 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

ITS ENVELOPE against t

Predictions from the supplied The LSZ reduction formula derived 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): ITS ENVELOPE

ITS ENVELOPE against tPredictions from the supplied The LSZ reduction formula derived 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: t: Time mt after the probe state. Vertical axis: ITS ENVELOPE. Exact values are available in the data table below.0.001680.031420.061160.090910.12061.1214.3117.510.6913.88xy

Horizontal axis (x): t: Time mt after the probe state

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

Reading 1t: Time mt after the probe state: 1.594ITS ENVELOPE: 0.1162
Read the exact data as a table
ITS ENVELOPE against t — supplied values
t: Time mt after the probe stateITS ENVELOPE
1.593750.11624257289941799
4.5468750.029865230221047704
7.4999999999999990.014414169889135717
10.4531249999999980.00881626542198758
13.4062499999999980.006086295069239051

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(t) = 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.59375) = 0.116243 and Y_out(13.4062) = 0.0060863. The target lies between them. On the selected branch the response is strictly decreasing.
  2. Endpoint interpolation gives t_line = 1.59375 + (0.0113832 − 0.116243)(13.4062 − 1.59375)/(0.0060863 − 0.116243) = 12.8382. 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 t = 8.79937; substituting back gives Y_out = 0.011383184.
  4. The interpolation error in the input is |t_line − t| = 4.03887. This is larger than the requested numerical tolerance.
  5. Invert both limits Y_out = 0.010588654 and 0.012177714. Sort the roots even if the response decreases: 8.40723 ≤ t ≤ 9.23971.
  6. The lower and upper input deviations are 0.392144 and 0.440333. 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 t ≈ 8.79937; compatible interval [8.40723, 9.23971]. Endpoint interpolation gives 12.8382.

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

Choose a topic