Friday, August 28, 2015

Art Winfree and the Bidomain Model of Cardiac Tissue

Art Winfree was a pioneer in applying physics and mathematics to cardiac electrophysiology. Russ Hobbie and I cite him often in the 5th edition of Intermediate Physics for Medicine and Biology. After his untimely death in 2002, I was asked to write an article for a special issue of the Journal of Theoretical Biology published in his honor. My paper, “Art Winfree and the Bidomain Model of Cardiac Tissue,” appeared in 2004.

My original submission for the special issue was a personal tribute to Art. It began
“Spiral waves have become so popular in Tucson they are even sold in hair styling salons (Figure 1)”
A photograph in a preprint from Art Winfree, with the caption "Spiral waves have become so popular in Tucson they are even sold in hair styling salons (Figure 1)"
Figure 1.
I had to laugh as I read the above quote in a preprint Art Winfree sent me. It was to be the opening sentence of a chapter appearing in a prestigious textbook on cardiac electrophysiology. Unfortunately, the sentence and the picture were deleted before the book's publication, although the picture (Fig. 1) did appear eventually in the second edition of Art’s The Geometry of Biological Time. For me, the quote captures the essence of Art: his humor, his irreverence, and his uncanny ability to find science in the world around him. I only met Art in person once, but we corresponded often by email, exchanging ideas, frustrations, and gossip. Of all the scientists who have influenced my research career, only my PhD advisor John Wikswo had a greater impact than Art Winfree did. In this paper, I describe several instances where my path crossed Art’s as we each attacked related problems in cardiac electrophysiology. In addition, I hope to show that Art made important contributions to what is known as the “bidomain model” of cardiac tissue.
Later in the article is one of my favorite passages.
I recall vividly a sunny day in April, soon after my second daughter Katherine was born. I was sitting on a rocking chair in the living room of our house in Kensington, Maryland, holding the sleeping infant in one arm and Art’s book When Time Breaks Down in the other. Outside I could see our dogwood tree in full blossom. As I read page after page, I remember thinking “life doesn’t get any better than this.” The book (and the daughter) changed my life.
Unfortunately, the editors of the special issue didn’t like my paper, saying they wanted a more traditional review article. In particular, they objected to my quoting Art’s emails he had sent me. So, I gave the paper a lobotomy and published a harmless but lifeless review. When the issue came out, I found a wonderful article by George Oster about Winfree, full of personal insights and even the text of one of Art’s emails. I wish now I had pushed harder to get my article published in its original form. The best article in the special issue was “Art Winfree, Artist of Science” by his daughter Rachael Winfree.

In the acknowledgments of my paper is the line “I would like to thank Jesse Malouf for his help editing this paper.” Jesse was a student in my honors college course about Pacemakers and Defibrillators. At Oakland University, Honors College has many of the best students in the university, but they are from all backgrounds and often have weak math skills. Jesse was a mathaphobe, but a wonderful writer. On one of my exams I had a mixture of questions, some requiring mathematical analysis and others needing an essay. Jesse skipped the math questions, but to make up for it he not only answered all the essay questions elegantly but also wrote a “bonus essay”. I never had a student hand in a bonus essay before! The next semester, I hired him to help me write the Winfree article. I fear many of his contributions to the original version were not included in the published one.

In the “olden days” the original draft of my Winfree article would be lost forever, or maybe would sit in some file cabinet unseen for decades. But nowadays, you can find anything on the internet (how did we live without it?). I have posted the original submission on my ResearchGate page. You can find it here.

Friday, August 21, 2015

The Coulter Counter

In Intermediate Physics for Medicine and Biology, Russ Hobbie and I often include applications of important topics in the homework problems. One such problem, new in Chapter 6 of the 5th edition, is an analysis of a Coulter counter.
Problem 23. The Coulter counter or resistive pulse technique is used to count and size particles in a wide variety of applications (Kubitschek 1969; DeBlois and Bean 1970), including the automated counting of blood cells. The cells being counted are assumed to be nonconducting and immersed in a conducting fluid. The fluid is made to flow through a narrow channel. When a suspended particle enters the channel there is a change in resistance. Assume a long channel of radius b with no end effects.
(a) What is the resistance of pure fluid of resistivity ρ = 1/σ in a segment of channel of length 2a?
(b) A cylindrical non-conducting cell of radius a and length 2a is in the channel. Its axis and the axis of the channel coincide. What is the resistance of a segment of channel of length 2a? Ignore end effects.
(c) Show that the resistance change (the difference between these two results) is proportional to the volume of the cell, V=2πa3, and inversely proportional to b4.
In the August issue of Physics Today is an article about extending the Coulter counter to sequencing DNA. Murugappan Muthukumar, Calin Plesa, and Cees Dekker write
In the 1940s Wallace Coulter set about finding a way to quickly count blood cells, which at the time was a slow and inefficient process. His approach was to pass cells, one by one, through a small hole connecting to compartments filled with electrolyte solution. Simultaneously, he applied a voltage across the compartment and measured the ionic current passed through the hole. As a cell passed through the hole, it would partially block the flow of electric charges, and the current would drop by an amount proportional to the volume of the cell….Coulter’s technique worked out wonderfully and revolutionized cell counting.
Then, the authors describe how this method can be used to sequence DNA.
The last two decades have seen a renaissance of the Coulter counter concept. The principle remains essentially the same, but nanopores—holes with a diameter of merely a few nanometers—have shrunk the length scale from that of single cells to that of single molecules. When DNA molecules are added to one side of the pore and an electric field is applied, the resulting electrophoretic force on the negatively charged DNA can pull the molecule through the pore in a head-to-tail fashion, leading to an observable blockade in the ionic current…

In the 1990s several research groups … began probing whether the different bases on a DNA strand might block measurably different amounts of ionic current as they pass through a nanopore. If so, the pattern of current generated by a DNA strand threaded through a nanopore might provide a linear readout of the strand’s base sequence… Although significant challenges remain to turn that vision into a practical reality, the goal appears to be within reach.
The authors then describe more details about the technique. Some use transmembrane proteins like the membrane channels described in Chapter 9 of IPMB. Others use tiny holes drilled into sheets of silicon nitride. Still others use a hybrid of these two.

Clearly the method will not work unless the DNA is a single strand. Wanunu (2012) discusses the molecular dynamics involved in unzipping a double strand to obtain two single strands, one of which can then be threaded through the pore to do the sequencing. The nanopores must be very narrow if you are to have any chance of distinguishing different bases attached to the DNA backbone.

Russ and I had no idea about these modern uses of the Coulter counter when we added the homework problem. This new application of the Coulter idea shows how a strong understanding of the fundamentals of physics as applied to medicine and biology can allow one to quickly move to the forefront of cutting-edge new technologies.

Friday, August 14, 2015

The Psychic Probe

Foundation,  by Isaac Asimov, superimposed on Intermediate Physics for Medicine and Biology.
Foundation,
by Isaac Asimov.
In the summer my wife and I sometimes take long car trips, and I often listen to audiobooks while I drive to keep me awake and alert. During a recent trip I listened to Isaac Asimov’s Foundation trilogy. Regular readers of this blog know that I’m a huge Asimov fan. I first read the Foundation series about forty years ago, and this was my third or fourth time through these delightful books.

In brief, the Foundation series tells the history of the decaying galactic empire, and describes the work of the psychohistorian Hari Seldon who has calculated mathematically how to reduce the duration of the dark ages following the empire’s fall from 30,000 years to merely 1000. All goes according to plan until the Mule, a mutant who can control other people’s emotions, causes all to go awry.

Foundation and Empire,  by Isaac Asimov, superimposed on Intermeidate Physics for Medicine and Biology.
Foundation and Empire,
by Isaac Asimov.
One of Asimov’s inventions in this future history is a device that can read minds, called the Psychic Probe. He writes in Foundation and Empire,
The general threw away his shredded, never-lit cigarette, lit another, and shrugged. “Well, it is beside the immediate point, this lack of first-class tech-men. Except that I might have made more progress with my prisoner were my Psychic Probe in proper order.”

The secretary’s eyebrows lifted. “You have a Probe?”

“An old one. A superannuated one which fails me the one time I needed it. I set it up during the prisoner’s sleep, and received nothing. So much for the Probe. I have tried it on my own men and the reaction is quite proper, but again there is not one among my staff of tech-men who can tell me why it fails upon the prisoner. Ducem Barr, who is a theoretician of parts, though no mechanic, says the psychic structure of the prisoner may be unaffected by the Probe since from childhood he has been subjected to alien environments and neural stimuli. I don’t know. But he may yet be useful. I save him in that hope.” 
Second Foundation,  by Isaac Asimov, superimposed on Intermediate Physics for Medicine and BIology.
Second Foundation,
by Isaac Asimov.
Russ Hobbie and I don’t mention the Psychic Probe in the 5th edition of Intermediate Physics for Medicine and Biology … or do we? Asimov didn’t explain the physical mechanism behind the Probe, but I can speculate. Four candidates are:
Asimov's Foundation Trilogy, superimposed on Intermediate Physcs for Medicine and Biology.
Asimov's Foundation Trilogy.
PET and fMRI are too slow to accurately follow rapid brain activity. PET detects brain metabolism and fMRI detects blood flow, both of which are only indirectly related to neuron firing. My best guess for the Psychic Probe is some combination of MEG and TMS. Apparently the probe can damage the brain when used aggressively, which suggests TMS. But it can also read minds when used more gently, which points toward MEG. A combo TMS/MEG unit could therefore both detect and alter brain function.

While working at NIH in the 1990s, I studied both magnetoencephalography and transcranial magnetic stimulation. Yikes! I may be partially responsible for the invention of the Psychic Probe!

Friday, August 7, 2015

Kramers’ Law

When preparing the 5th edition of Intermediate Physics for Medicine and Biology, Russ Hobbie and I added a homework problem about Kramers’ law. (We spelled it Kramer’s, but his name is Kramers with an s, so we should have written Kramers’.) Kramers’ law is Eq. 16.3a, the photon energy fluence dΨ/d() as a function of frequency ν for bremsstrahlung radiation
Kramers' law.
where Z is the atomic number, h is Planck’s constant, νo is the frequency of a photon having the same energy as the incident electrons, and C is a constant. In his paper “On the Theory of X-ray Absorption and of the Continuous X-ray Spectrum” (Philosophical Magazine, Volume 46, Pages 836–871, 1923), Kramers writes
The continuous x-ray spectrum has in the course of the last years been investigated by a number of physicists. The problem is here to determine how, for a given tension [voltage] on the tube and a given anticathode material [typically tungsten], the energy in the continuous spectrum is distributed among different frequencies…

The object of the present paper is to show how it is possible to account theoretically for the main features of the phenomena of x-ray absorption and continuous x-ray emission discussed above. The explanation of these phenomena may be traced back to the determination of the radiation processes which may occur when a free electron of given velocity approaches a positive nucleus with given charge.
Who was Kramers? According to the Dictionary of Scientific Biography, Hendrik Anthony (Hans) Kramers was born in Rotterdam, the Netherlands in 1894. He joined Niels Bohr’s Institute of Theoretical Physics, and in 1934 he moved to Leiden University, where he remained until his death in 1952. He’s known for many contributions to physics, including the Kramers-Kronig relations. The Dictionary of Scientific Biography article concludes
Kramers’ work, which covers almost the entire field of theoretical physics, is characterized both by outstanding mathematical skill and by careful analysis of physical principles. It also leaves us with the impression that he tackled problems because he found them challenging, not primarily because they afforded chances of easy success. As a consequence his work is somewhat lacking in spectacular results that can be easily explained to a layman; but among fellow theoreticians he was universally recognized as one of the great masters.
A Tale of Two Continents, by Abraham Pais, superimposed on Intermediate Physics for Medicine and Biology.
A Tale of Two Continents,
by Abraham Pais.
Here is my favorite Kramers story. Jewish physicist Abraham Pais described in his autobiography A Tale of Two Continents how he spent much of World War II in Holland hiding from the Gestapo. Kramers was one of the few people who knew of his hiding place, and would visit him weekly to talk physics. One day when Kramers was there, Gestapo agents knocked at the door and Pais had to hide in a small enclosure behind a panel in the wall. Pais writes
I kept sitting in the tiny space, practically bent over double, holding onto the panel, when I heard the door to my room, which lay at the other side of my hiding spot, open softly. Someone entered, I did not at first know who. Then that person sat down on a small bench that stood right at the wall behind which I was folded up. The person began to read, not loud but quite softly. It was Kramers. Earlier he had lent me a volume of Bradley’s Lectures on Shakespeare. What this good man was doing now was reading to me from that book, in order to calm my nerves.

Friday, July 31, 2015

The Divergence Theorem and Stokes’ Theorem

When preparing the 5th edition of Intermediate Physics for Medicine and Biology, Russ Hobbie and I added a new homework problem to Chapter 4.
Problem 4. Integrate Eq. 4.8 over a volume and subtract the result from Eq. 4.4. The resulting relationship is called the divergence theorem.
For those of you who don’t keep a copy of IPMB always close at hand (what’s the matter with you?), Eq. 4.8 is
The differential form of the continuity equation.
where “div j” is the divergence of the vector j, and Eq. 4.4 is
The integral form of the continuity equation.
When you put the two equations together, you get (spoiler alert) ,
The divergence theorem.
also known as the divergence theorem. It is one of the fundamental results of vector calculus.

Div, Grad, Curl, and All That, by H. M. Schey, superimposed on Intermediate Physics for Medicine and Biology.
Div, Grad, Curl, and All That,
by H. M. Schey.
If you want to learn more about the divergence theorem, I recommend H. M. Schey’s book Div, Grad, Curl and All That. He writes
For the remainder of this chapter we digress from the mainstream of our narrative to discuss a famous theorem that asserts a remarkable connection between surface integrals and volume integrals. Although this relation may be suggested by the work we have done in electrostatics, the theorem is a mathematical statement holding under quite general circumstances. It is independent of any physics and is applicable in many different places. It is called the divergence theorem, and sometimes Gauss’ theorem… It says that the flux of a vector function through some closed surface equals the triple integral of the divergence of that function over the volume enclosed by the surface.
Besides the divergence theorem, another basic tenet of vector calculus is Stokes’ theorem. Can we make a similar homework problem demonstrating that? Yes! Here is a new problem for Chapter 8.
Problem 23 ½. Integrate Eq. 8.22 over a surface and subtract the result from Eq. 8.21. The resulting relationship is called Stokes’ theorem.
If you don’t have IPMB handy, Eq. 8.21 is
The integral form of Faraday's law.
and Eq. 8.22 is
The differential form of Faraday's law.
where “curl E” is the curl of the vector E. When you put the two equations together, you get
Stokes theorem.
 Schey writes
We [now] discuss another famous theorem, one strongly reminiscent of the divergence theorem and yet, as we’ll see, quite different from it. This theorem, named for the mathematician Stokes, relates a line integral around a closed path to a surface integral over what is called a capping surface of the path…In words, Stokes’ theorem says that the line integral of the tangential component of a vector function over some closed path equals the surface integral of the normal component of the curl of that function integrated over any capping surface of the path.
The divergence theorem and Stokes’ theorem are a bit too mathematical to develop in IPMB, with its emphasis on biological and medical applications. Yet there they are, implicit in our discussions of diffusion and of transcranial magnetic stimulation. If you want to learn more, start with Schey’s wonderful (and relatively inexpensive) book Div, Grad, Curl.

Friday, July 24, 2015

So You Don’t Like Error Functions?

In Chapter 4 of the 5th edition of Intermediate Physics for Medicine and Biology, Russ Hobbie and I introduce the diffusion equation (Equation 4.26)
The one dimensional diffusion equation.
where C is the concentration and D is the diffusion constant. We then study one-dimensional diffusion where initially (t = 0) the region to the left of x = 0 has a concentration Co, and the region to the right has a concentration of zero (Section 4.13). We show that the solution to the diffusion equation is (Equation 4.75)
A solution of the diffusion equation involving an error function.
where erf is the error function.

Some students don’t like error functions (really?). Moreover, often we can gain insight by solving a problem in several ways, obtaining solutions that are seemingly different yet equivalent. Let’s see if we can solve this problem in another way and avoid those pesky error functions. We will use a standard method for solving partial differential equations: separation of variables. Before I start, let’s agree to solve a slightly different problem: initially C(x,0) is Co/2 for x less than zero, and −Co/2 for x greater than zero. I do this so the solution will be odd in x: C(x,t) = −C(−x,t). At the end we can add the constant Co/2 and get back to our original problem. Now, let’s begin.

Assume the solution can be written as the product of a function of only t and a function of only x: C(x,t) = T(t) X(x). Plug this into the diffusion equation, simplify, and divide both sides by TX 
The method of separation of variables applied to the diffusion equation.
The only way that the left and right hand sides can be equal at all values of x and t is if both are equal to a constant, which I will call −k2. This gives two ordinary differential equations
The equation for the temporal part of the diffusion equation after separation of variables.
 and
The solution to the first equation is an exponential
The solution to the equation for the temporal part of the diffusion equation after separation of variables.
and the solution to the second is a sine
The solution to the equation for the spatial part of the diffusion equation after separation of variables.
There is no cosine term because of the odd symmetry. Unfortunately, we don’t know the value of k. In fact, our solution can be a superposition of infinitely many values of k
A solution to the diffusion equation in integral form.
where A(k) specifies the weighting.

To determine A(k), use the Fourier techniques developed in Chapter 11. The result is
The Fourier transform of the solution to the diffusion equation in integral form.
How did I get that? Let me outline the process, leaving you to fill in the missing steps. I should warn you that a mathematician would worry about the convergence of the integrals we evaluate, but you and I’ll brush those concerns under the rug.

At t = 0, our solution becomes
Except for a missing factor of 2π, this looks just like the Fourier transform from Section 11.9 of IPMB. Next, multiply each side of the equation by sin(ωx), and integrate over all x. Then, use Equation 11.66b to express the integral of the product of sines as a delta function. You get
Both C(x,0) and sin(kx) are odd, so their product is even, and for x greater than zero C(x,0) is −Co/2. Therefore,
You know how to integrate sine (I hope you do!), so
Here is where things get dicey. We don’t know what cosine equals at infinity, but if we say it averages to zero the first term goes away and we get our result
Plugging in this expression for A(k) gives our solution for C(x,t). If we want to go back to our original problem with an initial condition of Co on the left and zero on the right, we must add Co/2. Thus
A solution of the one dimensional diffusion equation without error functions.
Let’s compare this solution with the one in Equation 4.75 (given above). Our new solution does not contain the error function! Those of you who dislike that function can celebrate. Unfortunately, we traded the error function for an integral that we can’t evaluate in closed form. So, you can have a function that you may be unfamiliar with and that has a funny name, or you can have an expression with common functions like the sine and the exponential inside an integral. Pick your poison.

Friday, July 17, 2015

Boston

Boston, the debut album by the rock band Boston.
Boston, by Boston.
I did something unusual for me last evening: I went to a rock concert. As a birthday present, my daughter Stephanie took me to Freedom Hill Amphitheater to listen to the band Boston.

I was in high school in 1976 when I bought Boston’s famous debut album. That was a big year: it was the bicentennial of the United States, Jimmy Carter was elected president, Nadia Comaneci was earning 10s and then-Bruce Jenner won the decathlon in the Olympic games, much to my chagrin the Cincinnati Reds won the world series (but it wasn’t quite as exciting a series as the year before, which was the best world series ever), and the Apple Computer Company was formed by Steve Jobs and Steve Wozniak. My family moved from Fort Wayne, Indiana to Ashland, Ohio, and I spent the year playing tuba in the high school band, managing the high school baseball team, reading my first Isaac Asimov book, wondering if I should study physics in college, and listening to Chicago, the Eagles, Peter Frampton, the Wings, and Boston. The severe winter of 1976-1977 in northern Ohio and the simultaneous energy crisis resulted in my high school missing several weeks of classes, so some of my friends and I had time to establish our own garage band: “Hades.” We didn’t have a singer, and my role was to pick out the melody on an electric keyboard while the guitars and drums banged out behind me. Only a few years later disco music and the Bee Gees drove me from rock to country music, which I have listened to ever since.

My ears are still tingling a bit from the concert. How loud was it? Chapter 13 in the 5th edition of Intermediate Physics for Medicine and Biology discusses the decibel scale for measuring sound intensity, a logarithmic scale defined as log10(I/Io), where the sound intensity is
-->I and the reference Io is the minimum perceptible sound (10−12 W m−2). Table 13.1 in IPMB says 120 dB is the threshold for pain, and 130 dB is typical for the peak sound at a rock concert. Stephanie and I were sitting in the back of the amphitheater, so I doubt we ever experienced 130 dB, but we were up there pretty high on the decibel scale. I probably didn’t lose any hair cells in my cochlea (see Section 13.5), but I wonder how the band plays concerts night after night without suffering hearing loss. As people age, they tend to lose the ability to hear high tones: presbycusis. I may not have heard the music last night in the same way I heard it in 1976; some of those frequencies may be lost to me forever.
The leader of Boston is Tom Scholz, their 68-year-old guitar player and keyboardist. Scholz was educated as a mechanical engineer, and is one of the few rock musicians who might enjoy reading IPMB. I found that I could identify with Scholz in some respects: he is past his prime, no longer topping the charts or breaking new ground in rock and roll. But, after decades in the business, he’s still out there performing, playing his music, and even sometimes writing new songs. It makes me want to go write another paper!

My daughter Stephanie and I at a Boston concert at Freedom Hill Amphitheater in Sterling Heights, Michigan.
My daughter Stephanie and I at a Boston concert
at Freedom Hill Amphitheater in Sterling Heights, Michigan.

Friday, July 10, 2015

The Machinery of Life

The Machinery of Life,  by David Goodsell, superimposed on Intermediate Physics for Medicine and BIology.
The Machinery of Life,
by David Goodsell.
In the very first section of the 5th edition of Intermediate Physics for Medicine and Biology (Sec. 1.1), Russ Hobbie and I discuss “Distances and Sizes.”
In biology and medicine, we study objects that span a wide range of sizes: from giant redwood trees to individual molecules. Therefore, we begin with a brief discussion of length scales.
At the end of this section, we conclude
To examine the relative sizes of objects in more detail, see Morrison et al. (1994) or Goodsell (2009).
I have talked about the book Powers of Ten by Morrison et al. previously in this blog. I have also mentioned David Goodsell’s book The Machinery of Life several times, but until today I have never devoted an entire blog entry to it.

In the 4th edition of IPMB, Russ and I cited the first edition of The Machinery of Life (1998), and that is the edition that sits on my bookshelf. When preparing the 5th edition, we updated our references, so we now cite the second edition of Goodsell's book (2009). Is there much difference between the two? Yes! Like when Dorothy left Kansas to enter Oz, the first edition is all black and white but the second edition is in glorious color. And what a difference color makes in a book that is first and foremost visual. The second edition of The Machinery of Life is stunningly beautiful. It is not just a colorized version of the first edition; it is a whole new book. Goodsell writes in the preface
I created the illustrations in this book to help bridge this gulf and allow us to see the molecular structure of cells, if not directly, then in an artistic rendition. I have included two types of illustrations with this goal in mind: watercolor paintings which magnify a small portion of a living cell by one million times, showing the arrangement of molecules inside, and computer-generated pictures, which show the atomic details of individual molecules. In this second edition of The Machinery of Life, these illustrations are presented in full color, and they incorporate many of the exciting scientific advances of the 15 years since the first edition.

As with the first edition, I have used several themes to tie the pictures together. One is that of scale. Most of us do not have a good concept of the relative sizes of water molecules, proteins, ribosomes, bacteria, and people. To assist with this understanding, I have drawn the illustrations at a few consistent magnifications. The views showing the interiors of living cells, as in the Frontispiece and scattered through the last half of the book, are all drawn at one million times magnification. Because of this consistent scale, you can flip between pages in these chapters and compare the sizes of DNA, lipid membranes, nuclear pores, and all of the other molecular machinery of living cells. The computer-generated figures of individual molecules are also drawn at a few consistent scales to allow easy comparison.

I have also drawn the illustrations using a consistent style, again to allow easy comparison. A space-filling representation that shows each atom as a sphere is used for all the illustrations of molecules. The shapes of the molecules in the cellular pictures are simplified versions of these space-filling pictures, capturing the overall form of the molecule without showing the location of every atom. The colors, of course, are completely arbitrary since most of these molecules are colorless. I have chosen them to highlight the functional features of the molecules and cellular environments.
I have often wondered how much molecular biology a biological or medical physicist needs to know. I suppose it depends on their research specialty, but in general I believe a physicist who has read and understood The Machinery of Life has most of what you need to begin working at the interface of physics and biology: An understanding of the relative scale of biological objects, an overview of the different types of biological molecules and their structures and functions, and a visual sense of how these molecules fit together to form a cell. To the physicist wanting an introduction to biology on the molecular scale, I recommend starting with The Machinery of Life. That’s why it was included in my ideal bookshelf.

Goodsell fans might enjoy visiting his website: http://mgl.scripps.edu/people/goodsell. There you can download a beautiful poster of different proteins, all drawn to scale. There are many other illustrations and publications. Enjoy!

Friday, July 3, 2015

Fermi Problems and the Annual Background Radiation Dose from Potassium-40

In the very first section of the 5th edition of Intermediate Physics for Medicine and Biology, Russ Hobbie and I emphasize the importance of being able to estimate.
One valuable skill in physics is the ability to make order of magnitude estimates, meaning to calculate something approximately right.
Our first four homework problems at the end of Chapter 1 ask the student to estimate. These exercises are examples of Fermi Problems, after physicist Enrico Fermi who was a master of the skill.

Let us try another of these problems, based on some of the concepts from Chapter 15 about the natural background radiation dose. To be specific, let’s estimate the annual background dose from the radioactive isotope potassium-40 inside our bodies.
Problem 4 ½ Estimate the annual background dose (in mSv/year) from potassium-40 in our bodies. Look up or guess any information you need for this calculation, and clearly explain any assumptions you make.
Begin by considering a single cell. Cells are about 10 microns in size, so their volume is about (10−5 m)3 = 10−15 m3 (the volume we use is not important; it will cancel out in the end). In nerves, the intracellular potassium ion concentration is about 100 mM, but this may overestimate the amount of potassium inside all types of cells. Moreover, the concentration of potassium ions in the extracellular space is much less that in the intracellular space. Let’s guess 50 mM for the average concentration of potassium in our body, meaning 50 millimoles/liter, or 50 moles/m3. If we multiply by Avogadro’s number (6 x 1023 molecules/mole), we get about 3 x 1025 molecules/m3. So, the number of potassium atoms per cell is 3 x 1010. This intermediate result is already interesting; there are over ten billion potassium ions in just one cell.

The 40K isotope is radioactive. Its abundance is about 0.01% (abundance data can be found in any table of the isotopes, or even by looking at wikipedia; I don’t know how you could guess that value from first principles). So, this means there are 3 x 106 40K atoms/cell (a little over a million). How rapidly do these decay? 40K has a half life of 1.25 x 109 years (again, see wikipedia), implying a decay rate of 0.693/1.25 x 109 = 5.5 x 10−10 decays per atom per year. Multiplying by the number of atoms/cell, we get 0.0017 decays per cell per year. This is another interesting result: an average cell has less than a one-percent chance of experiencing a 40K decay in a year. But we have a lot of cells (Russ and I estimate 2 x 1014 in IPMB), so your body suffers from about 3.4 x 1011 decays per year, or about ten thousand per second. (According to Wikipedia, this estimate is a factor of two too high, but we are not much worried about factors of two in such order-of-magnitude Fermi problems.)

How much energy does each decay deposit in our tissue? Beta decay accounts for 90% of all decays of 40K, and each decay releases an energy of about 1.3 MeV (decay energy data is a little harder to find, but appears in any good table of the isotopes; if you had merely guessed that 1 MeV is the order of magnitude of nuclear decay energies, you would not be too far off). Some of that energy goes to a neutrino, which leaves the body. Let’s assume that on average about 30% of the beta decay energy goes to the electron and that no electrons escape the body, so each decay deposits 0.4 MeV into the cell, or 6 x 10−14 joules. A gray (the unit of dose) is a joule per kilogram, so to calculate dose we need the mass of a cell, which to a first approximation is the product of the density of water (103 kg/m3) times the volume of the cell, or 10−12 kg (I told you the volume of the cell would cancel out). So, the dose is the number of decays per year (0.0017) times the energy per decay (6 x 10−14 J) divided by the mass (10−12 kg), or 0.0001 gray/year. Another unit of dose is the sievert, which accounts for biological damage in addition to energy deposition. A gray and a sievert are the same for electrons, so the annual background dose is 0.0001 Sv, or 0.1 mSv.

In Table 16.6 of IPMB, Russ and I estimate the annual background dose from all internal sources is about 0.3 mSv. Because 40K is not the only isotope in our body that is decaying (for example, carbon-14 is another), we seem to have gotten our order-of-magnitude estimate pretty close. One goal of Chapters 16 and 17 in IPMB is to refine such calculations. For medical purposes you need more accuracy; for a Fermi problem we did okay.

Friday, June 26, 2015

The Electric Potential of a Rectangular Sheet of Charge - Resolved

In the January 9 post of this blog, I challenged readers to find the electrical potential V(z) that will give you the electric field E(z) of Eq. 6.10 in the 5th edition of Intermediate Physics for Medicine and Biology
An expression for the electrical field produced by a rectangular sheet of charge.
In other words, the goal is to find V(z) such that E = − dV/dz produces Eq. 6.10. In the comments at the bottom of the post, a genius named Adam Taylor made a suggestion for V(z) (I love it when people leave comments in this blog). When I tried his expression for the potential, it almost worked, but not quite (of course, there is always a chance I have made a mistake, so check it yourself). But I was able to fix it up with a slight modification. I now present to you, dear reader, the potential:

A complicated expression for the electrical potential of a rectangular sheet of charge.

How do you interpret this ugly beast? The key is the last term, z times the inverse tangent. When you take the z derivative of V(z), you must use the product rule on this term. One derivative in the product rule eliminates the leading z and gives you exactly the inverse tangent you need in the expression for the electric field. The other gives z times a derivative of the inverse tangent, which is complicated. The two terms containing the logarithms are needed to cancel the mess that arises from differentiating tan−1.

I don’t know what there is to gain from having this expression for the potential, but somehow it comforts me to know that if there is an analytic equation for E there is also an analytic equation for V.