Friday, February 22, 2013

The Response of a Spherical Heart to a Uniform Electric Field

In Chapter 7 of the 4th edition of Intermediate Physics for Medicine and Biology, Russ Hobbie and I discuss the bidomain model of cardiac tissue.
Myocardial cells are typically about 10 μm in diameter and 100 μm long. They have the added complication that they are connected to one another by gap junctions, as shown schematically in Fig. 7.27. This allows currents to flow directly from one cell to another without flowing in the extracellular medium. The bidomain (two-domain) model is often used to model this situation [Henriquez (1993)]. It considers a region, small compared to the size of the heart, that contains many cells and their surrounding extracellular fluid.
The citation is to the 20-year-old-but-still-useful review article by Craig Henriquez of Duke University.
Henriquez, C. S. (1993) “Simulating the electrical behavior of cardiac tissue using the bidomain model,” Crit. Rev. Biomed. Eng., Volume 21, Pages 1–77.
According to Google Scholar, this landmark paper has been cited over 450 times (including a citation on page 202 of IPMB).

During the early 1990s I collaborated with another researcher from Duke, Natalia Trayanova. Our goal was to apply the bidomain model to the study of defibrillation of the heart. In the same year that Craig’s review appeared, Trayanova, her student Lisa Malden, and I published an article in the IEEE Transactions on Biomedical Engineering titled “The Response of the Spherical Heart to a Uniform Electric Field: A Bidomain Analysis of Cardiac Stimulation” (Volume 40, Pages 899–908). I’m fond of this paper for several reasons:
  • Like most physicists, I like simple models that highlight and clarify basic mechanisms. Our spherical heart model had that simplicity.
  • The article was the first to show that fiber curvature provides a mechanism for polarization of cardiac tissue in response to an electrical shock. Since our paper, researchers have appreciated the importance of the fiber geometry in the heart when modeling electrical stimulation.
  • The model emphasizes the role of unequal anisotropy ratios in the bidomain model. In cardiac tissue, both the intracellular and extracellular spaces are anisotropic (the electrical conductivity is different parallel to the myocardial fibers then perpendicular to them), but the intracellular space is more anisotropic than the extracellular space. Fiber curvature will only result in polarization deep in the heart wall if the tissue has unequal anisotropy ratios.
  • The calculation has important clinical implications. Fibrillation of the heart is a leading cause of death in the United States, and the only way to treat a fibrillating heart is to apply a strong electric shock: defibrillation. I’ve performed a lot of numerical simulations in my career, but none have the potential impact for medicine as my work on defibrillation.
  • The IEEE TBME publishes brief bios of the authors. Back in those days I published in this journal often, and my goal was to have my entire CV included, bit by bit, in these small bios. The one in this paper read “Bradley J Roth was raised in Morrison, Illinois. He received the BS degree in physics from the University of Kansas in 1982, and the PhD in physics from Vanderbilt University in 1987. His PhD dissertation research was performed in the Living State Physics Laboratory under the direction of Dr. J. WIkswo. He is now a Senior Staff Fellow with the Biomedical Engineering and Instrumentation Program, National Institutes of Health, Bethesda, MD. One of this research interests is the mathematical modeling of the electrical behavior of the heart. He is also interested in the production and interactions of magnetic fields with biological tissue, e.g. biomagnetism, magnetic stimulation, and magnetoacoustic imaging.”
  • The acknowledgments state “the authors thank B. Bowman for his assistance in editing the manuscript.” Barry was a great help to me in improving my writing skills during my years at NIH, and I’m glad that we mentioned him.
  • The paper cites several of my favorite books, including When Time Breaks Down by Art Winfree, Classical Electrodynamics by John David Jackson, and Handbook of Mathematical Functions with Formulas, Graphs, and and Mathematical Tables, by Abramowitz and Stegun.
  • The paper has been fairly influential. It’s been cited 97 times, which is small potatoes compared to Henriquez’s review, but not too shabby nevertheless; an average of almost five citations a year for 20 years.
  • It was a pleasure to collaborate with Natalia Trayanova, who I was to work with again seven years later on another study of cardiac electrical behavior (Lindblom, Roth, and Trayanova, Journal of Cardiovascular Electrophysiology, Volume 11, Pages 274–285, 2000).
  • The paper led to subsequent simulations of defibrillation that are much more realistic and sophisticated than our simple spherical model of twenty years ago. Trayanova has led the way in this research, first at Duke, then at Tulane, and now at Johns Hopkins. You can listen to her discuss her research here. If you have a subscription to the Journal of Visualized Experiments you can hear more here. For a recent review, see Trayanova et al. (2012). Also, see this article recently put out by Johns Hopkins University. 
Listen to Natalia Trayanova discuss developing computer simulations to improve arrhythmia treatments.
Cardiac Bioelectric Therapy: Mechanisms and Practical Implications with Intermediate Physics for Medicine and Biology.
Cardiac Bioelectric Therapy:
Mechanisms and Practical Implications.
To learn more about how physics and engineering can help us understand defibrillation, consult the book Cardiac Bioelectric Therapy: Mechanisms and Practical Implications, which has chapters by Trayanova and many of the other leading researchers in the field (including yours truly).

Friday, February 15, 2013

The Joy of X

The Joy of X,  by Steven Strogatz, superimposed on Intermediate Physics for Medicine and Biology.
The Joy of X,
by Steven Strogatz.
Steven Strogatz’s latest book is The Joy of X: A Guided Tour of Math, From One to Infinity. I have discussed books by Strogatz in previous entries of this blog, here and here. The preface defines the purpose of The Joy of X.
The Joy of X is an introduction to math’s most compelling and far-reaching ideas. The chapters—some from the original Times series [a series of articles about math that Strogatz wrote for the New York Times]—are bite-size and largely independent, so feel free to snack wherever you like. If you want to wade deeper into anything, the notes at the end of the book provide additional details, and suggestions for further reading.
My favorite chapter in The Joy of X was “Twist and Shout” about Mobius strips. Strogatz’s discussion was fine, but what I really enjoyed was the lovely video he called my attention to: “Wind and Mr. Ug”. Go watch it right now; it’s less than 8 minutes long. It is the most endearing mathematical story since Flatland.

Wind and Mr. Ug.

Of course, I’m always on the lookout for medical and biological physics, and I found it in Strogatz’s chapter called “Analyze This!,” in which he describes the Gibbs phenomenon. I have written about the Gibbs phenomenon in this blog before, but not so eloquently. Russ Hobbie and I introduce the Gibbs phenomenon in Chapter 11 of the 4th edition of Intermediate Physics for Medicine and Biology. When talking about the fit of a Fourier series to a square wave, we write
As the number of terms in the fit is increased, the value of Q [a measure of the goodness of the fit] decreases. However, spikes of constant height (about 18% of the amplitude of the square wave or 9% of the discontinuity in y) remain…These spikes appear whenever there is a discontinuity in y and are called the Gibbs phenomenon.
It turns out that the Gibbs phenomenon is related to the alternating harmonic series. Strogatz writes
Consider this series, known in the trade as the alternating harmonic series:
1 – 1/2 + 1/3 – 1/4 + 1/5 – 1/6 + … .
[…] The partial sums in this case are
S1 = 1
S2 = 1 – 1/2 = 0.500
S3 = 1 – 1/2 + 1/3 = 0.833 …
S4 = 1 – 1/2 + 1/3 – 1/4 = 0.583…

And if you go far enough, you’ll find that they home in on a number close to 0.69. The series can be proven to converge. Its limiting value is the natural logarithm of 2, denoted ln2 and approximately equal to 0.693147. […]

Let’s look at a particularly simple rearrangement whose sum is easy to calculate. Supposed we add two of the negative terms in the alternating harmonic series for every one of its positive terms, as follows:

[1 – 1/2 – 1/4] + [1/3 – 1/6 – 1/8] + [1/5 – 1/10 – 1/12] + …

Next, simplify each of the bracketed expressions by subtracting the second term from the first while leaving the third term untouched. Then the series reduces to

[1/2 – 1/4] + [1/6 – 1/8] + [1/10 – 1/12] + …

After factoring out ½ from all the fractions above and collecting terms, this becomes

½ [ 1 – 1/2 + 1/3 – 1/4 + 1/5 – 1/6 + …].

Look who’s back: the beast inside the brackets is the alternating harmonic series itself. By rearranging it, we’ve somehow made it half as big as it was originally—even though it contains all the same terms!”
Strogatz then relates this to a Fourier series

f(x) = sinx – 1/2 sin 2x + 1/3 sin 3x – 1/4 sin 4x + …

This series approaches a sawtooth curve. But when he examines its behavior with different numbers of terms in the sum, he finds the Gibbs phenomenon.
Something goes wrong near the edges of the teeth. The sine waves overshot the mark there and produce a strange finger that isn’t in the sawtooth wave itself… The blame can be laid at the doorstep of the alternating harmonic series. Its pathologies discussed earlier now contaminate the associated Fourier series. They’re responsible for that annoying finger that just won’t go away.
In the notes about the Gibbs phenomenon at the end of the book, Strogatz points us to a fascinating paper on the history of this topic
Hewitt, E. and Hewitt, R. E. (1979) “The Gibbs-Wilbraham Phenomenon: An Episode in Fourier Analysis,” Archive for the History of Exact Sciences, Volume 21, Pages129–160.
He concludes his chapter
This effect, commonly called the Gibbs phenomenon, is more than a mathematical curiosity. Known since the mid-1800s, it now turns up in our digital photographs and on MRI scans. The unwanted oscillations caused by the Gibbs phenomenon can produce blurring, shimmering, and other artifacts at sharp edges in the image. In a medical context, these can be mistaken for damaged tissue, or they can obscure lesions that are actually present.

Friday, February 8, 2013

Photodynamic Therapy

I am currently teaching Medical Physics (PHY 326) at Oakland University, and for our textbook I am using (surprise!) the 4th edition of Intermediate Physics for Medicine and Biology. In class, we recently finished Chapter 14 on Atoms and Light, which “describes some of the biologically important properties of infrared, visible, and ultraviolet light.”

Once a week, class ends with a brief discussion of a recent Point/Counterpoint article from the journal Medical Physics (see here and here for my previous discussion of Point/Counterpoint articles). I find these articles to be useful for introducing students to cutting-edge questions in modern medical physics. The title of each article contains a proposition that two leading medical physicists debate, one for it and one against it. This week, we discussed an article about photodynamic therapy (PDT) by Timothy C. Zhu (University of Pennsylvania, for the proposition) and E. Ishmael Parsai (University of Toledo, against the proposition):
Zhu, T. C., and E. I Parsai (2011) “PDT is Better than Alternative Therapies Such as Brachytherapy, Electron Beams, or Low-Energy X Rays for the Treatment of Skin Cancers,” Medical Physics, Volume 38, Pages 1133–1135.
When reading through the article, I thought I would check how extensively we discuss of PDT in IPMB. I found that we say nothing about it! A search for the term “photodynamic” or “PDT” comes back empty. So, this week (with an eye toward the 5th edition) I am preparing a very short new section in Chapter 14 about PDT.
14.8 ½ Photodynamic Therapy

Photodynamic therapy (PDT) uses a drug called a photosensitizer that is activated by light [Zhu and Finlay (2008), Wilson and Patterson (2008)]. PDT can treat accessible solid tumors such as basal cell carcinoma, a type of skin cancer [see Sec. 14.9.4]. An example of PDT is the surface application of 5-aminolevulinic acid, which is absorbed by the tumor cells and is transformed metabolically into the photosensitizer protoporphyrin IX. When this molecule interacts with light in the 600-800 nm range (red and near infrared), often delivered with a diode laser, it converts molecular oxygen into a highly reactive singlet state that causes necrosis, apoptosis (programmed cell death), or damage to the vasculature that can make the tumor ischemic. Some internal tumors can be treated using light carried by optical fibers introduced through an endoscope.
The two citations are to the articles 
Wilson, B. C. and M. S. Patterson (2008) “The Physics, Biophysics and Technology of Photodynamic Therapy,” Physics in Medicine and Biology, Volume 53, Pages R61–R109.

Zhu, T. C. and J. C. Finlay (2008) “The Role of Photodynamic Therapy (PDT) Physics,” Medical Physics, Volume 35, Pages 3127–3136.
The first PhD dissertation from the Oakland University Medical Physics graduate program dealt with photodynamic therapy: In Vivo Experimental Investigation on the Interaction Between Photodynamic Therapy and Hyperthermia, by James Mattiello (1987).

You can learn more about photodynamic therapy here and here. Please don’t confuse PDT with the alternative medicine (bogus) treatment “Sono Photo Dynamic Therapy.”

Friday, February 1, 2013

The Page 99 Test

English editor Ford Madox Ford advised people who are debating if they should read a particular book to “open the book to page ninety-nine and read, and the quality of the whole will be revealed to you.” This approach is now called the Page 99 Test. Although arbitrary, it provides a way to decide quickly if a book will interest you. Let’s try the Page 99 Test with the 4th edition of Intermediate Physics for Medicine and Biology. Section 4.12 comparing drift and diffusion ends on Page 99, and Section 4.13 about the solution to the diffusion equation begins. The page contains five displayed equations (four of them numbered, Eqs. 4.70 to 4.73) and three figures (Figs. 4.17 to 4.19). An example of the text of page 99 is the opening paragraph at the start of Sec. 4.13.
If C(x, 0) is known for t = 0, it is possible to use the result of Sec. 4.8 to determine C(x,t) at any later time. The key to doing this is that if C(x,t) dx is the number of particles in the region between x and x+dx at time t, it may be be interpreted as the probability of finding a particle in the interval (x, dx) multiplied by the total number of particles. (Recall the discussion on p. 91 about the interpretation of C(x,t).) The spreading Gaussian then represents the spread of probability that a particle is between x and x + dx.
Page 99 appears in the Table of Contents:
4.13 A General Solution for the Particle Concentration as a Function of Time . . . . . 99
and the title of this section appears as the running title at the top of the page. Page 99 appears three times in the index, under 1) Diffusion equation, general solution, 2) Fick’s law (frankly, I'm not sure why page 99 is listed for Fick’s law, as I don't see it mentioned explicitly anywhere on that page), and 3) Gaussian distribution. According to the Symbol List at the end of Chapter 4, the first use of the Greek symbol xi for position was on page 99. Somewhat unusually, no references are cited on page 99 (there are citations on the page before and the page after). No corrections to page 99 appear in the errata, and no words are emphasized using italics.

Does Intermediate Physics for Medicine and Biology pass the page 99 test? I think so. The topic—diffusion—is a physical phenomenon that is crucial for understanding biology. The mix of equations and figures is similar to the remainder of the book. Calculus is used without apology. If you like page 99, I think you will enjoy the rest of the book. And if you like page 99, you are going to love page 100, which contains more equations and figures, plus error functions, Green’s functions, random walks, and citations to classic texts such as Benedek and Villars (2000), Carslaw and Jaeger (1959), and Crank (1975). And if you liked page 100, on page 101 you find......

Friday, January 25, 2013

Aliasing

In Chapter 11 of the 4th edition of Intermediate Physics for Medicine and Biology, Russ Hobbie and I discuss aliasing.
If a component [in the Fourier spectrum] is present whose frequency is more than half the sampling frequency, it will appear in the analysis at a lower frequency. This is the familiar stroboscopic effect in which the wheels of the stagecoach appear to rotate backward because the samples (movie frames) are not made rapidly enough. In signal analysis, this is called aliasing. It can be seen in Fig. 11.15, which shows a sine wave sampled at regularly spaced intervals that are longer than half a period.
First of all, what is all this business about a stagecoach? Fifty years ago, when westerns were all the rage in movies and on TV, aliasing often occurred if the frame rate (typically 24 frames per second for old movies) was lower than the rotation rate of the wheel (if all the spokes of the wheel are equivalent, then you can take the “period of rotation” as the time it takes for one spoke to rotate to the position of the adjacent one, which may be much shorter than the time for the wheel to make one complete rotation). You can see an example of this in the John Wayne movie Winds of the Wasteland (1936), especially in the climactic scene of the stagecoach race. In this video of the movie, you can see aliasing of the stagecoach wheel briefly at time 55:40. For those of you who are more discriminating in your movie tastes, you can see another example of aliasing 14 minutes and 15 seconds into Stagecoach, a John Wayne classic from 1939 directed by John Ford. In my opinion, the greatest western is the John Ford masterpiece The Man Who Shot Liberty Valance. What more could you ask for than both John Wayne and Jimmy Stewart in the same production? You can see aliasing briefly when Stewart drives his buckboard out of Shinbone to practice his pistol shooting (without much success). Another time when you see a wheel rotate backwards in this movie does not involve aliasing; it is (Spoiler Alert!) after Stewart Wayne kills Valance (Lee Marvin), when Pompey (Woody Strode) takes the drunken Wayne to his ranch house where he backs up the buckboard (that was a joke….).

But I digress. Aliasing can happen in space as well as time, and can therefore affect images. If spatial frequencies in the structure of an object correspond to wavelengths smaller than the twice the pixel size, low spatial frequency artifacts, such as Moire patterns, can appear in the image, shown nicely in this figure. One can minimize aliasing by first filtering (anti-aliasing) before sampling. Some rather extreme cases of aliasing can been seen in Figs. 11.41 and 12.11 of Intermediate Physics for Medicine and Biology.

 Stagecoach, with John Wayne.

Friday, January 18, 2013

The Magic Angle

I recently found another error in the 4th edition of Intermediate Physics for Medicine and Biology. In Chapter 18 about magnetic resonance imaging, Homework Problem 18 reads
Problem 18  Suppose the two dipoles of the water molecule shown below point in the z direction while the line between them makes an angle θ with the x axis. Determine the angle θ for which the magnetic field of one dipole is perpendicular to the dipole moment of the other. For this angle the interaction energy is zero. This θ is called the “magic angle” and is used when studying anisotropic tissue such as cartilage [Xia (1998)].
Technically there is nothing wrong with this problem. However, if I were doing it over I would have the angle θ measured from the z axis, not the x axis. One reason is that this is the way θ is defined most often in the literature. Another is that in the solution manual we solve the problem as if θ were relative to the z axis, so the book and the solution manual are not consistent on the definition of θ. I should add, this problem was not present in the 3rd edition of Intermediate Physics for Medicine and Biology. It is a new problem I wrote for the 4th edition, so I can’t blame Russ Hobbie for this one (rats).

The citation in the homework problem is to the paper
Xia, Y. (1998) “Relaxation Anisotropy in Cartilage by NMR Microscopy (μMRI) at 14-μm Resolution,” Magnetic Resonance in Medicine, Volume 39, Pages 941–949.
The author, Yang Xia, is a good friend of mine, and a colleague here in the Department of Physics at Oakland University. He is well-known around OU because over the last decade he had the most grant money from the National Institutes of Health of anyone on campus. He uses a variety of techniques, including micro-magnetic resonance imaging (μMRI), to study cartilage and osteoarthritis. The abstract of his highly-cited paper reads
To study the structural anisotropy and the magic-angle effect in articular cartilage, T1 and T2 images were constructed at a series of orientations of cartilage specimens in the magnetic field by using NMR microscopy (μMRI). An isotropic T1, and a strong anisotropic T2 were observed across the cartilage tissue thickness. Three distinct regions in the microscopic MR images corresponded approximately to the superficial, transitional, and radial histological zones in the cartilage. The percentage decrease of T2 follows the pattern of the curve of (3cos2θ - 1)2 at the radial zone, where the collagen fibrils are perpendicular to the articular surface. In contrast, little orientational dependence of T2 was observed at the transitional zone, where the collagen fibrils are more randomly oriented. The result suggests that the interactions between water molecules and proteoglycans have a directional nature, which is somehow influenced by collagen fibril orientation. Hence, T2 anisotropy could serve as a sensitive and noninvasive marker for molecular-level orientations in articular cartilage.
Perhaps a better reference for our homework problem is another paper of Xia’s.
Xia, Y. (2000) “Magic Angle Effect in MRI of Articular Cartilage: A Review,” Investigative Radiology, Volume 35, Pages 602–621.
There in Fig. 3 of Xia’s review is a picture almost identical to the figure that immediately follows Homework Problem 18 in our book, except the angle θ is measured from the direction of the static magnetic field rather than perpendicular to it. Xia writes
T2 corresponds to the decay in phase coherence (dephasing) between the individual nuclear spins in a sample (protons in our case). Because each proton has a magnetic moment, it generates a small local dipolar magnetic field that impinges on its neighbor’s space (Fig. 3).43 This local field fluctuates constantly because the molecule is tumbling randomly. The T2 process can occur under the influence of this fluctuating magnetic field. At the end of signal excitation during an MRI experiment, the net magnetization (which produces the MRI signal) is coherent and points along a certain direction in space in the rotating frame of reference. This coherent magnetization vector soon becomes dephased because the local magnetic fields associated with the magnetic properties of neighboring nuclei cause the precessing nuclei to acquire a range of slightly different precessional frequencies. The time scale of this signal dephasing is reported as T2 and is characteristic of the molecular environment in the sample.43,44

For simple liquids or samples containing simple liquidlike molecules, the molecules tumble rapidly. The dipolar spin Hamiltonian (HD) that describes the dipolar interaction is averaged to zero, and its contribution to the spin relaxation vanishes. Relaxation characteristics exhibit a simple exponential decay that is well described by the Bloch equations.45 For samples containing molecules that are less mobile, HD is no longer averaged to zero and makes a significant contribution to the relaxation, resulting in a shorter T2. When HD is not zero, it is dominated by a geometric factor, (3cos2θ - 1), where θ is the angle between the position vector joining the two spins and the external magnetic field (see Fig. 3). A useful feature of this geometric factor is that it approaches zero as θ approaches 54.74° (Fig. 4). Therefore, even when HD is not zero, the contribution of HD to spin relaxation can be minimized if θ is set to 54.74°. This angle is called the magic angle in NMR.46
So, in the errata you will now find this:
Page 539: In Chapter 18, Homework Problem 18, “while the line between them makes an angle θ with the x axis” should be “while the line between them makes an angle θ with the z axis”. Also, in the accompanying figure following the homework problem, the angle θ should be measured from the z (vertical) axis, not the x (horizontal) axis. Corrected 1-18-13.
Is this the last error that we’ll find in our book? I doubt it; there are sure to be more we haven’t found yet. If you find any, please let us know.

Friday, January 11, 2013

5th Edition of Intermediate Physics for Medicine and Biology

Russ Hobbie and I are starting to talk about a 5th edition of Intermediate Physics for Medicine and Biology, and we need your help. We would like suggestions and advice about what changes/additions/deletions to make in the new edition.

We have prepared a survey to send to faculty members who we know have used IPMB as the textbook for a class they taught. However, our list may be incomplete, and input from any teacher, student, or reader would be useful. So, below is a copy of the survey. Please send responses to any or all of the questions to Russ (hobbie@umn.edu).

Thanks!
  1. What chapters did you cover when teaching from the 4th edition of IPMB?
  2. Were the homework problems appropriate?
  3. In the 4th edition we added a chapter on Sound and Ultrasound to IPMB. If we were to add one new chapter to the 5th edition, what should the topic be?
  4. Would color significantly improve the book for your purposes? How much extra money would you be willing to pay if the 5th edition contained many color pictures?
  5. What is the best feature of IPMB? What is the worst?
  6. Is the end-of-chapter list of symbols useful?
  7. Do your students use the Appendices? Suppose to save space one Appendix had to be deleted: which one should go?
  8. Did you have access to the solution manual? Was it useful? The solution manual we prepared using different software than the book itself. Did you see a noticeable difference in the quality of the book and the solution manual?
  9. Would you like students to have access to the solution manual?
  10. Did you use any information on the book website, such as the errata or text from previous editions?
  11. Are you aware of the book blog? Did you find it useful when teaching from IPMB? Do you find it interesting?
  12. How important is having a paperback version of the book?
  13. What textbooks did you consider other than IPMB? If IPMB did not exist, what book would you use for your class?

Friday, January 4, 2013

Non-Dynamical Stochastic Resonance: Theory and Experiments with White and Arbitrarily Coloured Noise

Section 11.18 of the 4th edition of Intermediate Physics for Medicine and Biology contains a discussion of stochastic resonance. This is a new section that Russ Hobbie and I added to the 4th edition, and features a discussion of a paper by Zoltan Gingl, Laszlo Kish (formerly “Kiss”), and Frank Moss.
Gingl, Z., L. B. Kiss, and F. Moss (1995) “Non-Dynamical Stochastic Resonance: Theory and Experiments with White and Arbitrarily Coloured Noise,” Europhysics Letters, Volume 29, Pages 191–196.
The paper is interesting (despite the annoying British spelling), and I reproduce part of the introduction below.
In the last decade’s physics literature, stochastic-resonance (SR) effect has been one of the most interesting phenomena taking place in noisy non-linear dynamical systems (see, e.g., [l-14]. The input of stochastic resonators [12] (non-linear systems showing SR) is fed by a Gaussian noise and a sinusoidal signal with frequency f0, that is, a random excitation and a periodic one are acting on the system. There is an optimal strength of the input noise, such that the system’s output power spectral density, at the signal frequency f0, has a maximal value. This effect is called SR. It can be viewed as: the transfer of the input sinusoidal signal through the system shows a “resonance vs. the strength of the input noise. It is a very interesting, and somewhat paradoxial effect, because it indicates that in these systems the existence of a certain amount of “indeterministic excitation is necessary to obtain the optimal “deterministic response. There are certain indications [2,13,14] that the principle of SR may be applied by nature in biological systems in order to optimise the transfer of neural signals.
Until last year, it was a common belief that SR phenomena occur only in (bistable, sometimes monostable [10] or multistable) dynamical systems [1-14]. Very recently, Wiesenfeld et al. [15] have proposed that certain systems with threshold-like properties should also show SR effects.
We present here an extremely simple system, invented by Moss, which displays SR. It consists only of a threshold and a subthreshold coherent signal plus noise as shown in fig. la). It is not a dynamical system, instead there is a single rule: whenever the signal plus the noise crosses the threshold unidirectionally, a narrow pulse of standard shape is written to a time series, as shown in fig. lb). The power spectrum of this series of pulses is shown in fig. 1c). It shows all the familiar features of SR systems previously studied [l, 2, 7, 16], in particular, the narrow, delta-like signal features riding on a broad-band noise background from which the signal-to-noise ratio (SNR) can be extracted. This system can be easily realized electronically as a level-crossing detector (LCD). There is a simple and very physically motivated theory of this phenomenon (due to Kiss), see below. Other, more detailed studies of various aspects of threshold-crossing dynamics have been made by Fox et al. [17], Jung [18] and Bulsara et al. [19].
We have experimentally realised and developed this simple SR system and carried out extensive analog and computer simulations on it. The theory of Kiss has been verified for the case of white and several sorts of coloured noises. Until now, the description of this new SR system, its physical realisation and the original theory have not appeared in the open literature, so in this letter we shall describe the new system and its developments made by us, present the outline and the main results of the theory and finally show some interesting experimental results…
Figure 1 in their paper is our Figure 11.50. It is an excellent figure, although I don’t know why they didn’t adjust the time axes so that the pulses in b) are aligned precisely with the signal crossings in a). The axes are almost correct, but are off just enough to be confusing, like when the video and audio signals are off by a fraction of a second in a movie or TV show.

The Gingl et al. paper is short and highly cited (over 200 citations to date, according to the Web of Science). However, it is not cited nearly as often as another paper published by Kurt Wisenfeld and Moss that same year:
Wisenfeld K. and F. Moss (1995) “Stochastic Resonance and the Benefits of Noise: From Ice Ages to Crayfish and SQUIDs,” Nature, Volume 373, Pages 33–36.
This paper, with over 1000 citations, reviews many applications of stochastic resonance.
Noise in dynamical systems is usually considered a nuisance. But in certain nonlinear systems, including electronic circuits and biological sensory apparatus, the presence of noise can in fact enhance the detection of weak signals. This phenomenon, called stochastic resonance, may find useful application in physical, technological and biomedical contexts.
Wisenfeld and Moss discuss how the crayfish may use stochastic resonance to detect weak signals with their mechanoreceptor hair cells.

Frank Moss (1934-2011) was the founding director of the Center for Neurodynamics at the University of Missouri at St Louis. Click here to read his obituary (he died two years ago today) in Physics Today, and click here to read a tribute to him in a focus issue of the journal Chaos.

Friday, December 28, 2012

Determining the Site of Stimulation During Magnetic Stimulation of a Peripheral Nerve

As December draws to a close and I reflect on all that’s happened over the last twelve months, I conclude that 2012 has been a good year. For me, it has also marked some important anniversaries. Thirty years ago (1982) I graduated from the University of Kansas with a bachelors degree in physics. Twenty-five years ago (1987) I obtained my PhD from Vanderbilt University. And twenty years ago (1992) I was at the National Institutes of Health in Bethesda, Maryland working on magnetic stimulation of nerves.

Today I want to focus on one particular paper published in 1992 that examined magnetic stimulation of a peripheral nerve: “Determining the Site of Stimulation During Magnetic Stimulation of a Peripheral Nerve” (Electroencephalography and Clinical Neurophysiology, Volume 85, Pages 253–264). To understand this article, we must first examine Frank Rattay’s analysis of electrical stimulation. Rattay showed that excitation along a nerve axon occurs where the “activating function” –λ2 d2Ve/dx2 is largest, with λ the length constant, Ve the extracellular potential produced by a stimulating electrode, and x the distance along the axon. Homework Problem 38 in Chapter 7 of the 4th edition of Intermediate Physics for Medicine and Biology guides you through Rattay’s derivation. In 1990, Peter Basser and I showed that this result also holds during magnetic stimulation. What is magnetic stimulation? In Chapter 8 of Intermediate Physics for Medicine and Biology, Russ Hobbie and I discuss Faraday’s law of induction, and then write
Since a changing magnetic field generates an induced electric field, it is possible to stimulate a nerve or muscle cells with out using electrodes…One of the earliest investigations was reported by Barker et al. (1985) who used a solenoid in which the magnetic field changed by 2 T in 110 μs to apply a stimulus to different points on a subject’s arm and skull.
The main difference between Rattay’s analysis of electrical stimulation and our analysis of magnetic stimulation was that Rattay expressed his activation function in terms of the electric potential produced by the stimulus electrode, whereas Basser and I considered the induced electric field along the axon, Ex, and wrote the activating function as λ2 dEx/dx. The most interesting feature of this result is that stimulation does not occur where the electric field is strongest, but instead where its gradient along the axon is greatest. In the early 1990s, this result was surprising (in retrospect, it seems obvious), so we set out to test it experimentally.

Basser and I both worked in NIH’s Biomedical Engineering and Instrumentation Program, and we had neither the expertise nor facilities to perform the needed experiments, but we knew who did. Since arriving at NIH in 1988, I had been working with Mark Hallett and Leo Cohen to develop clinical applications of magnetic stimulation. Also collaborating with Hallett was a delightful couple visiting from Italy, Jan Nilsson and his wife Marcela Panizza. Under Hallett’s overall leadership, with Nilsson and Panizza making the measurements, and with me occasionally making suggestions and cheerleading, we carried out the key experiments that confirmed Basser’s and my prediction about where excitation occurs. These studies were performed on human volunteers (at that time there were people who made their living as paid normal volunteers in clinical studies at NIH) and were carried out in the NIH clinical center. The abstract of our now 20-year old paper said
Magnetic stimulation has not been routinely used for studies of peripheral nerve conduction primarily because of uncertainty about the location of the stimulation site. We performed several experiments to locate the site of nerve stimulation. Uniform latency shifts, similar to those that can be obtained during electrical stimulation, were observed when a magnetic coil was moved along the median nerve in the region of the elbow, thereby ensuring that the properties of the nerve and surrounding volume conductor were uniform. By evoking muscle responses both electrically and magnetically and matching their latencies, amplitudes and shapes, the site of stimulation was determined to be 3.0 ± 0.5 cm from the center of an 8-shaped coil toward the coil handle. When the polarity of the current was reversed by rotating the coil, the latency of the evoked response shifted by 0.65 ± 0.05 msec, which implies that the site of stimulation was displaced 4.1 ± 0.5 cm. Additional evidence of cathode- and anode-like behavior during magnetic stimulation comes from observations of preferential activation of motor responses over H-reflexes with stimulation of a distal site, and of preferential activation of H-reflexes over motor responses with stimulation of a proximal site. Analogous behavior is observed with electrical stimulation. These experiments were motivated by, and are qualitatively consistent with, a mathematical model of magnetic stimulation of an axon.
Rather than describe this experiment in detail, I will let you analyze it yourself in a new homework problem (your three-days-late Christmas present). It is similar to a problem from an exam I gave to my biological physics (PHY 325) students.
Section 8.7

Problem 26 ½ (a) Rederive the cable equation for the transmembrane potential v (Eq. 6.55) using one crucial modification: generalize Eq. 6.48 to account for part of the intracellular electric field that arises from Faraday induction and therefore cannot be written as the gradient of a potential,
 Assume you measure v relative to the resting potential so Eq. 6.53 becomes jm = gm v, and let the extracellular potential be small so vi = v. Identify the new source term in the cable equation (the “activating function” for magnetic stimulation), analogous to vr in Eq. 6.55.
(b) Let
Calculate the activating function and plot both the electric field and the activating function versus x.
(c) Suppose you stimulate a nerve using this activating function, first with one polarity of the current pulse and then the other. What additional delay in the response of the nerve (as measured by the arrival time of the action potential at the far end) will changing polarity cause because of the extra distance the action potential must travel? Assume a = 4 cm and the conduction speed is 60 m/s.
At about the same time as we were doing this study, Paul Maccabee and his colleagues at the SUNY Health Science Center in Brooklyn were carrying out similar experiments using an in-vitro pig nerve model (a nerve in a dish), and came to similar conclusions (“Magnetic Coil Stimulation of Straight and Bent Amphibian and Mammalian Peripheral Nerve In Vitro: Locus of Excitation,” Journal of Physiology, Volume 460, Pages 201–219, 1993). Our paper was published first (Yes!!!) but their results were cleaner and more elegant, in part because they didn’t have the complication of the nerve being surrounded by irregularly shaped muscles and bones. Our paper has been fairly influential (53 citations to date in the Web of Science), but theirs has had an even greater impact (147 citations). A year later Maccabee and I together published a study of a new magnetic stimulation coil design.

What has happened to this cast of characters in the last 20 years? Hallett and Cohen remain at NIH, still doing great work. Nilsson is a biomedical engineer and Panizza is a neurophysiologist in Italy. Basser is at NIH, but is now with the Eunice Kennedy Shriver National Institute of Child Health and Human Development, where he works on MRI diffusion tensor imaging. Paul Maccabee is a neurologist and the Director of the EMG Laboratory at SUNY Brooklyn. I left NIH in 1995, and am now at Oakland University, where I teach, do research, and write a blog so I can wish readers of Intermediate Physics for Medicine and Biology a Happy New Year!

Friday, December 21, 2012

Royal Institution Christmas Lectures

With Christmas approaching, my attention naturally turns to the Royal Institution Christmas Lectures. The Royal Institution (Ri) website states
The Ri is an independent charity dedicated to connecting people with the world of science. We’re about discovery, innovation, inspiration and imagination. You can explore over 200 years of history making science in our Faraday Museum as well as engage with the latest research, ideas and debates in our public science events.

We run science programmes for young people at our Young Scientist Centre, present exciting, demonstration-packed events for schools and run mathematics masterclasses across the UK.

We are most famous for our Christmas Lectures which were started by Michael Faraday in 1825. Check out the 2011 Lectures here and don’t miss them this Christmas on BBC Four.

Anyone can join the Ri. If you’re interested in how the world works, or how to make it work better through science, the Ri is the place for you.
The 2012 Christmas Lectures, “The Modern Alchemist,” will be broadcast on BBC Four on December 26, 27, and 28 at 8pm. Don’t get BBC Four? Neither do I. But that’s OK, because you can watch the Christmas Lectures at the Ri website. In fact, you can watch the Christmas Lectures from past years too. You will have to open an account, which means you will need to give them your email address and other information, but you don’t need to pay anything; it’s free. Kind of like a Christmas present.

My favorite lecture is from 2010. Mark Miodownik stars in “Why Elephants Can’t Dance but Hamsters Can Skydive.” He talks about an issue discussed in Homework Problem 28 in Chapter 2 of the Fourth edition of Intermediate Physics for Medicine and Biology. Russ Hobbie and I ask the reader to analyze how fast animals of different sizes fall. In “Why Elephants Can’t Dance,” Miodownik performs a brilliant demonstration using two spherical animals—one about the size of a hamster, and the other about the size of a dog—made of some sort of jello-like gel. Suffice to say, the hamster-sized blob of gel does just fine when it hits the ground after a fall, but the dog-sized blob has some problems. The audience for the lecture is mostly children, but as Dickens wrote “it is good to be children sometimes, and never better than at Christmas.” The entire lecture is about why size matters in the animal kingdom.

Miodownik then talks about another topic in animal scaling that Russ and I don’t mention in our book, although I often bring it up when I teach Biological Physics at Oakland University. In two animals with the same shape but different sizes, their weight increases as the cube of their linear size, but the cross-sectional area of their legs increases as the square of the size. Therefore, a large animal has a harder time supporting its weight than a small animal does. Miodownik demonstrates this with two rubber pig-like spheres with rubber legs attached. The small sphere easily stands on its legs, while the large sphere just collapses. As the video says, size really does matter. Of course, elephants solve this problem by making their legs thick, which is why they can’t dance.

I recommend watching “Why Elephants Can’t Dance” while reading Chapter 2 of Intermediate Physics for Medicine and Biology. It will help you understand animal scaling.

Enjoy the Royal Institution Christmas Lectures, and have a Merry Christmas.