This is a blog article in progress, written by Jacob Biamonte. See also the quantum network theory series. To see discussions of the article as it was being written, visit the Azimuth Forum. For the final polished version, go to the Azimuth Blog.
NOTE: Here’s Part I
If you want to write your own article, please read the directions on How to blog.
Last time we talked about how defects in diamonds exhibit controllable quantum degrees of freedom. We showed that the energy levels of an elementary model of such a system resulted in a cubic equation. Today we are going to use a few oldtime tricks to visualize these energy levels.
Cubic equations embrace a long and twisted tale stretching back to the early days of mathematics. They lead to the eventual discovery of complex numbers, and even group theory.
Sometimes on Babylonian cuneiform tablets from the 20th to 16th centuries BC, people have noticed tables for calculating cubes and cube roots. The Babylonians could have used the tables to solve cubic equations, no one can provide evidence that they did.
The tale started to get more lively in the 16th century when del Ferro found a method for solving a class of cubic equations, namely those of the form $x^3 + mx = n$.
As we now know, all cubic equations can be reduced to this form if we allow m and n to be negative, but negative numbers were not known to del Ferro at that time!
Del Ferro kept his result to himself, until just before his death. He showed it to his student Fiore.
In 1530, Tartaglia entered the game. In fact, he was bold enough to let everyone know that he was able to solve two problems in cubic equations.
At that time, it seems the thing to do was to offer a challenge, and this is exactly what Fiore did. Each contestant put up money and proposed a number of problems for his rival. The winner would solve the most problems within 30 days. Sounds pretty reasonable.
Tartaglia received questions in the form x3 + mx = n, for which he had discovered a general method. Fiore questions came as x3 + mx2 = n,
Initially the system has a degenerate energy space configuration. This is called the zero-field splitting. In this case, the Hamiltonian is simply
With a static external magnetic field applied at angle $\theta$ with respect to the NV-symmetry axis and with strength $\beta$, one then arrives at the Hamiltonian together with a so called Zeeman splitting.
When probing the optical properties of the NV-center by driving the system with microwave pulses and detecting optical emissions, one can determine two resonant peaks.
We can consider some relationships of the eigenvalues of the Hamiltonian. The first is that
We know that the first transition occurs at the energy gap between the first excited energy state $\lambda_1$, and the ground state $\lambda_0$. In other words we know that the magic frequency is
The system is initially in the low-energy state corresponding to $\lambda_0$. So the transitions occur with respect to this state. The second transition is then
Combining this information into one equation gives us a value for $\lambda_0$ provided we know the zero-field splitting $\Delta$ and the transition frequencies from the ground state into the first and also second excited states.
We then consider the characteristic equation of the Hamiltonian, written in variable $\lambda$
There exists eigenvalues (roots) $\lambda_0$, $\lambda_1$, $\lambda_2$ such that the above polynomial can be expressed as
The coefficients of indeterminate $x$ are of interest to us (as opposed to e.g. $x^2$ or $x^0$). By the way, the expression of these terms into traces as above, are a special case of the (Newton–Girard identies)[], that’s how old this sort of trick is. If you really wanted to have fun, you might notice that they look just list a dispersion and write them like this …
The point is that by equating coefficients from the equation
we arrive at an expression for the magnitude of the magnetic field $\beta$ in terms of the measured transition frequencies $f_1$ and $f_2$ as
In other words, we can can use this NV to probe the external magnetic field! So it’s really a sensor, as we stated from the get go.
We then shift to the so called, center of mass of the polynomial. This can be thought of as adding a term proportional to the identity to make the Hamiltonian traceless. This amounts to the assignment
for $i = 0, 1, 2$ and from this we arrive at the depressed cubic
with
We’ve already solved for $\lambda_0$ in terms of the resonant frequencies $f_1$ and $f_2$. Changing this into an equation for $\kappa_0$ we arrive at
which we know is a root of $\kappa^3 + a \kappa + b = 0$. In other words, there exists $q, r, p$ such that
The cubic has a neat feature, if the roots are real, they are known to be related to the triangle as follows
Let’s first understand that, and then equate the quantities we just determine with the length of the sides, and the orientation (that is the rotation) of the triangle.
So there you have it. NV-centers in a nut shell, or at least in a rather long blog article. I certainly find this stuff interesting. As I mentioned at the start, the content here was inspired mainly from
There are a few more related references, including
For those that want to read another quick introduction, with more of an experimental focus, maybe check out
As I said, we’ve been working with Florian Dolde and others. So we have two of our own papers on NV-centers:
Quantum Simulation of Helium Hydride in a Solid-State Spin Register, Ya Wang, Florian Dolde, Jacob Biamonte, Ryan Babbush, Ville Bergholm, Sen Yang, Ingmar Jakobi, Philipp Neumann, Alán Aspuru-Guzik, James D. Whitfield, Jörg Wrachtrup, ArXiv preprints arXiv:1405.2696 (2014).
High fidelity spin entanglement using optimal control, Florian Dolde, Ville Bergholm, Ya Wang, Ingmar Jakobi, Sebastien Pezzagna, Jan Meijer, Philipp Neumann, T. Schulte-Herbrueggen, Jacob Biamonte, Jörg Wrachtrup, Nature Communications 5, 3371 (2014).
The sort of things I hear from experts in NV, about what will happen next, are mainly about scaling up the numbers of interacting NV-centers. Right now, the record is an impressive, but small handful of interacting particles (like 4 or 5), but with new the new methods being developed to create samples, people tell me this will likely double in the next two year. This would certainly open the door up to a lot more interesting applications. For example, say they are able to get 5 NV-centers to interact, this means we would then need to consider 59,049 energy levels. Modelling such a system using existing computing technologies will be challenging. Turning this on its head: such a system offers a computational resource, to at least, model itself.
Among other things, quantum physics explains how structures, such as crystals, form. In principle, using quantum physics, we could understand many more things than we actually do about the nature of mater.
The problem is that these equations seem impossible to solve using any computer not governed by quantum physics! The sad fact is that quantum computers which exist today as far to small to do anything practical.
One of the fun aspects of this is that we’re actually going to give a geometric interpretation to some things my friends Florian and Ya are actually measuring in their Stuttgart lab.