Mar 012012
 

Maxima is an open-source computer algebra system (CAS) and a damn good one at that if I may say so myself, being one of Maxima’s developers.

Among other things, Maxima has top-notch tensor algebra capabilities, which can be used, among other things, to work with Lagrangian field theories.

This week, I am pleased to report, SourgeForge chose Maxima as one of the featured open-source projects on their front page. No, it won’t make us rich and famous (not even rich or famous) but it is nice to be recognized.

 Posted by at 9:35 am
Feb 272012
 

The cover story in a recent issue of New Scientist was titled Seven equations that rule your world, written by Ian Stewart.

I like Ian Stewart; I have several of his books on my bookshelf, including a 1978 Hungarian edition of his textbook, Catastrophe Theory and its Applications.

However, I disagree with his choice of equations. Stewart picked the four Maxwell equations, Schrödinger’s equation, the Fourier transform, and the wave equation:

\begin{align}
\nabla\cdot E&=0,\\
\nabla\times E&=-\frac{1}{c}\frac{\partial H}{\partial t},\\
\nabla\cdot H&=0,\\
\nabla\times H&=\frac{1}{c}\frac{\partial E}{\partial t},\\
i\hbar\frac{\partial}{\partial t}\psi&=\hat{H}\psi,\\
\hat{f}(\xi)&=\int\limits_{-\infty}^{\infty}f(x)e^{-2\pi ix\xi}dx,\\
\frac{\partial^2u}{\partial t^2}&=c^2\frac{\partial^2u}{\partial x^2}.
\end{align}

But these equations really aren’t that fundamental… and some rather fundamental equations are missing.

For starters, the four Maxwell equations really should just be two equations: given a smooth (or at least three times differentiable) vector field \(A\) in 4-dimensional spacetime, we define the electromagnetic field tensor \(F\) and current \(J\) as

\begin{align}
F&={\rm d}A,\\
J&=\star{\rm d}{\star{F}},
\end{align}

where the symbol \(\rm d\) denotes the exterior derivative and \(\star\) represents the Hodge dual. OK, these are not really trivial concepts from high school physics, but the main point is, we end up with a set of four Maxwell equations only because we (unnecessarily) split the equations into a three-dimensional and a one-dimensional part. Doing so also obscures some fundamental truths: notably that once the electromagnetic field is defined this way, its properties are inevitable mathematical identities, not equations imposed on the theoretician’s whim.

Moreover, the wave equation really is just a solution of the Maxwell equations, and conveys no new information. It is not something you invent, but something you derive.

I really have no nit to pick with Schrödinger’s equation, but before moving on to quantum physics, I would have written down the Euler-Lagrange equation first. For a generic theory with positions \(q\) and time \(t\), this could be written as

$$\frac{\partial{\cal L}}{\partial q}-\frac{d}{dt}\frac{\partial{\cal L}}{\partial\dot{q}}=0,$$

where \({\cal L}\) is the Lagrangian, or Lagrange function (of \(q\) and \(\dot{q}\), and possibly \(t\)) that describes this particular physical system. The significance of this equation is that it can be derived from the principle of least action, and tells us everything about the evolution of a system. Once you know the generic positions \(q\) and their time derivatives (i.e., velocities) \(\dot{q}\) at some time \(t=t_0\), you can calculate them at any other time \(t\). This is why physics can be used to make predictions: for instance, if you know the initial position and velocity of a cannonball, you can predict its trajectory. The beauty of the Euler-Lagrange equation is that it works equally well for particles and for fields and can be readily generalized to relativistic theories; moreover, the principle of least action is an absolutely universal one, unifying, in a sense, classical mechanics, electromagnetism, nuclear physics, and even gravity. All these theories can be described by simply stating the corresponding Lagrangian. Even more astonishingly, the basic mathematical properties of the Lagrangian can be used to deduce fundamental physical laws: for instance, a Lagrangian that remains invariant under time translation leads to the law of energy conservation.

The Euler-Lagrange equation remains valid in quantum physics, too. The big difference is that the quantities \(q\) are no longer simple numbers; they are non-commuting quantities, so-called “q-numbers”. These q-numbers sometimes coincide with ordinary numbers but more often, they do not. Most importantly, if \(q\) happens to be an ordinary number, \(\dot{q}\) cannot be, and vice versa. So the initial position and momentum of a quantum system cannot both be represented by numbers at the same time. Exact predictions are no longer possible.

We can still make approximate predictions though, by replacing the exact form of the Euler-Lagrange equation with a probabilistic prediction:

$$\xi(A\rightarrow B)=k\sum\limits_A^B\exp\left(\frac{i}{\hbar}\int_A^B{\cal L}\right),$$

where \(\xi(A\rightarrow B)\) is a complex number called the probability amplitude, the squared modulus of which tells us the likelihood of the system changing from state \(A\) to state \(B\) and the summation is meant to take place over “all possible paths” from \(A\) to \(B\). Schrödinger’s equation can be derived from this, as indeed most of quantum mechanics. So this, then, would be my fourth equation.

Would I include the Fourier transform? Probably not. It offers a different way of looking at the same problem, but no new information content. Whether I investigate a signal in the time domain or the frequency domain, it is still the same signal; arguably, it is simply a matter of convenience as to which representation I choose.

However, Stewart left out at least one extremely important equation:

$$dU=TdS-pdV.$$

This is the fundamental equation of thermodynamics, connecting quantities such as the internal energy \(U\), the temperature \(T\), the entropy \(S\), and the medium’s equation of state (here represented by the pressure \(p\) and volume \(V\).) Whether one derives it from the first principles of axiomatic thermodynamics or from the postulates of statistical physics, the end result is the same: this is the equation that defines the arrow of time, for instance, as all the other fundamental equations of physics work the same even if the arrow of time is reversed.

Well, that’s five equations. What else would I include in my list? The choices, I think, are obvious. First, the definition of the Lagrangian for gravity:

$${\cal L}_\mathrm{grav}=R+2\Lambda,$$

where \(R\) is the Ricci curvature scalar that characterizes the geometry of spacetime and \(\Lambda\) is the cosmological constant.

Finally, the last equation would be, for the time being, the “standard model” Lagrangian that describes all forms of matter and energy other than gravity:

$${\cal L}_\mathrm{SM}=…$$

Its actual form is too unwieldy to reproduce here (as it combines the electromagnetic, weak, and strong nuclear fields, all the known quarks and leptons, and their interactions) and in all likelihood, it’s not the final version anyway: the existence of the Higgs-boson is still an open question, and without the Higgs, the standard model would need to be modified.

The Holy Grail of fundamental physics, of course, is unification of these final two equations into a single, consistent framework, a true “theory of everything”.

 Posted by at 1:18 pm
Feb 222012
 

Why exactly do we believe that stars and more importantly, gas in the outer regions of spiral galaxies move in circular orbits? This assumption lies at the heart of the infamous galaxy rotation curve problem, as the circular orbital velocity for a spiral galaxy (whose visible mass is concentrated in the central bulge) should be proportional to the inverse square root of the distance from the center; instead, observed rotation curves are “flat”, meaning that the velocity remains approximately the same at various distances from the center.

So why do we assume that stars and gas move in circular orbits? Well, it turns out that one key bit of evidence is in a 32-year old paper that was published by two Indian physicists: Radhakrishnan and Sarma (A&A 85, 1980) made observations of hydrogen gas in the direction of the center of the Milky Way, and found that the bulk of gas between the solar system and the central bulge has no appreciable radial velocity.

However, more recent observations may be contradicting this result. Just two years ago, the Radial Velocity Experiment (RAVE) survey (Siebert et al, MNRAS 412, 2010) found, using a sample of several hundred thousand relatively nearby stars, that a significant radial velocity exists, putting into question the simple model that assumes that circular orbits dominate.

 Posted by at 10:03 pm
Feb 222012
 

So maybe neutrinos don’t travel faster than light after all.

Instead, if rumors are to be believed, it was a simple instrumentation problem. There is no official confirmation yet, but according to a statement that also appears on Nature’s news blog, the OPERA team is indeed investigating two problems related to a timer oscillator and an optical fiber connection.

A while back, I wrote that I could identify four possible broad categories for conventional explanations of the OPERA result:

  1. Incorrectly synchronized clocks;
  2. Incorrectly measured distance;
  3. Unaccounted-for delays in the apparatus;
  4. Statistical uncertainties.

Of these, #4 was already out, as the OPERA team verified their result using short duration proton bunches that avoided the use of potentially controversial statistical methods. I never considered #2 a serious possibility, as highly accurate geographic localization is a well established art. Having read and re-read the OPERA team’s description of how they synchronized clocks, I was prepared to discount #1 as well, but then again, incorrect synchronization can arise as a result of equipment failure, so would that fall under #1 or #3?

In any case, it looks like #3, with a dash of #1 perhaps. Once again, conventional physics prevails.

That is, if we can believe these latest rumors.

 Posted by at 8:08 pm
Feb 212012
 

Some thirty thousand years ago, homo sapiens was busy perfecting techniques to produce primitive stone tools. They may have already invented nets, the bow and arrow, and perhaps even ceramics, but they were still a long way away from inventing civilization.

Around the same time, an arctic squirrel in north-eastern Siberia took the fruit of a narrow-leafed campion, a small arctic flower, and hid it in its burrow, never to be touched again. The fruit froze and remained frozen for over three hundred centuries.

It is frozen no longer; rather, it is blooming, thanks to the efforts of a research team led by Svetlana Yashina and David Gilichinsky of the Russian Academy of Sciences. Against all odds, the genetic material in the seed appears to have survived. I say “appears” because such an extraordinary claim will be subject to extraordinary scrutiny, but what I have been reading suggests that this is indeed real: the age of the fruit is confirmed by radioactive dating.

 Posted by at 9:21 am
Feb 162012
 

I always find these numbers astonishing.

The solar constant, the amount of energy received by a 1 square meter surface at 1 astronomical unit (AU) from the Sun is roughly s = 1.37 kW/m2. Given that 1 AU is approximately 150 million kilometers, or r = 1.5 × 1011 m, the surface area of a 1 AU sphere surrounding the Sun would be A = 4πr2 = 2.8 × 1023 m2. Multiplied by the solar constant, we get P = sA = 3.9 × 1026 W, or the energy E = sA = 3.9 × 1026 J every second. Using Einstein’s infamous mass-energy formula E = mc2, where c = 3 × 108 m/s, we can easily calculate how much mass is converted into energy: m = E/c2 = 4.3 × 109 kg. Close to four and a half million tons.

The dominant fusion process in the Sun is the proton-proton chain reaction, in which approximately 0.7% of the total mass of hydrogen is converted into energy. Thus 4.3 million tons of pure energy is equivalent to over 600 millon tons of hydrogen fuel burned every second. (For comparison, the largest ever nuclear device, the Soviet Tsar Bomba, burned no more than a few hundred kilograms of hydrogen to produce a 50 megaton explosion.)

Fortunately, there is plenty where that came from. The total mass of the Sun is 2 × 1030 kg, so if the Sun was made entirely of hydrogen, it could burn for 100 billion years before running out of fuel. Now the Sun is not made entirely of hydrogen, and the fusion reaction slows down and eventually stops long before all the hydrogen is consumed, but we still have a few billion years of useful life left in our middle-aged star. A much bigger (pun intended) problem is that as our Sun ages, it will grow in size; in a mere billion years, the Earth may well become uninhabitable as a result, with the oceans boiling away. I wonder if it’s too early to start worrying about it just yet.

 Posted by at 12:24 pm
Feb 162012
 

Other countries have launched satellites to observe the Earth; observe the Sun; observe the stars; perform physical, chemical, or biological experiments in space; or even for military purposes. But here is a first: trust a Swiss team to propose a microsatellite specifically designed to capture orbital junk and drag it back to the atmosphere to burn it up.

 Posted by at 11:28 am
Jan 272012
 

Normally, I would get tremendously excited to hear about a serious proposal to establish a permanent lunar colony. (Where do I sign up?)

Unfortunately, when Newt Gingrich floated this idea while campaigning in Florida, I did not feel excited at all. That is because I have very little doubt that this was simply an exercise in transparent political opportunism. Mr. Gingrich is hoping to gain some votes in the Space Coast, but I suspect that even residents there, whose livelihood for a long time has depended on a healthy space program, will see through his blatant pandering.

 Posted by at 1:52 pm
Jan 262012
 

NASA’s week of mourning begins tomorrow. The three deadly accidents in NASA’s history all happened in late January/early February. Apollo 1 caught fire 45 years ago on January 27, 1967, killing Grissom, White and Chaffee. Challenger exploded 26 years ago, on January 28, 1986, killing all seven on board. And Columbia broke up during reentry on February 1, 2003, just nine years ago, killing another seven people. Why these accidents all happened during the same calendar week remains a mystery.

 Posted by at 11:49 am
Jan 242012
 

When I write about things like precision orbit determination, I often have to discuss the difference between ephemeris time (ET) and coordinated universal time (UTC). ET is a “clean” time scale: it is essentially the time coordinate of an inertial coordinate frame that is attached to the barycenter of the solar system. On the other hand, UTC is “messy”: it is the time kept by noninertial clocks sitting here on the surface of the Earth. But the fact that terrestrial clocks sit inside the Earth’s gravity well and are subject to acceleration is only part of the picture. There are also those blasted leap seconds. It is because of leap seconds that terrestrial atomic time (TAI) and UTC differ.

Leap seconds arise because we insist on using an inherently wobbly planet as our time standard. The Earth wobbles, sometimes unpredictably (for instance, after a major earthquake) and we mess with our clocks. Quite pointlessly, as a matter of fact. And now, we missed another chance to get rid of this abomination: the International Telecommunication Union failed to achieve consensus, and any decision is postponed until 2015.

For the curious, an approximate formula to convert between TAI and ET is given by ET – TAI = 32.184 + 1.657×10–3 sin E, where E = M + 0.01671 sin M, M = 6.239996 + 1.99096871×10–7 t and t is the time in seconds since J2000 (that is, noon, January 1, 2000, TAI). To convert TAI to UTC, additional leap seconds must be added: 10 seconds for all dates prior to 1972, and then additional leap seconds depending on the date. Most inelegant.

Speaking of leap this and that, I think it’s also high time to get rid of daylight savings time. Its benefits are dubious at best, and I find the practice unnecessarily disruptive.

 Posted by at 12:23 pm
Jan 222012
 

A couple of weeks ago, somewhere I saw a blog comment that mentioned a book, Rad Decision, written by nuclear engineer James Aach.

Back in the late 1970s, when I got my hands on The Prometheus Crisis by Scortia and Robinson, I just couldn’t put the damn thing down; I read through the night and I finished the book by the morning. So naturally, I couldn’t resist the temptation to buy the last in-stock copy of Aach’s book on Amazon.ca.

And I am glad I did. My concerns that it would be a trashy, amateurishly written novel quickly dissipated. Indeed, in a sense it is a lot better than The Prometheus Crisis: the crisis in Aach’s book is far less dramatic, but the story is believable, the characters perhaps more credible.

My only concern: while this book teaches a lot about nuclear power (and why we should not fear it), its likely audience already knows. Those who would benefit the most from reading it, well, won’t.

 Posted by at 7:39 pm
Jan 152012
 

Microsoft’s Windows 7 weather widget tells me that the temperature is -30 degrees Centigrade this morning in Ottawa. I know that this particular reading is an outlier (I don’t know where MSN get their reading from, but often it’s several degrees below that of others) but it’s still darn cold outside… even on our balcony it was -23 this morning. Welcome to Canada in January, I guess…

 

 Posted by at 9:52 am
Jan 082012
 

Neutrinos recently observed by CERN’s OPERA experiment may have been traveling faster than light. Or may have not. I have been discussing with physicists a number of possibilities: the role of statistics, errors in time or distance measurements, comparisons to SN 1987A, Cherenkov radiation, or the necessity for a Lorentz-violating theoretical framework.

Fortunately, there is one thing I did not need to discuss: How faster-than-light neutrinos relate to the Koran. Physics educators in Pakistan, such as Pervez Hoodbhoy writing for the Express Tribune, are not this lucky: they regularly face criticisms from fundamentalists, and if they choose to confront these head-on, they provoke ominous reader comments that call on all Muslims to “reject this evil experiment”.

Yet, there is a glimpse of hope: a Pakistani reader mentions Carl Sagan’s The Demon-Haunted World, one of Sagan’s last books, and a superb one about rational thinking versus superstition. I don’t know how popular Sagan’s book is in Pakistan, but I am glad it’s not forgotten.

 Posted by at 5:29 pm
Jan 082012
 

Has NASA nothing better to do than harass aging astronauts such as Jim Lovell who, some forty years after having survived a near-fatal accident in deep space (caused by NASA’s negligent storage and handling of an oxygen tank), is auctioning off a checklist containing his handwritten notes? A checklist that, had it remained in NASA’s possession, would likely have ended up in a dumpster decades ago?

This is so not kosher. Let Lovell sell his memorabilia in peace. If anyone has a right to do it, the survivors of Apollo 13 certainly do.

 Posted by at 1:12 pm
Dec 242011
 

In 1968, the crew of Apollo 8, for the first time in the history of humanity, disappeared behind another celestial body. When they re-emerged on the other side and saw the Earth rise over the lunar landscape, on much of the Earth it was Christmas Day.

And this is when they sent us Earthlings a Christmas message, which ended with the words, “And from the crew of Apollo 8, we close with good night, good luck, a merry Christmas, and God bless all of you, all of you on the good Earth.”

You don’t need to be religious to find this moment awe-inspiring.

 Posted by at 9:14 am
Dec 142011
 

So I am reading details about the on-going search for the Higgs boson at the LHC. The media hailed the announcements this week as evidence that the hunt is nearing its goal… however, this is by no means conclusive, and instinctively, I’d be inclined to come to the opposite conclusion.

The Higgs boson, if exists as predicted, can decay into many things. It can decay into two photons. Just such a decay, consistent with a Higgs particle that is about 130 times heavier than a proton, was in fact observed by two of the LHC’s detectors, CMS:

and Atlas:

So far so good, but these signals are weak, far from conclusive. Never mind, both CMS and Atlas observed another slight peak. A Higgs particle can, in principle, also decay into two Z-bosons. Indeed, such a decay may be indicated by CMS (that ever so small bump near the extreme left of the plot):

and again, Atlas:

And on top of that, there is yet another decay mode, the Higgs particle decaying into a pair of W-bosons, but it is very difficult to see if anything exists at the extreme left of this plot:

So why does this leave me skeptical? Simple. First, we know that the ZZ and WW decay modes are far more likely than the diphoton (γγ) decay.

So naively, I would expect that if the signal is strong enough to produce noticeable bumps in the diphoton plot, very strong peaks should have been observed already in the ZZ and WW graphs. Instead, we see signals there that are even weaker than the bumps in the diphoton plots. While this is by no means rock solid proof that the Higgs does not exist, it makes me feel suspicious. Second… well, suppose that the Higgs does not exist. We always knew that it is the low energy region, namely the region that is still under consideration (the possibility of a Higgs that is heavier than 130 GeV is essentially excluded) where the Higgs search is the most difficult. So if no Higgs exist, this is precisely how we would expect the search to unfold: narrowing down the search window towards lower energies, just as the data becomes noisier and more and more bumps appear that could be misread as a Higgs that’s just not there.

Then again, I could just be whistling in the dark. We won’t know until we know… and that “until” is at least another year’s worth of data that is to be collected at the LHC. Patience, I guess, is a virtue.

 Posted by at 9:02 pm
Nov 182011
 

The latest OPERA results are in and they are very interesting. They used extremely tight bunches of protons this time, with a pulse width of only a few nanoseconds:

These bunches allowed the team to correlate individual neutrino events with the bunches that originated them. This is what they saw:

Numerically, the result is 62.1 ± 3.7 ns, consistent with their previously claimed result.

In my view, there are four possible categories of things that could have gone wrong with the OPERA experiment:

  1. Incorrectly synchronized clocks;
  2. Incorrectly measured distance;
  3. Unaccounted-for delays in the apparatus;
  4. Statistical uncertainties.

Because this new result does not rely on the statistical averaging of a large number of events, item 4 is basically out. One down, three to go.

 Posted by at 8:45 pm