{"Title": "Full spectral zeta values of homogeneous oscillators", "Definition": "The full spectral zeta values are $Z_N(n)=\\sum_{k=0}^{\\infty}E_{N,k}^{-n}$, where $E_{N,0}<E_{N,1}<\\cdots$ are the eigenvalues of $H_N=-d^2/dx^2+|x|^N$ on $L^2(\\mathbb R)$, each counted once. CITE{Watkins}", "Parameters": {"N": {"type": "Z", "display": "$N$", "title": "degree of the potential", "constraints": "$N\\geq1$"}, "n": {"type": "Z", "display": "$n$", "title": "inverse-power order", "constraints": "$n>(N+2)/(2N)$"}}, "Tags": ["special values", "analysis", "physics", "special functions"], "References": {"Watkins": {"bib": "J. Watkins, Spectral zeta functions of a 1D Schr\u00f6dinger problem, Journal of Nonlinear Mathematical Physics 19 (2012), 1250026, Section 2.", "arxiv": "1110.2004", "doi": "10.1142/S140292511250026X"}, "Voros1983": {"bib": "A. Voros, The return of the quartic oscillator. The complex WKB method, Annales de l\u2019Institut Henri Poincar\u00e9, Physique th\u00e9orique 39 (1983), 211\u2013338, Appendix C, especially equations (C.22), (C.23), and (C.28).", "mr": "729194", "zbl": "0526.34046"}, "Voros2000": {"bib": "A. Voros, Exercises in exact quantization, Journal of Physics A: Mathematical and General 33 (2000), 7423\u20137450, Appendix A. The corrected preprint is version 2.", "arxiv": "math-ph/0005029", "doi": "10.1088/0305-4470/33/41/311"}}, "Links": {"DLMF": {"title": "NIST DLMF: Airy differential equation and initial values", "url": "https://dlmf.nist.gov/9.2"}, "VorosArchive": {"title": "Numdam: The return of the quartic oscillator, Appendix C", "url": "https://www.numdam.org/item/AIHPA_1983__39_3_211_0/"}}, "Data properties": {"type": "R", "digits": "100", "rigour": "heuristic (agreement-checked)", "complete": "no", "complete-note": "it holds order $n=2$ for every degree $1\\leq N\\leq12$ and order $n=1$ for every degree $3\\leq N\\leq12$; the degree cutoff is a noncanonical practical prefix covering linear and harmonic controls and several steeper wells of each parity, while the first two inverse-operator moments give breadth across the degrees with explicit Gamma and hypergeometric formulas; the omitted first moments diverge", "rigour details": "Every value was computed at 150 and 200 decimal working digits. The first moments use the Gamma formula in Watkins equation (2.5a), including its factor 1/4. The second moments use equation (2.5b), with a Gamma-ratio asymptotic tail, at two cutoffs and two correction counts. The minimum agreement was 129 significant digits, and 100 digits are stored. The tail is not certified, so the rigour is heuristic (agreement-checked).\n\nAll entries were compared with separate full-trace formulas from CITE{Voros1983} and CITE{Voros2000}, with a third cutoff and working precision, agreeing to at least 156 digits. That check shares the generic summation helper. A further comparison using mpmath\u2019s own hypergeometric implementation agreed to at least 24 digits. Airy log derivatives and the harmonic odd-integer spectrum supply independent controls for degrees one and two. A known nonzero hypergeometric control, the first two quantum-Wronskian identities for every degree from 3 through 12, and a deliberate numerical perturbation were tested. The attached generator, checks, argument-selection note and numerical evidence record the methods, tolerances and limitations."}, "Display properties": {"number-header": "$Z_N(n)$", "group parameters": [["N"], ["n"]]}, "Numbers": {"1": {"2": "1.903178396159447800139065168030632969360008134672767094447247485590632327183550118473575762270501069"}, "2": {"2": "1.233700550136169827354311374984518891914212425905098828301668672027505602802400655375221675464819029"}, "3": {"1": "3.319386965493754693169476319946829002023560536117475592219264782156101837521626859818290822847120283", "2": "1.098003370680573465880414073379908567918383166689158210595599027328569244936287063325952316107760017"}, "4": {"1": "2.289908804320418350487700158170872460090165699629997424857846878235910597350966604608246696939283055", "2": "0.9963208276787715060724700465937150050135078031028254076131607718349949712416468399359828502501163302"}, "5": {"1": "1.922070491443430806944571952871298785373765682468328621508611869298918480339198212887921576327099698", "2": "0.9108397549065565510669343408855582684034704619311029148439157228421820402300559971170081979305731005"}, "6": {"1": "1.721346194668113966341474187559077603311862291076695736298545207880918722674511398540029647836713598", "2": "0.8387491545184126620672059374086251999303201267569848904935731531423690903595456176220196457332217374"}, "7": {"1": "1.589600430045033617309990840297260220894966166112886046442566142856265279624410293902798433481082535", "2": "0.7779567496716873852938031451573028143153605950969387636150390100708468296456187824033540073685081913"}, "8": {"1": "1.493995970626283323606449483976215692057478061489339631271175745340327534965228986029209139976675529", "2": "0.7264905365509530523276104496548539836689873482459224335917337074144197417298946322935003813605340313"}, "9": {"1": "1.420218357421152822820077921486188836704305652818475236813235253633683930709727361091504138335855076", "2": "0.6826375530927969955090308776399444955919809295822688294987029186446595921535004248801228803211452467"}, "10": {"1": "1.360906080874816610219754929814874604665365898681591710881641854039578548249194519520555462377414348", "2": "0.6449883867695178073696874156734389647644144971529971375192277476286622972307630513117654251004663053"}, "11": {"1": "1.311818983165571593593679839008201567953562395209460172667092562748054778746015293732248740728401489", "2": "0.6124120914440976069306457015509485417142925844630748902282096534948009200846493142461269316146061478"}, "12": {"1": "1.270306261143288234259549964910521536697501322389553580887760188612699158305650179350427774969227126", "2": "0.5840093424641228311308246685518228896715669841124739760982547285404538520450949108939539694109215394"}}, "Comments": {"normalization": "The kinetic coefficient and potential coefficient are both $1$, the additive energy shift is $0$, and the ground-state index is $k=0$. The potential is $|x|^N$ for odd as well as even degrees, and all eigenvalues are positive.", "parity": "The even-indexed eigenvalues form the half-line Neumann spectrum, and the odd-indexed eigenvalues form the half-line Dirichlet spectrum. Their inverse-power sums are added here, with no parity sign. The associated regularized determinants are in HREF{T540}[the Neumann determinant table] and HREF{T539}[the Dirichlet determinant table]. CITE{Watkins}", "convergence": "The asymptotic growth $E_{N,k}\\asymp k^{2N/(N+2)}$ gives absolute convergence for $n>(N+2)/(2N)$. No analytic continuation or finite-part value is used in these entries. CITE{Watkins}", "moments": "Integer orders are traces of inverse powers of the operator. They also give the coefficients of the expansion of the logarithmic spectral determinant at zero energy, starting at the first convergent order. CITE{Watkins}, equation (2.9); CITE{Voros2000}, equations (9)\u2013(11).", "linear": "For $N=1$, the two parity sectors give the HREF{Airy_zeta_values}[Airy zeta values] and the HREF{Airy_derivative_zeta_values}[Airy derivative zeta values]. The full trace is their sum."}, "Formulas": {"first-moment": "With $\\nu=1/(N+2)$ and $N>2$, $\\displaystyle Z_N(1)=\\frac{\\tan(2\\pi\\nu)}{\\tan(\\pi\\nu)}\\frac{\\sin(\\pi\\nu)}{2\\sqrt{\\pi}}(2\\nu)^{2-4\\nu}\\frac{\\Gamma(\\nu)\\Gamma(2\\nu)\\Gamma(3\\nu)}{\\Gamma(2\\nu+1/2)}$, where $\\Gamma$ is the HREF{Values_of_the_Gamma_function_at_rational_numbers}[Gamma function]. CITE{Voros2000}, equation (140).", "linear": "$Z_1(2)=\\kappa^2+\\kappa^{-1}$, where $\\kappa=-\\operatorname{Ai}^{\\prime}(0)/\\operatorname{Ai}(0)=3^{1/3}\\Gamma(2/3)/\\Gamma(1/3)$. The two terms are respectively the HREF{Airy_zeta_values#2}[Dirichlet] and HREF{Airy_derivative_zeta_values#2}[Neumann] second moments. CITE{DLMF}", "harmonic": "$E_{2,k}=2k+1$ and $Z_2(n)=(1-2^{-n})\\zeta(n)$ for integers $n\\geq2$, so $Z_2(2)=\\pi^2/8$. Here $\\zeta$ is the HREF{Values_of_the_Riemann_zeta_function_at_rational_numbers}[Riemann zeta function]. CITE{Voros2000}, Appendix A.2.3."}, "Similar tables": [{"table": "HREF{Eigenvalues_of_the_pure_quartic_sextic_octic_and_decic_oscillators}[Eigenvalues of pure anharmonic oscillators]", "relation": "gives individual full-line levels for four even degrees, while each trace here sums the entire infinite spectrum"}, {"table": "HREF{Airy_zeta_values}[Airy zeta values]", "relation": "gives the Dirichlet contribution at degree one"}, {"table": "HREF{Airy_derivative_zeta_values}[Airy derivative zeta values]", "relation": "gives the Neumann contribution at degree one"}, {"table": "HREF{T539}[Dirichlet spectral determinants of homogeneous oscillators]", "relation": "uses the odd-parity spectrum of the same potentials"}, {"table": "HREF{T540}[Neumann spectral determinants of homogeneous oscillators]", "relation": "uses the even-parity spectrum of the same potentials"}], "Keywords": ["homogeneous oscillator", "anharmonic oscillator", "Schrodinger", "Schr\u00f6dinger", "spectral zeta", "full spectral trace", "inverse-power trace", "spectral sum rule", "linear oscillator", "harmonic oscillator", "quartic oscillator"], "Programs": {"program-python": {"code": "# First moment for any integer degree N > 2.\nfrom mpmath import mp\nmp.dps = 120\nN = 4\nv = mp.mpf(1)/(N+2)\nA = mp.sin(mp.pi*v)/(2*mp.sqrt(mp.pi))*(2*v)**(2-4*v)\nA *= mp.gamma(v)*mp.gamma(2*v)*mp.gamma(3*v)/mp.gamma(2*v+mp.mpf(1)/2)\nprint(A*mp.tan(2*mp.pi*v)/mp.tan(mp.pi*v))", "language": "Python"}}}