research30 min studyLesson 1000 in Physics

Writing a reproducible physics research report

Build a complete quantitative claim from pendulum timing, units, cycle counts, uncertainty and a controlled amplitude correction.

On this page

What you’ll understand

  • Reconstruct a gravity estimate from raw elapsed time and cycle count.
  • Compare propagated measurement uncertainty with a finite-amplitude model correction.
  • Identify metadata errors that ordinary uncertainty propagation does not represent.

A reproducible report permits another researcher to reconstruct the claim from its inputs and assumptions, including a route to finding mistakes. For a pendulum measurement, publishing only a fitted value of gravitational acceleration is insufficient. The elapsed time, number of full cycles, effective length, release amplitude, unit conventions and timing procedure determine what was actually inferred. A report should connect those records to the model and then to the uncertainty and conclusion. This is a physics argument with an auditable evidence chain, rather than a list of writing preferences.

Preserve the quantities that determine the estimator

Use an ideal simple pendulum: a point bob on a massless inextensible support, a fixed frictionless pivot and uniform gravitational acceleration. The length runs from pivot to the bob’s center of mass. Let tau be the elapsed time for an explicitly counted integer number of complete cycles. The period is tau divided by that count. In the small-angle approximation, gravity is proportional to length and to the square of count divided by elapsed time.

A raw-data record should preserve those inputs rather than only a computed period. Count a cycle between equivalent phase crossings so that a half swing cannot be mistaken for a full oscillation. Record the reference for length and any timebase calibration. Under independent small input uncertainties, logarithmic differentiation gives the leading relative uncertainty of the gravity estimate. A wrong cycle count is a discrete metadata error, not a tiny continuous fluctuation to include in that formula.

T=τ/N,g0=4π2L(N/τ)2,u(g0)g0≃[u(L)L]2+4[u(τ)τ]2T=\tau/N,\qquad g_0=4\pi^2L(N/\tau)^2,\qquad \frac{u(g_0)}{g_0}\simeq\sqrt{\left[\frac{u(L)}{L}\right]^2+4\left[\frac{u(\tau)}{\tau}\right]^2}
L is in m, τ and T in s, and g0 in m/s². N is an exactly known full-cycle count. The uncertainty expression assumes independent length and elapsed-time inputs, small uncertainties and a valid small-angle model; it excludes amplitude and other model corrections.↔ Scroll wide equations horizontally.

Quantify a model correction before calling the result precise

The exact undamped pendulum period follows from energy conservation and a one-dimensional integral. Substituting sine of half the turning angle converts it into a complete elliptic integral. Expanding the integrand at small turning angle and averaging its sine-squared factor gives a period increase of angle squared over sixteen. This derives the leading finite-amplitude effect from the same ideal pendulum mechanics rather than assigning an arbitrary systematic allowance.

Because the measured finite-amplitude period is longer than the small-angle period, inserting it directly into the small-angle gravity formula underestimates gravity. At fixed measured period, the leading corrected estimate therefore increases by angle squared over eight in fractional terms. The sign matters. A corrected report should state the retained order and account for amplitude uncertainty if relevant. Air drag, a finite bob, support elasticity and pivot motion are other potential model changes with different observable signatures.

T=2πL/g[1+θ0216+O(θ04)],g=g0[1+θ028+O(θ04)]T=2\pi\sqrt{L/g}\left[1+\frac{\theta_0^2}{16}+O(\theta_0^4)\right],\qquad g=g_0\left[1+\frac{\theta_0^2}{8}+O(\theta_0^4)\right]
θ0 is the turning-angle amplitude in radians. The expansion concerns an ideal undamped pendulum at small amplitude. The second expression infers g at fixed observed T; it is not a claim that increasing amplitude physically changes gravity.↔ Scroll wide equations horizontally.

Make the conclusion traceable to the reported calculation

An auditable result states the estimated quantity, its units, uncertainty convention, applied corrections and validity range. Separate illustrative data from actual observations and name the primary source of the physical model. If code performs the reduction, preserve its version, constants, unit conversion and input data together with sufficient arithmetic that a short case can be checked by hand. A reproducible numerical answer need not be physically accurate if the measurement model or metadata is wrong.

The final claim should distinguish a derivation of the amplitude correction from an empirical gravity estimate and from a proposal to improve the apparatus. A comparison run at another length tests the predicted square-root period scaling and can reveal count or length errors. A small residual then has a concrete interpretation rather than serving as decoration. The worked example shows a correction larger than the quoted measurement uncertainty, while the practice example demonstrates a metadata mistake that no increase in numerical precision could repair.

Report measurement uncertainty and amplitude bias separately

For an illustrative run, L=1.000 m with standard uncertainty 0.002 m, and 20 full cycles take τ=40.14 s with standard uncertainty 0.04 s. Treat these inputs as independent. The turning amplitude is θ0=0.20 rad, known sufficiently well for this exercise. Use 4π²=39.4784. Calculate the small-angle g0, its propagated standard uncertainty, and the leading amplitude correction.

  1. The period is 40.14/20=2.007 s. Thus g0=39.4784/(2.007)²≈9.8009 m/s². The result follows from the stated effective length and full-cycle count.

  2. The relative propagated standard uncertainty is √[(0.002/1.000)²+(2×0.04/40.14)²]≈0.0028235. Multiplying gives u(g0)≈0.0277 m/s², so the small-angle reduction would report about 9.801 with standard uncertainty 0.028 m/s².

  3. The fractional leading amplitude correction is θ0²/8=0.005. It raises the inferred gravity by about 0.0490 m/s², giving 9.850 m/s² to the retained order. The correction exceeds the propagated measurement uncertainty and should be included with its model assumptions.

The result

The toy run gives g0≈9.801±0.028 m/s² as a standard-uncertainty small-angle result; the leading amplitude-corrected estimate is about 9.850 m/s². The uncertainty quoted first does not by itself cover omitted model effects.

Explore the model

Interactive lab

Interactive lab · function

Check the report’s period-length scaling

At fixed gravity and amplitude, normalized length is proportional to period squared. Reset and set both controls to Coefficient a=1 and Position x=1, interpreting x as normalized period; predict y=1. For a fresh comparison, Reset and set both controls to a=1 and x=0.5; predict y=0.25, corresponding to one quarter of the length.

Check the report’s period-length scalingThis lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Current readings: y at chosen x: 0.608; Tangent slope · 2ax: 1.35; Signed area · ax³/3: 0.182. All parameters are adjustable with the labeled controls below.−2.4−1.201.22.4−9−4.504.59xy(x, ax²)
y = ax²TangentSigned integral
y at chosen x
0.608
Tangent slope · 2ax
1.35
Signed area · ax³/3
0.182
t = 0 s

This lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forThis quadratic consistency check has no timed crossings, cycle counter, finite-amplitude pendulum dynamics, gravity inference or uncertainty propagation. The source records and physical corrections must still be supplied in the report.

Check your understanding

A run’s true full-cycle count is 10, but its elapsed time is divided by 20. What happens to the small-angle gravity estimate?

Try it yourself

A second illustrative run has L=0.2500 m and τ=10.035 s for 10 full cycles, at the same amplitude as the first. Calculate its small-angle g0 and compare its period with the first run. Then calculate the erroneous g0 if its count were recorded as 20. Use 4π²=39.4784.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

The length is one quarter of the first length. Keep elapsed time and cycle count separate before forming the period.

Compare with the solution

The correct period is 10.035/10=1.0035 s, exactly half of 2.007 s. Since L is one quarter and T² is one quarter, g0 is unchanged: 39.4784×0.2500/(1.0035)²≈9.8009 m/s². This is consistent with the fixed-amplitude square-root length scaling. Recording 20 cycles would instead give T=0.50175 s and g0≈39.2035 m/s², four times the correct estimate. The shared amplitude correction would affect both correct runs similarly; it cannot explain a factor-four discrepancy caused by cycle metadata.

Take it with you

  • A reproducible claim links raw inputs, metadata, units, model and inference.
  • A model correction can exceed propagated instrument uncertainty and must be stated separately.
  • Reconstruction and limiting checks can reveal mistakes that more numerical digits cannot fix.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.