Pith. sign in

REVIEW 3 major objections 5 minor 79 references

Time-dependent Gaussian basis sets for many-body systems using Rothe's method: A mean-field study

T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read This paper shows that Rothe's method reformulates TDHF and TDDFT orbital propagation as an optimization problem, and that 30–100 thawed complex Gaussians reproduce grid results for one-dimensional molecules in strong fields.

desk verdict A credible proof-of-principle for Rothe's method on mean-field orbital equations, with a real error-metric concern and an overstated 'grid-free' abstract. read the letter →

arxiv 2506.10701 v2 pith:3ZU3576Y submitted 2025-06-12 physics.chem-ph physics.comp-ph

Simon Elias Schrader , H{aa}kon Emil Kristiansen , Thomas Bondo Pedersen , Simen Kvaal This is my paper · ORCID
classification physics.chem-phphysics.comp-ph
keywords Rothe'smethodtime-dependentHartree-FockTDDFTthawedGaussianbasishigh-harmonicgenerationstrong-fielddynamicsadaptivesetsone-dimensionalmolecularmodels
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Strong-field processes such as high-harmonic generation push electrons into continuum states that ordinary atom-centered Gaussian basis sets cannot describe. This paper argues that Rothe's method, which replaces time stepping by a sequence of optimization problems, can be applied directly to the time-dependent Hartree-Fock (TDHF) and time-dependent density functional theory (TDDFT) orbital equations. In that formulation, thawed complex Gaussians with moving centers, momenta, and complex widths can follow unbound electrons, and the basis can be enlarged or pruned adaptively. For one-dimensional LiH and (LiH)$_2$ in 750 nm pulses, a few such Gaussians give qualitatively correct dipole moments and harmonic spectra, while 30–100 Gaussians reproduce grid reference calculations quantitatively at intensities up to $4\times10^{14}$ W/cm$^2$. If this holds in higher dimensions, it offers a compact and stable alternative to grid-based strong-field many-body simulations.

What carries the argument

The load-bearing object is the orbital Rothe error, Eqs. (27)–(28): $r_{i+1}(c,\alpha) = \sqrt{\sum_j \lVert \sum_m c_{j,m}\tilde{A}_i g_m(\alpha) - \tilde{A}_i^\dagger \varphi_j(t_i)\rVert^2}$, the $L^2$ mismatch between the Crank-Nicolson target and a Gaussian expansion of each orbital at the new time. The basis functions are thawed complex Gaussians $g_m(\alpha)=d_m\exp[-(a_m^2+ib_m)(x-\mu_m)^2+ip_m(x-\mu_m)]$, whose position, momentum, and complex width are optimized. Variable projection removes the linear coefficients analytically, adaptive addition and replacement are triggered when the optimized error exceeds $\varepsilon_{\Delta t}/N_T$, and a masking function absorbs outgoing flux. These pieces turn time propagation into a sequence of small nonlinear least-squares problems, avoiding the stiff ODEs and Gramian regularization of Gaussian wave-packet dynamics.

What would settle it

Rerun the LiH TDDFT case at $4\times10^{14}$ W/cm$^2$ and the high-accuracy LiH TDHF case with the Rothe error evaluated by analytic Gaussian integrals instead of grid quadrature, or on a grid so large that the width bounds are irrelevant, and compare the adaptive Gaussian trajectories and final densities with the sinc-DVR reference; if the basis evolution changes materially or the final densities diverge, the grid-based error metric is steering the method rather than merely measuring it.

Watch

Extended reading notes

Core claim

The central claim is that the orbital equations of TDHF and TDDFT can be propagated with a time-dependent Gaussian basis using Rothe's method. Over each time step the orbital-dependent mean field is frozen at its midpoint value and Crank-Nicolson propagation defines a target orbital; the new orbitals are then found by minimizing the orbital Rothe error $r_{i+1}$, the sum over occupied orbitals of the squared $L^2$ distance between the target and the best available Gaussian expansion. Because the linear coefficients are eliminated analytically via variable projection, only the nonlinear Gaussian parameters (position, momentum, real and imaginary width) are optimized at each step. The method monitors this error and adaptively adds, removes, or replaces Gaussians when the basis becomes insufficient or overcomplete, and a masking function absorbs outgoing density before re-orthonormalization. On the one-dimensional model systems tested, the result is systematic convergence to the grid reference: a handful of thawed Gaussians already captures the qualitative continuum dynamics, and 30–100 reproduce the grid dipole moments and high-harmonic spectra up to $4\times10^{14}$ W/cm$^2$.

Load-bearing premise

The adaptive strategy rests on the premise that the residual measured on the finite quadrature grid, after masking and with Gaussian widths bounded, honestly reflects how far the Gaussian orbitals are from the true propagated orbitals; if this metric misleads, the basis will grow in the wrong places and the agreement with grid calculations would not transfer to other systems or observables.

Editorial extensions

If this is right

  • Strong-field high-harmonic and ionization simulations for small molecules can be run with tens of time-adaptive Gaussians instead of large grids, with accuracy controlled by a single Rothe-error tolerance.
  • A handful of thawed Gaussians suffices for qualitatively correct continuum dynamics, so cheap exploratory calculations are possible before investing in larger bases.
  • Basis growth is systematic: adding Gaussians and lowering the tolerance reduces the discrepancy with grid densities by about two orders of magnitude in the tested cases.
  • Because the method works at the level of orbital equations, the same machinery can be carried over to correlated orbital-based methods such as time-dependent coupled cluster or MCTDH, as the paper itself anticipates.
  • The numerical difficulties of Gaussian wave-packet propagation based on equations of motion, namely stiff ODEs, non-invertible Gramians, and regularization bias, are bypassed because Gaussian parameters are found by optimization rather than by integration.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the Rothe error is instead evaluated with analytic Gaussian integrals, the current grid bottleneck disappears and the per-step cost should scale quadratically in the number of Gaussians, since the paper identifies the grid quadrature as the dominant computational cost; that is what would make the three-dimensional extension practical.
  • The masking step breaks unitarity before re-orthonormalization, so the present convergence evidence is tied to one-body observables such as dipole moments, HHG spectra, and densities; an obvious test is whether ionization yields or momentum distributions converge at the same Gaussian counts.
  • Comparing the error metric's behavior between TDHF and TDDFT could separate functional-driven complexity from basis-driven complexity in strong fields, since the paper finds TDDFT needs more Gaussians than TDHF for the same systems and intensity.
  • A direct test of the 'few Gaussians' claim would be to use four thawed Gaussians as a warm start for a correlated or configuration-interaction calculation, checking whether the subspace they span actually contains the dominant ionized configurations.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The manuscript applies Rothe's method, previously developed for the time-dependent Schrödinger equation, to the time-dependent Hartree-Fock (TDHF) and TDDFT orbital equations. The propagation step is reformulated as a variational optimization problem over thawed, complex-valued Gaussian basis functions with time-dependent widths, positions, and momenta. The authors demonstrate the approach on one-dimensional LiH and (LiH)2 model systems in 750 nm laser pulses at intensities 10^14 and 4×10^14 W/cm^2, comparing dipole moments, HHG spectra, and final densities against sinc-DVR grid reference calculations. They report that a small number of thawed Gaussians gives qualitatively correct spectra, and that increasing the number of Gaussians up to about 100 yields quantitative agreement for most systems, with systematic convergence of observables. The paper also includes a discussion of the role of the cumulative Rothe error and acknowledges that it is not a global error bound.

Significance. If the approach is validated, it offers a potentially attractive route to strong-field many-body dynamics without the exponential scaling of grid methods in higher dimensions, building on the authors' prior Rothe-method work. The paper is a proof-of-principle with clear derivations, well-documented numerical settings, and careful benchmarking against DVR references; these strengths are substantial. The main contribution is the extension of Rothe's method to orbital-based mean-field theories, which is a nontrivial step toward correlated methods such as TDCC and MCTDH. However, the quantitative claims in the abstract are stronger than the reported convergence data, and the fidelity monitor that drives basis adaptation is evaluated before masking and orthonormalization, so the error-control statement needs qualification.

major comments (3)
  1. [Abstract; Tables IV-V] The abstract states that grid calculations can be reproduced quantitatively using 30–100 Gaussians for intensities up to 4×10^14 W/cm^2. This range is not supported by the reported runs: Table V lists Mmax = 107 for LiH-DFT at 4×10^14 and Mmax = 107 and 118 for (LiH)2-DFT at 10^14 and 4×10^14, respectively; Table IV also contains Mmax = 90 and 107. Moreover, §V states that the (LiH)2 densities at 4×10^14 are not quite converged even at Mmax = 118, with ∫|Δρ|dx = 5.5×10^-1. The abstract should either be revised to reflect the actual basis sizes and the partial convergence, or additional calculations should be supplied that fit the claimed 30–100 range.
  2. [§IVA–C, §IVE] The quantity that triggers basis addition and replacement is the Rothe error r_{i+1} of Eqs. (27)–(28), evaluated before the masking function of Eq. (51) is applied and before the Löwdin orthonormalization described in §IVC. The error of the state that actually enters the next time step is therefore not the monitored quantity. Because the per-step Rothe error is the only fidelity monitor used in the adaptive procedure (§IVE), the statement that the basis is adapted according to the actual representation error is not fully established. The manuscript's own Table IV gives a case, (LiH)2 TDDFT at I0 = 10^14 W/cm^2 with εΔt = ∞, where the cumulative Rothe error decreases from 20.4 to 9.6 while the final density error increases from 0.740 to 0.769. The Discussion in §VI explains why rtot is not a global bound, but the local monitor is still the pre-mask, pre-orthonormalization quantity. I ask the authors to quantify, for a representative trajectory, how much the post-mask, post-orthonormalization error differs from the monitored error, or to state explicitly that the adaptation is heuristic rather than rigorously error-controlled.
  3. [Abstract; §IVA, §IVD] The abstract says 'removing the need for grids,' but the current implementation relies on grids in several essential ways: the initial ground-state orbitals are obtained from a grid calculation (§IID), the Rothe error is evaluated by quadrature on a grid (§IVA), and the width bounds [amin, amax] are imposed because the Gaussians must be numerically zero at the grid boundaries (§IVD), including a manual change of amin at t = 219 a.u. for the LiH-TDHF high-accuracy run. The paper appropriately labels the study as a proof of principle, but the unqualified abstract statement is stronger than what is demonstrated. The text should separate the grid-assisted implementation used here from the future analytical-integral implementation, and the abstract should be qualified accordingly.
minor comments (5)
  1. [Data Availability; Ref. 80] The Data Availability statement cites Ref. 80, but that Zenodo DOI corresponds to the authors' earlier paper on explicitly correlated Gaussian wave packets (Ref. 46), not to the present study. Please provide the correct DOI for the code and data supporting this manuscript.
  2. [§IVD; Tables IV-V] The 'frozen?' column in Tables IV and V is ambiguous, because even in the 'no' rows the nonlinear coefficients of the ground-state Gaussians are kept frozen; only the four additional Gaussians and all linear coefficients move. Please clarify the caption or the table notation.
  3. [§IIC, Eq. (53)] In Eq. (53), the line-search parameter δ is introduced without specifying whether it is a scalar or a vector, and the text does not state how its line search is performed. Please clarify.
  4. [§V] The sentence in §V referring to TDHF using '43 and 58 additional Gaussians' appears inconsistent with Table V if 'additional' is counted relative to the initial M0 (for LiH, Mmax = 66 gives 42 additional beyond the initial 24; for (LiH)2, Mmax = 95 gives 57 beyond 38). Please reconcile the counting or rephrase.
  5. [§IVA] The quadrature description for evaluating the Rothe error ('239 points' in (−17,17), trapezoidal with Δx = 0.4 outside) should state explicitly that this grid is different from the uniform Δx = 0.25 sinc-DVR grid used for the reference propagation, and comment on the effect of this mismatch on the reported comparisons.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: Rothe propagation is independently benchmarked against DVR grid references; self-citations are contextual only.

full rationale

The paper's central claim—that TDHF/TDDFT orbital propagation can be recast as a sequence of least-squares problems over thawed complex Gaussians—is not equivalent to its inputs. The Rothe functional r_{i+1} (Eqs. 27–28) is the Crank-Nicolson residual for the orbitals, and minimizing it is a bona fide reduced-basis approximation to the same linear system solved by the sinc-DVR reference; there is no fitted parameter that is later renamed as a prediction. The adaptive basis additions are driven by this residual, and although the residual is evaluated on a grid (Sec. IVA) and the paper concedes that r_tot is not a global error bound, this is a numerical-implementation limitation, not a circular reduction: the final observables (dipole moments, HHG spectra, densities) are compared directly to independent DVR calculations, and the paper reports cases in Tables IV/V where r_tot improves while the density error worsens, demonstrating that the monitor and the target are not conflated by construction. The citations to Refs. 44–46 are to the authors' own prior Rothe-method formulations, but the present paper gives the orbital working equations itself and validates them against external grid benchmarks; no uniqueness claim or unverified theorem is imported. Hence the derivation chain is self-contained at the tested level, and only minor, non-load-bearing self-citation is present.

Assumptions & free parameters 9 free parameters · 7 assumptions · 0 invented entities

The method introduces no new physical entities. The free parameters are numerical hyperparameters and basis-set sizes chosen by hand or by tolerance; the key load-bearing assumptions are the fidelity of the grid-based error metric and the expressivity of adaptive Gaussians for ionized dynamics.

free parameters (9)
  • Rothe error tolerance epsilon_delta_t = 0.1 to 30 depending on system, method, and intensity
    Adaptive basis enlargement threshold; chosen per run to balance accuracy and cost (Tables IV, V).
  • Penalty strength epsilon_p = 1e-4 (I0=1e14 W/cm2) or 1e-2 (I0=4e14 W/cm2)
    Penalty for near-linear dependence and small widths in basis refitting (Sec. IVF).
  • Overlap threshold smax = 0.99
    Triggers removal of Gaussians whose overlap exceeds this value (Sec. IVF).
  • Eigenvalue threshold lambda_min = 1e-10 or 1e-9 depending on intensity
    Detects global linear dependency in the overlap matrix (Sec. IVF).
  • Redundancy factor kappa = 1.1
    Gaussians whose removal increases the Rothe error by less than this factor are replaced (Sec. IVE).
  • Overlap penalty parameter x = 0.95 (0.7 for single-Gaussian step)
    Controls penalty onset for near-linear dependence in refitting objective (Sec. IVF).
  • Parameter update bounds s and q = s=0.1 (0.5 for new Gaussians), q=0.05 for a, q=0.1 for b, mu, p
    Restrict parameter changes per time step to stabilize optimization (Sec. IVD).
  • Width bounds amin, amax = amin=0.1 (0.04 at t=219 a.u. for LiH HF, I0=4e14 W/cm2), amax=2
    Grid-imposed limits on Gaussian widths; artifacts of the grid-based error evaluation (Sec. IVD).
  • Initial Gaussian count M0 = 20 for LiH, 34 for (LiH)2, plus 4 field Gaussians
    Chosen manually to represent ground-state orbitals; affects starting flexibility (Sec. VA).
assumptions (7)
  • standard math The TDHF and TDDFT orbital equations (Eq. 3) with the Fock or Kohn-Sham operator correctly describe the mean-field dynamics.
    Starting point of the method; standard theory from Refs. 13-17.
  • domain assumption The adiabatic 1D LDA functional (Eqs. 44-48, parameters from Ref. 49) is sufficiently accurate for the model systems.
    Adopted from the 1D DFT literature; its accuracy in strong fields for these models is assumed.
  • domain assumption The constant mean-field approximation (Eq. 21) and Crank-Nicolson propagator (Eq. 22) give a stable reference with global error O(Delta t).
    Defines the target propagation; first-order accuracy is accepted by the authors.
  • ad hoc to paper Adaptively added thawed complex Gaussians (Eq. 32) can represent the evolving orbitals, including ionized density, to the required accuracy with manageable basis sizes.
    Core expressivity premise; supported only by numerical evidence in this paper.
  • ad hoc to paper The grid-based Rothe error evaluation (Sec. IVA) with masking is a faithful proxy for the Hilbert-space representation error.
    Basis adaptation and convergence claims rest on this proxy; the method is grid-free only in the propagation step.
  • domain assumption Refitting masked orbitals as Gaussians (Sec. IVB) introduces negligible error in the physical region.
    Standard absorbing-boundary practice; not separately validated against unmasked propagation.
  • domain assumption Findings in 1D transfer to 3D molecules with similar Gaussian counts.
    Motivational extrapolation stated in Secs. I and VI; not demonstrated.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Time-dependent Gaussian basis sets for many-body systems using Rothe's method: A mean-field study." pith.science (2026). https://pith.science/paper/3ZU3576Y

@misc{pith2026250610701,
  author       = {Pith},
  title        = {Pith review of: Time-dependent Gaussian basis sets for many-body systems using Rothe's method: A mean-field study},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3ZU3576Y}},
  note         = {Machine review of arXiv:2506.10701}
}
abstract

A challenge in modeling time-dependent strong-field processes such as high-harmonic generation for many-body systems, is how to effectively represent the electronic continuum. We apply Rothe's method to the time-dependent Hartree-Fock (TDHF) and density functional theory (TDDFT) equations of motion for the orbitals, which reformulates them as an optimization problem. We show that thawed, complex-valued Gaussian basis sets can be propagated efficiently for these orbital-based approaches, removing the need for grids. In particular, we illustrate that qualitatively correct results can often be obtained by using just a few fully flexible Gaussians that describe the unbound dynamics for both TDHF and TDDFT. Grid calculations can be reproduced quantitatively using $30$--$100$ Gaussians for intensities up to $4\times10^{14}$ W/cm$^2$ for the one-dimensional molecular systems considered in this work.

Figures

Figures reproduced from arXiv: 2506.10701 by the authors.

Figure 1
Figure 1. FIG. 1: An illustration of the workflow for Rothe’s [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Time-dependent dipole moment (a) and HHG [PITH_FULL_IMAGE:figures/full_fig_p014_2.png] view at source ↗
Figure 4
Figure 4. FIG. 4: Time-dependent dipole moment (a) and HHG [PITH_FULL_IMAGE:figures/full_fig_p015_4.png] view at source ↗
Figures from the paper (4 more)
Figure 6
Figure 6. Figure 6: FIG. 6: Time-dependent dipole moment (a) and HHG [PITH_FULL_IMAGE:figures/full_fig_p016_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7: Time-dependent dipole moment (a) and HHG [PITH_FULL_IMAGE:figures/full_fig_p017_7.png]
Figure 9
Figure 9. Figure 9: FIG. 9: Time-dependent dipole moment (a) and HHG [PITH_FULL_IMAGE:figures/full_fig_p018_9.png]
Figure 10
Figure 10. Figure 10: FIG. 10: Electronic density [PITH_FULL_IMAGE:figures/full_fig_p019_10.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

79 extracted references · 31 canonical work pages

  1. [1]

    S.; McIntyre, I

    McPherson, A.; Gibson, G.; Jara, H.; Johann, U.; Luk, T. S.; McIntyre, I. A.; Boyer, K.; Rhodes, C. K. Studies of multiphoton production of vacuum-ultraviolet radiation in the rare gases. J. Opt. Soc. Am. B 1987, 4, 595--601, DOI: 10.1364/JOSAB.4.000595 http://dx.doi.org/10.1364/JOSAB.4.000595

    doi:10.1364/josab.4.000595 1987

  2. [2]

    F.; Lompre, L

    Ferray, M.; L'Huillier, A.; Li, X. F.; Lompre, L. A.; Mainfray, G.; Manus, C. Multiple-harmonic conversion of 1064 nm radiation in rare gases. J. Phys. B At. Mol. Opt. Phys. 1988, 21, L31, DOI: 10.1088/0953-4075/21/3/001 http://dx.doi.org/10.1088/0953-4075/21/3/001

    doi:10.1088/0953-4075/21/3/001 1988

  3. [3]

    Attosecond physics

    Krausz, F.; Ivanov, M. Attosecond physics. Rev. Mod. Phys. 2009, 81, 163--234, DOI: 10.1103/RevModPhys.81.163 http://dx.doi.org/10.1103/RevModPhys.81.163

    doi:10.1103/revmodphys.81.163 2009

  4. [4]

    B.; Krausz, F

    Corkum, P. B.; Krausz, F. Attosecond science. Nat. Phys. 2007, 3, 381--387, DOI: 10.1038/nphys620 http://dx.doi.org/10.1038/nphys620

    work page doi:10.1038/nphys620 2007

  5. [5]

    Attosecond imaging of molecular electronic wavepackets

    Haessler, S.; Caillat, J.; Boutu, W.; Giovanetti-Teixeira, C.; Ruchon, T.; Auguste, T.; Diveki, Z.; Breger, P.; Maquet, A.; Carr \'e , B.; Ta \"i eb, R.; Sali \`e res, P. Attosecond imaging of molecular electronic wavepackets. Nat. Phys. 2010, 6, 200--206, DOI: 10.1038/nphys1511 http://dx.doi.org/10.1038/nphys1511

    work page doi:10.1038/nphys1511 2010

  6. [6]

    Attosecond Electron Dynamics in Molecules

    Nisoli, M.; Decleva, P.; Calegari, F.; Palacios, A.; Mart \' n, F. Attosecond Electron Dynamics in Molecules . Chem. Rev. 2017, 117, 10760--10825, DOI: 10.1021/acs.chemrev.6b00453 http://dx.doi.org/10.1021/acs.chemrev.6b00453

    doi:10.1021/acs.chemrev.6b00453 2017

  7. [7]

    J.; Pabst, S.; Kamarchik, E.; Mazziotti, D

    Greenman, L.; Ho, P. J.; Pabst, S.; Kamarchik, E.; Mazziotti, D. A.; Santra, R. Implementation of the time-dependent configuration-interaction singles method for atomic strong-field processes. Phys. Rev. A 2010, 82, 023406, DOI: 10.1103/PhysRevA.82.023406 http://dx.doi.org/10.1103/PhysRevA.82.023406

    work page doi:10.1103/physreva.82.023406 2010

  8. [8]

    J.; McCurdy, C

    Haxton, D. J.; McCurdy, C. W. Two methods for restricted configuration spaces within the multiconfiguration time-dependent Hartree-Fock method. Phys. Rev. A 2015, 91, 012509, DOI: 10.1103/PhysRevA.91.012509 http://dx.doi.org/10.1103/PhysRevA.91.012509

    work page doi:10.1103/physreva.91.012509 2015

Show all 79 references
  1. [9]

    J.; Lawler, K

    Haxton, D. J.; Lawler, K. V.; McCurdy, C. W. Multiconfiguration time-dependent Hartree-Fock treatment of electronic and nuclear dynamics in diatomic molecules. Phys. Rev. A 2011, 83, 063416, DOI: 10.1103/PhysRevA.83.063416 http://dx.doi.org/10.1103/PhysRevA.83.063416

    2011 doi

  2. [10]

    Sato, T.; Pathak, H.; Orimo, Y.; Ishikawa, K. L. Communication: Time-dependent optimized coupled-cluster method for multielectron dynamics. J. Chem. Phys. 2018, 148

    2018

  3. [11]

    Ab initio quantum dynamics using coupled-cluster

    Kvaal, S. Ab initio quantum dynamics using coupled-cluster. J. Chem. Phys. 2012, 136

    2012

  4. [12]

    S.; Aurbakken, E.; Sch yen,

    Ofstad, B. S.; Aurbakken, E.; Sch yen, . S.; Kristiansen, H. E.; Kvaal, S.; Pedersen, T. B. Time-dependent coupled-cluster theory. WIREs Comput. Mol. Sci. 2023, 13, e1666, DOI: 10.1002/wcms.1666 http://dx.doi.org/10.1002/wcms.1666

    2023 doi

  5. [13]

    Runge, E.; Gross, E. K. U. Density-Functional Theory for Time-Dependent Systems. Phys. Rev. Lett. 1984, 52, 997--1000, DOI: 10.1103/PhysRevLett.52.997 http://dx.doi.org/10.1103/PhysRevLett.52.997

    1984 doi

  6. [14]

    Ullrich, C. A. Time-Dependent Density-Functional Theory: Concepts and Applications ; Oxford University Press: Oxford, 2011; DOI: 10.1093/acprof:oso/9780199563029.001.0001 http://dx.doi.org/10.1093/acprof:oso/9780199563029.001.0001

    2011

  7. [15]

    Single-Reference ab Initio Methods for the Calculation of Excited States of Large Molecules

    Dreuw, A.; Head-Gordon, M. Single-Reference ab Initio Methods for the Calculation of Excited States of Large Molecules. Chem. Rev. 2005, 105, 4009--4037, DOI: 10.1021/cr0505627 http://dx.doi.org/10.1021/cr0505627

    2005 doi

  8. [16]

    Maitra, N. T. Perspective: Fundamental aspects of time-dependent density functional theory. J. Chem. Phys. 2016, 144, 220901, DOI: 10.1063/1.4953039 http://dx.doi.org/10.1063/1.4953039

    2016 doi

  9. [17]

    Tancogne-Dejean, N. et al. Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems. J. Chem. Phys. 2020, 152, 124119, DOI: 10.1063/1.5142502 http://dx.doi.org/10.1063/1.5142502

    2020 doi

  10. [18]

    Molecular Electronic-Structure Theory ; Wiley: Chichester, 2000

    Helgaker, T.; J rgensen, P.; Olsen, J. Molecular Electronic-Structure Theory ; Wiley: Chichester, 2000

    2000

  11. [19]

    Universal Gaussian basis sets for an optimum representation of Rydberg and continuum wavefunctions

    Kaufmann, K.; Baumeister, W.; Jungen, M. Universal Gaussian basis sets for an optimum representation of Rydberg and continuum wavefunctions . J. Phys. B At. Mol. Opt. Phys. 1989, 22, 2223, DOI: 10.1088/0953-4075/22/14/007 http://dx.doi.org/10.1088/0953-4075/22/14/007

    1989 doi

  12. [20]

    The role of Rydberg and continuum levels in computing high harmonic generation spectra of the hydrogen atom using time-dependent configuration interaction

    Luppi, E.; Head-Gordon, M. The role of Rydberg and continuum levels in computing high harmonic generation spectra of the hydrogen atom using time-dependent configuration interaction . J. Chem. Phys. 2013, 139, 164121, DOI: 10.1063/1.4824482 http://dx.doi.org/10.1063/1.4824482

    2013 doi

  13. [21]

    P.; Lesiuk, M.; Przybytek, M.; Efimov, D

    Wo \'z niak, A. P.; Lesiuk, M.; Przybytek, M.; Efimov, D. K.; Prauzner-Bechcicki, J. S.; Mandrysz, M.; Ciappina, M.; Pisanty, E.; Zakrzewski, J.; Lewenstein, M.; Moszyński, R. A systematic construction of Gaussian basis sets for the description of laser field ionization and hi...

    2021 doi

  14. [22]

    Gaussian continuum basis functions for calculating high-harmonic generation spectra

    Coccia, E.; Mussard, B.; Labeye, M.; Caillat, J.; Ta \" eb, R.; Toulouse, J.; Luppi, E. Gaussian continuum basis functions for calculating high-harmonic generation spectra. Int. J. Quantum Chem. 2016, 116, 1120--1131, DOI: 10.1021/acs.jctc.8b00656 http://dx.doi.org/10.1021/acs...

    2016 doi

  15. [23]

    S.; Schlegel, H

    Durden, A. S.; Schlegel, H. B. Evaluation of Diffuse Basis Sets for Simulations of Strong Field Ionization Using Time-Dependent Configuration Interaction with a Complex Absorbing Potential . J. Phys. Chem. A 2025, 129, 3353--3367, DOI: 10.1021/acs.jpca.5c00195 http://dx.doi.or...

    2025 doi

  16. [24]

    Laser-induced electron dynamics including photoionization: A heuristic model within time-dependent configuration interaction theory

    Klinkusch, S.; Saalfrank, P.; Klamroth, T. Laser-induced electron dynamics including photoionization: A heuristic model within time-dependent configuration interaction theory . J. Chem. Phys. 2009, 131, 114304, DOI: 10.1063/1.3218847 http://dx.doi.org/10.1063/1.3218847

    2009 doi

  17. [25]

    Ab initio lifetime correction to scattering states for time-dependent electronic-structure calculations with incomplete basis sets

    Coccia, E.; Assaraf, R.; Luppi, E.; Toulouse, J. Ab initio lifetime correction to scattering states for time-dependent electronic-structure calculations with incomplete basis sets. J. Chem. Phys. 2017, 147, 014106, DOI: 10.1063/1.4991563 http://dx.doi.org/10.1063/1.4991563

    2017 doi

  18. [26]

    Burghardt, I.; Giri, K.; Worth, G. A. Multimode quantum dynamics using Gaussian wavepackets: The Gaussian-based multiconfiguration time-dependent Hartree (G-MCTDH) method applied to the absorption spectrum of pyrazine . J. Chem. Phys. 2008, 129, 174104, DOI: 10.1063/1.2996349 ...

    2008 doi

  19. [27]

    Quantum Dynamics Simulations Using Gaussian Wavepackets: The vMCG Method

    Richings, G.; Polyak, I.; Spinlove, K.; Worth, G.; Burghardt, I.; Lasorne, B. Quantum Dynamics Simulations Using Gaussian Wavepackets: The vMCG Method. Int. Rev. Phys. Chem. 2015, 34, 269--308, DOI: 10.1080/0144235X.2015.1051354 http://dx.doi.org/10.1080/0144235X.2015.1051354

    2015

  20. [28]

    A Strategy for Time Dependent Quantum Mechanical Calculations Using a Gaussian Wave Packet Representation of the Wave Function

    Sawada, S.-I.; Heather, R.; Jackson, B.; Metiu, H. A Strategy for Time Dependent Quantum Mechanical Calculations Using a Gaussian Wave Packet Representation of the Wave Function. J. Chem. Phys. 1985, 83, 3009--3027, DOI: 10.1063/1.449204 http://dx.doi.org/10.1063/1.449204

    1985 doi

  21. [29]

    Kay, K. G. The matrix singularity problem in the time-dependent variational method. Chem. Phys. 1989, 137, 165--175, DOI: 10.1016/0301-0104(89)87102-2 http://dx.doi.org/10.1016/0301-0104(89)87102-2

    1989 doi

  22. [30]

    H.; Byun, C

    Lee, M. H.; Byun, C. W.; Choi, N. N.; Kim, D. S. Solving Time -dependent Schrödinger Equation Using Gaussian Wave Packet Dynamics . J. Korean Phys. Soc. 2018, 73, 1269--1278, DOI: 10.3938/jkps.73.1269 http://dx.doi.org/10.3938/jkps.73.1269

    2018 doi

  23. [31]

    Kristiansen, H. E. Coupled cluster theory for electron dynamics. Ph.D.\ thesis, University of Oslo, 2022

    2022

  24. [32]

    Heller, E. J. Frozen Gaussians: A very simple semiclassical approximation . J. Chem. Phys. 1981, 75, 2923--2931, DOI: 10.1063/1.442382 http://dx.doi.org/10.1063/1.442382

    1981 doi

  25. [33]

    Van \' c ek, J.; Begu s i \'c , T. In Molecular Spectroscopy and Quantum Dynamics ; Marquardt, R., Quack, M., Eds.; Elsevier, 2020; pp 199--229, DOI: 10.1016/B978-0-12-817234-6.00011-8 http://dx.doi.org/10.1016/B978-0-12-817234-6.00011-8

    2020 doi

  26. [34]

    Regularized dynamical parametric approximation

    Feischl, M.; Lasser, C.; Lubich, C.; Nick, J. Regularized dynamical parametric approximation. 2024

    2024

  27. [35]

    V.; Child, M

    Shalashilin, D. V.; Child, M. S. Multidimensional Quantum Propagation with the Help of Coupled Coherent States. J. Chem. Phys. 2001, 115, 5367--5375, DOI: 10.1063/1.1394939 http://dx.doi.org/10.1063/1.1394939

    2001 doi

  28. [36]

    Kirrander, A.; Shalashilin, D. V. Quantum dynamics with fermion coupled coherent states: Theory and application to electron dynamics in laser fields. Phys. Rev. A 2011, 84, 033406, DOI: 10.1103/PhysRevA.84.033406 http://dx.doi.org/10.1103/PhysRevA.84.033406

    2011 doi

  29. [37]

    Symonds, C.; Wu, J.; Ronto, M.; Zagoya, C.; Figueira de Morisson Faria, C.; Shalashilin, D. V. Coupled-coherent-states approach for high-order harmonic generation. Phys. Rev. A 2015, 91, 023427, DOI: 10.1103/PhysRevA.91.023427 http://dx.doi.org/10.1103/PhysRevA.91.023427

    2015 doi

  30. [38]

    Sato, T.; Ishikawa, K. L. Time-dependent complete-active-space self-consistent-field method for multielectron dynamics in intense laser fields. Phys. Rev. A 2013, 88, 023402, DOI: 10.1103/PhysRevA.88.023402 http://dx.doi.org/10.1103/PhysRevA.88.023402

    2013 doi

  31. [39]

    L.; B r r ezinov\'a, I.; Lackner, F.; Nagele, S.; Burgd\"orfer, J

    Sato, T.; Ishikawa, K. L.; B r r ezinov\'a, I.; Lackner, F.; Nagele, S.; Burgd\"orfer, J. Time-dependent complete-active-space self-consistent-field method for atoms: Application to high-order harmonic generation. Phys. Rev. A 2016, 94, 023405, DOI: 10.1103/PhysRevA.94.023405 ...

    2016 doi

  32. [40]

    Sato, T.; Pathak, H.; Orimo, Y.; Ishikawa, K. L. Communication: Time-dependent optimized coupled-cluster method for multielectron dynamics. J. Chem. Phys. 2018, 148, 051101, DOI: 10.1063/1.5020633 http://dx.doi.org/10.1063/1.5020633

    2018 doi

  33. [41]

    Hybrid Gaussian-- B -spline basis for the electronic continuum: Photoionization of atomic hydrogen

    Marante, C.; Argenti, L.; Mart\' n, F. Hybrid Gaussian-- B -spline basis for the electronic continuum: Photoionization of atomic hydrogen. Phys. Rev. A 2014, 90, 012506, DOI: 10.1103/PhysRevA.90.012506 http://dx.doi.org/10.1103/PhysRevA.90.012506

    2014 doi

  34. [42]

    I.; Feist, J.; Noble, C

    Guan, X.; Zatsarinny, O.; Bartschat, K.; Schneider, B. I.; Feist, J.; Noble, C. J. General approach to few-cycle intense laser interactions with complex atoms. Phys. Rev. A 2007, 76, 053411, DOI: 10.1103/PhysRevA.76.053411 http://dx.doi.org/10.1103/PhysRevA.76.053411

    2007 doi

  35. [43]

    E.; Martín, F

    Bachau, H.; Cormier, E.; Decleva, P.; Hansen, J. E.; Martín, F. Applications of B-splines in atomic and molecular physics. Rep. Prog. Phys. 2001, 64, 1815, DOI: 10.1088/0034-4885/64/12/205 http://dx.doi.org/10.1088/0034-4885/64/12/205

    2001 doi

  36. [44]

    B.; Adamowicz, L

    Kvaal, S.; Lasser, C.; Pedersen, T. B.; Adamowicz, L. No need for a grid: Adaptive fully-flexible gaussians for the time-dependent Schr\"odinger equation . 2023

    2023

  37. [45]

    E.; Kristiansen, H

    Schrader, S. E.; Kristiansen, H. E.; Pedersen, T. B.; Kvaal, S. Time evolution as an optimization problem: The hydrogen atom in strong laser fields in a basis of time-dependent Gaussian wave packets . J. Chem. Phys. 2024, 161, 044105, DOI: 10.1063/5.0213576 http://dx.doi.org/1...

    2024 doi

  38. [46]

    E.; Pedersen, T

    Schrader, S. E.; Pedersen, T. B.; Kvaal, S. Multidimensional quantum dynamics with explicitly correlated Gaussian wave packets using Rothe’s method. J. Chem. Phys 2025, 162, 024109, DOI: 10.1063/5.0247732 http://dx.doi.org/10.1063/5.0247732

    2025 doi

  39. [47]

    P.; Adamowicz, L.; Pedersen, T

    Wo \'z niak, A. P.; Adamowicz, L.; Pedersen, T. B.; Kvaal, S. Rothe Time Propagation for Coupled Electronic and Rovibrational Quantum Dynamics. 2025; https://arxiv.org/abs/2503.09813

    2025 arXiv

  40. [48]

    Local-in-Time Error in Variational Quantum Dynamics

    Martinazzo, R.; Burghardt, I. Local-in-Time Error in Variational Quantum Dynamics. Phys. Rev. Lett. 2020, 124, 150601, DOI: 10.1103/PhysRevLett.124.150601 http://dx.doi.org/10.1103/PhysRevLett.124.150601

    2020 doi

  41. [49]

    I.; Casula, M.; Verstraete, M

    Helbig, N.; Fuks, J. I.; Casula, M.; Verstraete, M. J.; Marques, M. A. L.; Tokatly, I. V.; Rubio, A. Density functional theory beyond the linear regime: Validating an adiabatic local density approximation. Phys. Rev. A 2011, 83, 032503, DOI: 10.1103/PhysRevA.83.032503 http://d...

    2011 doi

  42. [50]

    O.; Stoudenmire, E

    Wagner, L. O.; Stoudenmire, E. M.; Burke, K.; White, S. R. Reference electronic structure calculations in one dimension. Phys. Chem. Chem. Phys. 2012, 14, 8581--8590, DOI: 10.1039/C2CP24118H http://dx.doi.org/10.1039/C2CP24118H

    2012 doi

  43. [51]

    Silva, J. W. O.; Vieira, D. Construction of Exchange-Correlation Potentials for Strongly Interacting One-Dimensional Systems. Braz. J. Phys. 2017, 47, 393--399, DOI: 10.1007/s13538-017-0508-x http://dx.doi.org/10.1007/s13538-017-0508-x

    2017 doi

  44. [52]

    E.; Stoudenmire, E

    Baker, T. E.; Stoudenmire, E. M.; Wagner, L. O.; Burke, K.; White, S. R. One-dimensional mimicking of electronic structure: The case for exponentials. Phys. Rev. B 2015, 91, 235141, DOI: 10.1103/PhysRevB.91.235141 http://dx.doi.org/10.1103/PhysRevB.91.235141

    2015 doi

  45. [53]

    V.; Rubio, A

    Helbig, N.; Tokatly, I. V.; Rubio, A. Exact Kohn–Sham potential of strongly correlated finite systems. J. Chem. Phys. 2009, 131, 224105, DOI: 10.1063/1.3271392 http://dx.doi.org/10.1063/1.3271392

    2009 doi

  46. [54]

    Szabo, A.; Ostlund, N. S. Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory ; Dover: New York, NY, USA, 1982

    1982

  47. [55]

    Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133--A1138, DOI: 10.1103/PhysRev.140.A1133 http://dx.doi.org/10.1103/PhysRev.140.A1133

    1965 doi

  48. [56]

    Dirac, P. A. M. Note on Exchange Phenomena in the Thomas Atom . Math. Proc. Camb. Philos. Soc. 1930, 26, 376--385, DOI: 10.1017/S0305004100016108 http://dx.doi.org/10.1017/S0305004100016108

    1930 doi

  49. [57]

    Gross, E. K. U.; Kohn, W. Local density-functional theory of frequency-dependent linear response. Phys. Rev. Lett. 1985, 55, 2850--2852, DOI: 10.1103/PhysRevLett.55.2850 http://dx.doi.org/10.1103/PhysRevLett.55.2850

    1985 doi

  50. [58]

    H.; Pereyra, V

    Golub, G. H.; Pereyra, V. The Differentiation of Pseudo-Inverses and Nonlinear Least Squares Problems Whose Variables Separate. SIAM J. Numer. Anal. 1973, 10, 413--432, DOI: 10.1137/0710036 http://dx.doi.org/10.1137/0710036

    1973 doi

  51. [59]

    P.; Rust, B

    O'Leary, D. P.; Rust, B. W. Variable projection for nonlinear least squares problems. Comput. Optim. Appl. 2013, 54, 579--593, DOI: 10.1007/s10589-012-9492-9 http://dx.doi.org/10.1007/s10589-012-9492-9

    2013 doi

  52. [60]

    Meyer, H.-D.; Gatti, F.; Worth, G. A. Multidimensional quantum dynamics: MCTDH theory and applications; John Wiley & Sons, 2009

    2009

  53. [61]

    Nocedal, J.; Wright, S. J. Numerical optimization; Springer: New York, 1999; DOI: 10.1007/b98874 http://dx.doi.org/10.1007/b98874

    1999 doi

  54. [62]

    On the equivalence of time-dependent variational principles

    Broeckhove, J.; Lathouwers, L.; Kesteloot, E.; Van Leuven, P. On the equivalence of time-dependent variational principles. Chem. Phys. Lett. 1988, 149, 547--550, DOI: 10.1016/0009-2614(88)80380-4 http://dx.doi.org/10.1016/0009-2614(88)80380-4

    1988 doi

  55. [63]

    Various variational approximations of quantum dynamics

    Lasser, C.; Su, C. Various variational approximations of quantum dynamics . J. Math. Phys. 2022, 63, 072107, DOI: 10.1063/5.0088265 http://dx.doi.org/10.1063/5.0088265

    2022 doi

  56. [64]

    Simple Theorems , Proofs , and Derivations in Quantum Chemistry ; Springer: New York, NY, USA, 2003; DOI: 10.1007/978-1-4757-6519-9 http://dx.doi.org/10.1007/978-1-4757-6519-9

    Mayer, I. Simple Theorems , Proofs , and Derivations in Quantum Chemistry ; Springer: New York, NY, USA, 2003; DOI: 10.1007/978-1-4757-6519-9 http://dx.doi.org/10.1007/978-1-4757-6519-9

    2003 doi

  57. [65]

    Thom, A. J. W.; Head-Gordon, M. Hartree–Fock solutions as a quasidiabatic basis for nonorthogonal configuration interaction . J. Chem. Phys. 2009, 131, 124113, DOI: 10.1063/1.3236841 http://dx.doi.org/10.1063/1.3236841

    2009 doi

  58. [66]

    Trigonometric pulse envelopes for laser-induced quantum dynamics

    Barth, I.; Lasser, C. Trigonometric pulse envelopes for laser-induced quantum dynamics. J. Phys. B At. Mol. Opt. Phys. 2009, 42, 235101, DOI: 10.1088/0953-4075/42/23/235101 http://dx.doi.org/10.1088/0953-4075/42/23/235101

    2009 doi

  59. [67]

    B.; Tukey, J

    Blackman, R. B.; Tukey, J. W. The measurement of power spectra from the point of view of communications engineering — Part I. Bell Syst. Tech. J. 1958, 37, 185--282, DOI: 10.1002/j.1538-7305.1958.tb03874.x http://dx.doi.org/10.1002/j.1538-7305.1958.tb03874.x

    1958

  60. [68]

    H.; Rubio, A

    De Giovannini, U.; Larsen, A. H.; Rubio, A. Modeling electron dynamics coupled to continuum states in finite volumes with absorbing boundaries. Eur. Phys. J. B 2015, 88, 56, DOI: 10.1140/epjb/e2015-50808-0 http://dx.doi.org/10.1140/epjb/e2015-50808-0

    2015 doi

  61. [69]

    On Löwdin's method of symmetric orthogonalization

    Mayer, I. On Löwdin's method of symmetric orthogonalization. Int. J. Quantum Chem. 2002, 90, 63--65, DOI: https://doi.org/10.1002/qua.981 http://dx.doi.org/https://doi.org/10.1002/qua.981

    2002 doi

  62. [70]

    Virtanen, P. et al. SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python . Nat. Methods 2020, 17, 261--272, DOI: 10.1038/s41592-019-0686-2 http://dx.doi.org/10.1038/s41592-019-0686-2

    2020 doi

  63. [71]

    L.; Adamowicz, L

    Bubin, S.; Pavanello, M.; Tung, W.-C.; Sharkey, K. L.; Adamowicz, L. Born–Oppenheimer and Non-Born–Oppenheimer, Atomic and Molecular Calculations with Explicitly Correlated Gaussians . Chem. Rev. 2013, 113, 36--79, DOI: 10.1021/cr200419d http://dx.doi.org/10.1021/cr200419d

    2013 doi

  64. [72]

    Discrete variable representations and sudden models in quantum scattering theory

    Lill, J.; Parker, G.; Light, J. Discrete variable representations and sudden models in quantum scattering theory. Chem. Phys. Lett. 1982, 89, 483--489, DOI: 10.1016/0009-2614(82)83051-0 http://dx.doi.org/10.1016/0009-2614(82)83051-0

    1982 doi

  65. [73]

    T.; Miller, W

    Colbert, D. T.; Miller, W. H. A novel discrete variable representation for quantum mechanical reactive scattering via the S-matrix Kohn method. J. Chem. Phys. 1992, 96, 1982--1991

    1992

  66. [74]

    Boyd, J. P. Chebyshev and Fourier Spectral Methods , 2nd ed.; Dover: New York, 2000

    2000

  67. [76]

    F.; Demmel, J.; Donato, J.; Dongarra, J.; Eijkhout, V.; Pozo, R.; Romine, C.; Van der Vorst, H

    Barrett, R.; Berry, M.; Chan, T. F.; Demmel, J.; Donato, J.; Dongarra, J.; Eijkhout, V.; Pozo, R.; Romine, C.; Van der Vorst, H. Templates for the solution of linear systems: building blocks for iterative methods; SIAM: Philadelphia, PA, USA, 1994; DOI: 10.1137/1.9781611971538...

    1994 doi

  68. [77]

    W.; Pfau, D.; Schaul, T.; Shillingford, B.; de Freitas, N

    Andrychowicz, M.; Denil, M.; Gomez, S.; Hoffman, M. W.; Pfau, D.; Schaul, T.; Shillingford, B.; de Freitas, N. Learning to learn by gradient descent by gradient descent. 2016; https://arxiv.org/abs/1606.04474

    2016 arXiv

  69. [78]

    Learning to Optimize : A Primer and A Benchmark

    Chen, T.; Chen, X.; Chen, W.; Heaton, H.; Liu, J.; Wang, Z.; Yin, W. Learning to Optimize : A Primer and A Benchmark . 2021; http://arxiv.org/abs/2103.12828, arXiv:2103.12828 [math]

    2021 arXiv

  70. [79]

    L2O- g^ : Learning to Optimize Parameterized Quantum Circuits with Fubini-Study Metric Tensor

    Huang, Y.-C.; Goan, H.-S. L2O- g^ : Learning to Optimize Parameterized Quantum Circuits with Fubini-Study Metric Tensor . 2024; http://arxiv.org/abs/2407.14761, arXiv:2407.14761 [quant-ph]

    2024 arXiv

  71. [80]

    Multidimensional quantum dynamics with explicitly correlated Gaussian wave packets using Rothe's method

    Schrader, S. E. Code and Data for "Multidimensional quantum dynamics with explicitly correlated Gaussian wave packets using Rothe's method" . 2025; https://doi.org/10.5281/zenodo.15648611 mcitethebibliography main.bib0000664000000000000000000023474615046310457011205 0ustar roo...

    2025

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.