May 312009
 

I’ve been learning a lot about Web development these days: Dojo and Ajax, in particular. It’s incredible what you can do in Javascript nowadays, sophisticated desktop applications running inside a Web browser. I am spending a lot of time building a complex prototype application that has many features associated with desktop programs, including graphics, pop-up dialogs, menus, and more.

I’ve also been learning a lot about the intricacies Brans-Dicke gravity and about the parameterized post-Newtonian (PPN) formalism. Brans-Dicke theory is perhaps the simplest modified gravity theory that there is, and I have to explain to someone why the gravity theory that I spend time working on doesn’t quite behave like Brans-Dicke theory. In the process, I find out things about Brans-Dicke theory that I never knew.

And, I’ve also been doing a fair bit of SCPI programming this month. SCPI is a standardized way for computers to talk to measurement instrumentation, and an old program I wrote used to use a non-standard way… not anymore.

Meanwhile, in all the spare time that I’ve left, I’ve been learning Brook+, a supercomputer programming language based on C… that is because my new test machine is a supercomputer, sort of, with its graphics card that doubles as a numeric vector processor capable in theory of up to a trillion single precision floating point instructions per second… and nearly as many in practice, in the test programs that I threw at it.

I’m also learning a little more about the infamous cosmological constant problem (why is the cosmological constant at least over 50 orders magnitude too small but not exactly zero?) and about quantum gravity.

As I said in the subject… busy days. Much more fun though than following the news. Still, I did catch in the news that Susan Boyle lost in Britains Got Talent… only because an amazing dance group won:

 Posted by at 3:07 am
May 142009
 

Given the less than perfect record of the Ariane 5 launch vehicle, there was reason for concern given that two new great observatories, Herschell and Planck, were launched on the same rocket this morning. Fortunately, the launch was successful, and both spacecraft are now on their merry way. Herschell is an infrared/submillimeter wavelength telescope, while Planck is “WMAP on steroids”, expected to provide much higher resolution views of the cosmic microwave background than its predecessor.

 Posted by at 4:01 pm
May 112009
 

The space shuttle Atlantis is on its way to the Hubble Space Telescope. If the planned repairs are successful, Hubble may get another five years or more before it has to be decommissioned. Originally, the plan was to return it to the Earth, but in the wake of Columbia, that has been deemed to risky… now, it will be deorbited, to ensure that its large mirror (which is not expected to burn up in the atmosphere) doesn’t fall on an inhabited area. But hopefully, that is still many years away.

I saw the package that is now in the belly of Atlantis when I was at Goddard last summer. Well, maybe not quite the same package as, after the repair mission was postponed, I believe they added a few bits to it, but still, it’s largely the same.

Sadly, this mission may also be the very last really useful mission of the shuttle fleet. That is not to say that they won’t fly several more missions… but they are all to the International Space Station, and while I am an enthusiastic supporter of manned space exploration, flying meaningless circles in low Earth orbit just isn’t it… it’s a waste of money, and a pointless risk to the astronauts’ lives.

 Posted by at 6:30 pm
Apr 272009
 

It’s April 27, in Ottawa, supposedly the second (third? sixth?) coldest capital city in the world. The temperature outside is presently 30.5°C outside (30°C according to the Weather Network) and still rising. Weren’t we wondering this time last year (okay, maybe a little earlier, but just a little) whether or not we were going to break the all-time snowfall record?

 Posted by at 8:15 pm
Apr 272009
 

I’ve run the first realistic tests of the kind of computation that I am planning to perform on my new machine with the GPU “supercomputer” card. Here is a “before” picture:

Self-gravitating star cluster on the CPU

Self-gravitating star cluster on the CPU

And now, the exact same program running on the GPU:

Self-gravitating star cluster on the GPU

Self-gravitating star cluster on the GPU

I’d say that’s quite an improvement. To say the least.

The calculation in this case computed the self-gravitational forces in a cluster of 10,000 stars… it seems that the GPU can perform this computation at least 20 times a second. That’s quite remarkable.

 Posted by at 6:13 pm
Apr 252009
 

Watching the outrage over the DHS memos that purportedly target all Americans on the political right as potential enemies of the state, I have come to the realization that a great many political conspiracy theories are based on a trivial error in formal logic: namely, that the implication operator is not commutative.

The implication operator, AB (A implies B) is true if A is false (B can be anything) or if both A and B are true. In other words, it is only false if A is true but B is false. However, AB does not imply BA; the former is true when A is false but B is true, but the latter isn’t.

Yet this is what is at the heart of many conspiracy theories. For instance, a DHS report might say, that those on the fringe of the political right are motivated by the Obama government’s more permissive stance on stem cell research. Some draw the conclusion that this report implies that all who are troubled by Obama’s stance on this issue must be right-wing extremists. I could write this symbolically as follows: we have

member(e, s) → prop(e, p)

where member(e, s) means that e is a member of set s, and prop(e, p) means that e has property p. This symbolic equation cannot be reversed: it does not follow that prop(e, p) → member(e, s).

A closely related mistake is the confusion of the universal and existential operators. The existential operator (usually denoted with an inverted E, but I don’t have an inverted E on my keyboard, so I’ll just use a regular E), E(s, p) says that the set s has at least one member to which property p applies. The universal operator (denoted with an inverted A; I’ll just use a plain A), A(s, p) says that all members of set s have property p. Clearly, the two do not mean the same. Yet all too often, people (on both sides of the political aisle, indeed a lot of the politically correct outrage happens because of this) make this error and assume that once it has been asserted that E(s, p), it is implied that A(s, p). (E.g., a logically flawless statement such as “some blacks are criminals” is assumed to imply the racist generalization that all blacks are criminals.)

One might wonder why formal logic is not taught to would be politicians. I fear that in actuality, the situation is far worse: that they do know formal logic, and use it to their best advantage assuming that you don’t.

 Posted by at 12:27 pm
Apr 202009
 

I am watchin Deep Impact tonight, a ten-year old film about a comet impacting the Earth. Why the Canadian History Channel is showing this film is a good question. Future history? Imagined history?

But putting that question aside, the movie made me go to Wikipedia again, and I ended up (re-)reading several articles there relating to the issue of global warming and controversies surrounding it.

One thing that struck me (and not for the first time) is this: criticism of global warming theories are often dismissed by the assertion that these go “against the mainstream” or are “not supported by scientific consensus.”

And global warming is by no means the only area of science where such arguments are frequently invoked. Take two topics that I have become involved with. There is scientific consensus that the inadequacy of Einstein’s theory of gravitation to explain the rotation of galaxies and large scale features of the universe is due to “dark matter” and “dark energy”. Even though no one knows what dark matter (or dark energy) is made of, and no one actually detected any dark matter or dark energy ever, the idea is treated as fact. True, dark matter theory can explain a few things and even made a few minor (but nonetheless impressive) predictions, but that doesn’t necessarily make it true, and it certainly doesn’t make the theory the only kid on the block worth considering. Still, try proposing an alternative gravity theory: no matter how firmly rooted in real physics it is, you will be fighting an uphill battle.

Or take the Higgs boson. This hypothetical particle (often along with the graviton) is often portrayed as if it has already been detected. It hasn’t. Indeed, the only thing experiments have accomplished to date is that they excluded the possibility that the Higgs boson exist at nearly the two-σ level. There are also significant unresolved issues with the Higgs boson that put the theoretical validity of the idea into question. Yet the “scientific consensus” is that the Higgs boson exists, and if you try to propose a quantum field theory without the Higgs, well, good luck!

Just to be clear about it, I am not saying that the climate skeptics got it right, and for all I know, maybe there is dark matter out there in abundant quantities, along with Higgs bosons behind every corner. But not because this is what the “scientific consensus” says but because the theory is supported by facts and by successful predictions. Otherwise, the theory remains “just a theory”, as the creationist crowd likes to say… neglecting the inconvenient fact that, of course, the theory of evolution is supported by an abundance of facts and successful predictions.

 Posted by at 2:55 am
Apr 042009
 

The gravitational theory that I’ve been working on for some time with John Moffat is called STVG, or Scalar-Tensor-Vector Gravity. It grew out of Moffat’s investigation of Nonsymmetric Gravity.

There is also a phenomenological formula called MOND (MOdified Newtonian Dynamics) that effectively flattens out the acceleration curve at high radii from a point source. MOND is nothing more than a formula designed by its creator, Mordechai Milgrom, to solve a specific problem, namely the rotation curves of galaxies. It is not rooted in any theory, and in fact, it is known to contradict some; for instance, it violates the law of conservation of energy and momentum. This is why Jacob Bekenstein endeavored to create a relativistic theory called TeVeS, which fixes MONDs problems, while still gives the approximate MOND acceleration formula in the case of weak fields.

Both STVG and TeVeS are gravity theories, and both happen to incorporate tensor, vector, and scalar fields. Beyond that, however, there’s nothing in common between the two theories.

Unfortunately, many Wikipedians don’t know this, and try from time to time to merge the STVG and TeVeS articles. Hopefully not any longer… I just posted a long, fairly complete description of STVG on Wikipedia.

 Posted by at 11:48 pm
Mar 312009
 

Folks working on quantum computers are busy trying to make sure that entangled states remain entangled, because decoherence is death for a quantum computation. But now, Gross et al. showed that too much entanglement may not be a good thing: it can result in quantum computers that offer no improvements in efficiency over conventional computers.

 Posted by at 12:04 am
Mar 202009
 

In the Futurama movie, The Beast with a Billion Backs, one scene features a blackboard with two different proofs of the Goldbach conjecture. The Goldback conjecture is one of the oldest unsolved problems in mathematics. One of those problems, like Fermat’s last theorem or the 4-color problem of maps that is deceivingly easy to state and fiendishly hard to prove: that every even number greater than 4 can be expressed as the sum of two primes.

Just how intriguing this problem is, it’s well illustrated by the following plot that shows the number of ways an even number between 4 and 1 million can be split into a sum of two primes:

Golbach to 1000000

Golbach partitioning to 1000000

This plot, taken from Wikipedia, clearly shows that the results cluster along curves (asymptotes? attractors?) that follow some kind of a power law with an exponent between ~0.68 and ~0.77, and there may also be some fractal splitting involved, too. This plot is known as Goldbach’s comet. All I have to do is look at it to understand why many people find number theory endlessly fascinating.

 Posted by at 5:54 pm
Mar 152009
 

I often get questions about the Pioneer anomaly and our on-going research. All too often, the questions boil down to this: what percentage of the anomaly can <insert theory here> account for?

This is a very bad way to think about the anomaly. It completely misses the fact that the Pioneer anomaly is not an observed sunward acceleration of the Pioneer spacecraft. The DATA is not a measured acceleration; what is measured is the frequency of the spacecraft’s radio signal.

I already capitalized the word DATA in the previous paragraph; let me also capitalize the words MODEL and RESIDUAL, as these are the right terms to use when thinking about the Pioneer anomaly.

As I said above, the DATA is the Doppler measurement of the spacecraft’s radio frequency.

The MODEL is a model of all forces acting on the spacecraft including gravity, on-board forces, solar pressure, etc; all effects acting on the spacecraft’s radio signal, including the Shapiro delay, solar plasma, the Earth’s atmosphere; and all effects governing the motion of the ground stations participating in the communication.

The RESIDUAL is the error, the difference between the MODEL’s prediction of the Doppler measurement vs. the actual measurement. This RESIDUAL basically appears as noise, but with characteristic signatures (a diurnal and an annual sinusoid along with discontinuous jumps at the time of maneuvers) that suggest mismodeling.

The goal is to make this RESIDUAL “vanish”; by that, we mean that only random noise remains, any diurnal, annual, or maneuver-related signatures are reduced to the level of background noise.

The RESIDUAL can be made to vanish (or at least, can be greatly reduced) by incorporating new contributions into the MODEL. These contributions may or may not be rooted in physics; indeed, orbit determination codes typically have the ability to add “unmodeled” effects (basically, mathematical formulae, such as a term that is a quadratic or exponential function of time) to the MODEL, without regard to the physical origin (if any) of these effects.

Anderson et al. found that if they add an unmodeled constant sunward acceleration to the MODEL, they can make the RESIDUAL vanish. This is the result that has been published as the Pioneer anomaly.

If one has a physical theory that predicts a constant sunward acceleration, it is meaningful to talk in terms of percentages. For instance, one may have a physical theory that predicts a constant sunward acceleration with magnitude cH where c is the speed of light and H is Hubble’s constant at the present epoch; it then makes sense to say that, “using the widely accepted value of H ~= 71 km/s/Mpc, the theory explains 79% of the Pioneer anomaly,” since we’re comparing two numbers that represent the same physical quantity, a constant sunward acceleration.

However, note (very important!) that the fact that a constant sunward acceleration fits the data does not exclude alternatives with forces that are not constant or sunward pointing; the DATA admits many different MODELs.

Now let’s talk about the thermal recoil force. It is NOT constant and it is NOT sunward pointing. As we recompute this force, incorporating  the best thermal model that we can compute into the MODEL and re-evaluate it, we obtain a new RESIDUAL. There are, then, the following possibilities:

  1. Suppose that the new RESIDUAL is as free of a mismodeling signature as the constant acceleration model and that its magnitude cannot be reduced by adding any unmodeled effects (i.e., we reached the level of our basic measurement noise.) Does it then make sense to speak of percentages? OK, so the thermal recoil force is 30%, 70%, 130%, you name it, of the constant sunward acceleration. But the thermal recoil force is neither constant nor sunward, and by incorporating it into the MODEL, we got a different trajectory than the constant sunward acceleration cas. Yet the RESIDUAL vanishes, so the MODEL fits the DATA just as well.
  2. Suppose that the new RESIDUAL is half the original RESIDUAL at least insofar as the apparent mismodeling is concerned. What does this mean? Does this mean that the thermal recoil force and the resulting acceleration is half that of the constant sunward value? Most certainly not. Say it’s 65%. Now did we explain 50% of the anomaly (by reducing the RESIDUAL to 50%) or did we explain 65% of the anomaly (by producing a thermal recoil acceleration that’s 65% of the published constant sunward value?)

Instead of playing with percentages, it makes a lot more sense to do this: after applying our best present understanding insofar as thermal recoil forces are concerned, we re-evaluate the MODEL. We compute the RESIDUAL. We check if this residual contains any signatures of mismodeling. If it doesn’t, we have no anomaly. If it does, we characterize this mismodeling by applying various unmodeled effects (e.g., a constant sunward force, exponential decay, etc.) to check if any of these can characterize the RESIDUAL. We then report on the existence of a (revised) anomaly with the formula for the unmodeled effect as a means to consisely characterize the RESIDUAL. If this revised anomaly is still well described by a constant sunward term, we may use a percentage figure to describe it… otherwise, it’s probably not helpful to do so.

 Posted by at 3:56 am
Mar 102009
 

Einstein was right and Silberstein was wrong, but it’s beautifully subtle why.

Einstein’s point was that in a region free of singularities, the circumference of an infinitesimal circle should be 2π times its radius. He then showed that an infinitesimal circle perpendicular to, and centered around, the line connecting the two mass points in Silberstein’s solution has the wrong radius, hence the line connecting the two masses must be singular.

Silberstein countered by pointing out an error in Einstein’s derivation and then showing that a particular quantity rigorously vanishes along this line, implying that yes, the circumference of the infinitesimal circles in question do, in fact, have the right radii. This solution appeared in his paper in Physical Review in 1936.

Einstein then, in his published Letter that he wrote in reply to Silberstein’s paper, provided a technical argument discussing square roots and derivatives that, while correct, is not very enlightening. (It’s one of those comments that make perfect sense once you know what the devil he is talking about, but there’s no possible way you can actually understand what he’s talking about from the comment alone. Not that I haven’t been guilty of committing the same sin, even writing things that I myself wasn’t able to understand a few months later without going back to my notes or calculations. But I digress.)

By solving for the metric and obtaining an explicit formula, I was able to better understand this mystery.

First, it should be noted that the axis connecting the two mass points is also the origin of the radial coordinate, so the coordinate system itself is singular here even if there is no physical singularity. To verify that there is no physical singularity at this location, we can draw tiny, infinitesimal circles centered around this axis and check if they obey the rules of Euclidean geometry (one fundamental rule in curved spacetime is that in a small enough region, things always look Euclidean.) This was the basis of Einstein’s argument.

We can express the ratio of the angular and radial components of the metric in the form of an exponential expression, say exp(N) where N is some expression. Clearly, we must have N(r → 0) = 0 in order for the circumference of an infinitesimal circle centered around the axis to be 2π its radius.

As it turns out, N appears in the form [] + C where [] is just some complicated expression. It turns out that this part [] is constant along the axis, but its constant value is different between the two mass points vs. outside the mass points. I can draw the axis and the two mass points like this:

––––––––––––––––X––––––––––––––––X––––––––––––––––

The bracketed expression is constant along this line (not constant outside the line, such as above or below it but that’s not relevant now) but has one constant value between the two points (the two X’s) and another constant value on the two sides. The two constant values are not the same.

We can choose the constant C to make N vanish either between the two mass points or outside the two mass points but not both at the same time. When N vanishes outside, but not between, the two mass points, the singularity that remains serves as a “strut” holding the two mass points apart:

––––––––––––––––X================X––––––––––––––––

Or, when N vanishes outside, it’s as if the two mass points were attached by “ropes” to the “sphere at infinity”, hanging from there in static equilibrium as gravity pulls them together:

================X––––––––––––––––X================

So how come Silberstein was able to dismiss Einstein’s criticism? He was able to do so by making a particular choice of the sign of a square root that made N vanish both between and outside the two mass points. But this is where it helps to write N in the form that I obtained, symbolically as [] + C, with the bracketed part having one constant value along the axis between the two masses and another value outside. Sure we can make N vanish everywhere along the axis… by allowing C to have different values between vs. outside the masses.

This is fine along the axis… the masses themselves are singularities, so they represent a discontinuity anyway, so we are free to choose different integration constants on two sides of a discontinuity. This is indeed what Silberstein has done by choosing a particular sign of a square root in his formulation of the metric.

But, as Einstein pointed out, such a choice of sign leads to a discontinuity in derivatives. More explicitly, what happens is that the function N is defined not just along the axis but everywhere in spacetime… to get from a point on the axis between the two masses to a point outside, we need not go along the axis, we can leave the axis, go around the mass, and then return to the axis. As we do this, we encounter no singularity but the value of C must jump from one constant to another at some point. In other words, by removing the “strut” or “rope” singularity from the axis, we introduced a much worse singular “membrane” that separates regions of space.

One morale of this story is that when we use a coordinate system like polar coordinates that is not well defined at the origin, we must be ultra careful about that spot… since the coordinate system is singular here, this is where things can go wrong even if they appear perfect everywhere else.

 Posted by at 1:14 pm
Mar 102009
 

I’m reading about the debate in the 1930s between Einstein and Silberstein about the (non-)existence of a static two-body vacuum solution of the Einstein field equations.

Silberstein claimed to have found just such a solution as a special case of Weyl’s metric. However, he then concluded that the existence of an unphysical solution implies that Einstein’s gravitational theory has to be modified.

Meanwhile, Einstein dismissed Silberstein’s solution on two grounds. First, he claimed that there are additional singularities; second, he claimed that a solution that yields singularities is in any case not a proper solution of a field theory, so it certainly cannot be used to discredit that theory.

I disagree with Silberstein… just because there exist solutions that are unphysical does not unmake a theory. The equations of ballistics also yield unphysical solutions, such as cannonballs going underground or flying backwards in time… but it simply means that we chose unphysical initial conditions, not that the theory is wrong.

I also disagree with Einstein’s second argument though… field theory or not, some singularities can be quite useful and physically meaningful, be it, say, the “point mass” in Newton’s theory, the “point source” in electromagnetism, or, well, singularities in general relativity representing compact (point) masses.

But both these issues are more philosophy than physics. I am more interested in Einstein’s first argument… is it really true that Silberstein’s solution yields more than two singularities?

That is because when I actually calculate with Silberstein’s metric, I find regular behavior everywhere except at the two singular points. I see no sign whatsoever of the supposed singular line between them. What am I missing?

 Posted by at 12:49 am
Feb 162009
 

I received some sad news yesterday from Hungary: my high school math teacher, Gusztáv Reményi, died last week, at the age of 88. He was a very kind teacher. Our class was a specialized mathematics class, and we were supposed to be the best in the country. In this class, being good at math didn’t just mean that, say, you got sent to national math competitions; you were expected to win them. Perhaps this made Mr. Reményi’s job easier, but I suspect that he would have done well with less talented pupils, too, if not because of his teaching style then due to his personality. If you met him and remembered nothing else, you’d have remembered his smile. I last met him a few years ago, at our high school reunion. He was old, he was frail, but the huge smile was still there, just as I remembered.

 Posted by at 3:54 pm
Feb 122009
 

Once again, I am reading an interesting paper on ArXiv.org (doesn’t matter which one, it wasn’t that interesting) and I notice that the author is a physicist from some Iranian university. ArXiv.org has many papers from Iran. No wonder that nation was able to launch a satellite and is working on a nuclear (weapons?) program, apparently with every hope of success. I am not sympathetic towards the regime of the ayatollahs, but the fact that Iran is not as black-and-white as some would like us to believe must be recognized. I also suspect that another fact, itself somewhat hard to reconcile with the picture of a monolithic, intellectually repressive theocracy, namely that as of 2007, 23 million out of Iran’s 66 million inhabitants had Internet access (according to the CIA World Factbook), has a great deal to do with the success and competence of Iranian physicists.

 Posted by at 2:45 am
Feb 072009
 

I remain troubled by this business with black holes.

In particular, the zeroth law. Many authors, such as Wald, say that the zeroth law states that a body’s temperature is constant at equilibrium. I find this formulation less than satisfactory. Thermodynamics is about equilibrium systems to begin with, so it’s not like you have a choice to measure temperatures in a non-equilibrium system; temperature is not even defined there! A proper formulation for the zeroth law is between systems: the idea that an equilibrium exists between systems 1 and 2 expressed in the form of a function f(p1, V1, p2, V2) being zero. Between systems 2 and 3, we have g(p2, V2, p3, V3) = 0, and between systems 3 and 1, we have h(p3, V3, p1, V1) = 0. The zeroth law says that if f(p1, V1, p2, V2) = 0 and g(p2, V2, p3, V3) = 0, then h(p3, V3, p1, V1) = 0. From this, the concept of empirical temperature can be obtained. I don’t see the analog of this for black holes… can we compare two black holes on the basis of J and Ω (which take the place of V and p) and say that they are in “equilibrium”? That makes no sense to me.

On the other hand, if you have a Pfaffian in the form of dA + B dC, there always exists an integrating denominator X (in this simple case, one doesn’t even need Carathéodory’s principle and assume the existence of irreversible processes) such that X dY = dA + B dC. So simply writing down dM – Ω dJ already gives rise to an equation in the form X dY = dM – Ω dJ. That κ and A serve nicely as X and Y may be no more than an interesting coincidence.

But then there is the area theorem such that dA > 0 (just like dS > 0). Is that another coincidence?

And then there is Hawking radiation. The temperature of which is proportional to the surface gravity, T = κ/2π, which is what leads to the identification S = A/4. Too many coincidences?

I don’t know. I can see why this black hole thermodynamics business is not outright stupid, but I remain troubled.

 Posted by at 9:50 pm