Dimensional analysis is a powerful tool based on a simple idea: look at the
kinds of quantity involved in a system, as reflected by their units of
measurement, and make sure that your equations and analysis are consisistent
about this simple information about the quantities. The kind

of a
quantity is also called its dimension

; but this leads to confusion with
dimension in the sense involved in linear spaces, so I'll call it a kind here.

There are two primary uses for dimensional (kindly ?) analysis: sanity-checking some formula one has derived (are the quantities to be added together, or declared equal to one another, of the same kind ?) and working out what formulae could be relevant – which generally comes down to determining the ways one can combine the known data of the system to yield a dimensionless quantity – that is, one of the same kind as pure numbers; a ratio of quantities of like kind. Another variant on the second approach is to identify combinations of data of the system that are of some simple kinds, that let us reduce our analysis to these combinations rather than dealing independently with all the data they combine. For example, the system may offer us a natural charge-to-mass ratio that we can use to eliminate the difference between charges and masses in the system – or, even to eliminate charges and masses altogether (see the second example below).

It is usual to express the kinds of quantities in terms of four primitive kinds:

- M, standing for mass
- L, standing for length
- T, standing for time, and
- Q, standing for electrical charge.

I'll write ~ for is of kind

, so any mass ~ M, speed ~ L/T, and so on.
Thus energy ~ L.L.M/T/T; electrostatic potential (voltage) is energy per charge
so ~ L.L.M/Q/T/T; and electrical impedence is potential over current so ~
L.L.M/T/Q/Q.

To illustrate the approach, let's consider the case of a simple pendulum. This comprises a massive body attached to one end of a rod (typically light enough that we can ignore its mass relative to the main one) whose other end is fixed but free to turn, with negligible friction, so that the rod and mass can swing back and forth. It is observed that such a pendulum swings to and fro, repeating its earlier motion cyclically with a definite period which doesn't seem to depend much on the size of the movement, as long as this is tiny compared to the length of the rod. The known data of our system are:

- m = mass
- How heavy is the object on the end of the rod ? Of kind M.
- h = length
- How far is this mass from the pivot point ? Of kind L.
- r = range
- How far apart are the extremes of the object's movement ? Of kind L.
- t = period
- How long does the pendulum take to complete one cycle of its movement ? Of kind T.
- g = gravitational field strength
- How rapidly do unsupported massive objects accellerate downwards ? Of kind L/T/T.

Now, we're given that t doesn't seem to depend much on r, at least as long as r/h is small, so let's ignore r for now. Only m has M in its kind, so we can immediately see that it can't play any part in the relationship among the other quantities: we can predict that the period of a pendulum doesn't depend on its mass (as long as this is large enough to dwarf the mass of its supporting rod, and large enough that its weight dwarfs the forces we're ignoring, such as air resistance). We can make this physical prediction based only on looking at the kinds of quantity involved !

Furthermore, when we look at h, t and g we can see that the only dimensionless quantity we can form from them is g.t.t/h. We can expect any such dimensionless quantity to be a constant determined by the abstract geometry (without reference to absolute scale) of our system; it may have stray factors of π and other pure numbers in it, but our system is simple enough that we can expect it to have little more. We can thus predict that the period, t, should be proportional to √(h/g).

Now, of course, such predictions should always be tested experimentally:
there may be some other quantities involved which we haven't noticed, which
might be of kinds that mix up the quantities we were taking into account,
thereby complicating the story. There might be some universal constant (or
parameter of our context, cf. gravitational field strength) which we didn't take
into account, that isn't dimensionless. However, it commonly happens that one
does identify the correct quantities and, from a knowledge of which forces are
involved, which universal constants are in play. From these, dimensional
analysis yields the most obvious

hypotheses to test; and these often are
correct, at least to a good approximation. This can save a lot of effort on
the way to discovering a full and accurate description of the system.

While I'll leave the main discussion of this elsewhere, Planck's units are an example of dimensional
analysis. In 1904, Max Planck noticed that certain theoretical complications in
the relationships between electrodynamics and thermodynamics could be resolved
by assuming that light's interaction with matter only happens in discrete
chunks, transferring an amount of energy proportional to the frequency of the
light involved. The constant of proportionality, h, is known as Planck's
constant; and h ~ energy / frequency = L.L.M/T. This only involves length mass
and time, but Planck noticed that it couldn't be expressed in terms of the kinds
of the other known universal constants: the speed of light, c ~ L/T, and
Newton's constant, G ~ L.L.L/M/T/T. (The electromagnetic constants
μ_{0}, ε_{0} and Z_{0} can all be determined,
given c, from any one; and each involves charge; so could not be relevant.)
This meant that c, h and G are three independent kinds expressible purely in
terms of the three base kinds M, L and T; from which he could infer that it must
be possible (because there are three of each) to invert the system and express
M, L and T in terms of c, h and G. The interested reader is encouraged to work
out how to do that, before reading on.

We can use c to eliminate either L or T from the others: G/c/c ~ L/M and h/c
~ L.M. Multiplying these together we get L.L while dividing the second by the
first yields M.M; the square roots of these are then √(h.G/c)/c ~ L and
√(c.h/G) ~ M. Applying c to the first we get √(h.G/c)/c/c ~ T.
When we work out the values of these, we find that the length and time are
spectacularly tiny (even when compared to, say, the diameter of an atomic
nucleus, although Planck didn't know about those at the time); while the mass is
merely small – a droplet of water with this mass is about a third of a
millimetre across and may have fairly complex life-forms doing interesting
things inside it. It is common to use these units in analysing problems in
theoretical physics, so as to make c, h and G all take the value 1 and thus save
the need to keep track of them: I, however, prefer to keep track of our arguably
spurious units, since doing so leaves me able to use dimensional analysis to
check that I haven't made stupid mistakes (or, often enough, makes it easier to
*find* my stupid mistake when I *have* made one).

When it comes to looking at a smooth manifold, we can distinguish displacements within the manifold from lengths and times: the metric (and the speed of light) perform conversions between these, but displacements themselves are not lengths or times. So add a kind S to describe displacements. At the same time, relativity definitely tells us to treat length and time as equivalent; so I'll suppress T in favour of L but keep track of the factors of c (the speed of light) that this implies.

The metric of space-time, g (not to be confused with the simple pendulum's
g), is combined multiplicatively with two displacements to yield the square of a
length (or time, but I'll stick with length). So g ~ L.L/S/S. Now, g is a
linear map from tangents, intrinsically of kind S, to gradients, intrinsically
of kind 1/S; but it's throwing in a dimensional factor of L.L in the process.
In practice, we can't measure displacements, *per se*, only lengths and the
like: so we can't directly examine any tensor rank unless its intrinsic factors
of S all cancel. Thus we can look at linear maps from tangents to tangents
(because their tensor rank is a tangent as output times a gradient that consumes
the input, so the implicit factors of S cancel) but we can't look at tangents
themselves. Equally, we can't measure g; but we can use it to convert tensor
ranks to ones we can measure; doing so amounts to converting every factor of S
to a factor of L.

The differential operator on space-time (with respect to which g is
constant) maps any function to one which, when multiplied by a displacement,
yields a change in the original function; so the differential operator's kind is
1/S. The Riemann tensor is a linear map from gradient fields to tensor fields
of third rank, all factors in which are gradient fields; it determines the
result of antisymmetrising the second derivative of a gradient field in the two
factors of gradient field brought into play by the two uses of the differential
operator. Thus the Riemann tensor's kind is just that of the square of the
differential operator, 1/S/S. The Ricci tensor is a trace of the Riemann tensor
which contracts out the tensor factors associated with the input

gradient
field, so is consequently of the same kind.

We have a gravitational field equation, linking the Ricci tensor to the energy-momentum-stress tensor, an electromagnetic source equation linking the relativistic electromagnetic field (which unifies the electric and magnetic fields) to the 4-current and various equations of consistency. These last all equate some single expression to zero – tacitly the zero of whatever kind the relevant thing is – so don't give us any scope for dimensional analysis (although they would if they equated some sum of expressions to zero). The Bianchi identities, constraining the Ricci tensor, are of this kind; as are the conservation-law for current, d^(μ(J)) = zero, and the consistency equation of the electromagnetic field, d^F = zero (see below for what μ, J and F are).

The electromagnetic field tensor, F, encodes the electric and magnetic fields; it's of the same kind as E or B.c, so it's a force per unit charge. However, it's also a tensor field of rank gradient times gradient, so the physical-ranked thing derived from it that really has this kind is F/g; so F/g ~ M.L/Q/T/T = c.c.M/Q/L and F ~ c.c.M.L/Q/S/S. The 4-current density, J, is a tangent quantity so needs a factor of S/L times its intrinsic current per unit area or speed times charge per unit volume, c.Q/L/L/L, yielding J ~ c.Q.S/L/L/L/L. It's linked to F by a field equation

- d^(μ(g\F/g)) = μ(J).Z
_{0}

wherein Z_{0} is the impedance of free space (a universal
constant), g\F/g = Alt(2, g, F) is a tangent times tangent quantity obtained
from F by applying g's inverse to each (gradient) tensor rank of F and
μ×μ = −det(g) defines (up to a constant sign) μ as an
alternating form, of kind Alt(dim, {gradients}) when the manifold has dimension
dim, which can be used as a linear map which contracts away tangent ranks of a
tensor (antisymmetrically); when it contracts away n of them, it contributes
dim−n remaining gradient ranks to the result. Thus μ's rank is L/S
raised to power dim; but it appears on both sides of the equation, so we can
ignore it in our dimensional analysis. We can now sanity-check our equation:

- g\F/g
- ~ (L.L/S/S) \ c.c.M/Q/L
- = c.c.S.S.M/Q/L/L/L;
- J.Z
_{0} - ~ (c.Q.S/L/L/L/L).(L.L.M/T/Q/Q)
- = c.S.M/Q/T/L/L
- = c.c.S.M/Q/L/L/L,

agreeing with what's left of the equation's other side after d^ (the alternating differential operator) has knocked out one factor of S.

The remaining contribution of F to the (dimensionally non-fatuous) field
equations is via its contribution to the energy-momentum-stress tensor, which is
(F/g\F −g.trace(F/g\F/g)/dim)/c/Z_{0} in which the second term
trivially has the same dimensions as the first (it divides by a factor of g and
multiplies by a factor of g) so I'll ignore it. Fitting together the parts we
have yields F/g\F/c/Z_{0} ~ c.c.M/L/S/S, in which c.c.M is indeed an
energy; the tensor is of rank gradient times gradient, the same as g, so
dividing it by g would turn the /S/S into /L/L and give us an energy per unit
volume, or energy density, as we should expect.

The field equation for gravity multiplies the energy-momentum-stress tensor
by κ = 8.π.G/c/c/c/c and identifies it with Ricci
−g.(trace(Ricci/g)/2 +Λ); as before, the g.trace term is
transparently of the same kind as the first, 1/S/S, so I'll ignore it, aside
from noting that the cosmological constant, Λ, has the same kind as
Ricci/g, so Λ ~ 1/L/L is an inverse area. Now, G is Newton's constant,
so ~ force.L.L/M/M = c.c.L/M, so κ ~ L/M/c/c and
κ.F/g\F/c/Z_{0} ~ 1/S/S matches Ricci.

In the last, we scaled F/g\F, which has two factors of F in it, by
κ/c/Z_{0} = 8.π.G/Z_{0}/c/c/c/c/c ~ Q.Q/M/M/c/c/c/c,
which is the square of a Q/M/c/c or charge-to-energy ratio. If we scale F ~
c.c.M.L/Q/S/S by a Q/M/c/c, the result would be a simple L/S/S, whose physical
expression is obtained by dividing by g to obtain a 1/L. The natural
charge-to-energy ratio derived from the universal constants of the field
equations is ρ = √(4.π.G/Z_{0}/c)/c/c – in which the
factor of 4.π comes from G's rôle in Newton's equation of gravitation
being equivalent to that of 1/(4.π.ε_{0}) =
c.Z_{0}/4/π in Coulomb's equivalent force law for electrostatics
– with a value of about 0.9587e-27 Coulomb / Joule. So let f = ρ.F ~
L/S/S and the equation of gravitation becomes

- Ricci −g.(trace(Ricci/g)/2 −Λ) = f/g\f −g.trace(f/g\f/g)/dim +κ.W

wherein W is the energy-momentum-stress contribution from everything but the electromagnetic field (which is presumably zero in the free-field case). We can re-arrange this as

- Ricci −f/g\f −g.(trace((Ricci/2 −f/g\f/dim)/g) +Λ) = κ.W

in which the left side is expressed in terms of fields which only involve length and displacement as kinds. Naturally, we must now replace F with f throughout our other field equations; d^F = zero simply becomes d^f = zero, which is easy enough. The interesting one is then

- d^(μ(g\f/g)) = μ(J).ρ.Z
_{0}

naturally encouraging us to introduce j = J.ρ.Z_{0} ~
S/L/L/L, a displacement density, whose physical manifestation would be an
inverse area. Its conservation law would simply be d^(μ(j)) = zero as for J.
The Lorentz force law would then give the (proper time) rate of transfer, per
unit volume, of 4-momentum between a charged medium carrying current
density

j and the electromagnetic field as

4-force density

- = g\F·J/c
in which the g\ is present to turn the covector (gradient-valued) F·J into a vector (tangent-valued) quantity and the /c is present to turn the current density J into a charge density to be acted on by the force per unit charge, F.

- = g\f·j/ρ/ρ/Z
_{0}/c

The constant ρ.ρ.Z_{0}.c = 4.π.G/c/c/c/c ~ L/M/c/c,
so its inverse (by which we're multiplying g\f·j) has dimensions of
energy over length, i.e. force; while g\f·j ~ S/L/L/L/L, an S/L per unit
volume, so our 4-force density's physical manifestation would indeed have
dimensions of force per unit volume. Note that the force constant, a kind of
Planck force, is the only physical constant

remaining in the free-field
theory (i.e. when W is zero); and its value is

- c.c.c.c/4/π/G ≈ 9.634e42 Newtons

which, when divided by Earth's surface gravity to get the mass whose weight it would equal, comes out at about 1.6 times the mass of our Milky Way galaxy.

Written by Eddy.