Lemniscate constant

From HandWiki
Short description: Ratio of the perimeter of Bernoulli's lemniscate to its diameter
Lemniscate of Bernoulli

In mathematics, the lemniscate constant ϖ is a transcendental mathematical constant that is the ratio of the perimeter of Bernoulli's lemniscate to its diameter, analogous to the definition of π for the circle.[1] Equivalently, the perimeter of the lemniscate (x2+y2)2=x2−y2 is 2ϖ. The lemniscate constant is closely related to the lemniscate elliptic functions and is approximately equal to 2.62205755.[2] It also appears in evaluation of the gamma and beta function at certain rational values. The symbol ϖ is a cursive variant of π known as variant pi represented in Unicode by the character U+03D6 ϖ .

Sometimes the quantities 2ϖ or ϖ/2 are referred to as the lemniscate constant.[3][4]

History

Gauss's constant, denoted by G, is equal to ϖ / π ≈ 0.8346268[5] and named after Carl Friedrich Gauss, who calculated it via the arithmetic–geometric mean as 1/M(1,2).[6] By 1799, Gauss had two proofs of the theorem that M(1,2)=π/ϖ where ϖ is the lemniscate constant.[7]

John Todd named two more lemniscate constants, the first lemniscate constant A = ϖ/2 ≈ 1.3110287771 and the second lemniscate constant B = π/(2ϖ) ≈ 0.5990701173.[8][9][10]

The lemniscate constant ϖ and Todd's first lemniscate constant A were proven transcendental by Carl Ludwig Siegel in 1932 and later by Theodor Schneider in 1937 and Todd's second lemniscate constant B and Gauss's constant G were proven transcendental by Theodor Schneider in 1941.[8][11][12] In 1975, Gregory Chudnovsky proved that π and ϖ are algebraically independent over ℚ, which implies that A and B are algebraically independent as well.[13][14] But the set {π,M(1,1/2),M′(1,1/2)} (where the prime denotes the derivative with respect to the second variable) is not algebraically independent over ℚ.[15] In fact, π=22M3(1,12)M′(1,12)=1G3M′(1,12). In 1996, Yuri Nesterenko proved that eπ, π, and ϖ are algebraically independent over ℚ.[16]

As of 2025 over 2 trillion digits of this constant have been calculated using y-cruncher.[17]

Forms

Usually, ϖ is defined by the first equality below, but it has many equivalent forms:[18]

ϖ=2∫01dt1−t4=2∫0∞dt1+t4=∫01dtt−t3=∫1∞dtt3−t=4∫0∞(1+t44−t)dt=22∫011−t44dt=3∫011−t4dt=2K(i)=12B(14,12)=122B(14,14)=Γ(1/4)222π=2−24ζ(3/4)2ζ(1/4)2=2.62205755429211981046483958989111941…,

where K is the complete elliptic integral of the first kind with modulus k, Β is the beta function, Γ is the gamma function and ζ is the Riemann zeta function.

The lemniscate constant can also be computed by the arithmetic–geometric mean M,

ϖ=πM(1,2).

Gauss's constant is typically defined as the reciprocal of the arithmetic–geometric mean of 1 and the square root of 2, after his calculation of M(1,2) published in 1800:[19]G=1M(1,2)John Todd's lemniscate constants may be given in terms of the beta function B: A=ϖ2=14B(14,12),B=π2ϖ=14B(12,34).

As a special value of L-functions

β′(0)=log⁡ϖπ

which is analogous to

ζ′(0)=log⁡12π

where β is the Dirichlet beta function and ζ is the Riemann zeta function.[20]

Analogously to the Leibniz formula for π, β(1)=∑n=1∞χ(n)n=π4, we have[21][22][23][24][25] L(E,1)=∑n=1∞ν(n)n=ϖ4 where L is the L-function of the elliptic curve E:y2=x3−x over ℚ; this means that ν is the multiplicative function given by ν(pn)={p−𝒩p,p∈ℙ,n=10,p=2,n≥2ν(p)ν(pn−1)−pν(pn−2),p∈ℙ∖{2},n≥2 where 𝒩p is the number of solutions of the congruence a3−a≡b2(mod⁡p),p∈ℙ in variables a,b that are non-negative integers (ℙ is the set of all primes). Equivalently, ν is given by F(τ)=η(4τ)2η(8τ)2=∑n=1∞ν(n)qn,q=e2πiτ where τ∈ℂ such that ℑ⁡τ>0 and η is the eta function.[26][27][28] The above result can be equivalently written as ∑n=1∞ν(n)ne−2πn/32=ϖ8 (the number 32 is the conductor of E) and also tells us that the BSD conjecture is true for the above E.[29] The first few values of ν are given by the following table; if 1≤n≤113 such that n doesn't appear in the table, then ν(n)=0: nν(n)nν(n)1153145−261−109−365−1213673−617281925−185−429−10891037−297184110101−2456109649−7113−14

As a special value of other functions

Let Δ be the minimal weight level 1 new form. Then[30] Δ(i)=164(ϖπ)12. The q-coefficient of Δ is the Ramanujan tau function.

Series

Viète's formula for π can be written:

2π=12⋅12+1212⋅12+1212+1212⋯

An analogous formula for ϖ is:[31]

2ϖ=12⋅12+12/12⋅12+12/12+12/12⋯

The Wallis product for π is:

π2=∏n=1∞(1+1n)(−1)n+1=∏n=1∞(2n2n−1⋅2n2n+1)=(21⋅23)(43⋅45)(65⋅67)⋯

An analogous formula for ϖ is:[32]

ϖ2=∏n=1∞(1+12n)(−1)n+1=∏n=1∞(4n−14n−2⋅4n4n+1)=(32⋅45)(76⋅89)(1110⋅1213)⋯

A related result for Gauss's constant (G=ϖ/π) is:[33]

ϖπ=∏n=1∞(4n−14n⋅4n+24n+1)=(34⋅65)(78⋅109)(1112⋅1413)⋯

An infinite series discovered by Gauss is:[34]

ϖπ=∑n=0∞(−1)n∏k=1n(2k−1)2(2k)2=1−1222+12⋅3222⋅42−12⋅32⋅5222⋅42⋅62+⋯

The Machin formula for π is 14π=4arctan⁡15−arctan⁡1239, and several similar formulas for π can be developed using trigonometric angle sum identities, e.g. Euler's formula 14π=arctan⁡12+arctan⁡13. Analogous formulas can be developed for ϖ, including the following found by Gauss: 12ϖ=2arcsl⁡12+arcsl⁡723, where arcsl is the lemniscate arcsine.[35]

The lemniscate constant can be rapidly computed by the series[36][37]

ϖ=2−1/2π(∑n∈ℤe−πn2)2=21/4πe−π/12(∑n∈ℤ(−1)ne−πpn)2

where pn=12(3n2−n) (these are the generalized pentagonal numbers). Also[38]

∑m,n∈ℤe−2π(m2+mn+n2)=1+3ϖ121/8π.

In a spirit similar to that of the Basel problem,

∑z∈ℤ[i]∖{0}1z4=G4(i)=ϖ415

where ℤ[i] are the Gaussian integers and G4 is the Eisenstein series of weight 4 (see Lemniscate elliptic functions § Hurwitz numbers for a more general result).[39]

A related result is

∑n=1∞σ3(n)e−2πn=ϖ480π4−1240

where σ3 is the sum of positive divisors function.[40]

In 1842, Malmsten found

β′(1)=∑n=1∞(−1)n+1log⁡(2n+1)2n+1=π4(γ+2log⁡πϖ2)

where γ is Euler's constant and β(s) is the Dirichlet-Beta function.

The lemniscate constant is given by the rapidly converging series

ϖ=π324e−π3(∑n=−∞∞(−1)ne−2nπ(3n+1))2.

The constant is also given by the infinite product

ϖ=π∏m=1∞tanh2(πm2).

Also[41]

∑n=0∞(−1)n6635520n(4n)!n!4=2457/4ϖ2π2.

Continued fractions

A (generalized) continued fraction for π is π2=1+11+1⋅21+2⋅31+3⋅41+⋱ An analogous formula for ϖ is[9] ϖ2=1+12+2⋅32+4⋅52+6⋅72+⋱

Define Brouncker's continued fraction by[42] b(s)=s+122s+322s+522s+⋱,s>0. Let n≥0 except for the first equality where n≥1. Then[43][44] b(4n)=(4n+1)∏k=1n(4k−1)2(4k−3)(4k+1)πϖ2b(4n+1)=(2n+1)∏k=1n(2k)2(2k−1)(2k+1)4πb(4n+2)=(4n+1)∏k=1n(4k−3)(4k+1)(4k−1)2ϖ2πb(4n+3)=(2n+1)∏k=1n(2k−1)(2k+1)(2k)2π. For example, b(1)=4π,b(2)=ϖ2π,b(3)=π,b(4)=9πϖ2.

In fact, the values of b(1) and b(2), coupled with the functional equation b(s+2)=(s+1)2b(s), determine the values of b(n) for all n.

Simple continued fractions

Simple continued fractions for the lemniscate constant and related constants include[45][46] ϖ=[2,1,1,1,1,1,4,1,2,…],2ϖ=[5,4,10,2,1,2,3,29,…],ϖ2=[1,3,4,1,1,1,5,2,…],ϖπ=[0,1,5,21,3,4,14,…].

Integrals

A geometric representation of ϖ/2 and ϖ/2

The lemniscate constant ϖ is related to the area under the curve x4+y4=1. Defining πn:=B(1n,1n), twice the area in the positive quadrant under the curve xn+yn=1 is 2∫011−xnndx=1nπn. In the quartic case, 14π4=12ϖ.

In 1842, Malmsten discovered that[47]

∫01log⁡(−log⁡x)1+x2dx=π2log⁡πϖ2.

Furthermore, ∫0∞tanh⁡xxe−xdx=log⁡ϖ2π

and[48]

∫0∞e−x4dx=2ϖ2π4,analogous to∫0∞e−x2dx=π2, a form of Gaussian integral.

The lemniscate constant appears in the evaluation of the integrals

πϖ=∫0π2sin⁡(x)dx=∫0π2cos⁡(x)dx

ϖπ=∫0∞dxcosh⁡(πx)

John Todd's lemniscate constants are defined by integrals:[8]

A=∫01dx1−x4

B=∫01x2dx1−x4

Circumference of an ellipse

The lemniscate constant satisfies the equation[49]

πϖ=2∫01x2dx1−x4

Euler discovered in 1738 that for the rectangular elastica (first and second lemniscate constants)[50][49]

arc length⋅height=A⋅B=∫01dx1−x4⋅∫01x2dx1−x4=ϖ2⋅π2ϖ=π4

Now considering the circumference C of the ellipse with axes 2 and 1, satisfying 2x2+4y2=1, Stirling noted that[51]

C2=∫01dx1−x4+∫01x2dx1−x4

Hence the full circumference is

C=πϖ+ϖ=3.820197789…

This is also the arc length of the sine curve on half a period:[52]

C=∫0π1+cos2(x)dx

Other limits

Analogously to 2π=limn→∞|(2n)!B2n|12n where Bn are Bernoulli numbers, we have 2ϖ=limn→∞((4n)!H4n)14n where Hn are Hurwitz numbers.

Notes

  1. ↑ See:
  2. ↑ See:
  3. ↑ "A064853 - Oeis". http://oeis.org/A064853. 
  4. ↑ "Lemniscate Constant". http://www.numberworld.org/digits/Lemniscate/. 
  5. ↑ "A014549 - Oeis". http://oeis.org/A014549. 
  6. ↑ Finch 2003, p. 420.
  7. ↑ Neither of these proofs was rigorous from the modern point of view. See Cox 1984, p. 281
  8. ↑ 8.0 8.1 8.2 Todd, John (January 1975). "The lemniscate constants". Communications of the ACM 18 (1): 14–19. doi:10.1145/360569.360580. 
  9. ↑ 9.0 9.1 "A085565 - Oeis". http://oeis.org/A085565.  and "A076390 - Oeis". http://oeis.org/A076390. 
  10. ↑ Carlson, B. C. (2010), "Elliptic Integrals", in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F. et al., NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0-521-19225-5, http://dlmf.nist.gov/19.20.E2 
  11. ↑ In particular, Siegel proved that if G4(ω1,ω2) and G6(ω1,ω2) with Im⁡(ω2/ω1)>0 are algebraic, then ω1 or ω2 is transcendental. Here, G4 and G6 are Eisenstein series. The fact that ϖ is transcendental follows from G4(ϖ,ϖi)=1/15 and G6(ϖ,ϖi)=0.
    Apostol, T. M. (1990). Modular Functions and Dirichlet Series in Number Theory (Second ed.). Springer. p. 12. ISBN 0-387-97127-0. 
    Siegel, C. L. (1932). "Über die Perioden elliptischer Funktionen." (in German). Journal für die reine und angewandte Mathematik 167: 62–69. doi:10.1515/crll.1932.167.62. https://eudml.org/doc/149791. 
  12. ↑ In particular, Schneider proved that the beta function B(a,b) is transcendental for all a,b∈ℚ∖ℤ such that a+b∉ℤ0−. The fact that ϖ is transcendental follows from ϖ=12B(14,12) and similarly for B and G from B(12,34).
    Schneider, Theodor (1941). "Zur Theorie der Abelschen Funktionen und Integrale". Journal für die reine und angewandte Mathematik 183 (19): 110–128. doi:10.1515/crll.1941.183.110. https://www.deepdyve.com/lp/de-gruyter/zur-theorie-der-abelschen-funktionen-und-integrale-mn0U50bvkB. 
  13. ↑ G. V. Choodnovsky: Algebraic independence of constants connected with the functions of analysis, Notices of the AMS 22, 1975, p. A-486
  14. ↑ G. V. Chudnovsky: Contributions to The Theory of Transcendental Numbers, American Mathematical Society, 1984, p. 6
  15. ↑ Borwein, Jonathan M.; Borwein, Peter B. (1987). Pi and the AGM: A Study in Analytic Number Theory and Computational Complexity (First ed.). Wiley-Interscience. ISBN 0-471-83138-7.  p. 45
  16. ↑ Nesterenko, Y. V.; Philippon, P. (2001). Introduction to Algebraic Independence Theory. Springer. p. 27. ISBN 3-540-41496-7. 
  17. ↑ Yee, Alexander J. (May 18, 2025). "Records set by y-cruncher". http://numberworld.org/y-cruncher/records.html. 
  18. ↑ See:
  19. ↑ Cox 1984, p. 277.
  20. ↑ "A113847 - Oeis". http://oeis.org/A113847. 
  21. ↑ Cremona, J. E. (1997). Algorithms for Modular Elliptic Curves (2nd ed.). Cambridge University Press. ISBN 0521598206. https://books.google.com/books?id=MtM8AAAAIAAJ.  p. 31, formula (2.8.10)
  22. ↑ In fact, the series ∑n=1∞ν(n)ns converges for ℜ⁡s>5/6.
  23. ↑ Murty, Vijaya Kumar (1995). Seminar on Fermat's Last Theorem. American Mathematical Society. p. 16. ISBN 9780821803134. 
  24. ↑ Cohen, Henri (1993). A Course in Computational Algebraic Number Theory. Springer-Verlag. pp. 382–406. ISBN 978-3-642-08142-2. 
  25. ↑ "Elliptic curve with LMFDB label 32.a3 (Cremona label 32a2)". https://www.lmfdb.org/EllipticCurve/Q/32/a/3. 
  26. ↑ The function F is the unique weight 2 level 32 new form and it satisfies the functional equation
    F(−1τ)=−τ232F(τ132).
  27. ↑ The ν function is closely related to the ξ function which is the multiplicative function defined by
    ξ(pn)={𝒩p′,p∈ℙ,n=1ξ(pn−1)+χ(p)n,p∈ℙ,n≥2
    where 𝒩p′ is the number of solutions of the equation
    a2+b2=p,p∈ℙ
    in variables a,b that are non-negative integers (see Fermat's theorem on sums of two squares) and χ is the Dirichlet character from the Leibniz formula for π; also
    ∑d|nχ(d)=ξ(n)
    for any positive integer n where the sum extends only over positive divisors; the relation between ν and ξ is
    ∑k=0n(−1)kξ(4k+1)ξ(4n−4k+1)=ν(2n+1)
    where n is any non-negative integer.
  28. ↑ The ν function also appears in
    ∑z∈𝔾;zz‾=nz=ν(n)
    where n is any positive integer and 𝔾 is the set of all Gaussian integers of the form
    (−1)a±b−12(a±bi)
    where a is odd and b is even. The ξ function from the previous note satisfies
    |{z:z∈𝔾∧zz‾=n}|=ξ(n)
    where n is positive odd.
  29. ↑ Rubin, Karl (1987). "Tate-Shafarevich groups and L-functions of elliptic curves with complex multiplication". Inventiones Mathematicae 89 (3): 528. doi:10.1007/BF01388984. Bibcode: 1987InMat..89..527R. https://eudml.org/doc/143493. 
  30. ↑ "Newform orbit 1.12.a.a". https://www.lmfdb.org/ModularForm/GL2/Q/holomorphic/1/12/a/a/. 
  31. ↑ Levin (2006)
  32. ↑ Hyde (2014) proves the validity of a more general Wallis-like formula for clover curves; here the special case of the lemniscate is slightly transformed, for clarity.
  33. ↑ Hyde, Trevor (2014). "A Wallis product on clovers". The American Mathematical Monthly 121 (3): 237–243. doi:10.4169/amer.math.monthly.121.03.237. https://math.uchicago.edu/~tghyde/Hyde%20--%20A%20Wallis%20product%20on%20clovers.pdf. Retrieved 2021-10-29. 
  34. ↑ Bottazzini, Umberto; Gray, Jeremy (2013). Hidden Harmony – Geometric Fantasies: The Rise of Complex Function Theory. Springer. doi:10.1007/978-1-4614-5725-1. ISBN 978-1-4614-5724-4.  p. 60
  35. ↑ Todd (1975)
  36. ↑ Cox 1984, p. 307, eq. 2.21 for the first equality. The second equality can be proved by using the pentagonal number theorem.
  37. ↑ Berndt, Bruce C. (1998). Ramanujan's Notebooks Part V. Springer. ISBN 978-1-4612-7221-2.  p. 326
  38. ↑ This formula can be proved by hypergeometric inversion: Let
    a⁡(q)=∑m,n∈ℤqm2+mn+n2
    where q∈ℂ with |q|<1. Then
    a⁡(q)=2F1(13,23,1,z)
    where
    q=exp⁡(−2π32F1(1/3,2/3,1,1−z)2F1(1/3,2/3,1,z))
    where z∈ℂ∖{0,1}. The formula in question follows from setting z=14(33−5).
  39. ↑ Eymard, Pierre; Lafon, Jean-Pierre (2004). The Number Pi. American Mathematical Society. ISBN 0-8218-3246-8.  p. 232
  40. ↑ Garrett, Paul. "Level-one elliptic modular forms". http://www-users.math.umn.edu/~garrett/m/mfms/notes_2015-16/10_level_one.pdf.  p. 11—13
  41. ↑ The formula follows from the hypergeometric transformation
    3F2(14,12,34,1,1,16z(1−z)2(1+z)4)=(1+z)2F1(12,12,1,z)2
    where z=λ(1+5i) and λ is the modular lambda function.
  42. ↑ Khrushchev, Sergey (2008). Orthogonal Polynomials and Continued Fractions (First ed.). Cambridge University Press. ISBN 978-0-521-85419-1.  p. 140 (eq. 3.34), p. 153. There's an error on p. 153: 4[Γ(3+s/4)/Γ(1+s/4)]2 should be 4[Γ((3+s)/4)/Γ((1+s)/4)]2.
  43. ↑ Khrushchev, Sergey (2008). Orthogonal Polynomials and Continued Fractions (First ed.). Cambridge University Press. ISBN 978-0-521-85419-1.  p. 146, 155
  44. ↑ Perron, Oskar (1957) (in German). Die Lehre von den Kettenbrüchen: Band II (Third ed.). B. G. Teubner.  p. 36, eq. 24
  45. ↑ "A062540 - OEIS". http://oeis.org/A062540. 
  46. ↑ "A053002 - OEIS". http://oeis.org/A053002. 
  47. ↑ Blagouchine, Iaroslav V. (2014). "Rediscovery of Malmsten's integrals, their evaluation by contour integration methods and some related results". The Ramanujan Journal 35 (1): 21–110. doi:10.1007/s11139-013-9528-5. https://www.researchgate.net/publication/257381156. 
  48. ↑ "A068467 - Oeis". https://oeis.org/A068467. 
  49. ↑ 49.0 49.1 Cox 1984, p. 313.
  50. ↑ Levien (2008)
  51. ↑ Cox 1984, p. 312.
  52. ↑ Adlaj, Semjon (2012). "An Eloquent Formula for the Perimeter of an Ellipse". p. 1097. https://www.ams.org/notices/201208/rtx120801094p.pdf. "One might also observe that the length of the “sine” curve over half a period, that is, the length of the graph of the function sin⁡(t) from the point where t=0 to the point where t=π , is 2l(1/2)=L+M."  In this paper M=1/G=π/ϖ and L=π/M=Gπ=ϖ.

References