\(\boldsymbol{\vec{K}} = yz \,\boldsymbol{\hat{x}} + xz \,\boldsymbol{\hat{y}}\)
Main ideas
Finding potential functions.
Students love this activity. Some groups will finish in 10 minutes or less;
few will require as much as 30 minutes.
*
Prerequisites
Fundamental Theorem for line integrals
The Murder Mystery Method
Warmup
none
Props
whiteboards and pens
Wrapup
Revisit integrating conservative vector fields along various paths, including
reversing the orientation and integrating around closed paths.
Details
In the Classroom
We recommend having the students work in groups of 2 on this activity, and not
having them turn anything in.
Most students will treat the last example as 2-dimensional, giving the answer
\(xyz\). Ask these students to check their work by taking the gradient; most
will include a \(\boldsymbol{\hat{z}}\) term. Let them think this through. The correct answer
of course depends on whether one assumes that \(z\) is constant; we have
deliberately left this ambiguous.
It is good and proper that students want to add together multivariable terms. Keep returning to the gradient, something they know well. It is better to discover the guidelines themselves.
Subsidiary ideas
3-d vector fields do not necessarily have a \(\boldsymbol{\hat{z}}\)-component!
Homework
A challenging question to ponder is why a surface fails to exist for nonconservative fields. Using an example such as \(y\,\boldsymbol{\hat{x}}+\boldsymbol{\hat{y}}\), prompt students to plot the field and examine its magnitude at various locations. Suggest piecing together level sets. There is serious geometry lurking that entails smoothness. Wrestling with this is healthy.
Essay questions
Write 3-5 sentences describing the connection between derivatives and integrals in the single-variable case. In other words, what is the one-dimensional version of MMM? Emphasize that much of vector calculus is generalizing concepts from single-variable theory.
Enrichment
The derivative check for conservative vector fields can be described using the
same type of diagrams as used in the Murder Mystery Method; this is just
moving down the diagram (via differentiation) from the row containing the
components of the vector field, rather than moving up (via integration). We
believe this should not be mentioned until after this lab.
When done in 3-d, this makes a nice introduction to curl --- which
however is not needed until one is ready to do Stokes' Theorem. We would
therefore recommend delaying this entire discussion, including the 2-d case,
until then.
Work out the Murder Mystery Method using polar basis vectors, by reversing the
process of taking the gradient in this basis.
Revisit the example in the Ampère's Law lab, using the Fundamental Theorem
to explain the results. This can be done without reference to a basis, but
it is worth computing \(\boldsymbol{\vec\nabla}\phi\) in a polar basis.
to perform a magnetic vector potential calculation using the superposition principle;
to decide which form of the superposition principle to use, depending on the dimensions of the current density;
how to find current from total charge \(Q\), period \(T\), and the geometry of the problem, radius \(R\);
to write the distance formula \(\vec{r}-\vec{r'}\) in both the numerator and denominator of the superposition principle in an appropriate mix of cylindrical coordinates and rectangular basis vectors;
Student discuss how many paths can be found on a map of the vector fields \(\vec{F}\) for which the integral
\(\int \vec{F}\cdot d\vec{r}\) is positive, negative, or zero. \(\vec{F}\) is conservative.
They do a similar activity for the vector field \(\vec{G}\) which is not conservative.
Students examine a plastic “surface” graph of the gravitational potential energy of an Earth-satellite system to explore the properties of gravitational potential energy for a spherically symmetric system.
Students hold rulers and meter sticks to represent a vector field. The instructor holds a hula hoop to represent a small area element. Students are asked to describe the flux of the vector field through the area element.
Students work in small groups to use the Biot-Savart law
\[\vec{B}(\vec{r}) =\frac{\mu_0}{4\pi}\int\frac{\vec{J}(\vec{r}^{\,\prime})\times \left(\vec{r}-\vec{r}^{\,\prime}\right)}{\vert \vec{r}-\vec{r}^{\,\prime}\vert^3} \, d\tau^{\prime}\]
to find an integral expression for the magnetic field, \(\vec{B}(\vec{r})\), due to a spinning ring of charge.
In an optional extension, students find a series expansion for \(\vec{B}(\vec{r})\) either on the axis or in the plane of the ring, for either small or large values of the relevant geometric variable.
Add an extra half hour or more to the time estimate for the optional extension.
Students are shown a topographic map of an oval hill and imagine that the classroom is on the hill. They are asked to point in the direction of the gradient vector appropriate to the point on the hill where they are "standing".
In this introductory lecture/SWBQ, students are given a picture as a guide. They then write down an algebraic expression for the vector differential in rectangular coordinates for coordinate equals constant paths.
This activity can be done as a mini-lecture/SWBQ as an introduction to Vector Differential--Curvilinear where students find the vector differential in cylindrical and spherical coordinates..
vector differentialrectangular coordinatesmath Found in: Static Fields, AIMS Maxwell, Vector Calculus II, Surfaces/Bridge Workshop, Problem-Solving course(s)Found in: Integration Sequence sequence(s)
Find the general form for \(d\vec{r}\) in cylindrical coordinates by determining \(d\vec{r}\) along the specific paths below.
Path 1 from \((s,\phi,z)\) to \((s+ds,\phi,z)\): \[d\vec{r}=\hspace{35em}\]
Path 2 from \((s,\phi,z)\) to \((s,\phi,z+dz)\): \[d\vec{r}=\hspace{35em}\]
Path 3 from \((s,\phi,z)\) to \((s,\phi+d\phi,z)\): \[d\vec{r}=\hspace{35em}\]
If all three coordinates are allowed to change simultaneously, by an
infinitesimal amount, we could write this \(d\vec{r}\) for any path as:
\[d\vec{r}=\hspace{35em}\]
This is the general line element in cylindrical coordinates.
Figure 1: \(d\vec{r}\) in cylindrical coordinates
Spherical Coordinates:
Find the general form for \(d\vec{r}\) in spherical coordinates by determining \(d\vec{r}\) along the specific paths below.
Path 1 from \((r,\theta,\phi)\) to \((r+dr,\theta,\phi)\):
\[d\vec{r}=\hspace{35em}\]
Path 2 from \((r,\theta,\phi)\) to \((r,\theta+d\theta,\phi)\):
\[d\vec{r}=\hspace{35em}\]
Path 3 from \((r,\theta,\phi)\) to \((r,\theta,\phi+d\phi)\): (Be careful, this is a tricky one!) \[d\vec{r}=\hspace{35em}\]
If all three coordinates are allowed to change simultaneously, by an
infinitesimal amount, we could write this \(d\vec{r}\) for any path as:
\[d\vec{r}=\hspace{35em}\]
This is the general line element in spherical coordinates.
Figure 2: \(d\vec{r}\) in spherical coordinates
Instructor's Guide
Main Ideas
This activity allows students to derive formulas for \(d\vec{r}\) in cylindrical, and spherical coordinates, using purely geometric reasoning. These formulas form the basis of our unified view of all of vector calculus, so this activity is essential. For more information on this unified view, see our publications, especially:
Using differentials to bridge the vector calculus gap
Students' Task
Using a picture as a guide, students write down an algebraic expression for the vector differential in different coordinate systems (cylindrical, spherical).
Introduction
Begin by drawing a curve (like a particle trajectory, but avoid "time" in the language) and an origin on the board. Show the position vector \(\vec{r}\) that points from the origin to a point on the curve and the position vector \(\vec{r}+d\vec{r}\) to a nearby point. Show the vector \(d\vec{r}\) and explain that it is tangent to the curve.
For the case of cylindrical coordinates, students who are pattern-matching will write
\(d\vec{r} = dr\, \hat{r} + d\phi\, \hat{\phi} + dz\, \hat{z}\). Point out that \(\phi\) is dimensionless and that path two is an arc with arclength \(r\, d\phi\).
Some students will remember the formula for arclength, but many will not. The following sequence of prompts can be helpful.
What is the circumference of a circle?
What is the arclength for a half circle?
What is the arclength for the angle \(\pi\over 2\)?
What is the arclength for the angle \(\phi\)?
What is the arclength for the angle \(d\phi\)?
For the spherical case, students who are pattern matching will now write
\(d\vec{r} = dr\, \hat{r} + d\phi\, \hat{\phi} + d\theta\, \hat{\theta}\). It helps to draw a picture in cross-section so that they can see that the circle whose arclength gives the coefficient of \(\hat{\theta}\) has radius \(r\sin\theta\). It can also help to carry around a basketball to write on to talk about the three dimensional geometry of this problem.
Wrap-up
The only wrap-up needed is to make sure that all students have (and understand the geometry of!) the correct formulas for \(d\vec{r}\).
Students write python programs to compute the potential due to a square of surface charge, and then to visualize the result. This activity can be used to introduce students to the process of integrating numerically.
Students work in small groups to use the superposition principle
\[V(\vec{r}) =\frac{1}{4\pi\epsilon_0}\int\frac{\rho(\vec{r}^{\,\prime})}{\vert \vec{r}-\vec{r}^{\,\prime}\vert} \, d\tau^{\prime}\]
to find an integral expression for the electrostatic potential, \(V(\vec{r})\), everywhere in space, due to a ring of charge.
In an optional extension, students find a series expansion for \(V(\vec{r})\) either on the axis or in the plane of the ring, for either small or large values of the relevant geometric variable. Add an extra half hour or more to the time estimate for the optional extension.
Students solve numerically for the potential due to a spherical shell of charge. Although this potential is straightforward to compute using Gauss's Law, it serves as a nice example for numerically integrating in spherical coordinates because the correct answer is easy to recognize.
Consider a column of atoms each of mass \(M\) at temperature \(T\) in
a uniform gravitational field \(g\). Find the thermal average
potential energy per atom. The thermal average kinetic energy is
independent of height. Find the total heat capacity per atom. The
total heat capacity is the sum of contributions from the kinetic
energy and from the potential energy. Take the zero of the
gravitational energy at the bottom \(h=0\) of the column. Integrate
from \(h=0\) to \(h=\infty\). You may assume the gas is ideal.
Consider a system consisting of four point charges arranged on the corners of a square in 3D Cartesian space of coordinates \((x,y,z)\).
Write a python function that returns the potential at any point in space caused by four equal point charges forming a square. Make the sides of the square parallel to the \(x\) and \(y\) axes and on the \(z=0\) plane.
To do this you will need the expression for a the potential due to a single point charge \(V= \frac{k_Cq}{r}\) where \(r\) is the distance from the point charge. You will also need to use the fact that the total potential is the sum of the potentials due to each individual point charge.
It is important that we ask students first to create a function for the potential, and only then try to visualize the potential. This allows students to reason about the computation for a single point in space (defined in their choice of coordinate systems).
Once you have written the above function, use it to plot the electrostatic potential versus position along the three cartesian axes.
Since the students have already written a function for their potential, they can create a plot by creating an array for \(x\) (or \(y\), or \(z\)), and then passing that array to their function, along with scalars for the other two coordinates. Many students will discover this simply by modifying an example script they find on the web, replacing \(\sin(x)\) or similar with their function. It is well worth showing this easier approach to students who attempt who attempt to write a loop in order to compute the potential at each point in space.
We ask students to explicitly plot the potential along axes because students seldom spontaneously think to create a 1D plot such as this.
Label your axes.
Work out the first non-zero term in a power series approximation for the potential at large \(x\), small \(x\), etc. Plot these approximations along with your computed potential, and verify that they agree in the range that you expect. Useful 1\(^{st}\) order Taylor expansions are:
\begin{eqnarray}
\sqrt{1+\epsilon} &\sim& 1+\frac{\epsilon}{2} \\
\frac{1}{1+\epsilon} &=& 1-\epsilon
\end{eqnarray}
where \(\epsilon\) is a small quantity.
This may need to be omitted on the first Tuesday of class, since students probably will not yet have seen power series approximations. It may work in this case to at least talk about what is expected at large distance, since "it looks like a point charge" is reasoning students do make.
Students struggle with the \(x\) approximations (assuming the square is in the xy plane). Each pair will probably need to have a little lecture on grouping terms according to the power of \(x\), and keeping only those terms for which they have every instance.
Extra fun
Create one or more different visualizations of the electrostatic potential. For example a 2D representation in the \(z=0\) plane.
More extra fun
Create a plot of the potential along a straight line that is not one of the axes. Hint: start from a line on the \(z=0\) plane, then try a random straight line. You can use your browser for help.
Even more extra fun
Move the charges around (e.g., off the \(z=0\) plane) and see what happens to your graphs
Dipole fun
Repeat the above (especially the limiting cases!) for four point charges in which half are positive and half negative, with the positive charges neighbors.
Common visualizations for 2D slices of space include contour plots, color plots, and "3D plots". Another option (less easy) would be to visualize an equipotential surface in 3 dimensions. It is worth reminding students to consider other planes than those at \(x=0\), \(y=0\), and \(z=0\).
Quadrupole fun
Repeat the above (especially the limiting cases!) for four point charges in which half are positive and half negative, with the positive charges diagonal from one another. It will help in this case to place the charges on the axes (rotating the square by 45 degrees), since otherwise the potential on each axis will be zero.
electrostatic potentialpython Found in: Computational Physics Lab II course(s)Found in: Computational integrating charge distributions sequence(s)
Students are prompted to consider the scalar superposition of the electric potential due to multiple point charges. First a single point charge is discussed, then four positive charges, then an electric quadrupole. Students draw the equipotential curves in the plane of the charges, while also considering the 3D nature of equipotentials.
Students consider how changing the volume of a system changes the internal energy of the system. Students use plastic graph models to explore these functions.
Students learn about the geometric meaning of the amplitude and period parameters in the sine function. They also practice sketching the sum of two functions by hand.
Many students do not know what it means to add two functions graphically. Students are shown graphs of two simple functions and asked to sketch the sum.
Students draw the 3D graphs of equations using three variables. They make choices for drawing a stack of curves in parallel planes and a curve in a perpendicular plane (e.g. substituting in values for \(x\), \(y\), or \(z\). )
One way to write volume charge densities without using piecewise functions is to use step \((\Theta)\) or \(\delta\) functions.
Consider a spherical shell with charge density
\[\rho (\vec{r})=\alpha3e^{(k r)^3} \]
between the inner radius \(a\) and the outer radius \(b\). The charge density is zero everywhere else.
(2 pts) What are the dimensions of the constants \(\alpha\) and \(k\)?
(2 pts) By hand, sketch a graph the charge density as a function of \(r\) for \(\alpha > 0\) and \(k>0\) .
(2 pts) Use step functions to write this charge density as a single function valid everywhere in space.
Using either this
Geogebra applet
or this
Mathematica notebook, explore the wave functions on a ring. (Note: The Geogebra applet may be a little easier to use and understand and is accessible if you don't have access to Mathematica, but it is more limited in the wave functions that you can represent. Also, the animation is pretty jumpy in some browsers, especially Firefox. Imagine that the motion is smooth.)
Look at graphs of the following states
\begin{align}
\Phi_1(\phi)&=\frac{1}{\sqrt{2}}(\left|{2}\right\rangle +\left|{-2}\right\rangle )\\
\Phi_2(\phi)&=\frac{1}{\sqrt{2}}(\left|{2}\right\rangle -\left|{-2}\right\rangle )\\
\Phi_3(\phi)&=\frac{1}{\sqrt{2}}(\left|{2}\right\rangle +i\left|{-2}\right\rangle )
\end{align}
Write a short description of how these states differ from each other.
Find a state for which the probability density does not depend on time. Write the state in both ket and wave function notation. These are called stationary states. Generalize your result to give a characterization of the set of all possible states that are stationary states.
Find a state that is right-moving. Write the state in both ket and wave function notation. Generalize your result to give a characterization of the set of all possible states that are right-moving.
Find a state that is a standing wave. Write the state in both ket and wave function notation. Generalize your result to give a characterization of the set of all possible states that are standing waves.
Students use known algebraic expressions for vector line elements \(d\boldsymbol{\vec{r}}\) to
determine all simple vector area \(d\boldsymbol{\vec{A}}\) and volume elements \(d\tau\) in cylindrical and spherical coordinates.
This activity is identical to
Scalar Surface and Volume Elements except uses a vector approach to find directed surface and volume elements.
A short improvisational role-playing skit based on the Star Trek series in which students explore the definition and notation for position vectors, the importance of choosing an origin, and the geometric nature of the distance formula.
\[\vert\vec{r}-\vec{r}^\prime\vert=\sqrt{(x-x^\prime)^2+(y-y^\prime)^2-(z-z^\prime)^2}\]
A circular cylinder of radius \(R\)
rotates about the long axis with angular velocity \(\omega\). The
cylinder contains an ideal gas of atoms of mass \(M\) at temperature
\(T\). Find an expression for the dependence of the concentration
\(n(r)\) on the radial distance \(r\) from the axis, in terms of
\(n(0)\) on the axis. Take \(\mu\) as for an ideal gas.
(4pts)
Find the electric field around an infinite, uniformly charged,
straight wire, starting from the following expression for the electrostatic
potential:
\begin{equation*}
V(\vec r)=\frac{2\lambda}{4\pi\epsilon_0}\, \ln\left( \frac{ s_0}{s} \right)
\end{equation*}
The concentration of potassium
\(\text{K}^+\) ions in the internal sap of a plant cell (for example,
a fresh water alga) may exceed by a factor of \(10^4\) the
concentration of \(\text{K}^+\) ions in the pond water in which the
cell is growing. The chemical potential of the \(\text{K}^+\) ions is
higher in the sap because their concentration \(n\) is higher there.
Estimate the difference in chemical potential at \(300\text{K}\) and
show that it is equivalent to a voltage of \(0.24\text{V}\) across the
cell wall. Take \(\mu\) as for an ideal gas. Because the values of the
chemical potential are different, the ions in the cell and in the pond
are not in diffusive equilibrium. The plant cell membrane is highly
impermeable to the passive leakage of ions through it. Important
questions in cell physics include these: How is the high concentration
of ions built up within the cell? How is metabolic energy applied to
energize the active ion transport?
David adds
You might wonder why it is even remotely plausible to consider the
ions in solution as an ideal gas. The key idea here is that the ideal
gas entropy incorporates the entropy due to position dependence, and
thus due to concentration. Since concentration is what differs between
the cell and the pond, the ideal gas entropy describes this pretty
effectively. In contrast to the concentration dependence, the
temperature-dependence of the ideal gas chemical potential will not be
so great.
Write the equation for the electrostatic potential due to a point charge.
Instructor's Guide
Prerequisite Knowledge
Students will usually have seen the electrostatic potential due to a point charge in their introductory course, but may have trouble recalling it.
Whole-Class Conversations
As students try to remember the formula, many will conflate potential, potential energy, force, and electric field. Their answers may have some aspects of each of these. We use this question to get the iconic equation into the students' working memory in preparation for subsequent activities. This question also be used to help student disambiguate these different physical quantities.
Correct answers you're likely to see
\[V=\frac{kq}{r}\]
\[V=\frac{1}{4\pi\epsilon_0}\frac{q}{r}\]
You may want to discuss which constants to use in which contexts, e.g. \(k\) is short and easy to write, but may be conflated with other uses of \(k\) in a give problem whereas \(\frac{1}{4\pi\epsilon_0}\) assumes you are working in a particular system of units.
Incorrect answers you're likely to see
Two charges instead of one
\[\cancel{V=\frac{kq_{1}q_{2}}{r}}\]
Distance squared in the denominator
\[\cancel{V=\frac{kq}{r^2}}\]
Possible follow-up questions to help with the disambiguation:
Relationship between potential and potential energy \(U = qV\)
Which function is the derivative of the other: \(1/r\) or \(1/r^2\)?
Which physical quantity (potential or electric field, potential energy or force) is the derivative of the other?
What is the electrostatic potential conceptually?
Which function falls off faster: \(1/r\) or \(1/r^2\)?
What are the dimensions of potential? Units?
Where is the zero of potential?
Wrap-up
This could be a good time to refer to the (correct) expression for the potential as an iconic equation, which will need to be further interpreted (”unpacked”) in particular physical situations. This is where the course is going next.
This SWBQ can also serve to help students learn about recall as a cognitive activity. While parts of the equations that students write may be incorrect, many other parts will be correct. Let the way in which you manage the class discussion model for the students how a professional goes about quickly disambiguating several different choices. And TELL the students that this is what you are doing. Deliberately invoke their metacognition.
Many students may not know that the electrostatic potential that we are talking about in this activity is the same quantity as what a voltmeter reads, in principle, but not in practice. You may need to talk about how a voltmeter actually works, rather than idealizing it. It helps to have a voltmeter with leads as a prop. Students often want to know about the “ground” lead. We often tie a long string to it (to symbolize making a really long wire) and send the TA out of the room with the string, “headed off to infinity” while discussing the importance of setting the zero of potential. The extra minute or two of humerous byplay gives the importance of the zero of potential a chance to sink in.
We use this small whiteboard question as a transition between The Distance Formula (Star Trek) activity, where students are learning about how to describe (algebraically) the geometric distance between two points, and the Electrostatic Potential Due to a Pair of Charges (with Series) activity, where students are using these results and the superposition principle to find the electrostatic potential due to two point charges.
This activity is the initial activity in the sequence Visualizing Scalar Fields addressing the representations of scalar fields in the context of electrostatics.
Found in: Static Fields, None course(s)Found in: Warm-Up, E&M Ring Cycle Sequence sequence(s)
These notes from week 6 of https://paradigms.oregonstate.edu/courses/ph441 cover the ideal gas from a grand canonical standpoint starting with the solutions to a particle in a three-dimensional box. They include a number of small group activities.
Students use a plastic surface representing the potential due to a charged sphere to explore the electrostatic potential, equipotential lines, and the relationship between potential and electric field.
A student is invited to “act out” motion corresponding to a plot of effective potential vs. distance. The student plays the role of the “Earth” while the instructor plays the “Sun”.
Students use a pre-written Mathematica notebook or a Geogebra applet to explore how the shape of the effective potential function changes as the various parameters (angular momentum, force constant, reduced mass) are varied.
Students observe the motion of a puck tethered to the center of the airtable. Then they plot the potential energy for the puck on their small whiteboards. A class discussion follows based on what students have written on their whiteboards.
Students examine a plastic "surface" graph of the electric potential due to two charged plates (near the center of the plates) and explore the properties of the electric potential.
Students examine a plastic "surface" graph of the gravitational potential energy of a Earth-satellite system to make connections between gravitational force and gravitational potential energy.