{"id":"53c72a4a-8560-41d5-bd7c-db7407f2af8d","arxiv_id":"1908.02605","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A velocity-gauge perfectly matched layer in a finite-difference time-domain Schrödinger solver accurately simulates strong-field ionization with a simulation domain tens of times smaller than a plain grid.","lead":"This paper tests a numerical trick called perfectly matched layers for simulating laser-driven electron ionization in small one-dimensional model potentials. When implemented in the right gauge, the method shrinks the simulation box by about a factor of 70 and beats an older absorbing-boundary technique for short-range potentials.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Velocity-gauge PML matching is exact for spatially uniform A(t); the load-bearing weakness is the first-order ECS stencil used in the head-to-head comparison.","rationale":"The reader's weakest assumption, the formal PML justification for the time-dependent velocity-gauge interaction, is not the main threat. I verified the continuum mapping for the p.A(t) term: because A(t) has no spatial dependence, the complex-coordinate substitution transforms the first-order derivative term exactly, and in the short-range case the exterior potential vanishes. Thus the PML's reflectionless property does not rely on a questionable extension of the modal analysis; even if the paper should state this one-line derivation explicitly, the numerical claim is not resting on a formal gap.\n\nThe real soft spot is the ECS comparison. The abstract's comparative claim is broad, but the ECS implementation is handicapped by first-order matching-point stencils, so the reported several-orders advantage may reflect an order mismatch rather than an intrinsic superiority of PML over ECS. This is a fairness and interpretation issue, not a refutation of the PML domain-reduction result. The 1e-15 accuracy, the analytic polarizability check, and the adiabatic ionization-rate comparison independently support the central PML claim.\n\nFor these reasons the conditional verdict remains appropriate: the numerical method is sound, but the ECS benchmark and the 'several orders of magnitude' phrasing need tightening. A single re-running of the ECS comparison with higher-order stencils would settle whether the comparative claim survives.","tokens_in":12733,"tokens_out":18181,"duration_ms":211123,"concrete_test":"Recompute the ECS error curves in Figs. 3-4 with p=2 (and p=3) Lagrange finite-difference formulas from Eqs. A5-A6 at x = +/-R0, keeping dx = 0.01, dt = 0.001, R0 = 20, d = 40, E0 = 0.1, omega = 0.52; evaluate sigma(20) at t = 200 and compare with the PML velocity-gauge value. If the ECS error drops by several orders of magnitude toward the PML value near 1e-15, the comparative claim is an artifact of the first-order stencil; if it stays near the current value around 1e-9, the claim stands.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The strongest numerical claim survives the formal criticism in the reader's weakest_assumption. Since A(t) in Eq. (3) is spatially uniform, setting phi(x,t) = psi(xi(x),t) with xi from Eq. (12) gives, by the chain rule, -iA(t) d_xi psi = -iA(t) c(x) d_x phi and -1/2 d_xi^2 psi = -1/2 c(x) d_x c(x) d_x phi. For V=0 in the exterior (true for Eq. (16) with R0 >= b), the PML equation (14) is satisfied exactly by the analytic continuation of the continuum solution at every t; time dependence of A(t) is not an obstacle. The reported 1e-15 error is therefore a meaningful demonstration of the PML concept, not a numerical accident.\n\nThe genuinely load-bearing weakness is the advertised comparison 'PML outperforms ECS by several orders of magnitude.' The ECS discretization in Appendix A uses p=1 Lagrange interpolation, making the spatial stencils only O(dx) at x = +/-R0 (Eqs. A10-A11), while the PML uses uniformly second-order stencils (Eqs. A3-A4). Figures 3-4 thus compare a globally second-order scheme against a first-order scheme at the matching point. The paper's own comment that the gap 'might be due to the poor performance of ECS when implemented in finite-difference schemes' concedes the uncontrolled variable. A higher-order ECS stencil could reduce or remove the claimed several-orders advantage. Separately, the phrase 'several orders of magnitude' for the domain reduction is 8000/120 about 67 in the 1D demonstration, less than two orders; this is wording overreach, not a correctness flaw.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a finite-difference time-domain (FDTD) Crank-Nicolson solver for the one-dimensional time-dependent Schrödinger equation with perfectly matched layers (PMLs), targeting strong-field ionization of an electron bound by a short-range potential. The PML is implemented by transforming spatial derivatives while leaving the potential untransformed, and the scheme is tested in both velocity and length gauges. The paper compares absorption functions, measures the interior error against a converged large-domain reference, compares PML with exterior complex scaling (ECS) in both gauges, and validates observables against analytical results: the frequency-dependent polarizability, the adiabatic Airy-function ionization rate, an asymptotic low-field expression, and the PPT formula. The central claims are that velocity-gauge PML reproduces the exact wave function in the interior to near machine precision with a domain of about 120 a.u. instead of 8000 a.u. without absorbers, that PML outperforms ECS for short-range potentials when both are implemented in finite differences, and that ECS remains preferable for long-range potentials.","tokens_in":13181,"tokens_out":8817,"duration_ms":98418,"significance":"The core numerical demonstration is strong and, if it holds, makes a useful contribution: a simple FDTD-PML scheme for short-range potentials that gives interior wave functions accurate to ~1e-15 and reduces the computational domain by nearly two orders of magnitude in the 1D test case. The paper earns credit for benchmarking every central observable against independent analytical or numerical results: the polarizability against Eq. (17), the adiabatic rate against Eqs. (22)-(24), and the strong-field rate against the PPT expression (25). The formal concern raised by the time-dependent velocity-gauge term is not actually fatal: since A(t) is spatially uniform, the chain rule gives c(x)∂_x[ψ_ex(ξ(x),t)] = ∂_ξ ψ_ex(...), so the transformed velocity-gauge operator exactly matches the exterior Hamiltonian when V=0. The main weakness is the advertised PML-versus-ECS comparison, which is confounded by the different discretization orders used for the two methods; this issue is load-bearing for the paper's comparative claim. The PML method itself appears sound, but the superiority over ECS is not yet established on the evidence presented.","major_comments":[{"comment":"The head-to-head comparison that yields the claim that PML outperforms ECS by several orders of magnitude is not a like-for-like comparison of the two methods. The ECS finite-difference stencils at x = ±R0, given in Eqs. (A10)-(A11), are only O(Δx) at the scaling radius, whereas the PML stencils in Eqs. (A3)-(A4) are uniformly O(Δx^2). With Δx = 1e-2 a.u. used throughout, the ECS error may be dominated by the first-order matching point. The sentence in Section V stating that the gap 'might be due to the poor performance of ECS when implemented in finite-difference schemes' explicitly concedes this uncontrolled variable but does not remove it. The authors should either implement ECS with a higher-order stencil at the scaling radius or provide a convergence study in Δx for both methods at fixed absorption parameters before claiming that PML is intrinsically superior to ECS for short-range potentials.","section":"Appendix A, Eqs. (A10)-(A11); Figs. 3-4 and Section V"}],"minor_comments":[{"comment":"The abstract and Section V state that PMLs reduce the computational domain 'by several orders of magnitude.' In the 1D demonstration the reduction is from 8000 a.u. to 120 a.u., a factor of about 67, which is less than two orders of magnitude. Please rephrase to 'nearly two orders of magnitude' or give the actual ratio.","section":"Abstract and Section V"},{"comment":"The derivation of the PML transformation is presented for spatially and temporally invariant potentials, and its application to the velocity-gauge term p·A(t) is stated without proof. Because A(t) is spatially uniform, a short chain-rule argument shows that the transformation is exact on the exterior for V=0; adding this argument would remove the apparent gap and strengthen the paper's formal basis.","section":"Section IV, Eq. (14)"},{"comment":"The text refers to the 'square absorption function' in Fig. 2(b), but Eq. (15) lists four functions: the nearly singular function, y^2, y^3, and a tanh-based function. Please clarify explicitly which entry in Eq. (15) each panel of Fig. 2 corresponds to, to avoid ambiguity.","section":"Section V, Fig. 2"},{"comment":"When comparing PML and ECS errors as a function of time, both methods use the same absorption width d=40 a.u.; for the R0=10 case the total domain is 100 a.u., which is smaller than the PML/ECS width used in the main comparisons. This is fine, but the caption or text should state explicitly that the total domain changes with R0 rather than remaining fixed.","section":"Section V, Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the core PML method and its numerical validation are convincing, and the formal concern about the time-dependent velocity-gauge term is not a real obstacle. The main substantive issue is the ECS comparison, which uses a first-order ECS stencil at the scaling radius and therefore does not establish the claimed several-orders-of-magnitude superiority. This is fixable within the scope of a revision, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid methods paper with one over-claimed comparison. The core result—velocity-gauge PML in a second-order FDTD scheme reproduces the exact wave function to ~1e-15 for a short-range potential on a 120 a.u. domain instead of 8000 a.u.—holds up. The stress-test note resolves the formal worry about the time-dependent velocity-gauge term: since A(t) is spatially uniform, the PML coordinate transformation acts like a time-dependent change of variables, and the matching is exact for the continuum solution when the potential vanishes outside the interior. That is real evidence, not a numerical accident.\n\nWhat is new: applying PML to the TDSE for strong-field ionization in the velocity gauge, demonstrating a large domain reduction and giving clear warnings that length-gauge PML is unreliable. The benchmarks are honest and extensive: analytical polarizability, Airy-function ionization rates, the PPT formula, and a converged 10000 a.u. reference. The paper also documents known limitations for long-range potentials and the low-frequency divergence, which is useful guidance for practitioners.\n\nSoft spots: the head-to-head with ECS is not a fair test. Appendix A shows the ECS stencil is only first-order at the matching point (p=1 Lagrange interpolation), while the PML stencils are second-order everywhere. So the claim that PML outperforms ECS by several orders of magnitude may largely reflect a weak ECS implementation, not an inherent PML advantage. The paper even hints at this with the phrasing \"might be due to the poor performance of ECS when implemented in finite-difference schemes,\" but the abstract and conclusion state the comparison without that caveat. Also, \"several orders of magnitude\" for the domain reduction is loose—8000/120 is a factor of 67, not two orders. Minor, but worth tightening.\n\nBottom line: the central methodological claim is believable and well tested; the ECS comparison needs a higher-order or corrected ECS stencil before the superiority claim is solid. This deserves a serious referee, and with modest revision—fix the comparison and soften the domain-reduction wording—it would be a useful contribution to computational atomic and condensed-phase dynamics.","headline":"Solid PML-for-Schrodinger paper with a believable core result and an overstated ECS comparison; worth refereeing after the ECS comparison is fixed.","tokens_in":13627,"tokens_out":1750,"would_cite":true,"duration_ms":17307,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A perfectly matched layer in the velocity gauge reproduces the exact interior wave function for short-range strong-field ionization, cutting the simulation domain from 8000 to 120 atomic units.","keywords":["finite-difference time-domain","perfectly matched layer","time-dependent Schrödinger equation","strong-field ionization","exterior complex scaling","velocity gauge","delta potential","absorption function"],"falsifier":"Propagate the velocity-gauge FDTD-PML for the same $\\delta$-potential at a field amplitude outside the tested range (e.g., $E_0=0.5$ a.u. at $\\omega=0.52$ a.u., or a few-cycle pulse with a nonzero carrier-envelope phase), compare the interior wave function to a reference solution on a domain of several thousand atomic units, and check whether the interior error at $t=200$ a.u. remains near $10^{-15}$; a jump of several orders of magnitude would show that the PML's absorption of the time-dependent interaction is parameter-dependent.","tokens_in":2057,"feed_emoji":"⚛️","tokens_out":2421,"duration_ms":78664,"temperature":0.7,"pith_summary":"The paper tests whether a perfectly matched layer (PML) can replace huge simulation boxes when solving the time-dependent Schrödinger equation for strong-field ionization on a finite-difference grid. For an electron bound to a short-range potential and driven by a monochromatic laser field, it claims the velocity-gauge PML reproduces the exact wave function inside the physical region with errors down to about $10^{-15}$ at $t=200$ a.u., using a total domain of 120 a.u. instead of 8000 a.u. The same scheme is shown to outperform exterior complex scaling (ECS) by several orders of magnitude when both are implemented with the same finite-difference formulas, though ECS remains superior for potentials that extend into the absorbing layer. The paper also reports that the length-gauge PML introduces much larger interior errors, and that at very low frequencies the length-gauge PML is the only one of the four methods that stays numerically stable.","feed_headline":"Absorber shrinks quantum grid from 8000 to 120 atomic units","feed_subtitle":"Velocity-gauge PML beats exterior complex scaling in finite-difference time-domain ionization runs.","key_machinery":"The central object is the PML coordinate transformation $x \\to x + i\\sigma_0 \\int_{R_0}^{x} f(x') dx'$ outside the interior, acting only on spatial derivatives: $\\partial_x \\to c(x)\\partial_x$ with $c(x)=1/[1+i\\sigma_0 f(x)]$. Inside $|x|\\le R_0$, $c=1$ and the original Schrödinger equation is recovered; outside, the complex stretch turns outgoing waves into decaying exponentials while ideally causing no reflection at the interface, so a Dirichlet boundary deep in the layer is harmless. The choice of absorption function $f$ (square, cubic, near-singular, or tanh-based) controls how gently the absorption turns on, and the paper shows that slowly growing $f$ yields smaller reflection errors at the cost of slower convergence.","core_discovery":"The central claim is that a perfectly matched layer implemented in the velocity gauge within a second-order finite-difference time-domain scheme reproduces the exact wave function in the interior region for an electron bound to a short-range potential driven by a strong monochromatic field. Using the square absorption function $f(y)=y^2$ with absorption strength $\\sigma_0=0.001$, the error inside $|x|\\le 20$ a.u. reaches $\\sim 10^{-15}$ at $t=200$ a.u. with an absorbing layer of width $d=40$ a.u., i.e., a total domain of 120 a.u. instead of the 8000 a.u. needed with Dirichlet boundaries. The same comparison shows the velocity-gauge PML beats exterior complex scaling implemented with the same finite-difference formulas by several orders of magnitude. The paper further claims that the length-gauge PML introduces much larger interior errors, that slowly growing absorption functions are preferable, and that at very low frequencies the length-gauge PML is the only one of the four methods that remains numerically stable.","pith_inferences":["The formal gap between the time-independent derivation and the time-dependent $\\mathbf p\\cdot\\mathbf A(t)$ term suggests the velocity-gauge PML advantage may not survive for arbitrary pulse shapes; scanning field amplitudes, frequencies, and carrier-envelope phases against a large-domain reference would map where the $10^{-15}$ accuracy holds.","Because slowly growing absorption functions reduce reflection, a non-uniform grid in the absorption layer, analogous to what the paper notes for ECS, could make very wide PML layers cheap and might improve accuracy for long-range potentials where PML is currently inferior.","The low-frequency stability of only the length-gauge PML suggests a practical hybrid: use velocity-gauge PML at moderate and high frequencies and switch to length-gauge PML in the adiabatic regime, combining the best accuracy with stability.","Extending the scheme to two and three dimensions would open a direct path to the exciton-dissociation problems the paper motivates; the long-range potential results caution that the PML would need a tail-aware treatment or a hybrid with ECS."],"forward_implications":["Strong-field ionization of short-range potentials can be simulated accurately on domains two orders of magnitude smaller than with plain Dirichlet boundaries, reducing memory and computation time.","For short-range potentials, velocity-gauge PML is the preferred absorber over ECS within finite-difference schemes, since ECS errors are several orders of magnitude larger.","Length-gauge PML must be avoided at moderate and high frequencies, where it introduces large interior errors, but it becomes the only stable option at very low frequencies.","Absorption functions that grow slowly reduce reflection errors at the cost of needing wider layers, and a square absorption function already reaches $10^{-15}$ with $d\\approx 40$ a.u.","The scheme reproduces known analytic results for the polarizability and for adiabatic ionization rates of the $\\delta$-potential, indicating that physical observables remain reliable even though the wave function is absorbed."],"supporting_citations":[{"why":"Berenger's original perfectly matched layer construction for Maxwell's equations, which the paper adapts to the Schrödinger equation.","marker":"[35]"},{"why":"Nissen and Kreiss's optimized PML for the Schrödinger equation with time-independent potentials, supplying the PML-Schrödinger formulation used here.","marker":"[43]"},{"why":"Scrinzi's error measure (Eq. 7), which defines the accuracy metric used throughout the paper.","marker":"[34]"},{"why":"He, Ruiz, and Becker's ECS treatment of the TDSE in laser fields, the ECS implementation the paper compares against and the source of the LG/VG exponential-behavior analysis.","marker":"[31]"},{"why":"McCurdy et al.'s review of ECS in finite-difference schemes, cited for the poor performance of ECS when implemented with finite differences.","marker":"[59]"},{"why":"Postma's analytic polarizability of the delta potential, the target for the weak-field verification.","marker":"[51]"},{"why":"Fernández and Castro's Airy-function condition for the DC ionization rate, the basis for the adiabatic-rate comparison.","marker":"[49]"},{"why":"Perelomov, Popov, and Terent'ev's adiabatic ionization rate, one of the analytic limits compared at low frequency.","marker":"[14]"},{"why":"Scharf et al.'s three-photon ionization line shape, used as a comparison for the frequency-dependent ionization rate.","marker":"[48]"}],"fun_headline_variants":["PML reduces quantum grid from 8000 to 120 a.u.","Velocity-gauge PML beats exterior complex scaling in FDTD","Absorbing layers cut ionization simulation grid 66-fold","PML outperforms exterior scaling for short-range potentials"],"cache_read_input_tokens":15616,"weakest_assumption_plain":"The PML is derived from a modal analysis for a potential that is constant in space and time, but the paper applies the same coordinate stretch to the explicitly time-dependent velocity-gauge interaction $\\mathbf p\\cdot\\mathbf A(t)$, and the perfect-absorption property for that term is not formally justified—only numerical agreement with a large-domain reference is offered as evidence.","fun_headline_variants_meta":{"raw":{"variants":["PML reduces quantum grid from 8000 to 120 a.u.","Velocity-gauge PML beats exterior complex scaling in FDTD","Absorbing layers cut ionization simulation grid 66-fold","PML outperforms exterior scaling for short-range potentials"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001232,"raw_usage":{"total_tokens":5060,"prompt_tokens":946,"completion_tokens":4114,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":562,"completion_tokens_details":{"reasoning_tokens":4043}},"tokens_in":562,"tokens_out":4114,"duration_ms":29744,"temperature":1.0,"reasoning_tokens":4043,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:39:39.011585+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Propagate the velocity-gauge FDTD-PML for the same $\\delta$-potential at a field amplitude outside the tested range (e.g., $E_0=0.5$ a.u. at $\\omega=0.52$ a.u., or a few-cycle pulse with a nonzero carrier-envelope phase), compare the interior wave function to a reference solution on a domain of several thousand atomic units, and check whether the interior error at $t=200$ a.u. remains near $10^{-15}$; a jump of several orders of magnitude would show that the PML's absorption of the time-dependent interaction is parameter-dependent.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Berenger's original perfectly matched layer construction for Maxwell's equations, which the paper adapts to the Schrödinger equation."},{"cited_title":"Zuo and Z","cited_arxiv_id":null,"evidence_quote":"Nissen and Kreiss's optimized PML for the Schrödinger equation with time-independent potentials, supplying the PML-Schrödinger formulation used here."},{"cited_title":"Scrinzi, Phys","cited_arxiv_id":null,"evidence_quote":"Scrinzi's error measure (Eq. 7), which defines the accuracy metric used throughout the paper."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"He, Ruiz, and Becker's ECS treatment of the TDSE in laser fields, the ECS implementation the paper compares against and the source of the LG/VG exponential-behavior analysis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"McCurdy et al.'s review of ECS in finite-difference schemes, cited for the poor performance of ECS when implemented with finite differences."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Postma's analytic polarizability of the delta potential, the target for the weak-field verification."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Fernández and Castro's Airy-function condition for the DC ionization rate, the basis for the adiabatic-rate comparison."},{"cited_title":"Perelomov, V","cited_arxiv_id":null,"evidence_quote":"Perelomov, Popov, and Terent'ev's adiabatic ionization rate, one of the analytic limits compared at low frequency."},{"cited_title":"Scharf, K","cited_arxiv_id":null,"evidence_quote":"Scharf et al.'s three-photon ionization line shape, used as a comparison for the frequency-dependent ionization rate."}],"review_version":1}