Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Magnetic Force & the Determination of μ₀

Background

● Background Overview

By measuring the magnetic force between two parallel current-carrying conductors, the permeability of free space, μ0\mu_0, will be experimentally determined. μ0\mu_0 is a fundamental constant of nature, and, with its electric equivalent ε0\varepsilon_0 that we determined in an earlier experiment, we will now be able to determine the magnitude of the velocity of propagation of an electromagnetic wave --- the speed of light. These constants are related to the speed of light, c, by the following relation, derived from Maxwell’s equations which describe the behavior of electric and magnetic fields:

c=1ε0μ0c=\sqrt{\frac{1}{\varepsilon_0\mu_0}}

You will measure μ0\mu_0 here, and, with your previous measurement of ε0\varepsilon_0, determine the speed of light cc.

The 2019 redefinition of SI units defined exact values of fundamental constants including the electron charge ee, the speed of light cc, and Planck’s constant hh. The second is defined in terms of the frequency of a Cesium atomic clock. As a result, the value of μ0\mu_0, the magnetic constant, must now be experimentally determined.

Note: Until 2019, μ0\mu_0 was defined to be exactly 4π×1074\pi \times 10^{-7} T⋅m/A. This 2019 value is so close to the current best measurement, that you can use μ0=4π×107\mu_0=4\pi \times 10^{-7} T⋅m/A as the accepted value.

The magnetic field strength BB a distance rr from the center of a very long, straight conductor carrying a current II is given by:

B=μ0I2πrB=\frac{\mu_0 I}{2\pi r}

A conductor of length LL [m] carrying a current II [A] in a magnetic field of strength BB [T] experiences a force FBF_B [N], given by:

FB=IL×B=ILBsinθ\vec{F_B}=I \vec{L} \times \vec{B} =I L B \sin{\theta}

where θ\theta is the angle between the vectors B\vec{B} and L\vec{L}. If the magnetic field is produced by the current in a second conductor, the currents in the two conductors exert equal magnitude and opposite forces on each other.

For the case of two long, parallel conductors, each carrying the same current II and separated, center-to-center, by the distance dcenter-to-centerd_\text{center-to-center}, the force between the two conductors is the force on one in the magnetic field of the second. Thus if (2) is used in (3), we obtain for the force:

FB=μ02πLI2dcenter-to-center\boxed{ F_B=\frac{\mu_0}{2\pi}\frac{L I^{2}}{d_\text{center-to-center}} }

By measuring the force between two such conductors, the value of μ0\mu_0 can be determined. By considering the parameters of the apparatus to be used, it can be seen that the magnitude of the force between the two conductors is quite small and will require some care to accurately measure. If LL is about 30 cm, II about 10 A, and dcenter-to-centerd_\text{center-to-center} about 2 mm, the force would be in the range of 10-4 N, or the weight of a few milligrams of mass. Since this force is well below the weight of a reasonable conductor to be used for the experiment, a counter weighted balance system will be used to provide the mechanical support for the movable conductor. This is essentially the same apparatus as Electric Force & ε₀ Lab. The movable conductor can then be loaded by a known mass whose weight can then be matched by the small magnetic force. The apparatus is schematically illustrated from a top-down perspective in Figure 1.

Apparatus dimensions from a top-down perspective, similar to the apparatus setup from Electric Force & the Determination of ε₀.

Figure 1:Apparatus dimensions from a top-down perspective, similar to the apparatus setup from Electric Force & the Determination of ε₀.

As before, similar to the setup and measurement methods from the Electric Force & ε₀ Lab, an optical lever, telescope and scale are used to observe the deflection of the balance, which carries the movable, upper conductor. The magnitude of the magnetic force can be measured when the magnetic force has restored the loaded system to its original position. In this case we make the same magnitude current flow in opposite directions in the two conductors so the force between the two conductors is repulsive. Thus on the upper conductor, there is a gravitational force down (the weight of the applied mass) and a magnetic force upward which we will adjust by controlling the current equal to the gravitational force. More conveniently than the electric force version of the balance beam, this equilibrium is stable!

With the required balancing current determined along with the other dimensions of the apparatus, we can now determine the value of μ0\mu_0.

● Equipment

Table 1:Equipment

CategoryItems
Optical Measurement• Telescope with crosshair
• Centimeter-scale ruler mounted on vertical pole
Power Supply• AC power from wall outlet controlled by large cylindrical transformer/potentiometer (0–100% of wall power, ~0–20 A) Dial markings are not accurate; ignore scale on knob
• Voltage transformer (small metal cube) reducing AC voltage to ~6 V
Current Measurement• Fluke multimeter used as AC ammeter (A~\tilde{\text{A}}, alt text), connected in series with circuit
• One lead in 10 A port alt text, other lead to COM port alt text
Masses & Tools• Small masses (5, 20, 50, 100, 200 mg)
• Plastic tweezers for handling masses
Parallel-Conductor Apparatus• Bottom conductor fixed in place
• Top conductor free to swing vertically when current or mass changes
• Adjust so conductors are parallel when 0 mg is applied
• Mirror used with telescope to read ruler and determine balance beam angle
• Beam lift knobs to reset top conductor position
• Leveling screws to level apparatus
Electrical Connections• 4x banana plug wires (12 AWG)
 • 2x 1 ft
 • 2x 3 ft
• Connect transformer \rightarrow ammeter \rightarrow parallel-conductor apparatus
Additional Equipment• Protective box with mirror-to-scale distance written on it

● Adjustment of apparatus

The apparatus has been carefully adjusted before your lab and should not require further significant adjustments. This section describes how the apparatus was prepared. If something seems to need adjusting, see the lab instructor. The apparatus is very similar to the apparatus setup from Electric Force & ε₀ Lab.

  1. The beam lift provides a definite location for the beam and thus guarantees continued alignment of the parallel conductors. Use the beam lift each time you change weights or relocate the counterweight.

  2. The fixed conductor can be adjusted vertically, the movable conductor horizontally, so that the conductors are parallel. If the conductors are NOT parallel with no mass in place, seek the instructor’s help.

  3. The counterweight behind the mirror can be used to change the equilibrium separation of the conductors.

  4. Level the apparatus with the adjusting screws so that it sits securely on the table.

  5. Mirror-to-scale distance was measured from rear of mirror (reflective surface) to roughly the middle of S1S_1 and S0S_0 on the scale (ruler).

Demo Video: Setup & Procedure

Some clarifications, additions, or corrections since this video is slightly outdated:

Demonstration video. Reviews overall setup and procedure. *Yes audio*

If embedding is broken, follow: https://www.youtube.com/watch?v=EpYTWOcaFSU

Experimental Procedure

● PRECAUTIONS

● Preview

● Preliminary Setup

If bb is the mirror-to-scale distance and aa is the length of the frame (see Figure 1), then a straightforward geometrical analysis (like in Electric Force & ε₀ Lab) will show that the vertical displacement of the conductor, yy, is given by:

y=Da2by=\frac{D a}{2 b}

where DD is the change in scale reading from the equilibrium value and the value when the two conductors are in contact (the contact reading is determined by adding a mass to the pan which depresses the beam until contact occurs). The factor of 2 results from the fact that the optical path reflected off of the mirror to the scale rotates through an angle twice that of the beam holding the movable conductor.

However, the change in scale reading doesn’t account for the conductors’ thickness. Since the B-field generated by current flowing in the cylindrical conductors is organized about their center, we must add in two radii to get the center-center separation between the two conductors. This gives us:

dcenter-to-center=Da2b+2rconductor\boxed{ d_\text{center-to-center}=\frac{D a}{2 b}+2r_\text{conductor} }

where rconductorr_\text{conductor} is the radius of the conductor (multiplied by two to deal with the fact that we’re dealing with two conductors).

  1. Investigate the use of the telescope and scale so that the rotation of the frame can be measured in terms of scale divisions (see Figure 2).

Schematic of the measuring apparatus, similar to the setup and measurement methods from Electric Force & ε₀ Lab.

Figure 2:Schematic of the measuring apparatus, similar to the setup and measurement methods from Electric Force & ε₀ Lab.

  1. Do not turn on the power supply until requested in the procedure later in ● Experimental Data Collection.

  2. Create a common data table with:

    • aa: length of frame, 0.215m0.215\,\text{m}

    • bb: mirror-to-scale distance (on top of protective box)

    • LL: wire/bar length 0.265m0.265\,\text{m}

    • rconductorr_\text{conductor}: wire radius 1.6mm1.6\,\text{mm}

    • S0S_0: equilibrium crosshair position

    • S1S_1: bars touching crosshair position

    • DD: ruler (a.k.a. scale) distance between S0S_0 and S1S_1

    • dcenter-to-centerd_\text{center-to-center}: length of frame

    • gg: acceleration due to gravity for Fairfield, CT, 9.803m/s29.803\,\text{m/s}^2

    • μ0,accepted\mu_{0\text{,accepted}}: accepted value of the permeability of free space, 4π×107T m/A or N/A24\pi \times 10^{-7}\,\text{T m/A or N/}\,\text{A}^2

    • cacceptedc_{\text{accepted}}: accepted value of the speed of light, 2.9979×108m/s2.9979\times10^8\,\text{m/s}

    • Also add in the plotting-derived (i.e. from LINEST) ε0\varepsilon_0 from the first lab for each member of your lab group

      • Of your group members’ values, take the average (ε0\overline \varepsilon_0), and use the standard deviation (σε0\sigma \varepsilon_0) of your groups’ values as the uncertainty for later speed of light use.

      • If you are completing the lab individually, use your single plotting-derived ε0\varepsilon_0 and its associated uncertainty δε0\delta \varepsilon_0

  1. In order to measure the equilibrium separation distance, dcenter-to-centerd_\text{center-to-center}, of the conductors, two readings are made.

    • The first is the scale reading, S0S_0, at the equilibrium position, i.e. when no mass has been applied to the pan on the movable conductor.

    • The second scale reading, S1S_1, is made when the top conductor contacts the lower conductor (Figure 2). Place a sufficient mass in the pan on the top conductor to make it contact the lower, stationary conductor. Record S1S_1.

      • ⚠️ POWER MUST BE OFF FOR THIS STEP ⚠️

      • You will periodically verify that S0S_0 and S1S_1 are not changed throughout the experiment.

    • Determine D=S1S0D = | S_1 - S_0 | from two scale readings (reminders: absolute value there is to represent the total distance between S0S_0 and S1S_1. If your crosshair crosses 0, make sure to consider the negative).

    • Calculate and record the conductor separation distance, dcenter-to-centerd_\text{center-to-center}. Refer to Figure 1, Figure 2, (6).

  2. Prepare a data table with columns for the variables below. You’ll have rows for each trial as well as for averages. Reminder, you will conduct 27 trials total by conducting three rounds of sequentially changing the applied masses:

    • Trial number

    • Lab member’s initials (person looking through telescope)

    • mappliedm_{\text{applied}}: Applied mass

    • FBF_B: Applied magnetic force What does this equal?

    • IminI_\text{min}: Minimum current II required to return to the equilibrium position

    • ImaxI_\text{max}: Maximum current II required to return to the equilibrium position

    • II: Current II required to return to the equilibrium position

    • δI\delta I: Estimated current uncertainty (to be assumed as majority source of uncertainty for today)

    • I2I^2: Current squared

    • δI2\delta I^2: Current squared

    • μ0,experimental\mu_{0\text{,experimental}}: Experimental value of μ0\mu_0

    • μ0,experimental,max\mu_{0\text{,experimental,max}}: Maximized experimental value based on δI2\delta I^2

    • δμ0,experimental\delta \mu_{0\text{,experimental}}: Uncertainty in experimental value

    • Δμ0,experimental vs. accepted\Delta \mu_{0\text{,experimental vs. accepted}}: Magntitude of difference between experimental and accepted values of μ0\mu_0

● Experimental Data Collection

  1. With the power off 🟥, make sure the conductor beam is still at equilibrium, then place the current trial’s mass mappliedm_\text{applied} on the upper conductor (see full list in Table 2). The top bar will swing downwards. In following steps, you will apply a current to counteract the gravitational force. Record the trial’s mass value and calculate FBF_\text{B} (reminder, use SI units; these masses are in units of milligrams; hint: what does the magnetic force balance?).

    Table 2:Experimental Trials Breakdown

    Student on telescopeTrial Applied Masses (mg)POWER ON 🟢 / Off 🟥
    Student #1200, 175, 150, 125, 100, 75, 50, 25ON 🟢
    Student #1 (and #2)0 (ensure S0S_0 has not changed)Off 🟥
    Student #2200, 175, 150, 125, 100, 75, 50, 25ON 🟢
    Student #2 (and #3)0 (ensure S0S_0 has not changed)Off 🟥
    Student #3 (or #1
    again for a group of two)
    200, 175, 150, 125, 100, 75, 50, 25ON 🟢
    Student #3 (and #1)0 (ensure S0S_0 has not changed)Off 🟥
  2. Calculate FBF_B (hint: what does this equal?) and determine current II required to balance the weight of the applied mass and return to equilibrium by finding a current range. Make sure the transformer setting is zero before turning on the power. Turn on 🟢 the power supply and slowly increase the current until the bars return to parallel (i.e. crosshair in telescope back to S0S_0, as it was with no applied mass on the top). Determine this by having an observer watching the scale reading with the telescope during this process. The telescope observer should be calling out instructions to the power supply operator to slowly approach the S0S_0 value.

    • When the S0S_0 value is approximately reached, call out to the power supply operator.

    • Power supply operator shall decrease the current until the telescope observer is no longer confident the crosshair is at the S0S_0 position, at which point the telescope observer shall call for a minimum current reading from the multimeter and record the lower end of the range of current IminI_\text{min}.

    • Power supply operator shall then increase the current until the telescope observer is no longer confident the crosshair is at the S0S_0 position, at which point the telescope observer shall call for a maximum current reading from the multimeter and record the upper end of the range of current ImaxI_\text{max}.

    • As soon as you make these readings, reduce the transformer setting to zero and turn off 🟥 the power switch. Under no circumstances should the current exceed 20 A.

    • Calculate applied current II as the average of the min and max currents, I=(Imax+Imin)/2I\,=\,(\,I_\text{max}\, +\,I_\text{min}\,)\,/\,2

    • Calculate applied current uncertainty δI\delta I as half the difference between max and min current, δI=(ImaxImin)/2\delta I\,=\,(\,I_\text{max}\,-\,I_\text{min}\,)\,/\,2

    • Calculate I2I^2 as well as δI2\delta I^2 by maximizing and taking the difference (i.e. δI2=(I+δI)2I2\delta I^2 = (I + \delta I)^2 - I^2)

  3. For the current trial, calculate μ0,experimental\mu_{0\text{,experimental}} from (4) and determine its uncertainty by maximizing by your measurement uncertainties (primarily due to current today) and taking the difference (i.e. δμ0,experimental=μ0,experimental,maxμ0,experimental\delta \mu_{0\text{,experimental}} = \mu_{0\text{,experimental,max}} - \mu_{0\text{,experimental}}, hint: should be a positive value).

  4. Also calculate the difference Δμ0,experimental vs. accepted\Delta \mu_{0\text{,experimental vs. accepted}} between your experimental and accepted permeability of free space values for quick comparisons.

  5. With power off, replace the applied mass for the next trial, and repeat the previous steps in ● Experimental Data Collection to find the current required to balance the system again, and subsequently calculate μ0±δμ0\mu_0 \pm \delta \mu_0 again.

● Data Analysis

○ Averages

  1. Calculate your average values from everyone’s trials (ignoring those that give a divide-by-zero error):

    • Average experimental value μ0,experimental,avg\mu_{0\text{,experimental,avg}}

    • Average experimental uncertainty δμ0,experimental,avg\delta \mu_{0\text{,experimental,avg}}

    • Δμ0,experimental vs. accepted\overline{\Delta \mu_{0\text{,experimental vs. accepted}}}, the average difference between experimental and accepted μ0\mu_0

  2. If, after performing the graphical data analysis below, you find some or all of the data unacceptable, repeat from Step 4, checking that the telescope and scale were not disturbed during your measurements.

○ Graphical

  1. The graphical display of data permits the comparison of all the values and associated errors at once. Points that depart markedly from the general trend of the data are quickly detected. We expect from theory (4) that magnetic force and current are related like:

    FB=kI2F_B=k I^{2}

    Using all of your data points from all rounds and trials, including those at equilibrium (0mg0\,\text{mg}) that you should have ignored earlier when you calculated μ0\mu_0 for each individual trial and the overall average:

    A) SCATTER PLOT 1 --- FBF_B vs. I2I^{2} (i.e. FBF_B as ordinate (yy) and I2I^{2} as abscissa (xx)). Fit a trend line through the all these points, display the trendline equation on the chart, and confirm the slope of the line kk matches what you found with the LINEST() function (for review, see Plotting in Excel). From (4) and (7), the slope value of kk is given by (8).

    k=μ0,slope-derived2πLdcenter-to-centerk=\frac{\mu_{0\text{,slope-derived}}}{2\pi}\frac{L}{d_\text{center-to-center}}
    • Calculate the plotting or slope-derived μ0,slope-derived\mu_{0\text{,slope-derived}} by rearranging (8) and using this experimentally determined value kk.

    • Similar to earlier, calculate δμ0,slope-derived\delta \mu_{0\text{,slope-derived}} by maximizing by the slope uncertainty (output from the LINEST() function, Plotting in Excel), and subsequently taking the difference: δμ0,slope-derived=μ0,slope-derived,maxμ0,slope-derived\delta \mu_{0\text{,slope-derived}} = \mu_{0\text{,slope-derived,max}} - \mu_{0\text{,slope-derived}}

    • Also calculate the difference between the slope-derived and accepted permeability of free space values to quickly compare whether your results agree with the accepted value.

    B) SCATTER PLOT 2 --- FBF_B vs. II. Fit a quadratic (a.k.a. polynomial order of 2) trend line, and display the equation.

○ Speed of Light

  1. Using the average of your group’s experimentally-determined plotting-derived values of the electric constant (ε0\overline \varepsilon_0) from the Electric Force & ε₀ Lab, and the plotting-derived value of μ0,slope-derived\mu_{0\text{,slope-derived}} determined above, calculate your estimate for the speed of light cc with (1).

    • As should be in your common data table, and if not yet completed, average the group members’ plotting-derived values of the electric constant (ε0\overline \varepsilon_0), and use the standard deviation σε0\sigma \varepsilon_0 of your groups’ values as its uncertainty.

      • If you are completing the lab individually, use your single plotting-derived ε0\varepsilon_0 and its associated uncertainty δε0\delta \varepsilon_0

  2. Estimate your uncertainty in cc by MINIMIZING and taking the difference. Reminder, for δμ0\delta \mu_0 from today’s error propagation, and δε0=σε0\delta \varepsilon_0 = \sigma \varepsilon_0. (Note: we’re minimizing instead of maximizing since there would be a chance of taking the square root of a negative number):

    δc=c1(ε0+δε0)×(μ0,slope-derived+δμ0,slope-derived)\delta c = c - \sqrt{\frac{1}{(\overline{\varepsilon_0}+\delta \varepsilon_0) \times (\mu_{0\text{,slope-derived}} + \delta \mu_{0\text{,slope-derived}})}}

● Summary and Cleanup

  1. Create a summary table of your data (e.g. average μ0\mu_0, slope-derived μ0\mu_0, and speed of light values with their uncertainties, as well as the differences between experimental and accepted values).

  2. When you are finished, reset your experimental setup before leaving.

Post-Lab Submission --- Interpretation of Results

● Finalized Spreadsheets

● Post-lab Writeup

The Whiteboard

2026 updates/notes on whiteboard summaries:

Overview. LINEST() function.

Figure 3:Overview. LINEST() function.

Separation distance equation; multimeter settings.

Figure 4:Separation distance equation; multimeter settings.

Plot notes; Speed of Light calculations.

Figure 5:Plot notes; Speed of Light calculations.