Friday, August 30, 2013

The Ascent of Sap in Trees

In Chapter 1 of the 4th edition of Intermediate Physics for Medicine and Biology, Russ Hobbie and I included a homework problem about moving water up trees.
Problem 34 Sap flows up a tree at a speed of about 1 mm s−1 through its vascular system (xylem), which consists of cylindrical pores of 20 μm radius. Assume the viscosity of sap is the same as the viscosity of water. What pressure difference between the bottom and top of a 100 m tall tree is needed to generate this flow? How does it compare to the hydrostatic pressure difference caused by gravity?
When you calculate the pressure needed to push water (that is, sap) up the tree through the xylem, you get (Spoiler Alert!) twenty atmospheres to overcome the viscous resistance of the pores, and ten atmospheres to overcome gravity. How does the tree generate all this pressure? That is a famous old problem known as the “ascent of sap.”

Now I admit that a 100-meter tree is, indeed, very tall; taller than even the Statue of Liberty. But it is not an unrealistic example. The majestic sequoias in California reach this height. The tallest known tree, named Hyperion, is a sequoia (coast redwood) in northern California’s Redwood National and State Parks that reaches a height of 115 m. The leaves at the top of that redwood need water to carry out photosynthesis. How do they get it?

First, let us consider some mechanisms that do not work. The tree cannot suck the water up, as if it were a gigantic drinking straw. Even if the tree could produce a true vacuum at its peak it could only create a pressure difference of one atmosphere, which corresponds to a rise of water of 10 m. Another idea is that the water rises by capillary action, like a giant wick. But the height that can be reached by climbing up a tube via surface tension is inversely proportional to the tube radius, and for xylem’s 20 micron radius tubes water will rise only a tiny fraction of the tree’s height (in Sec. 12.2 of his book Air and Water, Mark Denny estimates that water would rise in xylem by capillary action to only a height of three-fourths of a meter). Osmotic pressure won’t work either, for any realistic concentration gradient. So what is the answer?

There is still some controversy, but the generally accepted mechanism for the ascent of sap is called the cohesion-tension theory. In the leaves, capillary action through very tiny channels helps pull water upwards to replace that which evaporates from the leaf surface. In the larger pores of the xylem, the water is pulled by tension (negative pressure), somewhat like a steel cable pulling an elevator up its shaft. But can water support such a tension? It can, but there is a problem. If any air is present, the system will fail. Think of a piston filled half with water and half air. If you pull on the piston, you will just expand the air as its pressure is reduced. Now, consider the piston with only 1/4th air and 3/4th water; the air still expands when you pull. In fact, if there is even one bubble present in the water, pulling on the piston will cause it to expand. Only if the piston contains no air at all will the water be able to exert a tension force. In other words, water under negative pressure is susceptible to cavitation; the formation of bubbles. Fortunately, the structure of xylem is such that bubbles cannot grow indefinitely, but get trapped in one compartment.

For more details, see “The Cohesion-Tension Mechanism and the Acquisition of Water by Plant Roots,” by Ernst Steudle (Annual Review of Plant Physiology, Volume 52, Pages 847–875, 2001). Below I reproduce his summary of cohesion-tension theory. Note that 100 MPa is 1000 atm!
  • Water has high cohesive forces. It can be subjected to from some ten to several hundred MPa before columns break. When subjected to tensions, water is in a metastable state, i.e. pressure in xylem vessels is much smaller than the equilibrium water vapor pressure at the given temperature. 
  •  Walls of vessels represent the weak part of the system. They may contain air or seeds of water vapor. When a critical tension is reached in the lumen of xylem vessels, pits in vessel walls allow the passage of air through them, resulting in cavitation (embolism). 
  • Water in vessels of higher plants forms a continuous system from evaporating surfaces in the leaves to absorbing surfaces of the roots
and into the soil (soil-plant-air-continuum; SPAC). With few exceptions, water flow within the SPAC is hydraulic in nature, and the system can be described as a network of resistors arranged in series and in parallel. 
  • Evaporation from leaves lowers their water potential and causes water to move from the xylem to evaporating cells across leaf tissue. This reduces the pressure in the xylem, often to values well below zero (vacuum). 
  • Gradients in pressure (water potential) are established along transpiring plants; this causes an inflow of water from the soil into the roots and to the transpiring surfaces in the leaves.
Here is an animation that nicely summarizes this process.

I find the idea of water being hoisted up a tree by tens of atmospheres of tension to be fascinating, if a bit disconcerting. This phenomenon offers a fine example of the important role of physics in biology.

Friday, August 23, 2013

Stealth Nanoparticles Boost Radiotherapy

I hope, dear readers, that you all have been regularly browsing through http://medicalphysicsweb.org, the website from the Institute of Physics dedicated to medical physics news. I was particularly taken by the article published there this week titled “Stealth Nanoparticles Boost Radiotherapy.” Russ Hobbie and I don’t talk about nanoparticles in the 4th edition of Intermediate Physics for Medicine and Biology, but they are a hot topic in biomedical research these days. The article by freelance journalist Cynthia Keen begins
Imagine a microscopic bomb precisely positioned inside a cancer tumour cell that explodes when ignited by a dose of precision-targeted radiotherapy. The cancerous tumour is destroyed. The healthy tissue surrounding it survives.

This scenario may become reality within a decade if research by Massachusetts scientists on using nanoparticles to deliver cancer-fighting drugs proceeds smoothly. Wilfred F Ngwa, a medical physicist in the department of radiation oncology at Brigham and Women's Hospital and Dana Farber Cancer Institute in Boston, described the latest initiative at the AAPM annual meeting, held earlier this month in Indianapolis, IN. 
We discuss radiation therapy in Chapter 16 of IPMB. The trick of radiotherapy is to selectively kill cancer cells while sparing normal tissue. The nanoparticles are designed to target tumors
by applying tumour vasculature-targeted cisplatin, Oxaliplatin or carboplatin [three widely used, platinum-based chemotherapy drugs] nanoparticles during external-beam radiotherapy, a substantial photon-induced boost to tumour endothelial cells can be achieved. This would substantially increase damage to the tumour’s blood vessels, as well as cells that cause cancer to recur, while also delivering chemotherapy with fewer toxicities.
In general, nanoparticles typically have a size on the order of 10 to 100 nm. This size passes easily through the smallest blood vessels, but is too big to pass through ion channels in the cell membrane. It is about the size of a large biomolecule or a small virus. Nanoparticles are used in imaging and therapy. For an overview, see the review by Shashi Murthy (International Journal of Nanomedicine, Volume 2, Pages 129–141, 2007).

The medicalphysicsweb article concludes
“The promising result of using approved platinum-based nanoparticles combined with experimental results of the past two years convince us that our new RAID [radiotherapy application with in situ dose-painting] approach to cancer provides a number of possibilities for customizing and significantly improving radiotherapy,” Ngwa said at the press conference. This research is still in its early stages, with laboratory testing of the new approach in mice ongoing. If tests continue to prove successful, and a grant or private funding is available, it will lead to clinical trials in humans. The researchers are hopeful that they will be able to continue their work without any disruption and to move their novel treatment from laboratory to clinical use. 
Another news story about this research can be found here

Friday, August 16, 2013

We Need Theoretical Physics Approaches to Study Living Systems

An editorial titled “We Need Theoretical Physics Approaches to Study Living Systems,” which was published recently in the journal Physical Biology (Volume 10, Article number 040201), has resonated with me. Krastan Blagoev, Kamal Shukla and Herbert Levine discuss the importance of using simply physical models to understand complicated biological problems. The debate about how much detail to include in mathematical models is a constant source of tension between physicists and biologists, and even between physicists and biomedical engineers. I agree with the editorial’s authors: simple models are vitally important. Biologists (and even more so, medical doctors) put great emphasis on the complexity of their systems. But the value of a simple model is that it highlights the fundamental behavior of a system that is often not obvious from experiments. If we build realistic models including all the complexity, they will be just as difficult to understand as are the experiments themselves. Blagoev, Shukla and Levine say much the same (my italics).
In this editorial, we propose that theoretical physics can play an essential role in making sense of living matter. When faced with a highly complex system, a physicist builds simplified models. Quoting Philip W Anderson’s Nobel prize address, “the art of model-building is the exclusion of real but irrelevant parts of the problem and entails hazards for the builder and the reader. The builder may leave out something genuinely relevant and the reader, armed with too sophisticated an experimental probe, may take literally a schematized model. Very often such a simplified model throws more light on the real working of nature... ” In his formulation, the job of a theorist is to get at the crux of the system by ignoring details and yet to find a testable consequence of the resulting simple picture. This is rather different than the predilection of the applied mathematician who wants to include all the known details in the hope of a quantitative simulacrum of reality. These efforts may be practically useful, but do not usually lead to increased understanding.
In my own research, the best example of simple model building is the prediction of adjacent regions of depolarization and hyperpolarization during electrical stimulation of the heart. Nestor Sepulveda, John Wikswo, and I used the “bidomain model,” which accounts for essential properties of cardiac tissue such as the tissue anisotropy and the relative electrical conductivity of the intracellular and extracellular spaces (Biophysical Journal, Volume 55, Pages 987–999, 1989; I have discussed this study in this blog before). Yet, this model was an enormous simplification. We ignored the opening and closing of ion channels, the membrane capacitance, the curvature of the myocardial fibers, the cellular structure of the tissue, the details of the electrode-tissue interface, the three-dimensional volume of the tissue, and much more. Nevertheless, the model made a nonintuitive qualitative prediction that was subsequently confirmed by experiments. I think the reason this research has made an impact (over 200 citations to the paper so far) is that we were able to strip our model of all the unnecessary details except those key ones underlying the qualitative behavior. The gist of this idea can be found in a quote usually attributed to Einstein: Everything should be made as simple as possible, but no simpler. I must admit, sometimes it pays to be lucky when deciding which features of a model to keep and which to throw out. But it is not all luck; model building is a skill that needs to be learned.

The editorial continues (again, my italics)
A leading biologist once remarked to one of us that a calculation of in vivo cytoskeletal dynamics that did not take into account the fact that the particular cell in question had more than ten isoforms of actin could not possibly be correct. We need to counter that any calculation which takes into account all these isoforms is overwhelmingly likely to be vastly under-constrained and ultimately not useful. Adding more details can often bring us further from reality. Of course, the challenge for models is then falsification, i.e., finding robust predictions which can be directly tested experimentally.
How does one learn and practice model building? One place to start—regular readers of this blog will have already guessed my answer—is the 4th edition of Intermediate Physics for Medicine and Biology. This book, and especially the homework problems at the end of each chapter, provide plenty of examples of model building (for simple models applied to the study of the heart, see Chapter 10, Problems 37–40). I think that this aspect of the book sets it apart from many others texts, which cover the biology in more detail.

Krastan Blagoev is the director of the Physics of Living Systems program at the National Science Foundation. According to the NSF website
The program “Physics of Living Systems” (PoLS) in the Physics Division at the National Science Foundation targets theoretical and experimental research exploring the most fundamental physical processes that living systems utilize to perform their functions in dynamic and diverse environments. The focus should be on understanding basic physical principles that underlie biological function. Proposals that use physics only as a tool to study biological questions are of low priority.
Because I might someday apply for a grant from the PoLS program, let me note that Dr. Blagoev is a gentleman and a scholar, who has done much to advance the application of physics to biology. To learn more about Blagoev, see the April 2008 issue of The Biological Physicist, the newsletter for the Division of Biological Physics of the American Physical Society. Shukla is the director for the “Biomolecular Dynamics, Structure and Function” program at NSF, which I am unlikely ever to seek funding from, so I’ll just say he is probably a good guy too. Levine is the Director of the Center for Theoretical Biological Physics at Rice University.

Friday, August 9, 2013

Martha Chase (1927-2003)

Ten years ago yesterday, the American biologist Martha Chase passed away. Chase is famous for her participation in a fundamental genetics experiment. In collaboration with Alfred Hershey, she performed this experiment in 1952 at Cold Spring Harbor Laboratory (see last week's blog entry).  Their results supported the hypothesis that DNA is the biological molecule that carries genetic information. They showed that the DNA, not the protein, of the bacteriophage T2 (a virus that infects bacteria) entered E. coli upon infection.

The Eighth Day of Creation: The Makers of the Revolution in Biology, by Horance Freeland Judson, suuperimposed on Intermediate Physics for Medicine and Biology.
The Eighth Day of Creation:
The Makers of the Revolution in Biology,
by Horace Freeland Judson.
To describe this experiment, I quote from Horace Freeland Judson’s wonderful book The Eighth Day of Creation: The Makers of the Revolution in Biology.
Hershey and Chase decided to see if they could strip off the empty phage ghosts from the bacteria and find out what they were and where their contents had gone. DNA contains no sulphur; phage protein has no phosphorus. Accordingly, they began by growing phage in a bacterial culture with a radioactive isotope as the only phosphorus in the soup [P32], which was taken up in all the phosphate groups as the DNA of the phage progeny was assembled, or, in the parallel experiment, by growing phage whose coat protein was labelled with hot sulphur [S35]. They used the phage to infect fresh bacteria in broths that were not radioactive, and a few minutes after infection tried to separate the bacteria from the emptied phage coats. “We tried various grinding arrangements, with results that weren’t very encouraging,” Hershey wrote later. Then they made a technological breakthrough, in the best Delbruck fashion of homely improvisation. “When Margaret McDonald loaned us her blender the experiment promptly succeeded.”
This ordinary kitchen blender provided just the right shear forces to strip the empty bacteriophage coats off the bacteria. When tested, those bacteria infected by phages containing radioactive phosphorus were themselves radioactive, but those infected by phages containing radioactive sulphur were not. Thus, the DNA and not the protein is the genetic material responsible for infection. This was truly an elegant experiment. They key was the use of radioactive tracers. Russ Hobbie and I discuss nuclear physics and nuclear medicine in Chapter 17 of the 4th edition of Intermediate Physics for Medicine and Biology. We focus on medical applications of radioactive isotopes, but we should remember that these tracers also have played a crucial role in experiments in basic biology.

Hershey and Chase’s experiment, often called the Warring Blender experiment, is a classic studied in introductory biology classes. It was the high point of Chase’s career. She obtained her bachelor’s degree from the College of Wooster and was then hired by Hershey to work in his Cold Spring Harbor laboratory. She stayed at Cold Spring Harbor only three years, but in that time she and Hershey performed their famous experiment. In 1964 she obtained her PhD from the University of Southern California. Unfortunately, things did not go so well from Chase after that. Writer Milly Dawson tells the story.
In the late 1950s in California, she had met and married a fellow scientist, Richard Epstein, but they soon divorced… Chase suffered several other personal setbacks, including a job loss, in the late 1960s, a period that saw the end of her scientific career. Later, she experienced decades of dementia, with long-term but no short-term memory. [Waclow] Szybalski [a colleague at Cold Spring Harbor Laboratory in the 1950s] remembered his friend as “a remarkable but tragic person.”
A good description of the Hershey-Chase experiment can be found here. You can learn more about life of Martha Chase in obituaries here and here.  Szybalski’s reminiscences are recording in a Cold Spring Harbor oral history available here. Dawson’s tribute can be found here. And most importantly, the 1952 Hershey-Chase paper can be found here.

Friday, August 2, 2013

Cold Spring Harbor Laboratory

A photograph of me standing next to the entrance of Cold Spring Harbor Laboratory.
Me standing next to the entrance of
Cold Spring Harbor Laboratory.
Last week my wife, my mother-in-law, and I made a brief trip to Long Island, New York, where we made a quick stop at the Cold Spring Harbor Laboratory. What a lovely setting for a research center. We drove around the grounds, looking at the various labs. It sits right on a bay off the Long Island Sound, and looks more like a resort than a scientific laboratory. James Watson, of DNA fame, was the long-time director of Cold Spring Harbor Lab.

In the last few years, the lab has begun a thrust into “Quantitative Biology.” This area of research has much overlap with the 4th edition of Intermediate Physics for Medicine and Biology. I view this development as evidence that science is going in “our direction,” toward a larger role for physics and math in medicine and biology. The Cold Spring Harbor website describes the new Simons Center for Quantitative Biology.
Cold Spring Harbor Laboratory (CSHL) has recently opened the Simons Center for Quantitative Biology (SCQB). The areas of expertise in the SCQB include applied mathematics, computer science, theoretical physics, and engineering. Members of the SCQB will interact closely with other CSHL researchers and will apply their approaches to research areas including genomic analysis, population genetics, neurobiology, evolutionary biology, and signal and image processing.
We passed by CSHL during a trip that included stops at Sagamore Hill National Historic Site in Oyster Bay (President Theodore Roosevelt’s home), Planting Fields Arboretum, and the Montauk Point Lighthouse.

Friday, July 26, 2013

Bioelectromagnetism: Principles and Applications of Bioelectric and Biomagnetic Fields

Bioelectromagnetism: Principles and Applications of Bioelectric and Biomagnetic Fields, by Jakko Malmivuo and Robert Plonsey, superimposed on Intermediate Physics for Medicine and Biology.
Bioelectromagnetism,
by Malmivuo and Plonsey.
A good textbook about bioelectricity and biomagnetism is Bioelectromagnetism: Principles and Applications of Bioelectric and Biomagnetic Fields by Jaakko Malmivuo and Robert Plonsey (Oxford University Press, 1995). One of the best features of the book is that it is available for free online at www.bem.fi/book/index.htm. The book covers many of the topics Russ Hobbie and I discuss in Chapters 6–9 of the 4th edition of Intermediate Physics for Medicine and Biology: the cable equation, the Hodgkin and Huxley model, patch-clamp recordings, the electrocardiogram, biomagnetism, the bidomain model, and magnetic stimulation. The book’s introduction outlines its eight parts:
Part I discusses the anatomical and physiological basis of bioelectromagnetism. From the anatomical perspective, for example, Part I considers bioelectric phenomena first on a cellular level (i.e., involving nerve and muscle cells) and then on an organ level (involving the nervous system (brain) and the heart).

Part II introduces the concepts of the volume source and volume conductor and the concept of modeling. It also introduces the concept of impressed current source and discusses general theoretical concepts of source-field models and the bidomain volume conductor. These discussions consider only electric concepts.

Part III explores theoretical methods and thus anatomical features are excluded from discussion. For practical (and historical) reasons, this discussion is first presented from an electric perspective in Chapter 11. Chapter 12 then relates most of these theoretical methods to magnetism and especially considers the difference between concepts in electricity and magnetism.

The rest of the book (i.e., Parts IV–IX) explores clinical applications. For this reason, bioelectromagnetism is first classified on an anatomical basis into bioelectric and bio(electro)magnetic constituents to point out the parallelism between them. Part IV describes electric and magnetic measurements of bioelectric sources of the nervous system, and Part V those of the heart.

In Part VI, Chapters 21 and 22 discuss electric and magnetic stimulation of neural and Part VII, Chapters 23 and 24, that of cardiac tissue. These subfields are also referred to as electrobiology and magnetobiology. Part VIII focuses on Subdivision III of bioelectromagnetism—that is, the measurement of the intrinsic electric properties of biological tissue. Chapters 25 and 26 examine the measurement and imaging of tissue impedance, and Chapter 27 the measurement of the electrodermal response.

In Part IX, Chapter 28 introduces the reader to a bioelectric signal that is not generated by excitable tissue: the electro-oculogram (EOG). The electroretinogram (ERG) also is discussed in this connection for anatomical reasons, although the signal is due to an excitable tissue, namely the retina.
Jaakko Malmivuo is a Professor in the School of Electrical Engineering at Aalto University in Helsinki, Finland. He is also the director of the Ragnar Granit Institute.

Robert Plonsey is the Pfizer-Pratt University Professor Emeritus of Biomedical Engineering at Duke University. This year, he received the IEEE Biomedical Engineering Award “for developing quantitative methods to characterize the electromagnetic fields in excitable tissue, leading to a better understanding of the electrophysiology of nerve, muscle, and brain.” Plonsey is cited on 16 pages of Intermediate Physics for Medicine and Biology, the most of any scientist or author.

Friday, July 19, 2013

Reinventing Physics For Life-Science Majors

The July issue of Physics Today contained an article by Dawn Meredith and Joe Redish titled “Reinventing Physics for Life-Science Majors.” Much in the article is relevant to the 4th edition of Intermediate Physics for Medicine and Biology. The main difference between the goals of their article and IPMB is that they discuss the introductory physics course, whereas Russ Hobbie and I wrote an intermediate-level text. Nevertheless, many of the aims remain the same. Meredith and Redish begin
Physics departments have long been providing service courses for premedical students and biology majors. But in the past few decades, the life sciences have grown explosively as new techniques, new instruments, and a growing understanding of biological mechanisms have enabled biologists to better understand the physiochemical processes of life at all scales, from the molecular to the ecological. Quantitative measurements and modeling are emerging as key biological tools. As a result, biologists are demanding more effective and relevant undergraduate service classes in math, chemistry, and physics to help prepare students for the new, more quantitative life sciences.
Their section on what skills should students learn reads like a list of goals for IPMB:
  • Drawing inferences from equations…. 
  • Building simple quantitative models…. 
  • Connecting equations to physical meaning…. 
  • Integrating multiple representations…. 
  • Understanding the implications of scaling and functional dependence…. 
  • Estimating….”
Meredith and Redish realize the importance of developing appropriate homework problems for life-science students, which is something Russ and I have spent an enormous amount of time on when revising IPMB. “We have spent a good deal of time in conversation with our biology colleagues and have created problems of relevance to them that are also doable by students in an introductory biology course.” They then offer a delightful problem about calculating how big a worm can grow (see their Box 4). They also include a photo of a “spherical cow”; you need to see it to understand. And they propose the Gauss gun (see a video here) as a model for exothermic reactions. They conclude
Teaching physics to biology students requires far more than watering down a course for engineers and adding in a few superficial biological applications. What is needed is for physicists to work closely with biologists to learn not only what physics topics and habits of mind are useful to biologists but also how the biologist’s work is fundamentally different from ours and how to bridge that gap. The problem is one of pedagogy, not just biology or physics, and solving it is essential to designing an IPLS [Introductory Physics for the Life Sciences] course that satisfies instructors and students in both disciplines.

Friday, July 12, 2013

The Bohr Model

One hundred years ago this month, Niels Bohr published his model of the atom (“On the Constitution of Atoms and Molecules,” Philosophical Magazine, Volume 26, Pages 1–25, 1913). In the May 2013 issue of Physics Today, Helge Kragh writes
Published in a series of three papers in the summer and fall of 1913, Niels Bohr’s seminal atomic theory revolutionized physicists’ conception of matter; to this day it is presented in high school and undergraduate-level textbooks.
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.
I find Bohr’s model fascinating for several reasons: 1) it was the first application of quantum ideas to atom structure, 2) it predicts the size of the atom, 3) it implies discrete atomic energy levels, 4) it explains the hydrogen spectrum in terms of transitions between energy levels, and 5) it provides an expression for the Rydberg constant in terms of fundamental parameters. In his book The Making of the Atomic Bomb, Richard Rhodes discusses the background leading to Bohr’s discovery.
Johann Balmer, a nineteenth-century Swiss mathematical physicist, identified in 1885 … a formula for calculating the wavelengths of the spectral lines of hydrogen… A Swedish spectroscopist, Johannes Rydberg, went Balmer one better and published in 1890 a general formula valid for a great many different line spectra. The Balmer formula then became a special case of the more general Rydberg equation, which was built around a number called the Rydberg constant [R]. That number, subsequently derived by experiment and one of the most accurately known of all universal constants, takes the precise modern value of 109,677 cm−1.

Bohr would have known these formulae and numbers from undergraduate physics, especially since Christensen [Bohr’s doctorate advisor] was an admirer of Rydberg and had thoroughly studied his work. But spectroscopy was far from Bohr’s field and he presumably had forgotten them. He sought out his old friend and classmate, Hans Hansen, a physicist and student of spectroscopy just returned from Gottingen. Hansen reviewed the regularity of the line spectra with him. Bohr looked up the numbers. “As soon as I saw Balmer’s formula,” he said afterward, “the whole thing was immediately clear to me.”

What was immediately clear was the relationship between his orbiting electrons and the lines of spectral light… The lines of the Balmer series turn out to be exactly the energies of the photons that the hydrogen electron emits when it jumps down from orbit to orbit to its ground state. Then, sensationally, with the simple formula R = 2π2me4/h3 (where m is the mass of the electron, e the electron charge and h Planck’s constant—all fundamental numbers, not arbitrary numbers Bohr made up) Bohr produced Rydberg’s constant, calculating it within 7 percent of its experimentally measured value!...
In chapter 14 of the 4th edition of Intermediate Physics for Medicine and Biology, Russ Hobbie and I discuss the Bohr model, but interestingly we do not attribute the model to Bohr. However, at other locations in the book, we casually refer to Bohr’s model by name: see Problem 33 of Chapter 15 where we mention “Bohr orbits,” and Sections 15.9 and 16.1.1 where we refer to the “Bohr formula.” I guess we assumed that everyone knows what the Bohr model is (a pretty safe assumption for readers of IPMB). In Problem 4 of Chapter 14 (one of the new homework problems in the 4th edition), the reader is asked to derive the expression for the Rydberg constant in terms of fundamental parameters (you don’t get exactly the same answer as in the quote above; presumably Rhodes didn’t use SI units).

Bohr would become one the principal figures in the development of modern quantum mechanics. He also made fundamental contributions to nuclear physics, and contributed to the Manhattan project. He was awarded the Nobel Prize in Physics in 1922 “for his services in the investigation of the structure of atoms and of the radiation emanating from them.” He is Denmark’s most famous scientist, and for years he led the Institute of Theoretical Physics at the University of Copenhagen. A famous play, titled Copenhagen is about his meeting with former collaborator Werner Heisenberg in then-Nazi-controlled Denmark in 1941. Here is a clip.

Bohr and Heisenberg discussing the uncertainty principle, in Copenhagen.

Physicists around the world are celebrating this 100-year anniversary; for instance here, here, here and here.

I end with Bohr’s own words: an excerpt from the introduction of his first 1913 paper (references removed).
In order to explain the results of experiments on scattering of α rays by matter Prof. Rutherford has given a theory of the structure of atoms. According to this theory, the atoms consist of a positively charged nucleus surrounded by a system of electrons kept together by attractive forces from the nucleus; the total negative charge of the electrons is equal to the positive charge of the nucleus. Further, the nucleus is assumed to be the seat of the essential part of the mass of the atom, and to have linear dimensions exceedingly small compared with the linear dimensions of the whole atom. The number of electrons in an atom is deduced to be approximately equal to half the atomic weight. Great interest is to be attributed to this atom-model; for, as Rutherford has shown, the assumption of the existence of nuclei, as those in question, seems to be necessary in order to account for the results of the experiments on large angle scattering of the α rays.

In an attempt to explain some of the properties of matter on the basis of this atom-model we meet however, with difficulties of a serious nature arising from the apparent instability of the system of electrons: difficulties purposely avoided in atom-models previously considered, for instance, in the one proposed by Sir J. J. Thomson. According to the theory of the latter the atom consists of a sphere of uniform positive electrification, inside which the electrons move in circular orbits. The principal difference between the atom-models proposed by Thomson and Rutherford consists in the circumstance the forces acting on the electrons in the atom-model of Thomson allow of certain configurations and motions of the electrons for which the system is in a stable equilibrium; such configurations, however, apparently do not exist for the second atom-model. The nature of the difference in question will perhaps be most clearly seen by noticing that among the quantities characterizing the first atom a quantity appears—the radius of the positive sphere—of dimensions of a length and of the same order of magnitude as the linear extension of the atom, while such a length does not appear among the quantities characterizing the second atom, viz. the charges and masses of the electrons and the positive nucleus; nor can it be determined solely by help of the latter quantities.

The way of considering a problem of this kind has, however, undergone essential alterations in recent years owing to the development of the theory of the energy radiation, and the direct affirmation of the new assumptions introduced in this theory, found by experiments on very different phenomena such as specific heats, photoelectric effect, Röntgen [etc]. The result of the discussion of these questions seems to be a general acknowledgment of the inadequacy of the classical electrodynamics in describing the behaviour of systems of atomic size. Whatever the alteration in the laws of motion of the electrons may be, it seems necessary to introduce in the laws in question a quantity foreign to the classical electrodynamics, i. e. Planck’s constant, or as it often is called the elementary quantum of action. By the introduction of this quantity the question of the stable configuration of the electrons in the atoms is essentially changed as this constant is of such dimensions and magnitude that it, together with the mass and charge of the particles, can determine a length of the order of magnitude required. This paper is an attempt to show that the application of the above ideas to Rutherford’s atom-model affords a basis for a theory of the constitution of atoms. It will further be shown that from this theory we are led to a theory of the constitution of molecules.

In the present first part of the paper the mechanism of the binding of electrons by a positive nucleus is discussed in relation to Planck’s theory. It will be shown that it is possible from the point of view taken to account in a simple way for the law of the line spectrum of hydrogen. Further, reasons are given for a principal hypothesis on which the considerations contained in the following parts are based.

I wish here to express my thanks to Prof. Rutherford his kind and encouraging interest in this work.

Friday, July 5, 2013

The Spark of Life

The Spark of Life: Electricity in the Human Body, by Frances Ashcroft, superimposed on Intermediate Physics for Medicine and Biology.
The Spark of Life:
Electricity in the Human Body,
by Frances Ashcroft.
This week I finished the book The Spark of Life: Electricity in the Human Body, by Frances Ashcroft. In the introduction, Ashcroft explains the goal of her book.
In essence, this book is a detective story about a special kind of protein—the ion channel—that takes us from Ancient Greece to the forefront of scientific research today. It is very much a tale for today, as although the effects of static electricity and lightning on the body have been known for centuries, it is only in the last few decades that ion channels have been discovered, their functions unravelled and their beautiful, delicate, intricate structures seen by scientists for the first time. It is also a personal panegyric for my favourite proteins, which captured me as a young scientist and never let me go; they have been a consuming passion throughout my life. In Walt Whitman’s wonderful words, my aim is to “sing the body electric.”
The book examines much of the history behind topics that Russ Hobbie and I discuss in the 4th edition of Intermediate Physics for Medicine and Biology, such as the work of Hodgkin and Huxley on the squid nerve axon, the electrocardiogram, and modern medical devices such as the pacemaker and cochlear implants. The book is definitely aimed at a general audience; having worked in the field of bioelectricity, I sometimes wish for more depth in the discussion. For instance, anyone wanting to know the history of pacemakers and defibrillators would probably prefer something like Machines in our Hearts, by Kirk Jeffrey. Nevertheless, it was useful to find the entire field of bioelectricity described in one relatively short and easily readable book. With its focus on ion channels, I consider this book as a popularization of Bertil Hille’s text Ion Channels of Excitable Membranes. Her book was also useful to me as a review of various drugs and neurotransmitters, which I don’t know nearly as much about as I should.

Here is a sample of Ashcroft’s writing, in which she tells about Rod MacKinnon’s determination of the structure of a potassium channel. Russ and I discuss MacKinnon’s work in Chapter 9 (Electricity and Magnetism at the Cellular Level) of IPMB.
A slight figure with an elfin face, MacKinnon is one of the most talented scientists I know. He was determined to solve the problem of how channels worked and he appreciated much earlier than others that the only way to do so was to look at the channel structure directly, atom by atom. This was not a project for the faint-hearted, for nobody had ever done it before, no one really knew how to do it and most people did not even believe it could be done in the near future. The technical challenges were almost insurmountable and at that time he was not even a crystallographer. But MacKinnon is not only a brilliant scientist; he is also fearless, highly focused and extraordinarily hard-working (he is famed for working around the clock, snatching just a few hours’ sleep between experiments). Undeterred by the difficulties, he simultaneously switched both his scientific field and his job, resigning his post at Harvard and moving to Rockefeller University because he felt the environment there was better. Some people in the field wondered if he was losing his mind. In retrospect, it was a wise decision. A mere two years later, MacKinnon received a standing ovation—an unprecedented event at a scientific meeting—when he revealed the first structure of a potassium channel. And ion channels went to Stockholm all over again.
Ashcroft is particularly good at telling the human interest stories behind the discoveries described in the book. There were several interesting tales about neurotoxins; not only the well-known tetrodotoxin, but also others such as saxitoxin, aconite, batrachotoxin, and grayanotoxin. The myotonic goats, who because of an ion channel disease fall over stiff whenever startled, are amazing. Von Humboldt’s description of natives using horses to capture electric eels is incredible. The debate in Parliament about August Waller’s demonstration of the electrocardiogram, using his dog Jimmie as the subject, was funny. If, like me, you enjoy such stories, read The Spark of Life.

I will end with Ashcroft’s description of how Hodgkin and Huxley developed their mathematical model of the action potential in the squid giant axon, a topic covered in detail in Chapter 6 of IPMB.
Having measured the amplitude and time course of the sodium and potassium currents, Hodgkin and Huxley needed to show that they were sufficient to generate the nerve impulse. They decided to do so by theoretically calculating the expected time course of the action potential, surmising that if it were possible to mathematically simulate the nerve impulse it was a fair bet that only the currents they had recorded were involved. Huxley had to solve the complex mathematical equations involved using a hand-cranked calculator because the Cambridge University computer was “off the air” for six months. Strange as it now seems, the university had only one computer at that time (indeed it was the first electronic one Cambridge had). It took Huxley about three weeks to compute an action potential: times have moved on—it takes my current computer just a few seconds to run the same simulation. What is perhaps equally remarkable is that we often still use the equations Hodgkin and Huxley formulated to describe the nerve impulse.

Three years after finishing their experiments, in 1952, Hodgkin and Huxley published their studies in a landmark series of five papers that transformed forever our ideas about how nerves work. The long time between completing their experiments and publication seems extraordinary to present-day scientists, who would be terrified of being scooped by their rivals. Not so in the 1950s—Huxley told me, “It never even entered my head.” In 1963, Hodgkin and Huxley were awarded the Nobel Prize. Deservedly so, for they got such beautiful results and analysed them so precisely that they revolutionized the field and provided the foundations for modern neuroscience.
For more about The Spark of Life, see here, here, here, and here. Listen to and watch Ashcroft being interviewed by Denis Noble here, and giving the Croonian Lecture here.

Francis Ashcroft being interviewed by Denis Noble.

Francis Ashcroft giving the Croonian Lecture.

Friday, June 28, 2013

Lotka-Volterra equations

Russ Hobbie and I don’t study population dynamics much in the 4th edition of Intermediate Physics for Medicine and Biology. To me, it’s more of a mathematical biology topic rather than an application of physics to biology. However, we do discuss one well-known model for population dynamics, the Lotka-Volterra equations, in a Homework Problem in Chapter 2.
Problem 34 Consider a classic predator-prey problem. Let the number rabbits be R and the number of foxes be F. The rabbits eat grass, which is plentiful. The foxes eat only rabbits. The number of rabbits and foxes can be modeled by the Lotka-Volterra equations
dR/dt = a R – b R F
dF/dt = - c F + d R F .
(a) Describe the physical meaning of each term on the right-hand side of each equation. What does each of the constants a, b, c, and d denote?
(b) Solve for the steady-state values of R and F.

These differential equations are difficult to solve because they are nonlinear (see Chapter 10). Typically, R and F oscillate about the steady-state solutions found in part (b). For more information, see Murray (2001).
There are two steady-state solutions. One is the trivial R = F = 0. The most interesting aspect of this solution is that it is not stable. If R and F are both small, the nonlinear terms in the Lotka-Volterra equations are negligible, and the number of foxes falls exponentially (they is no prey to eat) but the number of rabbits rises exponentially (there is no predator to eat them).

The other steady-state solution is (spoiler alert!) R = c/d and F = a/b. We claim in the problem that these equations are difficult to solve, and that is true in general, at least when searching for analytical solutions. However, if we focus on small deviations from this steady-state, we can solve the equations. Let  

R = c/d + r 
F = a/b + f ,

where r and f are small (much less than the steady state solutions). Plug these into the original differential equations, and ignore any terms containing r times f (these “doubly small” terms are negligible). The new equations for r and f are

dr/dt = - b (c/d) f 
df/dt =   d (a/b) r .

Now let’s use my favorite technique for solving differential equations: guess and check. I will guess

r = A sin(ωt)
f = B cos(ωt) .

If we plug these expressions into the differential equations, we get a solution only if ω2 = ac. In that case, B = -(d/b) √(a/c) A. You can’t get A in this way; it depends on the initial conditions.

A plot of the solution shows two oscillating populations, with the rabbits lagging 90 degrees behind the foxes. In words, suppose you start with foxes at their equilibrium value, but a surplus of rabbits above their equilibrium. In this case, there are lots of rabbits for the foxes to eat, so the foxes gorge themselves and their population grows. However, as the number of foxes rises, the number of rabbits starts to fall (they are being ravaged by all those foxes). After a while, the number of rabbits declines back to its equilibrium value, but by then the number of foxes has surged above its steady-state value. Foxes continue to devour rabbits, reducing the rabbit population below equilibrium. Now there are too many foxes competing for too few rabbits, so the fox population starts to shrink as some inevitably go hungry. During this difficult time, both populations are plummeting as a large but decreasing number of ravenous foxes hunt the rare and frightened rabbits. When the foxes finally fall back to their equilibrium value there is a shortage of rabbits, so the foxes continue to starve and their number keeps falling. With less foxes, the rabbits breed like…um…rabbits and begin to make a comeback. Once they climb to their equilibrium value, there are still relatively few foxes, so the rabbits prosper all the more. With the rabbit population surging, there is plenty of food for the foxes, and the fox population begins to increase. During these happy days, both populations thrive. Eventually, the foxes return to their equilibrium value, but by this time the rabbits are plentiful. But this is just where we started, so the process repeats, over and over again. I needed a lot of words to explain about those foxes and rabbits. I think you can begin to see the virtue of a succinct mathematical analysis, rather than a verbose nonmathematical description.

For larger oscillations, the nonlinear nature of the model becomes important. The populations still oscillate, but not sinusoidally. For some parameters, one population may rise slowly and then suddenly drop precipitously, only to gradually rise again. You can see some of those results here and here.

The Lotka-Volterra model is rather elementary. For instance, there is no damping; the oscillations never decay away but instead continue forever. Moreover, the oscillations do not approach some fixed amplitude (a limit cycle behavior). Instead, the amplitude depends entirely on the initial conditions. Many more realistic models have a threshold, above which oscillations occur but below which the systems returns to its steady state.

Mathematical Biology, by James Murray, superimposed on Intermediate Physics for Medicine and Biology.
Mathematical Biology,
by James Murray.
Population dynamics is a large field, of which we only scratch the surface. One place to learn more is James Murray’s book that we cite at the end of the homework problem:

Murray, J. D. (2001) Mathematical Biology. New York, Springer-Verlag.
The most recent (3rd) edition of Murray’s book is actually in two volumes:
Murray, J. D. (2002) Mathematical Biology: I. An Introduction. New York, Springer-Verlag. 

Murray, J. D. (2002) Mathematical Biology: II. Spatial Models and Biomedical Applications. New York, Springer-Verlag.
Hear Murray talk about his research here


Listen to James Murray talk about Mathematical Biology.
https://www.youtube.com/watch?v=6Yj5Nyb_VyU

Alfred Lotka (1880–1949) was an American scientist. In 1925 he published a book, Elements of Physical Biology, that is in some ways a precursor to Intermediate Physics for Medicine and Biology, or perhaps an early version of Murray’s Mathematical Biology. You can download a copy of the book here.