Skip to article frontmatterSkip to article content

2Newton’s Laws

Now we turn to one of the most profound breakthroughs in the history of science: the laws of motion formulated by Isaac Newton. These laws provide a systematic framework for understanding how and why objects move, and form the backbone of classical mechanics. Using these three laws we can predict the motion of a falling apple, a car accelerating down the road, or a satellite orbiting Earth (though some adjustments are required in this context to make use of e.g. GPS!). More than just equations, they express deep principles about the nature of force, mass, and interaction.

In this chapter, you will begin to develop the core physicist’s skill: building a simplified model of the real world, applying physical principles, and using mathematical tools to reach meaningful conclusions.

2.1Newton’s Three Laws

Much of physics, in particular Classical Mechanics, rests on three laws that carry Newton’s name:

N1 has, in fact, been formulated by Galileo Galilei. Newton has, in his N2, build upon it: N1 is included in N2, after all:
if F=0\vec{F} = 0, then dpdt=0p=constantv=constant\frac{d\vec{p}}{dt} = 0 \rightarrow \vec{p} = \text{constant} \rightarrow \vec{v} = \text{constant}, provided mm is a constant.

Most people know N2 as

F=ma\vec{F} = m \vec{a}

For particles of constant mass, the two are equivalent:
if m=constant m = \text{constant} , then

dpdt=mdvdt=ma\frac{d\vec{p}}{dt} = m\frac{d\vec{v}}{dt} = m\vec{a}

Nevertheless, in many cases using the momentum representation is beneficial. The reason is that momentum is one of the key quantities in physics. This is due to the underlying conservation law, that we will derive in a minute. Moreover, using momentum allows for a new interpretation of force: force is that quantity that - provided it is allowed to act for some time interval on an object - changes the momentum of that object. This can be formally written as:

dp=FdtΔp=Fdtd\vec{p} = \vec{F} dt \leftrightarrow \Delta \vec{p} = \int \vec{F} dt

The latter quantity IFdt\vec{I} \equiv \int \vec{F} dt is called the impulse.

In Newton’s Laws, velocity, acceleration and momentum are key quantities. We repeat here their formal definition.

Solution to Exercise 1
Interpret
Develop

First we make a sketch.

This is obviously a 1-dimensional problem. So, we can leave out the vector character of e.g. the force.

Solution to Exercise 3

The gravitational force acts from the earth on the jumper. Newton’s law states that the jumper thus acts a gravitational force on the earth. Hence, the earth accelerates towards the jumper!

Although this sounds silly, when comparing this idea to the sun and the planets, we must draw the conclusion that the sun is actually wobbling as it is pulled towards the various planets! See also this animated explanation

2.2Newton’s laws applied

2.2.1Force addition, subtraction and decomposition

Newton’s laws describe how forces affect motion, and applying them often requires combining multiple forces acting on an object, see Figure 3. This is done through vector addition, subtraction, and decomposition—allowing us to find the net force and analyze its components in different directions, see this chapter in the book on linear algebra for a full elaboration on vector addition and subtraction.

Three forces acting on a particle. In which direction will it accelerate?

Figure 3:Three forces acting on a particle. In which direction will it accelerate?

Solution to Exercise 4
Fnet=Fi=F1+F2+F3=(104)+(113)+(10.51)=(1+110+1+0.54+3+1)=(10.50)\begin{aligned} \vec{F}_{net} &= \sum\vec{F}_i=\vec{F}_1+\vec{F}_2+\vec{F}_3\\ &= \begin{pmatrix} 1 \\ 0 \\ -4 \end{pmatrix} + \begin{pmatrix} 1 \\ 1 \\3 \end{pmatrix} + \begin{pmatrix} -1 \\ -0.5 \\ 1 \end{pmatrix} &= \begin{pmatrix} 1 + 1 -1 \\ 0 + 1 + -0.5 \\ -4 + 3 + 1 \end{pmatrix} = \begin{pmatrix} 1 \\ 0.5 \\ 0 \end{pmatrix} \end{aligned}

Hence, the net force acting on the particle is 12+.52=1.1N\sqrt{1^2+.5^2}=1.1\text{N} and the particle will accelerate in the direction (10.5)\begin{pmatrix} 1 \\ 0.5 \end{pmatrix}, just like in the previous example. The magnitude of the acceleration is a=F/ma=F/m and can only be calculated when the mass of the particle is specified.


2.2.2Acceleration due to gravity

In most cases the forces acting on an object are not constant. However, there is a classical case that is treated in physics (already at secondary school level) where only one, constant force acts and other forces are neglected. Hence, according to Newton’s second law, the acceleration is constant.

When we first consider only the motion in the z-direction, we can derive:

a=Fm=const.a=\frac{F}{m}=\text{const.}

Hence, for the velocity:

v(t)=t0teadt=a(tet0)+v0v(t) = \int_{t_0}^{t_e} a dt = a(t_e-t_0) + v_0

assuming t0=0 and te=tv(t)=v0+att_0=0 \text{ and } t_e=t \Rightarrow v(t)=v_0+at the position is described by

s(t)=0tv(t)dt=0tat+v0dt=12at2+v0t+s0s(t) = \int_{0}^{t} v(t)dt = \int_{0}^{t} at + v_0 dt = \frac{1}{2}at^2 + v_0t + s_0

Rearranging:

s(t)=12at2+v0t+s0s(t) = \frac{1}{2}at^2 + v_0t + s_0
Solution to Exercise 5
Interpret
Develop
Evaluate
Assess
A free body diagram of the situation with all relevant quantities.

A free body diagram of the situation with all relevant quantities.

Only gravity acts on the stone (in the downward direction). We will call the position of the stone at time tt: s(t)s(t)
Initial conditions: t=0s(0)=s0=1.5m and s˙=v=v0=10m/st=0 \rightarrow s(0) = s_0 = 1.5\mathrm{ m} \text{ and } \dot{s} = v = v_0 = 10 \mathrm{ m/s}

NOTE: Some of these solutions can be derived more easily using the concept of conservation of energy which will be covered in one of the next chapters.

Solution to Exercise 7

The horizontal traveled distance is the same per time unit. For the vertical traveled distance it decreases until vy=0v_y=0 and then increases.

Interpret
Develop
Evaluate
Solution to Exercise 8

The acceleration of gravity is found by setting the gravitation force equal to mg-mg:

GmMr2=mg(r)g(r)=GMr2-G\frac{mM}{r^2} = -mg(r) \Rightarrow g(r) = G\frac{M}{r^2}

with MM the mass of the earth.

At the surface of the earth, r=Rer = R_e we have for the value of ge=9.81m/s2g_e = 9.81\mathrm{m/s}^2. We look for the height above the earth surface where gg has dropped to 9.80m/s29.80 \mathrm{m/s}^2. If we call this height HH, we write for the distance to the center of the earth r=Re+Hr = R_e +H.

Thus, we look for g(r)ge=9.80(m/s2)9.81(m/s2)=0.999\frac{g(r)}{g_e} = \frac{9.80 (\mathrm{m/s}^2)}{9.81 (\mathrm{m/s}^2)} = 0.999:

g(r)ge=GM/r2GM/Re2Re2r2=Re2(Re+H)2=9.809.81=0.999\frac{g(r)}{g_e} = \frac{GM/r^2}{GM/R_e^2} \rightarrow \frac{R_e^2}{r^2} = \frac{R_e^2}{(R_e+H)^2} = \frac{9.80}{9.81} = 0.999

If we solve HH from this equation we find: H=3.25kmH = 3.25 \mathrm{km} (we used Re=6378kmR_e = 6378 \mathrm{km}).

If we would have said: ‘significant change’ in means g9.819.71m/s2g \rightarrow 9.81 \to 9.71\mathrm{m/s}^2, we would have found H=32.8kmH = 32.8 \mathrm{km}.

Source
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.animation import FuncAnimation
from IPython.display import HTML
plt.rcParams['font.family'] = 'Segoe UI Emoji'

# Simulationparameters
dt = 0.05
t_max = 10
t_values = np.arange(0, t_max, dt)

# Physical parameters
vx = 1.0
Fy = 1.0
m = 1.0
ay = Fy / m

# Calculation trajectory
x = vx * t_values
y = np.zeros_like(t_values)

x_burn_start = 2.0
x_burn_end = 4.0
i_start = np.argmax(x >= x_burn_start)
i_end = np.argmax(x >= x_burn_end)

for i in range(i_start, i_end+1):
    t_burn = t_values[i] - t_values[i_start]
    y[i] = 0.5 * ay * t_burn**2

vy_final = ay * (t_values[i_end] - t_values[i_start])
y0 = y[i_end]
t0 = t_values[i_end]

for i in range(i_end, len(t_values)):
    y[i] = y0 + vy_final * (t_values[i] - t0)

# Plot
fig, ax = plt.subplots(figsize=(8, 4))
ax.set_xlim(0, np.max(x)+1)
ax.set_ylim(0, np.max(y)+1)
ax.set_xlabel("x")
ax.set_ylabel("y")
ax.set_title("Rocket with propulsion on between x=2 en x=4")

# Rocket 
rocket = ax.text(0, 0, "\U0001F680", fontsize=14)

# Trail
trail, = ax.plot([], [], 'r-', lw=1)

# Time
time_text = ax.text(0.98, 0.95, '', transform=ax.transAxes,
                    ha='right', va='top', fontsize=12)

# Init
def init():
    rocket.set_position((0, 0))
    trail.set_data([], [])
    time_text.set_text('')
    return rocket, trail, time_text

# Update
def update(frame):
    rocket.set_position((x[frame], y[frame]))
    trail.set_data(x[:frame+1], y[:frame+1])
    time_text.set_text(f"t = {t_values[frame]:.2f} s")
    return rocket, trail, time_text

# Animation
ani = FuncAnimation(fig, update, frames=len(t_values),
                    init_func=init, interval=dt*1000, blit=True)

plt.close()
HTML(ani.to_jshtml())
Loading...
Solution to Exercise 10

Two are easy: uniform motion (F=ma=0s=v0tF=ma=0 \rightarrow s=v_0t) and constant acceleration a=consts=1/2at2a=const \rightarrow s=1/2at^2, with a=2v02sa=\frac{_2v_0^2}{s}.

Consider the third being a harmonic oscillating force field: F(t)=Asin(2πft)F(t)=A\sin(2\pi ft) Then the equation of motion becomes:

a=F/m=Amsin(2πft)a = F/m = \frac{A}{m}\sin(2\pi ft)
v=adt=Am2πfcos(2πft)+C0v = \int{a}dt = \frac{A}{m2\pi f}\cos(2\pi ft) + C_0

Assuming v(0)=0C0=Am2πfv(0)=0 \rightarrow C_0 = -\frac{A}{m2\pi f}

And,

x=vdt=Am(2πf)2sin(2πft)+C0t+C1x = \int{v}dt = \frac{A}{m(2\pi f)^2}\sin(2\pi ft) + C_0t + C_1

Assuming x(0)=0C1=0x(0)=0 \rightarrow C_1 = 0

Hence:

x=Am(2πf)2sin(2πft)Am2πftx = \frac{A}{m(2\pi f)^2}\sin(2\pi ft) - \frac{A}{m2\pi f}t

Now, finding traveling the same distance in the same time AND the harmonic oscillation is complete (hence, f=1tef=\frac{1}{t_e}):

v0te=Ate2m(2π)2sin(2π)Atem2πtev_0t_e = \frac{A t_e^2}{m(2\pi)^2}\sin(2\pi) - \frac{A t_e}{m2\pi}t_e
v0te=Ate2m2πv_0t_e = - \frac{A t_e^2}{m2\pi}
v0=Atem2πv_0 = - \frac{A t_e}{m2\pi}
mA=tev02π\frac{m}{A} = - \frac{t_e}{v_02\pi}
Source
# Animatie van een deeltje met constante snelheid en een deeltje met constante versnelling
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.animation import FuncAnimation
from IPython.display import HTML

# Parameters
v = 1.0                         # snelheid in x-richting
dt = 0.05                       # tijdstap in seconden
t_max = 10 + dt                 # totale simulatie tijd
y = 0.5                         # constante y-positie
a = 2*v**2/(v*(t_max-dt))       # versnelling in x-richting
m = 1.0                         # massa van het deeltje
A = - v*2*np.pi / (t_max-dt)    # amplitude van de sinusgolf
f = 1 / (t_max-dt)              # frequentie van de sinusgolf

# Vooraf posities berekenen
t_values = np.arange(0, t_max, dt)
x_values = v * t_values
x_values_2 = 1/2 * a * t_values**2  # voor een andere beweging
x_values_3 = A / (m * (2 * np.pi * f)**2) * np.sin(2 * np.pi * f * t_values) - A / (m * 2 * np.pi * f) * t_values

y_values = np.full_like(x_values, y)  # constante y

# Setup plot
fig, ax = plt.subplots(figsize=(8, 4))
ax.set_xlim(0, x_values[-1] + 1)
ax.set_ylim(0, 1)
ax.set_xlabel("x")
ax.set_ylabel("y")


particle, = ax.plot([], [], 'ro', markersize=10, label='uniform motion')
particle_2, = ax.plot([], [], 'bo', markersize=10, label='uniform acceleration')
particle_3, = ax.plot([], [], 'go', markersize=10, label='harmonic oscillating force')

ax.legend(loc='lower right')
time_text = ax.text(0.98, 0.95, '', transform=ax.transAxes,
                    ha='right', va='top', fontsize=12)

# Initialisatie
def init():
    particle.set_data([], [])
    particle_2.set_data([], [])
    particle_3.set_data([], [])
    time_text.set_text('')
    return particle, time_text

# Update per frame
def update(frame):
    x = x_values[frame]
    x_2 = x_values_2[frame]
    x_3 = x_values_3[frame]
    y = y_values[frame]
    t = t_values[frame]
    particle.set_data([x],[y])
    particle_2.set_data([x_2], [.5*3*y])
    particle_3.set_data([x_3], [.5*y])
    time_text.set_text(f"t = {t:.2f} s")
    return particle, time_text

# Animatie
ani = FuncAnimation(fig, update, frames=len(t_values),
                    init_func=init, interval=dt*1000, blit=True)

plt.close()
HTML(ani.to_jshtml())
Loading...

2.2.3Frictional forces

There are two main types of frictional force:

  • Static friction prevents an object from starting to move. It adjusts in magnitude up to a maximum value, depending on how much force is trying to move the object. This maximum is given by
Fstatic,max=μsNF_{static,max}=\mu_sN

where μs\mu_s is the coefficient of static friction and NN is the normal force. If the applied force exceeds this maximum, the object begins to slide.

  • Kinetic (dynamic) friction opposes motion once the object is sliding. Its magnitude is generally constant and given by
Fkinetic=μkNF_{kinetic} = \mu_k N

where μk\mu_k is the coefficient of kinetic friction. This force does not depend on the velocity of the object, only on the normal force and surface characteristics.

Friction always acts opposite to the direction of intended or actual motion and is essential in both preventing and controlling movement.

Material PairStatic Friction (μs\mu_s)Kinetic Friction (μk\mu_k)
Rubber on dry concrete1.00.8
Steel on steel (dry)0.740.57
Wood on wood (dry)0.50.3
Aluminum on steel0.610.47
Ice on ice0.10.03
Glass on glass0.940.4
Copper on steel0.530.36
Teflon on Teflon0.040.04
Rubber on wet concrete0.60.5
Leather on wood0.560.4

Values are approximate and can vary depending on surface conditions.

Solution to Exercise 11
  1. There a two forces acting on mm parallel to the inclined plane: friction and gravity’s component parallel to the slope. These two determine the motion along the slope: if we tilt the plane the component of gravity parallel to the slope gets bigger. The particle will start moving when we pass: Fgx=Fsmgsin(θ)=mgμscos(θ)θmax=tan1(μs)F_{g_x} = F_s \rightarrow mg\sin(\theta) = mg\mu_s\cos(\theta) \Rightarrow \theta_{max} = \tan^{-1}(\mu_s)
  2. Once the particle is sliding downward, gravity and the kinetic friction determine how fast:
    Fnet=FgxFfma=mgsin(θ)mgμkcos(θ)F_{net} = F_{g_x} - F_f \rightarrow ma = mg\sin(\theta) - mg\mu_k\cos(\theta) \Rightarrow
    and
    a=g(sin(θ)μkcos(θ))a = g(\sin(\theta)-\mu_k\cos(\theta))

2.3Conservation of Momentum

From Newton’s 2nd^{\text{nd}} and 3rd^\text{rd} law we can easily derive the law of conservation of momentum.
Assume there are only two point-particle (i.e. particles with no size but with mass), that exert a force on each other. No other forces are present. From N2 we have:

dp1dt=F21dp2dt=F12\frac{d\vec{p}_1}{dt} = \vec{F}_{21} \\ \\ \frac{d\vec{p}_2}{dt} = \vec{F}_{12}

From N3 we know:

F21=F12\vec{F}_{21} = -\vec{F}_{12}

And, thus by adding the two momentum equations we get:

dp1dt=F21dp2dt=F12=F21}\left. \begin{array}{ll} \frac{d\vec{p}_1}{dt} &= \vec{F}_{21} \\ \\ \frac{d\vec{p}_2}{dt} &= \vec{F}_{12} = -\vec{F}_{21} \end{array} \right \} \Rightarrow
dp1dt+dp2dt=0ddt(p1+p2)=0\frac{d\vec{p}_1}{dt} + \frac{d\vec{p}_2}{dt} = 0 \rightarrow \frac{d}{dt} \left ( \vec{p}_1 + \vec{p}_2 \right ) = 0
p1+p2=const  i.e. does not depend on time\Rightarrow \vec{p}_1 + \vec{p}_2 = const ~\text{ i.e. does} \textit{ not } \text{depend on time}

Note the importance of the last conclusion: if objects interact via a mutual force then the total momentum of the objects can not change. No matter what the interaction is. It is easily extended to more interacting particles. The crux is that particles interact with one another via forces that obey N3. Thus for three interacting point particles we would have (with Fij \vec{F}_{ij} the force from particle i felt by particle j):

dp1dt=F21+F31dp2dt=F12+F32=F21+F32dp3dt=F13+F23=F31F32}\left. \begin{array}{ll} \frac{d\vec{p}_1}{dt} &= \vec{F}_{21} + \vec{F}_{31} \\ \\ \frac{d\vec{p}_2}{dt} &= \vec{F}_{12} + \vec{F}_{32}= -\vec{F}_{21} + \vec{F}_{32}\\ \\ \frac{d\vec{p}_3}{dt} &= \vec{F}_{13} + \vec{F}_{23}= -\vec{F}_{31} - \vec{F}_{32} \end{array} \right \}

Sum these three equations:

dp1dt+dp2dt+dp3dt=0ddt(p1+p2+p3)=0p1+p2+p3=const.  i.e. does not depend on time\frac{d\vec{p}_1}{dt} + \frac{d\vec{p}_2}{dt} + \frac{d\vec{p}_3}{dt} = 0 \rightarrow \frac{d}{dt} \left ( \vec{p}_1 + \vec{p}_2 + \vec{p}_3 \right ) = 0 \\ \\ \Rightarrow \vec{p}_1 + \vec{p}_2 + \vec{p}_3 = \text{const.}~\text{ i.e. does} \textit{ not } \text{depend on time}

For a system of NN particles, extension is straight forward.

2.3.1Momentum example

The above theoretical concept is simple in its ideas:

  • a particle changes its momentum whenever a force acts on it;
  • momentum is conserved;
  • action = - reaction.

But it is incredible powerful and so generic, that finding when and how to use it is much less straight forward. The beauty of physics is its relatively small set of fundamental laws. The difficulty of physics is these laws can be applied to almost anything. The trick is how to do that, how to start and get the machinery running. That can be very hard. Luckily there is a recipe to master it: it is called practice.

Solution to Exercise 12

Let’s do this one together. We follow the standard approach of IDEA: Interpret (and make your sketch!), develop (think ‘model’), evaluate (solve your model) and assess (does it make any sense?).

Interpret
Develop
Evaluate
Assess

First a sketch: draw what is needed, no more, no less.

2.4Forces & Inertia

Newton’s laws introduce the concept of force. Forces have distinct features:

  • forces are vectors, that is, they have magnitude and direction;
  • forces change the motion of an object:
    • they change the velocity, i.e. they accelerate the object
a=Fmdv=adt=Fdtm\vec{a} = \frac{\vec{F}}{m} \leftrightarrow d\vec{v} = \vec{a}dt = \frac{\vec{F}dt}{m}
  • or, equally true, they change the momentum of an object
dpdt=Fdp=Fdt\frac{d\vec{p}}{dt} = \vec{F} \leftrightarrow d\vec{p} = \vec{F}dt

Many physicists like the second bullet: forces change the momentum of an object, but for that they need time to act.

Momentum is a more fundamental concept in physics than acceleration. That is another reason why physicists prefer the second way of looking at forces.

2.4.1Inertia

Inertia is denoted by the letter m m for mass. And mass is that property of an object that characterizes its resistance to changing its velocity. Actually, we should have written something like mi m_i , with subscript i denoting inertia.

Why? There is another property of objects, also called mass, that is part of Newton’s Gravitational Law.

Two bodies of mass m1 m_1 and m2 m_2 that are separated by a distance r12r_{12} attract each other via the so-called gravitational force (r^12\hat{r}_{12} is a unit vector along the line connecting m1m_1 and m2m_2):

F12=Gm1m2r122r^12\vec{F}_{12} = - G \frac{m_1 m_2}{r^2_{12}}\hat{r}_{12}

Here, we should have used a different symbol, rather than m m . Something like mg m_g , as it is by no means obvious that the two ‘masses’ mi m_i and mg m_g refer to the same property. If you find that confusing, think about inertia and electric forces. Two particles with each an electric charge, q1q_1 and q2q_2, respectively exert a force on each other known as the Coulomb force:

FC,12=14πϵ0q1q2r122r^12\vec{F}_{C,12} = \frac{1}{4\pi\epsilon_0} \frac{q_1 q_2}{r^2_{12}}\hat{r}_{12}

We denote the property associated with electric forces by q q and call it charge. We have no problem writing

F=maFC=14πϵ0qQr2r^\vec{F} = m \vec{a} \\ \vec{F}_C = \frac{1}{4\pi \epsilon_0} \frac{q Q}{r^2} \hat{r}

We do not confuse qq by mm or vice versa. They are really different quantities: qq tells us that the particle has a property we call ‘charge’ and that it will respond to other charges, either being attracted to, or repelled from. How fast it will respond to this force of another charged particle depends on mm. If mm is big, the particle will only get a small acceleration; the strength of the force does not depend on mm at all. So far, so good. But what about mgm_g? That property of a particle that makes it being attracted to another particle with this same property, that we could have called ‘gravitational charge’. It is clearly different from ‘electrical charge’. But would it have been logical that it was also different from the property inertial mass, mim_i?

F=miaFg=GmgMgr2r^\begin{aligned} \vec{F} &= m_i \vec{a} \\ \vec{F}_g &= -G \frac{m_g M_g}{r^2} \hat{r} \end{aligned}

As far as we can tell (via experiments) mi m_i and mg m_g are the same. Actually, it was Einstein who postulated that the two are referring to the same property of an object: there is no difference.

2.4.1.1Force field

We have seen, forces like gravity and electrostatics act between objects. When you push a car, the force is applied locally, through direct contact. In contrast, gravitational and electrostatic forces act over a distance — they are present throughout space, though they still depend on the positions of the objects involved.

One powerful way to describe how a force acts at different locations in space is through the concept of a force field. A force field assigns a force vector (indicating both direction and magnitude) to every point in space, telling you what force an object would experience if placed there.

For example, the graph below at the left shows a gravitational field, described by Fg=GmMr2r^\vec{F}_g=G\frac{mM}{r^2}\hat{r}. Any object entering this field is attracted toward the central mass with a force that depends on its distance from that mass’s center.

Source
import numpy as np
import matplotlib.pyplot as plt

# Define the grid
r = np.linspace(1.4, 3, 3)
theta = np.linspace(0, 2*np.pi, 11)

R, T = np.meshgrid(r, theta)

X = np.cos(T) * R
Y = np.sin(T) * R

# Gravitational field at a point (r,theta)
def gravitational_field(R, T):
    # Magnitude 
    F = -1 / (R**2) 

    # X-component
    F_x = F * np.cos(T)
    # Y-component
    F_y = F * np.sin(T)

    return F_x, F_y

# Calculate the field
F_x, F_y = gravitational_field(R, T)

# Create a figure with a fixed aspect ratio
fig, axes = plt.subplots(1, 2, figsize=(10, 5))

# Plot gravitational field
axes[0].set_xlim((-3, 3))
axes[0].set_ylim((-3, 3))
axes[0].quiver(X, Y, F_x, F_y, color='#00a6d6ff', scale=3.5, width=0.004)
axes[0].scatter(0, 0, s=800, c='k')
axes[0].set_aspect('equal', 'box')
axes[0].set_xticks([])
axes[0].set_yticks([])
axes[0].set_xlabel('x')
axes[0].set_ylabel('y')

# Simulation parameters
electron1 = np.array([-1, 0])  # x, y position of the electron
electron2 = np.array([1, 0])

# Define the grid
x = np.linspace(-2, 2, 400)
y = np.linspace(-2, 2, 400)

X, Y = np.meshgrid(x, y)

# Electric field at a point (x, y)
def electric_field(x, y):
    r1 = np.sqrt((x - electron1[0])**2 + (y - electron1[1])**2)
    r2 = np.sqrt((x - electron2[0])**2 + (y - electron2[1])**2)

    # Electric field due to electron1
    E1_x = - (x - electron1[0]) / r1**3
    E1_y = - (y - electron1[1]) / r1**3

    # Electric field due to electron2
    E2_x = - (x - electron2[0]) / r2**3
    E2_y = - (y - electron2[1]) / r2**3

    # Total electric field
    E_x = E1_x + E2_x
    E_y = E1_y + E2_y

    return E_x, E_y

# Calculate the electric field
E_x, E_y = electric_field(X, Y)

# Plot electric field lines
axes[1].set_xlim((-2, 2))
axes[1].set_ylim((-2, 2))
axes[1].streamplot(X, Y, E_x, E_y, color='#00a6d6ff', density=1, linewidth=1, zorder=1)
axes[1].scatter(electron1[0], electron1[1], s=200, c='k', zorder=10)
axes[1].scatter(electron2[0], electron2[1], s=200, c='k', zorder=10)
axes[1].text(electron1[0], electron1[1], '-', color='white', ha='center', va='center', fontsize=20, zorder=11)
axes[1].text(electron2[0], electron2[1], '-', color='white', ha='center', va='center', fontsize=20, zorder=11)
axes[1].set_aspect('equal', 'box')
axes[1].set_xticks([])
axes[1].set_yticks([])
axes[1].set_xlabel('x')
axes[1].set_ylabel('y')

plt.tight_layout()
plt.show()
<Figure size 1000x500 with 2 Axes>

2.4.1.2Measuring mass or force

So far we did not address how to measure force. Neither did we discuss how to measure mass. This is less trivial than it looks at first side. Obviously, force and mass are coupled via N2: F=ma F = m a .

Can force be measured using a balance?

Figure 13:Can force be measured using a balance?

The acceleration can be measured when we have a ruler and a clock, i.e. once we have established how to measure distance and how to measure time intervals, we can measure position as a function of time and from that velocity and acceleration.

But how to find mass? We could agree upon a unit mass, an object that represents by definition 1kg. In fact we did. But that is only step one. The next question is: how do we compare an unknown mass to our standard. A first reaction might be: put them on a balance and see how many standard kilograms you need (including fractions of it) to balance the unknown mass. Sounds like a good idea, but is it? Unfortunately, the answer is not a ‘yes’.

As on second thought: the balance compares the pull of gravity. Hence, it ‘measures’ gravitational mass, rather than inertia. Luckily, Newton’s laws help. Suppose we let two objects, our standard mass and the unknown one, interact under their mutual interaction force. Every other force is excluded. Then, on account on N2 we have

{m1a1=F21m2a2=F12=F21\left\{ \begin{array}{l} m_1 a_1 = F_{21} \\ m_2 a_2 = F_{12} = -F_{21} \end{array} \right.

where we used N3 for the last equality. Clearly, if we take the ratio of these two equations we get:

m1m2=a2a1\frac{m_1}{m_2} = \left | \frac{a_2}{a_1} \right |

irrespective of the strength or nature of the forces involved. We can measure acceleration and thus with this rule express the unknown mass in terms of our standard.

Now that we know how to determine mass, we also have solved the problem of measuring force. We just measure the mass and the acceleration of an object and from N2 we can find the force. This allows us to develop ‘force measuring equipment’ that we can calibrate using the method discussed above.

2.4.2Eötvös experiment on mass

The question whether inertial mass and gravitational mass are the same has put experimentalists to work. It is by no means an easy question. Gravity is a very weak force. Moreover, determining that two properties are identical via an experiment is virtually impossible due to experimental uncertainty. Experimentalist can only tell the outcome is ‘identical’ within a margin. Newton already tried to establish experimentally that the two forms of mass are the same. However, in his days the inaccuracy of experiments was rather large. Dutch scientist Simon Stevin concluded in 1585 that the difference must be less than 5%. He used his famous ‘drop masses from the church’ experiments for this (they were primarily done to show that every mass falls with the same acceleration).

A couple of years later, Galilei used both fall experiments and pendula to improve this to: less than 2%. In 1686, Newton using pendula managed to bring it down to less than 1‰ .

An important step forward was set by the Hungarian physicist, Loránd Eötvös (1848-1918). We will here briefly introduce the experiment. For a full analysis, we need knowledge about angular momentum and centrifugal forces that we do not deal with in this book.

2.4.2.1The experiment

The essence of the Eötvös experiment is finding a set up in which both gravity (sensitive to the gravitational mass) and some inertial force (sensitive to the inertial mass) are present. Obviously, gravitational forces between two objects out of our daily life are extremely small. This will be very difficult to detect and thus introduce a large error if the experiment relies on measuring them. Eötvös came up with a different idea. He connected two different objects with different masses, m1m_1 and m2m_2, via a (almost) massless rod. Then, he attached a thin wire to the rod and let it hang down.

Torsion balance used by Eötvös.

Figure 15:Torsion balance used by Eötvös.

This is a sensitive device: any mismatch in forces or torques will have the setup either tilt or rotate a bit. Eötvös attached a tiny mirror to one of the arms of the rod. If you shine a light beam on the mirror and let it reflect and be projected on a wall, then the smallest deviation in position will be amplified to create a large motion of the light spot on the wall.

In Eötvös experiment two forces are acting on each of the masses: gravity, proportional to mgm_g, but also the centrifugal force Fc=miRω2F_c = m_i R \omega^2, the centrifugal force stemming from the fact that the experiment is done in a frame of reference rotating with the earth. This force is proportional to the inertial mass. The experiment is designed such that if the rod does not show any rotation around the vertical axis, then the gravitational mass and inertial mass must be equal. It can be done with great precision and Eötvös observed no measurable rotation of the rod. From this he could conclude that the ratio of the gravitational over inertial mass differed less from 1 than 5108 5 \cdot 10^{-8}. Currently, experimentalist have brought this down to 11015 1 \cdot 10^{-15}.