{"id":"9c9d4d9b-5ada-467a-9bb6-53ededfa8a8e","arxiv_id":"1908.02395","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A variational width ansatz yields 1D and 2D dipolar Gross-Pitaevskii equations that reproduce 3D ground-state densities more accurately than the fixed-width reductions.","lead":"This paper derives improved one- and two-dimensional equations for elongated and flattened dipolar Bose-Einstein condensates, allowing the cloud width in the tightly trapped direction to vary instead of staying fixed. The reduced equations match full three-dimensional simulations better than the standard fixed-width reduction at moderate trap aspect ratios.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The slow-variation approximation (∇²φ≈∇⊥²φ) is the load-bearing assumption of both derivations, and its validity in the relaxed-trap regime claimed by the abstract is unquantified.","rationale":"The reader's weakest_assumption identifies the separable Gaussian ansatz with slowly varying width. I agree, and I focus specifically on the slow-variation part because it is explicitly invoked in both derivations (Section III after Eq. 10–11; Section IV after Eq. 28) and is not justified quantitatively. The paper's numerical tests show good visual agreement at the tested parameters, but they cover only two trap ratios per geometry, one interaction strength, and two polarization angles, and do not establish the 'relaxed away' claim for γ close to 1. If the omitted width-derivative terms become important for γ near 1, the variational equations themselves are not derived from the original 3D action, and the agreement seen at larger γ may not persist. I do not find an internal algebraic inconsistency in the equations: in the non-dipolar limit the width equations reduce to the known NPSE forms, and the DDI terms appear to retain the correct structure. One minor reporting issue: in Section I, 'in the manner of Salasnich et al. [64]' should probably cite [68] (the NPSE paper), while [64] is the fixed-width dipolar reduction; this is a citation error and does not affect the argument. Because the reader already assigned CONDITIONAL, my read does not change the verdict; the concern strengthens the need for the proposed quantitative check of the slow-variation approximation.","tokens_in":14292,"tokens_out":13511,"duration_ms":143242,"concrete_test":"Using the converged 3D ground states for cigar traps at γ=20, 10, 5, and 2 (with ϵ_dd=0.9 and β̃1D=2, as in Figs. 3–4), fit a Gaussian transverse profile at each z to extract σ(z), then compute the dimensionless ratio R = ∫|f(z)|² (∂σ/∂z)² dz / ∫|f(z)|² dz. If R≪0.05 for all γ, the slow-variation approximation is supported; if R grows to ≳0.1 as γ→1, the central claim of accuracy in relaxed traps is not justified and the reduced equations should be compared against 3D at these γ quantitatively (e.g., via L2 density error) before trusting the abstract's scope.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The 1D and 2D reductions both depend on the explicit slow-variation assumption stated in Section III ('we assume that φ is slowly varying along the z-direction, allowing us to write ∇²φ≈∇⊥²φ') and Section IV (analogously for η). This drops all axial (or in-plane) derivative terms of the width, including (∂σ/∂z)² and ∂²σ/∂z², from the effective action. The paper's central claim is accuracy 'even as the trapping is relaxed away from the strict quasi-one- and quasi-two-dimensional regimes,' yet no quantitative estimate of the omitted terms is given, and the validation is visual agreement at γ=10 and 80 (1D) and γ=1/10 and 1/80 (2D) for a single interaction strength (ϵ_dd=0.9) and two polarization angles. As γ approaches 1, the axial density profile becomes more strongly modulated and σ(z) can vary on short length scales, so the neglected terms are not obviously small. If they are non-negligible, the reduced equations would lose accuracy exactly in the regime the abstract claims to cover. The paper should demonstrate that the omitted kinetic-energy terms are small over the claimed parameter range.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper derives effective one-dimensional and two-dimensional Gross-Pitaevskii equations for cigar- and pancake-shaped dipolar Bose-Einstein condensates, extending the non-polynomial Schrödinger equation approach to dipolar systems. The authors start from the 3D dipolar GPE action, assume a separable wavefunction whose transverse (1D case) or axial (2D case) profile is a Gaussian with variational width, integrate out the tightly confined direction(s), and obtain a coupled PDE for the low-dimensional wavefunction and algebraic equations for the width(s). They treat the polarization angle generically, but validate numerically only for α=0 and α=π/2, comparing ground-state densities and aspect ratios with full 3D dGPE simulations at trap ratios γ=10,80 (cigar) and γ=1/10,1/80 (pancake), with interaction strength ϵ_dd=0.9 for densities. The central claim is that the variable-width equations agree with the 3D results significantly better than the fixed-width reductions, including in moderately relaxed geometries.","tokens_in":14558,"tokens_out":19508,"duration_ms":202992,"significance":"If the equations are correct and the approximations are under control, this is a genuinely useful contribution: it provides computationally cheaper reduced models for dipolar BECs that capture magnetostriction and variable confinement without fitting any parameters to the 3D benchmarks. The derivation is systematic, based on an action functional, and the comparison is made against independent 3D dGPE solutions. The approach could benefit studies of solitons, rotons, and dynamical phenomena in elongated and flattened dipolar condensates. However, the claimed regime of validity is broader than the presented evidence: the slow-variation approximation is not quantified, the validation is mostly visual and covers a narrow parameter set, and one of the central equations appears to contain a typo. These issues prevent acceptance in the current form.","major_comments":[{"comment":"The width equation as printed appears to have a load-bearing typo. For ϵ_dd=0 it reduces to 1−γ²σ̃²+β̃1D|f̃|²=0, i.e. σ̃²=(1+β̃1D|f̃|²)/γ²; with the Fig. 3 parameters (γ=10, β̃1D=2) and a typical |f̃|²≈0.2 this gives σ̃≈0.12, whereas the plotted σ(z) is approximately 0.3–0.38. The standard NPSE-type result is σ̃⁴=(1+...)/γ², and indeed the anisotropic width equations (21)–(22) contain σ̃_x⁴ and σ̃_y⁴. The α=0 limit of Eqs. (21)–(22) therefore cannot reproduce Eq. (18) as written. Please correct Eq. (18) and verify that the numerical code actually solves the corrected equation.","section":"§III, Eq. (18)"},{"comment":"The reduction drops all derivatives of the variational width with respect to the free coordinate, including (∂_zσ)² and ∂_z²σ in the 1D case and the analogous in-plane terms for η in the 2D case. The abstract claims accuracy 'even as the trapping is relaxed away from the strict quasi-one- and quasi-two-dimensional regimes,' but no estimate of the omitted terms is given. At γ=10 the plotted σ(z) varies by roughly 10–20% across the cloud, so the slow-variation assumption is not obviously negligible. Please provide a quantitative estimate of the neglected kinetic-energy contributions over the claimed parameter range (e.g., the ratio of the omitted terms to the retained transverse kinetic energy), and ideally test additional trap ratios closer to γ=1.","section":"§III after Eq. (11); §IV after Eq. (28)"},{"comment":"The validation is qualitative and narrowly sampled. Density comparisons are shown only for ϵ_dd=0.9, α=0 and α=π/2, and two trap ratios per geometry; the conclusion states agreement 'over a range of trap ratios, polarization angles and interaction strengths,' but the interaction-strength dependence is demonstrated only for the aspect ratio, not for density profiles. The abstract's 'arbitrary polarization angle' is not tested (e.g., α=π/4 is absent). Please include quantitative error measures (relative L2 or H1 density differences, chemical potentials, or peak-density errors) and additional parameter points to substantiate the claims of 'strong agreement' and 'significant improvement' over the fixed-width models.","section":"§III and §IV, Figs. 3–8"},{"comment":"The reported grid parameters appear inconsistent with the physical sizes of the condensates shown in the figures. For pancake geometries, 256 points with ∆z=0.0008𝓁⊥ give a total z-extent of 0.205𝓁⊥, which is smaller than the axial width η≈0.31–0.35𝓁⊥ displayed in Fig. 6. For cigar geometries with γ=80, ∆z=0.04𝓁⊥ gives a z-extent of about 1.15𝓁z, which may be too small to contain the axial density shown in Fig. 3 if the axis is in units of 𝓁z. Please clarify the actual box sizes or correct the grid-spacing values; if the boxes are indeed that small, the 3D dGPE results used as benchmarks could be affected by truncation.","section":"§II, Numerical methods"}],"minor_comments":[{"comment":"The word 'intrego-differential' should be 'integro-differential'.","section":"§III B"},{"comment":"The word 'discrepany' should be 'discrepancy'.","section":"Fig. 6 caption"},{"comment":"There appears to be a missing '+' sign between '-ℏ²/(4mη²)' and '1/(4mωz²η²)'; as printed, the second term is dimensionally inconsistent and should likely read '+(1/4)mωz²η²'.","section":"Eq. (29)"},{"comment":"The definition of Φ2D_dd contains a trailing 'f' that seems spurious: the Fourier bracket should act on |f|² only, as in Eq. (32).","section":"Eq. (30)"},{"comment":"References [22] and [72] are the same work (H.-Y. Lu et al., Phys. Rev. A 82, 023622 (2010)); please merge them.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the general approach is sound, but the apparent typo in Eq. (18) and the unquantified slow-variation approximation are substantive issues that must be addressed. The numerical grid parameters also need clarification. I see no reason to doubt the authors' good faith; the figures presumably come from a corrected code, but the printed equations and the validation are not yet reliable enough for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nHere's my read on Knight et al. The core result is new and useful: they extend the Salasnich variable-width NPSE to dipolar BECs, with anisotropic widths and arbitrary polarization angle, for both cigar and pancake geometries. The derivation is clean—action functional, Euler-Lagrange equations, no fitted constants—and the comparison against the 3D dGPE shows real improvement over the fixed-width reduction at the tested parameters. The aspect-ratio versus epsilon_dd curves (Figs. 5 and 8) are a nice quantitative check, especially because the fixed-width model would be flat at unity.\n\nThe main soft spot is the narrow validation. The density comparisons are all at epsilon_dd = 0.9, two polarization angles, and two trap ratios per geometry. That's a thin slice of the parameter space. Agreement is shown only visually; there are no error metrics like L2 density differences or chemical potential differences. For a paper whose selling point is accuracy in the relaxed regime, that's a gap. The numerical setup is also under-specified: no code or data release, and no convergence checks on the 3D solver's grid and cutoff. Those are fixable but should be addressed.\n\nThe stress-test concern about the slow-variation approximation is real. In Section III they drop axial derivatives of the transverse profile, including terms like (d sigma/dz)^2 and d^2 sigma/dz^2. As gamma approaches 1, sigma(z) can modulate more, and those terms may not be negligible. They don't quantify the omitted terms, and they don't test near gamma = 1. I wouldn't call it a fatal flaw—the benchmarks at gamma = 10 and 1/10 are reasonably far from quasi-1D/2D and still agree well visually—but the abstract's 'even as the trapping is relaxed' claim exceeds the evidence. A proper revision should either scan gamma down toward 1 or estimate the size of the dropped kinetic terms.\n\nWho is this for? Anyone doing theory of dipolar BECs in elongated or flat traps who needs a cheap 1D/2D model that captures width variation. It deserves a serious referee; the derivation alone is worth checking, and the benchmark is a useful service. I'd send it to peer review with a request for quantitative errors, broader parameter scans, and code/data release.","headline":"Useful extension of variable-width NPSE to dipolar BECs, but validation is thin and the relaxed-regime claim outruns the evidence.","tokens_in":15056,"tokens_out":2165,"would_cite":true,"duration_ms":23311,"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":"This paper claims that keeping the width of a dipolar condensate's tightly-confined Gaussian profile as a variational field, instead of fixing it to the oscillator length, yields effective one- and two-dimensional Gross-Pitaevskii…","keywords":["dipolar Bose-Einstein condensate","Gross-Pitaevskii equation","dimension reduction","variational width","non-polynomial Schrödinger equation","magnetostriction","anisotropic trap","quasi-one-dimensional and quasi-two-dimensional regimes"],"falsifier":"Compute the 3D dipolar GPE ground state for a moderately elongated cigar ($\\gamma\\approx10$) with stronger dipoles ($\\epsilon_{dd}\\to1$) and larger $s$-wave nonlinearity, and compare the actual transverse density slices at the cloud's centre with the best-fit Gaussian of width $\\sigma(0)$; if the slices are visibly non-Gaussian (flat-topped, bimodal, or asymmetric in a way the two-width ansatz cannot represent) while the variable-width equations still force a Gaussian, the claim that the method tracks the 3D ground state in the relaxed regime is refuted.","tokens_in":14110,"feed_emoji":"🧲","tokens_out":12297,"duration_ms":115692,"temperature":0.7,"pith_summary":"This paper derives effective one- and two-dimensional equations for cigar-shaped and pancake-shaped dipolar Bose-Einstein condensates, ultracold atomic clouds whose magnetic dipoles are aligned by an external field. The key move is to let the width of the condensate in the tightly confined direction(s) vary with position and time, and to fix it by minimizing the energy, rather than locking it to the harmonic-oscillator length as standard low-dimensional reductions do. The authors show that the resulting coupled equations reproduce the ground-state densities of the full three-dimensional dipolar Gross-Pitaevskii equation much more accurately than the fixed-width reduction, and that the improvement persists when the trap is only moderately elongated or flattened. This matters because the strict quasi-one- and quasi-two-dimensional regimes are hard to reach experimentally, and cheaper but accurate low-dimensional models make studies of rotons, solitons, and supersolids far more tractable.","feed_headline":"Letting cloud width vary sharpens low-dim dipolar gas equations","feed_subtitle":"New 1D/2D equations match full 3D ground states even when traps relax away from strict quasi-low-dimensional limits.","key_machinery":"The load-bearing object is the separable Gaussian ansatz with a variational width. For a cigar-shaped condensate the wavefunction is written as a 1D field along $z$ times a Gaussian in the transverse plane whose width $\\sigma(z,t)$ is itself a field; for a pancake it is a 2D field in the plane times a Gaussian in $z$ with width $\\eta(x,y,t)$. Substituting this ansatz into the 3D action functional and integrating out the tightly confined coordinate(s), while approximating the Laplacian of the transverse profile by its transverse (or axial) part alone, produces a reduced action. The Euler-Lagrange equations then give a coupled system: a Gross-Pitaevskii-type equation for the low-dimensional wavefunction and an algebraic equation for the width(s), where the width equation inherits extra interaction terms from the derivative of the nonlocal dipolar potential. Setting the width to the fixed harmonic-oscillator length recovers the conventional quasi-1D/quasi-2D equations, so the variational width is exactly what carries the improvement.","core_discovery":"The central claim is that a variational-width Gaussian ansatz, $\\psi = f(z,t)\\phi(x,y;\\sigma(z,t))$ for cigars and $\\psi = f(x,y,t)\\phi(z;\\eta(x,y,t))$ for pancakes, yields effective 1D and 2D dipolar Gross-Pitaevskii equations that stay close to the full 3D solutions far outside the strict quasi-1D/quasi-2D limits. The reduced system is a coupled set: a partial differential equation for the low-dimensional wavefunction together with algebraic equations for the width(s), obtained by varying the action with respect to both. For cigar-shaped condensates with dipoles tilted away from the long axis, the transverse profile is allowed to be anisotropic with two widths $\\sigma_x$ and $\\sigma_y$, which lets the equations reproduce the dipole-induced magnetostriction; for pancakes the axial width $\\eta(x,y,t)$ varies across the plane. In ground-state comparisons against full 3D numerical solutions, the variable-width equations match the density profiles and the transverse aspect ratio as a function of dipole strength for trap ratios $\\gamma=80$ and $\\gamma=10$ (cigars) and $\\gamma=1/10$ and $\\gamma=1/80$ (pancakes), with $\\epsilon_{dd}=0.9$ and both $\\alpha=0$ and $\\alpha=\\pi/2$ polarizations, and they substantially outperform the fixed-width reduction in the moderately trapped cases.","pith_inferences":["The same reduction, applied to time-dependent problems, could serve as a fast diagnostic for dipole-driven instabilities: the algebraic width equations should develop non-convex or multiple solutions at the onset of collapse or roton softening, signalling instability without a full 3D simulation.","The two-width cigar equations are a natural candidate for resolving the known disagreement between 3D and standard 1D predictions for dark-soliton oscillations in dipolar gases; the authors flag this as future work, not a demonstrated result.","Adding beyond-mean-field (Lee-Huang-Yang) corrections to the action before carrying out the variational reduction should yield low-dimensional supersolid models; the paper notes the need for quantum fluctuations but does not implement them.","The transverse aspect-ratio curves the paper computes are directly measurable, so a focused experiment on in-situ aspect ratio versus dipole strength could test the reduced equations against 3D physics without resolving the full 3D profile."],"forward_implications":["For cigar-shaped dipolar BECs with dipoles aligned along the long axis, the variable-width 1D equations match 3D ground-state densities closely at trap ratios $\\gamma=80$ and $\\gamma=10$, where the fixed-width equations visibly deviate.","For dipoles polarized perpendicular to the long axis, allowing two anisotropic widths $\\sigma_x$ and $\\sigma_y$ reproduces the 3D transverse aspect ratio as a function of $\\epsilon_{dd}$, something a fixed isotropic width cannot do.","For pancake-shaped BECs, the variable-width 2D equations match 3D densities at $\\gamma=1/10$ and $1/80$, correcting the fixed-width underestimate at the moderate ratio.","Because the reduced equations are far cheaper to solve than the full 3D dGPE, they offer a practical route to studying rotons, solitons, and collapse in elongated or flattened dipolar gases while retaining key 3D effects."],"supporting_citations":[{"why":"Supplies the variational-width Gaussian ansatz and Euler-Lagrange reduction procedure that this paper extends from non-dipolar to dipolar condensates.","marker":"[68]"},{"why":"Provides the standard fixed-width 1D/2D dipolar GPEs that the paper compares against and improves upon.","marker":"[64]"},{"why":"Gives the 3D dipolar Gross-Pitaevskii model and the Fourier-space dipole-dipole interaction kernel used in the derivation.","marker":"[10]"},{"why":"Shows the predecessor non-polynomial Schrödinger equation captured 3D collapse physics, motivating the same strategy for dipolar gases.","marker":"[66]"},{"why":"Supplies the split-step Fourier imaginary-time method used to obtain the 3D reference ground states.","marker":"[70, 71]"},{"why":"Provides the semi-analytic Fourier transform of the truncated dipole potential used in the numerical comparisons.","marker":"[72]"},{"why":"Establishes (with [64]) that fixed-width reductions depart from 3D results outside the strict quasi-low-dimensional regime, the target of this paper.","marker":"[69]"}],"fun_headline_variants":["Variable-width ansatz sharpens 1D/2D dipolar BEC equations","Adaptive width makes low-dim dipolar BEC equations match 3D","Cigar and pancake BEC equations improve with variational width","Variable-width scheme handles tilted dipoles in low-dim BEC equations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results rely on the condensate's profile in the tightly confined direction(s) remaining close to a Gaussian that changes only slowly along the free direction; if the true profile is strongly non-Gaussian or the width varies too rapidly, the reduced equations will not reproduce the 3D behaviour.","fun_headline_variants_meta":{"raw":{"variants":["Variable-width ansatz sharpens 1D/2D dipolar BEC equations","Adaptive width makes low-dim dipolar BEC equations match 3D","Cigar and pancake BEC equations improve with variational width","Variable-width scheme handles tilted dipoles in low-dim BEC equations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000734,"raw_usage":{"total_tokens":3306,"prompt_tokens":990,"completion_tokens":2316,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":606,"completion_tokens_details":{"reasoning_tokens":2236}},"tokens_in":606,"tokens_out":2316,"duration_ms":17147,"temperature":1.0,"reasoning_tokens":2236,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:45:32.638353+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the 3D dipolar GPE ground state for a moderately elongated cigar ($\\gamma\\approx10$) with stronger dipoles ($\\epsilon_{dd}\\to1$) and larger $s$-wave nonlinearity, and compare the actual transverse density slices at the cloud's centre with the best-fit Gaussian of width $\\sigma(0)$; if the slices are visibly non-Gaussian (flat-topped, bimodal, or asymmetric in a way the two-width ansatz cannot represent) while the variable-width equations still force a Gaussian, the claim that the method tracks the 3D ground state in the relaxed regime is refuted.","supporting_citations":[{"cited_title":"Salasnich, A","cited_arxiv_id":null,"evidence_quote":"Shows the predecessor non-polynomial Schrödinger equation captured 3D collapse physics, motivating the same strategy for dipolar gases."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the semi-analytic Fourier transform of the truncated dipole potential used in the numerical comparisons."},{"cited_title":"Baillie and P","cited_arxiv_id":null,"evidence_quote":"Establishes (with [64]) that fixed-width reductions depart from 3D results outside the strict quasi-low-dimensional regime, the target of this paper."}],"review_version":1}