Friday, October 31, 2014

The Biological Risk of Ultraviolet Light From the Sun

In Section 14.9 (Blue and Ultraviolet Radiation) of the 4th edition of Intermediate Physics for Medicine and Biology, Russ Hobbie and I discuss the biological impact of ultraviolet radiation from the sun. Figure 14.28 in IPMB illustrates a remarkable fact about UV light: only a very narrow slice of the spectrum presents a risk for damaging our DNA. Why? For shorter wavelengths, the UV light incident upon the earth’s atmosphere is almost entirely absorbed by ozone and never reaches the earth’s surface. For longer wavelengths, the UV photons do not have enough energy to damage DNA. Only what is known as UVB radiation (wavelengths of 280-315 nm) poses a major risk.

This does not mean that other wavelengths of UV light are always harmless. If the source of UV radiation is, say, a tanning booth rather than the sun, then you are not protected by miles of ozone-containing atmosphere, and the amount of dangerous short-wavelength UV light depends on the details of the light source and the light’s ability to penetrate our body. Also, UVA radiation (315 – 400 nm) is not entirely harmless. UVA can penetrate into the dermis, and it can cause skin cancer by other mechanisms besides direct DNA damage, such as production of free radicals or suppression of the body’s immune system. Nevertheless, Fig. 14.28 shows that UVB light from the sun is particularly effective at harming our genes.

Russ and I obtained Fig. 14.28 from a book chapter by Sasha Madronich (The atmosphere and UV-B radiation at ground level. In: Environmental UV Photobiology, Plenum, New York, 1993). Madronich begins
“Ultraviolet (UV) radiation emanating from the sun travels unaltered until it enters the earth’s atmosphere. Here, absorption and scattering by various gases and particles modify the radiation profoundly, so that by the time it reaches the terrestrial and oceanic biospheres, the wavelengths which are most harmful to organisms have been largely filtered out. Human activities are now changing the composition of the atmosphere, raising serious concerns about how this will affect the wavelength distribution and quantity of ground-level UV radiation.”
Madronich wrote his article in the early 1990s, when scientists were concerned about the development of an “ozone hole” in the atmosphere over Antarctica. Laws limiting the release of chlorofluorocarbons, which catalyze the break down of ozone, have resulted in a slow improvement in the ozone situation. Yet, the risk of skin cancer continues to be quite sensitive to ozone concentration in the atmosphere.

Our exposure to ozone is also sensitive to the angle of the sun overhead. Figure 14.28 suggests that at noon in lower latitudes, when the sun is directly overhead, the ozone exposure is about three times greater than when the sun is at an angle of 45 degrees (here in Michigan, this would be late afternoon in June, or noon in September; we never make it to 45 degrees in December). The moral of this story: Beware of exposure to UV light when frolicking on the Copacabana beach at noon on a sunny day in January!

Friday, October 24, 2014

A Log-Log Plot of the Blackbody Spectrum

Section 14.7 (Thermal Radiation) of the 4th edition of Intermediate Physics for Medicine and Biology contains one of my favorite illustrations: Figure 14.24, which compares the blackbody spectrum as a function of wavelength λ and as a function of frequency ν. One interesting feature of the blackbody spectrum is that its peak (the wavelength or frequency for which the most thermal radiation is emitted) is different depending on if you plot it as a function of wavelength (Wλ(λ,T) in units of W m-3) or frequency (Wν(ν,T) in units of W s m-2). The units make more sense if we express the units of Wλ as W m-2 per m, and the units of Wν as W m-2 per Hz.

A few weeks ago I discussed the book The First Steps in Seeing, in which the blackbody spectrum was plotted using a log-log scale. This got me to thinking, “I wonder how Fig. 14.24 would look if all axes were logarithmic?” The answer is shown below.


The caption for Fig. 14.24 is “The transformation from Wλ(λ,T) to Wν(ν,T) is such that the same amount of power per unit area is emitted in wavelength interval (λ, dλ) and the corresponding frequency interval (ν, dν). The spectrum shown is for a blackbody at 3200 K.” I corrected the wrong temperature T in the caption as printed in the 4th edition.

The bottom right panel of the above figure is a plot of Wλ versus λ. For this temperature the spectrum peaks just a bit below λ = 1 μm. At longer wavelengths, it falls off approximately as λ-4 (shown as the dashed line, known as the Rayleigh-Jeans approximation). At short wavelengths, the spectrum rises abruptly and is exponential.

The top left panel contains a plot of Wν versus ν. The spectrum peaks at a frequency just below about 0.3 THz. At low frequencies it increases approximately as ν2 (again, the Rayleigh-Jeans approximation). At high frequencies the spectrum falls dramatically and exponentially.

The connection between these two plots is illustrated in the upper right panel, which plots the relationship ν = c/λ. This equation has nothing to do with blackbody radiation, but merely shows a general relationship between frequency, wavelength, and the speed of light for electromagnetic radiation.

Why is it useful to show these functions in a log-log plot? First, it reinforces the concepts Russ Hobbie and I introduced in Chapter 2 (Exponential Growth and Decay) of IPMB. In a log-log plot, power laws appear as straight lines. Thus, in the book’s version of Fig. 14.24 the equation ν = c/λ is a hyperbola, but in the log-log version this is a straight line with a slope of negative one. Furthermore, the Rayleigh-Jeans approximation implies a power-law relationship, which is nicely illustrated on a log-log plot by the dashed line. In the book’s version of the figure, Wλ falls off at both large and small wavelengths, and at first glance the rates of fall off look similar. You don’t really see the difference until you look at very small values of Wλ, which are difficult to see in a linear plot but are apparent in a logarithmic plot. The falloff at short wavelengths is very abrupt while the decay at long wavelengths is gradual. This difference is even more striking in the book’s plot of Wν. The curve doesn’t even go all the way to zero frequency in Fig. 14.24, making its limiting behavior difficult to judge. The log-log plot clearly shows that at low frequencies Wν rises as ν2.

Both the book’s version and the log-log version illustrate how the two functions peak at different regions of the electromagnetic spectrum, but for this point the book’s linear plot may be clearer. Another advantage of the linear plot is that I have an easier time estimating the area under the curve, which is important for determining the total power emitted by the blackbody and the Stefan-Boltzmann law. Perhaps there is some clever way to estimate areas under a curve on a log-log plot, but it seems to me the log plot exaggerates the area under the small frequency section of the curve and understates the area under the large frequencies (just as on a map the Mercator projection magnifies the area of Greenland and Antarctica). If you want to understand how these functions behave completely, look at both the linear and log plots.

Yet another way to plot these functions would be on a semilog plot. The advantage of semilog is that an exponential falloff shows up as a straight line. I will leave that plot as an exercise for the reader.

For those who want to learn about the derivation and history of the blackbody spectrum, I recommend Quantum Physics of Atoms, Molecules, Solids, Nuclei, and Particles (although any good modern physics book should discuss this topic). A less mathematical but very intuitive description of why Wλ and Wν peak at different parts of the spectrum is given in The Optics of Life. For a plot of photon number (rather than energy radiated) as a function of λ or ν, see The First Steps in Seeing.

Friday, October 17, 2014

A Theoretical Model of Magneto-Acoustic Imaging of Bioelectric Currents

Twenty years ago, I became interested in magneto-acoustic imaging, primarily influenced by the work of Bruce Towe that was called to my attention by my dissertation advisor and collaborator John Wikswo. (See, for example, Towe and Islam, A magneto-acoustic method for the noninvasive measurement of bioelectric currents, IEEE Trans. Biomed. Eng., 35:892-894, 1988). The result was a paper by Wikswo, Peter Basser, and myself (A theoretical model of magneto-acoustic imaging of bioelectric currents, IEEE Trans. Biomed. Eng., 41:723-728, 1994). This was my first foray into biomechanics, a subject that has become increasingly interesting to me, to the point where now it is the primary focus of my research (but that’s another story; for a preview look here).

I started learning about biomechanics mainly through my friend Peter Basser. We both worked at the National Institutes of Health in the early 1990s. Peter used continuum models in his research a lot, and owned a number of books on the subject. He also loved to travel, and would often use his leftover use-or-lose vacation days at the end of the year to take trips to exotic places like Kathmandu. When he was out of town on these adventures, he left me access to his personal library, and I spent many hours in his office reading classic references like Schlichting’s Boundary Layer Theory and Love’s A Treatise on the Mathematical Theory of Elasticity. Peter and I also would read each other’s papers, and I learned much continuum mechanics from his work. (NIH had a rule that someone had to sign a form saying they read and approved a paper before it could be submitted for publication, so I would give my papers to Peter to read and he would give his to me.) In this way, I became familiar enough with biomechanics to analyze magneto-acoustic imaging. Interestingly, we published our paper in the same year Basser began publishing his research on MRI diffusion tensor imaging, for which he is now famous (see here).

As with much of my research, our paper on magneto-acoustic imaging addressed a simple “toy model": an electric dipole in the center of an elastic, conducting sphere exposed to a uniform magnetic field. We were able to calculate the tissue displacement and pressure analytically for the cases of a magnetic field parallel and perpendicular to the dipole. One of my favorite results in the paper was that we found a close relationship between magneto-acoustic imaging and biomagnetism: “Magneto-acoustic pressure recordings and biomagnetic measurements image action currents in an equivalent way: they both have curl J [the curl of the current density] as their source”.

For about ten years, our paper had little impact. A few people cited it, including Amalric Montalibet and later Han Wen, who each developed a method to use ultrasound and the Lorentz force to generate electrical current in tissue. I’ve described this work before in a review article about the role of magnetic forces in medicine and biology, which I have mentioned before in this blog. Then, in 2005 Bin He began citing our work in a long list of papers about magnetoacoustic tomography with magnetic induction, which again I've written about previously. His work generated so much interest in our paper that in 2010 alone it was cited 19 times according to Google Scholar. Of course, it is gratifying to see your work have an impact.

But the story continues with a more recent study by Pol Grasland-Mongrain of INSERM in France. Building on Montalibet’s work, Grasland-Mongrain uses an ultrasonic pulse and the Lorentz force to induce a voltage that he can detect with electrodes. The resulting electrical data can then be analyzed to determine the distribution of electrical conductivity (see Ammari, Grasland-Mongrain, et al. for one way to do this mathematically). In many ways, their technique is in competition with Bin He’s MAT-MI as a method to image conductivity.

Grasland-Mongrain also publishes his own blog about medical imaging. (Warning: The website is in French, and I have to rely on Google Translate to read it. It is my experience that Google has a hard time translating technical writing). There he discusses his most recent paper about imaging shear waves using the Lorentz force. Interestingly, shear waves in tissue is one of the topics Russ Hobbie and I added to the 5th edition of Intermediate Physics for Medicine and Biology, due out next year. Grasland-Mongrain’s work has been highlighted in Physics World and Focus Physics, and a paper about it appeared this year in Physical Review Letters, the most prestigious of all physics journals (and one I’ve never published in, to my chagrin).

I am amazed by what can happen in twenty years.

As a postscript, let me add a plug for toy models. Russ and I use a lot of toy models in IPMB. Even though such simple models have their limitations, I believe they provide tremendous insight into physical phenomena. I recently reviewed a paper in which the authors had developed a very sophisticated and complex model of a phenomena, but examination of a toy model would have told them that the signal they calculated was far, far to small to be observable. Do the toy model first. Then, once you have the insight, make your model more complex.

Friday, October 10, 2014

John H Hubbell

In the references at the end of Chapter 15 (Interaction of Photons and Charged Particles with Matter) in the 4th edition of Intermediate Physics for Medicine and Biology, you will find a string of publications authored by John H. Hubbell (1925-2007), covering a 27-year period from 1969 until 1996. Data from his publications are plotted in Fig. 15.2 (Total cross section for the interactions of photons with carbon vs. photon energy), Fig. 15.3 (Cross sections for the photoelectric effect and incoherent and coherent scattering from lead), Fig. 15.8 (The coherent and incoherent differential cross sections as a function of angle for 100-keV photons scattering from carbon, calcium, and lead), Fig. 15.14 (Fluorescence yields from K-, L-, and M-shell vacancies as a function of atomic number Z), and Fig. 15.16 (Coherent and incoherent attenuation coefficients and the mass energy absorption coefficient for water).

Hubbell’s 1982 paper “Photon mass attenuation and energy-absorption coefficients from 1 keV to 20 MeV” (International Journal of Applied Radiation and Isotopes, volume 33, pages 1269-1290) has been cited 976 times according to the Web of Science. In fact, it has been cited so many times that it was selected as a citation classic, and Hubbell was invited to write a one-page reminiscence about the paper. It began modestly
“Some papers become highly cited due to the creativity, genius, and vision of the authors, presenting seminal work stimulating and opening up new and multiplicative lines of research. Another, more pedestrian class of papers is ‘house-by-the-side-of-the-road’ works, highly cited simply because these papers provide tools required by a substantial number of researchers in a single discipline or perhaps in several diverse disciplines, as is here the case."
At the time of his death, the International Radiation Physics Society Bulletin published the following obituary
"The International Radiation Physics Society (IRPS) lost one of its major founding members, and the field of radiation physics one of its advocates and contributors of greatest impact, with the death this spring of John Hubbell.

John was born in Michigan in 1925, served in Europe in World War II [he received a bronze star], and graduated from the University of Michigan with a BSE (physics) in 1949 and MS (physics) in 1950. He then joined the National Bureau of Standards (NBS), later NIST, where he worked during his entire career. He married Jean Norford in 1955, and they had three children. He became best known and cited for National Standards Reference Data Series Report 29 (l969), "Photon Cross Sections, Attenuation Coefficients, and Energy Absorption Coefficients from 10 keV to 100 GeV". He was one of the two leading founding members of the International Radiation Physics Society in 1985, and he served as its President 1994-97. While he retired from NIST in 1988, he remained active there and in the affairs of IRPS, until the stroke that led to his death this year.”
Learn more about John Hubbell here and here.

Friday, October 3, 2014

Update on the 5th edition of IPMB

A few weeks ago, Russ Hobbie and I submitted the 5th edition of Intermediate Physics for Medicine and Biology to our publisher. We are not done yet; page proofs should arrive in a few months. The publisher is predicting a March publication date. I suppose whether we meet that target will depend on how fast Russ and I can edit the page proofs, but I am nearly certain that the 5th edition will be available for fall 2015 classes (for summer 2015, I am hopeful but not so sure). In the meantime, use the 4th edition of IPMB.

What is in store for the 5th edition? No new chapters; the table of contents will look similar to the 4th edition. But there are hundreds—no, thousands—of small changes, additions, improvements, and upgrades. We’ve included many new up-to-date references, and lots of new homework problems. Regular readers of this blog may see some familiar additions, which were introduced here first. We tried to cut as well as add material to keep the book the same length. We won’t know for sure until we see the page proofs, but we think we did a good job keeping the size about constant.

We found several errors in the 4th edition when preparing the 5th. This week I updated the errata for the 4th edition, to include these mistakes. You can find the errata at the book website, www.oakland.edu/~roth/hobbie.htm. I won’t list here the many small typos we uncovered, and all the misspellings of names are just too embarrassing to mention in this blog. You can see the errata for those. But let me provide some important corrections that readers will want to know about, especially if using the book for a class this fall or next winter (here in Michigan we call the January-April semester "winter"; those in warmer climates often call it spring).
  • Page 78: In Problem 61, we dropped a key minus sign: “90 mV” should be “-90 mV”. This was correct in the 3rd edition (Hobbie), but somehow the error crept into the 4th (Hobbie and Roth). I can’t figure out what was different between the 3rd and 4th editions that could cause such mistakes to occur.
  • Page 137: The 4th edition claimed that at a synapse the neurotransmitter crosses the synaptic cleft and “enters the next cell”. Generally a neurotransmitter doesn’t “enter” the downstream cell, but is sensed by a receptor in the membrane that triggers some response.
  • Page 338: I have already told you about the mistake in the Bessel function identity in Problem 10 of Chapter 12. For me, this was THE MOST ANNOYING of all the errors we have found. 
  • Page 355: In Problem 12 about sound and hearing, I used an unrealistic value for the threshold of pain, 10-4 W m-2. Table 13.1 had it about right, 1 W m-2. The value varies between people, and sometimes I see it quoted as high as 10 W m-2. I suggest we use 1 W m-2 in the homework problem. Warning: the solution manual (available to instructors who contact Russ or me) is based on the problem as written in the 4th edition, not on what it would be with the corrected value.
  • Page 355: Same page, another screw up. Problem 16 is supposed to show how during ultrasound imaging a coupling medium between the transducer and the tissue can improve transmission. Unfortunately, in the problem I used a value for the acoustic impedance that is about a factor of a thousand lower than is typical for tissue. I should have used Ztissue = 1.5 × 106 Pa s m-1. This should have been obvious from the very low transmission coefficient that results from the impedance mismatch caused by my error. Somehow, the mistake didn’t sink in until recently. Again, the solution manual is based on the problem as written in the 4th edition.
  • Page 433: Problem 30 in Chapter 15 is messed up. It actually contains two problems, one about the Thomson scattering cross section, and another (parts a and b) about the fraction of energy due to the photoelectric effect. Really, the second problem should be Problem 31. But making that change would require renumbering all subsequent problems, which would be a nuisance. I suggest calling the second part of Problem 30 as Problem "30 ½". 
  • Page 523: When discussing a model for the T1 relaxation time in magnetic resonance imaging, we write “At long correlation times T1 is proportional to the Larmor frequency, as can be seen from Eq. 18.34.” Well, a simple inspection of Eq. 18.34 reveals that T1 is proportional to the SQUARE of the Larmor frequency in that limit. This is also obvious from Fig. 18.12, where a change in Larmor frequency of about a factor of three results in a change in T1 of nearly a factor of ten. 
  • Page 535: In Chapter 18 we discuss how the blood flow speed v, the repetition time TR, and the slice thickness Δz give rise to flow effects in MRI. Following Eq. 18.56, we take the limit when v is much greater than "TR/Δz". I always stress to my students that units are their friends. They can spot errors by analyzing if their equation is dimensionally correct. CHECK IF THE UNITS WORK! Clearly, I didn’t take my own advice in this case.
For all these errors, I humbly apologize. Russ and I redoubled our effort to remove mistakes from the 5th edition, and we will redouble again when the page proofs arrive. In the meantime, if you find still more errors in the 4th edition, please let us know. If the mistake is in the 4th edition, it could well carry over to the 5th edition if we don’t root it out immediately.