jump to navigation

Vector Spherical Harmonics May 26, 2010

Posted by Sarah in Uncategorized.
Tags: , , ,

Spherical harmonics are the angular part of the solution to Laplace’s equation on the sphere: they have the property that
\Delta Y_{lm} = \frac{l (l+1)}{r^2} Y_{lm}

They form an orthonormal basis in Hilbert space and have many nice properties. In particular, the space of spherical harmonics is an irreducible representation of the group of rotations SO(3).

Vector spherical harmonics are an extension of the concept for use with vector fields. We define three vector spherical harmonics,
\mathbf{Y_{lm}} = Y_{lm} \hat{\mathbf{r}}
\mathbf{\Psi_{lm}} = r \nabla Y_{lm}
\mathbf{\Phi_{lm}} = \vec{r} \times \nabla Y_{lm}
These are orthogonal just like the usual spherical harmonics, and they allow every vector field to be expanded in vector spherical spherical harmonics.

They also turn out to be useful in the study of magnetostatics (that is, the study of static magnetic fields) as in this paper by Barrerra, Estevez, and Giraldo (Here.)

The gist: expanding functions on the sphere in spherical harmonics makes it possible to replace the Poisson equation with an ordinary differential equation. This means that if we define a charge distribution and require it to be bounded in extent, we can determine the potential outside the charge distribution. The key is canceling out the angular component to simplify the partial differential equation.

Extending this notion to vector fields allows us to do the same thing with differential equations involving the gradient (and Laplacian and curl etc), expanding in vector spherical harmonics and letting the angular components cancel. This will allow us to calculate magnetic multipoles analogous to electrical multipoles. When there is no charge, the magnetic induction field is the gradient of some function, the magnetic scalar potential.
\mathbf{B} = \nabla \mathbf{\Phi}.
The potential can be expanded in spherical harmonics; the coefficients are called the magnetostatic multipole moments.


Representation Theory: Basics and Heisenberg Representation (3) May 25, 2010

Posted by Sarah in Uncategorized.
add a comment

We return to our regularly scheduled programming — more notes from the seminar.

The next example is the Heisenberg group. Recall that this is the group of 3-by-3 upper triangular matrices with ones on the diagonal. What are the irreducible representations of this group H?

Let (V, \omega) = (\mathbb{F}_p \times \mathbb{F}_p; \omega(v, v') = det(v v'). V is denoted a symplectic vector space, and \omega is its form. The Heisenberg group H is
H = V \times \mathbb{F}_p
endowed with the group law
(v, z) * (v', z') = (v + v', z + z' + 1/2 \omega (v, v'))

One-dimensional representation:
H \to V \to \mathbb{C}^\times
\rho_\theta(h) = \theta(Pr(h))
where \theta is the function from V to \mathbb{C}^\times. So this is just the projection.

Stone-von Neumann Theorem:

There exists a unique up to isomorphism irreducible representation $latex(\pi, H, \mathcal{H})$ such that
\pi |_z = \psi(z) Id|_{\mathcal{H}}

Lagrangian models.

\mathcal{H}_L = \mathbb{C}(L \ H, \psi) = \{ f: H \to \mathbb{C}, f(z l h) = \psi(z) f(h)
$dim \mathcal{H}_L = p$
We can think of L as a line through the space V.

\pi_L = H \to \mathcal{H}_L
[\pi_L(h') f](h) = f(h h')

Claim (\pi_L, H, \mathcal{H}) irreducible.
We obtain (p-1) \times p-dimensional irreducible representations
P^2 * 1^2 + (p-1)*p^2 = p^3 = |H|

F_{M, L} \in Hom_H(\mathcal{H}_L, \mathcal{H}_M) \simeq \mathbb{C}
F_{M, L}[f] (h) = \sum_{m \in M} f(mh)

Quote of the Day: William Vollmann May 25, 2010

Posted by Sarah in Uncategorized.
add a comment

I spent a 12-hour train ride through all the twilight mountains of upstate New York, reading William Vollmann’s wonderful book about riding freight trains. It is now time for me to read all his books — I believe the man is a kindred spirit.

Every time I surrender, even necessarily, to authority which disregardingly or contemptuously violates me, so I violate myself. Every time I break an unnecessary law, doing so for my own joy and to the detriment of no other human being, so I regain myself, and become strong in the parts of me that the security man can never see.

Even better — it seems Will Wilkinson agrees. And pretty much any friend of Will Wilkinson is a friend of mine.

(N.B. I am fighting a long-standing bad habit here — it is my aim not to mention politics even in the off-topic posts on this blog. Talking politics for me is the way tequila is for other people: it always seems like a good idea at the time and a bad idea in the morning. But this is literature, not politics, and therefore kosher — or so I tell myself.)

A Yardstick for Smell: Thinking in PCA May 22, 2010

Posted by Sarah in Uncategorized.
Tags: ,
add a comment

I mention this Discover article because it’s really an applied math story being popularized as a biology story. Noam Sobel and his team from Weitzmann made a database of thousands of smell-producing molecules, along with their properties (like size, or how tightly the molecules pack.) They then ran Principal Components Analysis on the data to find the critical properties, and observed that a single dimension or score encapsulated most of the variation — and even could predict which smells humans would find more and less appealing.

What I want to harp on is the use of PCA. To those of us who actually do math and science full time, it’s a pretty familiar, standard technique. It’s over a hundred years old. It’s very simple. It doesn’t even work in the difficult cases (a lot of modern random matrix theory is about identifying when PCA doesn’t work.) But I don’t think Principal Components Analysis (or multidimensional eigenvector analysis generally) gets the attention it deserves in the popular imagination. Here, given a glut of data, is a tool that can tell you what are the most useful ways to measure it. We’re used to science giving us measurements; but this is science telling us what measurements are important. Sobel’s research isn’t really a biological discovery at all; it’s a data analysis discovery.

And when you’re used to thinking in the language of dimensionality reduction, suddenly you see when it’s needed but missing. A gap, like “Oh, this really should be a dimensionality reduction problem.” For instance, I was reading some of Simon Baron-Cohen’s research about autism; for the unfamiliar, his theory is that people lie on a spectrum between “systematizing” and “empathizing” thoughts, with autistics on the “systematizing” extreme, much better at logical puzzles than interpersonal communication. I was reading his papers and I immediately thought “Something’s missing here.” If I wanted to know whether a single axis distinguished autistics from non-autistics, well, I’d look at all kinds of neurological and psychological properties, get a big cloud of data from autistics and non-autistics, and see if most of the variation was due to this empathizing/systematizing business.

Now, I’m not remotely a psychologist, but it looks like Baron-Cohen doesn’t do that. He gives out a survey, observes that autistics are more systematizing and less empathizing than non-autistics, and essentially says “Ta-Da!” And I go, “Huh?” How on earth does he know that this is the main difference between autistics and non-autistics? Where are all the correlations and variances? Now, maybe that’s standard practice in psychology. I’ve read a bit more in the social sciences, and it’s standard practice there. But to me, this kind of research is just crying out for a little data analysis. I think there are whole areas of research where people aren’t yet thinking in PCA. Maybe that should change.

Representation Theory: Basics and Heisenberg Representation (2) May 21, 2010

Posted by Sarah in Uncategorized.
add a comment

More notes from Shamgar Gurevich’s lectures bringing representation theory to the applied-math populace.

We start with an example. Let G = S_4 = Aut\{a, b, c, d\}. Denote the irreducible representations of G by G^v. There are 5 of these, one for each conjugacy class (in the case of the symmetric group, a conjugacy class is a partition.)

Constructions: S_4 goes into the complex numbers over the tetrahedron, which is V_0 \oplus V_{const}, where V_0 is the three-dimensional space of the front face of the tetrahedron. The representation is
\rho_3: S_4 \to GL(V_0).
Why is \rho_3 irreducible?

General tool: the intertwining number, \langle \rho, \pi \rangle.

Proposition: \rho is irreducible iff \langle \rho, \rho \rangle =1.
Proof: \rho = \oplus m_i \rho_i.
1 =  = \sum m_i^2
Since these are integers, the equation is true iff there is only one subrepresentation.

In our case,
\rho = \rho_3 + \rho_1
\langle \rho, \rho \rangle = \langle \rho_3, \rho_3 \rangle  + 2\langle \rho_3, \rho_1 \rangle  +\langle \rho_1, \rho_1 \rangle .
We have X, Y as G-sets: $ latex (\pi_X, G, \mathbb{C}(X))$.

\langle X, Y \rangle is the number of G orbits in X \times Y.
In our case, \langle \rho, \rho \rangle= 2 implies \langle \rho_3, \rho_3 \rangle= 1.

Proof of proposition:
K: Hom(\mathbb{C}(X), \mathbb{C}(Y)) \to \mathbb{X \times Y}
both groups are isomorphic to G. This implies
Hom_G (\mathbb{C}(X), \mathbb{C}(Y)) \to \mathbb{C}(X \times Y)^G.

Another 3-dimensional representation:
\pi_3=  sgn \otimes \rho_3: S_4 \to GL(V_0).
Clearly this is also irreducible.
How do we know it’s not equivalent to \rho_3?

Well, take three function f_1, f_2, f_3 defined on the triangle (all entries have to add to 1.) f_1 is 0 on both vertices a, b, f_2 is 1 on both, and f_3 is 1 on a and -1 on b.
These form a basis. The trace of \rho_3 in this basis is 1 + 1 + -1 = 1, and the trace of \pi_3 in this basis is sign(a, b) = -1. So they are not identical.

(To be continued — I’m posting from Montreal, so I don’t have long blocks of computer time.)

King of the Road May 18, 2010

Posted by Sarah in Uncategorized.

I’m all done!
Thesis defense seemed to go smoothly — I’m very happy with it. And this is the end of all my work in college! Now I pack and clean my room, because I’m going to Montreal tomorrow with some friends. It is literally the cheapest possible vacation, and it’s going to be awesome.
We start with a ten-hour train ride at the crack of dawn, and I have the perfect book for it: William Vollmann’s Riding Towards Everywhere, a book about … riding the rails.

Montreal! There will be sightseeing, and hostel adventures, and poutine. Then I go back to school to see Reunions and graduate; I’ll sing Old Nassau one more time and get all misty, and then it’s back home, where I will do nothing but run, lift, and read math. It’s going to be a good summer.

Representation Theory: Basics and Heisenberg Representation (1) May 17, 2010

Posted by Sarah in Uncategorized.
add a comment

Here’s a new installment of lecture notes from Shamgar Gurevich’s seminar on representation theory. (This was a small group of applied math people, professors and grad students and me, looking at applications of representation theory to the problem of cryo-electron microscopy. It’s very pretty stuff and well taught.)

The Discrete Fourier Transform (DFT) is defined as
1/\sqrt{P'} (e^{\frac{2 \pi i}{p} \omega}) : \mathbb{C}(\mathbb{F}_p \to \mathbb{C}(\mathbb{F}_p
We want to find a natural basis B_{DFT} of eigenvectors, and explain how to compute the DFT fast. Related is the following theorem:

Fix 1 \neq \psi: \mathbb{F}_p \to \mathbb{C}^\times an additive character. There is a unique up to isomorphism irreducible representation (\pi, H, \mathcal{H}) such that
\pi |_{\mathbb{Z}(H)} = \psi(\cdot) Id_{\mathcal{H}}.

Let’s show first:
R \to V = \mathbb{C}(G) \simeq \oplus_{\rho \in Irr(G)} dim \rho V_{\rho}.

In particular,

|Irr(G)| < \infty
|G| = \sum_{\rho \in Irr(G)} dim(\rho)^2
To show this, let us answer
G \to_\pi V =  \oplus_{\rho \in Irr(G)} m_\rho V_{\rho}.
How to compute m_\rho?

\langle G, \pi \rangle = dim Hom_G(G, \pi)

Properties of intertwining numbers \langle G, \pi \rangle:

1. G \simeq \hat{G}, \pi \simeq \hat{\pi}; \implies \langle G, \pi \rangle = \langle \hat{G}, \hat{\pi}; \rangle

2. \langle G_1 \oplus G_2, \pi \rangle = \langle G_1, \pi \rangle + \langle G_2, \pi \rangle

3. G, \pi \in Irr(G) \langle G, \pi \rangle = \delta_{G, \pi}

1 is clear; 3 is proven by Schur’s Lemma.
To prove 2, observe
Hom(V_1 \oplus V_2, W) = Hom(V_1, W) \oplus Hom(V_2, W)
Hom_G(V, W) = Hom(V, W)^G
(V \oplus W)^G = V^G \oplus W^G

The result is, if
\pi \simeq \oplus  m_\rho * \rho
m_\rho = \langle \rho, \pi \rangle

The application to the first theorem is that if
R: G \to \mathbb{C}(G) and given a representation (\rho, G, V) then
\langle R, \rho \rangle = dim V.

Proof: Hom(\mathbb{C}(G), V) = V

We are now ready to prove the next theorem: the number of irreducible representations of G equals the number of conjugacy classes of G. The idea is an isomorphism between the “geometric side” and the “spectral side.”

\mathbb{C}(G) \to \oplus_{\rho \in Irr(G)} Hom(V_\rho, V_\rho)
This is an isomorphism:
R(f) = \sum_{\rho \in Irr(G)} dim \rho * \rho(f)
= \sum_{g \in G} f(g) \rho(g)
Assume \rho(f) = 0 \forall \rho. Then R(f) = 0. But R(f) \delta_e = f.
The left hand side are the functions constant on conjugacy classes of G, while the right hand side are
Hom_G (V_\rho, V_\rho).
In particular, the number of conjugacy classes is the number of irreducible representations.
And G is abelian if every irreducible representation has dimension 1.

I have math friends! May 17, 2010

Posted by Sarah in Uncategorized.
add a comment

Apologies for light blogging — I’ve been in finals and thesis mode. (I defend today!) I have some nice representation theory notes lined up, I promise. Meanwhile I note that two very nice math people also have brand new blogs. Brooke and Jimmy went to school at UChicago, which makes them cool. Jimmy writes a lot of pedagogical stuff, which is pretty damn important — imagine if everyone who was supposed to learn multivariable calculus actually learned it and liked it! Brooke seems to have the same quals prep plan as me (at least for the algebra/topology stuff.)

Read ’em!

I Am Trying To Break Your Heart May 10, 2010

Posted by Sarah in Uncategorized.
add a comment

I am a Chicago girl, and so it is my duty to present you with this incredible soul cover of Wilco’s “I Am Trying To Break Your Heart” by JC and the Uptown Sound.

They also have a clever little blog here which is improving my musical taste as we speak.

We Fine May 7, 2010

Posted by Sarah in Uncategorized.
add a comment

Last get-together for the math majors, in Fine Hall, and a final group photo.

I’m in front in the blue top.
I’m going to miss this place, and I’m really going to miss the people.