{"id":"1dfdd47c-5057-43fc-8738-558396ff05fc","arxiv_id":"2502.07569","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A time-splitting multiscale finite element method for the random-potential semiclassical cubic Schrödinger equation is proven and tested with second-order spatial/temporal and near-first-order random-space convergence.","lead":"This paper introduces faster time-splitting finite element solvers for the semiclassical nonlinear Schrödinger equation with random, oscillatory potentials, using multiscale coarse meshes. If the proven error bounds hold, simulations of wave localization and delocalization in disordered media become substantially cheaper.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 4.3 inherits exponential localization from prior work without proof that it holds uniformly over random samples or with the modified constraint; the MsFEM spatial error bound lacks support in the stochastic setting.","rationale":"The paper's main deliverable is Theorem 4.3, an a priori bound for the full discretization in space, time, and randomness. Every spatial contribution to that bound flows through Lemma 4.3, which is an application of the MsFEM approximation theory summarized in Propositions 3.1 and 3.2. The paper proves Proposition 3.1 in a short argument, and that part is straightforward given the first-order optimality conditions of (3.14)-(3.15). However, Proposition 3.2 is not proved; it is quoted from [52], and the present work changes the optimization constraints via lambda(H) = (1, phi_H^p) and moves to sample-dependent potentials. For the theorem to hold, the decay constant beta must be uniformly bounded below 1 across the random samples, since the expectation in (4.23) averages over omega. Assumption 2.1 as stated controls only the L-infinity norm of each sample; it gives no direct control of beta(omega). This is the precise point where the analysis is thinnest. The numerical section does not test localization directly, and the random-potential experiments (e.g., Section 5.3.1 with beta = 0, m = 5, epsilon = 1/8, H approx 0.063) appear to have ||v||_infinity approx 8.66, exceeding epsilon^2/H^2 approx 3.96, so at face value they lie outside Assumption 2.1. I therefore agree with the reader that the key weakness is the imported MsFEM theory. Because this is exactly the concern that motivates the CONDITIONAL verdict, my stress-test does not change that verdict; it only sharpens the required condition: verify numerically that (3.22) holds uniformly over random samples, and state the H-epsilon resolution condition imported from [52]. Other parts of the paper, including the temporal splitting analysis (Theorem 4.1), the KL-truncation error (Lemma 2.2), and the qMC quadrature bound (Lemma 4.5), are standard and internally consistent, and the small-Delta-t convergence tables are consistent with second order, so there is no reason to move to REJECT.","tokens_in":29180,"tokens_out":13518,"duration_ms":115657,"concrete_test":"Take the potential (5.2) with the parameters of Section 5.3.1 (sigma = 1, beta = 0, m = 5, epsilon = 1/8, H = 6h). For at least 100 qMC samples, solve (3.14)-(3.15) for all NH coarse nodes on the same fine mesh as the paper. For each basis function compute rho_l(p, omega) = ||grad phi_p||_{L^2(D \\ D_l)} / ||grad phi_p|| for l = 1..4. If max_{p, omega} rho_l decays like beta^l with beta < 1 uniformly and, in addition, ||psi_H^n - psi_ref|| follows the H^2/epsilon^3 trend when H is varied, the imported Proposition 3.2 is credible in the stochastic setting; otherwise Theorem 4.3 lacks a verified foundation. Also record whether the Assumption 2.1 bound ||v||_infinity <= C epsilon^2/H^2 holds for these samples.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is Lemma 4.3 (Eq. 4.18): it produces the H^2/epsilon^3 spatial term used in Theorem 4.3 (Eq. 4.23). Its proof reduces the MsFEM error to Proposition 3.2's exponential decay (3.22) and the split Vh = Vms ⊕ Wh (Prop. 3.1). Prop. 3.2 is quoted from [52] for a deterministic setting, while here the basis is redefined by (3.14)-(3.15) with lambda(H) = (1, phi_H^p); the rescaling argument changes the minimizers and the transfer of [52] is not shown. Moreover, for random potentials Assumption 2.1 bounds only the sample-wise L-infinity norm; it does not imply a uniform-in-omega decay rate beta < 1 in (3.22). Losing uniformity in beta would make the constants in Lemma 4.3 and the expectation in (4.23) uncontrolled. The numerical random experiments (Section 5.3.1 and 5.3.3 with sigma = 1, beta = 0, m = 5, epsilon = 1/8, H approx 0.063) have ||v||_infinity approx 8.66 while epsilon^2/H^2 approx 3.96, so they may even fall outside Assumption 2.1. Thus the central stochastic MsFEM claim is not independently substantiated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes two time-splitting finite element methods (SI and SII) for the semiclassical cubic nonlinear Schrödinger equation with random potentials, and combines them with a multiscale finite element method (MsFEM) whose basis functions are constructed from equality-constrained optimization problems. The main theoretical result is Theorem 4.3, which gives a root-mean-square error bound for the expectation of a linear functional of the solution in terms of H^2/ε^3, Δt^2/ε^4, m^{-χ}/ε, and C_{γ,m}N^{-r}. The paper also presents numerical experiments in 1D and 2D that show second-order spatial and temporal convergence and a near-first-order qMC convergence rate, and it describes a POD-based multiscale reduced basis method in Appendix A.","tokens_in":29481,"tokens_out":11213,"duration_ms":102301,"significance":"If the central error analysis is correct, this would be a useful contribution: it extends MsFEM error estimates to the semiclassical NLSE with random potentials and provides a rigorous rate in the random sampling dimension. The paper has several strengths: it proposes a concrete basis construction via optimization, it gives explicit convergence rates that are stated in a checkable form, and it includes numerical experiments with clean second-order rates and a comparison of MC and qMC. The POD reduced-basis appendix also addresses a practical bottleneck in constructing basis functions for many random samples. However, the significance is currently limited by a mismatch between the theory and the numerics (linear functionals vs. quadratic density) and by gaps in the proof of the load-bearing MsFEM spatial error estimate, as detailed below.","major_comments":[{"comment":"Theorem 4.3 and Lemma 4.5 are stated only for continuous linear functionals G, yet all random-potential experiments in §5.3 compute the expected density E(|ψ|^2), which is a quadratic functional. The error analysis does not cover this quantity, and an extension to quadratic (or at least Lipschitz) functionals, with control on ∥ψ∥_{H^2}, would be needed before the numerical rates for E(|ψ|^2) can be cited as evidence for the theorem.","section":"§5.3, Theorem 4.3"},{"comment":"The SI scheme implements the nonlinear substep as pointwise multiplication of nodal coefficients by exp(-iλΔt/(2ε)|U^n|^2). This is not the Galerkin discretization of the nonlinear flow in the FEM space and introduces an interpolation error that is not estimated. The analysis in Lemma 4.2 and Theorem 4.1 concerns the continuous splitting operator, while the fully discrete scheme in (3.10) is what is actually used; the notation L in (4.9) conflates the continuous and discrete operators. Consequently, the claimed O(H^2/ε^3) spatial error is not proven for the implemented scheme.","section":"§3.1, Eq. (3.10)"},{"comment":"The proof first quotes a bound O(H^2/ε^2) from [52] for the linear Schrödinger step and then concludes (4.18) with O(H^2/ε^3); the additional factor 1/ε is not derived. The lemma is stated for the full NLSE, but the cubic nonlinear term is only treated through a stability factor, and the MsFEM error introduced by the nonlinear substep is not estimated. As a result, the central spatial error term used in Theorem 4.3 is not verified by the argument given.","section":"§4.2.1, Lemma 4.3"},{"comment":"The exponential decay estimate (3.22) is imported from the deterministic analysis [52] and is asserted for the present rescaled constraints and for each random potential sample, but the paper does not prove that the decay constant β is uniform in ω. Assumption 2.1 is not stated as a hypothesis of Theorem 4.3, and the numerical test in §5.3.1 appears to violate it: with σ=1, β=0, m=5, ε=1/8, and H≈0.063, one has ∥v∥∞≈8.66 while ε^2/H^2≈3.96. A uniform-in-ω version of Proposition 3.2, or at least a numerical verification of Assumption 2.1 for the experiments, is needed to support the stochastic MsFEM claim.","section":"§4.2, Prop. 3.2 and Theorem 4.3"}],"minor_comments":[{"comment":"The expression '∥Lnψin − Lnψin∥' should read '∥L^n ψ0_h − L^n ψin∥' to match the error split in (4.9).","section":"§4.1, proof of Theorem 4.2"},{"comment":"In the covariance kernel, the notation '|xi − yj|^2' appears to have an index error: the second index j should presumably be i in the sum over dimensions.","section":"§2, Eq. (2.7)"},{"comment":"The measured convergence slopes for the MC and qMC curves are not reported; stating the numerical slopes would substantiate the claim of 'almost first-order' convergence in the random space.","section":"§5.3.2, Fig. 7"},{"comment":"The proof of Proposition 3.1 is only a sketch; a fuller derivation or a precise quotation of [52, Lemma 3.2] including the required resolution condition would improve verifiability.","section":"§3.2, Prop. 3.1"},{"comment":"The formula for ϱ(α) appears to have a parenthesis imbalance; please check and correct the expression.","section":"§4.2, Lemma 4.5, Eq. (4.21)"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses an important problem and contains promising numerical results, but the proof of the main theorem has gaps related to the discrete nonlinear substep and to the transfer of localization estimates to the random setting. I would be willing to review a revised version if the authors close these gaps or restate the claims to match the proven results. The manuscript may be suitable for the journal after such revisions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a real extension of MsFEM to the semiclassical cubic NLSE with random potentials, with a clean second-order scheme and a genuine convergence theorem. The main thing to know before reading: the theory is for linear observables while the numerics use the quadratic mass density, and the spatial error analysis leans on an imported exponential-decay result. Both are fixable, but they matter for how much weight you put on the claims.\n\nWhat is genuinely new and good: the mesh-dependent normalization λ(H) = (1, φ^H_p) is a neat fix for the spurious scale that otherwise appears in the cubic term, and the authors show it does not change the basis space. The SI variant with eigendecomposition handles discontinuous potentials better than SII, which is a useful observation. The convergence analysis is an honest attempt: Theorem 4.3 gives a root-mean-square error with terms H²/ε³, Δt²/ε⁴, m^{-χ}/ε, and N^{-r}, and the numerical tables show the advertised second-order rates and faster qMC convergence. The paper also proposes a POD-based reduced basis method with timing comparisons, which is practically relevant.\n\nSoft spots, in proportion: First, Theorem 4.3 is stated for continuous linear functionals G, but the numerical validation in Section 5.3 measures E(|ψ|²), a quadratic functional. The proof splits the error using linearity of G; that step does not directly apply to density, and the paper does not supply an L^p error estimate to bridge it. This is not a fatal flaw, but it means the numerics do not directly verify the theorem as stated. Second, Proposition 3.2, the exponential decay of multiscale basis functions, is imported from the authors' prior work [52]. The paper gives only a one-sentence transfer argument. The rescaling via λ(H) should not destroy decay, since it just scales the minimizers, but the uniformity of β over random samples is not stated. A referee should ask for that to be made explicit. Third, the numerical examples with σ=1, β=0, m=5 may violate Assumption 2.1 (‖v‖∞ can reach about 8.7 while ε²/H² is around 4). That does not invalidate the method, but the paper should acknowledge or adjust the experiments. Fourth, no code or data is provided, which limits reproducibility.\n\nWho is this for? Numerical analysts working on multiscale finite elements, semiclassical Schrödinger equations, or uncertainty quantification. They will find the method and rates useful. It deserves a serious referee; I would send it to review with a request to address the observable mismatch, the transfer of Proposition 3.2, and the assumption check.","headline":"A credible extension of MsFEM to the semiclassical cubic NLSE with random potentials, but the theory covers linear observables while the numerics use quadratic density, and the spatial error analysis leans on an imported localization result.","tokens_in":30014,"tokens_out":3909,"would_cite":true,"duration_ms":36070,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q55","65M60","81Q05","47H40"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves an $L^2$ error bound for a time-splitting multiscale finite element method for the semiclassical nonlinear Schrödinger equation with random potentials, giving second-order accuracy in space and time and almost…","keywords":["semiclassical nonlinear Schrödinger equation","random potentials","multiscale finite element method","time-splitting methods","quasi-Monte Carlo","Karhunen-Loève expansion","error estimate","Anderson localization"],"falsifier":"Compute, for a single random sample satisfying Assumption 2.1, the multiscale basis from (3.14)--(3.15) and measure two things: the decay of $\\|\\nabla\\phi_p\\|_{L^2(D\\setminus D_\\ell)}$ as $\\ell$ grows, and the $L^2$ error of the projection of a known smooth function onto $V_{\\mathrm{ms}}$ as $H\\to0$. If the decay is not exponential with ratio $\\beta<1$ or the error does not scale like $H^2/\\varepsilon^3$, then the propositions behind Theorem 4.3 fail for that sample.","tokens_in":28985,"feed_emoji":"⚛️","tokens_out":14596,"duration_ms":114909,"temperature":0.7,"pith_summary":"The paper proposes two Strang time-splitting finite element methods for the semiclassical cubic nonlinear Schrödinger equation with a random potential, and couples them to multiscale finite element spaces to keep the spatial degrees of freedom low. Multiscale basis functions are generated by solving constrained energy-minimization problems, with a rescaling constant that prevents the nonlinear term from introducing a spurious mesh-size dependence. The central result is a root-mean-square error bound for the sampled expectation of any linear functional of the wave function: $C(H^2/\\varepsilon^3+\\Delta t^2/\\varepsilon^4+m^{-\\chi}/\\varepsilon+C_{\\gamma,m}N^{-r})$, which the paper interprets as second-order accuracy in space and time and almost first-order accuracy in the number of quasi-Monte Carlo samples. The paper also shows numerically that the scheme retains its order with discontinuous potentials and that it reproduces Anderson-style localization for the linear equation and delocalization as the nonlinearity strengthens.","feed_headline":"Random-potential Schrödinger solved on coarse grids at second order","feed_subtitle":"MsFEM plus Strang splitting and qMC keeps the error bound while cutting degrees of freedom.","key_machinery":"The multiscale basis functions are the central mechanism. They are obtained from the equality-constrained quadratic program $\\min a(\\phi,\\phi)$ subject to $\\int_D \\phi\\,\\varphi_q^H\\,dx=\\lambda(H)\\delta_{pq}$, where $a(\\phi,\\phi)=\\frac{\\epsilon^2}{2}\\|\\nabla\\phi\\|^2+(v\\phi,\\phi)$ and $\\lambda(H)=(1,\\varphi_q^H)$. This normalization is what suppresses the mesh-dependent scale in the nonlinear term. A weighted quasi-interpolation operator $I_H$ produces the splitting $V_h=V_{\\mathrm{ms}}\\oplus W_h$, and an iterative decay argument imported from earlier work gives the exponential localization of the basis. In time, Strang splitting with $L_1$ solved by eigendecomposition (SI) or by a Crank--Nicolson update (SII) avoids nonlinear iterations. In the random dimension, the truncated Karhunen--Loève expansion parameterizes the potential and randomly shifted lattice rules provide the sampling error term $C_{\\gamma,m}N^{-r}$.","core_discovery":"The paper's claim is that a time-splitting multiscale finite element method can solve the semiclassical nonlinear Schrödinger equation with random potentials at the rate stated in Theorem 4.3: if $\\psi_{\\mathrm{in}}\\in H^4$ and the random potential satisfies Assumptions 2.1 and 2.2, the root-mean-square error between $E[G(\\psi^\\varepsilon(t_n))]$ and its quasi-Monte Carlo approximation $Q_{m,N}[G(\\psi^{\\varepsilon,n}_{H,m})]$ is bounded by $C(H^2/\\varepsilon^3+\\Delta t^2/\\varepsilon^4+m^{-\\chi}/\\varepsilon+C_{\\gamma,m}N^{-r})$. The spatial component of this bound comes from the multiscale basis construction (3.14)--(3.15): the basis functions minimize the energy with an orthogonality constraint whose scaling $\\lambda(H)=(1,\\varphi_q^H)$ cancels a mesh-dependent factor that would otherwise break the cubic nonlinearity. The paper further claims that the first splitting variant, which solves the linear part by eigendecomposition, keeps second-order convergence even for discontinuous potentials, while the classical variant requires smoother potentials and time steps of order $\\epsilon$.","pith_inferences":["The scale-cancelling choice $\\lambda(H)=(1,\\varphi_q^H)$ is specific to the cubic term; for other nonlinearities such as a quintic term, the same optimization-constraint mechanism would likely need a different normalization, so checking whether the advertised rate survives is a direct test of the mechanism.","The POD reduced-basis construction of Appendix A is demonstrated only as a timing improvement in 1D and 2D; an immediate question is whether the same offline/online decomposition preserves the $H^2/\\varepsilon^3$ spatial error in 3D, where the constants in the regularity bounds are largest.","The analysis assumes the Karhunen--Loève eigenvalues decay polynomially and the potential is almost surely bounded; if the potential is lognormal or unbounded, the $m^{-\\chi}$ and $N^{-r}$ terms would need re-derivation, and the numerical evidence here does not cover that case.","The finding that nonlinearity destroys Anderson localization is consistent with earlier studies; a natural next step would be to compare the delocalization threshold in $\\lambda$ against predictions from lattice or random-matrix models, but the paper does not do that."],"forward_implications":["A coarse multiscale mesh with mesh size $H$ gives the same order of spatial accuracy as a fine finite element mesh, so solving on $V_{\\mathrm{ms}}$ cuts the number of unknowns while retaining the $\\Delta t^2$ temporal order.","The MsFEM spatial error depends on $\\varepsilon$ only through $H^2/\\varepsilon^3$ and needs only $\\partial_t\\psi^\\varepsilon\\in L^2$, whereas the classical FEM spatial estimate needs $H^2$ with $\\partial_t\\psi^\\varepsilon\\in H^2$; the authors use this to say the MsFEM is better adapted to the semiclassical regime.","With the SI splitting, the scheme is second-order in space even for discontinuous potentials, so the method can be used for rough random media.","In the random dimension, the quasi-Monte Carlo rate is $N^{-r}$ for any $r=1-\\delta$ with $0<\\delta<1/2$; combined with the KL truncation error $m^{-\\chi}/\\varepsilon$, this gives a full uncertainty-quantification estimate.","The numerical experiments predict a physical transition: mass density localizes for the linear equation and delocalizes when $\\lambda$ is large, tracked by the moment $A(t)$. If these rates hold, the method is a practical tool for long-time wave-in-random-media simulations."],"supporting_citations":[{"why":"supplies the space decomposition $V_h=V_{\\mathrm{ms}}\\oplus W_h$ and the exponential localization of the multiscale basis that Lemma 4.3 and Theorem 4.3 inherit.","marker":"[52]"},{"why":"introduces the multiscale reduced-basis optimization construction for the Schrödinger equation with multiscale and random potentials and the eigendecomposition-based time stepping.","marker":"[14]"},{"why":"provides the Strang splitting stability and order estimates used in Theorems 4.1 and 4.2.","marker":"[10]"},{"why":"gives the time-splitting spectral discretization in the semiclassical regime and the $\\Delta t=O(\\epsilon)$, $h=O(\\epsilon)$ conditions behind SII.","marker":"[8]"},{"why":"supplies the weighted quasi-interpolation estimate $\\|f-I_H f\\|\\le H^2\\|f\\|_{H^2}$ used in the spatial decomposition.","marker":"[27]"},{"why":"provides the KL-truncation and quasi-Monte Carlo error analysis for the NLSE with random potentials that Lemma 4.6 and Theorem 4.3 adapt.","marker":"[53]"},{"why":"gives the randomly shifted lattice rule and component-by-component construction with the $N^{-r}$ sampling error used in Lemma 4.5.","marker":"[17]"},{"why":"supplies the iterative decay argument behind the exponential localization of the multiscale basis functions.","marker":"[32]"}],"fun_headline_variants":["MsFEM solves random Schrödinger at second order on coarse grids","Time-splitting MsFEM: second order for random-potential Schrödinger","Second-order MsFEM cuts unknowns for random nonlinear Schrödinger","Random Schrödinger at second order with reduced degrees of freedom","Multiscale FEM: second-order accuracy for random Schrödinger"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The spatial half of the headline error bound assumes, from earlier work, that every optimized multiscale basis function stays exponentially localized and that the finite element space splits as $V_h=V_{\\mathrm{ms}}\\oplus W_h$ for every random potential sample with the assumed boundedness; the paper cites this instead of re-deriving or testing it.","fun_headline_variants_meta":{"raw":{"variants":["MsFEM solves random Schrödinger at second order on coarse grids","Time-splitting MsFEM: second order for random-potential Schrödinger","Second-order MsFEM cuts unknowns for random nonlinear Schrödinger","Random Schrödinger at second order with reduced degrees of freedom","Multiscale FEM: second-order accuracy for random Schrödinger"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001518,"raw_usage":{"total_tokens":6094,"prompt_tokens":970,"completion_tokens":5124,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":586,"completion_tokens_details":{"reasoning_tokens":5034}},"tokens_in":586,"tokens_out":5124,"duration_ms":30671,"temperature":1.0,"reasoning_tokens":5034,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T12:16:58.806339+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute, for a single random sample satisfying Assumption 2.1, the multiscale basis from (3.14)--(3.15) and measure two things: the decay of $\\|\\nabla\\phi_p\\|_{L^2(D\\setminus D_\\ell)}$ as $\\ell$ grows, and the $L^2$ error of the projection of a known smooth function onto $V_{\\mathrm{ms}}$ as $H\\to0$. If the decay is not exponential with ratio $\\beta<1$ or the error does not scale like $H^2/\\varepsilon^3$, then the propositions behind Theorem 4.3 fail for that sample.","supporting_citations":[{"cited_title":"Besse, B","cited_arxiv_id":null,"evidence_quote":"provides the Strang splitting stability and order estimates used in Theorems 4.1 and 4.2."},{"cited_title":"Wu and Z","cited_arxiv_id":null,"evidence_quote":"supplies the space decomposition $V_h=V_{\\mathrm{ms}}\\oplus W_h$ and the exponential localization of the multiscale basis that Lemma 4.3 and Theorem 4.3 inherit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"introduces the multiscale reduced-basis optimization construction for the Schrödinger equation with multiscale and random potentials and the eigendecomposition-based time stepping."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the time-splitting spectral discretization in the semiclassical regime and the $\\Delta t=O(\\epsilon)$, $h=O(\\epsilon)$ conditions behind SII."},{"cited_title":"Henning and A","cited_arxiv_id":null,"evidence_quote":"supplies the weighted quasi-interpolation estimate $\\|f-I_H f\\|\\le H^2\\|f\\|_{H^2}$ used in the spatial decomposition."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the KL-truncation and quasi-Monte Carlo error analysis for the NLSE with random potentials that Lemma 4.6 and Theorem 4.3 adapt."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the iterative decay argument behind the exponential localization of the multiscale basis functions."}],"review_version":1}