

Feature Articles: Adapting to the Changing Present and Creating a Sustainable Future Vol. 20, No. 10, pp. 67–76, Oct. 2022. https://doi.org/10.53829/ntr202210fa10 Number Theory and Quantum Physics Based on Symmetry—Themes from Quantum OpticsAbstractThe connection between mathematics and physics has produced many successes over the centuries. Modern mathematics is used, for example, in cosmology, high energy physics, and condensed matter theory. However, the connection between the theoretical study of quantum optics and modern mathematics, which is also necessary for the construction of quantum computers, is not always deep. In this article, I introduce the quantum Rabi model, which is the most fundamental theoretical model in quantum optics, and its relationship with modern mathematics, especially number theory, from the viewpoint of symmetry. I discuss the motivation for my research, how it has progressed, and future plans. I also introduce that mathematical conjectures arising from this study. Keywords: noncommutative harmonic oscillator, quantum Rabi model, zeta functions 1. Mathematics and mathematical researchThe purpose of mathematical research is to enrich mathematics by solving unsolved problems and discovering new problems, as well as by developing new mathematical visions and theories that can be applied to solutions of many problems. Mathematical research also gives rise to applications that are not immediately apparent. For instance, Riemannian geometry, developed in tandem with physicists, is largely independent from applications and is pure mathematics. It is now used, via Einstein’s theory of general relativity, in making GPS (Global Positioning System) highly effective. Using the example of “research on quantum models and number theory from the viewpoint of symmetry” I describe my research motivation, objectives and results, and future goals based on them. The models mentioned in this article are the noncommutative harmonic oscillator (NCHO) [1], quantum Rabi model (QRM), and asymmetric QRM (AQRM) [2, 3]. Mathematics is, in part, an endeavor to understand, in a unified manner, facts that appear to be far apart, which sometimes takes quite long time. For example, there is the application of the discovery in ancient Greece of the infinity of prime numbers and uniqueness of prime factorization to today’s public key cryptography. It took more than 2500 years to understand that mathematical discoveries can be used to maintain secrecy. Also, the structure of the universe is gradually being revealed through nonEuclidean geometry and group theory, which arose out of the human mind. We now know that economic phenomena and the movement of pollen can be understood in the same way since they are due to the mathematical description of Brownian motion, i.e., the development of stochastic analysis. In fact, there is a special value in the clarification of unquestionable affinities between subjects that have long been believed to be unrelated to each other. According to Henri Poincaré, mathematics is the art of identification. As ShinIchiro Tomonaga^{*1} wrote: “What has worked well is experiment, demonstration/proof, and abstraction through mathematics. Abstraction makes it very universal. These are the great strengths of physics.” Mathematical theorems, however, are not invented but discovered. This becomes clear from Kunihiko Kodaira’s view^{*2} that mathematical research is, using paper and pencil, to dig out mathematical truths that do exist independent from human minds but are hidden. The symmetry in the title refers to the actions of groups. We may recall the congruence transformation groups of Euclidean planes and spaces and the crystallographic groups familiar from physics and chemistry. Group theory was founded by Évariste Galois, who died in a duel at the age of 20. Before Galois, there were formulas for roots of quadratic, cubic, and quartic equations, and there were efforts to find those for quintic equations. However, what Galois did was to consider the entirety of the substitutions of the roots of an equation (the permutation group) and found from its structure that there is in general no root formula that uses the usual algebraic operations (addition, subtraction, multiplication, division) and application of radicals (square roots, cube roots, etc.) to coefficients of a polynomial when the equation is of fifth degree or higher^{*3}. This is a major paradigm shift. What is the structure? It is related to the strength of the noncommutativity of the group. There is a large difference between putting on a suit after putting on a shirt and putting on a shirt after putting on a suit. I hope this gives an idea of noncommutative operations (actions).
2. NCHOQuantum harmonic oscillators are fundamental in quantum theory. Their energies (eigenvalues) are given as halfintegers in the standard normalization, and the corresponding eigenstates are Hermitian functions, i.e., given essentially by the Hermite polynomials times . A clean description of this can be obtained from the representation theory^{*4} of the threedimensional Lie algebra consisting of 2 × 2 matrices with trace 0. The captures the infinitesimal action of the Lie group SL_{2}(R) formed by all real 2 × 2 matrices with determinant 1. Lie group theory originated from the work of Sophus Lie, who was impressed by the Galois theory of algebraic equations and wanted to construct a sort of Galois theory for algebraic differential equations. Lie groups are groups as the name implies and geometric objects called manifolds. In the geometric context, is considered as the tangent space of SL_{2}(R) at the identity element of the group. The Hermite functions form a basis of the Hilbert space L^{2}(R) of square integrable functions on the real line. Many special functions, such as Hermite functions, Bessel functions, and Jacobi polynomials, are obtained by specializing one of the parameters a, b, or c of the Gaussian hypergeometric function F(a, b, c; x) to integer values, etc. The differential equation (Gaussian ordinary differential equation (ODE)) satisfied by F(a, b, c; x) has three regular singular points (e.g., 0, 1, ∞). The family of Gaussian hypergeometric functions is indispensable to the number theory of elliptic curves, which has been the flower of mathematics since the 19th century. In representation theory, however, special functions derived from Gaussian hypergeometric functions and its multivariable versions can be understood essentially as matrix elements of representations of Lie groups. I will now turn our attention to the following. (1) The eigenvalues of (quantum) harmonic oscillators H = a^{†}a + , where a and a^{†} are respectively the creation and annihilation operators subjecting [a^{†}, a] = 1, are halfintegers. Therefore, the spectral zeta function^{*5} of the harmonic oscillator determined from the eigenvalues essentially coincides with the Riemann zeta function ζ(s)^{*6}. (2) The symmetry of SL_{2}(R) remains in the world of Gaussian hypergeometric functions. For instance, the matrix elements of representations of SL_{2}(R) and, almost equivalently, the solution of the SL_{2}(R)invariant differential equation all can be expressed essentially by Gaussian hypergeometric functions. The group SL_{2}(R) is behind (1) and (2) above. To explore the wider world, we need to weaken the symmetry. The only noncommutative property of harmonic oscillators is the canonical commutation relation (CCR) controlled by ^{*7}. The NCHO^{*8}, which defines a system of ODEs, with two parameters α and β, adds matrix noncommutativity. The name comes from this noncommutativity. Of course, it would not be rational to define it in an arbitrary manner. From expected physical applications and corresponding mathematically reasonable consideration, we thought that it should be at least required that the Hamiltonian of the NCHO have only real discrete eigenvalues and that it be kept a sufficiently weak symmetry to allow the harmonic oscillator to be recovered when α = β. The former requires that be a positive Hermitian operator (i.e., selfadjoint operator). The latter is that the spectral zeta function ζ(s) of be a good generalization of ζ(s) though neither a functional equation nor Euler product may be expected [4]. The convergence region of the defining series of ζ(s) is only in the halfplane, while ζ(s) has analytic continuation to the whole plane. The Riemann Hypothesis^{*9} asserts that the real part of an imaginary number s with ζ(s) = 0 is always . There are two analogues of the Riemann Hypothesis at the top (and among the bestknown achievements) of 20th century mathematics. One is the Weil Conjecture for the congruent zeta function and the other is the Selberg zeta function for the prime geodesic distribution. Both conjectures were settled by the fact that each zeta function has a determinant expression in terms of suitable Hermitian operators^{*10}. The study of a possible determinantal expression of ζ(s) is also important in the attack on the original Riemann Hypothesis.
3. QRM and NCHOIsidor Isaac Rabi described the interaction of classical light with atomic quantum states in 1936. In 1963, Edwin Jaynes and Fred Cummings quantized light as well and defined the QRM and its rotating wave approximation, now known as the JaynesCummings model. However, while the integrability of the QRM was unknown, the JaynesCummings model was useful because of the existence of continuous invariants, making the theoretical treatment easy, and suitability to experiments if the lightmatter coupling is very weak, roughly a fraction 10^{−6} of the mode frequency^{*11}. However, with the possibility of developing a quantum computer on circuit quantum electrodynamics (QED) platforms, the weak coupling approximation is no longer justified, and one must use the QRM in the ultrastrong coupling regime. In 2011, Daniel Braak [2] demonstrated the integrability of the QRM, and research on the periphery of the QRM has advanced rapidly since. In fact, the experimental measurement results [5] using superconducting artificial atoms were consistent with the theory of the QRM and its generalization, the AQRM (see also [6]). The NCHO has expanded the world from Gauss’s as expected. In fact, the eigenvalue problem of the NCHO was found to be equivalent to the existence of holomorphic solution of the Heun ODE^{*12} in the complex region containing 0, 1, and not including αβ. This fact was discovered by Hiroyuki Ochiai [7] and fully understood 15 years later [8]. The (A)QRM also has a Heun picture [2]. However, unlike the NCHO, it is a confluent Heun differential equation. Nevertheless, by merging the two singularities as αβ → ∞, a confluent Heun picture of the QRM can be obtained from the Heun picture of the NCHO [8]. In other words, the NCHO can be considered as a covering model of the QRM^{*13} (Fig. 1, the NCHO is a covering model of the QRM).
Quantization is discretization and passing from a commutative world to a noncommutative world. The former is the original arena of number theory. In fact, the special values of ζ(s) in the NCHO shows rich number theory as I explain below. For instance, as to the special values of the spectral zeta function of the harmonic oscillator, essentially given by ζ(s), there is the old Basel problem to determine the value of the infinite series ^{*14}, which was solved by Euler. The even integer values of ζ(s) are generally expressed in terms of powers of π and rational numbers. However, the value and properties of ζ(s) at odd points remained unknown for more than 300 years until Roger Apéry’s proof that ζ(3) is irrational in 1978. The best result on irrationality is still the theorem that at least one of ζ(5), ζ(7), ζ(9), or ζ(11) is an irrational number. However, the research on special values is a deep problem in modern mathematics, such as the geometric conjecture on special values of various zeta functions by Pierre Deligne^{*15} who solved the Weil Conjecture using Grothendieck’s arithmetic geometry, and the BirchSwinnertonDyer conjecture on special values at the central point of Lfunctions (relatives of zeta) defined from elliptic curves and its generalization by Bloch, Beilinson, and Kato. Apéry used a sequence of rational numbers, now called the Apéry numbers, to prove the impossibility that ζ(3) is rational. Frits Beukers took up the challenge of finding out why and further proving irrationality and developed a brilliant study of elliptic curves, K3 surfaces, and automorphic forms behind the proof of the irrationality [9]. As a result, he showed remarkable congruence relations of the Apéry numbers, which have become the subject of research in algebraic number theory. Surprisingly, the special values ζ(2) and ζ(3) for the NCHO also yield sequences of numbers with properties like those of the Apéry numbers [10, 11, 12]. Elliptic curves and automorphic forms underlie these properties, and when s = 4, even automorphic integrals^{*16}, which generalize the notion of automorphic forms, appear [10] via certain nearly holomorphic automorphic forms initially studied by Goro Shimura^{*17}. The next development will be the study of a wider range/class of automorphic forms (also related to the Langlands program, which is extremely important in mathematics [13]).
4. Spectra of the AQRM and hidden symmetriesThe importance of the partition function (weighted sum of states) in statistical mechanics is that it encodes the time evolution of the system. Knowing the exact behavior of a quantum system during a long period is of utmost importance to any task in quantum computing. Purely numerical approaches are severely limited due to the accumulation of rounding errors in the longtime limit. Therefore, it is necessary to derive analytical formulas for its heat kernel and partition function. In fact, the partition function and spectral zeta function have a onetoone correspondence through the Mellin integral transform. However, obtaining these analytical formulas was considered extremely difficult. Physicists use a highly intuitive form of expression called Feynman path integral to carry out various asymptotic or approximate (as well as precise) calculations to gain a better understanding of the physical system of interest. The translation invariant measure and integral on the space of paths cannot be treated mathematically in a rigorous manner, except in some special cases. Of course, various studies are still ongoing to understand it. I, together with ReyesBustos, have attempted to derive a heat kernel for the (A)QRM^{*18}. The method, which has been successfully applied [14], is based on the action of the symmetric group of the mth order on (m = 1,2,...), the Fourier transform on ^{*19}, and graph theory, where Z_{2} is Z/2Z or the field of two elements F_{2} = {0,1}. The resulting analytic formula is practically given by the series over nonnegative integers m for integrals over msimplex (an mfold integral) of elementary functions given by products and quotients of exponential functions. It can also be interpreted as the infinite sum (integral) of the orbit integrals of over the entire equivalence class by the infinite symmetric group of discrete paths connecting 0 and ∞ (Fig. 2, Discrete paths and infinite symmetry groups (Geometric Interpretation)). From a representation theoretic point of view, the action of on corresponds to a decomposition into the direct sum of irreducible representations (Fig. 3 (Algebraic Interpretation)). The formulas obtained are similar to formulas for partition functions obtained for the spinboson model or the Kondo model (describing magnetic impurities in nonmagnetic host metals, named after Jun Kondo, 1930–2022), a formally exact description of the matrix element of the heat kernel, and to approximate partition functions for various physical systems. For a general model, the (possible) existence of an infinitedimensional algebraic system acting on the noncountable paths of the Feynman integrals and the acquisition of discrete paths as their equivalence analogue should also be pursued (Fig. 4, Conjecture on the paths of the Feynman path integral).
From the analytical formula for the heat kernel, the partition function is obtained, and the contour integral expression of the spectral zeta function is gained^{*20}. Thus, the analytic connection to the whole plane is obtained, the asymptotic evaluation of the number of eigenvalues (Weyl’s law) follows, and we can start the special value study [15]. What is the number theory and geometric structure of the QRMs? This is one of the most important questions in this area. The AQRM is of major importance in circuit QED implementations of the QRM, thus for applications to quantum computing. Unlike the QRM, it has no visible Z_{2}symmetry (parity), if the bias is not zero. Nevertheless, it shows characteristic spectral degeneracies at halfinteger values of the bias parameter. In fact, the AQRM is degenerate if and only if the bias parameter (of the bias operator) of the Hamiltonian^{*21} is a halfinteger ℓ/2 (ℓ ∈ Z) [3, 16]. In physics, however, degeneracy is considered to indicate the existence of a symmetry [17]. Unlike classical mechanics, the concept of integrability is not clear in quantum mechanics [18]. However, the analogy with classical systems has led to the search for operators that commute with H_{AQRM} but are not functions of H_{AQRM}, indicating the presence of an additional conserved quantity [19]. Inspired by this study, we [20] conducted a thorough analysis to completely determine the operators that are commuting with the H_{AQRM}. Interestingly, the ring^{*22} of operators commuting with H_{AQRM}, has a unique natural generator J_{ℓ} up to a constant multiple. Moreover, the square of J_{ℓ} is expressed as J_{ℓ}^{2} = p_{ℓ} (H_{AQRM}, g, Δ), a polynomial in H_{AQRM}, where p_{ℓ} is an ℓth degree polynomial. The equation y^{2} = p_{ℓ} (x; g, Δ) defines a hyperelliptic curve [21]. On the other hand, the degenerate eigenvalues of the AQRM are of the form N ± + g^{2} with positive integer N. Fixing N and ℓ, (g^{2}, Δ^{2}) is obtained as zeros of the constraint polynomial (g^{2}, Δ^{2}) and (g^{2}, Δ^{2}) [12]. In fact, (g^{2}, Δ^{2}) divides (g^{2}, Δ^{2}), and the quotient polynomial is (g^{2}, Δ^{2}) (>0) [16]. It should also be noted that the AQRM degeneracies are, geometrically, conical intersection points in the threedimensional spectral graph (Left of Fig. 5).
From experimental calculations, we know that If this holds in general, the expectation “hidden symmetry behind the degeneracy” becomes explainable concretely provided we may think of the existence of a nontrivial operator commuting with the Hamiltonian (which is in fact an involution modulo the Hamiltonian) as evidence of quantum integrability^{*23} (Fig. 5 and Fig. 6). There is, however, no way to prove it, and it remains a conjecture [21]. However, this conjecture is very similar in form to a certain theorem^{*24} in the Diophantine geometry^{*25} [22] in its form. Therefore, a clue to the proof might be found there. We also expect that the AQRM has a finite number of degenerate eigenvalues for fixed parameters g and Δ. Although it is still unresolved, we can find some formal similarities with the Mordell Conjecture^{*26} in this case.
5. ConclusionLet us discuss a few phrases that characterize mathematics. First, “All things are numbers” by Pythagoras. This is reminiscent of the big data and artificial intelligence (AI) that dominate modern society. The genius Leonardo da Vinci said, “Engineering is the paradise of the mathematical sciences. After all, it is here that the fruits of mathematics are borne.” Galileo Galilei’s “The universe is written in mathematical language” is well known. Wigner’s “The unreasonable effectiveness of mathematics in the natural sciences,”^{*27} is a brilliant point. Although mathematical research is like musical composition in some respects, it has been difficult to find a mechanism to make mathematical papers play like music, such as reproducing a symphony from a score. However, the use of mathematics has emerged because of machine learning, AI, and other technologies. Although the motivations for research are diverse, including pure mathematics driven solely by curiosity, it now seems as if an opportunity has arrived for basic mathematics and applied mathematics to play a brilliant concerto. I close this article with the hope that a new “mathematical music” (in Japanese, we may express it as (sugaku) meaning enjoining mathematics and is a homonym with (sugaku) meaning the study of numbers, i.e., mathematics)^{*28}.
References
