Monday, May 11, 2020

Are Tens of Thousands Dying from Radon Each Year or Not?

Radon Action Month Proclamation
by Michigan Governor Gretchen Whitmer.
Governor Whitmer declared last January to be Radon Action Month in Michigan. Just how serious of a health hazard is radon?

Russ Hobbie and I discuss radon in Chapter 16 of Intermediate Physics for Medicine and Biology.
Radon is produced naturally in many types of rock. It is a noble gas, but its 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 soil into buildings and contributes a large fraction of the exposure to the general population.
The Environmental Protection Agency published A Citizen’s Guide to Radon: The Guide to Protecting Yourself and Your Family from Radon. It recommends that you fix your home if your radon level is greater than 4 pCi/L (a picocurie per liter is equal to 37 decays per second per cubic meter of air).

Let’s put that into perspective. According to EPA’s 2003 Assessment of Risks from Radon in Homes, radon causes about 21,000 lung cancer deaths per year.
Based on its analysis, EPA estimates that out of a total of 157,400 lung cancer deaths nationally in 1995, 21,100 (13.4%) were radon related. Among NS [non-smokers], an estimated 26% were radon related... The estimated risks from lifetime exposure at the 4 pCi/L action level are: 2.3% for the entire population, 4.1% for ES [smokers], and 0.73% for NS. A Monte Carlo uncertainty analysis that accounts for only those factors that can be quantified without relying too heavily on expert opinion indicates that estimates for the U.S. population and ES may be accurate to within factors of about 2 or 3.
The data listed above are 25 years old. The American Cancer Society estimates that in 2020 about 136,000 deaths in the US will be from lung cancer. The reduction compared to 1995 is probably because fewer people smoke. I don’t see any way to accurately estimate the current number of radon-related deaths in the USA, but 20,000 per year is a reasonable guess based on the EPA’s 2003 Assessment.

Does this number make sense? Assume a 5% excess risk of dying from cancer per 1 Sievert of radiation dose. If we use the background dose from radon given in IPMB of about 2 mSv/year, then the excess risk is 0.0001/year. The current population of the US is about 330 million. Multiplying 330,000,000 times 0.0001 gives 33,000. This is the same order of magnitude as our 20,000 ballpark guess. Both of these estimates are uncertain, but they suggest that a few tens of thousands of deaths in the US each year are caused by radon. While this is not as bad as the coronavirus (80,000 deaths in a couple months), it’s still worrisome.

Tens of thousands dead. Really? Such estimates are based on the controversial linear-no-threshold model. A 2016 open-access article “Rectifying Radon’s Record: An Open Challenge to the EPA” by Jeffry Siegel, Charles Pennington, Bill Sacks, and James Welsh (International Journal of Radiology and Imaging Technology, Volume 2, Article Number 014) states
The American Lung Association has recently led a national workgroup to develop The National Radon Action Plan: A Strategy for Saving Lives. The U.S. Environmental Protection Agency (EPA) is the lead governmental organization projected to implement this plan. The stated intent of the plan is to address the “radon problem” in the United States, with the aim of saving 3,200 lives by the year 2020 through preventing at least a portion of the lung cancer mortality that is assumed to arise from inhaling modest doses of radon in homes, offices, and buildings. The plan identifies a number of actions that government can take in the spirit of saving lives by avoiding the inhalation of radon and its progeny. We are among a growing number of investigators who recognize the substantial body of evidence demonstrating that the radiation doses associated with indoor radon inhalation are not harmful. Radon, at these doses, is unlikely to be a cause of lung cancer, and, on the contrary, may be beneficial in various ways, including its paradoxical tendency to protect against lung cancer. In the present paper, we review and critique the past policies of the EPA with respect to indoor radon and the very impetus for the plan. We indicate that the plan should not be implemented because a preponderance of the evidence indicates an unintended consequence: implementation of the plan is likely to increase, rather than decrease, the risk of lung cancer.
What are we to believe? I don’t know the real risk of radon exposure. One thing I do know is that we need to figure out whether or not the linear-no-threshold model is correct. Are tens of thousands of our citizens dying each year from radon exposure? It seems to me that with so many lives at stake, our nation needs to invest the time and money necessary to answer this question. Either the EPA overestimates the risk, in which case we can focus on other more pressing issues, or it accurately estimates the risk, in which case we have an epidemic on our hands.

Friday, May 8, 2020

Ping Pong

A ping pong paddle on top of Intermediate Physics for Medicine and Biology.
When I was in a teenager, I played a lot of ping pong. It was the sport that I played best (which isn’t saying much) and loved most (except for baseball). Dad bought a ping-pong table for the family when I was young, and he was my first opponent. As I grew up, I spent hours practicing with friends, perfecting our serves, slices, and slams. In high school, my friend Terry Fife and I played hundreds of games, and we were evenly matched. The psychological warfare during those battles was fierce.

I continued playing ping pong in college. I bought one of those paddles that has foam under the rubber so the ball stays in contact with the surface for a long time, letting you impart more spin. During my first two years at the University of Kansas, I studied but also spent a lot of time playing ping pong at the dorm. I was better than most of the guys, although my friend Son Do could beat me, much to my chagrin. Too many evenings were squandered playing cards or ping pong. Only during my last couple years at KU did I get serious about physics.

Life in Moving Fluids,
by Steven Vogel.
Understanding ping pong requires us to discuss the Magnus effect, which Steven Vogel illustrates by considering a spinning cylinder in a flowing fluid. In Life in Moving Fluids—a book often cited in Intermediate Physics for Medicine and Biology—Vogel writes
A look at the resulting streamlines (Figure 10.12 [a modified version of which is shown below]) clarifies what’s happening in this superposition of rotation and translation of a cylinder. On one side of the cylinder the two motions in the fluid oppose one another, so the velocities are lower and the streamlines are farther apart. On the other side, the motions are additive, velocities are increased, and the streamlines are closer together. By Bernoulli’s principle pressure will be elevated on the side where flow speeds are lower and will be reduced on the side where the speeds are higher. Thus a net pressure or force will act in a direction normal to the free-speed flow—in short, lift
Figure 10.12. If a solid body such as a cylinder rotates as it translates through a fluid, the resulting asymmetry of flow generates a force normal to the free-stream flow. We call this force lift. Adapted from Life in Moving Fluids by Steven Vogel.
This phenomenon, the lift of a rotating cylinder moving through a fluid, is called the “Magnus effect,” after H. G. Magnus (1802-1870).

The Magnus effect (at a little lower intensity) works for spheres as well as for cylinders. It’s a really big deal in sports in which spheres are thrown, hit, or otherwise put into motion since (except for a golf slice) a confusingly nonstraight course is distinctly meritorious. Two pleasant books on such contemporary compulsions are Sport Science, by P. J. Brancazio (1984), with good references, and The Physics of Baseball, by R. K. Adair (1990), with more on the Magnus effect specifically.
Me playing ping pong at the
University of Kansas, circa 1980.
Vogel’s description of the Magnus effect is clear and interesting, but I think it’s not very relevant. In ping pong, the goal of spin is not to make the ball follow a curved path through the air, but instead is to make the ball leave your opponent’s paddle traveling in the wrong direction! After putting backspin on the ball, I loved to see Terry Fife’s return dive into the net. And I enjoyed nothing more than giving the ball some sidespin, and then watching Son Do’s next hit miss the table on the right. Those were the days!

Thursday, May 7, 2020

Intensity-Modulated Radiation Therapy

Radiation Oncology: A Physicist's-Eye View,
by Michael Goitein.

In an earlier post, I mentioned Michael Goitein textbook Radiation Oncology: A Physicist’s Eye View. In Intermediate Physics for Medicine and Biology, Russ Hobbie and I cite Goitein’s book when we discuss intensity-modulated radiation therapy.
In classical radiotherapy, the beam was either of uniform fluence across the field, or it was shaped by an attenuating wedge placed in the field. Intensity-modulated radiation therapy (IMRT) is achieved by stepping the collimator leaves during exposure so that the fluence varies from square to square in Fig. 16.45 (Goitein 2008; Khan 2010, Ch. 20)

It was originally hoped that CT reconstruction techniques could be used to determine the collimator settings at different angles. This does not work because it is impossible to make the filtered radiation field negative, as the CT reconstruction would demand. IMRT with conventional treatment planning improves the dose pattern (Goitein (2008); Yu et al. (2008)), providing better sparing of adjacent normal tissue and allowing a boost in dose to the tumor.
Goitein’s analysis of IMRT is full of insight, and I quote an excerpt below. Enjoy.

INTENSITY-MODULATED RADIATION THERAPY (IMRT)

So far, we have implicitly assumed that each radiation field is near-uniform over its cross section; dose uniformity of a field within the target volume has, in fact, historically been an explicit goal of radiotherapy. However, the radical suggestion to allow the use of non-uniform fields was made some two decades ago, independently by Anders Brahme and Alan Cormack, fresh from co-inventing the CT scanner—and, in the context of π-meson therapy... (Cormack 1987; Brahme 1988; Pedroni 1981). Their idea was based on the judgment that, using mathematical techniques, an irradiation scheme using non-uniform beams could be found which would more closely achieve the ideal of delivering the desired dose to the target volume while limiting the dose to the normal tissues outside the target volume to some predefined value.

Brahme’s and Cormack’s approaches were motivated by the observation that, in CT reconstruction, one can deduce from the intensity reduction of X-rays traversing an object along a series of straight paths what the internal structure of the object is. By inverting the mathematics, one can deduce the intensities ( pencil beam weights) of a series of very small beams (pencil beams) that pass through the object and deliver dose within it. This procedure leads to highly nonuniform individual fields which, in combination, deliver the desired (usually, uniform) dose to the target volume.

There are two very substantial flaws to the original idea. The first is that, when the problem is posed to deliver zero dose outside the target volume as was initially proposed, many of the computed intensitiesare negative—a highly unphysical result. The second is that there is no a priori way of specifying a physically possible dose distribution to serve as the goal of the optimization.

However, the basic idea of using non-uniform beams has proven enormously fruitful. A workable computational solution is to use optimization algorithms to iteratively adjust the pencil beam weights such that the resulting dose distribution maximizes some score function. The search is computationally intensive and therefore poses interesting technical challenges. However, the still bigger challenge is to find score functions which give a viable measure of clinical goodness. Increasingly, biophysical models of the dose-response of both tumors and normal tissues are being investigated and are beginning to be used as elements of such score functions. These matters are discussed in Chapters 5 and 9.

Intensity-modulated radiation therapy (IMRT), as treatments featuring non-uniform beams are called, has been most intensely developed for X-ray therapy. However, it is equally appropriate for other radiation modalities—including protons. With charged particles one has an extra degree of freedom. One can vary the beam intensity as a function of lateral position and as a function of penetration (energy).
Who was Michael Goitein? Below is an excerpt from his 2017 obituary in the International Journal of Radiation Oncology Biology Physics.
Michael Goitein was a visionary thought provocateur and is rightly judged to have been an exceptionally innovative and creative physicist in radiation oncology (Fig. 1). He was a critical player in the development of proton radiation therapy, with many of his advances widely used in current proton and photon therapy. His highly important and numerous contributions to medical science have been well rewarded with many important awards: Fulbright Fellowship, 1961-65; US Research Career Development Award, 1976-81; Fellow, American Association of Physicists in Medicine, 2000; Gold Medal of the American Society for Radiation Oncology, 2003; and Lifetime Achievement Award of the European Society for Radiotherapy and Oncology, 2014. Significantly, he was a cofounder of the Proton Therapy Co-Operative Oncology Group and served as the second president. Furthermore, he participated in many of its subsequent functions. In addition, Michael was an invited lecturer at a long list of national and international conferences and medical centers. He published 163 articles and 3 books that have positively affected the practice of radiation oncology (1). It must also be mentioned that he has authored 4 books of a personal nature, one of which is a book of his poems (2, 3, 4, 5).
To learn more, read his ASTRO interview.

Wednesday, May 6, 2020

Never at Rest

Never at Rest, by Richard Westfall.
Isaac Newton’s name appears many times in Intermediate Physics for Medicine and Biology. You can learn more about him in Richard Westfall’s wonderful book Never at Rest: A Biography of Isaac Newton. As we all sit in quarantine because of the coronavirus pandemic, I thought you might like to read about Newton’s experience with the plague. Here is an excerpt from Never at Rest.
In the summer of 1665, a disaster descended on many parts of England including Cambridge. It had “pleased Almighty God in his just severity,” as Emmanuel College put it, “to visit this towne of Cambridge with the plague of pestilence.” Although Cambridge could not know it and did little in the following years to appease divine severity, the two-year visitation was the last time God would choose to chastise them in this manner [until 2020]. On 1 September, the city government canceled Sturbridge Fair [one of the largest fairs in Europe] and prohibited all public meetings. On 10 October, the senate of the university discontinued sermons at Great St. Mary’s and exercises in the public schools. In fact, the colleges had packed up and dispersed long before. Trinity [College, Cambridge] recorded a conclusion on 7 August that “all Fellows & Scholars which now go into the Country upon occasion of the Pestilence shall be allowed [the] usual Rates for their Commons for [the] space of [the] month following...” For eight months the university was nearly deserted…

Many of the students attempted to continue organized study by moving with their tutors to some neighboring village. Since Newton was entirely independent in his studies and had had his independence confirmed with a recent B.A. [Newton received his bachelors degree in August 1665],… he returned… to Woolsthorpe [the Newton family home]...

Much has been made of the plague years in Newton’s life. He mentioned them in his account of mathematics. The story of the apple [hitting him on the head, triggering the discovery of the universal law of gravity], set in the country, implies the stay in Woolsthorpe. In another much-quoted statement written in connection with the calculus controversy [a debate between Newton and Leibniz about who first invented calculus] about fifty years later, Newton mentioned the plague years again.
In the beginning of the year 1665 I found the Method of approximating series & the Rule for reducing any dignity of any Binomial into such a series [the binomial theorem]. The same year in May I found the method of Tangents of Gregory & Slusius [a way of finding the slope of a curve], & in November had the direct method of fluxions [diferential calculus] & the next year in January had the Theory of Colours [later published in Opticks] & in May following I had entrance into [the] inverse method of fluxions [integral calculus]. And the same year I began to think of gravity extending to [the] orb of the Moon & (having found out how to estimate the force with [which a] globe revolving within a sphere presses the surface of the sphere) from Keplers rule of the periodical times of the Planets being in sesquialterate proportion [to the 3/2 power] of their distances from the center of their Orbs, I deduced that the forces [which] keep the Planets in their Orbs must [be] reciprocally as the squares of their distances from the centers about [which] they revolve [the inverse square law]: & thereby compared the force requisite to keep the Moon in her Orb with the force of gravity at the surface of the earth and found them answer pretty nearly. All this was in the two Plauge years of 1665-1666. For in those days I was in the prime of my age for invention and minded Mathematicks & Philosophy [physics] more than at any time since.
 So what are you doing while stuck at home during the coronavirus pandemic?

Tuesday, May 5, 2020

A Toy Model For Tomography

In Chapter 12 of Intermediate Physics for Medicine and Biology, Russ Hobbie and I discuss tomography. The algorithms for solving the tomography problem involve calculus, Fourier analysis, and convolutions. I love all that mathematics, but some people don’t (😮). Is there a way to introduce tomography to students who don’t have a mathematical background?

Let’s start with a simple object divided into six pixels. Our goal is to determine the value of some property for each pixel. If this were positron emission tomography, this property would be the concentration of a radioactive substance. If this were a CT scan, this property would be the x-ray absorption. How we interpret the property doesn’t matter; we’re just going to assign a number to each pixel.

The Froward Problem

Suppose this is your object.


We assume what you can measure is the sum of the pixels along one direction: a projection.

You can take projections from different orientations; sum the pixels in that direction.
A tomography machine in the hospital measures projections.

The Inverse Problem

So far we have examined the forward problem: determine the projections from the object. Next, consider the inverse problem: determine the image from the projections. Here’s another example.

How do you obtain an image of the object from multiple projections? That is the fundamental problem of tomography. In other words, how do you figure out what number to put into each pixel so that it gives the projections shown above? Stop reading and try to guess the image. When you’re done, continue reading.

Perhaps you found this image.
Good job, but your friends working from the same projections may have found different images.

Check for yourself; they all have the same three projections. If you permit negative values, you can find even more images.
The values don’t have to be integers.
The number of solutions is infinite. What can we do to pick the correct image? Use more projections!

You’ll find only one image consistent with all six projections: the first one I listed (try it yourself). If you have enough projections, you can find a unique image. You’ve solved the problem of tomography.

If you add noise to the projections, you may have no solution. In that case, you would need to use something like the least squares method to estimate the image. But that’s another story.

If you have only six pixels (and no noise) you can compute the image using trial and error, and some logic, sort of like Sudoku. A medical image, however, might have tens of thousands of pixels! What do you do then? That’s exactly what Russ and I discuss in Chapter 12 of Intermediate Physics for Medicine and Biology.

I will leave you with a final example to solve on your own. Enjoy!

Monday, May 4, 2020

IPMB Scavenger Hunt!

Here in Michigan Big Gretch has us all stuck at home cause of rona, and some of you might be quarantined with a younger brother. What can you do with him that is both educational and fun? Try the Intermediate Physics for Medicine and Biology Scavenger Hunt!

In the first section of IPMB, Russ Hobbie and I talk about distances and sizes. We show a couple illustrations (Figures 1.1 and 1.2) containing biological objects drawn to scale. Our goal is to help you build intuition about the relative size of things. This is an important skill, both for a sophomore premed student and for your third-grade niece. So, grab your little rug-rat and let’s play. Start with a picture of severe acute respiratory syndrome coronavirus, the virus that causes covid-19.

The starting point for playing IPMB Scavenger Hunt!
The starting point for playing IPMB Scavenger Hunt!
This game is best played in Powerpoint, so I suggest downloading the .pptx version of the figure. Now, you and your little sister search the internet for pictures of biological objects, copy and paste them into the figure, and label them. Make sure you expand or contract the image so its size matches the scale bars in the figure. You should be able to fit eight to ten objects in the picture before it gets crowded. When done, you will have created your own version of Fig. 1.1 or 1.2. But yours will be cooler that ours, because it will be in color and will probably contain real images as opposed to poorly drawn cartoons. An example that I created is shown below.

An example of playing Intermediate Physics for Medicine and Biology Scavenger Hunt.
An example resulting from playing IPMB Scavenger Hunt! A couple good sources of pictures are the CDC Image Library, and the Molecule of the Month Website.
What are the goals of this activity?
If your young cousin doesn’t want to stop, you can create your own version of the game using a larger or smaller initial object (perhaps a nerve axon and a 10 micron scale bar).

Your tenderfoot friend may select objects that don’t fit into the figure. For instance, she may want to include a red blood cell, which is too large for the slide. Don’t say no. Find an object around the house (perhaps a pillow) that is about the size of a red blood cell on this scale, thereby reinforcing the relative size of things. Don’t worry if you don’t have Powerpoint. Instead, pick your own size scale (for example, one foot equals one micron) and go around the house collecting everyday objects to represent “viruses,” “bacteria,” and “cells.” The only requirement is that they are all to scale.

If the little fella can’t get enough biology, have him look at some of David Goodsell’s drawings (be sure to see his latest coronavirus painting, below). If you have plenty of time and money, buy him a plush covid-19 toy.

Enjoy!
The coronavirus.
Illustration by David S. Goodsell, RCSB Protein Data Bank;
doi: 10.2210/rcsb_pdb/goodsell-gallery-019

Friday, May 1, 2020

Faust in Copenhagen

Faust in Copenhagen, by Gino Segrè, superimposed on Intermediate Physics for Medicine and Biology.
Faust in Copenhagen,
by Gino Segrè.
My friend and fellow biological physicist Gene Surdutovich loaned me his copy of Faust in Copenhagen: A Struggle for the Soul of Physics, by Gino Segrè. It’s about seven physicists—Niels Bohr, Max Delbrück, Paul Dirac, Paul Ehrenfest, Werner Heisenberg, Lise Meitner, and Wolfgang Pauli—who played a key role in the development of quantum mechanics in the 1920s. Segrè uses a 1932 meeting in Copenhagen as the hinge about which his story revolves.
This is a book about seven physicists, six men and one woman, who attended a small annual gathering in Copenhagen in April 1932. To be honest, only six of them were actually there. The seventh, Wolfgang Pauli, had originally intended to go, as he had in earlier years and would do so again, but he decided that spring instead to take a vacation. He was there in spirit, as you will see.
Of the seven, the one who plays the least significant role in the story is Delbrück. Nevertheless, he is the one most closely related to Intermediate Physics for Medicine and Biology. One topic I have examined in this blog is how scientists make the transition from physics to biology. Below I provide brief excerpts from Faust in Copenhagen, explaining how Delbrück did it. It appears in the book’s epilogue; I will add links and fill in background that you would have known had you read the rest of the book.

How Delbrück Became a Biologist

While working as the “family-theorist” in Meitner’s laboratory [in Berlin, she was the only experimental physicist of the seven], Delbrück often returned to see Bohr [his former mentor] and the continually changing group at Blegdamsvej [a street in Copenhagen where Bohr’s institute was located]. After the April 1932 meeting he had gone back to Berlin but in August was once again in Copenhagen. Sleepily stepping out of the carriage of the overnight train from Berlin, Delbrück was surprised to see Bohr’s trusted collaborator Leon Rosenfeld waiting for him with a message from Bohr. Delbrück was to go straightaway to the great meeting hall of the Rigsdag, the Danish parliament, where Bohr was about to deliver the opening lecture at the International Congress of Light Therapists… Knowing Delbrück’s interest in biology, Bohr wanted him to attend. This turned out to be the talk that changed Delbrück’s life [my italics].

In 1932 some biologists still clung to the notion that special forces not described by the ordinary laws of physics and chemistry were responsible for the existence of life… [In his talk, Bohr] asked: Could complementarity apply to life itself? [Complementarity was the principle that two opposing properties could both be true, but could not be observed simultaneously.] Perhaps the distinction between living and nonliving was not so easy to understand. Might the act of measurement be a critical step in the assessment [as it is in the Copenhagen interpretation of quantum mechanics]?...

Bohr thought there might be some intangible quality that differentiates living from nonlinving, a quality that could not be precisely quantified. He was not alluding to a divine spark, but rather to something analogous to the impossibility of saying with certainty whether light is a photon or a wave. Bohr’s conjecture was provocative, as it was meant to be, but in the end it turned out to be wrong. DNA and RNA are the answer to life, not complementarity.

But the potential connection between physics and biology had surfaced. Would Delbrück have become a biologist if he had not thought he might find something like complementarity at the root of life’s existence? Would molecular biology have developed as it did without the structure imposed on it by the discipline Delbrück had learned from Bohr? Many who know the field think the answer to both questions is no. When James Watson wrote Delbrück in 1953 that he and Francis Crick had found their double helix model of DNA, Delbrück was struck by the elegance of the structure, but disappointed by its simplicity. Hearing that life did not require any basic new principles, he remarked that if felt to him as if the hydrogen atom had been fully explained in the 1920s without the need for quantum mechanics. [The Bohr model explained the spectrum of hydrogen in 1913, but the more modern theory of quantum mechanics refined and improved that model in the 1920s.]

By then Delbrück had become a Bohr-life figure in the new field of molecular biology [Bohr, one of the senior members of the group of seven along with Ehrenfest and Meitner, was much loved as a mentor and father-figure as well as a scientist]. He had created at Cold Spring Harbor Laboratory on Long Island and at the California Institute of Technology Copenhagen-like atmospheres for young biologists…

[Delbrück], together with Salvador Luria and Alfred Hershey, with whom he shared the 1969 Nobel Prize in Physiology or Medicine, had shown how bacteriophages, the viruses that attack bacteria, were the simplest and most efficient tool for studying genetics… In many ways the importance of these findings were comparable to Bohr’s 1913 discovery of the rules for the hydrogen atom’s radiation...

Perhaps these qualities [of personality that are similar to Bohr’s, of integrity and intellectual openness] are the basis of the deep lifelong link between Pauli and Debrück. Luria, in his memoirs, remembered his first meeting with Delbrück. At the time both were recent refugees from Europe. It was New York City in late December 1941, only a few weeks after Pearl Harbor:
Max took me to dinner with two other scientists, one of them the great physicist Wolfgang Pauli. I was properly intimidated, but Pauli simply asked me “Sprechen Sie Deutsch?” and without waiting for a reply proceeded to eat and speak German so prodigiously fast that I understood not a word [of the seven, only Bohr (Danish) and Dirac (English) did not speak German as their first language]. I would have been even more scared had I known of Pauli’s classic remark “So young and he has already contributed so little.” [Pauli was known and loved for his quick and biting wit.]
Fortunately young Luria soon contributed a great deal.

As for Pauli, his last significant letter, written two months before he died, was a twelve-page one to Max. Beginning with a discussion of the festivities for Lise Meitner’s eightieth birthday, it ends on a more personal note: “I cannot forget the dear manner, which goes back to old times, when you took leave of me. I have the impression that here something has been renewed between us that is important.”

Delbrück would hike and camp with young biologists in the Southern California desert rather than walk with then along the North Sea shore of Denmark [where Bohr had walked with young Heisenberg, discussing quantum mechanics], but the spirit was the same as Bohr’s had been. On the other hand, rather than Bohr’s polite, “I don’t mean to criticize, only to understand,” Max remained famous for his rudeness... Though the spirit was Bohr’s, the style was Pauli’s...

Delbrück died in 1981 in California.
Looking for a good read while cooped up in your house because of the coronavirus? I recommend Faust in Copenhagen. It’s wunderbar.

Thursday, April 30, 2020

Radiation Oncology: A Physicist's Eye View

I miss Oakland University’s library. It’s been locked up now for about six weeks because of the coronavirus pandemic, so I can’t checkout books. As the physics department’s liaison to the library, I’ve used our meager annual allocation to purchase many of the books cited in Intermediate Physics for Medicine and Biology. I enjoy browsing through the stacks more than I realized.

I now better appreciate online access to textbooks. IPMB can be downloaded through the OU library catalog at no charge for anyone logging on as an OU student, faculty,or staff. I have access to other books online, and the number seems to be increasing every year. 

Radiation Oncology: A Physicist's Eye View, by Michael Goitein, superimposed on Intermediate Physics for Medicine and Biology.
Radiation Oncology:
A Physicist's Eye View
,
by Michael Goitein.

For example, I can download a textbook cited in Chapter 16 (Medical Uses of X-Rays) of IPMB.
Goitein M (2008) Radiation Oncology: A Physicist’s Eye View. Springer, New York.
It’s an excellent book, and you may hear more about it from me in the coming weeks. The preface begins:
This book describes how radiation is used in the treatment of cancer. It is written from a physicist’s perspective, describing the physical basis for radiation therapy, and does not address the medical rationale or clinical aspects of such treatments. Although the physics of radiation therapy is a technical subject, I have used, to the extent possible, non-technical language. My intention is to give my readers an overview of the broad issues and to whet their appetite for more detailed information, such as is available in textbooks.
“Whetting their appetite” is also a goal of IPMB, and of this blog. In fact, education in general is a process of whetting the appetite of readers so they will go learn more on their own.

I particularly like Michael Goitein’s chapter on uncertainty. The lessons it provides apply far beyond radiation oncology. In fact, it’s a lesson that is never more relevant than during a pandemic. Enjoy!

UNCERTAINTY MUST BE MADE EXPLICIT

ISO [International Organization for Standardization] (1995) states that “the result of a measurement [or calculation] … is complete only when accompanied by a statement of uncertainty.” Put more strongly, a measured or computed value which is not accompanied by an uncertainty estimate is meaningless. One simply does not know what to make of it. For reasons which I do not understand, and vehemently disapprove of, the statement of uncertainty in the clinical setting is very often absent. And, when one is given, it is usually unaccompanied by the qualifying information as to the confidence associated with the stated uncertainty interval—which largely invalidates the statement of uncertainty.

The importance of first estimating and then providing an estimate of uncertainty has led me to promulgate the following law:
Goitein's First Law of Uncertainty.
[Goitein's First Law of Uncertainty.]
There is simply no excuse for violating either part of Law number 1. The uncertainty estimate may be generic, based on past experience with similar problems; it may be a rough “back-of-the-envelope” calculation; or it may be the result of a detailed analysis of the particular measurement. Sometimes it will be sufficient to provide an umbrella statement such as “all doses have an associated confidence interval of ±2% (SD) unless otherwise noted.” In any event, the uncertainty estimate should never be implicit; it should be stated.

In graphical displays such as that of a dose distribution in a two- dimensional plane, the display of uncertainty can be quite challenging. This is for two reasons. First, it imposes an additional dimension of information which must somehow be graphically presented. And second because, in the case, for example, of the value of the dose at a point, the uncertainty may be expressed as either a numerical uncertainty in the dose value, or as a positional uncertainty in terms of the distance of closest approach. One approach to the display of dose uncertainty is shown in Figure 6.4 of Chapter 6.

HOW TO DEAL WITH UNCERTAINTY

To act in the face of uncertainty is to accept risk. Of course, deciding not to act is also an action, and equally involves risk. One’s decision as to what action to take, or not to take, should be based on the probability of a given consequence of the action and the importance of that consequence. In medical practice, it is particularly important that the importance assigned to a particular consequence is that of the patient, and not his or her physician. I know a clinician who makes major changes in his therapeutic strategy because of what I consider to be a trivial cosmetic problem. Of course, some patients might not find it trivial at all. So, since he assumes that all patients share his concern, I judge that he does not reflect the individual patient’s opinion very well. Parenthetically, it is impressive how illogically most of us perform our risk analyses, accepting substantial risks such as driving to the airport while refusing other, much smaller ones, such as flying to Paris (Wilson and Crouch, 2001). (I hasten to add that I speak here of the risk of flying, not that of being in Paris.)

People are often puzzled as to how to proceed once they have analyzed and appreciated the full range of factors which make a given value uncertain. How should one act in the face of the uncertainty? Luckily there is a simple answer to this conundrum, which is tantamount to a tautology. Even though it may be uncertain, the value that you should use for some quantity as a basis for action is your best estimate of that quantity. It’s as simple as that. You should plunge ahead, using the measured or estimated value as though it were the “truth”. There is no more correct approach; one has to act in accordance with the probabilities. To reinforce this point, here is my second law:
Goitein's Second Law of Uncertainty.
[Goitein's Second Law of Uncertainty.]
It may seem irresponsible to promote gambling when there are life-or-death matters for a patient at stake; the word has bad connotations. But in life, since almost everything is uncertain, we in fact gamble all the time. We assess probabilities, take into account the risks, and then act. We have no choice. We could not walk through a doorway if it were otherwise. And that is what we must do in the clinic, too. We cannot be immobilized by uncertainty. We must accept its inevitability and make the best judgment we can, given the state [of] our knowledge.

Wednesday, April 29, 2020

The Toroid Illustration (Fig. 8.26)

In Chapter 8 (Biomagnetism) of Intermediate Physics for Medicine and Biology, Russ Hobbie and I show an illustration of a nerve axon threaded through a magnetic toroid to measure its magnetic field (Fig. 8.26).
Fig. 8.26. A nerve cell preparation is threaded through the magnetic
toroid to measure the magnetic field. The changing magnetic flux in
the toroid induces an electromotive force in the winding. Any external
current that flows through the hole in the toroid diminishes the magnetic field.
While this figure is clear and correct, I wonder if we could do better? I started with a figure of a toroidal coil from a paper I published with my PhD advisor John Wikswo and his postdoc Frans Gielen.
Gielen FLH, Roth BJ, Wikswo JP Jr (1986) Capabilities of a Toroid-Amplifier System for Magnetic Measurements of Current in Biological Tissue. IEEE Trans. Biomed. Eng. 33:910-921.
Starting with Figure 1 from that paper (you can find a copy of that figure in a previous post), I modified it to resemble Fig. 8.26, but with a three-dimensional appearance. I also added color. The result is shown below.

An axon (purple) is threaded through a toroid to measure the magnetic field.
The toroid has a ferrite core (green) that is wound with insulated copper wire
(blue). It is then sealed in a coating of epoxy (pink). The entire preparation
is submerged in a saline bath. The changing magnetic flux in the ferrite induces
an electromotive force in the winding. Any current in the bath that flows
through the hole in the toroid diminishes the magnetic field.
Do you like it?

To learn more about how I wound the wire and applied the epoxy coating, see my earlier post about The Magnetic Field of a Single Axon. The part about "any current in the bath that flows through the hole in the toroid diminishes the magnetic field" is described in more detail in my post about the Bubble Experiment.

Tuesday, April 28, 2020

Chernobyl Then and Now: A Global Perspective

Last year I was supposed to give a talk at Oakland University for a symposium about “Chernobyl Then and Now: A Global Perspective.” It was part of an exhibition at the OU Art Gallery titled “McMillan’s Chernobyl: An Intimation of the Way the World Would End.” My role at the symposium was to explain the factors that led to the explosion of the Chernobyl nuclear power plant. I was chosen by the organizer, OU Professor of Art History Claude Baillargeon, because I had taught a class about The Making of the Atomic Bomb in the Oakland’s Honors College.

Readers of Intermediate Physics for Medicine and Biology should become familiar with the Chernobyl disaster because it illustrates how exposure to radiation can affect people over different time scales, from short term acute radiation sickness to long-term radiation-induced cancer.

It turned out I could not attend the symposium. My friend Gene Surdutovich stepped in at the last minute to replace me, and because he is from Ukraine—where the disaster occurred—he provided more insight than I could have. However, I thought the readers of this blog might want to read a transcript of the talk I planned to present. It was supposed to be my “TED Talk,” aimed at a broad audience with limited scientific background. No Powerpoint, no blackboard; just a few balls and a pencil as props.
The nuclear reactor in Chernobyl had an inherently unstable design that led to the worst nuclear accident in history. To understand why the design was so unstable, we need to review some physics.

The nucleus of an atom contains protons and neutrons. The number of protons determines what element you have. For instance, a nucleus with 92 protons is uranium. The number of neutrons determines the isotope. If a nucleus has 92 protons and 146 neutrons it is uranium-238 (because 92 + 146 = 238). Uranium-238 is the most common isotope of uranium (about 99% of natural uranium is uranium-238). If the nucleus has three fewer neutrons, that is only 143 neutrons instead of 146, it’s uranium-235, a rare isotope of uranium (about 1% of natural uranium is uranium-235).

No stable isotopes of uranium exist, but both uranium-235 and uranium-238 have very long half-lives (a half-life is how long it takes for half the nuclei to decay). Their half-lives are several billion years, which is about the same as the age of the earth. So many of the atoms of uranium that originally formed with the earth have not decayed away yet, and still exist in our rocks. We can use them as nuclear fuel.

Although uranium-235 is the rarer of the two long-lived isotopes, it is the one that is the fuel for a nuclear reactor. The uranium-235 nucleus is “fissile” meaning that it is so close to being unstable that a single neutron can trigger it to break in two pieces, releasing energy and two additional neutrons. This is called nuclear fission.

A nuclear chain reaction can start with a lot of uranium-235 and a single neutron. The neutron causes a uranium-235 nucleus to fission, breaking into two pieces plus releasing two additional neutrons and energy. These two neutrons hit two other uranium-235 nuclei, causing each of them to fission, releasing a total of four neutrons plus more energy. These four neutrons hit four other uranium-235 nuclei, releasing eight neutrons….and so on. The atomic bomb dropped on Hiroshima at the end of World War Two was based on just such an uncontrolled uranium-235 chain reaction. Fortunately, there are ways to control the chain reaction, so it can be used for more peaceful purposes, such as a nuclear reactor.

One surprising feature of uranium-235 is that SLOW neutrons are more likely to split the nucleus than FAST neutrons. How this effect was discovered is an interesting story. Enrico Fermi, an Italian physicist, was studying nuclear reactions in the 1930s by bombarding different materials with neutrons. He observed more nuclear reactions if his apparatus sat on a wooden table top than if it sat on a marble table top! What? It turns out wood was better at slowing the neutrons than marble. Think how confusing this must have been for Fermi. He was so confused that he tried submerging the apparatus in a pond behind the physics building and the reactions increased even more!

A uranium-235 chain reaction triggered by neutrons works best with slow neutrons. Therefore, nuclear reactors need a “moderator”: a substance that slows the neutrons down. The moderator is the key to understanding what happened at Chernobyl.

The best moderators are materials whose nuclear mass is about the same as the mass of a neutron. If the nucleus was a lot heavier than the neutron, the neutron would not slow down after the collision. Imagine this tennis ball is the light neutron, and this big basketball is the heavy nucleus. When the neutron hits the nucleus, it just bounces off. It changes direction but doesn’t slow down. Now, imagine this neutron collides with a very light particle, represented by this ping pong ball. When the relatively heavy neutron hits the light particle, it will just push it out of the way like a train hitting a mosquito. The neutron itself won’t slow down much. To be effective at slowing the neutron down, the nucleus needs to be about the same mass as the neutron. What has a similar mass to a neutron? A proton. What nucleus contains a single proton? Hydrogen. Watch what happens when a neutron and a hydrogen nucleus collide? This ball is the neutron, and this ball is the proton: the hydrogen nucleus. Right after the collision, the neutron stops! It is like when a moving billiard ball slams into a stationary billiard ball; the one that was moving stops, and the one that was stationary starts moving. Interacting with hydrogen is a great way to slow down neutrons. Therefore, hydrogen is a great moderator. Where do you find a lot of hydrogen? Water (H2O). It was the hydrogen in the water of the wooden table top that was so effective at slowing Fermi’s neutrons. The water in the pond behind the physics building was even better; it had even more hydrogen.

Other elements that have relatively light nuclei are also good moderators, such as carbon (carbon’s nucleus has 6 protons and 6 neutrons). It’s somewhat heavier than you want in order to slow neutrons optimally, but it’s not bad, and its abundant, cheap, and dense. During the Manhattan Project, Fermi (who had fled fascist Italy and settled in the United States) built the first nuclear reactor in a squash court under the football stadium at the University of Chicago. His reactor was a “pile” of uranium balls, with each ball surrounded by blocks of graphite (almost pure carbon, like the lead in this pencil). The uranium was the fuel and the graphite was the moderator.

Before we talk more about moderators, you might be wondering why Fermi’s reactor didn’t explode, destroying Chicago? One reason was that his uranium was a mix of uranium-235 and uranium-238, and was in fact 99% uranium-238. The uranium-238 doesn’t contribute to the chain reaction; it’s not fissile. To make matters worse, uranium-238 can absorb a neutron and dampen the chain reaction. When uranium-238 captures a neutron to become uranium-239, it takes a neutron “out of action” so to speak. During the Manhattan Project the United States spent enormous amounts of time and money separating uranium-235 from uranium-238, so it could use almost pure uranium-235 in the atomic bomb. But Fermi didn’t have any such enriched uranium. Also, Fermi controlled his reactor using a super-duper neutron absorber, cadmium. Cadmium sucks up the neutrons, stopping the chain reaction. Fermi could push in or pull out cadmium control rods to keep the speed of the reaction “just right.” As an emergency backup he had one big cadmium control rod suspended over the reactor by a rope. One of Fermi’s assistants stood by with an axe. If things started to go out of control, his job was to cut the rope dropping the cadmium rod and stopping the reaction. Fortunately, Fermi took great pains to operate the reactor carefully, and no such problems occurred. Had things gone wrong, the reactor probably wouldn’t have exploded like a bomb. It would have just gotten very hot and melted, causing a “meltdown” with all sorts of radiation release, like at Chernobyl. It’s a scary thought because it was in the middle of Chicago, but we were at war against the Nazis, so people took some risks.

Now back to the moderator. Let’s consider three different moderators. First, “heavy water.” This is water containing a rare, heavy isotope of hydrogen, hydrogen-2 (its nucleus consists of one proton and one neutron). While it is not quite as good as hydrogen-1 at slowing down neutrons, it’s still very good, and it has one advantage. Hydrogen-1 (a single proton) can sometimes absorb a neutron to become hydrogen-2. It’s as if occasionally these two balls stick together when they hit. This capture of a neutron slows the chain reaction. Hydrogen-2, however, rarely absorbs a neutron to become hydrogen-3, so it’s a great moderator: it slows the neutrons without absorbing them. During World War Two, the Germans tried to construct a nuclear reactor using heavy water as the moderator. The problem was, heavy water is difficult and expensive to make. There was a plant in Norway that produced heavy water, and it was controlled by the Germans. The British sent in a commando raid that sabotaged the plant, causing all that precious heavy water to flow down the drain. Heavy water is so expensive it isn’t used nowadays in reactors, and we won’t discuss it anymore.

The second moderator we’ll consider is regular water made using hydrogen-1 (I’ll call it just “water” as opposed to “heavy water”). Nowadays most nuclear reactors in the United States use water as the moderator. They also use water as the coolant. You need a coolant to keep the reactor from getting too hot and melting. Also, the coolant is how you get the heat out of the reactor so you can use it to run a steam engine and generate power. So in the United States, water in a nuclear reactor has two purposes: it’s the moderator and the coolant. Suppose that the reactor, for some reason, gets too hot and the water starts boiling off. That will cause the moderator to boil away. No more moderator, no more slowing down the neutrons. No more slowing down the neutrons, no more chain reaction. This is a type of a negative feedback loop that makes the reactor inherently safe. It’s like the thermostat in your house: if the house gets too hot, the thermostat turns off the furnace, and the house cools down. Recall that hydrogen-1 can also absorb neutrons, and in theory that could cause the reactor to speed up when the water boils away because there is less neutron absorption. So neutron absorption and moderation are opposite effects. But a reduction of neutron absorption is less important than the disappearance of the moderator, so on the whole when water boils the reaction slows down. We say that the reactor has a “negative void coefficient.” The “void” means the water is boiling, forming bubbles. The “negative” means this negative feedback loop occurs, keeping the reaction from increasing out of control.

Now for the third moderator: carbon. The Russians built something called an RBMK reactor. This is a Russian acronym, so I won’t try to explain what the different letters mean. Suffice to say, an RBMK reactor is a nuclear reactor that uses carbon as the moderator. Chernobyl was an RBMK reactor. Like Fermi’s original reactor, the carbon was in the form of graphite. In addition, an RBMK reactor uses water as the coolant. Graphite is the moderator and water is the coolant. Now, suppose this type of reactor begins to heat up and the water starts to boil away. The hydrogen in the water is not the primary moderator, the carbon in the graphite is. So, the reaction doesn’t slow down when the water boils away; the carbon moderator is still there, slowing the neutrons. But remember, the hydrogen in water sometimes absorbs a neutron, taking it out of action. This neutron capture decreases as the water boils away, so the reaction increases. Increased heat causes water to boil, causing the reaction to speed up, causing increased heat, causing more water to boil, causing the reaction to speed up even more, causing yet more increased heat … This is a positive feedback loop; a vicious cycle. The reactor has a “positive void coefficient.” It’s as if the thermostat in your house was wired wrong, so when the house got hot the furnace turned ON, heating the house more. Normally the reactor is designed with all sorts of controls to prevent this positive feedback loop from taking off. For instance, control rods can be pushed in and out as needed. But, if for some reason these controls are not in place, the reactor will heat up dramatically and quickly, just as it did at Chernobyl.

Why do we have nuclear reactors? Nuclear reactors produce heat to power a steam engine, which in turn generates electricity. The steam needs to be at high pressure, so it can turn the turbine. Therefore, the reactor is in a pressure container. It’s like a pressure cooker. If the water boils too much, the pressure builds up until the container can’t handle it anymore and bursts, releasing steam. It’s a little like your whistling tea pot, except instead of whistling when the water boils, the reactor explodes. And unlike your teapot, the reactor releases radioactive elements along with the steam. You get a cloud of radioactivity.

Another problem with an RBMK reactor is that graphite burns. It’s pure carbon. It’s like coal. Once the pressure container bursts, oxygen can get in igniting the graphite, starting a fire. The graphite then spews radioactive smoke up into the atmosphere. Many of the people killed in the Chernobyl accident were firemen, trying to put out the fire.

Another issue, a little less important but worth mentioning, was the control rods. Chernobyl had control rods made out of boron, which like cadmium is an excellent neutron absorber. It vacuums up the neutrons and stops the chain reaction. The problem was, the control rods were tipped with graphite. As you push in a control rod, initially it would be like adding moderator, quickening the reaction. Only when the rod was completely pushed in would the boron absorb neutrons, slowing the reaction. So, the control rods would eventually suppress the chain reaction, but initially they made things worse. If, like at Chernobyl, a problem developed quickly, the control rods couldn’t keep up.

I won’t go in to all the comedy of errors that were the immediate cause of the accident at Chernobyl. The reactor was undergoing a test, and several of the controls were turned off. Some safeguards were still in place, but mistakes, poor communication, and ignorance prevented them from working. Whatever the immediate cause of the accident, the crucial point is that the reactor design itself was unstable. It’s like trying to balance this pencil on its tip. You can do it if you are careful and have some controls, but it’s inherently unstable. If you are not always vigilant, the pencil will fall over. The unstable design of the Chernobyl reactor made it a disaster waiting to happen.
If you would like to hear me give this talk (slightly modified), you can watch the YouTube video below. This winter I was teaching the second semester of introductory physics, and when the coronavirus pandemic arrived I had to switch to an online format. I recorded a lecture about Chernobyl when we were discussing nuclear energy.

My Chernobyl talk, given to my Introductory Physics class,
online from home because of Covid-19.