TY - JOUR

T1 - From divergent perturbation theory to an exponentially convergent self-consistent expansion

AU - Remez, Benjamin

AU - Goldstein, Moshe

N1 - Publisher Copyright: © 2018 American Physical Society.

PY - 2018/9/19

Y1 - 2018/9/19

N2 - For many nonlinear physical systems, approximate solutions are pursued by conventional perturbation theory in powers of the nonlinear terms. Unfortunately, this often produces divergent asymptotic series, collectively dismissed by Abel as "an invention of the devil." Although a lot of progress has been made on understanding the mathematics and physics behind this, new approaches are still called for. A related method, the self-consistent expansion (SCE), has been introduced by Schwartz and Edwards. Its basic idea is a rescaling of the zeroth-order system around which the solution is expanded, to achieve optimal results. While low-order SCEs have been remarkably successful in describing the dynamics of nonequilibrium many-body systems (e.g., the Kardar-Parisi-Zhang equation), its convergence properties have not been elucidated before. To address this issue we apply this technique to the canonical partition function of the classical harmonic oscillator with a quartic gx4 anharmonicity, for which perturbation theory's divergence is well-known. We obtain the Nth order SCE for the partition function, which is rigorously found to converge exponentially fast in N, and uniformly in g≥0. We use our results to elucidate the relation between the SCE and the class of approaches based on the so-called "order-dependent mapping." Moreover, we put the SCE to test against other methods that improve upon perturbation theory (Borel resummation, hyperasymptotics, Padé approximants, and the Lanczos τ-method), and find that it compares favorably with all of them for small g and dominates over them for large g. The SCE is shown to converge to the correct partition function for the double-well potential case, even when expanding around the local maximum. Our treatment is generalized to the case of many oscillators, as well as to any nonlinearity of the form g|x|q with q≥0 and complex g. These results allow us to treat the Airy function, and to see the fingerprints of Stokes lines in the SCE.

AB - For many nonlinear physical systems, approximate solutions are pursued by conventional perturbation theory in powers of the nonlinear terms. Unfortunately, this often produces divergent asymptotic series, collectively dismissed by Abel as "an invention of the devil." Although a lot of progress has been made on understanding the mathematics and physics behind this, new approaches are still called for. A related method, the self-consistent expansion (SCE), has been introduced by Schwartz and Edwards. Its basic idea is a rescaling of the zeroth-order system around which the solution is expanded, to achieve optimal results. While low-order SCEs have been remarkably successful in describing the dynamics of nonequilibrium many-body systems (e.g., the Kardar-Parisi-Zhang equation), its convergence properties have not been elucidated before. To address this issue we apply this technique to the canonical partition function of the classical harmonic oscillator with a quartic gx4 anharmonicity, for which perturbation theory's divergence is well-known. We obtain the Nth order SCE for the partition function, which is rigorously found to converge exponentially fast in N, and uniformly in g≥0. We use our results to elucidate the relation between the SCE and the class of approaches based on the so-called "order-dependent mapping." Moreover, we put the SCE to test against other methods that improve upon perturbation theory (Borel resummation, hyperasymptotics, Padé approximants, and the Lanczos τ-method), and find that it compares favorably with all of them for small g and dominates over them for large g. The SCE is shown to converge to the correct partition function for the double-well potential case, even when expanding around the local maximum. Our treatment is generalized to the case of many oscillators, as well as to any nonlinearity of the form g|x|q with q≥0 and complex g. These results allow us to treat the Airy function, and to see the fingerprints of Stokes lines in the SCE.

UR - http://www.scopus.com/inward/record.url?scp=85054501689&partnerID=8YFLogxK

U2 - https://doi.org/10.1103/PhysRevD.98.056017

DO - https://doi.org/10.1103/PhysRevD.98.056017

M3 - مقالة

SN - 2470-0010

VL - 98

JO - Physical Review D

JF - Physical Review D

IS - 5

M1 - 056017

ER -