Physics Lab Experiments (2016)
Physics Lab Experiments (2016)
EXPERIMENTS
LICENSE, DISCLAIMER OF LIABILITY, AND LIMITED WARRANTY
By purchasing or using this book (the “Work”), you agree that this license
grants permission to use the contents contained herein, but does not give you
the right of ownership to any of the textual content in the book or ownership
to any of the information or products contained in it. This license does not
permit uploading of the Work onto the Internet or on a network (of any kind)
without the written consent of the Publisher. Duplication or dissemination
of any text, code, simulations, images, etc. contained herein is limited to and
subject to licensing terms for the respective products, and permission must
be obtained from the Publisher or the owner of the content, etc., in order to
reproduce or network any portion of the textual material (in any media) that
is contained in the Work.
The sole remedy in the event of a claim of any kind is expressly limited to
replacement of the book, and only at the discretion of the Publisher. The
use of “implied warranty” and certain “exclusions” vary from state to state,
and might not apply to the purchaser of this product.
PHYSICS LAB
EXPERIMENTS
Matthew M. J. French
Original title and copyright: Finding Your Way in the Undergraduate Physics Labo-
ratory. Copyright ©2013 by M.M.J. French. ISBN 978-0-9549780-8-2. All rights
reserved. Published by The Pantaneto Press.
This publication, portions of it, or any accompanying software may not be repro-
duced in any way, stored in a retrieval system of any type, or transmitted by any
means, media, electronic display or mechanical display, including, but not limited
to, photocopy, recording, Internet postings, or scanning, without prior permission
in writing from the publisher.
The publisher recognizes and respects all marks used by companies, manufactur-
ers, and developers as a means to distinguish their products. All brand names and
product names mentioned in this book are trademarks or service marks of their
respective companies. Any omission or misuse (of any kind) of service marks or
trademarks, etc. is not an attempt to infringe on the property of others.
Our titles are available for adoption, license, or bulk purchase by institutions, corpo-
rations, etc. For additional information, please contact the Customer Service Dept.
at 800-232-0223(toll free).
All of our titles are available in digital format at authorcloudware.com and other
digital vendors. The sole obligation of MERCURY LEARNING AND INFORMATION to the
purchaser is to replace the book, based on defective materials or faulty workman-
ship, but not based on the operation or functionality of the product.
CONTENTS
Chapter 1 Introduction 1
1.1 Characteristics of the Laboratory 2
1.2 Demonstrating Undergraduate Physics Laboratory 4
Chapter 2 Safety 9
2.1 Hazard Symbols 10
2.2 Electrical 12
2.3 High Voltages 13
2.4 Gas Cylinders 15
2.5 Vacuum Pumps 19
2.5.1 Pressure Units 20
2.5.2 Hand Pumps and Vacu-vins 22
2.5.3 Diaphragm Pump 22
2.5.4 Rotary Pump 22
2.5.5 Diffusion Pump 23
2.5.6 Turbomolecular Pump 23
2.6 Cryogenics 25
2.6.1 Liquid Nitrogen 25
2.6.2 Dry Ice 30
2.6.3 Liquid Helium 32
2.6.4 Helium-3 33
2.6.5 Liquid Oxygen 35
2.6.6 Argon 36
2.6.7 Cryogen Free Cooling 37
2.7 Lasers 38
2.8 Ionizing Radiation 40
3.7 Ammeters 50
3.8 Voltmeters 51
3.9 Resistor Color Codes 52
3.10 Multimeters 52
3.11 Further Ideas on Measurement of Resistance 53
3.12 Constant Current Sources 55
3.13 Lock-in Amplifers 58
3.14 Further Ideas on Reducing Noise in Electrical Measurements 61
3.15 Signal Generator 64
3.15.1 TTL 64
3.16 Oscilloscopes 65
3.17 Light Gates 70
3.17.1 Velocity with One Light Gate 70
3.17.2 Velocity with Two Light Gates 71
3.17.3 Acceleration with Two Light Gates 71
3.17.4 Acceleration with One Light Gate 71
3.18 PASCO Science Workshop Interface 72
3.19 Data Harvest 74
3.20 Pressure Gauges 77
3.21 Temperature Measurement 78
3.21.1 Near Room Temperature 79
3.21.2 Low Temperatures 82
3.21.3 High Temperatures 84
3.22 PID and Temperature Control 84
3.23 Kaye and Laby Tables 86
3.24 Microscopes 87
3.24.1 Optical Microscopes 87
3.24.2 Transmission Electron Microscope 89
3.24.3 Other Microscopes 90
3.25 Spectroscopes 92
3.26 Computers 92
3.26.1 Plugs, Sockets, and Connectors 92
3.26.2 Computer Data Logging 97
3.26.3 Still Cameras and Scanners 105
3.27 Video Cameras 106
3.28 Other Researcher's Data 107
3.29 Wii Controllers 108
3.30 Traditional Photographs 109
Index 275
ABOUT THE AUTHOR
1
INTRODUCTION
In the first year, students will usually work in pairs (or very
small groups) undertaking a different experiment each week. It
is important students are able to work with other scientists: for
some experiments it is essential to have two people to perform the
experiment. By discussing experimental methods and the prob-
lems experienced with their partner students can challenge their
understanding and learn from each other. However, it is essential
to work together as a team and that both are engaged in each part
Introduction • 3
Figure 1.1: The undergraduate laboratories at the University of Oxford. Left: condensed
matter laboratory. Right: optics laboratory.
2
SAFETY
Table 2.1: Chemical Hazard Symbols. Usually the border is red and the diagram is black.
Table 2.2: Non-Chemical Hazard Symbols. Usually the text, border, and diagram is black
and the background is yellow.
2.2 ELECTRICAL
Figure 2.1: Examples of simple low EMF power supplies. Usually the yellow outputs give AC.
The black (negative) and red (positive) outputs give DC.
14 • Physics Lab Experiments
Figure 2.2: Left: Low (electric) potential difference connectors with exposed metal contacts.
Right: High (electric) potential difference connectors with shielded metal contacts. The
shielding automatically moves back as the plug is connected into a socket. Both leads can
be piggy-backed to connect multiple leads to the same place.
Figure 2.5: A cart used to move gas cylinder. This cart has an extra set of wheels so the user
doesn’t have to support part of the weight of the cylinder while moving the cylinder.
Safety • 17
Figure 2.6: A gas regulator. This one is for oxygen only: the knob is labeled and the crossed
out symbol of an oil can is present in the pressure gauges.
18 • Physics Lab Experiments
check the screw threads on the cylinder and regulator are in good
condition. With closed eyes, blow into the cylinder head to remove
any dust.
Then place the regulator above the thread and using a spanner,
turn the nut (clockwise from above for right handed threads) until
the nut will not turn any further. Do not use excessive force.
To obtain gas, ensure the knob is turned fully counterclockwise (the
off position). Use the cylinder key to turn the spindle (see Figure 2.7)
fully counterclockwise (the on position). The pressure gauge now gives
a reading, indicating the contents of the cylinder. As gas is used the pres-
sure reading will decrease: an empty cylinder will read approximately
zero when the spindle is open. A quick check for leaks in the regulator
can be made. Gas leaking rapidly will make a hissing noise. Slowly turn
the regulator knob clockwise (to the on position) to obtain the gas. The
second gauge indicates the delivery pressure of the gas. The further the
knob is turned clockwise the higher the delivery pressure.
To close down a gas cylinder when the experiment is finished
turn the regulator knob fully counterclockwise to the off position.
Use the cylinder key to move the spindle fully clockwise to the
off position.
Next turn the knob clockwise to remove any gas in the regulator.
This prevents damage to the diaphragm in the regulator. The pres-
sure gauge will now read zero. Finally turn the knob fully counter-
clockwise again.
To remove the regulator when the cylinder is empty and needs
changing follow the instructions in the paragraph above to close
down a gas cylinder. Then unscrew (counterclockwise from above
Figure 2.7: Left: A gas cylinder key. Right: The spindle on a gas bottle.
Safety • 19
for a right handed thread) the nut using a spanner. The regulator
should now be free. It is vital that the spindle is fully turned off
(turned fully clockwise) before starting to unscrew the nut. This pre-
vents gas escaping rapidly from the cylinder potentially causing the
torpedo effect described earlier.
Note carefully that the spindle is off when turned fully clockwise
where as the regulator is off when turned fully counterclockwise.
Compressed gases are graded by their purity which is usually
expressed as a percentage, for example: 99 % or 99.999 %. How-
ever this is often re-expressed in terms of the number of nines:
99 % becomes 2N (corresponding to 2 nines) and 99.999 % becomes
5N (corresponding to 5 nines). A further variation is 99.95 % which
becomes 3N5 (corresponding to 3 nines and then a five).
In the EU, the colors of cylinders of some gases are becoming
standardized following the new standard EN 1089-3. Colors for cyl-
inders of the most common gases are given in Table 2.3.
Rotary pumps are the work horse pump of the physics lab. They
provide a cheap, reliable way to achieve a vacuum down to a pres-
sure of around 0.1 mbar. A two stage rotary pump consists of two
rotary pumps in series and will produce a vacuum of around 0.001
mbar = 1×10−3 mbar. Figure 2.9 shows a rotary pump.
The pump consists of a vane projecting from a cylinder rotating
eccentrically which compresses air and forces it out through a valve.
Rotary pumps have oil inside. This is a special type of oil which has a
very low vapor pressure i.e., it produces very little vapor. If it did not
have a low vapor pressure it could not be used as the oil vapor would
pollute the vacuum the pump was trying to produce. Despite the
low vapor pressure tiny droplets (i.e., a mist) of oil from the pump
can be carried away by the air it exhausts. This is more of a concern
when the pumping flow rate is high or the pump is run for pro-
longed periods of time or in a confined space/small room. There are
Safety • 23
2.6 CRYOGENICS
Figure 2.11: Phase diagram for Nitrogen. Nitrogen is a gas at STP (standard temperature
and pressure).
Figure 2.13: Liquid Nitrogen Dewars. Left: a 25 liter “onion” Dewar with a cart Right: a 10
liter Dewar.
reaction experienced when hot objects are touched) and present the
danger severe cold burns.
Care should be taken to avoid splashes of liquid nitrogen enter-
ing the eyes especially when pouring liquid nitrogen from one con-
tainer to another and so wearing goggles may be appropriate. Open
shoes or sandals must be avoided by all those working in the same
room as liquid nitrogen is used.
A further danger presented by liquid nitrogen is asphyxiation.
When a liquid turns into a gas it expands significantly (remember
the diagrams showing the particles in a solid, liquid, gas drawn
Safety • 29
where Vr is the room volume, L is the gas release units of m3/h and
n is the number of air changes in the room per hour. Assuming a
typical value of n=0.4 and that L is twice the manufacturers quoted
value of 0.2 liquid liters per day (i.e. L=0.4 × 683 = 273 gas liters
per day = 11.4 gas liters per hour = 0.0114 m3/h where the gas
expansion factor is 683) gives the minimum room volume to not
fall below 19.5 % oxygen as 0.38 m3. This is insignificant compared
to the requirement concerning the complete spillage of a Dewar.
Following Annex 4 of the same code, accounting for the complete
spillage of the 25 liter dewar the minimum room volume for stor-
age/use will be given by
30 • Physics Lab Experiments
Vd f g
0.21 Vr –
1000 (2.2)
C
Vr
The primary danger from handling dry ice is severe frost bite,
therefore dry ice should be handled only while wearing thick gloves
(suede gardening gloves or special cryogenic gloves may be appro-
priate). Similar gloves should also be used while handing objects
cooled by dry ice. With extreme care, individual pieces of dry ice my
be handled without gloves provided they are not gripped tightly nor
held for more that a second or two.
As with liquid nitrogen, care must be taken over storage espe-
cially in small room or confined spaces. The typical carbon dioxide
content of the atmosphere is around 390 ppm (parts per million) [4].
The safe limit is considered 5000 ppm, above this level symptoms
consist of an increased breathing rate, tiredness, and headaches. If
any of the symptoms appear immediately remove the affected person
to the open air and follow up with artificial respiration if necessary.
Figure 2.14: Phase diagram for Carbon Dioxide. Carbon Dioxide is a gas at STP
(standard temperature and pressure).
32 • Physics Lab Experiments
2.6.4 HELIUM-3
Figure 2.16: Phase diagram for Oxygen. Oxygen is a gas at STP (standard temperature and
pressure).
2.6.6 ARGON
2.7 LASERS
The older classification system which has now been superseded may
still be in use on older devices. This is shown in Table 2.8.
Figure 2.18: Left panel: a lead lined wooden storage box and forceps. Right panel: inside
the box and the source inside a metal cup holder.
Table 2.9: List of common radioactive sources. Protactinium and Radon are listed with
radiation given out by the decay chain which generates them from Uranium-238 and Tho-
rium-232 respectively.
3
MAKING
MEASUREMENTS
3.1 SI UNITS
3.2 SI PREFIXES
Figure 3.1: A Vernier Scale. The top is the main scale and the bottom in the traveling scale.
The left panel shows a reading of 0.0 mm; The middle panel shows a reading of 0.1 mm;
The right panel shows a reading of 0.6 mm.
Figure 3.2: A digital calipers measuring; Left: the diameter of a wire and Right: the inside
diameter of a pipe.
3.5 MICROMETER
3.6 BALANCES
precisions are available and it should always be ensured that the bal-
ance used is appropriate for the measurement needed. For example:
it would not be appropriate to use a 0–2 kg balance which reads
to the nearest 1 g to try to measure 100 mg of metal, neither is it
appropriate to check the mass of a 1 kg standard mass with a 0 - 500
g balance which reads to the nearest 1mg.
Balances involve the use of one or more pressure sensors which
produce an electrical output dependent upon the force acting on
them. This electrical output has been calibrated by the manufac-
turer as corresponding to a particular force. This force is then con-
verted to a mass by dividing by an average value for the acceleration
due to gravity. According to the latest data from the NASA Grace
mission the acceleration due to gravity varies around the Earth by at
least 0.001 ms−2. Thus quoting masses to more than 4 significant fig-
ures is inappropriate on balances not calibrated for use in a specific
location, although this may not be true for mass differences.
The balance should be positioned directly on a flat surface: ensure
there are no wires or other objects under the balance. When mak-
ing measurements on a table, desk or worksurface ensure that no one
touches, leans or sits on the surface as this can sometimes significantly
alter the reading on the scale. It is useful to perform a quick experi-
ment to see the significance of this. Some balances have a cover to
prevent draughts (possibly from people walking past or convection
currents from heating or air conditioning) affecting the results: this
should be placed down and closed when any readings are taken.
To take a measurement first ensure that the scale pan is clean and
free from loose material and the balance reads zero. If not the “Tare”
or “Zero” button should be pressed. If a container is used often to
hold powder which is being measured out, this can be placed on the
balance before it is zeroed. Finally add the object to be measured to
the scale pan and take a reading.
3.7 AMMETERS
minimum (or no) change to the behavior of the rest of the circuit an ideal
ammeter will allow current to flow through with no resistance. In prac-
tice ammeters have a very low resistance, which is assumed to be small
compared to the circuit being measured. To picture this more clearly
consider the equation for the total resistance of resistors in series:
RT R Rammeter (3.1)
3.8 VOLTMETERS
minimum (or no) change to the behavior of the rest of the circuit an
ideal voltmeter will have an infinite resistance and allow no current
to flow through it. In practice, voltmeters have a very high resis-
tance, which is assumed to be large compared to the circuit being
measured. To picture this more clearly consider the equation for the
total resistance of resistors in parallel:
1 1 1
(3.2)
RT R Rvoltmeter
3.10 MULTIMETERS
Figure 3.5: Simple diagram of a four-wire method of measuring the sample resistance. R1,
R2, R3, and R4 are resistances which represent the contact and lead resistances.
V V
I (3.4)
R Rseries Rsample
Vin V
– B (3.5)
Ri Rf
Where Vin is the (electric) potential difference from the lock-in oscil-
lator output; Ri and Rf are the values of the resistors and VB is the
Making Measurements • 57
Figure 3.7: Left panel shows a circuit layout of a resistor (typically Rseries = 1 kΩ to 10 kΩ)
used as a constant current source. It is connected in series with the sample and the oscillator
output on a lock-in amplifier. Right panel shows a circuit diagram of a (electric) potential
difference controlled constant current source. The op-amp used is an Linear Technology
LTC1150CN8#PBF (RS order number 5455629).
VB V
Isample B (3.6)
Rg Rf
R f Rg
Isample – V (3.7)
RiRg in
The lock-in oscillator output is connected to Vin and the sample cur-
rent contacts are those attached to the op-amp output. The current
in the sample is given by Equation 3.7.
A final method of generating a constant current is to use a dedi-
cated off the shelf constant current source such as a Keithley 6221
AC and DC current source. This has two current output leads which
are connected to the sample.
58 • Physics Lab Experiments
Vpsdoutput V
f
f ,sig VLsin w f t q f ,sig sin wr t qref
(3.8)
Vpsdoutput 2V
1
sig
VL cos wr – w f t qsig – qref
f
(3.9)
–
1
V V cos wr w f t qsig – qref
2 sig L
For each component frequency, f, that is detected this output is two
AC signals, one at the difference frequency wr−wf and the other at the
sum frequency wr+wf, which are then passed through a low pass filter.
In the general case for wr and wf there is no signal from the low pass
filter as both signals are AC and are filtered out. If wr=wf then part of
the signal is no longer AC, but DC and proportional to the signal ampli-
tude Vlpfoutput 2 Vsig VL cos qsig – qref . This gives the X output of the
1
Figure 3.8: The front panel of a Stanford Research Systems (SRS) 830 lock-in amplifier.
60 • Physics Lab Experiments
Figure 3.9: Left panel shows the circuit setup for a normal 4-wire measurement configura-
tion. Right panel shows the circuit setup when setting the phase. Current is passed though
the sample and a series resistor, R. The (electric) potential difference drop is measured
across only the resistor.
Z X2 Y2 (3.10)
R 2 XC2 X L2 (3.11)
R
L (3.12)
Figure 3.10: Inductive cross-talk between wires in two closed circuits carrying AC signals.
Making Measurements • 63
m0 I sin wr t
B r, t (3.13)
2p r
The other two wires are part of a neighboring circuit and situated
parallel to the first wire a distance of r1 and r2 away. The separation
of these two wires (r2−r1) continues for length l. Ampere’s law can
be used to show the magnetic flux induced through this area of the
circuit is
B r, t dA
S
(3.14)
r2 m0 I sin wr t
l
r1 2p r
0
dzdr (3.15)
l m0 I sin wr t r2
ln (3.16)
2p r r1
3.15.1 TTL
3.16 OSCILLOSCOPES
To begin using an oscilloscope first find and operate the power but-
ton on the front panel. It is then usual for a power LED to turn on
and perhaps a number of other LEDs. If nothing happens, check for
an additional power switch on the back of the oscilloscope and that
the mains lead is properly connected and that the mains is turned
on.
Wait for the beam to be displayed. An analogue oscilloscope may
take a moment or two to turn on and a digital one will take longer as
the software starts up. If the trace is not displayed on the screen then
check the trigger level dial is set to the center and the horizontal and
vertical position controls are set to the center.
Next connect the signal via the BNC connector on the front
panel to at least one channel. The coupling should be set to DC in
most cases and the mode set to either CH1, CH2, or DUAL depend-
ing on where the input signal(s) are connected. Set the vertical gain
and timebase dials to an appropriate setting so that the waveform
occupies a large part
To measure the acceleration with just one light gate, two veloc-
ity measurements are still needed. These can be taken using a “U”
shaped card. The first velocity is measured using the first card to
break the beam and the second velocity is measured using the sec-
ond card to break the beam. Here it is important to measure the
width of each of the vertical sections of the “U”. The acceleration is
calculated in a similar way by diving the difference in the two veloci-
ties by the duration of time between them.
72 • Physics Lab Experiments
Figure 3.14: A PASCO Science WorkShop 500 Interface box connected to a photogate.
Making Measurements • 73
Figure 3.19. Alternatively, in a stand alone mode, the results are dis-
played on the built in screen or recorded and saved to the datalogger.
For example to setup a velocity measurement using a single light
gate, connect the light gate to the input labeled A. Press any button
to turn the device on. The arrow buttons can be used to navigate
through different menu options. Navigate to “Time & Motion” and
press the green arrow button. Navigate to “Speed” and press the
green arrow button. Navigate to “Speed at A” and press the green
arrow button. Choose the “1 Interrupt card” option and press the
green arrow button. Finally navigate to the length of the interrupt
card being used and press the green button. The datalogger will dis-
play the latest speed as well as an average speed.
To use a Light Level probe, connect this to input 1. Navigate to
“Meter” and press the green arrow button. The display will show a
reading updated a few times per second.
Use the red square button to go back to the previous menu. Nav-
igate to “Logging” and press the green arrow. Use the arrow keys
Making Measurements • 77
P1 – 2 r gh (3.19)
TF 1.8TC 32 (3.20)
TK TC 273.15 (3.21)
2h f 3 1
I( f , T) (3.22)
c2 ehf / kT – 1
P AesT 4. (3.23)
t d
x t Pe t I e t dt D dt e t
0
(3.24)
half the value which caused the oscillation. The present value will
then overshoot the setpoint. The I value can then be increased from
zero until the offset is corrected in a reasonable time. Finally the D
value can be increased so that any instability in the system is cor-
rected in a reasonable time.
Table 3.5 gives the effect of changing one of the PID constants
independently: An alternative to setting the PID gain constants
manually is to follow one of a number of pre-defined algorithms for
setting them such as the Ziegler–Nichols method [8]. As before the I
and D gain constants are set to zero. The P gain constant is increased
until it reaches the ultimate gain PU at which the output of the loop
oscillates with a constant amplitude. The values of PU and TU (where
TU is the time period of oscillation) are used to set the values of the
PID gain constants using the relationships shown in Table 3.6 [9].
3.24 MICROSCOPES
Figure 3.20: Illustration of chromatic aberation of a lens. Red light (the dashed line) is
refracted less than blue light (the dotted line).
Figure 3.22: Left: Electron diffraction pattern from Gold. Right: A dark field image of
islands of Gold.
3.25 SPECTROSCOPES
3.26 COMPUTERS
Mains
Mains
Mains
3 Phase
Mains
VGA
DVI
S-Video
94 • Physics Lab Experiments
Composite
HDMI
3.5mm
Audio
USB Type A
USB Type B
mini USB
Making Measurements • 95
USB 3
miniUSB 3
FireWire
(IEEE1394)
mini
FireWire
BNC
9pin Serial
(RS-232)
96 • Physics Lab Experiments
25pin Serial
(RS-232)
Parallel
Parallel
GPIB
(IEEE488)
PS-2
Ethernet
(RJ45)
Modem
(RJ11)
the GPIB bus in order to respond with the required data. These
sometimes complex commands can be written out once in Delphi in
an interface unit specific to that instrument. Each Delphi program
that is then written uses these interface units to communicate with
instruments. A standard interface unit (called gpib) is also required
in order for the Delphi program to communicate with the GPIB
card in the computer. All the steps described in the this section need
only be performed once. After this is all setup it provides a simple
and rapid method to write new programs.
As an example consider reading the X output from a SR830 lock-
in amplifier. A simple interface unit would be as follows below. The
unit name (DCSTAN830) is given at the top, then there are two main
sections “interface” and “implementation”. The interface section con-
tains definitions of class properties, global variables and all procedures
contained in the unit. The first line of the implementation section
contains a list of other units from which code is referenced rather
than duplicated. The code at the very end encapsulated with a “begin”
and “end“ statement creates the class and sets a text handle which
corresponds to the text handle linked to the correct GPIB address
set earlier in the CBCONF32 program. There are two procedures
also in the implementation section, the first initialises the instrument
and the second “GetVolts” sends a message over the GPIB interface
(“SNAP?1,2,3,4”) and waits for the response to come back before
interpreting it and saving part of it as the variable x. These messages
and the form of the response are detailed for each instrument in its
manual.
unit DCSTAN830;
interface
Procedure TLIA.Initialise;begin
Making Measurements • 99
ieeeHdl:=InitDevice(ieeeStr);
ibclr(ieeeHdl);
writeieee(ieeeHdl,’OUTX1’);
writeieee(ieeeHdl,’*CLS’);
end;
begin LIA[1]:=TLIA.Create;
LIA[1].ieeeStr:=’SR830’;
end.
As well as the interface unit the actual program is needed - the sim-
ple one shown below sets a label on the Delphi form to show the
lock-in X value when a button is pressed. At the top in the “uses” sec-
tion the name of the interface unit (DCSTAN830) has been added:
this automatically includes the code from the file above to save it
being typed out again. The “FormCreate” function is run automati-
cally when the program is opened and runs the instrument initialisa-
tion procedure. The “Button1Click” function is run when the button
on the form is clicked. This first runs the “GetVolts” procedure to
get the value from the instrument and then sets the label caption to
be the value returned and saved to the variable x.
unit Unit1;
interface
uses
Windows, Messages, SysUtils, Variants, Classes,
Graphics,
Controls, Forms,
Dialogs, StdCtrls, dcstan830, dmisc;
type
TForm1 = class(TForm)
Button1: TButton;
Label1: TLabel;
procedure Button1Click(Sender: TObject);
private
Private declarations
100 • Physics Lab Experiments
public
Public declarations
end;
var
Form1: TForm1;
implementation
$R *.dfm
LabVIEW
LabVIEW software gives a graphical point and click interface
which allows acquisition and processing of data and data logging. It
also has the ability to make decisions based on the measurements
and control equipment. LabVIEW programs can be run on Win-
dows, Mac and Linux.
A LabVIEW program is called a VI. When LabVIEW is started,
from the “Create Project’ window select the “Blank VI” template:
see Figure 3.23. This brings up two windows: one shows the “block
diagram” where the code is developed and the other shows the
“front panel” which is where the user interface can be customized to
include graphs and buttons: see Figure 3.24.
In the “block diagram” window, pictogram components can
be wired together to create a program. Each block has inputs and
outputs just like a function in a traditional programming language.
When a component executes it produces data which passes down
the wire to the next block. The movement of data determines the
order in which the components are executed in the program.
102 • Physics Lab Experiments
Figure 3.24: The LabVIEW program. On the left is the ‘block diagram’ window and on the
right is the ‘front panel’ window.
Figure 3.26: A LabVIEW program which reads data and produces a graph.
Making Measurements • 105
ticularly useful for measuring the speed of objects which are moving
or events which happen too rapidly to be seen with the human eye.
The experiment should be setup so that a ruler or scale is in camera
shot and in close proximity to the object. Alternatively a resolution cal-
ibration must be performed giving a value in units of meters per pixel.
Captured video can be opened in a program such as Apple
QuickTime. The video should be played back frame by frame. The
position of the object can be read from the ruler or by measuring
its position in terms of pixels from a fixed point. When the same
measurements are taken from the next frame the velocity can be
calculated using:
the y axis. For each point input the coordinates based on the scale
printed on the graph. This is shown in Figure 3.29.
The program allows the user to either identify points manually
or to automatically detect and trace a curve. Both methods produce
a table of x and y coordinates which can be exported to a graph draw-
ing program.
Over recent year there has been great interest in using Nintendo
Wii games console controllers as data loggers. Wii controllers have
a built in 3 axis accelerometer and an infrared camera which can be
used to give the distance from a pair of infrared sources. They can
communicate with computers via their built in bluetooth radio.
Early work [11],[12] involved a significant amount of computer
programming. Dedicated and easy to use “Wiimote” software has
now been written [13],[14] to allow the connection of Wii controller
with a computer and subsequent collection of data in real time.
The computer must have a Bluetooth radio adapter. Sometimes
these are built into computers (especially laptops), but a cheap
Making Measurements • 109
USB dongle can also be purchased. The more stable and reliable
ones tend to be those which use the Microsoft Bluetooth Stack
drivers. Once the computer is setup, the Wii controller needs to
be paired by pressing buttons 1 and 2 together and initiating the
pairing via the computer. If a passcode is requested in Windows
it can be left blank or skipped by pressing Alt-S on the computer
keyboard.
Once the Wii controller is paired, the Wiimote software can be
started. In accelerometer mode, the acceleration can be recorded
in the x, y and z axis. It measures in units of g between ± 3g with a
resolution of 0.04 g. When the controller is stationary lying flat on
the desk, one of the axes will measure 1 g. This can be used to mea-
sure the acceleration of a dynamics cart, a pendulum, a mass spring
system, a lift, a car or a rotating disk.
The software can also be set in “position sensor” mode. The bar
which usually sits on top of the TV is commonly referred to as the
“Sensor Bar”. In fact, it does not sense anything, but is two infrared
LEDs positioned 20cm apart. These two sources emit infrared light
which is detected by the 1024 pixel by 768 pixel infrared camera
in the Wii controller. Given the angular fields of view of the cam-
era are 41° horizontally and 31° vertically it is possible to find the
average angle subtended per pixel as 0.040°. The software in the
Wii is able to measure the number of pixels, n, between the two
light sources on its infrared camera and convert this into the angu-
lar separation. Since the sources are known to be 20 cm apart, the
distance, d, from the WiiMote to the Sensor Bar can be calculated
using triangulation:
0.1
d (3.27)
tan(0.020 n)
fixer solutions has a large effect on the final image and will depend
on their strength. The correct timings will usually be given, however
as a guide consider around 1 minute in the developer, 30 seconds in
the first wash, 2 minutes in the fixer, and 2 minutes in the final wash.
To create a paper photograph a second stage of development is
required. The developed film is inserted into a slot below the light
source on an enlarger. The light shining through the film creates an
image on the easel. Once the size and position of the image have
been adjusted the light on the enlarger is turned off. A sheet of pho-
tographic paper is placed onto the easel. The light source can then
be turned on using a timer, exposing the photographic paper to the
image. The final stage is to develop the photographic paper in a simi-
lar way to the development of the film. The developer solution is
usually specific to the paper, but the same two wash baths and fixer
bath as for the negative can usually be used.
Always carefully mop up any chemicals which are spilt and clear
away all chemicals at the end of each session/day as they can leave
residues behind when they evaporate. Chemicals should not be
returned to the original bottle as this will cause contamination.
CHAPTER
4
DATA ANALYSIS
AND ERRORS
l
T 2p (4.1)
g
Figure 4.1: Targets illustrating the differences between the terms accurate and precise.
The center of the target corresponds to the true value.
The stars are the individual measurements.
116 • Physics Lab Experiments
Figure 4.2: Graphs illustrating the differences between the terms accurate and precise.
The vertical line corresponds to the true value. The vertical axis is a probability distribution
(or the number of measurements) and horizontal axis is the magnitude of
the quantity being measured.
Data Analysis and Errors • 117
4.3 UNCERTAINTY
4.4 RESOLUTION
4.5 TOLERANCE
long with a tolerance of 1 mm. This means that the metal needs to be
cut somewhere between 0.999 m and 1.001 m. Another frequently
occurring example is a resistor. One of the color bands on the resis-
tor corresponds to the tolerance. For example a 100 resistor with
a gold band (which indicates a tolerance of ± 5 %) should have a
resistance between 95 and 105 .
4.6 SENSITIVITY
The response time is the time taken for the measurement device
to respond and reach its final reading after a dramatic change in
input. Digital meters typically have a faster response time than an
analogue meter with a pointer because the analogue meter needs a
comparatively long time to physically move. Changes in the quan-
tity being measured which happen faster than the response time
of the measurement device will be averaged out and not seen. A
good demonstration of this would be using an alternating current
from a signal generator as the input to a center-zero analog meter.
Starting from a low frequency (with period much longer than the
response time) where the pointer is seen to oscillate the frequency
can be gradually increased until the period is much shorter than the
120 • Physics Lab Experiments
response time where the pointer will simply twitch slightly around a
reading of zero.
m y
N
i 1 yi
(4.2)
N
y
1
– y
2
s2 i (4.3)
N
i 1
Table 4.1: Example data for calculation of the mean and variance.
122 • Physics Lab Experiments
1
yi – y .
N
would be possible to find the variance from s 2
N i 1
y i – y is
N
However, this does not work in practice because i 1
y
1
s
2 2
i – 2 yy i y 2 (4.4)
N
i 1
N
x i2 2 yy i y 2
N –
i 1
N
N
(4.5)
y 2N (4.6)
y 2 – 2y 2
N
y2 – y2 (4.7)
where
N
yi2
y 2
i 1
N (4.8)
and
N
N
yi
y (4.9)
i 1
Figure 4.2 shows two number lines with values marked upon
them: one has a small spread and the other a much larger spread.
Both have the same mean, however more confidence could be
placed in the calculated mean being equal to (or very close) to the
true value if the data followed the pattern with the smaller spread:
this reinforces the point that the mean does not give us the whole
picture when analyzing results.
As an alternative to the number line, a set of measurements may
be represented by a histogram. Figure 4.4 shows a histogram for
the 10 measurements of the current flowing through a lamp given
Figure 4.4: Histogram of the measurements of the current flowing through a bulb given in
the text.
Figure 4.5: Histogram of the measurements of a quantity when the number of measure-
ments tends to infinity.
Data Analysis and Errors • 125
–
p y dy 1 (4.10)
y – Y p y dy
2
s2 (4.12)
–
In practical terms, using a small set of real data the best estimate of
standard error on a single measurement is given by the equation:
y
1 N
– y s
2
s i (4.13)
N –1 N –1
i 1
s s
sm (4.14)
N N –1
6
s 9.74 10.7 (4.15)
5
10.7
sm 4.4 (4.16)
6
x 250 4 (4.17)
error for this individual reading, but the error in a group of read-
ings is random because it is equally likely to be positive as negative.
The error will range from nothing to half a division (0.5 cm), and so
the average error will be about a quarter of a division (0.25 cm). The
standard error (corresponding to 1 standard deviation if we assume
a normal distribution) will be about a third of a division (0.33 cm).
Care must be taken to give a reasonable estimate of the error in
measurements like this.
1 y0 y2
P y
s 2p – y0
exp – 2 dy
2s
(4.18)
y–m
z (4.19)
s
The “z” value for data being within one standard deviation of the
mean is 1, since the difference y−µ = . This corresponds to a prob-
ability of 0.68 or 68 % that a given measurement is within 1 standard
Data Analysis and Errors • 129
Table 4.2: Probability that a given measurement is within ‘z’ standard deviations of the
mean.
deviation of the mean. Table 4.2 shows the probabilities for higher
numbers of standard deviations.
4.15.1 LINEAR
a b c (4.20)
da db dc (4.21)
b – b c – c
2
(4.25)
b –b c – c 2 b –b c – c
2 2
(4.26)
sb2 sc2 2 b – b c – c (4.27)
The final term depends on whether the errors in b and c are corre-
lated in some way. If there is no correlation then the size of the error
in b has no effect on the size of the error in c and visa versa. This
means that the term b – b c – c will average to zero. This reduces
the error in a2 to the form given in Equation 4.22.
a b x c y (4.28)
Following the method applied in the linear case above, take the
mean square values added in quadrature:
2 2 2
sa 2 sb 2 sc
a x b y c (4.31)
A proof of this follows the same pattern as the formal proof given in
the linear case.
132 • Physics Lab Experiments
Again, Equations 4.31 and 4.32 are only valid if the errors in b and
c are uncorrelated.
It can be demonstrated that the overall error in not the same for
correlated and uncorrelated errors by considering the measurement
of the diameter, d, of a wire and then using this to calculate the cross
sectional area, A, of the wire. Taking the equation for the area as:
pd 2
A (4.33)
4
pd d
A (4.35)
4
This is only half the error found in Equation 4.34 and does not give
a valid estimate of the error in A.
Data Analysis and Errors • 133
f f y 1, y 2 , y 3,, y n (4.37)
f – f0
2
sf2 i (4.38)
f tan x (4.39)
f 0 tan(89) 57 (4.40)
f 1 tan(89.5) 115 (4.41)
Now calculate:
f – f 0 57 – 115 3364
2 2
sf2 i (4.42)
sf 58 (4.43)
f 57–19
58
(4.44)
Notice that a different upper and lower limit has been found and it
is entirely appropriate to quote the error like this.
A line of best fit may be drawn by hand onto a graph. The aim
should be to draw a single smooth line without using repeated short
strokes of the pencil. There should be no sharp kinks in the graph
unless the data obviously indicates this. The line need not go directly
through any of the points and should have approximately the same
number of points sitting above the line as below. Sometimes a flex-
ible ruler is useful for drawing curved lines or a ruler with a slit in the
middle (available from http://www.bestfitlineruler.com) is useful for
straight lines. The gradient of a straight line should be taken using as
large a triangle as reasonably practical. For a curved line, the gradi-
ent may be found at individual points by taking the gradient of the
tangent to the curve at this point.
Computer programs can be used to draw a line of best fit (some-
times called a trendline) based on a least squares fitting (described
later in more detail) to a mathematical function such as y=mx+c or
y=a+bx+cx2.
Error bars should be added to graphs for all x and y values where
possible. They give a measure of the uncertainty in the position of
the point on the axes. If a number of repeats are taken for each point
then the error bar should extend above the point by a distance corre-
sponding to m and also below the point by a distance corresponding
to m. Horizontal error bars are more usually related to the resolution
136 • Physics Lab Experiments
y mx c (4.45)
V – Ir (4.46)
y c mx (4.47)
l
T 2p (4.48)
g
If the square of the time period is plotted on the y axis and the length
is plotted on the x axis, the gradient of the graph should equal 42/g
and the acceleration due to gravity may be measured. The y axis
intercept in this case should be zero. If a significant y axis intercept
is found then a systematic error should be suspected in the data.
Log graphs are usually used because of one, or both, of two main
benefits. Using a logarithmic scale on one or both axes of a graph
allows data which span a very large range to be plotted in a more
meaningful way. Measurements in the range 1 to 10 take up the
same proportion of the axis as measurements in the range 10 to 100
or 100 to 1000. If a log scale is used for the x axis data then logx
rather than x must be plotted on the x axis. It is most convenient to
plot a log graph using a computer although it may be done on paper
using special logarithmic graph paper (which has non-equal spacing
between the grid lines).
A second benefit is that by plotting both the x and y values on a
log scale a power law in the data can found, since:
y xn (4.51)
log y n log x (4.52)
If logy is plotted on the y axis and logx is plotted on the x axis then
the data should form a straight line with gradient n.
df x f x d x – f x
(4.53)
dx dx
Step 3 of the “Chart wizard” allows a chart title and axis titles to
be set and gives a choice which grid lines to display (along with a few
other less note-worthy options).
A trend line may be added to the graph by right clicking the
mouse on one of the data points and selecting “Add Trendline”:
see Figure 4.8. An appropriate trend line should be selected (often
Figure 4.8: Adding a polynomial trend line of order 2 to the data and ensuring the equation
is displayed on the chart.
With Excel 2010 highlight the data with the mouse (holding the
Ctrl button while using the mouse allows selection of non-contig-
uous regions). Next, select the Insert ribbon and then the “Scat-
ter” option as shown in Figure 4.10. In most cases the “XY Scatter”
option will needed as this will plot one set of points verses the other.
Data Analysis and Errors • 141
Figure 4.9: Adding X axis error bars. They have a fixed size of ± 0.25.
This will insert a chart onto the spreadsheet. If Excel has plot-
ted the axes the wrong way round these can swaped by selecting the
“Select Data” button from the “Chart Tools, Design” ribbon. This
opens a window which allows the problem to be solved: see Figure
4.11.
To add titles to the x and y axes select the “Axis Titles” button
from the “Chart Tools, Layout” ribbon as shown in Figure 4.12. The
chart title can be set in the same way.
Figure 4.13: Adding a linear trend line to the data and ensuring the equation is displayed on
the chart.
Figure 4.14: Adding X axis error bars. They have a fixed size of ± 1.
144 • Physics Lab Experiments
4.17.3 KALEIDAGRAPH
By default the fit parameters box is added to the chart. This gives
the values of the coefficients m0, m1, and m2 where y=(m0)+(m1)
x+(m2)x2. As the data plotted follows a simple y=x2 format, it is no
surprise to see that m0 0, m1 0 and m2=1: see Figure 4.19.
Error bars can be added by right clicking on the graph and
choosing “Error Bars...” from the popup menu. In the “Error Bar
Variables” window insert a cross in the “X Err” box to add X error
4.17.4 MATLAB
To draw a graph in MATLAB, first load the data points into vec-
tors within the software. A very simple graph can be drawn using the
following code which can be entered into the MATLAB editor and
saved as a .m file:
i=0:1:100;x=i;y=i.*i;plot(x,y);
x2=j;
y2=200*j;
hold on plot(x,y);
plot(x2,y2);
hold off
The line styles can be changed to dashed and the colors altered to
red and blue and axis labels added by typing:
hold on
plot(x,y,’-b’,’LineWidth’,2);
plot(x2,y2,’-r’,’LineWidth’,2);
xlabel(‘x title’)
ylabel(‘y title’)
Data Analysis and Errors • 149
Figure 4.22: A MATLAB graph showing two lines and axes labels.
hold off
The properties of the graphs such as the line style, line thickness,
colors, labels, and fonts can be set in the code as illustrated above,
but they can also be modified by selecting the “Tools” menu from
150 • Physics Lab Experiments
the graph window and choosing “Edit Plot”. Then right click on the
graph, the axes or the lines to bring up menus to allow the properties
to be modified for the individual graph.
MATLAB does not have a built in way of simply adding a best fit
line to a graph. One of the least squares method described later in
Section 4.18 must be employed. Usually the necessary computer
code is already written (often called a numerical recipe) and
the function must me called and the results it gives interpreted.
A detailed explanation of this is illustrated in Section 4.18.5.
Only vertical errors bars can be added in MATLAB. The code
needed for this is:
i=0:10:100; x=i; y=i. *i; e=0.1. *y;
errorbar(x,y,e);
y a bx (4.54)
a y – bx
x i – x y i – y
N
S xy i 0
(4.55)
b
S xx
2
xi – x
N
i 0 (4.56)
152 • Physics Lab Experiments
This is the method that Microsoft Excel uses to add a linear trend-
line to a graph. The R2 value given by Excel provides a measure of
the correlation between the data and the trendline. It does not, as
is often thought, provide a measure of the goodness of the fit of the
line to the data. The R value is actually the product moment correla-
tion coefficient where:
x i – x y i – y
N
S xy i 0
R
S xx S yy
N
i 0
xi – x
2
N
i 0
yi – y
2
(4.57)
N
xy –
i
x y
i
N
i 0
i
N
i 0
i
i 0 N
N 2 N 2
xi yi
N i 0 N i 0 (4.58)
x i2 – y i2 –
i 0 N i 0 N
Figure 4.24: An Excel trendline, showing the equation of the straightline and the R2 value
are given by Equations 84, 85, and 86.
This is precisely what Excel finds when it plots a trendline for this
data as shown in Figure 4.24.
N
eˆ i (4.60)
i 1 x i
which is
xˆ yˆ zˆ
x y z (4.61)
Data Analysis and Errors • 155
f x, y 8x 3 – 3x 2 xˆ 2 y yˆ (4.62)
Taking the starting point as (a,b), the point (a,b) must be found such
that a= a – —f(a, b) and bb – — f (a,b). is the step size, if this
is small enough then f(a,b) > f(a,b) and so by finding f(a,b) a smaller
value for the function f(x, y) has been found. This procedure can be
reproduced iteratively until, on a given iteration f (a, b) and f(a, b)
differ by less than a set tolerance. Often the function f(a, b) would be
the sum of the square differences between a set of experimental data
points and a set of data points derived from a equation of a line of best
fit. The steepest descent method is simple, easy to apply computation-
ally and each step is very fast. It is also quite efficient at finding the
minimum point if just one exists. If more than one minimum point
exists, or if there is a local minimum, this method may struggle or reach
different answers depending upon the starting point. The downside of
this algorithm is that the progress towards the minimum gets slower
and slower and often a very large number of steps are required. Also,
this only works where the function is differentiable. If this differen-
tiation needs to applied numerically, the algorithm can become quite
slow. In the case of fitting a function to experimental data the function
f(x,y) would be the sum of the square differences between the experi-
mental data and the data generated by the fitting function.
4.18.3 GAUSS-NEWTON
ri yi – f xi, b (4.63)
156 • Physics Lab Experiments
r1 r1
(b) (b)
b1 bM
J(b) (4.65)
rN rN
b (b) b (b)
1 M
–1
b k 1 b k – J b k J b k J bk r bk
T T
(4.66)
– e b2 – b1 e b2
2b2
–e – 2b1 e 2b2
J (b) J b1,b2 – e 5b2 – 5b1 e 5b2 (4.78)
– e 7b2 – 7b1 e 7b2
– e 9b2 – 9b1 e 9b2
and then:
–1.2214 –2.1985
–1.4918 –5.3706
J (b0) –2.7183 –24.4645 (4.84)
–4.0552 –51.0955
–6.0497 –98.0043
So:
–1 (4.85)
b1 b0 – J b0 J b0 J b0 r b0
T T
The process is then repeated to find b2, b3 ... bz where there are
z iterations. If z is large enough it will eventually converge to
(1.800,0.200).
The Gauss-Newton algorithm may converge slowly or not at all
if the initial guess, b0, is far from the minimum and it can only be
used to minimize the sum of the squares differences of a function.
However, it can converge very quickly and no second derivatives
need to be calculated.
So:
N
y – f x , b J
2
S (b ) i i i (4.91)
i 1
2
S (b ) y – f (x, b) – J (4.92)
which can be solved for d by simply taking the inverse of JTJ and
applying it to both sides.
Levenberg modified this slightly to introduce the damping term
mentioned earlier:
J T
J ldiag J T J d J T y – f (x, b) (4.95)
Figure 4.25: The different triangles involved in a two parameter version of the Nelder-Mead
function minimization algorithm.
W and the simplex is shrunk to BMS (see Figure 4.25d). This process
repeats until a predefined convergence criterion is reached such as
S becoming less than a certain value, failing to change between suc-
cessive iterations or a maximum number of iterations is reached.
MATLAB implements this method using its fminsearch routine
in the optimization toolbox. The MATLAB function fminsearch can
be modified to calculate and return the covariance matrix i.e., the
errors on the fit parameters using the method outlined in the appen-
dix of [16]. The declaration at the top is modified to include an addi-
tional parameter N which is passed in the function call which gives
the total number of data points and also to return the error vector.
This additional code was added after line 382 in the code for
fminsearch and is based on the appendix of [16].
%calculate covariance matrix
%qmatrix
temp=size(v);
numberofparams=temp(1,1);
for i=1:numberofparams
qmatrix(i,:)=v(i,:)-v(i,1);
end
qmatrix(:,1)=[];
%bmatrix
for i=1:numberofparams
for j=1:numberofparams
if i==j
yi=fv(1,i+1);
y0=fv(1,1);
162 • Physics Lab Experiments
p0i=(v(:,1)+v(:,i+1))./2;
y0i=funfcn(p0i,varargin:);
bmatrix(i,i)=2*(yi+y0-2*y0i);
a(i)=2*y0i-(yi+3*y0)/2;
else
pij=(v(:,i+1)+v(:,j+1))./2;
yij=funfcn(pij,varargin:);
yi=fv(1,i+1);
y0=fv(1,1);
p0i=(v(:,1)+v(:,i+1))./2;
p0j=(v(:,1)+v(:,j+1))./2;
y0i=funfcn(p0i,varargin:);
y0j=funfcn(p0j,varargin:);
bmatrix(i,j)=2*(yij+y0-y0i-y0j);
a(i)=2*y0i-(yi+3*y0)/2;
end
end
end
covarmatrix=qmatrix*inv(bmatrix)*qmatrix’;
ymin=y0-a*inv(bmatrix)*a’;
normalizefactor=2*ymin/(N-numberofparams);
covarmatrix=covarmatrix.*normalizefactor;
%getstandarddeviation
for i=1:numberofparams
errorvector(i)=sqrt(covarmatrix(i,i));
end
expData = interp1(Data(:,1),Data(:,2),x,’spline’);
%interpolate data to values at simulation points
given by x
params(1)=40; %starting parameter
params(2)=5; %starting parameter
numberofdatapoints=length(x);
disp([‘initialguess=’num2str(params)])
estimates=fminsearch(@fit,params,numberofdatapoints,
options,x,expData);
164 • Physics Lab Experiments
disp(‘Estimates’)
disp(num2str(estimates))
fitprm = estimates;
vN – p (4.97)
Figure 4.26: A graph showing the maximum and minimum gradients consistent with the
error bars.
1 N ri2
gradient i 1
2
(4.98)
N xi – x N –2
2
i 1
and
1 x2
N ri2
N i 1
2
intercept 2 (4.99)
i 1
N xi – x N –2
line graphs. This package is not always loaded into Excel. In Excel
2003 click on “Tools” and then choose “Ad Ins”. Then select the
“Analysis ToolPak” checkbox. The option to select “Data Analysis”
is available from the “Tools” menu. Find and select the “Regres-
sion” option. This brings up the “Regression” window shown in
Figure 4.27.
Once the “Input X Range” and the “Input Y Range” have been
set and the “Output Range” selected, the calculations begin. The
output looks similar to Figure 4.28.
In Figure 4.28, the intercept and gradient are given in cells B27
and B28. Cells C27 and C28 give the standard error on these values.
In Excel 2010 click on “File” then “Options”. Next, choose
“Add-Ins” from the menu at the left. Find and click once on “Analy-
sis ToolPak” in the list of Add-Ins. Press the “Go...” button. Select
the “Analysis ToolPak” checkbox. Click on the “Data” ribbon, a new
option has appeared on the righthand side called “Data Analysis”.
Press this, and then find “Regression” from the list. The “Regres-
sion” window and the output table is the same as for Excel 2003.
For more complex graphs the uncertainty in the fitting param-
eters can be found by fitting the data to an equation using a least
squares fitting method. This is illustrated in Section 4.18.5.
168 • Physics Lab Experiments
Figure 4.29: Entering the data ranges into Fourier Analysis window.
N
f (4.100)
t
2pi
H (i) 0.5 1 – cos
N – 1
(4.101)
Figure 4.32: Fast Fourier transform of a sine wave with frequency 10 Hz. Left panel shows
leakage with no windowing function present. Right panel shows leakage reduced with a
windowing function.
CHAPTER
5
PRESENTING RESULTS
5.2 REPORTS
As well as keeping a neat lab book showing all the results and
analysis of each experiment, students will often be required to write
up two or three experiments formally as though they were scien-
tific articles. This is a vital part of becoming a good Physicist, as
once a cutting edge experiment is completed, the results analyzed
and a new discovery has potentially been made the results need to
be communicated with other Physicists. This has traditionally been
achieved by writing the experiment up as a scientific paper and then
submitting it to a scientific journal. The journal then ask a num-
ber of other Physicists (the identities of whom are unknown the the
authors) working in similar fields at other universities or research
laboratories to be a referee for the paper and to make critical com-
ments upon the work. Referee comments are passed back to the
authors for their comment or so that they make modifications to
their paper. Based upon these comments the journal editors will
make a decision to accept or reject the paper for publication.
The aim of writing formal reports is to develop skills in writing
papers: these will be vital in a future career in academia and almost
certainly useful for a future career in any other field. They provide a
good opportunity to gain a deeper understanding of an experiment
beyond that achieved in a one or two day laboratory session.
This section gives an overview of the typical structure of a formal
written report and provides some advice on how to include refer-
ences and write a report on a computer.
5.2.1 SOFTWARE
5.2.2 STRUCTURE
Theory
A theory section should explain the physics behind the problem
and show how the experiment carried out and the quantities mea-
sured can be related to the problem.
Experimental Methods
A description of the apparatus used, including a diagram. This
should include information on how the apparatus was built or
connected together and any special features of the equipment or
particular problems which had to be overcome. A description of
the measurements taken and any special techniques which were
required.
Results and Discussion
Raw data, especially large tables of numerical data should be
placed in an appendix at the end of the report. This section should
show how the data was analyzed to give a final result and contain any
important intermediate results. A small number of graphs or sum-
mary tables are appropriate, although most extensive data sets can
be put into an appendix. It is vital to include an estimate of uncer-
tainties and to consider all sources of experimental error in the final
results.
Conclusion
The results of the experiments, together with appropriate
uncertainties, should be summarized and compared with other
values. The results should be set in context and their significance
noted. Suggestions for possible future work can be made where
appropriate.
References
Sources of information such as other papers, articles, and books
should be listed here.
Appendix
Detailed tables, graphs, derivations, or computer programs
should be presented in the appendix.
In Microsoft Word or OpenOffice Writer, the headings
may be formatted to have a specific style: perhaps a larger font
and in bold. In LaTeX, simply use the “section” commands as
follows:
Presenting Results • 179
\section{Section Title}
Some text.
\subsection{Subsection Title}
Some more text.
\begin{center}
\begin{figure}
\includegraphics [width=4in] {folder/filname.png}
182 • Physics Lab Experiments
\begin{table}
\begin{tabular}{ll}
\hline x & y=x$^2$ \\
\hline 1 & 1 \\
2 & 4 \\
3 & 9 \\
4 & 16 \\
\hline
\end{tabular}\label{tablename}
\end{table}
These commands to insert images and tables mean that LaTeX will
determine their position in the document to best suit the flow of test
and figures across pages. If insertion at a fixed point in the text of
preferred, regardless of how the text flows around, omit the begin/
end figure, or begin/end table commands.
To create bullet point lists in Latex, use the “itemize” structure:
\begin{itemize}
\item The first item
\item The second item
\item The third item
\end{itemize}
Presenting Results • 183
\begin{enumerate}
\item The first item
\item The second item
\item The third item
\end{enumerate}
5.2.4 REFERENCING
When writing a report always make clear which parts are written
by the author and which parts are taken from, based upon, derived
from or copied from other sources. For instance it is not possible to
state “It is well known that the charge on an electron is 1.6×10−19 C”.
Statements like this should refer to a source which backs it up, for
example: “The data in the present report give a value of the charge
on an electron e=2.4±0.2×10−19 C, which compares well with the
study of Price and Arthur [18] who obtained e=2.2±0.1×10−19 C”.
Always check references carefully, especially the volume and
page numbers and the spelling of authors’ names/any unfamiliar
words in the title.
There are two accepted styles for labeling references: the
numerical system and the Harvard system.
Referencing: The Numerical System
The numerical system is generally in wider use in the sciences,
especially Physics. It is the method adopted by most journals for
the publication of original research. References are listed in a num-
bered list at the end of the paper. In the main text of the paper, the
contribution of a reference in indicated with the relevant number
either as a superscript34 or in square brackets [34]. It is good prac-
tice, although not essential, for the references to be numbered in the
order in which they are referred to in the text.
Various pieces of software can help with this and make the job of
ordering numerical references less tedious. Microsoft Word has the
option for endnotes. To insert a new reference choose “Insert, Ref-
erence, Footnote”, then choose the “Endnotes” bullet point. This
184 • Physics Lab Experiments
will give a list at the end of the document, rather than at the end of
each page (if the “Footnote” option is chosen).
In LaTeX, a separate bibliography file can be created with the
file extension “.bib”. In the “.tex” document, the bibliography file is
linked by including the following lines of code:
\bibliographystyle{unsrt}
\bibliography{file_name}
This “unsrt” style produces a list in which the references are listed
in the order in which they are referenced in the text. Changing the
style to “plain” gives a list in which the references are alphabetical by
first author surname. The “.bib” file is a text file with a list of refer-
ences. For example, the entry for one of the references used in this
publication is arranged as follows:
@Article{liquidoxygen,
title = Making Liquid Oxygen,
author = French, M. M. J. and Hibbert, Michael,
journal = Phys. Ed.,
volume = 45,
pages = 221,
year = 2010,
}
@MISC{winedt,
note = http://www.winedt.com
}
\bibliographystyle{unsrt}
\bibliography{references}
5.2.5 EQUATIONS
Microsoft Word
In Microsoft Word 97 - 2003, equations can be inserted using
the “Insert, Object, Microsoft Equation 3”. This brings up a tool
bar, from which Greek characters, mathematical expressions, and
layouts can be selected and inserted into an equation object. Figure
5.4 shows the tool bar being used to select a fraction arrangement.
\begin{equation}
\frac{a}{b}
\end{equation}
which produces
a
(5.1)
b
\begin{equation}
\left[ \left( x \right) \right]
\end{equation}
which produces
x (5.2)
\begin{equation}
\int_0^\pi A\times\cos^2 \theta \partial\theta = A\
times\int_0^\pi \left(\frac{e^{\imath\theta}+e^{-\
imath\theta}}{2}\right)^2
\partial\theta
\end{equation}
which produces
2
p p eiq e–iq
A cos qq A q
2
(5.3)
0 0 2
5.3 POSTERS
5.3.1 SOFTWARE
B n A n A n – 1 (5.4)
Tolerances for the paper size are ±1.5 mm for lengths up to 150mm,
±2mm for lengths between 150 and 600 mm and ± 3mm for any
dimension above 600 mm.
In Microsoft PowerPoint 2003, the page size is changed by select-
ing “Page Size” from the “File” menu and changing the “Width” and
“Height” settings. In OpenOffice Impress, the page size is changed
by selecting “Page” from the “Format” menu and changing the
“Width” and “Height” settings.
Finally, it should be ensured that images are of sufficient resolu-
tion that they will not appear “blocky” and distorted when they are
printed out on a poster. Ideally there should be at least 150 dots per
inch and probably less than 300 dots per inch to avoid excessively
large file sizes. So if a graph is required to appear about 6 inches
wide (15 cm) on the final poster, it should be between 900 and 1800
pixels wide.
5.4 PRESENTATIONS
Too many people stand behind a table holding notes and appear to
the audience to be “hiding” or seeking protection behind furniture.
Slides should not be blocks of text, perhaps a maximum of 30
words per slide as a general guide, nor should they have more than
one graph per slide. Elaborate slide transitions and flashing text
should always be avoided. Make sure text is large enough to read
(this includes graph axis labels) and the color of the text and back-
ground aid rather than hinder its legibility. For example, black on
white is easy to read, black on dark blue or yellow on white is very
hard to read. For ease of reading fonts should not be so complex
and fussy that the audience struggle to read them, although recent
research suggests [21] that if serif fonts such as Times New Roman
or Monotype Corsiva are used, the audience will remember more of
the presentation.
Start with a title slide. This should be followed with a slide giving
a brief outline of the talk. A good way to plan a short presentation
is to complete the outline slide and then turn each point from the
outline into a slide title. Finally, include an acknowledgments and
references slide unless this information is included elsewhere. As a
rough guide around one slide per minute is needed, so a 10 minute
presentation should have 10 slides in addition to slides for the title,
outline, references, and acknowledgments.
Figure 5.11 shows an example set of slides for a short presenta-
tion. The accompanying notes are as follows:
1. Originally a Sci-Fi creation, which was first proposed scien-
tifically by Carl Sagan in 1961.
2. It hasn’t been proven that humans can live in enclosed envi-
ronments for long periods of time i.e. space stations.
3. Correct Temperature - not too hot (velocity of particles
is greater than escape velocity of planet) and not too cold
(can’t heat it enough); Ability to retain an atmosphere - size
of planet - enough gravity; Location of planet - solar winds
remove atmosphere; magnetic field existence/strength;
Within easy traveling distance of earth - transport of equip-
ment, other resources if needed, people; Availability of
resources - things needed to terraform - lots of water and
CO2; Mars - close, once had water flowing on surface,
198 • Physics Lab Experiments
5.4.1 SOFTWARE
6
COMMON
EXPERIMENTS
This experiment is often the first which you will undertake in the
undergraduate laboratory. The aim of this experiment is to measure
the acceleration due to gravity to a precision of one part in a thou-
sand. It also has the objective to teach you to understand errors. You
will learn how to treat random error in measurements of the time
period of the pendulum and how to combine errors from different
sources to give an estimate of the error in the final value of the accel-
eration due to gravity.
You will use a rigid pendulum pivoted about a knife edge. For
small oscillations the period, T, of the pendulum is given by:
I
T 2p (6.1)
mgh
4p 2 I
g (6.2)
mhT 2
Common Experiments • 203
6
1. The sum of the six readings xi
i 1
You will find that you are probably some way off achieving an error
of 1 part in 2000 or 0.0005. You could, of course, continue to repeat
the 5 period method enough times as the standard error in one
period decreases as the number of repeats, n, increases. However
due to the 1 n dependence a very large number of repeats will be
required.
An alternative approach is to reduce m by reducing rather
than increasing n. You can do this by increasing the number of peri-
ods measured from 5: this will reduce the random error, but also the
systematic error in the period T.
If your reaction time is around ± 0.1s then to achieve 1 part in
2000, you will need to measure for around 200 seconds. To ensure
you reach the target of 1 part in 2000, try to measure a time of
around 300 – 400 seconds. However, it will be very easy to lose
count of the number of periods in such a long time so try the fol-
lowing method.
Using the stopwatch time the period of an integer number of
oscillations (it doesn’t matter exactly how many) which take around
60 seconds in total. The actual number of periods can be found by
using your previous value of T and knowing there was an integer
number of periods. You should now be able to fill in the table below:
You should repeat the procedure, but timing for around 300 -
400 seconds and fill in the table again:
You should now have a value for T with an error of less than 0.0005.
Finding h and m
These should be relatively straightforward. The mass m can
be determined by placing the rigid pendulum onto a balance.
The calibration of the balance should be checked using standard
masses and a number of repeats taken. The mass and the error in
the mass should be recorded. If the pendulum has mass around
1kg, then a measurement to the nearest 10th of a gram will give an
error of one part in 10,000. To measure the height h the center of
mass of the pendulum must be found by balancing the rigid pen-
dulum on a knife edge. The pendulum rod can then be marked at
the balance point. The measurement h should be taken from the
knife edge at the top of the pendulum rod to the mark you have
just made. Over a distance of around 500 mm, you will need to
measure to the nearest 0.5 mm to achieve an error of less than 1
part in 1000.
Finding the Moment of Inertia
The final measurement to be taken is the moment of inertia,
I, of the pendulum about the knife edge. With a compound object
like the rigid pendulum, the moment of inertia can be most eas-
ily found by using the “parallel axis” theorem. This states that the
moment of inertia, I, of an object about a given axis is related to
its moment of inertia, Ip, about a parallel axis through the center
of mass by:
I I p mh 2 (6.3)
206 • Physics Lab Experiments
where m and h have the same definitions we have already used. See
Figure 6.1.
You should find that Ip is only about a tenth of the magnitude of
mh2 and so only needs to be measured to a precision of about 1 part
in 100.
Ip, the moment of inertia of the pendulum about the axis
through its center of mass, can be found by comparing torsional
oscillations of the pendulum with those of a regularly shaped rect-
angular bar:
2
I p Tp
I b Tb (6.4)
where Ib is the moment of inertia of the bar about its center of mass,
Tp is the time period of the torsional oscillation for the pendulum
and Tb is the same for the bar.
Common Experiments • 207
12
a b
mb 2 2
Ib (6.5)
where a and b are the length and width of the bar and mb is the mass
of the bar. Since you are looking for Ip to 1 part in 100, mb is needed
to 1 part in 100 and a and b are needed to 1 part in 200.
You should suspended the pendulum from a wire about its cen-
ter of mass. The pendulum can then be displaced so that it under-
goes a torsional oscillation. This means that it twists from side to
side in the horizontal plane, with the center mass staying directly
beneath the point of suspension. The time period of the torsional
oscillation should be measured following a similar procedure to that
used previously. This is repeated for the bar. Since you are aiming
for an error in Ip of less than 1 part in 100, the time periods must
have an error of less than 1 part in 200.
To determine the errors in Ib and then Ip do the following. Write
a +b2 = X then:
2
X 2 a a 2 b b
2 2 2
(6.6)
and calculate X. You may like to try proving this using the equa-
tions given in Section 4.15 of this book. The error in Ib can now be
found from:
2 2 2
I b mb X
I m X (6.7)
b b
2 2
I p I b T p
2 2
Tb
I p I b 4 T p 4 Tb (6.8)
208 • Physics Lab Experiments
Since:
4p 2 I p mh 2
g (6.9)
mhT 2
(6.10)
so:
2 2 2
g h T
g h 2 T (6.11)
If Ip/Ip > 0.01 then a more complex equation must be derived for
the error in g.
Discussion
Finally, it is worth identifying the largest sources of error in the
experiment and writing a discussion focusing on whether the error
in g met the target of being less than 1 part in 1000. You should
consider how the value of g could be improved either by modifica-
tion of the experimental technique or by redesigning the method
completely. Systematic errors also need to be considered: is your
location on earth or height above sea level likely to influence the
result significantly: by how much?
nx
t t0 (6.12)
c
You may need to adjust the position of the ends of the fiber optic
cables to optimise the amplitude of the received signal, especially
with the longer cables where the detected signal is weaker. The cur-
sor controls on the oscilloscope can be used to measure the time
delay. Choose the most appropriate points to measure between and
try to justify your choice.
You should consider how reproducible your result is and whether
or not it should be repeated. Finally, repeat the measurement with a
range of optical fibers up to 100 m in length.
To analyze these results, plot a graph (including error bars for
both quantities) of the time delay verses the length of cable. Use the
graph to determine the delay time due to the electronics and the
ratio n/c. The refractive index, n, will be given to you, allowing you to
calculate c. Using multiple lines on the graph, estimate the error in
n/c and use this to express the error in c. How many standard devia-
tions is your result away from the defined value? Consider the effect
of systematic errors: where might they arise, what could be done to
minimize them and would they lead to an increased or decreased
value for c?
Foucault’s Method
A beam of laser light is focused using Lens 1 and then passed
through a beam splitter. Part of the light enters the microscope and
part makes its way to the rotating mirror. When light is incident on
the rotating mirror it is reflected towards a fixed, spherical mirror.
The fixed mirror always reflects the light back along the same path,
returning it to the rotating mirror and then the beam splitter. Some
will be directed towards the microscope where it can be viewed in
conjunction with the light which traveled the shorter route. A dia-
gram of the apparatus is shown in Figure 6.2.
First, lets consider the point of reflection on the fixed mirror as
the angle of the rotating mirror changes. Figure 6.3 shows the path of
a light ray which is reflected from the rotating mirror at angle 2. At
a slightly later time, the rotating mirror has moved and now the light
ray reflects at angle 2(+). This means the two light rays reflect
from the fixed mirror at different points. If the distance between the
rotating and fixed mirror is D then, the distance between the two
points of reflection on the fixed mirror is given using the general
equation for the perimeter of a circle r:
S 2D dq (6.13)
Figure 6.3: Reflection from the fixed mirror as the rotating mirror changes angle.
212 • Physics Lab Experiments
v X
dp S S (6.14)
u YD
where X is the distance between Lens 1 and Lens 2 minus the focal
length of Lens 1 and Y is the distance between Lens 2 and the rotat-
ing mirror. Substituting from Equation 6.13 gives:
2DX dq
dp (6.15)
YD
Common Experiments • 213
2D w
dq (6.17)
c
4XD 2w
dp (6.18)
cY D
4XD2w
c (6.19)
d p Y D
w
m (6.20)
dp
Common Experiments • 215
4XD 2
c (6.21)
mW
Abstract
Electrons emitted from an indirectly heated cathode were accel-
erated through known electric and magnetic fields and by varying
the solenoid current, focused into a narrow beam that then hit a
fluorescent screen. When the beam was focused on the screen the
electrons performed an integer number of helix rotations. For a
series of anode–cathode potential differences, measurements of the
solenoid current and number of helix rotations were taken. After
careful analysis of the experimental data acquired, the charge:mass
ratio, [e/m], of the electron was calculated to be 2.99 × 1011 Ckg−1
±0.55 × 1011 Ckg−1.
Introduction
The calculation of the charge:mass ratio of the electron was
important in the early investigation and discovery of the structure
of the atom and in proving the existence of units of matter smaller
than the atom.
J.J. Thompson first carried out an experiment similar to this
in 1897. He was attempting to disprove the widely believed idea
that the hydrogen atom was the lightest particle in existence. His
experiment, in which he measured the deflection of charged par-
ticles in a magnetic field perpendicular to the direction of motion
of the particles, allowed him to show that the electron was nearly
2000 times lighter than the hydrogen ion. This experiment is based
on an adaptation of Thompson’s, first carried out by H. Busch in
1922.
e
This experiment aims to calculate and to test Busch’s adapta-
m
tion of Thompson’s experiment.
Theory
The equation F=Ee + Bev gives the force (measured in New-
tons) on an electron of charge e (measured in Coulombs) and mass
m (measured in kg), moving with velocity v (measured in m/s), in
a magnetic field B (measured in Tesla) and electric field E (mea-
sured in volts/meter). Since the acceleration due to this force is
Common Experiments • 217
Bev t m Be
w
vt m
v t mv t
r
w Be
L LBe m
n
Tv 2pm 2eV
Left: Figure (6.5) Electron paths when both have an integer number of helical rotations pro-
ducing a focused spot. Right: Figure (6.6) Electron paths when both do not have an integer
number of helical rotations producing an unfocused spot.
Common Experiments • 219
So:
e 8p 2n 2V
2 2
m LB
e
Where is the charge:mass ratio of the electron, n is the number
m
of helix rotations, V is the anode-cathode potential difference, L is
the anode-screen distance and B is the magnetic field.
The field inside an infinitely long solenoid is given by B=µ0NI
where µ0 is the permeability of free space, N is the number of turns
per meter and I is the current flowing in the solenoid. A correction
factor of K is needed since the solenoid is only of a finite length.
Where
X A D Z U A D W
K ln ln
2D A Y 2 D A V
taken, using an ammeter, for points when the spot was in focus
(i.e., integer values of n). As the current was increased from zero,
the first time the spot was in focus was assumed to be n = 1, with n
increasing by one each time the spot was in focus.
The solenoid was measured to find its dimensions to get a value
e
for K and to calculate .
m
After B had been calculated for each current, a graph of B2
against n2 was plotted so as the gradient was
e L2
m 8p 2V
Figure (6.8) shows the center of the apparatus, including the anode,
cathode, and the fluorescent screen. This is surrounded by the sole-
noid.
Results
The value of K was calculated to be 0.87. L, the anode-screen
distance was measured to be 0.051 m, the total length of the sole-
noid, L0, was measured to be 0.13 m and the total number of turns
on the solenoid was 1500. This gave N as 11500 turns/meter.
e
The gradients of the graphs gave for each (electric) potential
m
difference as:
Discussion
After plotting the graphs for each of the (electric) potential dif-
ferences, their gradients were found not to agree with each other
exactly, but to vary slightly. An average value for the gradient was
taken as the final answer that was around 1.7 times the accepted
e
This gives an average value of as 2.99 ×1011 Ckg−1.
m
222 • Physics Lab Experiments
value. This means that the gradients of the graphs I plotted from my
experiment were a little too large.
The error bars on the graph show the errors with respect to the
uncertainty in the solenoid current. The error range was obtained by
using the extremes of these error bars to calculate the gradient and
e
then . They do not show, what is probably the most significant
m
source of error in this experiment, the error with respect to deciding
the exact time when the spot was focused on the screen. This was a
visual measurement which is prone to error, which cannot be mea-
sured accurately.
Even though the gradients are a little large, the line of best fit
still passes very close, but slightly to the positive side of the origin
in all cases. This suggests that there could be some systematic error
involved in these readings - which has a greater effect at larger val-
ues of n. Supporting this idea is the fact that the values of B should
be multiples of the value for n = 1. This does not quite fit the experi-
mental data. The values fall short of what this predicts them to be–
and increasingly more so as n increases.
One possible cause of this could be the partial focusing of the
beam of electrons before it emerges from the anode. To eliminate
Common Experiments • 223
this the deflection of the spot when the electrons acquire a large
transverse component of velocity in a known part of their path
can be measured. This transverse component may be produced
by applying a potential difference between one pair of deflector
plates. The spot comes into focus when the electrons have trav-
eled an integer number of turns in their helical path from the
center of the deflector plates to the screen. The distance L, in
the earlier equation, can be substituted for the distance between
e
the deflector plates and the screen, and calculated using the
m
same methods.
Also, using an AC (electric) potential difference between the
deflector plates could make it easier to see when the spot is in focus
since the spot is then drawn out into a line and a collapse of the line
into a spot indicates an integer number of helical rotations. Since
the points lie so close to the lines of best fit, there is little evidence
to suggest any significant random error in this experiment. As well
as the extension involving applying AC (electric) potential differ-
ences to the deflector plates, a larger CRT could be used - as this
would make it easier to see the spot and when it’s in focus/not in
focus.
Conclusions
From my experiment the charge:mass ratio of the electron is
2.99 ±0.55 ×1011 Ckg−1. Considering the inaccuracy in the visual
“measurement” of when the spot is in focus, I consider my answer to
be consistent with others” results, even though the accepted value of
1.76 ×1011 Ckg−1 is outside my error range, it is of the same order of
magnitude and only a factor of 1.4 out. I consider Busch’s adaptation
of Thompson’s experiment to valid - especially if the above improve-
ment of applying an AC (electric) potential difference to the deflec-
tor plates behaves as predicted.
References
[1] Haliday, Resnick, Krane, “Physics”, Fourth Edition, 1992,
pp738 – 740
224 • Physics Lab Experiments
Appendix
Common Experiments • 225
Abstract
Some of the properties of radioactive count rate and the
Geiger Muller tube were investigated. The binomial distribution
for independent events was approximated by the Normal and Pois-
son distribution. The Poisson was found to be a better approxima-
tion since the probability of decay was small. The dead time of
the Geiger Muller tube was found to be 495 ± 73 µs. The inverse
square law for the count rate was tested with a gamma and a beta
source, the results showed an almost perfect fit to theory. Gamma
rays were found, as expected, to be more penetrating than beta
particles. Exponential attenuation was also investigated with lead
plates. The theoretical equation M = M0 e−µt was found to fit the
experimental data quite well, but a different equation, M M 0 e – kt ,
was found to provide a much improved fit, especially where the
standard thickness lead plates were not combined to provide addi-
tional thicknesses to test.
Introduction
History of radioactivity
Antoine Henri Becquerel (1852 - 1908)
In March of 1896, a French physicist, Antoine Henri Becquerel,
placed some wrapped photographic plates into a draw along with
some crystals containing uranium. Much to his surprise, when he
returned to the draw sometime later he discovered the plates were
exposed. He concluded that the uranium crystals must have been
emitting something, since the exposure had occurred without the
presence of the usual initiating energy source - sunlight. Becquerel,
though, did not pursue his discovery of radioactivity any further.
This was left to the Curies.
226 • Physics Lab Experiments
Theory
Binomial Distribution
Radioactive decay is a random process, so the probability of any
particular atom decaying is independent of whether any other par-
ticular atom has decayed. This is characteristic of the Binomial dis-
tribution.
P r n Cr pr q n–r
F t dt F t . 1 – ldt : (6.22)
228 • Physics Lab Experiments
dF
F t d t ... F t . 1 – ld t : (6.23)
dt
As t 0, (6.23) becomes:
dF
– lF t : (6.24)
df
F t Ce – lt : (6.25)
Since F(0)=1 then C=1 i.e., the particle had not decayed at t=0.
Equation (6.25) can be multiplied by N0, the total number of atoms
present to give:
N N0 e – lt
: (6.26)
1
e
– r – np 2 npq
2 2
P r n Cr pr q n–r
2pnpq
np
r
P r n Cr pr q n–r e –np
r!
ytheoretical – yexperimental
2
x
2
ytheoretical
x2
x v2
v
1 1 N 2 N1 – N 3
t – –
N3 N 32 N 1N 2N 3
Experimental Methods
Sources Used
All these sources have very long half lives so as the count rate
will be effectively constant over the whole period of the experiment.
Common Experiments • 233
Diagram of Apparatus
Background Count
The background count (i.e., the count rate with no source
present) was measured and subtracted from all subsequent mea-
surements. The count rate was found to be 0.48727 counts sec−1.
Dead Time Measurement
In order to calculate the dead time, t, values of N1, N2 and N3
were obtained. A value for N1 was obtained by performing a series of
measurements of the count over a period of 10 s with just one 204Tl
source present. A value for N2 was obtained by the same method with
the other 204Tl source present. Finally, a value for N3 was obtained by
the same method with both the sources present. Putting these values
into the equation found in the earlier section called Dead Time gives
a value for the dead time. Later, we will see the results and final
value for dead time.
Dead Time and the Background Count Rate
The background count was adjusted for the dead time. The
adjusted background count rate was now found to be 0.48739 counts
sec−1. This new value was subtracted from all subsequent results.
The Inverse Square Law
Since the particles emitted by a radioactive source are random
in both time and direction, the count rate should fall with the square
of the distance between the source and the Geiger Muller tube.
The experiment was performed by placing the source on each of the
shelves in the chamber and then taking a reading of the count rate over
a 10 second period. The two sources tested were 90Sr and 22Na. The
intershelf distance was constant and found to be 12.8 mm. Graphs of
1 M , where M is the count, against distance could then be plotted.
Using the intershelf distance method meant that there was no
need to know the absolute distance between the source and the
Geiger Muller tube window. In order to plot the results the distance
of the first shelf was arbitrarily chosen to be 1m. This method still
allowed for comparison of the results and for the inverse square law
to be tested.
After the points were plotted and a line of best fit was added,
Excel was used to give the equation of the line and also a value of
236 • Physics Lab Experiments
R2. From the equation of the line the x-axis intercept was also calcu-
lated. The results are given in the Inverse Square Law section.
Exponential Attenuation
The absorption of particles by a certain material can be inves-
tigated by placing the source at a constant distance from the Gei-
ger Muller tube and inserting various thicknesses of the material
in between. Alpha and Beta particles lose energy progressively
and have a definite range in a certain material but gamma rays are
expected to be absorbed exponentially.
In order to test whether this was the case, the 22Na source and
various thicknesses of lead were used. These plates were inserted in
various combinations in between the source and the Geiger Muller
tube and the count rate measured. The results are provided in the
Appendix.
In M = In M 0 mt
y c mx
Common Experiments • 237
it can be seen that if a graph of lnM against t is plotted with the inter-
cept fixed at lnM0, the gradient of this graph would be the constant
−µ . Graphs are shown in the Exponential Attenuation section.
Calibration Graph
The threshold (electric) potential difference (A) can be seen to
be 650 V and the start of the plateau (B) is at 750 V. For the rest of
the experiment the Geiger Muller Tube was run with a EMF supply
of 750 V.
The bars show the experimental results and the line shows the
theoretical results. The mean is 109, the variance is 10.2, the prob(2)
is 0.47 and the 2 value is 73.1.
Graph of the Poisson Approximation
Again the bars show the experimental results and the line shows
the theoretical results. In order to sufficiently reduce the mean
count rate a piece of aluminium 0.032 inches thick was found to be
necessary. The mean is 3.5, the variance is 3.5, the prob(2) is 0.61
and the 2 value is 7.29.
The Nx values shown above are already corrected for the back-
ground count and then the dead time was calculated. The mean,
standard deviation and standard error on the mean were calculated
240 • Physics Lab Experiments
- see Lyons (1991). The mean is 510 µs and the standard error on
the mean is 70 µs.
Beta Source
With the beta source, the x-axis intercept was found to be at
-0.19 and the R2 value was found to be 0.9977.
Common Experiments • 241
Gamma Source
Exponential Attenuation
Figure 6.18
242 • Physics Lab Experiments
Figure 6.19
With the gamma source, the x-axis intercept was found to be at -0.52
and the R2 value was found to be 0.998.
Figure 6.20
Common Experiments • 243
Figure 6.21
Figure 6.22
m s
r Ma
244 • Physics Lab Experiments
Discussion
Conclusions
The start of the plateau for the Geiger Muller tube was found to
be 750 V. Both the Normal and Poisson distribution are valid approx-
imations to the Binomial in this case, although the Poisson provides
the closer fit between experimental data and theory. The dead time
of the Geiger Muller tube was calculated to be 490±70 µs. The count
rate for both the gamma and beta sources was found to accurately
obey the inverse square law. The data gathered while investigating
exponential attenuation supports the equation, M=M0 e−µt given in
Kaplan (1956), but the modified equation, M M 0e k t fits the data
more accurately. The cross section of a scattering center was found
to be 33.32 barns.
References
Appendix
Dead Time
If N counts are recorded in 1 second, and the tube is dead for a
time t after each count then in 1 second the tube will be dead for Nt
seconds. If an ideal counter has M counts in 1 seconds then
N M MNt
N M (1 Nt )
M 3 M1 M 2
N3 N1 N2
1 – N 2t 1 – N 1t 1 – N 2t
Simple rearranging yields a quadratic in t:
0 = N 1 N 2 N 3 t 2 2 N 1 N 2t (N 2 N 1 N 3)
N 2 N 1 N 3 2 N 1 N 2t
Abstract
The latent heat of vaporization of liquid nitrogen was calcu-
lated by measuring the boil-off from an insulated dewar flask when
a known amount of heat was supplied. It was found to be signifi-
cantly higher than the accepted value of 199 kJkg−1 at 249 kJkg−1.
The high value can be accounted for because the volume measuring
was stopped as soon as the heater was switched off, and not contin-
ued as the heater cooled. The specific heat capacity of a range of
materials was then calculated over two temperature ranges, 293 K
252 • Physics Lab Experiments
Introduction
Theory
Q
Lvap
m
Heat Capacity
The heat capacity, CT, is defined as the ratio of the energy input
by heating and the resulting temperature change (where Q is the
energy input and T is the temperature change):
Q
CT
T
dQ
CT
dT
CT
C
n
where n is the number of moles of material. The units used for are
JK−1mol−1.
Debye Theory
The classical theory explaining the specific heat capacity of a
solid assumes that in a crystal lattice there can only be vibrational
energy - there is no rotation or translational energy. Each atom on
a lattice site acts as an oscillator in three independent directions, so
there are 3N oscillators (where N is the number of atoms in the lat-
tice) each with two degrees of freedom (potential energy and kinetic
energy). Classically there is kT/2 of energy per degree of freedom
(where k is the Boltzman constant and T is the temperature) so:
1
3N 2 2 kT 3NkT
256 • Physics Lab Experiments
dQ
CV 3R
dT
hf
qE
k
and
2
q e qE T
CV 3Nk E
T
eq
T 2
E
–1
Common Experiments • 257
12pV f 3
D D ebye
c s3
where V is volume and cs is the speed of sound. Since the most vibra-
tions which can occur is 3N, the maximum number of states is also
3N:
f max 12pV f 2
0 c s3
df 3N
so:
3 3N c s3
f max
4pV
k qD hf max
13
3 hc s
qD
4p ka
1
n f hf kT
e –1
f max
E 0
hf D D ebye n f df
f max 12phV f3
E 0 c s3 e hf kT – 1
df
E
3
12p 4 T
CV Nk
T V 5 qD
W QV
Q It
So:
W It (V2 V1)
W VIt
PV nRT
a b c
then
and if: v
x
f xy or f
y
then:
2 2 2
sf sx s y
f x y
Experimental Methods
Figure 6.24: The Latent Heat Calorimeter. A diagram of the apparatus used in the Latent
Heat of Vaporization Liquid Nitrogen section.
the calorimeter were topped up with liquid nitrogen as the boil off
was considerable. When the boil off had reduced to a steady and low
rate, the experiment was started.
Firstly a background rate was taken, so as this contribution could
be subtracted from the results so that the contribution from the
heating coil alone could be calculated.
The heating coil was then turned on, and the current and (elec-
tric) potential difference were recorded using a multimeter. As the
liquid nitrogen in the calorimeter boiled off the volume was recorded
from the gas volume meter.
Using the current and (electric) potential difference supplied to
the heater and the time period over which the measurements were
taken, the heat dissipated by the heater was calculated. The volume
of N2 evolved as a result of this amount of heat was read from the gas
volume meter. This was converted into the number of moles as by
dividing by 24 (as outlined in the Number of Moles of a Gas section).
The mass of N2 produced was calculated using:
m
n
Mr
the next part could be performed more easily. The tube was then
removed from the calorimeter and a little liquid nitrogen was
added to it. The end of the rubber hose was corked so as to seal
the tube. As the liquid nitrogen boiled off, the bladder enlarged.
While the tube and hose were held upside down with the cork
removed the bladder was squeezed and the nitrogen gas expelled.
This was repeated until all the nitrogen had boiled off (i.e., the
bladder stopped enlarging when the hose was re-corked). Since
oxygen gas (1.43 kgm−3) is denser than nitrogen gas (1.25 kgm−3),
is sinks out of the tube, encouraged by the nitrogen, thus leav-
ing an oxygen-free atmosphere inside the stainless steel tube. The
hose was then immediately re-corked. This method proved very
successful at preventing the formation of liquid oxygen in the
stainless steel tube.
The stainless steel tube was then reinserted into the calorimeter
and any significant boil off from the liquid nitrogen in the calorim-
eter or the large dewar was topped up. When the boil off had settled
to a low steady value a background reading was started.
After the background reading was taken the cork was briefly
removed and the sample was added to the hose while squeezing the
bladder so as no air could enter the tube or hose. The hose was then
immediately re-corked. The sample dropped into the bottom of the
stainless steel tube and readings were taken from the gas volume
meter every 30 seconds until the rate of boil off had dropped below
that of the background rate.
The sample was not added to the hose earlier and left there
while the background reading was taken (with the hose bent over
to prevent the sample entering the calorimeter) since the ambient
air temperature inside the hose was not the same as the ambient air
temperature in the room. This would have meant the temperature
range over which the measurements were taken would not have been
known.
Each time a reading was taken a new background reading was
taken since it was more accurate than trying to keep the level of the
liquid nitrogen in the calorimeter and the dewar constant between
readings.
This part of the experiment was performed with samples of
copper, aluminium, graphite, zinc, lead, and rock salt.
264 • Physics Lab Experiments
does not build up. The possible formation of liquid oxygen (highly
flammable) was avoided by using the method outlined in the Specific
Heat Capacity Over Temperature Range 293 K → 77 K section.
Calculation of Results Over Temperature Range 293 K
77 K
The volume of N2 evolved over the previous 30 seconds was
calculated from the results. This was then plotted against time to
produce a graph showing the rate of N2 boil-off (see Table (6.1)
for an example of this graph). The background reading was then
added to the same graph. At the point where the two lines meet, it
is assumed that the sample is at 77 K (the same temperature as the
liquid nitrogen). A line of best fit was added to this graph using a
spreadsheet program to enable the intercept (the time when the
temperature of the sample reached 77 K) to be calculated exactly.
The original results from the gas volume meter were then plot-
ted against time, and a line of best fit added using the spreadsheet
program (see Table (6.2) for an example of this graph). When the
intercept from the first graph was substituted into the equation of
the second line of best fit, the reading on the gas volume meter at
the time when the temperature of the sample was 77 K was found.
The reading from the gas volume meter at t = 0 was then subtracted
from this to find the actual volume of N2 evolved.
This volume was then converted into the number of moles of N2
using the theory laid out in the Number of Moles of a Gas section
and into a number of grams.
Using the average value of latent heat of vaporization of N2 cal-
culated earlier this was converted into a number of Joules. Dividing
by the change in temperature (i.e. 293 K – 77 K = 216 K or 87 K
– 77 K = 10 K) gave a figure in Joules per Kelvin. Using the mass of
the sample and the Mr, the number of moles of the sample was cal-
culated and then the specific heat capacity in Joules per Kelvin per
mole (JK−1mol−1) was calculated.
Calculation of Results Over Temperature Range 86 K 77 K
The background reading observed during the first 600 seconds was
extrapolated to cover the time period of the rest of the experiment.
The difference between the actual volume and the extrapolated vol-
ume was then used in the same way as the change in volume was used
in the Calculation of Results Over Temperature Range 293 K → 77 K.
266 • Physics Lab Experiments
The average Lvap was found to be 249 kJkg−1 with a propagated error
calculated to be ± 4 kJkg−1.
Mass, Mr, Volume, and Density of Samples
Figure 6.25: Example of a graph used to calculate the time when the background level was
reached. This graph refers to the first result for copper over the range 293 K to 77 K. The
horizontal line shows the background rate and the curve is a line of best fit added to the
points showing the rate at which the nitrogen was evolved during the experiment.
The equation is the equation of the line of best fit.
268 • Physics Lab Experiments
Figure 6.26: Example of a graph used to calculate the volume of nitrogen evolved in the time
taken for the rate of production of nitrogen to reach background levels. This graph again
refers to the first result for copper over the temperature range 293 K to 77 K. The curve is the
line of best fit added to the points showing the reading on the gas volume meter. The equa-
tion is the equation of the line of best fit.
Accepted Results
Density
Discussion
Conclusion
Bibliography
REFERENCES
[1] http://www.cleapss.org.uk
[2] Mythbusters TV Program, Discovery Channel, Epi-
sode 63 (2006).
[3] British Compressed Gases Association Code of Prac-
tice CP30 “The Safe Use of Liquid Nitrogen Dewars
up to 50 liters’.
[4] http://www.esrl.noaa.gov/gmd/ccgg/trends/global.html.
[5] Bureau International des Poids et Mesures: www.
bipm.org/utils/common/pdf/its-90/SUPChapter4.pdf.
[6] M. M. J. French and Michael Hibbert. Making liquid
oxygen. Phys. Ed., 45:221, 2010.
[7] M. M. J. French. Mobile phone faraday cage. Phys.
Ed., 46:290, 2011.
[8] J. G. Ziegler and N. B. Nichols. Optimum settings for
automatic controllers. Transactions of the ASME.,
64:759, 1942.
[9] A. S. McCormack and K. R. Godfrey. Rule-based
autotuning based on frequency domain identification.
IEEE Transactions on Control Systems Technology,
6:1, 1998.
[10] http://www.amsterchem.com/onlinehelp.
php?page=welcome.htm.
[11] M. Vannoni and S. Straulino. Low-cost accelerom-
eters for physics experiments. Eur. J. Phys., 28:781,
2007.
[12] F. Rooney and W. Somers. Using the wiimote in intro-
ductory physics experiments. TCNJ J. Stud. Schol-
arsh., 12:1, 2010.
[13] M. D. Wheeler. Physics experiments with nintendo
wii controllers. Phys. Ed., 46:57, 2011.
[14] http://wiimotephysics.codeplex.com/.
[15] http://pdg.lbl.gov/2012/tables/rpp2012-sum-quarks.
pdf.
[16] J. A. Nelder and R. Mead. A simplex method for func-
tion minimization. Computer Journal, 7:308, 1965.
[17] http://www.openoffice.org.
274 • Physics Lab Experiments
[18] http://www.miktex.org.
[19] http://www.winedt.com.
[20] http://www.irfanview.com/.
[21] M. M. J. French et al. Changing fonts in education:
How the benefits vary with ability and dyslexia. The
Journal of Educational Research, 2013.
[22] http://www.prezi.com.
INDEX
W Y
“Wiimote” software, 108–109 Y error bars tab, 140
“position sensor” mode, 109
Window functions, 171–172 Z
Wii controllers, 108–109 Zero error, 117–118
Ziegler–Nichols method, 86
X Znamia Project, presentation, 198–199
X axis error bars, 143
X error bars tab, 140