Physics · Book 3 · Bachelor Year 1

University Physics — Year 1

University Physics — Year 1 · Bachelor Year 1

1Units, Dimensions, and Measurement

In a basement laboratory a pendulum swings, a stopwatch clicks, and a student writes g=9.77m/s2g = 9.77\,\mathrm{m}/\mathrm{s}^{2}. Is that the “right” value? The textbook says 9.819.81. Whether the two numbers agree, whether the experiment was well designed, which of its steps limited it — none of this can be decided from the digits alone. Physics measures, and a measurement without its uncertainty is a rumor. This chapter sets up the common language of every chapter to come: units, dimensions, orders of magnitude, and the honest arithmetic of uncertainties.

The NIST-4 Kibble balance, one of the instruments that now realize the kilogram from Planck’s constant: a mass is weighed against an electromagnetic force measured in electrical units. Photograph: J. L. Lee, NIST (public domain).
The NIST-4 Kibble balance, one of the instruments that now realize the kilogram from Planck’s constant: a mass is weighed against an electromagnetic force measured in electrical units. Photograph: J. L. Lee, NIST (public domain).

1.1 The International System of Units

Definition 1.1 (Physical quantity and unit)

A physical quantity is a property that can be measured: a length, a duration, a mass, a current. Measuring it means comparing it with a reference quantity of the same kind, the unit: the result is a number times a unit, =1.257m\ell = 1.257\,\mathrm{m}. The number alone means nothing; the unit alone measures nothing.

Definition 1.2 (Base units of the SI)

The International System of Units (SI) rests on seven base units: the second (s\mathrm{s}), the meter (m\mathrm{m}), the kilogram (kg\mathrm{kg}), the ampere (A\mathrm{A}), the kelvin (K\mathrm{K}), the mole (mol\mathrm{mol}) and the candela (cd\mathrm{cd}). Since 2019 each is defined by fixing the numerical value of a constant of nature: the cesium hyperfine frequency ΔνCs=9192631770Hz\Delta\nu_{\mathrm{Cs}} = 9\,192\,631\,770\,\mathrm{Hz} defines the second; the speed of light c=299792458m/sc = 299\,792\,458\,\mathrm{m}/\mathrm{s} then defines the meter; the Planck constant h=6.62607015×1034Jsh = 6.626\,070\,15 \times 10^{-34}\,\mathrm{J}\,\mathrm{s} defines the kilogram; the elementary charge e=1.602176634×1019Ce = 1.602\,176\,634 \times 10^{-19}\,\mathrm{C} the ampere; the Boltzmann constant kB=1.380649×1023J/Kk_B = 1.380\,649 \times 10^{-23}\,\mathrm{J}/\mathrm{K} the kelvin; the Avogadro constant NA=6.02214076×1023mol1N_A = 6.022\,140\,76 \times 10^{23}\,\mathrm{mol}^{-1} the mole; and a fixed luminous efficacy KcdK_{\mathrm{cd}} the candela. Every other unit is a derived unit, a product of powers of base units: 1N=1kgm/s21\,\mathrm{N} = 1\,\mathrm{kg}\,\mathrm{m}/\mathrm{s}^{2}, 1J=1Nm1\,\mathrm{J} = 1\,\mathrm{N}\,\mathrm{m}, 1W=1J/s1\,\mathrm{W} = 1\,\mathrm{J}/\mathrm{s}, 1V=1W/A1\,\mathrm{V} = 1\,\mathrm{W}/\mathrm{A}, 1Pa=1N/m21\,\mathrm{Pa} = 1\,\mathrm{N}/\mathrm{m}^{2}.

The seven base units (outer ring) and the seven defining constants (inner ring). A fixed constant defines its unit only together with units already defined (dashed): c turns seconds into meters, h needs the meter and the second to define the kilogram, e needs the second to define the ampere.
The seven base units (outer ring) and the seven defining constants (inner ring). A fixed constant defines its unit only together with units already defined (dashed): cc turns seconds into meters, hh needs the meter and the second to define the kilogram, ee needs the second to define the ampere.

Remark 1.3 (Names and symbols)

Units named after people are written lowercase in full (newton, joule, pascal) and capitalized as symbols (N\mathrm{N}, J\mathrm{J}, Pa\mathrm{Pa}); capitalized in full, the word is the person. Prefixes scale by powers of ten, from p (101210^{-12}) through n, μ\mu, m, k, M, G to T (101210^{12}); the kilogram is the one base unit carrying a prefix in its name.

Example 1.4 (Unpacking a derived unit)

The volt: V=W/A=J/(As)=kgm2/(As3)\mathrm{V} = \mathrm{W}/\mathrm{A} = \mathrm{J}/(\mathrm{A}\,\mathrm{s}) = \mathrm{kg}\,\mathrm{m}^{2}/(\mathrm{A}\,\mathrm{s}^{3}). The ohm: Ω=V/A=kgm2/(A2s3)\Omega = \mathrm{V}/\mathrm{A} = \mathrm{kg}\,\mathrm{m}^{2}/(\mathrm{A}^{2}\,\mathrm{s}^{3}). The farad: F=C/V=As/V=A2s4/(kgm2)\mathrm{F} = \mathrm{C}/\mathrm{V} = \mathrm{A}\,\mathrm{s}/\mathrm{V} = \mathrm{A}^{2}\,\mathrm{s}^{4}/(\mathrm{kg}\,\mathrm{m}^{2}). Such unpacking is the safest check that a formula has been remembered correctly — the next section makes it systematic.

1.2 Dimensional analysis

Definition 1.5 (Dimension)

The dimension of a quantity XX, written [X][X], records how XX is built from the seven base quantities, independently of the units chosen: length LL, mass MM, time TT, electric current II, temperature Θ\Theta, amount of substance NN, luminous intensity JJ. A speed has dimension [v]=LT1[v] = L\,T^{-1}, a force [F]=MLT2[F] = M\,L\,T^{-2}, an energy [E]=ML2T2[E] = M\,L^2\,T^{-2}. A quantity of dimension 11 (an angle, a ratio, a refractive index) is dimensionless.

Theorem 1.6 (Principle of dimensional homogeneity)

In a physical law, the two sides of an equality have the same dimension, and every term of a sum has the dimension of the whole. The argument of an exponential, a logarithm, a sine or a cosine is dimensionless.

Proof. The laws of physics do not depend on the units humans choose. Changing the unit of length by a factor λ\lambda multiplies every quantity of dimension LaL^a \cdots by λa\lambda^{-a}: an equality between terms of different dimensions would hold in one system of units and fail in another. For a function such as exp\exp, the series 1+x+x2/2+1 + x + x^2/2 + \dots adds powers of xx of different dimensions unless xx is dimensionless.

Method 1.7 (Checking a formula)

  1. Write the dimension of every symbol.
  2. Reduce each side (each term) to a product LaMbTcL^a M^b T^c \cdots.
  3. Equal exponents on both sides: the formula may be right. Unequal: it is certainly wrong — a missing factor gg, a squared quantity that should be plain, an exponent with dimensions.

A homogeneous formula can still be wrong by a dimensionless factor (22, π\pi, 12\tfrac12): dimensions never see those.

Example 1.8 (The period of a pendulum)

A pendulum of length \ell and mass mm swings in gravity gg. Suppose its period is τ=kαmβgγ\tau = k\,\ell^\alpha m^\beta g^\gamma with kk dimensionless. Then T=LαMβ(LT2)γ=Lα+γMβT2γT = L^\alpha M^\beta (L T^{-2})^\gamma = L^{\alpha+\gamma} M^\beta T^{-2\gamma}, so β=0\beta = 0, γ=12\gamma = -\tfrac12, α=12\alpha = \tfrac12:

τ=k/g.\tau = k\sqrt{\ell/g} .

Without solving a single equation of motion: the mass drops out, and doubling the length multiplies the period by 2\sqrt2. Dynamics (Chapter 14) will supply k=2πk = 2\pi for small swings.

Proposition 1.9 (Dimensionless groups)

If a quantity YY depends on nn quantities X1,,XnX_1, \dots, X_n that involve rr independent base dimensions, the relation can be rewritten between n+1rn + 1 - r dimensionless combinations of the XiX_i and YY. When n+1r=1n + 1 - r = 1, the single combination must be a constant, and dimensional analysis gives the law up to a numerical factor.

Proof. Admitted at this level.

Remark 1.10 (Using the proposition)

The general statement (the Vaschy–Buckingham theorem) is a piece of linear algebra on the exponent vectors. What matters here is the recipe: list the relevant quantities, count dimensions, form the dimensionless groups. In Example 1.8, n=3n = 3 (,m,g\ell, m, g), r=3r = 3 (L,M,TL, M, T), so a single group τg/\tau\sqrt{g/\ell} exists and must be constant. Drag on a sphere of radius RR moving at speed vv in a fluid of density ρ\rho and viscosity η\eta involves n=4n = 4 quantities and r=3r = 3 dimensions: two groups, F/(ρv2R2)F/(\rho v^2 R^2) and ρvR/η\rho v R/\eta (the Reynolds number), and dimensions alone cannot say how the first depends on the second — experiment must.

1.3 Orders of magnitude

Definition 1.11 (Order of magnitude)

The order of magnitude of a quantity is the power of ten nearest to it: 6.4×106m6.4 \times 10^{6}\,\mathrm{m} (the Earth’s radius) has order of magnitude 10710^7; 3×103m3 \times 10^{3}\,\mathrm{m} has order 10310^3; 3.2×103m3.2 \times 10^{3}\,\mathrm{m} is closer to 10310^3 than to 10410^4 on a logarithmic scale (the boundary is 103.16\sqrt{10} \approx 3.16). Two quantities are “of the same order” when their ratio is below about 33.

Method 1.12 (Estimating)

To estimate a quantity nobody has measured for you:

  1. break it into factors you can estimate to within a factor of 22 or 33 each;
  2. estimate each factor, keeping one significant figure;
  3. multiply, keeping the power of ten and one digit.

Errors in different factors partly cancel; the result is usually right to within a factor of 33 to 1010, which is often all one needs to rule a hypothesis in or out.

Example 1.13 (The mass of the atmosphere)

Atmospheric pressure P01×105PaP_0 \approx 1 \times 10^{5}\,\mathrm{Pa} is the weight of the air column above one square meter, P0=mcolgP_0 = m_{\text{col}}\,g, so mcol105/10=1×104kgm_{\text{col}} \approx 10^5/10 = 1 \times 10^{4}\,\mathrm{kg} per square meter. The Earth’s surface is 4πR24×3×(6×106)25×1014m24\pi R^2 \approx 4 \times 3 \times (6\times10^6)^2 \approx 5 \times 10^{14}\,\mathrm{m}^{2}; the atmosphere weighs about 5×1018kg5 \times 10^{18}\,\mathrm{kg}. The accepted value is 5.1×1018kg5.1 \times 10^{18}\,\mathrm{kg} — the estimate needed only P0P_0 and RR.

Definition 1.14 (Significant figures)

The significant figures of a written number are its digits from the first non-zero one: 0.08200.0820 has three, 8.20×1028.20 \times 10^{-2} the same three, 820820 is ambiguous (two or three) unless written 8.2×1028.2 \times 10^{2} or 8.20×1028.20 \times 10^{2}. A result is written with the significant figures its uncertainty allows (Method 1.20 below): not “g=9.7723m/s2g = 9.7723\,\mathrm{m}/\mathrm{s}^{2}” when the fourth digit is unknown.

1.4 Measurement uncertainty

Definition 1.15 (Measurand, error, uncertainty)

The measurand is the quantity one intends to measure. Repeating a measurement gives scattered values: the measurement is variable. The error of one result is its unknown difference to the true value; the standard uncertainty u(x)u(x) of a result xx is the standard deviation of the values one could reasonably attribute to the measurand — a quantity one can evaluate. Its sources: the operator (reaction time, parallax), the environment (temperature, vibrations), the instrument (resolution, calibration), the method (a model that is only approximately true).

Proposition 1.16 (Type A evaluation)

From nn independent repeated values x1,,xnx_1, \dots, x_n of mean xˉ=1nxi\bar x = \frac1n \sum x_i and experimental standard deviation

s=1n1i=1n(xixˉ)2,s = \sqrt{\frac{1}{n-1}\sum_{i=1}^n (x_i - \bar x)^2},

the best estimate of the measurand is xˉ\bar x and its standard uncertainty is

u(xˉ)=sn.u(\bar x) = \frac{s}{\sqrt n} .

Partial proof. That ss estimates the spread of a single measurement is the definition of the experimental standard deviation (n1n - 1 rather than nn because the mean itself was estimated from the data). That the mean of nn independent values scatters n\sqrt n times less than one value is the addition of variances: the variance of a sum of nn independent variables is nn times the variance of one, so the variance of the mean is s2n/n2=s2/ns^2 n / n^2 = s^2/n. The probability course of this year proves both statements.

Fifty timings of the same pendulum period, binned by 0.01\, s. Their spread is s 0.016\, s (one timing); the mean of the fifty is known to s/√50 0.002\, s. The bell curve is the Gaussian of the same mean and spread.
Fifty timings of the same pendulum period, binned by 0.01s0.01\,\mathrm{s}. Their spread is s0.016ss \approx 0.016\,\mathrm{s} (one timing); the mean of the fifty is known to s/500.002ss/\sqrt{50} \approx 0.002\,\mathrm{s}. The bell curve is the Gaussian of the same mean and spread.

Proposition 1.17 (Type B evaluation)

When a value is read once, its uncertainty is evaluated from what is known of the instrument. If the true value can lie anywhere between xax - a and x+ax + a with no preference (a uniform distribution of half-width aa),

u(x)=a3.u(x) = \frac{a}{\sqrt 3} .

For a scale graduated in steps δ\delta (half-width a=δ/2a = \delta/2), u=δ/(23)=δ/12u = \delta/(2\sqrt3) = \delta/\sqrt{12}; the same applies to the last digit of a digital display. A manufacturer’s tolerance quoted as “±a\pm a” is treated the same way unless the data sheet says otherwise.

Proof. The variance of a uniform distribution on [xa,x+a][x - a, x + a] is 12aaat2 ⁣dt=12a2a33=a23\frac{1}{2a}\int_{-a}^{a} t^2\,\dd t = \frac{1}{2a}\cdot\frac{2a^3}{3} = \frac{a^2}{3}; its standard deviation is a/3a/\sqrt3.

A reading known only to lie within ± a: uniform density 1/(2a), standard deviation a/√3 0.58\,a — smaller than a, because the extremes are no likelier than the center.
A reading known only to lie within ±a\pm a: uniform density 1/(2a)1/(2a), standard deviation a/30.58aa/\sqrt3 \approx 0.58\,a — smaller than aa, because the extremes are no likelier than the center.

Theorem 1.18 (Combined uncertainty)

Let y=f(x1,,xn)y = f(x_1, \dots, x_n) be computed from independent measured values with standard uncertainties u(xi)u(x_i). Then

u(y)2=i=1n(fxi)2u(xi)2.u(y)^2 = \sum_{i=1}^n \left(\frac{\partial f}{\partial x_i}\right)^2 u(x_i)^2 .

In particular:

  • sum or difference y=x1±x2y = x_1 \pm x_2: u(y)2=u(x1)2+u(x2)2u(y)^2 = u(x_1)^2 + u(x_2)^2;
  • product or quotient y=x1x2y = x_1 x_2 or x1/x2x_1/x_2: (u(y)y)2=(u(x1)x1)2+(u(x2)x2)2\left(\dfrac{u(y)}{y}\right)^2 = \left(\dfrac{u(x_1)}{x_1}\right)^2 + \left(\dfrac{u(x_2)}{x_2}\right)^2;
  • power y=kxαy = k\,x^\alpha: u(y)y=αu(x)x\dfrac{u(y)}{\abs y} = \abs\alpha\,\dfrac{u(x)}{\abs x}.

Proof. For small deviations δxi\delta x_i, the first-order Taylor expansion gives δy=i(f/xi)δxi\delta y = \sum_i (\partial f/\partial x_i)\,\delta x_i. The deviations are independent and of zero mean, so the variance of the sum is the sum of the variances, each scaled by the square of its coefficient — the displayed formula. The special cases follow: (x1x2)/x1=x2\partial(x_1 x_2)/\partial x_1 = x_2 gives u(y)2=x22u(x1)2+x12u(x2)2u(y)^2 = x_2^2 u(x_1)^2 + x_1^2 u(x_2)^2, divide by y2=x12x22y^2 = x_1^2 x_2^2; for y=kxαy = kx^\alpha, y/x=αy/x\partial y/\partial x = \alpha y/x.

Example 1.19 (The density of a cylinder)

A cylinder: m=47.32gm = 47.32\,\mathrm{g} (u=0.01gu = 0.01\,\mathrm{g}), d=12.00mmd = 12.00\,\mathrm{mm} (u=0.02mmu = 0.02\,\mathrm{mm}), h=50.2mmh = 50.2\,\mathrm{mm} (u=0.1mmu = 0.1\,\mathrm{mm}). ρ=4m/(πd2h)\rho = 4m/(\pi d^2 h), so the relative uncertainty squared adds (0.01/47.32)2=4.5×108(0.01/47.32)^2 = 4.5 \times 10^{-8}, (2×0.02/12.00)2=1.1×105(2 \times 0.02/12.00)^2 = 1.1 \times 10^{-5} and (0.1/50.2)2=4.0×106(0.1/50.2)^2 = 4.0 \times 10^{-6}: the diameter dominates, because it is squared and measured to 0.17%0.17\%. u(ρ)/ρ=1.5×105=0.39%u(\rho)/\rho = \sqrt{1.5 \times 10^{-5}} = 0.39\%; ρ=8334kg/m3\rho = 8334\,\mathrm{kg}/\mathrm{m}^{3}, so u(ρ)=33kg/m3u(\rho) = 33\,\mathrm{kg}/\mathrm{m}^{3}: ρ=(8.33±0.03)×103 kg/m3\rho = (8.33 \pm 0.03)\times10^3\ \mathrm{kg}/\mathrm{m}^{3} — brass.

Method 1.20 (Writing a result)

  1. Round the uncertainty to one significant figure (two if the first is 11 or 22).
  2. Round the value to the same decimal place.
  3. Write value and uncertainty with the unit: g=(9.77±0.03)m/s2g = (9.77 \pm 0.03)\,\mathrm{m}/\mathrm{s}^{2}, or g=9.77m/s2g = 9.77\,\mathrm{m}/\mathrm{s}^{2}, u(g)=0.03m/s2u(g) = 0.03\,\mathrm{m}/\mathrm{s}^{2}.

Keep extra digits during the calculation; round only at the end.

Definition 1.21 (Normalized deviation)

Two results x1±u1x_1 \pm u_1 and x2±u2x_2 \pm u_2 of the same measurand are compared through their normalized deviation

z=x1x2u12+u22.z = \frac{\abs{x_1 - x_2}}{\sqrt{u_1^2 + u_2^2}} .

They are compatible when z2z \leq 2; a larger zz points to an unaccounted effect — or an underestimated uncertainty. When x2x_2 is a reference value with no quoted uncertainty, take u2=0u_2 = 0.

Remark 1.22 (Why 2)

If the two results scatter as Gaussians, their difference has standard deviation u12+u22\sqrt{u_1^2 + u_2^2}, and a Gaussian exceeds twice its standard deviation in about 5%5\% of cases. “z2z \leq 2” is a convention, not a theorem: it accepts roughly one good measurement in twenty as a false disagreement.

1.5 Fitting a model to data

Proposition 1.23 (Least-squares line)

Given nn points (xi,yi)(x_i, y_i) expected to obey y=ax+by = a x + b, the straight line minimizing the sum of squared vertical distances i(yiaxib)2\sum_i (y_i - a x_i - b)^2 has

a=i(xixˉ)(yiyˉ)i(xixˉ)2,b=yˉaxˉ.a = \frac{\sum_i (x_i - \bar x)(y_i - \bar y)}{\sum_i (x_i - \bar x)^2}, \qquad b = \bar y - a\,\bar x .

Proof. The sum is a quadratic function of (a,b)(a, b); setting both partial derivatives to zero gives two linear equations: i(yiaxib)=0\sum_i (y_i - ax_i - b) = 0, which is yˉ=axˉ+b\bar y = a \bar x + b, and ixi(yiaxib)=0\sum_i x_i (y_i - ax_i - b) = 0. Substituting bb into the second and centering the variables yields aa. The minimization itself is carried out in the mathematics course on functions of two variables.

Method 1.24 (Validating a model graphically)

  1. Choose variables that make the expected law a straight line: τ2\tau^2 against \ell for the pendulum, lny\ln y against lnx\ln x for a power law y=Cxαy = Cx^\alpha (slope α\alpha).
  2. Plot the points with their uncertainty bars.
  3. Fit the line; the model is acceptable if the bars straddle the line without a systematic trend of the residuals, and if the intercept is what the model predicts (zero, here).
  4. Read the physics off the slope, with its uncertainty.
2 against  for five pendulum lengths, with uncertainty bars; the least-squares line passes through the origin within uncertainty and its slope 4π2/g gives g.
τ2\tau^2 against \ell for five pendulum lengths, with uncertainty bars; the least-squares line passes through the origin within uncertainty and its slope 4π2/g4\pi^2/g gives gg.

Example 1.25 (Reading gg from a slope)

The five points of the figure have ˉ=0.800m\bar\ell = 0.800\,\mathrm{m}, τ2=3.232s2\overline{\tau^2} = 3.232\,\mathrm{s}^{2}, and the sums (iˉ)(τi2τ2)=1.614ms2\sum (\ell_i - \bar\ell)(\tau_i^2 - \overline{\tau^2}) = 1.614\,\mathrm{m}\,\mathrm{s}^{2}, (iˉ)2=0.400m2\sum (\ell_i - \bar\ell)^2 = 0.400\,\mathrm{m}^{2}: slope a=4.035s2/ma = 4.035\,\mathrm{s}^{2}/\mathrm{m}, intercept b=0.004s2b = 0.004\,\mathrm{s}^{2}, negligible. With τ2=(4π2/g)\tau^2 = (4\pi^2/g)\,\ell, g=4π2/a=9.78m/s2g = 4\pi^2/a = 9.78\,\mathrm{m}/\mathrm{s}^{2}.

1.6 Exercises

Exercise 1.1

Express the joule, the watt, the pascal and the coulomb in base units. Deduce that Pam3\mathrm{Pa}\,\mathrm{m}^{3} is a unit of energy.

Solution

Solution of Exercise 1.1.

J=kgm2s2\mathrm{J} = \mathrm{kg}\,\mathrm{m}^{2}\,\mathrm{s}^{-2}; W=kgm2s3\mathrm{W} = \mathrm{kg}\,\mathrm{m}^{2}\,\mathrm{s}^{-3}; Pa=kgm1s2\mathrm{Pa} = \mathrm{kg}\,\mathrm{m}^{-1}\,\mathrm{s}^{-2}; C=As\mathrm{C} = \mathrm{A}\,\mathrm{s}. Pam3=kgm1s2m3=kgm2s2=J\mathrm{Pa}\,\mathrm{m}^{3} = \mathrm{kg}\,\mathrm{m}^{-1}\,\mathrm{s}^{-2}\,\mathrm{m}^{3} = \mathrm{kg}\,\mathrm{m}^{2}\,\mathrm{s}^{-2} = \mathrm{J} (the work of pressure forces, P ⁣dVP\,\dd V).

Exercise 1.2

Give the dimensions of a pressure, a power, an electric charge, a voltage, and of the constants GG (in F=Gm1m2/r2F = Gm_1m_2/r^2), hh (in E=hνE = h\nu) and kBk_B (in E=kBTE = k_B T).

Solution

Solution of Exercise 1.2.

Pressure ML1T2M L^{-1} T^{-2}; power ML2T3M L^2 T^{-3}; charge ITI\,T; voltage == power/current =ML2T3I1= M L^2 T^{-3} I^{-1}. [G]=[Fr2/m2]=MLT2L2M2=L3M1T2[G] = [F r^2/m^2] = M L T^{-2} L^2 M^{-2} = L^3 M^{-1} T^{-2}; [h]=[E/ν]=ML2T2T=ML2T1[h] = [E/\nu] = M L^2 T^{-2} \cdot T = M L^2 T^{-1}; [kB]=[E/T]=ML2T2Θ1[k_B] = [E/T] = M L^2 T^{-2} \Theta^{-1}.

Exercise 1.3

Are these formulas homogeneous? v=2ghv = \sqrt{2gh}; E=12mvE = \tfrac12 m v; P=ρghP = \rho g h for a pressure; x=v0t+12gt2x = v_0 t + \tfrac12 g t^2; T=2πg/T = 2\pi\sqrt{g/\ell}; I=I0et/RCI = I_0 \eu^{-t/RC}.

Solution

Solution of Exercise 1.3.

2gh\sqrt{2gh}: LT2L=LT1\sqrt{L T^{-2} \cdot L} = L T^{-1}, homogeneous. 12mv\tfrac12 mv: MLT1ML2T2M L T^{-1} \neq M L^2 T^{-2}, not an energy. ρgh\rho g h: ML3LT2L=ML1T2M L^{-3} \cdot L T^{-2} \cdot L = M L^{-1} T^{-2}, a pressure. v0t+12gt2v_0 t + \tfrac12 g t^2: both terms LL, homogeneous. 2πg/2\pi\sqrt{g/\ell}: T2=T1\sqrt{T^{-2}} = T^{-1} — a frequency, not a period. I0et/RCI_0\eu^{-t/RC}: RCRC in ΩF=s\Omega\,\mathrm{F} = \mathrm{s}, argument dimensionless, homogeneous.

Exercise 1.4

A ruler is graduated in millimeters; a digital balance displays 0.01g0.01\,\mathrm{g} steps; a voltmeter’s data sheet promises “±0.5%\pm 0.5\% of the reading”. Give the type B standard uncertainty of a 12.7cm12.7\,\mathrm{cm} length, a 5.43g5.43\,\mathrm{g} mass and a 4.80V4.80\,\mathrm{V} reading.

Solution

Solution of Exercise 1.4.

Ruler: δ=1mm\delta = 1\,\mathrm{mm}, u=1/12=0.29mmu = 1/\sqrt{12} = 0.29\,\mathrm{mm}. Balance: u=0.01/12=0.003gu = 0.01/\sqrt{12} = 0.003\,\mathrm{g}. Voltmeter: a=0.005×4.80=0.024Va = 0.005 \times 4.80 = 0.024\,\mathrm{V}, u=0.024/3=0.014Vu = 0.024/\sqrt3 = 0.014\,\mathrm{V}.

Exercise 1.5 ★★

The speed vv of waves on a stretched string depends on its tension FF (a force) and its mass per unit length μ\mu. Find vv up to a dimensionless factor. Same question for the speed of sound in a gas of pressure PP and density ρ\rho.

Solution

Solution of Exercise 1.5.

v=kFαμβv = k F^\alpha \mu^\beta: LT1=(MLT2)α(ML1)βL T^{-1} = (M L T^{-2})^\alpha (M L^{-1})^\beta, so α+β=0\alpha + \beta = 0, αβ=1\alpha - \beta = 1, 2α=1-2\alpha = -1: α=12\alpha = \tfrac12, β=12\beta = -\tfrac12, v=kF/μv = k\sqrt{F/\mu} (in fact k=1k = 1). Sound: P/ρP/\rho has dimension ML1T2/(ML3)=L2T2M L^{-1} T^{-2} / (M L^{-3}) = L^2 T^{-2}, so v=kP/ρv = k\sqrt{P/\rho} (k=γ1.2k = \sqrt\gamma \approx 1.2 for air).

Exercise 1.6 ★★

Eight timings of the same fall, in milliseconds: 452452, 447447, 461461, 455455, 449449, 458458, 444444, 454454. Compute the mean, the experimental standard deviation and the standard uncertainty of the mean; write the result.

Solution

Solution of Exercise 1.6.

Sum =3620ms= 3620\,\mathrm{ms}, tˉ=452.5ms\bar t = 452.5\,\mathrm{ms}. Deviations 0.5,5.5,8.5,2.5,3.5,5.5,8.5,1.5-0.5, -5.5, 8.5, 2.5, -3.5, 5.5, -8.5, 1.5; squares sum to 226226; s=226/7=5.7mss = \sqrt{226/7} = 5.7\,\mathrm{ms}; u(tˉ)=5.7/8=2.0msu(\bar t) = 5.7/\sqrt8 = 2.0\,\mathrm{ms}. Result: t=(452.5±2.0)mst = (452.5 \pm 2.0)\,\mathrm{ms}.

Exercise 1.7 ★★

A resistance is measured as R=U/IR = U/I with U=4.80VU = 4.80\,\mathrm{V}, u(U)=0.03Vu(U) = 0.03\,\mathrm{V} and I=12.4mAI = 12.4\,\mathrm{mA}, u(I)=0.2mAu(I) = 0.2\,\mathrm{mA}. Compute RR and u(R)u(R); which measurement should be improved first?

Solution

Solution of Exercise 1.7.

R=4.80/0.0124=387ΩR = 4.80/0.0124 = 387\,\Omega. Relative uncertainties 0.03/4.80=0.63%0.03/4.80 = 0.63\% and 0.2/12.4=1.6%0.2/12.4 = 1.6\%; combined 0.632+1.62=1.7%\sqrt{0.63^2 + 1.6^2} = 1.7\%, u(R)=7Ωu(R) = 7\,\Omega: R=(387±7)ΩR = (387 \pm 7)\,\Omega. The current dominates — improve the ammeter (or its range) first.

Exercise 1.8 ★★

Two groups measure the speed of sound: (343±4)m/s(343 \pm 4)\,\mathrm{m}/\mathrm{s} and (336±3)m/s(336 \pm 3)\,\mathrm{m}/\mathrm{s}. Are they compatible? A third group finds (352±2)m/s(352 \pm 2)\,\mathrm{m}/\mathrm{s}: compatible with the first? with the second?

Solution

Solution of Exercise 1.8.

z=7/16+9=1.42z = 7/\sqrt{16 + 9} = 1.4 \leq 2: compatible. Third vs first: z=9/4+16=2.0z = 9/\sqrt{4 + 16} = 2.0: at the limit, barely compatible. Third vs second: z=16/4+9=4.4z = 16/\sqrt{4 + 9} = 4.4: incompatible — at least one of the two has a systematic error or an underestimated uncertainty.

Exercise 1.9 ★★

The volume of a sphere is computed from its diameter d=25.40mmd = 25.40\,\mathrm{mm} measured with a caliper, u(d)=0.02mmu(d) = 0.02\,\mathrm{mm}. Compute VV, its relative and absolute uncertainty, and write the result.

Solution

Solution of Exercise 1.9.

V=πd3/6=π×25.403/6=8580mm3V = \pi d^3/6 = \pi \times 25.40^3/6 = 8580\,\mathrm{mm}^{3}. u(V)/V=3u(d)/d=3×0.02/25.40=0.24%u(V)/V = 3\,u(d)/d = 3 \times 0.02/25.40 = 0.24\%, u(V)=20mm3u(V) = 20\,\mathrm{mm}^{3}: V=(8.58±0.02)cm3V = (8.58 \pm 0.02)\,\mathrm{cm}^{3}.

Exercise 1.10 ★★★

A small sphere of radius RR falls slowly in a viscous liquid; the drag depends on RR, on the speed vv and on the viscosity η\eta, of dimension ML1T1M L^{-1} T^{-1}. Show that F=kηRvF = k\,\eta R v. At high speed the drag depends instead on RR, vv and the fluid density ρ\rho: find the new law. Estimate, for a raindrop (R=1mmR = 1\,\mathrm{mm}, v=5m/sv = 5\,\mathrm{m}/\mathrm{s}, air: ρ=1.2kg/m3\rho = 1.2\,\mathrm{kg}/\mathrm{m}^{3}, η=1.8×105Pas\eta = 1.8 \times 10^{-5}\,\mathrm{Pa}\,\mathrm{s}), which regime applies by comparing the two estimates.

Solution

Solution of Exercise 1.10.

F=kηaRbvcF = k\eta^a R^b v^c: MM: a=1a = 1; TT: ac=2-a - c = -2, c=1c = 1; LL: a+b+c=1-a + b + c = 1, b=1b = 1: F=kηRvF = k\eta R v (Stokes, k=6πk = 6\pi). With ρ\rho instead: MM: a=1a = 1; TT: c=2-c = -2, c=2c = 2; LL: 3+b+2=1-3 + b + 2 = 1, b=2b = 2: F=kρR2v2F = k\rho R^2 v^2. Raindrop: ηRv=1.8×105×103×51×107N\eta R v = 1.8\times10^{-5} \times 10^{-3} \times 5 \approx 1 \times 10^{-7}\,\mathrm{N}; ρR2v2=1.2×106×25=3×105N\rho R^2 v^2 = 1.2 \times 10^{-6} \times 25 = 3 \times 10^{-5}\,\mathrm{N}, three hundred times larger: the inertial (high-speed) regime applies — the ratio ρvR/η330\rho v R/\eta \approx 330 is the Reynolds number.

Exercise 1.11 ★★★

Four measurements of the voltage across a resistor for increasing currents: (I,U)=(2.0mA,0.41V)(I, U) = (2.0\,\mathrm{mA}, 0.41\,\mathrm{V}), (4.0mA,0.78V)(4.0\,\mathrm{mA}, 0.78\,\mathrm{V}), (6.0mA,1.22V)(6.0\,\mathrm{mA}, 1.22\,\mathrm{V}), (8.0mA,1.59V)(8.0\,\mathrm{mA}, 1.59\,\mathrm{V}). Compute the least-squares slope and intercept; deduce RR; is the intercept compatible with zero if each UU carries u=0.02Vu = 0.02\,\mathrm{V}?

Solution

Solution of Exercise 1.11.

Iˉ=5.0mA\bar I = 5.0\,\mathrm{mA}, Uˉ=1.000V\bar U = 1.000\,\mathrm{V}; deviations in II: 3,1,1,3-3, -1, 1, 3; in UU: 0.59,0.22,0.22,0.59-0.59, -0.22, 0.22, 0.59. (IiIˉ)(UiUˉ)=3.98\sum (I_i - \bar I)(U_i - \bar U) = 3.98, (IiIˉ)2=20\sum (I_i - \bar I)^2 = 20: a=0.199V/mA=199Ωa = 0.199\,\mathrm{V}/\mathrm{mA} = 199\,\Omega, b=1.0000.199×5=0.005Vb = 1.000 - 0.199 \times 5 = 0.005\,\mathrm{V}. With u(U)=0.02Vu(U) = 0.02\,\mathrm{V} per point, an intercept of 0.005V0.005\,\mathrm{V} is well within uncertainty: compatible with zero, the resistor is ohmic, R199ΩR \approx 199\,\Omega.

Exercise 1.12 ★★★

From G=6.67×1011m3/(kgs2)G = 6.67 \times 10^{-11}\,\mathrm{m}^{3}/(\mathrm{kg}\,\mathrm{s}^{2}), h=6.63×1034Jsh = 6.63 \times 10^{-34}\,\mathrm{J}\,\mathrm{s} and c=3.00×108m/sc = 3.00 \times 10^{8}\,\mathrm{m}/\mathrm{s}, build by dimensional analysis a length, a time and a mass (the Planck units). Compute them and compare with the size of a proton (1×1015m1 \times 10^{-15}\,\mathrm{m}) and with a human mass.

Solution

Solution of Exercise 1.12.

=Gahbcc\ell = G^a h^b c^c: MM: a+b=0-a + b = 0; TT: 2abc=0-2a - b - c = 0; LL: 3a+2b+c=13a + 2b + c = 1. So b=ab = a, c=3ac = -3a, 2a=12a = 1: P=Gh/c3=4.42×1044/2.7×1025=4.1×1035m\ell_P = \sqrt{Gh/c^3} = \sqrt{4.42\times10^{-44}/2.7\times10^{25}} = 4.1 \times 10^{-35}\,\mathrm{m}; tP=P/c=1.4×1043st_P = \ell_P/c = 1.4 \times 10^{-43}\,\mathrm{s}; mP=hc/G=1.99×1025/6.67×1011=5.5×108kgm_P = \sqrt{hc/G} = \sqrt{1.99\times10^{-25}/6.67\times10^{-11}} = 5.5 \times 10^{-8}\,\mathrm{kg}. The Planck length is twenty orders of magnitude below a proton; the Planck mass, 55µg55\,\text{µ}\mathrm{g}, is a speck of dust — huge for a particle, tiny for a person (10910^9 times lighter).

The raw material of the chapter: a caliper, a stopwatch, a ball, a rule, and a column of repeated readings that never quite agree.
The raw material of the chapter: a caliper, a stopwatch, a ball, a rule, and a column of repeated readings that never quite agree.

1.7 Problem: Measuring gg with a pendulum

Problem 1.1

Weekend problem — a string, a bob, a stopwatch and a tape: how well can a kitchen table measure the gravitational field of the Earth, and which of its four numbers is the weak one

A pendulum is a steel ball hung from a fixed point by a light string; small swings have period τ=2π/g\tau = 2\pi\sqrt{\ell/g} where \ell is the distance from the suspension point to the ball’s center. The aim is a value of gg with a defensible uncertainty.

Part I — The model.

  1. Check the homogeneity of τ=2π/g\tau = 2\pi\sqrt{\ell/g}.
  2. The ball has mass mm. Argue by dimensions alone that the period cannot depend on mm if it depends only on \ell, gg and mm.
  3. The period also depends, slightly, on the amplitude θ0\theta_0 of the swing. Why does dimensional analysis allow this, and what must one do experimentally to make the formula applicable?
  4. Express gg as a function of \ell and τ\tau.

Part II — Measuring the length. The tape measure is graduated in millimeters; the string is 99.0cm99.0\,\mathrm{cm} from suspension to the top of the ball, the ball’s diameter is 2.00cm2.00\,\mathrm{cm}, both read once.

  1. Compute \ell.
  2. Evaluate the type B uncertainty of a single reading on the tape.
  3. Locating the suspension point and the ball’s center adds an estimated ±5mm\pm5\,\mathrm{mm} of possible error (uniform). Combine with the previous item into u()u(\ell), and give the relative uncertainty u()/u(\ell)/\ell.

Part III — Measuring the period. To reduce the effect of reaction time, one times ten periods. Eight such timings, in seconds: 20.1220.12, 20.0520.05, 20.1820.18, 20.0920.09, 20.1520.15, 20.0220.02, 20.1120.11, 20.0820.08.

  1. Why time ten periods rather than one? Why repeat the timing eight times?
  2. Compute the mean t10\overline{t_{10}} and the experimental standard deviation ss of the eight timings.
  3. Deduce the standard uncertainty of the mean, then the period τ\tau and its uncertainty u(τ)u(\tau).
  4. Compare u(τ)/τu(\tau)/\tau with u()/u(\ell)/\ell.

Part IV — The result and its weak point.

  1. Compute gg.
  2. Write u(g)/gu(g)/g in terms of u()/u(\ell)/\ell and u(τ)/τu(\tau)/\tau, and compute it.
  3. Write the result in the standard form.
  4. Compare with the reference value 9.81m/s29.81\,\mathrm{m}/\mathrm{s}^{2} through the normalized deviation. Verdict?
  5. Which of the two measured quantities limits the precision of gg? Propose two concrete improvements and say by how much each would reduce the relative uncertainty.
  6. Suppose that, by a careless alignment, every length was overestimated by 5mm5\,\mathrm{mm}. Is this accounted for in the uncertainty budget? What name does this kind of error carry, and how would one detect it?

Part V — Five lengths and a straight line. The experiment is repeated for =0.400\ell = 0.400, 0.6000.600, 0.8000.800, 1.0001.000 and 1.200m1.200\,\mathrm{m}, giving τ2=1.62\tau^2 = 1.62, 2.412.41, 3.253.25, 4.044.04 and 4.84s24.84\,\mathrm{s}^{2}, each with u(τ2)1%u(\tau^2) \approx 1\%.

  1. Why plot τ2\tau^2 against \ell rather than τ\tau against \ell?
  2. Compute ˉ\bar\ell and τ2\overline{\tau^2}.
  3. Compute the least-squares slope aa and intercept bb.
  4. Deduce gg from the slope.
  5. Interpret the intercept: what would a clearly non-zero bb reveal?
  6. The computer’s fit reports u(a)=0.05s2/mu(a) = 0.05\,\mathrm{s}^{2}/\mathrm{m}. Deduce u(g)u(g) and write the result.
  7. Compare the single-length result of Part IV with this one by the normalized deviation.
  8. A student proposes to measure instead one swing of a 10m10\,\mathrm{m} pendulum in a stairwell with the same tape and stopwatch. Estimate the relative uncertainty on gg that this would give, and conclude on the best strategy.
Solution

Solution of Problem 1.1.

1. L/(LT2)=T2=T\sqrt{L/(L T^{-2})} = \sqrt{T^2} = T.

2. τ=kαmβgγ\tau = k\ell^\alpha m^\beta g^\gamma gives T=Lα+γMβT2γT = L^{\alpha+\gamma} M^\beta T^{-2\gamma}, so β=0\beta = 0: no combination of \ell and gg can cancel a mass.

3. θ0\theta_0 is dimensionless, so any factor f(θ0)f(\theta_0) is invisible to dimensions. Keep the amplitude small (below about 1010^\circ), where f1f \approx 1 to better than 0.2%0.2\%.

4. g=4π2/τ2g = 4\pi^2\ell/\tau^2.

5. =99.0+1.00=100.0cm=1.000m\ell = 99.0 + 1.00 = 100.0\,\mathrm{cm} = 1.000\,\mathrm{m}.

6. u=1mm/12=0.29mmu = 1\,\mathrm{mm}/\sqrt{12} = 0.29\,\mathrm{mm}.

7. Alignment: 5/3=2.9mm5/\sqrt3 = 2.9\,\mathrm{mm}; combined 0.292+2.92=2.9mm\sqrt{0.29^2 + 2.9^2} = 2.9\,\mathrm{mm} — the tape reading is negligible. u()/=0.29%u(\ell)/\ell = 0.29\%.

8. The reaction-time error (about 0.1s0.1\,\mathrm{s}) is the same whether one or ten periods are timed, so its effect on one period is ten times smaller. Repeating reveals the scatter (hence ss) and divides the uncertainty of the mean by 8\sqrt8.

9. Sum =160.80s= 160.80\,\mathrm{s}, t10=20.100s\overline{t_{10}} = 20.100\,\mathrm{s}. Deviations 0.02,0.05,0.08,0.01,0.05,0.08,0.01,0.020.02, -0.05, 0.08, -0.01, 0.05, -0.08, 0.01, -0.02; squares sum to 188×104188 \times 10^{-4}; s=188×104/7=0.052ss = \sqrt{188\times10^{-4}/7} = 0.052\,\mathrm{s}.

10. u(t10)=0.052/8=0.018su(\overline{t_{10}}) = 0.052/\sqrt8 = 0.018\,\mathrm{s}; τ=2.0100s\tau = 2.0100\,\mathrm{s}, u(τ)=0.0018su(\tau) = 0.0018\,\mathrm{s}.

11. u(τ)/τ=0.09%u(\tau)/\tau = 0.09\%, three times smaller than u()/=0.29%u(\ell)/\ell = 0.29\%.

12. g=4π2×1.000/2.01002=39.48/4.040=9.772m/s2g = 4\pi^2 \times 1.000/2.0100^2 = 39.48/4.040 = 9.772\,\mathrm{m}/\mathrm{s}^{2}.

13. gτ2g \propto \ell\,\tau^{-2}: u(g)g=(u())2+(2u(τ)τ)2=0.292+0.182%=0.34%\dfrac{u(g)}{g} = \sqrt{\left(\dfrac{u(\ell)}{\ell}\right)^2 + \left(2\dfrac{u(\tau)}{\tau}\right)^2} = \sqrt{0.29^2 + 0.18^2}\,\% = 0.34\%.

14. u(g)=0.0034×9.772=0.033m/s2u(g) = 0.0034 \times 9.772 = 0.033\,\mathrm{m}/\mathrm{s}^{2}: g=(9.77±0.03)m/s2g = (9.77 \pm 0.03)\,\mathrm{m}/\mathrm{s}^{2}.

15. z=9.7729.81/0.033=1.22z = \abs{9.772 - 9.81}/0.033 = 1.2 \leq 2: compatible.

16. The length (0.29%0.29\% against 0.18%0.18\% for the period term). Locating both ends to ±1mm\pm1\,\mathrm{mm} brings u()/u(\ell)/\ell to 0.06%0.06\% and u(g)/gu(g)/g to 0.062+0.182=0.19%\sqrt{0.06^2 + 0.18^2} = 0.19\%; doubling the length halves u()/u(\ell)/\ell and lengthens τ\tau, helping both terms. Timing twenty periods instead of ten halves the period term only: u(g)/g=0.292+0.092=0.30%u(g)/g = \sqrt{0.29^2 + 0.09^2} = 0.30\% — a small gain while the length is the bottleneck.

17. No: a shift common to every reading is a systematic error (a bias); the budget covers random scatter and declared tolerances, not a consistent offset. It is detected by changing method or instrument, or — Part V — by a non-zero intercept of the τ2\tau^2\ell line.

18. τ2=(4π2/g)\tau^2 = (4\pi^2/g)\,\ell is a straight line through the origin: linearity and the zero intercept are both testable by eye and by least squares; τ\tau \propto \sqrt\ell is not.

19. ˉ=0.800m\bar\ell = 0.800\,\mathrm{m}; τ2=16.16/5=3.232s2\overline{\tau^2} = 16.16/5 = 3.232\,\mathrm{s}^{2}.

20. Deviations in \ell: 0.4,0.2,0,0.2,0.4-0.4, -0.2, 0, 0.2, 0.4; in τ2\tau^2: 1.612,0.822,0.018,0.808,1.608-1.612, -0.822, 0.018, 0.808, 1.608. Cross sum =1.614= 1.614, (iˉ)2=0.400\sum(\ell_i - \bar\ell)^2 = 0.400: a=4.035s2/ma = 4.035\,\mathrm{s}^{2}/\mathrm{m}, b=3.2324.035×0.800=0.004s2b = 3.232 - 4.035 \times 0.800 = 0.004\,\mathrm{s}^{2}.

21. g=4π2/a=39.48/4.035=9.78m/s2g = 4\pi^2/a = 39.48/4.035 = 9.78\,\mathrm{m}/\mathrm{s}^{2}.

22. b0b \approx 0, as the model predicts. A clearly non-zero bb would betray a systematic offset δ\delta in every length (τ2=a(measδ)\tau^2 = a(\ell_{\text{meas}} - \delta) gives b=aδb = -a\delta), e.g. a mislocated ball center — or an amplitude effect.

23. u(g)/g=u(a)/a=0.05/4.035=1.2%u(g)/g = u(a)/a = 0.05/4.035 = 1.2\%, u(g)=0.12m/s2u(g) = 0.12\,\mathrm{m}/\mathrm{s}^{2}: g=(9.78±0.12)m/s2g = (9.78 \pm 0.12)\,\mathrm{m}/\mathrm{s}^{2}.

24. z=9.7729.784/0.0332+0.122=0.1z = \abs{9.772 - 9.784}/\sqrt{0.033^2 + 0.12^2} = 0.1: fully compatible. The single careful measurement is more precise; the line validates the model and guards against a length offset.

25. τ=2π10/9.8=6.3s\tau = 2\pi\sqrt{10/9.8} = 6.3\,\mathrm{s}; one swing timed with a reaction time of about 0.2s0.2\,\mathrm{s} gives u(τ)/τ3%u(\tau)/\tau \approx 3\%, hence u(g)/g6%u(g)/g \approx 6\%, while u()/u(\ell)/\ell drops to 0.03%0.03\%: twenty times worse than Part IV. The length gain is wasted by timing a single period. Best strategy: long pendulum and many periods, repeated — ten swings of the 10m10\,\mathrm{m} pendulum, eight times, gives u(τ)/τ0.03%u(\tau)/\tau \approx 0.03\% and u(g)/g0.07%u(g)/g \approx 0.07\%.

Terms defined in this chapter

See all 393 terms in the glossary