Blog

How many index can you cram into an equation

Index notation, where a summation over some set is indicated by subscripts and superscripts, can be a curse in physics, as things can go crazy pretty fast for some systems. This usually occurs either when the equations require a large amount of them due to the sheer number of terms or just their variety, in the case of the most esoteric physics.

For instance, he's an example from my own master thesis, where some calculation required me to compute Taylor terms up to the 4th order :

\begin{align} \psi (x_2,t + \epsilon ) &=\sqrt{\frac{m}{2\hbar\pi i\varepsilon}} \int_{-\infty}^{\infty} d\eta\ e^{\frac{im}{2\hbar} g_{ab}(x_2) \frac{\eta^a\eta^b}{\varepsilon}}\nonumber\\ &(1 - \frac{im}{2\hbar\varepsilon}(\lambda \eta^\mu \eta^\nu \eta^\rho g_{\mu\nu,\rho})- \frac{m^2}{8\hbar^2\varepsilon^2}(\lambda^2 \eta^\mu \eta^\nu \eta^\rho \eta^\sigma ( g_{\mu\nu,\rho\sigma}+ \eta^\tau \eta^\chi g_{\sigma\tau,\chi}g_{\mu\nu,\rho})) +\mathcal{O}(\eta^7))\nonumber\\ &(1 +\frac{im}{2\hbar\epsilon} g_{ab}(\tilde{x}) C^{ab}- \frac{m^2}{8\hbar^2\epsilon^2} g_{ab}(\tilde{x})g_{pq}(\tilde{x})C^{ab}C^{pq}+\mathcal{O}(C^3)) \nonumber\\ &(\psi(x_2,t)+ \eta^\alpha \psi(x_2,t)_{,\alpha}+ \frac 12 \eta^\alpha \eta^\beta \psi(x_2,t)_{,\alpha\beta} +\mathcal{O}(\eta^3)) \nonumber\\ &(\sqrt{g(x_2)} + \eta^\aleph \sqrt{g(x_2)}_{,\aleph}+ \frac 12 \eta^\aleph \eta^\beth \sqrt{g(x_2)}_{,\aleph\beth}+ \mathcal{O}(\eta^3)) \end{align}

with some hebrew letters added for fun. Index notation is a topic of some humor, for instance from this parody paper :

We use the following notation: Greek letters for vile spinor indices, Roman letters for isospin indices, Cyrillic letters for vector indices, Hebrew letters for 10-dimensional vector indices, Roman numerals for Dirac spinor indices, and radiation therapy for Hodgkin's disease.

Due to this, indexes tend to be dropped, assumed to be implicit by the choice of objects used in the equation. This happens a lot for spinors, in particular.

Let's now see what we could do with the opposite approach, to use as many indexes as possible. Take for example the Dirac equation, in pure unindexed form :

$$(i\hbar \partial\!\!\!/ - mc) \psi = 0$$

The Feynman slash notation $\partial\!\!\!/$ implies here a summation over two spacetime indexes with Dirac's gamma matrices, $\gamma^\mu \partial_\mu$ (this isn't strictly speaking two spacetime indexes since $\gamma$ isn't a tensor but a soldering form, but it is close enough to it). The equation then becomes

$$(i\hbar \gamma^\mu \partial_\mu - mc) \psi = 0$$

For a total of a single index. From there, we'll have to add the actual spinor indices. There are four types of spinor indices, all noted with upper case letters : the basic spinor $\psi^A$, the dual space spinor $\psi_A$, the complex conjugate spinor $\psi^{\dot A}$ and the complex conjugate dual spinor $\psi_{\dot A}$. The spinor of the Dirac equation is actually none of those, but a vector of spinors called the Dirac bispinor.

$$\psi = \begin{pmatrix}\xi^{A}\\\chi^{\dot A}\end{pmatrix}$$

$\xi$ is the left-handed spinor while $\chi$ is the right-handed spinor. In that notation, the gamma matrices become a matrix of the solder forms on those spinors.

$$\gamma^\mu = \begin{pmatrix} 0 && \sigma^\mu_{A\dot A}\\ \bar \sigma^\mu_{A\dot A} && 0 \end{pmatrix}$$

This allows for instance the usual formula to map a spinor to a vector via $\psi^{\dot A} \sigma^\mu_{A\dot A} \psi^A$. The Dirac equation then becomes

\begin{eqnarray} i\hbar (\sigma^\mu_{A\dot A} \partial_\mu - mc I_{A\dot A})\chi^A &=& 0\\ i\hbar (\bar \sigma^\mu_{A\dot A} \partial_\mu - mc I_{A\dot A}) \xi^A &=& 0 \end{eqnarray}

$I$ simply being the identity matrix. We can set the mass to $0$ and end up with the Weyl equation for a single right-handed spinor with no loss of indexes.

$$i\hbar (\sigma^\mu_{A\dot A} \partial_\mu - mc I_{A\dot A})\chi^A = 0$$

Up to three indexes. Now we go onto gauge fields : we replace the partial derivative $\partial_\mu$ with the gauge derivative $\nabla_\mu$. Gauge derivatives are Kozsul connection on a principal bundle and as such are fairly complicated, something of the form (up to some physical factors)

$$\nabla_\mu = \partial_\mu + \tau_\alpha^{ij} A_\mu^\alpha$$

Where $A$ is the gauge field and $\tau$ is the representation of the gauge group. In our case, we'll want to see the gauge derivative of a product of gauge groups, such as $\mathrm{SU}(3) \times \mathrm{SU}(2) \times \mathrm{U}(1)$, for which it will be a sum of the form

$$\nabla_\mu = \mathrm{I}^{ij} \mathrm{I}^{kl} \mathrm{I}^{mn} \partial_\mu + \tau_\alpha^{ij} G_\mu^\alpha \mathrm{I}^{kl} \mathrm{I}^{mn} + \mathrm{I}^{ij} \sigma_\alpha^{kl} W_\mu^\alpha \mathrm{I}^{mn} + \mathrm{I}^{ij} \mathrm{I}^{kl} u^{mn}_\alpha B_\mu^\alpha $$

$G$ is the gluon field, with $\tau$ the Gell-Mann matrices, $W$ is the weak isospin field, with $\sigma$ the Pauli matrices, and $B$ is the weak hypercharge field, with $u$ simply being some one dimensional matrix (we could remove its indices but let's keep it). The spinor field itself is a product of the original spinor space $\mathbb C^2$, $\chi^A_{ikm}$, each index corresponding to a different charge of the standard model. If we drop the mass term (to avoid getting too long), the new equation then becomes something like

$$i\hbar (\sigma^\mu_{A\dot A} (\mathrm{I}^{ij} \mathrm{I}^{kl} \mathrm{I}^{mn} \partial_\mu + \tau_\alpha^{ij} G_\mu^\alpha \mathrm{I}^{kl} \mathrm{I}^{mn} + \mathrm{I}^{ij} \sigma_\alpha^{kl} W_\mu^\alpha \mathrm{I}^{mn} + \mathrm{I}^{ij} \mathrm{I}^{kl} u^{mn}_\alpha B_\mu^\alpha) ) \chi^A_{ikm} = 0$$

Although we could have arbitrarily many indexes here, for a reasonable theory we are up to 10. We can then add the various flavors of particles. Considering every particle as a variation of the triplet, doublet and singlet of the standard model, there are six of them : the three couples of leptons and their neutrinos $(e, \nu_e)$, $(\mu, \nu_\mu)$ and $(\tau, \nu_\tau)$, and the three generations of quarks, $(u,d)$, $(c,s)$ and $(b,t)$. Without the mass term, we can simply add the flavor index (otherwise, the mass term would depend on both the flavor and the charges)

$$i\hbar (\sigma^\mu_{A\dot A} (\mathrm{I}^{ij} \mathrm{I}^{kl} \mathrm{I}^{mn} \partial_\mu + \tau_\alpha^{ij} G_\mu^\alpha \mathrm{I}^{kl} \mathrm{I}^{mn} + \mathrm{I}^{ij} \sigma_\alpha^{kl} W_\mu^\alpha \mathrm{I}^{mn} + \mathrm{I}^{ij} \mathrm{I}^{kl} u^{mn}_\alpha B_\mu^\alpha) ) \chi^A_{ikmf} = 0$$

Now we can add curved spacetime to the mix. This changes two things : we replace spacetime indices with frame indices (and add the frame field), and we also need another connection, this time the $\mathrm{SO}(1,3)$ connection of the frame bundle. $$i\hbar (e^\mu_a \sigma^a_{A\dot A} (\mathrm{I}^{ij} \mathrm{I}^{kl} \mathrm{I}^{mn} \partial_\mu + \tau_\alpha^{ij} G_\mu^\alpha \mathrm{I}^{kl} \mathrm{I}^{mn} + \mathrm{I}^{ij} \sigma_\alpha^{kl} W_\mu^\alpha \mathrm{I}^{mn} + \mathrm{I}^{ij} \mathrm{I}^{kl} u^{mn}_\alpha B_\mu^\alpha) ) \chi^A_{ikmf} = 0$$ [...]

While rarely done, it is entirely possible to write quantum states in index notation, since they are vectors themselves. For instance, take the two wavefunctions $\Psi$ and $\Phi$. We can write their product as

$$\langle \Psi, \Phi \rangle = \Psi_{\alpha} \Phi^\alpha$$

We can then add the wavefunction to the equation of motion via the Schwinger-Dyson equation

$$\left\langle \frac{\delta S}{\delta \phi(x)} \right\rangle = \Psi_\alpha \left(\frac{\delta S}{\delta \phi(x)}\right)^\alpha_\beta \Psi^\beta = 0$$
View article

What's going on at Black Mesa

The original Half-Life had rather low resolution whiteboards, making it difficult to fit much legible text on it. All the whiteboards at Black Mesa were the same texture

Fig 1. Whiteboards in Half Life 1

With the same formulas, which were mostly easily recognizable equations, such as $E = mc^2$ and the Newtonian gravity

$$F_g = G\frac{Mm}{r^2}$$

For the fanmade recreation of Half-Life 1 in the Source engine, Black Mesa, the whiteboards were spruced up by making them higher resolution, with relevant content and unique for each room, although still fairly unrealistic in how well-written and small the text is on a laboratory whiteboard. A good compilation of all the whiteboards is available here. I'll look into the content of the theoretical physics-related ones, as I'm not particularly good with chemistry, crystallography or alien biology.

The first few whiteboards are about wormholes and faster-than-light travel, as befits Black Mesa's research on teleportation (although they appear in Sector C, which is more material science, not in Sector F, the Lambda Complex, where the actual teleportation research happens). The first one is on wormholes :

Black Mesa whiteboard 1

Let's break this down a bit. The first important line is this

$$ds^2 = -e^{2\Phi} dt^2 + \frac{1}{1 + \frac{b}{r}} dr^2 + r^2 \left[ d\theta^2 + \sin^2 \theta d\varphi^2 \right]$$

This is the metric tensor of the Morris-Thorne wormhole, expressed in the Schwarzdchild coordinates. From then on, the rest of the whiteboard roughly follows the same procedure as Morris and Thorne's original paper, Wormholes in spacetime and their use for interstellar travel: A tool for teaching general relativity. To the right are the properties of the functions making up the metric tensor.

Defining $\Phi$ as the redshift function, and putting the throat of the wormhole as the minimum value [of the shape function $b(r)$]. A drawing of the embedding of the surface defined by $(t= \text{const}, r, \theta = \pi/2, \varphi)$ is done, embedded in the three-dimensional space with cylindrical coordinates $(r, \varphi, z)$.

The proper radial distance $l$ is then defined with

$$l(r) = \pm \int_{r_0}^r \frac{dr'}{\sqrt{1 - \frac{b(r)}{r}}}$$

It's a bit hard to make out exactly what they ask but from context it seems to be the usual requirement that $dr/dl = 0$ at the throat, or equivalently $dl/dr \to \infty$ [this implies that at the throat ($r = r_0$), we have $b(r_0) = r_0$].

Once the metric is dealt with, they define a tetrad (a set of orthonormal basis vectors)

\begin{eqnarray} e_{\hat{t}} &=& e^{-\Phi} e_t\\ e_{\hat{r}} &=& (1 - \frac{b}{r}) e_r\\ e_{\hat{\theta}} &=& \frac{1}{r} e_\theta\\ e_{\hat{\varphi}} &=& \frac{1}{r\sin \theta} e_\varphi \end{eqnarray}

Then they compute the Einstein tensor along that tetrad

\begin{eqnarray} G_{\hat{t}\hat{t}} &=& \frac{b'}{r^2}\\ G_{\hat{r}\hat{r}} &=& - \frac{b}{r^3} + 2 (1 - \frac{b}{r}) \frac{\Phi'}{r}\\ G_{\hat{\theta}\hat{\theta}} &=& G_{\hat{\varphi}\hat{\varphi}} = \left[ 1 - \frac{b}{r} \right] \left[ \Phi'' + \Phi' (\Phi' + \frac{1}{r}) \right] - \frac{1}{2 r^2} \left[ b'r - b\right] \left[ \Phi' + \frac{1}{r} \right] \end{eqnarray}

And from there, using the Einstein field equations $G_{\hat{\alpha}\hat{\beta}} = 8\pi G T_{\hat{\alpha}\hat{\beta}}$, they compute the stress-energy tensor, and associated to it, the energy density $\rho$, the radial tension $\tau$ and the tangential pressure $p$.

\begin{eqnarray} T_{\hat{t}\hat{t}} &=& \rho = \frac{b'}{8 \pi G r^2}\\ - T_{\hat{r}\hat{r}} &=& \tau =\frac{1}{8 \pi G} \left[ - \frac{b}{r^3} + 2 (1 - \frac{b}{r}) \frac{\Phi'}{r}\right]\\ T_{\hat{\theta}\hat{\theta}} &=& T_{\hat{\varphi}\hat{\varphi}} = p = \frac{1}{8\pi G} \left\{ \left[ 1 - \frac{b}{r} \right] \left[ \Phi'' + \Phi' (\Phi' + \frac{1}{r}) \right] - \frac{1}{2 r^2} \left[ b'r - b\right] \left[ \Phi' + \frac{1}{r} \right] \right \} \end{eqnarray}

They then compute the quanity $\rho - \tau$, which is the product of the stress-energy tensor with a null curve of tangent $k^{\hat{\mu}} = (1, 1, 0, 0)$ in the tetrad basis, that is

$$T_{\hat{\alpha}\hat{\beta}} k^{\hat{\alpha}} k^{\hat{\beta}} = \rho - \tau$$

$ 8 \pi G(\rho - \tau)$ corresponding to the Einstein tensor $G_{\hat{\alpha}\hat{\beta}} k^{\hat{\alpha}} k^{\hat{\beta}}$. This quantity is then found to be equal to

$$ 8 \pi G(\rho - \tau) = - \frac{e^{2\Phi}}{r} \left[ e^{-2\Phi} (1 - \frac{b}{r}) \right]' $$

Since $b(r)$ has a minima at the throat of the wormhole ($r = r_0$) where $b(r_0) = r_0$, and any point beyond the throat will have $b(r) < r$, this leads to

$$e^{-2\Phi} (1 - \frac{b}{r}) > 0$$

Meaning that there is at least one point superior to $r_0$ where the derivative is positive, in particular this will be true at the throat itself, hence

$$ 8 \pi G(\rho(r_0) - \tau(r_0)) \leq 0$$

which means that the throat requires exotic matter, as the comment indicates. The details of why this is so is in the following section

This starts first with writing the stress-energy tensor in matrix form. Since the metric leads to a Type I stress energy tensor, it is of the form

$$T_{\hat{\alpha}\hat{\beta}} = \begin{bmatrix} \rho &0&0&0 \\ 0 & p_1 & 0 & 0 \\ 0 & 0 & p_2 & 0 \\ 0 & 0 & 0 & p_3 \end{bmatrix}$$

Then the Null Energy Condition (or NEC) is written down. The NEC states that given any null vector, the product of the stress-energy tensor with those vectors is positive or null,

$$T_{\alpha \beta} k^{\alpha} k^{\beta} \geq 0$$

For Type I stress-energy tensors, this will correspond to the following condition

$$\forall j, \rho + p_j \geq 0$$

that is, the sum of the energy density and momentum density in any direction is positive or zero, which is precisely the condition violated by the metric, for the radial momentum. The next equations concern other energy conditions, starting with the Weak Energy Condition (or WEC), which states that for any timelike vector $V$, we have

$$T_{\alpha \beta} V^{\alpha} V^{\beta} \geq 0$$

meaning that, for a Type 1 stress-energy tensor, we get

$$\rho \geq 0 \wedge \forall j, \rho + p_j \geq 0$$

The energy density is positive or zero, as well as its sum with the momentum density in any direction, which is also rewritten as $p_i \in [-\rho, \rho]$ later.

Then there is the Strong Energy Condition (or SEC)

$$(T_{\alpha \beta} + \frac{T}{2} g_{\alpha \beta}) V^{\alpha} V^{\beta} \geq 0$$

with $T = T_{\alpha \beta} g^{\alpha \beta}$, which can once again be rewritten as

$$\forall j, \rho + p_j \geq 0 \wedge \rho + \sum_j p_j \geq 0$$

The next three lines are the averaged versions of those energy conditions, the averaged null energy condition (ANEC), averaged weak energy condition (AWEC) and averaged strong energy condition (ASEC). For any curve $\gamma$ with tangent vector $k$ or $V$,

\begin{eqnarray} \int T_{\alpha \beta} k^\alpha k^\beta d\lambda &\geq& 0\\ \int T_{\alpha \beta} V^\alpha V^\beta d\tau &\geq& 0\\ \int (T_{\alpha \beta} k^\alpha k^\beta + \frac{1}{2} T) d\tau&\geq& 0 \end{eqnarray}

For a tangent vector $k^\alpha = \xi (1; \cos \psi_i)$ and $V^\alpha = \gamma (1; \beta \cos \varphi_i)$, this gives the relations \begin{eqnarray} \int (\rho + \sum_j \cos^2 \psi_j p_j) \xi^2 d\lambda &\geq& 0\\ \int (\rho + \beta^2 \sum_j \cos^2 \varphi_j p_j) \gamma^2 &\geq& 0\\ \int \left[ \gamma^2 (\rho + \beta^2 \sum_j \cos^2 \varphi_j p_j) \frac{1}{2} \rho + \frac{1}{2} \sum_j p_j \right] d\tau&\geq& 0 \end{eqnarray}

The follow-up isn't too clear, outside of the bit of text saying that violation of the WEC, SEC and DEC implies a violation of the ANEC (which is true).

Another whiteboard contains a few different things concerning faster-than-light travel.

First there is this metric tensor

$$ds^2 = -c^2 dt^2 + dl^2 + (k^2 + l^2) (d\theta^2 + \sin^2 \theta d\varphi^2)$$

This is actually the same metric as previously, except that it is in the proper distance coordinate system, using the relation $l(r)$ seen above, giving the metric $$ds^2 = -c^2 dt^2 + dl^2 + r^2(l) (d\theta^2 + \sin^2 \theta d\varphi^2)$$

with the radius function $r = \sqrt{l^2 + k^2}$, corresponding to the Ellis-Bronnikov wormhole.

The diagram

isn't actually the corresponding spacetime, but it is the diagram of the Kranikov tube, showing the light cones at various points, $c^2 dt^2 > dx^2 + dy^2 + dz^2$ simply indicating timelike curves.

The final part of the whiteboard concerns the Casimir effect

While realistically not ideal for propping up wormholes or Krasnikov tubes, the Casimir effect is a simple example of violation of most energy conditions, which wormholes do require. It's a fairly standard treatment (once again seemingly taken from the wikipedia article on the topic), first defining the eigenfunctions of a scalar field between two conducting plates separated by a length $a$ along the $z$ axis

$$\psi_n(x,y,z;t) = e^{-i\omega_n t} e^{ik_x x + i k_y y} \sin(k_n z)$$

with the momentum in the $z$ direction quantized as $k_n = \pi n / a$, and correspondingly the dispersion relation for the frequency

$$\omega_n = c \sqrt{k_x^2 + k_y^2 + \frac{\pi^2 n^2}{a^2}}$$

The expectation value of the energy of that field being

$$\langle E \rangle = \frac{\hbar}{2} 2 \int \frac{dk_x dk_y}{(2\pi)^2} \sum_{i= 1}^\infty A \omega_n$$

If we define $q^2 = k_x^2 + k_y^2$ and switch to polar coordinates, the integral becomes

$$\langle E \rangle = \frac{\hbar}{2} 2 \int \frac{1}{(2\pi)} r dr \sum_{i= 1}^\infty A \omega_n$$

or

$$\langle E \rangle = \frac{\hbar}{2} 2 \int \frac{1}{(2\pi)} r dr \sum_{i= 1}^\infty A c \sqrt{q^2 + \frac{\pi^2 n^2}{a^2}}$$

This quantity is fairly divergent, just being an integral of an unbounded positive function over $\mathbb R$, as well as an infinite sum of a strictly increasing function, so we need to perform a regularization, in this case with the use of a zeta regulator. Instead of summing the quantity $\omega_n$, we'll sum $\omega_n |\omega_n|^{-s}$, with the limit $s \to 0$. This gives the appropriate forumla

$$\frac{\langle E(s) \rangle}{A} = \frac{\hbar c^{1 - s}}{4\pi^{2}} \sum_n \int_0^\infty 2\pi q dq |q^2 + \frac{\pi^2 n^2}{a^2}|^{\frac{1-s}{2}} $$

with the integral evaluating to

$$\frac{\langle E(s) \rangle}{A} = - \frac{\hbar c^{1 - s} \pi^{2 - s}}{2 a^{3-s}} \frac{1}{3-s} \sum_n |n|^{3-s}$$

Which is obviously divergent for $s = 0$. This corresponds to the limit of $s \to 0$ of the Riemann zeta function $\zeta(s - 3)$. While the Riemann zeta function isn't well-defined around that point, it is possible to analytically continue it in that region. We can then obtain

$$\frac{\langle E \rangle}{A} = \lim_{s \to 0} \frac{\langle E(s) \rangle}{A} = - \frac{\hbar c \pi^2}{6 a^3} \zeta(-3)$$

And the analytical continuation of $\zeta$ gives

$$\zeta(-3) = \frac{1}{120}$$

For a total energy of

$$\frac{\langle E \rangle}{A} = -\frac{\hbar c \pi^2}{720 a^3}$$

It is followed by a simple calculation of the Casimir force generated by moving the plates,

$$\frac{F_c}{A} = - \frac{d}{da} \frac{\langle E \rangle}{A} = -\frac{\hbar c \pi^2}{240a^4}$$

The following whiteboard is a calculation of evaporation time for Schwarzschild black holes via Hawking radiation (the derivation of which is possibly taken from the wikipedia article on the topic) :

Consider a black hole of mass $M$, corresponding to a total energy of $E = Mc^2$. The total power radiated will be equal to the change in its energy

$$P = -\frac{dE}{dt} = -c^2 \frac{dM}{dt}$$

The power radiated $P$ is known from the Hawking derivation of the radiation

$$P = \frac{\hbar c^6}{15360 \pi G^2 M^2}$$

of which we will split out the mass to obtain

$$P = \frac{K_{ev}}{M^2},\ K_{ev} = \frac{\hbar c^6}{15360 \pi G^2}$$

Then we can simply solve the evaporation time as a simple ordinary differential equation of the mass, via

$$-c^2 \frac{dM}{dt} = \frac{K_{ev}}{M^2}$$

meaning that the energy radiated away is lost to the total mass. Using the common physicist trick of treating differentials as real objects, we get

$$M^2 dM= -\frac{K_{ev}}{c^2} dt$$

The black hole starts at time $t = 0$ at mass $M = M_0$ and finishes at time $t = t_{ev}$ at mass $M = 0$, with the solution

$$\int_{M_0}^0 M^2 dM= -\frac{K_{ev}}{c^2} \int_{0}^{t_ev} dt \Rightarrow - \frac{M_0^3}{3} = -\frac{K_{ev}}{c^2} t_{ev}$$

which gives us the evaporation time as

\begin{eqnarray} t_{ev} &=& \frac{c^2 M_0^3}{3 K_{ev}} = (\frac{c^2 M_0^3}{3} ) (\frac{15360 \pi G^2}{\hbar c^6})\\ &=& \frac{5120 \pi G^2 M_0^3}{\hbar c^4} \end{eqnarray}

The rest of the calculations are just the evaporation times for two masses : the Planck mass and one solar mass. For the Planck mass

$$M_P = \sqrt{\frac{\hbar c}{G}}$$

the corresponding evaporation time is

$$t_{ev} = \frac{5120 \pi G^2 M_p^3}{\hbar c^4} = 5120 \pi t_p = 5120 \pi \sqrt{\hbar G}{c^5} \approx 8.671 \times 10^{-40} s$$

while for one solar mass,

$$t_{ev} = \frac{5120 \pi G^2 M_\odot^3}{\hbar c^4} = 6.617 \times 10^{74} s$$

While the first part is some sort of mix of chemical and alchemical notation, the second is related to radioactive decay. The diagram is the Feynman diagram of a beta netative ($\beta^-$) decay, as seen from the perspective of quantum chromodynamics

with a neutron (composed of three quarks $udd$) emitting a $W^-$ boson, via the process $d \to u + W^-$, turning it into a proton ($uud$). The $W^-$ boson then decays into an electron and anti-electronic neutrino, $W^- \to e^- + \bar \nu_e$.

The text itself is a standard derivation of the eponymous half life formula for radioactive decay.

$$t_{1/2} = \frac{\ln(2)}{\lambda} = \tau \ln(2)$$ $$N_t = N_0 (\frac{1}{2})^{\frac{t}{t_{\frac{1}{2}}}}$$
View article