pith. sign in

arxiv: 2604.14812 · v1 · submitted 2026-04-16 · 🪐 quant-ph · math-ph· math.AP· math.MP

Time-Dependent Logarithmic Perturbation Theory for Quantum Dynamics: Formulation and Applications

Juan Carlos del Valle , Paul Bergold , Karolina Kropielnicka This is my paper

Pith reviewed 2026-05-10 11:07 UTC · model grok-4.3

classification 🪐 quant-ph math-phmath.APmath.MP
keywords time-dependent perturbation theorylogarithmic perturbation theoryquantum dynamicsSchrödinger equationenergy shiftsAC Stark shiftharmonic oscillatorhydrogen atom
0
0 comments X

The pith

Time-dependent logarithmic perturbation theory extends the log-wavefunction expansion to driven quantum systems and yields closed integral expressions for corrections and energy shifts.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a time-dependent version of logarithmic perturbation theory for the Schrödinger equation. The logarithm of the wave function is expanded in powers of a coupling constant, producing a recursive hierarchy of equations. This hierarchy is controlled by a gauge-rotated Hamiltonian of the unperturbed problem and solved by Duhamel's integral formula, preserving the closed-form character known from the time-independent case. The resulting expressions directly supply the instantaneous energy shift and time-averaged dynamic shifts that relate to AC Stark effects and polarizabilities. The approach is tested on a laser-driven harmonic oscillator, where it recovers the exact solution, and on the hydrogen atom, where it produces accurate observables for multi-photon processes.

Core claim

By expanding the logarithm of the time-dependent wave function in powers of a coupling constant, the Schrödinger equation generates a hierarchy of linear equations governed by a gauge-rotated unperturbed Hamiltonian; each order admits an explicit integral representation via Duhamel's formula, furnishing computable expressions for the instantaneous energy shift and for dynamic energy shifts expressed as time averages of pseudopotential expectation values.

What carries the argument

The gauge-rotated Hamiltonian of the unperturbed system, which replaces the original Hamiltonian in the perturbative hierarchy and permits closed integral solutions for the logarithmic corrections through Duhamel's formula.

If this is right

  • Instantaneous energy shifts become directly computable from the perturbative corrections without solving the full time-dependent problem.
  • Dynamic energy shifts appear naturally as time averages of pseudopotential expectation values and connect to AC Stark shifts and electric polarizabilities.
  • The same closed-integral structure applies to multi-photon processes in atoms under time-dependent laser fields.
  • Exact recovery of the driven harmonic oscillator solution confirms that the hierarchy reproduces known analytic results order by order.
  • Numerical tests on the hydrogen atom demonstrate that relevant observables can be obtained with high accuracy in the perturbative regime.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The method could reduce secular divergences that appear in conventional time-dependent perturbation theory when applied to periodic driving.
  • Because the corrections remain integrals over the unperturbed propagator, the framework may extend straightforwardly to Floquet-type periodic problems without additional averaging approximations.
  • For systems where the gauge rotation is analytically tractable, the approach supplies a practical route to high-order analytic expressions that standard perturbative expansions cannot achieve.

Load-bearing premise

The perturbative series in the coupling constant converges for the chosen systems and time scales, and the gauge rotation of the unperturbed Hamiltonian remains well-defined and tractable.

What would settle it

A numerical comparison in which the perturbative energy shifts or wave-function corrections for the driven harmonic oscillator deviate from the known exact solution at moderate coupling strengths or longer times would show the method fails to converge as claimed.

Figures

Figures reproduced from arXiv: 2604.14812 by Juan Carlos del Valle, Karolina Kropielnicka, Paul Bergold.

Figure 1
Figure 1. Figure 1: Real and imaginary parts of the instantaneous energy shift (II.35) as functions of [PITH_FULL_IMAGE:figures/full_fig_p019_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Time-dependent induced dipole moment calculated using first-order TDLPT and [PITH_FULL_IMAGE:figures/full_fig_p020_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Diagrammatic representation of the terms contributing to the correction Φ [PITH_FULL_IMAGE:figures/full_fig_p021_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Graphical representation of the hybrid coordinates [PITH_FULL_IMAGE:figures/full_fig_p031_4.png] view at source ↗
read the original abstract

We present a time-dependent extension of logarithmic perturbation theory for nonrelativistic quantum dynamics governed by the Schr\"odinger equation, in which the logarithm of the wave function is expanded in powers of a coupling constant. The resulting hierarchy of equations defining the perturbative corrections is governed by a gauge-rotated Hamiltonian of the unperturbed system and leads to closed-integral expressions for the time-dependent corrections based on Duhamel's formula. This closed-integral structure of corrections is a hallmark of time-independent logarithmic perturbation theory and is preserved in the present extension. This structure, in particular, provides a computable expression for the instantaneous energy shift. Furthermore, dynamic energy shifts arise naturally within this framework in the form of time-averaged expectation values of pseudopotentials and can be related, for example, to AC Stark shifts and electric polarizabilities. As an illustration, we apply the method to the harmonic oscillator and the hydrogen atom, both driven by a time-dependent laser field. The harmonic oscillator provides a proof of principle for which the exact solution is recovered, while the hydrogen atom illustrates the method applied to atomic systems. Supported by numerical simulations, we demonstrate the applicability to obtain relevant physical observables with high accuracy. The present approach offers a promising alternative for analytical studies of time-dependent multi-photon processes in the perturbative regime.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 3 minor

Summary. The manuscript develops a time-dependent extension of logarithmic perturbation theory for the Schrödinger equation. The logarithm of the wave function is expanded in powers of a coupling constant, yielding a hierarchy of equations governed by a gauge-rotated unperturbed Hamiltonian. Perturbative corrections are obtained as closed integrals via Duhamel's formula, preserving the structure of the time-independent case and furnishing explicit expressions for the instantaneous energy shift as well as time-averaged dynamic shifts (e.g., AC Stark shifts). The formalism is applied to a laser-driven harmonic oscillator, where the exact solution is recovered, and to the laser-driven hydrogen atom, where numerical results indicate high accuracy for relevant observables.

Significance. If the derivation and numerical checks hold, the work supplies a computable perturbative framework that retains the closed-integral character of logarithmic perturbation theory while extending it to driven systems. The exact recovery for the harmonic oscillator constitutes a strong internal verification of the gauge-rotation and integral-construction steps. The hydrogen-atom application illustrates utility for atomic multi-photon processes and polarizability calculations, offering a potential alternative to conventional time-dependent perturbation theory when only energy shifts are required.

major comments (2)
  1. [§2.2, Eq. (8)] §2.2, Eq. (8): The gauge-rotated Hamiltonian is introduced without an explicit demonstration that the rotation remains well-defined and invertible for the time-dependent driving term; a short appendix verifying that the operator exponential commutes appropriately with the time-dependent perturbation would remove any ambiguity in the hierarchy derivation.
  2. [§4.2, Fig. 4] §4.2, Fig. 4: The reported accuracy for the hydrogen-atom energy shift is stated as 'high' but no order-by-order convergence table or residual norm versus perturbative order is supplied; without this, it is difficult to judge the radius of convergence for the chosen laser parameters and time scales.
minor comments (3)
  1. The abstract claims the method 'provides a computable expression for the instantaneous energy shift,' yet the main text does not explicitly contrast this expression with the usual time-dependent perturbation theory formula for the same quantity.
  2. Notation: the symbol for the gauge-rotation operator is introduced in §2 but reused without redefinition in §3; a single consistent definition or a notation table would improve readability.
  3. The hydrogen-atom section would benefit from a brief statement of the basis-set size and time-stepping method used in the numerical reference calculations against which the perturbative results are compared.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of our work and the recommendation for minor revision. The comments are constructive and help improve the clarity of the derivation and the presentation of numerical results. We address each major comment below.

read point-by-point responses

  1. Referee: [§2.2, Eq. (8)] §2.2, Eq. (8): The gauge-rotated Hamiltonian is introduced without an explicit demonstration that the rotation remains well-defined and invertible for the time-dependent driving term; a short appendix verifying that the operator exponential commutes appropriately with the time-dependent perturbation would remove any ambiguity in the hierarchy derivation.

    Authors: We agree that an explicit verification strengthens the presentation. Although the gauge rotation follows the standard construction from the time-independent case and extends directly via the time-ordered exponential, we will add a short appendix in the revised manuscript. This appendix will explicitly demonstrate the commutation of the operator exponential with the time-dependent perturbation, confirming that the rotation is well-defined and invertible for the driving term under consideration. revision: yes

  2. Referee: [§4.2, Fig. 4] §4.2, Fig. 4: The reported accuracy for the hydrogen-atom energy shift is stated as 'high' but no order-by-order convergence table or residual norm versus perturbative order is supplied; without this, it is difficult to judge the radius of convergence for the chosen laser parameters and time scales.

    Authors: We appreciate this observation. To allow a clearer assessment of convergence, we will include in the revised manuscript an order-by-order table of the hydrogen-atom energy shift together with the corresponding residual norms (or L2 norms of the correction terms) as a function of perturbative order. These data will be provided for the specific laser parameters and time scales used in Figure 4, enabling readers to evaluate the practical radius of convergence. revision: yes

Circularity Check

0 steps flagged

Derivation self-contained via standard Duhamel and gauge rotation; no reduction to inputs

full rationale

The paper extends logarithmic perturbation theory to the time-dependent Schrödinger equation by expanding log(ψ) in the coupling constant, deriving a hierarchy governed by a gauge-rotated unperturbed Hamiltonian, and obtaining corrections via Duhamel's formula. These are standard independent techniques (Duhamel's formula is a classical integral equation solution; gauge rotation is a unitary transformation). The harmonic-oscillator application recovers the exact solution, confirming the structure without circularity. No self-definitional steps, fitted inputs renamed as predictions, or load-bearing self-citations appear in the derivation chain. The central claims remain independent of the target results.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The approach rests on standard assumptions of non-relativistic quantum mechanics and perturbation theory; no free parameters, invented entities, or ad-hoc axioms beyond the usual Schrödinger framework are indicated in the abstract.

axioms (2)
  • domain assumption The wave function admits a perturbative expansion of its logarithm in powers of the coupling constant for the driven systems considered.
    This is the foundational assumption of logarithmic perturbation theory extended to time dependence.
  • standard math Duhamel's formula applies to the gauge-rotated unperturbed Hamiltonian to yield closed integral expressions for corrections.
    Standard integral equation technique for time-dependent linear evolution equations.

pith-pipeline@v0.9.0 · 5543 in / 1381 out tokens · 35128 ms · 2026-05-10T11:07:48.752726+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

51 extracted references · 51 canonical work pages

  1. [1]

    However, we can resort to the PDE (II.9), rewritten in hybrid coordinates, which defines the corrections and solve it by either numerical or analytical means

    First-Order Correction In contrast to the harmonic oscillator application, the action of the semigroup e −it ˆL on Q1 =V int cannot be evaluated directly. However, we can resort to the PDE (II.9), rewritten in hybrid coordinates, which defines the corrections and solve it by either numerical or analytical means. In particular, the first-order correction s...

    work page 2000

  2. [2]

    (III.35) Equation (III.34) can be identified with the TDSE governing the evolution of ansstate, while (III.35) describes the evolution of adstate

    Second-Order Correction and the Structure of Higher-Order Terms Similarly to (III.20), we assume the following natural factorization for the second correc- tion Φ2(r, z, t) = er r Φ2,0(r, t) + z2 r2 Φ2,2(r, t) ,(III.33) where the functionsΦ 2,0 andΦ 2,2 obey the following equations i∂ t Φ2,0 =− 1 2 ∂2 r Φ2,0 + −1 r + 1 2 Φ2,0 − Φ2,2 r2 − erΦ2 1,1 2r3 ,(II...

    work page 2022

  3. [3]

    P. M. Paul, E. S. Toma, P. Breger, G. Mullot, F. Aug´ e, P. Balcou, H. G. Muller, and P. Agos- tini, Observation of a train of attosecond pulses from high harmonic generation, Science292, 1689 (2001)

    work page 2001

  4. [4]

    V´ eniard, R

    V. V´ eniard, R. Ta¨ ıeb, and A. Maquet, Phase dependence of (n+ 1)-color (n >1) ir-uv photoionization of atoms with higher harmonics, Phys. Rev. A54, 721 (1996)

    work page 1996

  5. [5]

    Bharti, D

    D. Bharti, D. Atri-Schuller, G. Menning, K. R. Hamilton, R. Moshammer, T. Pfeifer, N. Douguet, K. Bartschat, and A. Harth, Decomposition of the transition phase in multi- sideband schemes for reconstruction of attosecond beating by interference of two-photon tran- sitions, Phys. Rev. A103, 022834 (2021)

    work page 2021

  6. [6]

    J. J. Sakurai and J. Napolitano,Modern Quantum Mechanics; 2nd ed.(Addison-Wesley, San Francisco, CA, 2011)

    work page 2011

  7. [7]

    L. I. Schiff,Quantum Mechanics, 3rd ed. (McGraw-Hill, New York, 1968)

    work page 1968

  8. [8]

    L. D. Landau and L. M. Lifshitz,Quantum Mechanics Non-Relativistic Theory, Third Edition: Volume 3, 3rd ed. (Butterworth-Heinemann, 1981)

    work page 1981

  9. [9]

    K. T. Chung, New formulation of time-dependent perturbation theory, Phys. Rev.163, 1343 (1967)

    work page 1967

  10. [10]

    Simon, Resonances inn-body quantum systems with dilatation analytic potentials and the foundations of time-dependent perturbation theory, Ann

    B. Simon, Resonances inn-body quantum systems with dilatation analytic potentials and the foundations of time-dependent perturbation theory, Ann. Math.97, 247 (1973)

    work page 1973

  11. [11]

    Kutzelnigg, Stationary perturbation theory, Theor

    W. Kutzelnigg, Stationary perturbation theory, Theor. Chem. Acc.83, 263 (1992)

    work page 1992

  12. [12]

    D. W. Hone, R. Ketzmerick, and W. Kohn, Time-dependent floquet theory and absence of an adiabatic limit, Phys. Rev. A56, 4045 (1997)

    work page 1997

  13. [13]

    Rodriguez-Vega, M

    M. Rodriguez-Vega, M. Lentz, and B. Seradjeh, Floquet perturbation theory: formalism and application to low-frequency limit, New J. Phys.20, 093022 (2018)

    work page 2018

  14. [14]

    Dalgarno and A

    A. Dalgarno and A. L. Stewart, On the perturbation theory of small disturbances, Proc. R. Soc. Lond. A238, 269 (1956). 23

    work page 1956

  15. [15]

    F. M. Fern´ andez and E. A. Castro, A new implementation of stationary perturbation theory for separable quantum-mechanical problems, J. Chem. Phys.98, 6392 (1993)

    work page 1993

  16. [16]

    R. J. Swenson and S. H. Danforth, Hypervirial and Hellmann-Feynman theorems applied to anharmonic oscillators, J. Chem. Phys.57, 1734 (1972)

    work page 1972

  17. [17]

    F. M. Fernandez and J. A. Morales, Perturbation theory without wave functions for the Zeeman effect in hydrogen, Phys. Rev. A46, 318 (1992)

    work page 1992

  18. [18]

    Aharonov and C

    Y. Aharonov and C. K. Au, New approach to perturbation theory, Phys. Rev. Lett.42, 1582 (1979)

    work page 1979

  19. [19]

    C. K. Au and Y. Aharonov, Logarithmic perturbation expansions, Phys. Rev. A20, 2245 (1979)

    work page 1979

  20. [20]

    A. V. Turbiner, The eigenvalue spectrum in quantum mechanics and the nonlinearization procedure, Sov. Phys., Usp.27, 668 (1984)

    work page 1984

  21. [21]

    F. M. Fern´ andez,Introduction to Perturbation Theory in Quantum Mechanics, 1st ed. (CRC Press, Boca Raton, 2000)

    work page 2000

  22. [22]

    A. D. Dolgov and V. S. Popov, Higher orders and structure of the perturbation-theory series for the anharmonic oscillator, Zh. Eksp. Teor. Fiz.75, 2010 (1978)

    work page 2010

  23. [23]

    Eletsky, V

    V. Eletsky, V. Popov, and V. Weinberg, Logarithmic perturbation theory for the screened coulomb potential, Phys. Lett. A84, 235 (1981)

    work page 1981

  24. [24]

    Privman, New method of perturbation-theory calculation of the stark effect for the ground state of hydrogen, Phys

    V. Privman, New method of perturbation-theory calculation of the stark effect for the ground state of hydrogen, Phys. Rev. A22, 1833 (1980)

    work page 1980

  25. [25]

    J. C. del Valle, A. V. Turbiner, and A. M. E. Ruiz, Two-body neutral coulomb system in a magnetic field at rest: From hydrogen atom to positronium, Phys. Rev. A103, 032820 (2021)

    work page 2021

  26. [26]

    S. K. Bandyopadhyay and K. Bhattacharyya, Logarithmic perturbation theory: A reappraisal, Int. J. Quantum Chem.90, 27 (2002)

    work page 2002

  27. [27]

    Reed and B

    M. Reed and B. Simon,Methods of Modern Mathematical Physics IV: Analysis of Operators (Academic Press, 1978)

    work page 1978

  28. [28]

    Bader, A

    P. Bader, A. Iserles, K. Kropielnicka, and P. Singh, Efficient methods for linear Schr¨ odinger equation in the semiclassical regime with time-dependent potential, Proc. A.472, 20150733, 18 (2016)

    work page 2016

  29. [29]

    Kropielnicka and J

    K. Kropielnicka and J. C. del Valle, Second–order exponential splittings in the presence of unbounded and time-dependent operators: construction and convergence, SIAM J. Numer. 24 Anal.63, 288 (2025)

    work page 2025

  30. [30]

    Madelung, Quantentheorie in hydrodynamischer form, Z

    E. Madelung, Quantentheorie in hydrodynamischer form, Z. Phys.40, 322 (1927)

    work page 1927

  31. [31]

    H. E. Wilhelm, Hydrodynamic model of quantum mechanics, Phys. Rev. D1, 2278 (1970)

    work page 1970

  32. [32]

    F. D. Mera, S. A. Fulling, J. D. Bouas, and K. Thapa, WKB approximation to the power wall, J. Phys. A46, 195302 (2013)

    work page 2013

  33. [33]

    Goldfarb, J

    Y. Goldfarb, J. Schiff, and D. J. Tannor, Complex trajectory method in time-dependent WKB, J. Chem. Phys.128, 164114 (2008)

    work page 2008

  34. [34]

    Lasser and C

    C. Lasser and C. Lubich, Computing quantum dynamics in the semiclassical regime, Acta Numer.29, 229–401 (2020)

    work page 2020

  35. [35]

    Leung, J

    S. Leung, J. Qian, and S. Serna, Fast Huygens sweeping methods for Schr¨ odinger equations in the semi-classical regime, Methods Appl. Anal.21, 31 (2014)

    work page 2014

  36. [36]

    Hochbruck and A

    M. Hochbruck and A. Ostermann, Exponential integrators, Acta Numer.19, 209 (2010)

    work page 2010

  37. [37]

    Kato,Perturbation Theory for Linear Operators, 2nd ed., Classics in Mathematics (Springer, Berlin, Heidelberg, 1995) pp

    T. Kato,Perturbation Theory for Linear Operators, 2nd ed., Classics in Mathematics (Springer, Berlin, Heidelberg, 1995) pp. XXI + 623, see p. 488

    work page 1995

  38. [38]

    Pazy,Semigroups of Linear Operators and Applications to Partial Differential Equations, 1st ed., Applied Mathematical Sciences, Vol

    A. Pazy,Semigroups of Linear Operators and Applications to Partial Differential Equations, 1st ed., Applied Mathematical Sciences, Vol. 44 (Springer, New York, NY, 1983) p. 282

    work page 1983

  39. [39]

    A. V. Turbiner, A new approach to the eigenvalue problem in quantum mechanics: convergent perturbation theory for rising potentials, J. Phys. A14, 1641 (1981)

    work page 1981

  40. [40]

    P. W. Langhoff, S. T. Epstein, and M. Karplus, Aspects of time-dependent perturbation theory, Rev. Mod. Phys.44, 602 (1972)

    work page 1972

  41. [41]

    E. H. Kerner, Note on the forced and damped oscillator in quantum mechanics, Can. J. Phys. 36, 371 (1958)

    work page 1958

  42. [42]

    C. J. Joachain, N. J. Kylstra, and R. M. Potvliege,Atoms in Intense Laser Fields(Cambridge University Press, 2011)

    work page 2011

  43. [43]

    A. V. Turbiner, The hydrogen atom in weak electric and magnetic fields, Zh. Eksp. Teor. Fiz. 84, 1329 (1983)

    work page 1983

  44. [44]

    Douguet, M

    N. Douguet, M. Guchkov, K. Bartschat, and S. F. d. Santos, Efficient time-dependent method for strong-field ionization of atoms with smoothly varying radial steps, Atoms12, 10.3390/atoms12070034 (2024)

    work page doi:10.3390/atoms12070034 2024

  45. [45]

    Krausz and M

    F. Krausz and M. Ivanov, Attosecond physics, Rev. Mod. Phys.81, 163 (2009). 25

    work page 2009

  46. [46]

    Li, Y.-L

    P.-C. Li, Y.-L. Sheu, C. Laughlin, and S.-I. Chu, Dynamical origin of near- and below-threshold harmonic generation of Cs in an intense mid-infrared laser field, Nat. Commun.6, 7178 (2015)

    work page 2015

  47. [47]

    Reiff, T

    R. Reiff, T. Joyce, A. Jaro´ n-Becker, and A. Becker, Single-active electron calculations of high- order harmonic generation from valence shells in atoms for quantitative comparison with tddft calculations, J. Phys. Commun.4, 065011 (2020)

    work page 2020

  48. [48]

    Iserles, H

    A. Iserles, H. Z. Munthe-Kaas, S. P. Nørsett, and A. Zanna, Lie-group methods, inActa Numer., 2000, Acta Numer., Vol. 9 (Cambridge Univ. Press, Cambridge, 2000) pp. 215–365

    work page 2000

  49. [49]

    Zanna, The Fer expansion and time-symmetry: a Strang-type approach (2001) pp

    A. Zanna, The Fer expansion and time-symmetry: a Strang-type approach (2001) pp. 435–459, special issue: Themes in geometric integration

    work page 2001

  50. [50]

    Blanes, F

    S. Blanes, F. Casas, J. A. Oteo, and J. Ros, The Magnus expansion and some of its applications, Phys. Rep.470, 151 (2009)

    work page 2009

  51. [51]

    Urban, A Dyson-type series solution to some time-dependent Schr¨ odinger equations on finite adeles, J

    R. Urban, A Dyson-type series solution to some time-dependent Schr¨ odinger equations on finite adeles, J. Pseudo-Differ. Oper. Appl.16, Paper No. 24, 24 (2025). 26 Appendix A: Gauge-Rotated Representation of the Semigroup In this section, we will show that ˆLis the generator of aC 0-semigroup e−it ˆL on a weighted L2 space and admits a gauge-rotated repr...

    work page 2025