Friday, May 31, 2013

Rounding Off the Cow

In the October 2012 issue of the American Journal of Physics, Dawn Meredith and Jessica Bolker published an article about “Rounding Off the Cow: Challenges and Successes in an Interdisciplinary Physics Course for Life Science Students” (Volume 80, Pages 913–922). The article is interesting, and much of the motivation for their work is nearly identical to that of Russ Hobbie and I in writing the 4th edition of Intermediate Physics for Medicine and Biology. They focus on an introductory physics class, whereas Russ and I wrote an intermediate level textbook. Nevertheless, many of the ideas and challenges are the same. Here, I want to focus on their Table 1, in which they list topics that are emphasized and deemphasized compared to standard introductory classes.

Table I. Changes in topic emphasis compared to standard course
Semester 1Semester 2
Included/stressedKinematicsHeat transfer
DynamicsKinetic theory of gases
Static torqueEntropy
EnergyDiffusion, convection, conduction
Stress/strain and fractureSimple harmonic motion
Fluids (far more)Waves (sound, optics)
Omitted/de-emphasizedProjectile motionHeat engines
Relative motionMagnetism (less)
Rotational motionInduction (qualitatively)
StaticsAtomic physics (instrumentation)
CollisionsRelativity
Newton’s law of gravitation
Kepler’s laws

How does this list compare with the content of IPMB? We don’t stress kinematics and dynamics much. In fact, most of our mechanics discussion centers on static equilibrium. Interestingly, Meredith and Bolker emphasize static torque, which is absolutely central to our analysis of biomechanics in Chapter 1. Rotational equilibrium and torque is what explains why bones, muscles and tendons often experience forces far larger than the weight of the body. It also underlies our rather extensive discussion of the role of a cane. We discuss mechanical energy in Chapter 1, but energy doesn’t become an essential topic until our Chapter 3 on thermodynamics. We agree completely with Meredith and Bolker’s listing of “stress/strain and fracture” and “Fluids (far more)”, and I second the “far more”. Our Chapter 1 contains a lot of fluid dynamics, including the biologically-important concept of buoyancy, the idea of high and low Reynolds number, and applications of fluid dynamics to the circulatory system.

The time allotted to an introductory physics class is limited, so something must get deemphasized to free up time for topics like fluids. Meredith and Bolker mention projectile motion (we agree, it is nowhere in IPMB), relative motion (not crucial if not covering relativity), and rotational motion (we don’t emphasize this either, except when analyzing the centrifuge). I don’t really understand the omission of statics, because as I said earlier static mechanical equilibrium is crucial for biomechanics. They deemphasize collisions, and so do we, although we do discuss the collision of an electron with a photon when analyzing Compton Scattering in Chapter 15. Newton’s law of gravity and Kepler’s laws of planetary motion are absent from both our book and their class.

In the second semester, Meredith and Bolker stress heat transfer (convection and conduction), the kinetic theory of gases, and entropy. Russ and I discuss all these topics in our Chapter 3. Diffusion is a topic they emphasize, and rightly so. It is a topic that is typically absent from an introductory physics class, but is crucial for biology. We discuss it in detail in Chapter 4 of IPMB. Meredith and Bolker list simple harmonic motion among the topics they stress. We talk about harmonic motion in Chapter 10, but mainly as a springboard for the study of nonlinear dynamics. Much of the analysis of linear harmonic motion is found in IPMB in an appendix. Finally, they stress waves (sound and optics). We do too, mainly in our Chapter 13 about sound and ultrasound; a new chapter in the 4th edition.

Topics they omit or deemphasize in the second semester include heat engines. We barely mention heat engines at the end of Chapter 3, and the well-known Carnot heat engine is never analyzed in our book. Meredith and Bolker deemphasize magnetism and magnetic induction. As a researcher in biomagnetism, I would hate to see these topics go. Russ and I analyze biomagnetism in Chapter 8. However, I could see how one might be tempted to deemphasize these topics; biomagnetic fields are very weak and do not play a large role in either biology or medicine. I personally would keep them in, and they remain an important part of IPMB. They do not stress “Atomic Physics (Instrumentation)” and I am not sure exactly what they mean, especially with their parenthetical comment about instruments. We talk a lot about atomic physics in Chapter 14 on Atoms and Light. Finally, Meredith and Bolker omit relativity, and so do Russ and I, except as needed to understand photons. We never discuss the more traditional phenomena of relativity, such as the Lorentz contraction, time dilation, or simultaneity.

Some topics should get about the same amount of attention as in a traditional class, but with slight changes in emphasis. For instance, I would cover geometrical optics, including lenses (when discussing the eye and eyeglasses) but I would skip mirrors. I would cover nuclear physics, but I would skip fission and fusion, and focus on radioactive decay, including positron decay (positron emission tomography).

I think that Meredith and Bolker provide some useful guidance on how to construct an introductory physics class for students interested in the life sciences. Russ and I aim at an intermediate class for students who have taken a traditional introductory class and want to explore applications to biology and medicine in more detail. Our book is clearly at a higher mathematical level: we use calculus while most introductory physics classes for life science majors are algebra based. But for the most part, we agree with Meredith and Bolker about what physics topics are central for biology majors and pre-med students.

Friday, May 24, 2013

Eleanor Adair (1926-2013)

Eleanor Adair, who studied the health risks of microwave radiation, died on April 20 in Hamden, Connecticut at the age of 86. A 2001 interview with Adair, published in the New York Times, began
Eleanor R. Adair wants to tell the world what she sees as the truth about microwave radiation.

New widely reported studies have failed to find that cellular phones, which use microwaves to transmit signals, cause cancer. And most academic scientists say the microwave radiation that people are exposed to with devices like cell phones is harmless. But still, Dr. Adair knows that many people deeply fear these invisible rays.

She knows that many people hear the word “radiation” and assume that all radiation is dangerous, equating microwaves to the very different X-rays.

Microwaves, she points out, are at the other end of the electromagnetic spectrum from high energy radiation like X-rays and gamma rays. And unlike gamma rays and X-rays, which can break chemical bonds and injure cells, even causing cancer, microwaves, she says, can only heat cells. Of course, if cells get hot enough, they can die, but the heat level has to be closer to that in an oven than the extremely low level from cell phones.
The interview ends with this exchange:
Q. If I were to say to people, “Hey there’s this really cool idea: Why heat your whole house when you could use microwaves to heat yourself?” they would say: “You’ve got to be kidding. Don’t you know that microwaves are dangerous? They can even cause cancer.” What do you say to people who respond like that?

A. I try to educate them in exactly what these fields are. That they are part of the full electromagnetic spectrum that goes all the way from the radio frequency and microwave bands, through infrared, ultraviolet, the gamma rays and all that.

And the difference between the ionizing X-ray, gamma ray region and the microwave frequency is in the quantum energy. The lower you get in frequency the lower you get in quantum energy and the less it can do to the cells in your body.

If you have a really high quantum energy such as your X-rays and ionizing-radiation region of the spectrum, this energy is high enough that it can bump electrons out of the orbit in your cells and it can create serious changes in the cells of your body such that they can turn into cancers and various other things that are not good for you.

But down where we are working, in the microwave band, you are millions of times lower in frequency and there the quantum energy is so low that they can’t do any damage to the cells whatsoever. And most people don’t realize this.

Somehow, something is missing in their basic science education, which is something I keep trying to push. Learn the spectrum. Learn that you’re in far worse shape if you lie out on the beach in the middle of summer and you soak up that ultraviolet radiation than you are if you use your cell phone.

Q. Some people say that with the ever-increasing exposure of the population to microwaves—cell phones have really taken off in the past few years—we need to redouble our research efforts to look for dangerous effects of microwaves on cells and human tissues. Do you agree?

A. No. All the emphasis that we need more research on power line fields, cell phones, police radar—this involves billions of dollars that could be much better spent on other health problems. Because there is really nothing there.
We don’t cite Adair’s research in the 4th edition of Intermediate Physics for Medicine and Biology, but we do cover the interaction of electromagnetic fields with tissue in Chapter 9. Much of our discussion is about powerline (60 Hz) fields, but many of the same considerations apply to microwaves. In our discussion, we do cite Robert Adair, Eleanor’s husband and an emeritus professor of physics at Yale, who shares his wife’s interest in the health effect of microwave radiation.

Adair won the d’Arsonval Award, presented by the Bioelectromagnetics Society, to recognize her accomplishments in the field of bioelectromagnetics. In an editorial announcing the award, Ben Greenebaum writes (Bioelectromagnetics, Volume 29, Page 585, 2008)
It gives me great pleasure to introduce Dr. Eleanor R. Adair, the recipient of the Bioelectromagnetics Society’s 2007 D’Arsonval Award, as she presents her Award Lecture (Fig. 1). Dr. Adair is being honored by the Society for her body of work investigating physiological thermoregulatory responses to radio frequency and microwave fields. Her bioelectromagnetic career began with extensive experimental studies of electromagnetic radiation-induced thermophysiological responses in monkeys and concluded with experiments that accomplished the critical extrapolation of the earlier findings to humans. I believe that this body of work constitutes a majority of the literature on the latter topic.

She spent most of her career as a research scientist at the John B. Pierce Foundation Laboratory at Yale University, but finished it as a scientist at the US Air Force’s Brooks City Base in San Antonio, Texas. As she notes in her D’Arsonval address [Adair, 2008], she took her undergraduate degree at Mount Holyoke College in 1948 and her doctorate in psychology at the University of Wisconsin-Madison in 1955. Interspersed among her academic accomplishments in Madison were others—marriage to Robert Adair and children. We should not forget that combining a research career and family at that time was much rarer and required overcoming greater difficulties than those still encountered today. Those of us who have interacted with Dr. Adair over the years know that she has determination in plenty.

Dr. Adair was a charter member of the Society and was its Secretary-Treasurer (1983–1986) during a difficult time, when the Society decided to replace its first Executive Director with Bill Wisecup. She has also been active outside the Society, both with groups concerned with research into bioelectromagnetic effects and with groups concerned with the implications of these results.

However, it is for her overall scientific contributions to bioelectromagnetics that she is being presented the D’Arsonval Award. The criteria for the Award state that “. . . the D’Arsonval Medal is to recognize outstanding achievement in research in the field of Bioelectromagnetics.” And that is the topic that she will address today in her presentation entitled, “Reminiscences of a Journeyman Scientist.”
For those who want to read Adair's own words, you can find her presentation at:
Adair. E. R. (2008) “Reminiscences of a journeyman scientist: Studies of thermoregulation in non-human primates and humans,”  Bioelectromagnetics  Volume 29, Pages 586–597.

Friday, May 17, 2013

The Lorenz equations and chaos

Fifty years ago, Edward Lorenz (1917–2008) published an analysis of Rayleigh-Benard convection that began the study of a field of mathematics called chaos theory. In the 4th edition of Intermediate Physics for Medicine and Biology, Russ Hobbie and I introduce chaos by analyzing the logistic map, which
is an example of chaotic behavior or deterministic chaos. Deterministic chaos has four important characteristics: 1. The system is deterministic, governed by a set of equations that define the evolution of the system. 2. The behavior is bounded. It does not go off to infinity. 3. The behavior of the variables is aperiodic in the chaotic regime. The values never repeat. 4. The behavior depends very sensitively on the initial conditions.
The sensitivity to initial conditions is sometimes called the “butterfly effect,” a term coined by Lorenz. His model is a simplified description of the atmosphere, and has implications for weather prediction.

The mathematical model that Lorenz analyzed consists of three first-order coupled nonlinear ordinary differential equations. Because of their historical importance, I have written a new homework problem that introduces Lorenz’s equations. These particular equations don’t have any biological applications, but the general idea of chaos and nonlinear dynamics certainly does (see, for example, Glass and Mackey’s book From Clock’s to Chaos.
Section 10.7

Problem 33 1/2. Edward Lorenz (1963) published a simple, three-variable (x, y, z) model of Rayleigh-Benard convection
dx/dt = σ (y – x)
dy/dt = x (ρ – z) – y
dz/dt = xy – β z
where σ=10, ρ=28, and β=8/3.
(a) Which terms are nonlinear?
(b) Find the three equilibrium points for this system of equations.
(c) Write a simple program to solve these equations on the computer (see Sec. 6.14 for some guidance on how to solve differential equations numerically). Calculate and plot x(t) as a function of t for different initial conditions. Consider two initial equations that are very similar, and compute how the solutions diverge as time goes by.
(d) Plot z(t) versus x(t), with t acting as a parameter of the curve.

Lorenz, E. N. (1963) “Deterministic nonperiodic flow,” Journal of the Atmospheric Sciences, Volume 20, Pages 130–141.
If you want to examine chaos in more detail, see Steven Strogatz’s excellent book Nonlinear Dynamics and Chaos. He has an entire chapter (his Chapter 9) dedicated to the Lorenz equations.

The story of how Lorenz stumbled upon the sensitivity of initial conditions is a fascinating tale. Here is one version in a National Academy of Sciences Biographical Memoir about Lorenz written by Kerry Emanuel.
At one point, in 1961, Ed had wanted to examine one of the solutions [to a preliminary version of his model that contained 12 equations] in greater detail, so he stopped the computer and typed in the 12 numbers from a row that the computer had printed earlier in the integration. He started the machine again and stepped out for a cup of coffee. When he returned about an hour later, he found that the new solution did not agree with the original one. At first he suspected trouble with the machine, a common occurrence, but on closer examination of the output, he noticed that the new solution was the same as the original for the first few time steps, but then gradually diverged until ultimately the two solutions differed by as much as any two randomly chosen states of the system. He saw that the divergence originated in the fact that he had printed the output to three decimal places, whereas the internal numbers were accurate to six decimal places. His typed-in new initial conditions were inaccurate to less than one part in a thousand.

“At this point, I became rather excited,” Ed relates. He realized at once that if the atmosphere behaved the same way, long-range weather prediction would be impossible owing to extreme sensitivity to initial conditions. During the following months, he persuaded himself that this sensitivity to initial conditions and the nonperiodic nature of the solutions were somehow related, and was eventually able to prove this under fairly general conditions. Thus was born the modern theory of chaos.
To learn more about the life of Edward Lorenz, see his obituary here and here. I have not read Chaos: Making a New Science by James Gleick, but I understand that he tells Lorenz’s story there.

Friday, May 10, 2013

Graduation

Today, my wife Shirley and I will attend the graduation of our daughter Kathy from Vanderbilt University. She is getting her undergraduate degree, with a double major in biology and history.

Kathy spent part of her time in college working with John Wikswo in the Department of Physics. As regular readers of this blog may know, Wikswo was my PhD advisor when I was a graduate student at Vanderbilt in the 1980s. Russ Hobbie and I often cite Wikswo’s work in the 4th edition of Intermediate Physics for Medicine and Biology, for his contributions to both cardiac electrophysiology and biomagnetism. You can see a picture of Kathy and John in an article about undergraduate research in Arts and Science, the magazine of Vanderbilt University’s College of Arts and Science. Interestingly, there are now publications out there with “Roth and Wikswo” among the authors that I had nothing to do with; for example, see poster A53 at the 6th q-bio conference (Santa Fe, New Mexico, 2012). You can watch and listen to Wikswo give his TEDxNashville talk here. Kathy also worked with Todd Graham of the Department of Cell and Developmental Biology at Vanderbilt. This fall, she plans to attend graduate school studying biology at Michigan State University.

Friday, May 3, 2013

Biodamage Via Shock Waves Initiated by Irradiation With Ions

Just two doors down the hall from my office in Hannah Hall of Science at Oakland University is the office of my friend Gene Surdutovich. Gene teaches physics at OU, and is an international expert on the interaction of heavy ions with biological tissue. Recently he and his collaborators have published a review describing their work, titled “Biodamage Via Shock Waves Initiated by Irradiation With Ions,” in the journal Scientific Reports (Volume 3, Article Number 1289). Their abstract states
Radiation damage following the ionising radiation of tissue has different scenarios and mechanisms depending on the projectiles or radiation modality. We investigate the radiation damage effects due to shock waves produced by ions. We analyse the strength of the shock wave capable of directly producing DNA strand breaks and, depending on the ion’s linear energy transfer, estimate the radius from the ion’s path, within which DNA damage by the shock wave mechanism is dominant. At much smaller values of linear energy transfer, the shock waves turn out to be instrumental in propagating reactive species formed close to the ion’s path to large distances, successfully competing with diffusion.
Except for the British spelling, I enjoyed reading this paper very much.

Russ Hobbie and I discuss the interaction of ions with tissue in the 4th edition of Intermediate Physics for Medicine and Biology, in the context of proton therapy.
Protons are also used to treat tumors. Their advantage is the increase of stopping power at low energies. It is possible to make them come to rest in the tissue to be destroyed, with an enhanced dose relative to intervening tissue and almost no dose distally (“downstream”) as shown by the Bragg peak in Fig. 16.51.
Surdutovich and his colleagues note that the Bragg peak occurs for heavier ions too, such as carbon. What is really fascinating about their work is one of the mechanisms they propose. They note that
while the Bragg peak location answers a question of where most of the damage occurs, we are raising a question of how the damage happens… We investigate the effects that stem from a large inhomogeneity of the dose distribution in the vicinity of the Bragg peak on biological damage.
Their hypothesis is that when a carbon ion interacts with tissue, energy is abruptly deposited
within a cylinder of about one nm radius, which is so small that the temperature within this cylinder increases by over 1000 K by 10−13 s (we will refer to it as the ‘‘hot cylinder’’). This increase of temperature brings about a rapid increase of pressure (up to 1 GPa) compared to the atmospheric pressure outside the cylinder. Such circumstances cause the onset of a cylindrical shock wave described by the strong explosion scenario. The pressure rapidly increases on the wave front and then decreases in the wake.
I find this mechanism to be fascinating. The theory is noteworthy for its multiscale approach: analyzing events spanning a wide range of time, space, and energy scales. Interestingly, their theory also spans multiple chapters in Intermediate Physics for Medicine and Biology: Chapter 1 about pressure, Chapter 3 on temperature and heat, Chapter 4 on diffusion, Chapter 13 on sound, Chapter 15 about the interaction of radiation with tissue, and Chapter 16 about the medical use of radiation. They conclude
The notion of thermomechanical effects represents a paradigm shift in our understanding of radiation damage due to ions and requires re-evaluation of relative biological effectiveness because of collective transport effects for all ions and direct covalent bond breaking by shock waves for ions heavier than argon. These effects will also have to be considered for high-density ion beams, irradiation with intensive laser fields, and other conditions prone to causing high gradients of temperature and pressure on a nanometre scale.
This paper appears in an open access online journal, so you don’t need a subscription to read it. Be sure to look at the paper’s supplementary information, especially the movie showing a molecular dynamics simulation of the shock wave distorting a DNA molecule. And if you read very closely, you will find this nugget appearing in the acknowledgments: “We are grateful to … Bradley Roth who critically read the manuscript.”

By the way, this fall Gene is scheduled to teach PHY 325, Biological Physics, using the textbook….you guessed it….Intermediate Physics for Medicine and Biology.

Friday, April 26, 2013

Molecular Structure of Nucleic Acids: A Structure for Deoxyribose Nucleic Acid

Sixty years ago yesterday the journal Nature published a letter by James Watson and Francis Crick titled “Molecular Structure of Nucleic Acids: A Structure for Deoxyribose Nucleic Acid” (Volume 171, Pages 737–738), announcing the discovery of the double helix structure of DNA.

The 4th edition of Intermediate Physics for Medicine and Biology doesn’t discuss the structure of DNA much. As Russ Hobbie and I say in the preface, “Molecular biophysics has been almost completely ignored: excellent texts already exist, and this is not our area of expertise.” Yet, we do mention DNA occasionally. In the first section of our book, about distances and sizes, we say
Genetic information is stored in long, helical strands of deoxyribonucleic acid (DNA). DNA is about 2.5 nm wide, and the helix completes a turn every 3.4 nm along its length.
Problem 2 in the first chapter analyzes the volume of DNA
Problem 2 Our genetic information or genome is stored in the parts of the DNA molecule called base pairs. Our genome contains about 3 billion (3×109) base pairs, and there are two copies in each cell. Along the DNA molecule, there is one base pair every one-third of a nanometer. How long would the DNA helix from one cell be if it were stretched out in a line? If the entire DNA molecule were wrapped up into a sphere, what would be the diameter of that sphere?
A problem in Chapter 3 considers errors in DNA in the context of the Boltzmann factor.
Problem 25 The DNA molecule consists of two intertwined linear chains. Sticking out from each monomer (link in the chain) is one of four bases: adenine (A), guanine (G), thymine (T), or cytosine (C). In the double helix, each base from one strand bonds to a base in the other strand. The correct matches, A–T and G–C, are more tightly bound than are the improper matches. The chain looks something like this, where the last bond shown is an “error.

A  T  G  C  G
T  A  C  G  A (error)

The probability of an error at 300 K is about 10−9 per base pair. Assume that this probability is determined by a Boltzmann factor e−U/kBT, where U is the additional energy required for a mismatch.
(a) Estimate this excess energy.
(b) If such mismatches are the sole cause of mutations in an organism, what would the mutation rate be if the temperature were raised 20° C?
We discuss DNA again in Chapter 16, when considering radiation damage to tissue.
Cellular DNA is organized into chromosomes. In order to understand radiation damage to DNA, we must recognize that there are four phases in the cycle of cell division

Figure 16.33 shows, at different magnifications, a strand of DNA, various intermediate structures which we will not discuss, and a chromosome as seen during the M phase of the cell cycle. The size goes from 2 nm for the DNA double helix to 1400 nm for the chromosome. In addition to cell survival curves one can directly measure chromosome damage. There is strong evidence that radiation, directly or indirectly, breaks a DNA strand. If only one strand is broken, there are efficient mechanisms that repair it over the course of a few hours using the other strand as a template. If both strands are broken, permanent damage results, and the next cell division produces an abnormal chromosome.19 Several forms of abnormal chromosomes are known, depending on where along the strand the damage occurred and how the damaged pieces connected or failed to connect to other chromosome fragments. Many of these chromosomal abnormalities are lethal: the cell either fails to complete its next mitosis, or it fails within the next few divisions. Other abnormalities allow the cell to continue to divide, but they may contribute to a multistep process that sometimes leads to cancer many cell generations later.
The Double Helix, by James Watson, superimposed on Intermeidate Physics for Medicine and Biology.
The Double Helix,
by James Watson.
The story of how the structure of DNA was discovered is nearly as fascinating as the structure itself. James Watson provides a first-person account in his book The Double Helix. It begins
I have never seen Francis Crick in a modest mood. Perhaps in other company he is that way, but I have never had reason so to judge him. It has nothing to do with his present fame. Already he is much talked about, usually with reverence, and someday he may be considered in the category of Rutherford or Bohr. But this was not true when, in the fall of 1951, I came to the Cavendish Laboratory of Cambridge University to join a small group of physicists and chemists working on the three-dimensional structures of proteins. At that time he was thirty-five, yet almost totally unknown. Although some of his closest colleagues realized the value of his quick, penetrating mind and frequently sought his advice, he was often not appreciated, and most people thought he talked too much.
The Eighth Day of Creation, by Horace Freeland Judson, superimposed on Intermediate Physics for Medicine and Biology.
The Eighth Day of Creation,
by Horace Freeland Judson.
The Double Helix is one of those iconic books that everyone should read for the insights it provides into how science is done, and for what is simply a fascinating story. The tale is also told from a more unbiased perspective in Horace Freeland Judson’s book The Eighth Day of Creation: The Makers of the Revolution in Biology. Let us end with Judson’s discussion of Watson and Crick’s now 60-year-old letter.
The letter to Nature appeared in the April 25th issue. To those of its readers who were close to the questions, and who had not already heard the news, the letter must have gone off like a string of depth charges in a calm sea. “We wish to suggest a structure for the salt of deoxyribose nucleic acid (D.N.A.). This structure has novel features which are of considerable biological interest,” the letter began; at the end, “It has not escaped our notice that the specific pairing we have postulated immediately suggests a possible copying mechanism for the genetic material.” The last sentence has been called one of the most coy statements in the literature of science.

Friday, April 19, 2013

Hyperpolarized 129Xe MRI of the Human Lung

Chapter 18 of the 4th edition of Intermediate Physics for Medicine and Biology is devoted to magnetic resonance imaging. Russ Hobbie and I discuss many aspects of MRI, including functional MRI and diffusion tensor imaging. One topic we do not cover is using hyperpolarized spins to study lung function. Fortunately, a clearly written review article, “Hyperpolarized 129Xe MRI of the Human Lung,” by John Mugler and Talissa Altes, appeared recently in the Journal of Magnetic Resonance Imaging (Volume 37, Pages 313–331, 2013). Below I reproduce excerpts from their introduction, with citations removed.
CLINICAL MAGNETIC RESONANCE IMAGING (MRI) of the lung is challenging due to the lung’s low proton density, which is roughly one-third that of muscle, and the inhomogeneous magnetic environment within the lung created by numerous air–tissue interfaces, which lead to a T2* value on the order of 1 msec at 1.5 T. Although advances continue with techniques such as ultrashort echo time (UTE) imaging of the lung parenchyma, conventional proton-based MRI is at a fundamental disadvantage for pulmonary applications because it cannot directly image the lung airspaces. This disadvantage of proton MRI can be overcome by turning to a gaseous contrast agent, such as the noble gas helium-3 (3He) or xenon-129 (129Xe), which upon inhalation permits direct visualization of lung airspaces in an MR image. With these agents, the low density of gas compared to that of solid tissue can be compensated by using a standalone, laser-based polarization device to increase the nuclear polarization to be roughly five orders of magnitude (10,000 times) higher than the corresponding thermal equilibrium polarization would be in the magnet of a clinical MR scanner. As a result, the MR signal from these hyperpolarized noble gases is increased by a proportionate amount and is easily detected using an MR scanner tuned to the appropriate resonance frequency. MRI using hyperpolarized gases has led to the development of numerous unique strategies for evaluating the structure and function of the human lung that provide advantages relative to current clinically available methods. For example, as compared with nuclear-medicine ventilation scintigraphy scans using 133Xe or aerosolized technetium-99m DTPA, hyperpolarized-gas ventilation MR images provide improved temporal and spatial resolution, expose the patient to no ionizing radiation, and can be repeated multiple times in a single day if desired. Although inhaled xenon has also been used as a contrast agent with computed tomography (CT), which can provide high spatial and temporal resolution, the high radiation dose and low contrast on the resulting ventilation images has dampened enthusiasm for the CT-based technique.

Although the first hyperpolarized-gas MR images were obtained using hyperpolarized 129Xe, and images of the human lung were acquired with hyperpolarized 129Xe only a few years later, the vast majority of work in humans has been performed using hyperpolarized 3He instead. This occurred primarily because 3He provided a stronger MR signal, due to its larger nuclear magnetic moment (and hence larger gyromagnetic ratio) compared to 129Xe and historically high levels of polarization (greater than 30%) achieved for 3He, and because there are no significant safety concerns associated with inhaled helium. However, in the years following the terrorist attacks of 9/11 there was a surge in demand for 3He for use in neutron detectors for port and border security, and this demand far exceeded the replenishment rate from the primary source, the decay of tritium used in nuclear warheads. As a result, 3He prices skyrocketed and availability plummeted. Currently, the U.S. government is regulating the supply of 3He, allocating this precious resource among users whose research or applications depend on 3He’s unique physical properties. This includes an annual allocation for medical imaging, which allows research on hyperpolarized 3He MRI of the lung to continue. Nonetheless, unless a new source for 3He is found it is clear that insufficient 3He is available to permit hyperpolarized 3He MRI of the lung to translate from the research community to a clinical tool.

In contrast to 3He, 129Xe is naturally abundant on Earth and its cost is relatively low. Thus, 129Xe is the obvious potential alternative to 3He as an inhaled contrast agent for MRI of the lung. While the 3He availability crisis has accelerated efforts to develop and evaluate hyperpolarized 129Xe for human applications, it is important to understand that 129Xe is not just a lower-signal alternative to 3He, forced upon us by practical necessity. In particular, the relatively high solubility of xenon in biological tissues and an exquisite sensitivity to its environment, which results in an enormous range of chemical shifts upon solution, make hyperpolarized 129Xe particularly attractive for exploring certain characteristics of lung function, such as gas exchange and uptake, that cannot be accessed using hyperpolarized 3He. The quantitative characteristics of gas exchange and uptake are determined by parameters of physiologic relevance, including the thickness of the blood–gas barrier, and thus measurements that quantify this process offer a potential wealth of information on the functional status of the healthy and diseased lung.

Historically, polarization levels for liter-quantities of hyperpolarized 129Xe have been roughly 10%, while those for similar quantities of hyperpolarized 3He have been greater than 30%. (Recall that the thermal equilibrium polarization of water protons at 1.5T is 0.0005%—four to five orders of magnitude lower.) Given 129Xe’s lower nuclear magnetic moment, this situation has put hyperpolarized 129Xe at a distinct disadvantage relative to 3He. A recent, key advance for 129Xe is the development of systems that can deliver liter quantities of hyperpolarized 129Xe with polarization on the order 50%. This now puts 129Xe on a competitive footing with 3He, positioning MRI of the human lung using hyperpolarized 129Xe to advance quickly in the immediate future, and making hyperpolarized 129Xe MRI of interest to the broader radiology and medical-imaging communities.
This idea of gas hyperpolarization is fascinating. How does one hyperpolarize the gas? Mugler and Altes explain:
Although it is possible to image either 129Xe or 3He by simply placing the gas (in a suitable container) in the magnet of an MR scanner, the low density of gas compared to that of solid tissue results in a signal that is too low to be of practical use for imaging the human lung… Nonetheless, the nuclear polarization can be increased dramatically compared to that produced by the magnet of the MR scanner by using a method called opticalpumping and spin exchange (OPSE), which was originally developed for nuclear-physics experiments many years before being applied to medical imaging.

As its name implies, OPSE is, in concept, a two-step process. The first step, optical pumping, involves using a laser to generate electron-spin polarization in a vapor of an alkali metal. This process takes place within a glass container, called an optical cell…positioned within a magnetic field... A small amount of the alkali metal, typically rubidium, is placed in the cell, which is heated…during the polarization process to create rubidium vapor. The optical cell is illuminated with circularly polarized laser light…at a specific wavelength (795 nm) to optically pump the rubidium atoms. This pumping preferentially populates one of the two spin states for the valence electron, thereby polarizing the associated electron spins and resulting in electron-spin polarization approaching 100%. In the second step of OPSE, collisions between spin-polarized rubidium atoms and noble-gas (129Xe or 3He) atoms within the cell result in spin exchange—the transfer of polarization from rubidium electrons to noble-gas nuclei...
To learn more, you can hear John Mugler discuss hyperpolarized gas MRI in the lung on youtube.

John Mugler discusses hyperpolarized gas MRI in the lung.

Friday, April 12, 2013

Radon

The largest source of natural background radiation is radon gas. In Chapter 17 of the 4th edition of Intermediate Physics for Medicine and Biology, Russ Hobbie and I discuss radon.
The naturally occurring radioactive nuclei are either produced continuously by cosmic γ ray bombardment, or they are the products in a decay chain from a nucleus whose half-life is comparable to the age of the earth. Otherwise they would have already decayed. There are three naturally occurring radioactive decay chains near the high-Z end of the periodic table. One of these is the decay products from 238U, shown in Fig. 17.27. The halflife of 238U is 4.5 × 109 yr, which is about the same as the age of the earth. A series of α and β decays lead to radium, 226Ra, which undergoes α decay with a half-life of 1620 yr to radon, 222Rn.

Uranium, and therefore radium and radon, are present in most rocks and soil. Radon, a noble gas, percolates through grainy rocks and soil and enters the air and water in different concentrations. Although radon is a noble gas, its decay products have different chemical properties and attach to dust or aerosol droplets which can collect in the lungs. High levels of radon products in the lungs have been shown by both epidemiological studies of uranium miners and by animal studies to cause lung cancer.
In Chapter 16 we consider radon in the context of the risk of the general population to low levels of background radiation.
The question of a hormetic effect or a threshold effect [as opposed to the linear no-threshold model of radiation exposure] has received a great deal of attention for the case of radon, where remediation at fairly low radon levels has been proposed. Radon is produced naturally in many types of rock. It is a noble gas, but its radioactive decay products can become lodged in the lung. An excess of lung cancer has been well documented in uranium miners, who have been exposed to fairly high radon concentrations as well as high dust levels and tobacco smoke. Radon at lower concentrations seeps from the soil into buildings and contributes up to 55% of the exposure to the general population.
Building Blocks of the Universe, by Isaac Asimov, superimposed on Intermediate Physics for Medicine and Biology.
Building Blocks of the Universe,
by Isaac Asimov.
Given the high profile of radon in our book, I thought readers might be interested in a bit of the history of this element. A brief discussion of the discovery of radon can be found in Isaac Asimov’s book Building Blocks of the Universe. After Asimov describes the Curies’ discoveries of radium and polonium from uranium ore in the late 1890s, he writes
When the radium atom breaks up, it forms an atom of radon, element No. 86. Radon is a gas, a radioactive gas! It fits into the inert gas column of the periodic table, right under xenon, and has all the chemical characteristics of the other inert gases.

Radon was first discovered in 1900 by a chemist named F. E. Dorn, and he called it radium emanation because it emanated from (that is, was given off by) radium. [William] Ramsay and R. Whytlaw-Gray collected the gas in 1908, and they called it niton from a Greek word meaning “shining.” In 1923, though, the official name became “radon” to show that the gas arose from radium…

Other gases arise from the breakdown of thorium and actinium … and have been called thoron and actinon, respectively. These are, as it turns out, varieties [isotopes] of radon. However, there have been suggestions that the element be named emanon (from “emanation”) since it does not arise from the breakdown of radium only, as “radon” implies.
Asimov's Biographical Encyclopedia of Science and Technology, by Isaac Asimov, superimposed on Intermediate Physics for Medicine and Biology.
Asimov's Biographical Encyclopedia
of Science and Technology,
by Isaac Asimov.
Asimov’s Biographical Encyclopedia of Science and Technology describes the scientist who discovered radon in more detail.
Dorn, Friedrich Ernst, German physicist
Born: Guttstadt (now Dobre, Miasto, Poland), East Prussia, July 27, 1848
Died: Halle, June 13, 1916

Dorn was educated at the University of Konigsberg and taught physics at the universities of Darmstadt and Halle. He turned to the study of radioactivity in the wake of Madame Curie’s discoveries and in 1900 showed that radium not only produced radioactive radiations, but also gave off a gas that was itself radioactive.

This gas eventually received the name radon and turned out to be the final member of Ramsay’s family of inert gases. It was the first clear-cut demonstration that in the process of giving off radioactive radiation, one element was transmuted (shades of the alchemists!) to another. This concept was carried further by Boltwood and Soddy.

Friday, April 5, 2013

Leon Glass wins Winfree Prize

From Clocks to Chaos, by Glass and Mackey, superimposed on Intermediate Physics for Medicine and Biology.
From Clocks to Chaos,
by Glass and Mackey.
Leon Glass was honored recently by the Society for Mathematical Biology. Their website states
The Society for Mathematical Biology is pleased to announce that this year’s recipient of the Arthur T. Winfree prize is Prof. Leon Glass of McGill University. Awarded every two years to a scientist whose work has “led to significant new biological understanding affecting observation/experiments,” this prize commemorates the creativity, imagination and intellectual breadth of Arthur T. Winfree.

Beginning with simple and brilliantly chosen experiments, Leon launched the study of chaos in biology. Among the applications he and his many collaborators and students pursued was the novel idea of “dynamical disease” and the better understanding of pathologies like Parkinson’s disease and cardiac arrhythmias. His elegant work (with Michael Guevara and Alvin Shrier) on periodic stimulation of heart cells demonstrated and explained how the interaction of nonlinearities with oscillations create complex dynamics and chaos.

The book From Clocks to Chaos, which he co-authored with Michael Mackey, was an instant classic that illuminated this difficult subject for a whole generation of mathematical biologists. His combination of imagination, experimental and mathematical insight, and ability to communicate fundamental principles has launched new fields of research and inspired researchers ranging from applied mathematicians to medical researchers.
Leon Glass is the Isadore Rosenfeld Chair in Cardiology at McGill University. Russ Hobbie and I cite From Clocks to Chaos (discussed previously in this blog) in the 4th edition of Intermediate Physics for Medicine and Biology, especially in Chapter 10 when discussing nonlinear dynamics. According to Google Scholar, the book has been cited 1800 times. Even more highly cited (over 2600 times) is Mackey and Glass’s paper “Oscillation and Chaos in Physiological Control Systems” (Science, Volume 197, Pages 287–289, 1977), which Russ and I also cite.

Other books and papers mentioned in IPMB include
Bub, G., A. Shrier, and L. Glass (2002) “Spiral wavegeneration in heterogeneous excitable media,” Phys. Rev. Lett., Volume 88, Article Number 058101.

Glass, L., Y. Nagai, K. Hall, M. Talajic, and S. Nattel (2002) “Predicting the entrainment of reentrant cardiacwaves using phase resetting curves,” Phys. Rev. E, Volume 65, Article Number 021908.

Guevara, M. R., L. Glass, and A. Shrier (1981) “Phaselocking,period-doubling bifurcations and irregular dynamicsin periodically stimulated cardiac cells,” Science Volume 214, Pages 1350–1353.

Glass, L. (2001) “Synchronization and rhythmicprocesses in physiology,” Nature, Volume 410, Pages 277–284.

Kaplan, D., and L. Glass (1995) Understanding NonlinearDynamics. New York, Springer-Verlag.
You can listen to Glass talk about cardiac arrhythmias below.

Leon Glass talks about cardiac arrhythmias.

Friday, March 29, 2013

1932: A Watershed Year in Nuclear Physics

I should have posted this article last year, to mark the 80th anniversary of the annus mirabilis of nuclear physics. Unfortunately, I didn’t think of it until I read “1932: A Watershed Year in Nuclear Physics” by Joseph Reader and Charles Clark, which appeared in the March 2013 issue of Physics Today. The article describes four major discoveries that changed nuclear physics forever.

 

Deuterium

The first landmark result was published by Harold Urey on January 1, 1932, in which he reported his discovery of deuterium, the isotope 2H. Russ Hobbie and I mention deuterium in Homework Problem 45 of Chapter 4 in the 4th edition of Intermediate Physics for Medicine and Biology.
Problem 45 Using the definitions in Problem 44, write the diffusion constant in terms of λ and vrms. By how much do you expect the diffusion constant for “heavy water” (water in which the two hydrogen atoms are deuterium, 2H) to differ from the diffusion constant for water? Assume the mean free path is independent of mass.
Unlike many elements, for which the various stable isotopes differ in mass be only a few percent, deuterium has twice the mass of normal hydrogen. Even when deuterium is incorporated into water, the H2O molecule’s mass increases by a significant 11%. Heavy water has been used as a non-radioactive biological tracer.

The Neutron

A second advance of 1932, and in my opinion the most important, is the discovery of the neutron by James Chadwick. Of course, the idea of a neutron is central to nuclear physics, and you cannot make sense of isotopes without it (I wonder how Urey interpreted the deuterium before the neutron was discovered). Russ and I discuss neutrons throughout Chapter 17 on Nuclear Medicine, and in particular we describe boron neutron capture therapy in Chapter 16
Boron neutron capture therapy (BNCT) is based on a nuclear reaction which occurs when the stable isotope 10B is irradiated with neutrons, leading to the nuclear reaction (in the notation of Chapter 17)
105B + 10n → 42α + 73Li
... Both the alpha particle and lithium are heavily ionizing and travel only about one cell diameter. BNCT has been tried since the 1950s; success requires boron-containing drugs that accumulate in the tumor.

The Positron

Discovery number three is the positron, the first example of antimatter. Carl Anderson found evidence of this positive electron in cosmic ray tracks in cloud chambers. Positrons appear in IPMB in two key places. In Chapter 15 (The Interaction of Photons and Charged Particles with Matter) positrons are key to pair production.
A photon with energy above 1.02 MeV can produce a particle–antiparticle pair: a negative electron and a positive electron or positron… Since the rest energy (mec2) of an electron or positron is 0.51 MeV, pair production is energetically impossible for photons below 2mec2 = 1.02 MeV.
The positron appears again in our discussion of β+ decay in Chapter 17.
Two modes of decay allow a nucleus to approach the stable line. In beta or electron) decay, a neutron is converted into a proton. This keeps A [mass number] constant, lowering N [neutron number] by one and raising Z [atomic number] by one. In positron (β+) decay, a proton is converted into a neutron. Again A remains unchanged, Z decreases and N increases by 1. We find β+ decay for nuclei above the line of stability and β- decay for nuclei below the line.
Isotopes that undergo β+ decay are used in positron emission tomography.
If a positron emitter is used as the radionuclide, the positron comes to rest and annihilates an electron, emitting two annihilation photons back to back. In positron emission tomography (PET) these are detected in coincidence….

PET can provide a functional image with information about metabolic activity. A very common positron agent is 18F fluorodeoxyglucoseglucose in which a hydroxyl group has been replaced with 18F. The PET signal is largest in those cells that have taken up the 18F because they are actively metabolizing glucose. PET has become particularly important in studies of brain function, where active neurons are detected by an increase in their metabolism, and in locating metastatic cancer.

Accelerators

The last of the four great developments of 1932 is the first use of accelerators to study nuclear reactions. John Cockcroft and Ernest Walton built an accelerator to produce high energy protons, which smashed into 7Li to produce two alpha particles. Their work was soon followed by the invention of the cyclotron by Ernest Lawrence, which is now the main tool for producing the unstable isotopes used in PET. Russ and I explain that
Positron emitters are short-lived, and it is necessary to have a cyclotron for producing them in or near the hospital. This is proving to be less of a problem than initially imagined. Commercial cyclotron facilities deliver isotopes to a number of nearby hospitals.
The Making of the Atomic Bomb, by Richard Rhodes. superimposed on Intermediate Physics for Medicine and Biology.
The Making of the Atomic Bomb,
by Richard Rhodes.

Soon after the miraculous year of 1932 Hitler came to power in Germany, and nuclear physics became much more than a scientific curiosity. The story of how the discoveries of Urey, Chadwick, Anderson, Cockcroft and Walton led relentlessly to the Manhattan Project is told masterfully in Richard Rhodes’ book The Making of the Atomic Bomb.

I have a few personal connections to this watershed year. My father Ron Roth, now retired and living in Lenexa Kansas, was born in 1932, proving that we are not so far removed from that historic time. In addition, my academic genealogy goes back to James Chadwick and Ernest Rutherford (whose lab Cockcroft and Walton worked in). Finally, Carl Anderson worked under the supervision of Robert Millikan, who was born in Morrison, Illinois, the small town I grew up in.