{"id":"b60427ce-ba8b-4974-9d05-13485063aa36","arxiv_id":"2606.23033","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Derives closed-form asymptotics for the anisotropic Rytova-Keldysh interaction and a closed-form expression for the 2D exciton ground-state binding energy using an anisotropic Gaussian variational ansatz.","lead":"This paper derives closed-form asymptotic expressions for the strongly anisotropic Rytova-Keldysh potential using steepest descent and momentum-space scaling. These enable an analytical expression for the ground-state binding energy of 2D excitons that accounts for competition between anisotropic mass and dielectric tensors.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Steepest-descent asymptotics for the 1D integral may lose uniform validity across anisotropy ratios and crossover distances","rationale":"The reader’s weakest_assumption directly names the two approximation steps that underwrite the headline analytical result. Because the full text was not previously available, the present check isolates exactly those steps; agreement follows because no other internal inconsistency appears more load-bearing once the derivations are examined.","tokens_in":1751,"tokens_out":356,"duration_ms":11866,"concrete_test":"Fix three representative anisotropy ratios (e.g., 1, 5, 20) and numerically quadrature the original 1D integral at 20–30 distances spanning short, intermediate and long regimes; compare pointwise to the paper’s closed-form asymptotics. If relative deviation exceeds ~5 % anywhere in the intermediate window used to justify the log-well approximation, recompute the variational binding energy with the numerical potential and check whether the analytic expression shifts by more than the claimed precision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on mapping the full spatial landscape via steepest descent plus momentum-space scaling of the exact 1D integral representation, then feeding the resulting short-range anisotropic log well into an anisotropic Gaussian variational solution of the Wannier equation. Steepest descent on a direction-dependent saddle can produce direction-specific error terms that are not controlled uniformly when the anisotropy ratio becomes large; the intermediate-range amplitude scaling (driven by the weakest polarizability) is especially sensitive to sub-leading contributions that the leading saddle may omit. If those corrections alter the effective short-range confinement felt by the variational wavefunction, the closed-form binding-energy expression no longer faithfully captures the mass-dielectric competition.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript maps the full spatial landscape of the strongly anisotropic Rytova-Keldysh potential by applying steepest-descent and momentum-space coordinate scaling to its exact one-dimensional integral representation, yielding closed-form asymptotics in all regimes. It then feeds the resulting short-range anisotropic logarithmic well into an anisotropic Gaussian variational ansatz for the Wannier equation, producing a closed-form analytical expression for the 2D exciton ground-state binding energy that encodes the competition between the effective-mass and dielectric tensors.","tokens_in":1883,"tokens_out":428,"duration_ms":23896,"significance":"If the asymptotics are uniformly controlled and the variational ansatz is adequate, the work supplies an analytical handle on anisotropy effects in 2D excitons that is currently missing from the literature; the closed-form binding-energy expression is a concrete strength that could be tested against numerics or experiment.","major_comments":[{"comment":"The steepest-descent analysis of the direction-dependent saddle (used to obtain the short-range logarithmic well that is then inserted into the variational ansatz) does not appear to supply uniform error estimates with respect to the anisotropy ratio. When the ratio becomes large, sub-leading corrections to the intermediate-range amplitude scaling (driven by the weakest polarizability) can alter the effective confinement felt by the Gaussian trial function, undermining the claim that the final binding-energy formula faithfully captures the mass-dielectric competition.","section":"asymptotic analysis of the 1D integral (short-range limit)"}],"minor_comments":[{"comment":"The abstract states that numerical studies already suggest robustness of an effective isotropic approximation; the manuscript should explicitly compare its new closed-form asymptotics against those numerical benchmarks in at least one figure or table.","section":null},{"comment":"Notation for the principal polarizabilities and the geometric-mean short-range coefficient should be introduced once and used consistently; occasional re-definition of symbols in different sections reduces readability.","section":null}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the positive summary and the detailed comment. We address the concern regarding the asymptotic analysis below.","responses":[{"response":"We thank the referee for highlighting the absence of explicit uniform error estimates with respect to the anisotropy ratio. The manuscript applies steepest descent to the exact 1D integral representation to extract the leading short-range anisotropic logarithmic asymptotics (governed by the geometric mean of the principal polarizabilities) and the intermediate-range amplitude scaling. While the derivation is performed for fixed anisotropy parameters and does not include bounds that are uniform as the ratio tends to infinity, the leading logarithmic term remains the dominant contribution in the short-distance regime relevant to the exciton ground state. The variational Gaussian ansatz is constructed precisely from this leading short-range form; the self-consistent optimization of the trial-function width is insensitive to sub-leading intermediate-range corrections at the level of the leading binding-energy expression. In the revised manuscript we will add a dedicated paragraph providing explicit error estimates for the steepest-descent approximation and demonstrating that, for anisotropy ratios of physical interest (up to order 10), the omitted corrections enter only at higher order in the binding energy and do not alter the leading mass-dielectric competition encoded in the closed-form formula. We therefore maintain that the central analytical result remains valid, while acknowledging that a more rigorous uniform-control analysis would strengthen the presentation.","revision_made":"yes","referee_comment":"[asymptotic analysis of the 1D integral (short-range limit)] The steepest-descent analysis of the direction-dependent saddle (used to obtain the short-range logarithmic well that is then inserted into the variational ansatz) does not appear to supply uniform error estimates with respect to the anisotropy ratio. When the ratio becomes large, sub-leading corrections to the intermediate-range amplitude scaling (driven by the weakest polarizability) can alter the effective confinement felt by the Gaussian trial function, undermining the claim that the final binding-energy formula faithfully captures the mass-dielectric competition."}],"tokens_in":1316,"tokens_out":424,"duration_ms":34515,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The core advance is the set of explicit asymptotic expressions for the potential in all regimes, plus the closed-form binding-energy formula that keeps the competition between the mass and dielectric tensors. The intermediate-range scaling tied to the weakest polarizability and the short-range log well set by the geometric mean of the polarizabilities are the concrete new pieces.\n\nThe work is straightforward: start from the known one-dimensional integral, apply standard steepest-descent and rescaling steps, extract the limiting forms, and plug the short-range piece into an anisotropic Gaussian trial function for the Wannier equation. That produces an analytical result without free parameters beyond the input tensors.\n\nThe main limitation is the absence of any shown numerical checks or error estimates on the asymptotics. The stress-test concern about uniform validity of the saddles when the anisotropy ratio grows large is reasonable on its face; if sub-leading terms matter in the crossover region, they could shift the effective confinement felt by the variational wavefunction. Without the explicit derivations or comparisons in the text, it is hard to judge how large that effect is.\n\nThis is aimed at theorists who calculate excitons in anisotropic 2D materials and want usable analytical expressions rather than repeated numerical integrals. A reader already working in that niche can test the formulas directly against their own numerics.\n\nI would send it to peer review so the derivations and any internal checks can be examined.","headline":"The paper turns the integral form of the anisotropic Rytova-Keldysh potential into closed-form asymptotics via steepest descent and scaling, then feeds the short-range part into a Gaussian variational solution for the exciton binding energy.","tokens_in":2365,"tokens_out":374,"would_cite":false,"duration_ms":18762,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"The anisotropic Rytova-Keldysh potential reduces to a closed-form anisotropic logarithmic well at short range that yields an explicit variational expression for 2D exciton ground-state binding energies.","keywords":["Rytova-Keldysh potential","anisotropic screening","2D excitons","Wannier equation","exciton binding energy","variational ansatz","dielectric anisotropy"],"falsifier":"Direct numerical solution of the Wannier equation with the exact integral potential for chosen anisotropy ratios, compared against the variational binding-energy formula, or spectroscopic measurement of exciton binding energies in a strongly anisotropic monolayer such as black phosphorus.","tokens_in":2628,"feed_emoji":"","tokens_out":758,"duration_ms":23393,"temperature":0.7,"pith_summary":"This paper maps the full spatial dependence of the Rytova-Keldysh interaction when in-plane anisotropy is strong by converting its one-dimensional integral representation into closed-form asymptotic expressions through steepest-descent and momentum-space coordinate scaling. The resulting short-range limit is an anisotropic logarithmic confinement whose strength is set by the geometric mean of the two principal polarizabilities. An anisotropic Gaussian variational trial function is then inserted into the Wannier equation, producing a compact analytical formula for the exciton binding energy that directly encodes the competition between the effective-mass tensor and the dielectric tensor. A sympathetic reader cares because the result supplies a practical, non-numerical route to exciton energies in directionally dependent atomically thin materials.","feed_headline":"Anisotropic Rytova-Keldysh potential yields closed-form exciton binding energy","feed_subtitle":"Short-range logarithmic well and Gaussian ansatz give explicit energy capturing mass-dielectric competition","key_machinery":"The one-dimensional integral representation of the anisotropic Rytova-Keldysh potential, reduced via steepest descent and coordinate scaling to closed-form asymptotics, together with the anisotropic Gaussian variational ansatz applied to the Wannier equation.","core_discovery":"By applying the method of steepest descent and momentum-space coordinate scaling to the exact one-dimensional integral representation, closed-form asymptotic expressions are obtained for the strongly anisotropic Rytova-Keldysh potential in every spatial regime; the short-range limit is an anisotropic logarithmic well governed by the geometric mean of the principal polarizabilities, which then permits an anisotropic Gaussian variational ansatz to furnish a closed-form analytical expression for the ground-state binding energy of 2D excitons that explicitly captures the competition between the effective mass and dielectric tensors.","pith_inferences":["The closed-form binding-energy expression could be inserted directly into models of optical response in van der Waals heterostructures without requiring repeated numerical integration.","The same steepest-descent reduction may apply to other nonlocal potentials written as one-dimensional integrals, offering a general route to analytic screening in anisotropic media.","Predicted binding energies can be tested against angle-resolved photoluminescence data on strained or naturally anisotropic 2D semiconductors."],"forward_implications":["Intermediate-range screening amplitude scales with the direction of weakest polarizability.","Short-range confinement is an anisotropic logarithmic well set by the geometric mean of the principal polarizabilities.","Exciton binding energy depends explicitly on both the effective-mass tensor and the dielectric tensor.","The effective isotropic approximation remains robust for highly directional systems because the derived asymptotics justify it."],"fun_headline_variants":["Anisotropic Rytova-Keldysh gives closed-form exciton binding","Short-range log well yields analytic 2D exciton energy","Gaussian ansatz solves anisotropic exciton binding equation","Steepest descent derives full anisotropic Rytova-Keldysh map"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The steepest-descent and coordinate-scaling procedures produce accurate closed-form asymptotics in every spatial regime and the anisotropic Gaussian variational ansatz is adequate to capture the ground-state solution of the Wannier equation.","fun_headline_variants_meta":{"raw":{"variants":["Anisotropic Rytova-Keldysh gives closed-form exciton binding","Short-range log well yields analytic 2D exciton energy","Gaussian ansatz solves anisotropic exciton binding equation","Steepest descent derives full anisotropic Rytova-Keldysh map"]},"model":"grok-4.3","cost_usd":0.004896,"raw_usage":{"total_tokens":2405,"prompt_tokens":678,"num_sources_used":0,"completion_tokens":61,"cost_in_usd_ticks":48962000,"prompt_tokens_details":{"text_tokens":678,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":1666,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":678,"tokens_out":61,"duration_ms":10851,"temperature":1.0,"reasoning_tokens":1666,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-26T07:24:42.386595+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Direct numerical solution of the Wannier equation with the exact integral potential for chosen anisotropy ratios, compared against the variational binding-energy formula, or spectroscopic measurement of exciton binding energies in a strongly anisotropic monolayer such as black phosphorus.","supporting_citations":[],"review_version":1}