^{*}

With the advent of Computer Algebra System (CAS) such as Mathematica [1], challenging symbolic longhand calcula-tions can effectively be performed free of error and at ease. Mathematica’s integrated features allow the investigator to combine the needed symbolic, numeric and graphic modules all in one interactive environment. This assists the author to focus on interpreting the output rather than exerting the efforts of relating the scattered separate modules. In this note the author, utilizing these three features, explores the magneto-static field and its associated vector potential of a steady looping current. In particular by deploying the numeric features of Mathematica the exact value of the vector potential of the looping current conducive to its 3D graph is presented.

In the course of quantifying the flux of a magnetodynamic field of a permanent mobile magnet through a certain planar conducting loop, the author needed to evaluate its companion static magnetic field. The evaluation of the latter entails modeling the current distribution within the permanent magnet. One such model deals with a current confined to the rim of a nut-shaped magnet. Other models extensively have been considered as well; for a comprehensive detailed analysis see Sarafian [

This note is composed of three sections. In addition to Introduction and Motivation, in Section 2, Analysis, we extend the current classic theoretical status of evaluating the vector potential and its associated magneto-static field of a steady looping current conducive to its 2D and 3D visual displays. In the last section we close the note with a few remarks.

In general it is assumed that the electric current is the source of the magnetic field. The issue of interest is, “for a given current distribution, to evaluate its associated magnetic field.” Although in general one may apply Biot-Savart law [

The integral is evaluated on the volume of the current distribution denoted with the volume element; with being the observation point vector. The permeability of the free space in SI units is. We apply this formulation to a steady current i looping in a circle of radius a positioned on a horizontal plane centered about the origin of a coordinate system, see

If one chooses a spherical coordinate system, the current density distribution can be written as,

Here d is Dirac delta function, and is the unit vector along the polar direction of the spherical coordinate system denoted with the vector. The denominator of (1) is

where we utilized

and

The volume element in (1) is

Upon substituting these pieces in (1) and integrating over and we arrive at. The f dependent of the potential is traced to

However, utilizing the symmetry of the current yields a f independent potential; that is to say irrespective the polar angle observation point the potential has the same value.

Setting f = 0 yields,

Equation (3) is evaluated analytically in [

In step 2, introduce a new variable. By manipulating the integrand we get,

where. In step 3, we recall the definition of the Complete Elliptical integral of the First and the Second kind, namely,

and

respectively. Then it is a matter of exercise to show these integrals are invariant after replacing with. Hence, the second term of (4) yields 2 K(k). The integrand of the first term of (4) can be also manipulated as and the integration yields. Putting these pieces together we arrive at,

With this analytic expression at hand, we pursue evaluating its associated magnetic field. Applying the fundamental relationship between the vector potential and the field, namely, , in a spherical coordinate system we write,

where the curvilinear scale factors are , [

According to (5) the vector potential is an explicit function of elliptic integrals. Hence, it is not hard to envision that evaluating (7) is mathematically challenging. A thorough search reveals no such calculation has been reported in scientific literature. The classic textbook in electromagnetism, “Classical Electrodynamics”, Jackson [

Utilizing (5 & 7) deploying Mathematica we symbolically evaluate the magnetic field. What follows are the needed Mathematica codes. In these codes EllipticK and EllipticE are the Mathematica library functions for the Complete Elliptic Integrals of the First and the Second kind, respectively.

The radial and the azimuthal components of the magnetic field are denoted by. These are:

These expressions are the components of the magnetic filed at the position vector. For instance, for the given looping current, one intuitively expects that, on the horizontal plane, the radial component of the field to vanish. Since the azimuthal component of the field on the horizontal plane practically is along the vertical axis it should intuitively direct itself in two opposite vertical directions, depending on the point of interest being inside or outside of the loop. The latter observation conforms to the basis for one of the man made “Right Hand Rules”. For numeric evaluation and for simplicity we consider a loop with a one unit radius. Setting the azimuthal angle q = p/2 conforms to the horizontal plane. The next two code lines give the value of the radial component of the field and the plot of the azimuthal B_{q} (vertical B_{z}) field.

Plot[Evaluate[]_{, }{r, 0.1, 2}, PlotStyle –> {Black,Thick}_{,}_{ }AxesLabel –> {“r”, “~B = Bz”}, GridLines ® Automatic]

According to

Graphic display of the potential is missing in the scientific literature; this hinders ones physical insight. Utilizing Mathematica graphics, we fill in the gap. To accomplish this, it is helpful to plot the contours of the potential about the ring. Utilizing the symmetry of the ring we plot the contours only on the xz-plane. First we change the coordinate system from spherical to its associated Cartesian coordinate. This requires replacing x = rsinq. This

also entails reformatting the argument of the elliptical integrals. The code reads,

contourplotExactzoomed = ContourPlot[Afxz/.a®1, {x, 0, 2}, {z, –1, 1}, ContourShading ® False, ContourStyle ® Thick, GridLines ® Automatic, FrameLabel ® {“x”, “z”}]

disk=Graphics[Disk[{1,0},0.02]];

Show[{contourplotExactzoomed,disk}]

This graph assists in forming an opinion about the polar component of the magnetic vector potential. In other words, it shows for the points close to the rim of the ring the potential contours are closed loops asymmetrically wrapped around the ring. For the points near by the symmetry axis, the contours are open; they are almost straight lines parallel to the symmetry axis. For the outside far points the contours are open curves; the further the points are from symmetry axis the curves are less bent.

tableContourN = Table[Show[{ContourPlot[Evaluate [Afxz/.a ® 1] == n, {x, 0, 2}, {z, –1, 1}, ContourShading ® False, ContourStyle ® Thick, GridLines ® Automatic, Frame Label ® {“x”, “z”}](*{n, 2.9, 3.1, 0.1}*), disk}], {n, 3.15, 2.9, –0.05}]

In other words, plots shown in

Applying we may now visualize the magnetic field. The fields are tangent to the potential curves, i.e. the contours and the fields are on the same plane. The fields are plotted applying Mathematica StreamPlot. In doing so first we replace the spherical coordinate system with its associated Cartesian coordinate. We then zoom in three different areas: the entire region, the region close to the center and the region far from the loop. These are shown sequentially in

These plots show the relationship between the potential and the field. For instance the first graph shows how the field in order to sustain its tangential alignment reorients itself; compare this feature vs. the field for the points close to the rim and the ones close the center. One readily may observe useful features analyzing the second and the third graphs as well.

We further our analysis by displaying three 3D graphs for the vector potential

tableSpherical3DN = Table[SphericalPlot3D[Afrq/.a ® 1.,{q, 0, p}, {f, 0, p}, Mesh ® False, PlotStyle ® Directive[Hue[0.2r], Opacity[0.7], Specularity[White, 10], ImageSize®600,PlotRange ® {{–1, 1}, {–1, 1}, {–5, 5}}]], {r, 0.01, 3, 0.5}]

In many words, the magnetic vector potential is a 3D function. For the best scenario its 2D profile utilizing Mathematica contour plot is depicted in Figures 3 and 4. Here in

In

Traditional methodology of mathematical physics with the advent of Computer Algebra System (CAS) has gone through a quantum transition. CAS can assist and minimize the tedious longhand symbolic calculations requiring time consuming efforts that once hindered the exploration of new horizons. Numeric computation of issues of interest as well as their accompanied graphs utilizing CAS is a norm of mathematical physics. Consequently, the authors of scientific academic textbooks are encouraged to utilize CAS to present results. An example of how these features are culminated conducive to a comprehensive research project is presented in this note. The author systematically following the traditional steps of mathematical physics has tackled a topic of interest utilizing a CAS such as Mathematica bringing the project to fruition. The details discussed in this project are an example demonstrating how Mathematica plays an essential role exploring new features not yet reported in scientific literature. The entire project including the text is completed utilizing the latest version of Mathematica, V 8.0.1 [

The author would like to thank Mrs. Nenette Sarafian Hickey for carefully reading over the manuscript and making editorial comments.