Speculative idea: the genus of an eigenfunction and Kac’s “can one hear the shape of a drum?”

This post will be of a different nature to my other posts: I want to talk about a speculative idea that: (i) I had while sitting in a talk; (ii) I haven’t worked through at all; (iii) involves higher power mathematics than what I have at my disposal right now; and (iv) may be stupid, wrong, too naive, or already well-known.

These kind of ideas occur to theorists all the time and, due to resource constraints (i.e., time, energy, money, etc.), seldom get mentioned, except to, say, a conference speaker, and are quickly forgotten by both parties by the time the conference is over.

However, with access to this new medium I hope to broadcast these ideas to more people, some of whom may be experts on the maths required, and maybe, just maybe, someone will be inspired to check it out.

Dumb questions are especially encouraged in these posts because, after all, I’m posting about a dumb question!

Please read on for my first speculative idea post.

1. Nodal domains

Recently Uzi Smilansky gave a talk to the quantum dynamics group at Royal Holloway entitled “Can one count the shape of a drum?” (see here for the paper the talk was based on). This is heady stuff! Uzy is a real physicist of the highest caliber and he is able to understand subtle problems with powerful physical arguments that, as yet, cannot be studied at a mathematical level of rigour.

Uzy is responsible for a kind of revolution in our understanding of the simple (Schrödinger) wave equation: whereas a lot of research energy has been spent on understanding the eigenvalue distribution of the 2D time-independent Schrödinger wave equation

\displaystyle  -\frac{\hbar^2}{2m}\left(\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}\right)\psi(x,y) +V(x,y)\psi(x,y) - E\psi(x,y) = 0, \ \ \ \ \ (1)

a great deal of interesting open questions remain!

Uzy has a new approach to the question, first put by Victor Kac: “can one hear the shape of a drum?” What does this mean? Suppose you knew the all of the eigenvalues {\{E_n\}_{n=0}^\infty} of the Schrödinger equation (or, as originally posed, the simple wave equation) for a particle in an arbitrarily shaped 2D box. Could you reconstruct the shape of the box knowing just the eigenvalues? The answer, is no. Although, the simplest counterexamples are nonconvex. In fact, one of the reasons the question is interesting is precisely because it could so easily have been yes: there are large families of domains which can be uniquely reconstructed from knowing the spectrum of the Schrödinger equation in that domain.

To understand Uzy’s contribution we need to define a nodal domain. Let {\psi_n(x,y)} be an eigenfunction of (1). A nodal domain is a connected region where {\psi(x,y)} is positive or negative. The boundaries between nodal domains are given by the zeros of {\psi(x,y)}, i.e., by the (non simply) connected curve

\displaystyle  \psi(x,y) = 0. \ \ \ \ \ (2)

Write {\nu_n} for the number of nodal domains for the {n}th eigenfunction. Uzy’s idea is as follows: suppose you know both {E_n} and {\nu_n}, then you can reconstruct the domain. He presents a great deal of evidence for this conjecture and has gone on to study the properties of nodal domains in great detail.

In his talk, Uzy mentioned a problem: how do you count nodal domains? This seems like a stupid question when you first encounter it; after all, don’t you just write down the eigenfunction {\psi_n} and count the nodal domains? Well, yes and no: the problem is that you need to count 2D connected regions and there is no elegant mathematical way to represent this count short of running a 2D search algorithm. Contrast this with the situation in 1D where nodal domains can be easily counted by instead counting the zeros of {\psi_n} and these, in turn, can be studied via the poles of {1/\psi_n}, upon which the machinery of complex analysis can be brought to bear. Additionally, nodal domains are pretty sensitive to boundary conditions: they easily split if the slightest perturbation is added. Thus they are not very robust. Arguably eigenvalues are similarly non-robust but they aren’t that sensitive to the shape of the perturbation: under a perturbation all that can happen is that degeneracies get lifted. Of course, you might well argue that this holds true for the distribution {\nu_n}, and to this I have no real reply other than we’d like to have a robust way to count them which somehow accounts for perturbations. This obsession with perturbations is because when you solve these things numerically, or approximate them theoretically, the errors induced can be modelled as perturbations. One would like to have a way to track how errors propagate into {\nu_n}.

2. Introducing the genus

I think it was Michael Nielsen who once said to me: “if the answer is ugly, then maybe the question was wrong” (the original source of this quote, and accuracy, is unknown to me). I think that the difficulty in counting nodal domains that I mentioned in the previous section could be a symptom of the fact the question asked is slightly wrong. I really don’t want to sound like I’m criticising the whole nodal domain approach, which is a truly wonderful idea, not least because Uzy has gotten actual results. Rather, I would like to suggest a modification of the question which might be easier to study but which might lose some of the expressive power of nodal domains.

The idea is based one of the fundamental ideas involved in the study of algebraic curves such as the elliptic curve

\displaystyle  y^2 = x^3 + ax + b \ \ \ \ \ (3)

which can be written implicitly as

\displaystyle  f(x,y) = 0. \ \ \ \ \ (4)

The idea is that when {x} and {y} are promoted to complex numbers then (4) becomes an equation over {\mathbb{C}^2}. But when we regard {\mathbb{C}^2} as {\mathbb{R}^4} in the natural way then (4) gives rise to two real equations over {\mathbb{R}^4}:

\displaystyle  \mbox{Re}(f)(\mbox{Re}(x), \mbox{Im}(x), \mbox{Re}(y), \mbox{Im}(y)) = 0 \ \ \ \ \ (5)


\displaystyle  \mbox{Im}(f)(\mbox{Re}(x), \mbox{Im}(x), \mbox{Re}(y), \mbox{Im}(y)) = 0 \ \ \ \ \ (6)

and represents a two-dimensional surface in {\mathbb{R}^4}. (See here for some example visualisations.)

If {\psi(x,y)=0} is the equation of the boundary between nodal domains for an eigenfunction {\psi(x,y)} then we can, as above, regard {\psi(x,y)=0} as defining an equation over {\mathbb{C}^2}, or defining a two-dimensional surface in {\mathbb{R}^4}. The actual boundary is then the intersection of this surface with the two-dimensional plane defined by {\mbox{Im}(x) = 0} and {\mbox{Im}(y) = 0}.

Given that we can interpret nodal boundaries in this way, i.e. as 2D surfaces, then we can try ask questions about these surfaces which should be topologically robust. The first question which comes to mind, inspired by the fact that all closed 2D surfaces have been classified up to homeomorphism, is to ask what is the genus {g} of the (projective version of the) equation defined by {\psi(x,y)=0}. (By projectively equivalent I mean the equation when interpreted over the complex numbers with infinity adjoined, so the space becomes closed.) This quantity arises from the fact that every 2D surface is homeomorphic to a sphere with {g} handles.

So that’s it: the idea is to study not the sequence {\{E_n, \nu_n\}_{n=0}^\infty} of eigenvalues and numbers of nodal domains, but rather to study {\{E_n, g_n\}_{n=0}^\infty} of eigenvalues and corresponding genuses of the surface defined by {\psi_n(x,y) = 0}. (Note that it’s probably important that (1) be defined over all {\mathbb{R}^2} to begin with, and we mimic say, boxes, with infinite potentials {V}.)

Presumably (I guess!) the genus is easier to understand in the presence of perturbations than the number of nodal domains. The problem is I have no idea how to compute the genus in terms of the original differential operator, so it’s not really any gain whatsoever. (Indeed, it may be impossible to calculate practically or even approximate!)

Any comments or suggestions or contributations, as always, are most welcome!


2 Responses to Speculative idea: the genus of an eigenfunction and Kac’s “can one hear the shape of a drum?”

  1. Steve Flammia says:

    Hi Tobias,

    Interesting post. I look forward to reading Uzi’s paper at some point. One way which might be useful in studying nodal domains more directly as level curves in R^2 would be using Morse theory. Morse theory works quite well for “generic” situations, and of course most potentials that we deal with in physics are not “generic”. So probably this approach has plenty of limitations. But perhaps it is worth thinking about.

    In particular, it is a standard calculation in Morse theory to find the Euler characteristic of a surface. I don’t know anything about elliptic curves, so I don’t know if this is related to the genus as a surface in R^4, but it does seem to yield potentially useful information.

  2. tobiasosborne says:

    Dear Steve,

    Many thanks for your comment!

    That sounds like a great idea (and, really, a more direct answer to the question of how to count nodal domains). I know even less about morse theory than algebraic curves (and I know very little about them)… But skimming through Hirsch’s evergreen book (http://books.google.co.uk/books?id=iSvnvOodWl8C) certainly does inspire hope that this is a better approach!

    I didn’t realise Morse theory was such a powerful computational tool.

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

%d bloggers like this: