Extremal Determinants of Laplace-Beltrami Operators for Rectangular Tori

In this work we study the determinant of the Laplace-Beltrami operator on rectangular tori of unit area. We will see that the square torus gives the extremal determinant within this class of tori. The result is established by studying properties of the Dedekind eta function for special arguments and refined logarithmic convexity and concavity results of the classical Jacobi theta functions of one real variable are deeply involved.


Introduction
The search for extremal geometries is a popular topic in many branches of mathematics and mathematical physics. In this work, we pick up a result by Osgood, Phillips and Sarnak [26] on extremals of determinants of Laplace-Beltrami operators on tori and restrict the assumptions, excluding their solution of the following problem.
Among all 2-dimensional tori of area 1, which torus maximizes the determinant of the Laplace-Beltrami operator?
The answer in [26] is that the torus identified with the plane modulo a hexagonal (sometimes called triangular or equilateral) lattice gives the unique solution. However, if we only consider rectangular lattices, this solution is not possible and the natural assumption is that the square lattice will lead to the optimal solution. We will prove that this is indeed the case. Both problems, the one for general and the one for rectangular lattices, are closely related to the study of extremal values of the heat kernel on the torus [3], [5], finding extremal bounds of Gaussian Gabor frames of given density [17], [18], [19], as well as the study of certain theta functions [25].
It is worth noting that in all cases the extremal solutions are the same as for the classical sphere packing and covering problem in the plane. This immediately raises the question about extremal solutions for the above problems in higher dimensions. An interesting aspect is that in higher dimensions the optimal arrangements for the sphere packing and covering problem differ from each other. Hence, in order to guess what the right solution could be, one first has to decide whether one deals with a packing or a covering problem. However, we will not discuss higher dimensions in this work.
Another common theme is that we deal with theta functions in one way or another, which take a prominent role in several branches of mathematics. They appear in the studies on energy minimization [6], the study of Riemann's zeta and xi function [9], [12] or the theory of sphere packing and covering [11], including the recent breakthrough for sphere packings in dimension 8 by Viazovska [31] and in dimension 24 [10]. They also appear in the field of time-frequency analysis and the study of Gaussian Gabor frames [18], [19], [21] or the study of the heat kernel [30] to name just a few.
An open question concerns how the above problems can be linked to old, unsolved problems in geometric function theory, namely finding the exact values of Bloch's constant (1925) [7] and Landau's constant (1929) [24]. The correct solutions are conjectured to be given in the work of Ahlfors and Grunsky [1] and in the work of 1 Rademacher [27], respectively.
Baernstein repeatedly suggested that a better understanding of the behavior of extremal values of the heat kernel on the torus might lead to new insights for the mentioned constants [2], [3], [5] and, after some research, the author shares this opinion. In fact, this work is a result of the author's study on a conjecture of Baernstein, Eremenko, Frytnov and Solynin [4] related to Landau's constant, which Eremenko posed again as an open problem in an unpublished preprint in 2011 [13]. However, besides the fact that in both cases varying metrics on rectangular tori are involved, it is not clear to the author how deep the connection between this work and the mentioned conjecture in [4] and [13] truly is.
Let us return to the question posed at the beginning. In [26] we find the following result, which is 2 Corollary 1.3(b) in that work. For a lattice Λ, let ∆ Λ be the Laplace-Beltrami operator for the torus T Λ = C/Λ of unit area, then 2 )| 4 ≈ 0.35575 . . . with equality if and only if Λ is hexagonal. The result was established by first showing that the hexagonal lattice gives a local maximum by exploiting general facts about modular forms. Then, a numerical check gave the result that this local maximum is indeed global.
In another work, Sarnak [29] mentions that 3 Karnaukh has shown that the square lattice gives the only other critical point and that it is a saddle point. In particular, this implies that there exists a one-parameter family of lattices, within which the square lattice yields a local maximum of the determinant. Actually, we will see that among all rectangular tori, the square torus gives the global maximum. Our main result is as follows.
Theorem 1.1 (Main Result). For α ∈ R + , we denote the rectangular torus of unit area by T α = C/(α −1/2 Z × i α 1/2 Z) and the Laplace-Beltrami operator by ∆ α . Then with equality if and only if α = 1.  In [26] there is a typo concerning the spanning vectors of the extremal lattice, it should read Λ = 2 √ . 3 Peter Sarnak does not mention any reference for this claim. Anton Karnaukh was his PhD student around that time, but the author also could not find any result in that direction in Karnaukh's thesis [23].
Our proof of Theorem 1.1 will at first be parallel to the proof in [26], in particular, we will show that the problem about the determinant can be transferred to a problem of finding the maximum of the Dedekind eta function on a ray in the upper half plane. After that point, the proof will differ greatly from the methods in [26]. Not only will we show that the square lattice yields the global maximum, we will also give a precise behavior of the determinant as the lattice parameter varies. The key in our proof is to exploit the fact that the eta function can be decomposed into a product of Jacobi's theta functions. The result will follow from certain logarithmic convexity and concavity results, partially established already in [19].
This work is structured as follows; • In Section 2 we recall the definitions of Laplace-Beltrami operators on tori and their determinants as well as the results from the work of Osgood, Phillips and Sarnak [26]. Also, we will see how the determinant connects with the Dedekind eta function. • In Section 3 we define Jacobi's classical theta functions of one real variable and show how they can be used to express the Dedekind eta function. The proof of the main result will follow from refined logarithmic convexity and concavity statements related to Jacobi's theta functions as described by Faulhuber and Steinerberger [19].

The Laplace-Beltrami Operator on the Torus
In this section we recall the results established in [26] and how the problem of finding extremal surfaces for determinants of Laplace-Beltrami operators connects with the Dedekind eta function. Osgood, Phillips and Sarnak studied the determinant of the Laplace-Beltrami operator on a surface with varying metric as a function of the metric. In the case of the torus with flat metric, the varying of the metric can be interpreted as varying the lattice associated to the torus (and keeping the standard metric). Furthermore, we will introduce the heat kernel of a Laplace-Beltrami operator and its determinant.
We denote the Laplace-Beltrami operator on the 2-dimensional torus T Λ = C/Λ by ∆ Λ . The torus is identified with the fundamental domain of the lattice Λ which is a discrete subgroup of C. A lattice Λ is generated by integer linear combinations of two nonzero complex numbers z 1 and z 2 with the property that z 1 z 2 / ∈ R; The area of the torus is then defined to be the area of a fundamental domain, i.e., area(Λ) = |Im(z 1 z 2 )| = |x 1 y 2 − x 2 y 1 |, z k = x k + i y k , k = 1, 2.
In this work we solely deal with rectangular lattices, i.e., we can choose a basis of the lattice z 1 , z 2 Z with the property that the ratio of z 1 and z 2 is purely imaginary; We note that any rectangular lattice can be identified with αZ × i βZ where α, β ∈ R + . In this work, it will be no restriction to assume that the lattice has unit area, i.e., αβ = 1.
2.1. The Laplace-Beltrami Operator on Manifolds and its Heat Kernel. We will now introduce Laplace-Beltrami operators on connected Riemannian manifolds as well as the associated heat kernel and the determinant. For further reading on heat kernels we refer to [22]. Also, the procedure of introducing determinants of Laplace-Beltrami operators is described in [5], [26] or [28], and we will follow these references.
Just for the moment, let us change to a more general notation. Let M be a connected Riemannian manifold and ∆ M the corresponding Laplace-Beltrami operator. The heat semigroup is defined as The action of the group P t on a function f ∈ L 2 (M ) is (by abuse of notation) given by the integral operator where µ(y) is the 4 Lebesgue measure on M . The integral kernel p t is called the heat kernel of the Laplace-Beltrami operator. For any y ∈ M , the heat kernel fulfills the heat equation and for any y ∈ M p t ( . , y) → δ y , t → 0.
If the spectrum of ∆ M is discrete and consists of eigenvalues {λ k } ∞ k=1 with {φ k } ∞ k=1 being the sequence of corresponding eigenfunctions, constituting an orthonormal basis for L 2 (M ), then the heat kernel can be expanded as Also, in this case the trace of the heat kernel is given by The left-hand side of the above equation defines the trace of the heat kernel in general. To the Laplace-Beltrami operator one associates a zeta function in the following way This function can be continued analytically in C\{1}. Formally, the determinant of ∆ M is given by This product is not necessarily meaningful and the proper definition of the determinant is to use the zeta regularization . A closely related function is the height function of a Riemannian manifold [28], . 4 (M, µ) can also be a weighted manifold with µ being any measure with smooth positive density with respect to the Riemannian measure.
This function is an isospectral invariant and it is clear that problems about det ′ ∆ M can be transferred to problems about h M and vice versa.

2.2.
The Determinant on the Torus. We return to the case of the torus C/Λ where Λ is a lattice of unit area. As we will see, it is no restriction to assume that Λ = c 1, z Z , where c is chosen such that the area of the lattice is 1, i.e., c = Im(z) −1/2 . In other words, the problem under consideration is invariant under rotation and scaling, just like the classical sphere packing and covering problems. In this case, the eigenvalues of the Laplace-Beltrami operator ∆ Λ are (2π|z λ |) 2 , where z λ ∈ Λ. At this point, we mention that the eigenvalues are actually given by (2π|z ⊥ λ |) 2 , z ⊥ λ ∈ Λ ⊥ , the dual lattice. However, for 2-dimensional lattices, the relation between a lattice and its dual lattice is simply given by This means that the dual lattice is a 90 degrees rotated, scaled version of the original lattice. However, as mentioned we deal with lattices of unit area, hence the scaling factor is irrelevant, as is the rotation. Thus, in our concrete situation, there is no need to distinguish between the lattice Λ and its dual Λ ⊥ (this only results in a re-labeling of the eigenvalues). The zeta function is, consequently, given by where the prime indicates that the sum does not include the origin. Using the definition of the lattice Λ = z 1 , z 2 Z , we re-write the zeta function as For the last equality we set z = z 2 z 1 and y = Im(z) > 0, where the second condition is imposed by the fact that the lattice has unit area. Osgood, Phillips and Sarnak [26] now use the fact that the last series is a multiple of the Eistenstein series The final step in order to compute det ′ ∆ Λ is now to differentiate Z Λ with respect to s and evaluate at 0. We compute that Now, the following result is a consequence of Kronecker's limit formula.
The Dedekind eta function η(τ ) is defined in the upper half plane H by the infinite product Hence, It follows that the (zeta regularized) determinant of the Laplace-Beltrami operator ∆ Λ , with Λ = y −1/2 1, z Z , z = x + iy, y > 0, is given by

Rectangular Tori.
We have seen that maximizing the determinant det ′ ∆ Λ , where the lattice is given by Λ = y −1/2 1, z Z corresponds to maximizing y|η(z)| 4 . In the case of a rectangular torus C/(y −1/2 Z × i y 1/2 Z), the problem in focus is We could end this work now with the fact, already observed by Baernstein and Vinson [5], that maximizing det ′ ∆ Λ is implied by minimizing tr(e t∆ Λ ). Montgomery showed that the unique minimizer of tr(e t∆ Λ ) is the hexagonal lattice [25]. In a recent work on Gaussian Gabor frames, the results of Faulhuber and Steinerberger [19] imply that the square lattice is the unique minimizer of tr(e t∆ Λ ) within the class of rectangluar lattices. Hence, Theorem 1.1 is implied by the results in [19]. However, we will come up with a proof independent from the results in [19] and [25], but, however, the techniques are similar. The rest of this work is meant to give a deeper insight into the mentioned problems and maybe these insights can be valuable for some of the related problems mentioned in the introduction.

Jacobi's Theta Functions
In this section we study properties of Jacobi's theta functions which will lead to a deeper understanding of Theorem 1.1. We start by defining the theta function in accordance with the textbook of Stein and Shakarchi [30].
Definition 3.1. For z ∈ C and τ ∈ H (the upper half plane) we define the theta function as This function is an entire function with respect to z and holomorphic with respect to τ . As stated in [30], the function arises in many different fields of mathematics, such as the theory of elliptic functions, the theory of modular functions, as a fundamental solution of the heat equation on the torus as well as in the study of Riemann's zeta function. Also, it is used to prove results in combinatorics and number theory.
The function can also be expressed as an infinite product.
Proposition 3.2 (Jacobi triple product). For z ∈ C and τ ∈ H we have It also fulfills the following identity.
. For a proof of Theorem 3.3 and more details on theta functions of two complex variables as well as the product representation we refer to the textbook of Stein and Shakarchi [30].
Whereas derivatives of Θ with respect to z are often studied, it seems that studies of its derivatives with respect to τ are less common. The following Lemma contains a symmetry result for the logarithmic derivative of Θ with respect to τ . Lemma 3.4. For z ∈ C and τ ∈ H we have In particular, for z = 0 we get Proof. We start by taking the logarithm on both sides of the identity in Theorem 3.3 Differentiating with respect to τ on both sides and a multiplication by τ yields, after rearranging the terms, the desired result.
As a next step, we define Jacobi's theta functions of two variables which we will then restrict to certain domains. For z ∈ C and τ ∈ H we define We note that ϑ 3 (z, τ ) = Θ(z, τ ) and that any ϑ j can be expressed via Θ from equation (3.1). The functions which we will study in the rest of this work are restrictions of the above functions with (z, τ ) = (0, ix), x ∈ R + . In particular, all these functions are real-valued.
Definition 3.5. For x ∈ R + we define the real-valued theta functions in the following way.
It does not make much sense to study properties of θ 1 (x) if defined as above, because ϑ 1 (0, ix) = 0 for all x ∈ R + . However, we will see that θ 1 is involved in the proof of our main result in some sense.
All of the above functions can also be expressed by infinite products.
These representations can be quite useful when studying the logarithmic derivatives of these functions. We note that for z ∈ R and purely imaginary τ = ix, x ∈ R + the theta function Θ(z, ix) is maximal for z ∈ Z and minimal for z ∈ Z + 1 2 . These special cases correspond to Jacobi's θ 3 and θ 4 function θ 3 (x) = Θ (0, ix) , θ 4 (x) = Θ 1 2 , ix . For x ∈ R + , the θ 2 -function can be expressed via Θ in the following way.
A recurring theme will be the frequent use of the differential operator x d dx . We will use the notation x d dx n for its iterated repetition, i.e., In particular, we will study properties of the logarithmic derivative of a certain function f : R + → R + , on a logarithmic scale. Let us explain how this statement should be interpreted. By using the variable transformation y = log(x) we extend the domain from R + to R and consider the new function log (f (e y )). Therefore, we get an extra exponential factor each time we take a (logarithmic) derivative. We have d dy log (f (e y )) = e y f ′ (e y ) f (e y ) .
By reversing the transformation of variables, we come back to the original scale, but the factor x stays. Without being explicitly mentioned, these methods were used in [18], [19], [25] to establish uniqueness results about extremal theta functions on lattices. We will now provide a new point of view on as well as new results related to the mentioned articles.
For what follows it is necessary to clarify the notation of the Fourier transform and Poisson's summation formula which are given by respectively. Both formulas certainly hold for Schwartz functions and since we will apply both only on Gaussians we do not have to worry about more general properties for the formulas to hold. As a consequence of the Poisson summation formula we find the following, well-known identities x . We note the common theme in the last identities which leads to the following lemma.
Lemma 3.6. Let r ∈ R and suppose that f, g ∈ C 1 (R + , R) do not possess zeros. If f and g satisfy the (generalized Jacobi) identity Proof. Both, f and g are either positive or negative on R + since they are real-valued, continuous and do not contain zeros. As they also fulfill the identity x r f (x) = g 1 x both of them possess the same sign. Without loss of generality we may therefore assume that both functions are positive because otherwise we change the sign which does not affect the results. Therefore, we have log (x r f (x)) = log g 1 x .
Using the differential operator x d dx on both sides directly yields We remark that an alternative proof of Lemma 3.6 is given in [19] for the special case f = g = θ 3 and r = 1 2 (a generalization of the proof in [19] is possible with the arguments given there). With this proof, the result was derived. However, it seemed to have gone unnoticed that the result can be put into a more general context. With the previous lemma we also get the results which were not given in [19]. Also, the author just learned that the above formulas can also be found in the monograph of Borwein and Borwein [8, chap. 2.3]. All of the above results actually contain symmetry statements about logarithmic derivatives of Jacobi's theta functions on a logarithmic scale. We note that the identities (3.7) and (3.8) are special cases of Lemma 3.4 when Θ is restricted to certain rays in C × H.

3.1.
Properties of Theta-2 and Theta-4. We start with monotonicity properties of the logarithmic derivatives of Jacobi's θ 2 and θ 4 functions on a logarithmic scale. The following result was already given in [19] and the proof can be established by using the product representation of θ 4 .
is strictly decreasing on R + . As a consequence of Proposition 3.7 we derive the following result.
Proposition 3.8. The Jacobi theta function θ 4 is strictly logarithmically concave or, equivalently, Proof. The statement that x θ ′ 4 (x) θ 4 (x) is strictly decreasing on R + was already proved in [19]. We observe that lim x→∞ x θ 4 (x) > 0 and, therefore, < 0, and it follows that which proves that is strictly decreasing, hence, θ 4 is logarithmically concave. As a consequence of Proposition 3.7 and Lemma 3.6 we obtain the following property of θ 2 , which was already claimed, but not proven, in [19].
Proof. Proposition 3.7 tells us that x θ ′ 4 (x) θ 4 (x) is strictly decreasing on R + . By using Lemma 3.6 or the first identity in (3.8) we obtain Strict logarithmic convexity of θ 2 is now equivalent to We set .
We proceed with proving the result involving θ 2 . We already know that θ 4 is logarithmically concave, which is equivalent to the fact that θ 4 (x) is monotonically decreasing. By using identity (3.8) we get . Therefore, we have We note that a similar simple argument as used for the expression involving θ 2 does not work for θ 4 . Due to identity (3.8) it also seems plausible to assume that for a monomial weight x r , the second power, i.e. x 2 , is the limit for the monotonicity properties of the logarithmic derivative of θ 2 and θ 4 as stated in Proposition 3.11 and numerical investigations point in that direction.
Using the fact that is strictly decreasing and Proposition 3.10 we conclude that Moreover, it seems to be true that x 2 d 2 dx 2 log(θ 2 (x)) is strictly decreasing, which would in return imply that x 2 θ ′ 4 (x) θ 4 (x) is strictly decreasing. However, we leave this problem open and close this section with symmetry properties involving θ 2 , θ 4 and repeated applications of x d dx . Proposition 3.12. For x ∈ R + and 2 ≤ n ∈ N the following holds.
Proof. The statement is an obvious consequence of identity (3.8) and follows by induction.

3.2.
Properties of Theta-3. We will now study properties of the logarithm of θ 3 on a logarithmic scale.
This result was proved in [19]. As a consequence of Lemma 3.6 we get the following results which were already proved in the author's doctoral thesis [17].
Proposition 3.14. For x ∈ R + and 2 ≤ n ∈ N the following holds.
Proof. The statement is an obvious consequence of identity (3.7) and follows by induction.
The following statement and its proof were also already given in [17]. The techniques are similar to the techniques used by Montgomery in [25]. The computer algebra system Mathematica [32] was used at some points in the proof in order to compute explicit values or closed expressions for geometric series, but in principle all computations can also be checked by hand. The proof was adjusted to a level which should be quite accessible with no or only little help of computer algebra software. Proposition 3.15. For x ∈ R + , the function is positive for x ∈ (0, 1) and negative for x > 1. Also, the function is anti-symmetric in the following sense Proof. To simplify notation we set .
Proposition 3.14 settles the part concerning the anti-symmetry of expression (3.10). Therefore, it is also clear that we only need to prove the statement about the negativity for x > 1, the result for x ∈ (0, 1) then follows immediately. We start with the following calculation; In particular, Proposition 3.14 implies that ψ(1) + 3ψ ′ (1) + ψ ′′ (1) = 0. We proceed in two steps: (i) We will show that x d dx 3 log (θ 3 (x)) is negative for x > 1.1.
First step: We will use a combination of the asymptotic and the local behavior to establish the result claimed in (i).
Since ψ is the logarithmic derivative of θ 3 , we use Jacobi's triple product formula for θ 3 to obtain a series representation for ψ. In order to control the expression in equation (3.11) we compute the derivatives of ψ up to order 2.
It is easy to verify that ψ (n) (x) = O(e −πx ), n = 0, 1, 2 (this actually holds for any n ∈ N). However, due to the monomial term in k, which is of the order of the derivative plus 1, we will also consider the contributions of terms involving e −2πx to gain more control on the local behavior close to 1.
The techniques are standard and the proof is inspired by the proofs in Montgomery's article [25]. In particular, we will estimate parts of the series by the leading term(s), by appropriate geometric series or by their values at x = 1 (or use combinations of the mentioned methods). We will try to indicate which kind of estimates are used at which point in the proof.
We start with an estimate for ψ for x > 1. We will first estimate the denominators and after that estimate the positive part by a geometric series and the negative part by the leading term. This will also be the common theme throughout the proof.
We proceed with an upper bound for ψ ′ for x > 1.
In the second to last line we used the fact that (for |q| < 1) The closed expressions for the according geometric series are obtained by taking derivatives of the classical geometric series and for the series involving only even term, q is substituted by q 2 .
With the same techniques we bound ψ ′′ from above for x > 1.
This time we truncated the negative series over the odd integers after 2 terms, as the monomial already has degree 3 and, therefore, the second term contributes some weight to the series evaluated near x = 1. Now, for x > 1 we have the following estimate; It is obvious that the last expression asymptotically tends to zero from below, hence, there exists some x 0 > 1 such that We will now determine x 0 . By the above estimates, it is enough to show that the following inequality holds for x > x 0 > 1; It is not hard to check that −6 + 60.3 x − 55.1 x 2 is strictly decreasing for x > 1. Also, we note that expressions of the form a + b x + c x 2 e −πx are strictly decreasing for For a = 6.35, b = −118.2 and c = 255.5, we find out that for x > 1 the expression is indeed strictly decreasing. Therefore, it suffices to find a value x 0 such that (3.12) is true, as it then holds for all x > x 0 . For x = 1.1 the left-hand side of (3.12) is smaller than −0.4 and we found the desired value x 0 = 1.1.

Second step:
We will now show that x d dx 3 log (θ 3 (x)) is strictly decreasing on an interval containing (1, x 0 ). To do so, we apply the differential operator x d dx on (3.10), i.e., In order to control the last expression for x ∈ (1, x 0 ), we need good estimates on ψ ′′′ on this interval. We start in the same manner as for the lower order derivatives.
In total, it follows that for x ∈ (1, x 0 ) the following inequality holds; We need to verify that . This is equivalent to showing that In order to establish this result, we note that Hence, we get the estimate for x ∈ (1, x 0 ). Therefore, by showing the truth of the (stronger) inequality (3.14) −5.7 + 132 x − 282 x 2 + 144 x 3 < 0, x ∈ (1, x 0 ), we can finish the proof. It is not hard to show that the left-hand side of the above inequality has critical points at 11 36 / ∈ (1, x 0 ) (local maximum) and 1 (local minimum). Finally, by checking the values on the boundary of the interval (1, x 0 ) we see that (3.14) indeed holds.
By combining the two steps we see that the function x d dx 3 log (θ 3 (x)) is strictly decreasing at least on the interval (1, 1.1) and negative for x > 1.1. As the value at x = 1 is zero, we can finally conclude that the expression given in equation (3.10) is negative for x > 1. Due to the already mentioned anti-symmetry with respect to the point (1, 0) the function has to be positive for 0 < x < 1.
Vesalainen and Ernvall-Hytönen [14], [16] proved Proposition 3.15 independently from this work. They used the result to prove a conjecture by Hernandez and Sethuraman 5 [20] (see also the references in [16]). The following result by Vesalainen and Ernvall-Hytönen is an immediate consequence of Proposition 3.15.
Theorem 3.16 (Ernvall-Hytönen, Vesalainen). For α, β ∈ R + with 1 ≤ α < β and x ∈ R + we define the function . 5 The author does not claim any credit for solving the conjecture by Hernandez and Sethuraman since the author was not aware of the conjecture prior to the appearance of [16].
Then g(x) = g( 1 x ) and g assumes its global maximum only for x = 1. Also, the function is strictly increasing for x ∈ (0, 1) and strictly decreasing for x ∈ (1, ∞).

3.3.
A Result Involving Theta-1 and a Proof of Theorem 1.1. In contrast to the fact that ϑ 1 (0, τ ) = 0 for all τ ∈ H, its derivative with respect to z evaluated at (0, ix), x ∈ R + is not the zero function, in particular The above function is often denoted by θ ′ 1 (x) (see e.g., [11]), where the prime indicates differentiation with respect to z. The representation as an infinite product is also well-known and follows from the Jacobi triple product representation of Θ (Proposition 3.2). Since the notation θ ′ 1 could be misleading in this work, we rather use the notation from [8] and set θ + 1 (x) = ∂ z ϑ 1 (0, ix), x ∈ R + . Now we have the following identity (see e.g., [8, chap. 3] or [11,Chap. 4]) (3.16) θ 2 (x)θ 3 (x)θ 4 (x) = θ + 1 (x). From this identity we derive the following result. The function is positive and assumes its global maximum only for x = 1 and x ψ ′ 1 (x) ψ 1 (x) is strictly decreasing. Furthermore, ψ 1 (x) = ψ 1 ( 1 x ) and, hence, Proof. The different representations of ψ 1 follow immediately from (3.15) and (3.16). It is obvious from the product representation that ψ 1 is positive. The symmetry of ψ 1 follows easily from (3.5) and (3.6). The (anti-)symmetry for the repeated application of the differential operator x d dx immediately follows by induction. Differentiating log(ψ 1 ) on a logarithmic scale, i.e., applying the differential operator x d dx , yields .
In particular, this means that ψ 1 has a critical point at x = 1 since which follows from (3.7) and (3.8). The fact that ψ 1 assumes its global maximum only for x = 1 will follow from the fact that x ψ ′ 1 (x) ψ 1 (x) is strictly decreasing because, since x > 0 and ψ 1 (x) > 0 we have that Since ψ ′ 1 (1) ψ 1 (1) = 0, the negativity of ψ ′ 1 for x > 1 already follows if we can show that x ψ ′ 1 (x) ψ 1 (x) is decreasing. We use the infinite product representation to establish this result. It is quickly verified that, except for the constant term, all terms in this expression are decreasing and the proof is finished. Figure 1. The function ψ 1 (x) and its derivative x d dx ψ 1 (x) plotted on a logarithmic scale. Figure 1 shows the function ψ 1 (x) and its derivative with logarithmic scaling on the x-axis, revealing the symmetry properties. In Figure 2 we see the logarithmic derivative of ψ 1 on a logarithmic scale, i.e., x d dx log(ψ 1 (x)). The behavior of this function is strongly influenced by the above established properties of Jacobi's theta functions.