{"id":"e096f3f9-2643-48f8-b97c-7d7c349d94e2","arxiv_id":"2501.09641","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For quasi-2D dipolar Bose condensates, an in-plane polarization component turns the perpendicular-polarization critical point into three critical lines and makes the hexagonal and honeycomb supersolids axially deformed with anisotropic superfluid response.","lead":"Tilting the magnetic polarization of a flat dipolar Bose gas splits the single critical point of the known phase diagram into three separate critical lines, producing compressed hexagonal, stripe, and stretched honeycomb supersolid phases. The work maps this new phase diagram and shows the tilt makes superfluidity anisotropic.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lattice-symmetry restriction in Eqs. (8)-(9) is the load-bearing assumption; a competing oblique or square phase could change the reported critical-line topology.","rationale":"The paper's strongest claim is that tilting the polarization replaces the alpha=0 critical point with three critical lines separating homogeneous-stripes-compressed-hexagonal and stripes-stretched-honeycomb, with first- and second-order parts and triple points. The argument's load-bearing step is the energy comparison among the allowed states. The reader identified the restricted ansatz as the weakest assumption; I agree. The alpha=0 benchmark is a good check, and the single- vs many-mode comparison at alpha=30 degrees shows quantitative shifts but no topological change for the allowed geometries. However, neither check explores other Bravais lattices. In a system with an in-plane polarization, the interaction is anisotropic and could favor a lattice with shear that is not representable by Eq. (9); the one-parameter theta controls only uniaxial compression/stretching and fixes the y-components of the reciprocal-lattice vectors. The superfluid anisotropy results and the first/second-order assignments inherit this limitation. There is no reason to suspect the authors of oversight beyond a standard variational restriction; the paper is honest about the ansatz choices. A single unrestricted minimization at a few phase-diagram points would settle the concern. Because the reader's CONDITIONAL verdict already reflects this uncertainty, no verdict change is needed.","tokens_in":12893,"tokens_out":5261,"duration_ms":59870,"concrete_test":"Run a fully unconstrained 2D minimization of the functional Eq. (5) for alpha=30 degrees at several representative points, e.g., (rho, as/add) = (120, 0.755), (350, 0.779), and points along each critical line, using imaginary-time evolution or a plane-wave basis with no imposed lattice symmetry and starting from both random noise and the predicted compressed-hexagonal/stretched-honeycomb states. If a converged pattern with lower energy has an oblique, square, or otherwise excluded symmetry, the phase diagram and possibly the three-critical-line topology change; if the energy minimum always lies in the assumed family, the ansatz restriction is validated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central phase diagram is built by comparing energies of only three variational families: homogeneous, stripes with wavevector along y (Eq. 8), and a one-parameter set of centered-rectangular lattices obtained by uniaxial deformation of the hexagonal/honeycomb structure (Eq. 9). This excludes oblique (sheared) lattices, square lattices, and any pattern whose reflection symmetries differ from those assumed in Sec. II C. Because tilting the polarization breaks continuous rotational symmetry, the true ground-state lattice need not belong to this restricted family; if a sheared triangular or other oblique structure has lower energy in a region of the (rho, as/add) plane, the reported critical lines, their first- and second-order segments, and the two triple points would shift or even change topology. The many-mode calculation in Sec. IV validates convergence within the same lattice-geometry class only; it does not test alternative symmetries. The transverse-width approximation (sigma taken from the homogeneous solution for all phases) is a secondary quantitative concern. Thus the main claim, while plausible and consistent with the alpha=0 benchmark, is conditional on an untested lattice-symmetry assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a quasi-2D dipolar Bose gas harmonically confined along z, with the polarization direction tilted by an angle α from the normal. The authors minimize an extended Gross-Pitaevskii energy functional including LHY quantum fluctuations, using a Thomas-Fermi ansatz in the transverse direction and three planar variational families: homogeneous, stripes with wavevector along y, and centered-rectangular lattices obtained by uniaxial compression/stretching of hexagonal/honeycomb states. Their central claim is that for α>0 the single critical point of the α=0 phase diagram is replaced by three critical lines, two triple points, and phase boundaries containing both first- and second-order segments. They further report that tilting deforms the hexagonal and honeycomb patterns along the in-plane polarization direction and makes the superfluid fraction tensor anisotropic. A many-mode spectral calculation at α=30° is presented as a check of the single-mode phase diagram.","tokens_in":13136,"tokens_out":5513,"duration_ms":63116,"significance":"If the central topology claim holds, the paper significantly extends the known phase diagram of planar dipolar supersolids and makes concrete, experimentally testable predictions for the lattice spacing and superfluid anisotropy of 162Dy samples. The work has genuine strengths: the α=0 limit reproduces the previously established hexagonal/stripes/honeycomb diagram; the energy functional is minimized from first principles with no experimental input; the superfluid fraction is computed from the ground-state density via a well-defined auxiliary problem; and the many-mode check at α=30° is a good-faith convergence test. The main reservation is that the phase diagram is obtained by comparing only three symmetry-restricted variational families, so the claimed topology is conditional on an untested lattice-symmetry assumption.","major_comments":[{"comment":"The ground-state search is restricted to homogeneous states, stripes with wavevector along y, and centered-rectangular lattices with basis vectors e1=(0,1) and e2=(1/2 cot θ, -1/2). This family has no shear (oblique) degree of freedom and forces reflection symmetry about the x and y axes. Because tilting the polarization breaks continuous rotational symmetry, competing oblique, rotated-stripe, square, or other lattice symmetries are not excluded a priori. All phase boundaries in Fig. 2 are obtained by comparing energies within this restricted set, so the central claim of three critical lines and two triple points is conditional. I ask for a concrete test, e.g., a variational calculation with a general oblique lattice (allowing e2=(a,b)) and with a rotated stripe orientation at α=30°, or an imaginary-time real-space evolution; if an excluded symmetry wins anywhere in the (ρ, as/add) plane, the reported topology must be revised.","section":"Sec. II.C, Eqs. (8)-(9)"},{"comment":"The phase diagrams for all α values in Fig. 2 are computed in the single-mode approximation, retaining only the first Fourier harmonic of each ansatz. The many-mode comparison in Fig. 3 is performed only at α=30°. While this check supports the topology at that angle, it does not validate the evolution of the critical lines and first-order boundaries as α is varied, in particular the convergence back to a single critical point as α→0. I request many-mode calculations for at least one additional tilt angle, or a quantitative error estimate for the single-mode boundaries over the full α range shown.","section":"Sec. III, Fig. 2 and Sec. IV, Fig. 3"},{"comment":"The transverse Thomas-Fermi width σ is minimized only for the homogeneous solution and then reused for the stripes, hexagonal, and honeycomb ansätze. Since σ enters the LHY term as σ^{-3/2} and multiplies all interaction terms, this approximation can bias the relative energies of the modulated phases and shift the reported boundaries. The statement that σ 'does not differ significantly' is not quantified. I ask for a test in which σ is an independent variational parameter for each modulated phase, or at least a numerical comparison of the optimized vs. homogeneous σ at representative points in the phase diagram.","section":"Sec. II.A, text after Eq. (7)"}],"minor_comments":[{"comment":"The phrase 'all transition lines contain first- and second-order regions' is not defined precisely; please clarify whether each phase boundary changes order at an isolated critical point or whether the statement refers only to the three critical lines.","section":"Abstract and Sec. III"},{"comment":"The claim that the single-mode approximation becomes exact for the stripes-homogeneous boundary is asserted rather than derived; a short argument showing that higher harmonics vanish at that boundary would strengthen the statement.","section":"Sec. IV"},{"comment":"The derivation of the effective potential F(qx,qy,α) is compressed; in particular, the scaling of qx and qy by σ/cos(α) and its relation to the rotated coordinate system would benefit from an explicit intermediate step.","section":"Eq. (6)"},{"comment":"The pink stripe marking the transition region is said to have width 1°, but the criterion used to locate the transition from the computed data is not specified; please state how the transition point and the width were determined.","section":"Sec. IV, Fig. 4"},{"comment":"References [60] and [64] are the same paper by Blakie; please merge them to avoid duplicate citations.","section":"References"},{"comment":"There are minor typographical issues, including 'anzats' for 'ansatz' in Sec. II.C and inconsistent use of 'C.P.'/'T.P.' acronyms in the figure captions; these should be cleaned up.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for a quantum-gas journal and the α=0 benchmarking is solid. My main concern is the restricted ansatz family, which is a testable assumption rather than a fatal flaw; I would like the revision to include an explicit search over oblique/rotated lattice symmetries at least at α=30°, since the central topology claim depends on it. I also note that a duplicate reference and a few wording issues should be corrected during revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Daniel, Grossklags, and collaborators present a mean-field variational study of a quasi-2D dipolar BEC with polarization tilted at angle α relative to the plane normal. The genuinely new result is the α>0 phase diagram: the single critical point known for perpendicular polarization splits into three critical lines with two triple points, and the hexagonal/honeycomb lattices deform along the in-plane polarization axis, producing anisotropic superfluidity. The α=0 limit correctly reproduces the established phase diagram including the stripes phase from Ripley et al., which gives me confidence the machinery is sound.\n\nThe paper does several things well. The derivation in Sec. II is standard and internally consistent. The many-mode computation at α=30° supports the single-mode critical-line topology, and the superfluid-tensor method is carefully checked by recovering isotropy at α=0. The experimental estimates with 162Dy parameters are concrete and useful.\n\nThe main soft spot is exactly what the stress-test note identifies. The ground-state search is restricted to homogeneous, stripes with wavevector perpendicular to the in-plane polarization, and a one-parameter family of axially deformed hexagonal/honeycomb lattices. No oblique, sheared, square, or rotated-stripe patterns are tested. Since the tilt breaks continuous rotational symmetry, the true ground state need not belong to this family. If any competing lattice wins anywhere in the (ρ, as/add) plane, the reported critical lines and triple points could shift or change topology. The many-mode calculation does not address this: it only tests convergence within the same lattice-geometry class. So the central claim is conditional on an untested symmetry assumption. That is a genuine limitation, and the paper would be stronger if it acknowledged it explicitly or tested at least a few alternative structures.\n\nSecondary concerns are milder: the transverse width σ is taken from the homogeneous solution for all modulated phases, which could shift boundaries quantitatively, and no code, data, or convergence details are provided, limiting independent verification. These are standard issues for variational studies and not disqualifying.\n\nOverall, this is a competent and honest mean-field study of an unexplored regime, with a plausible central claim and a good benchmark. It deserves a serious referee. The referee should ask for either a wider variational search at representative points or a clear statement of the restricted-symmetry caveat. I would read the follow-up before citing this as established topology, but the predictions are well-motivated and worth having on record.","headline":"A solid, well-benchmarked mean-field study whose main topology claim rests on a restricted lattice-symmetry ansatz that the paper does not test; worth refereeing with that caveat pressed.","tokens_in":13642,"tokens_out":3597,"would_cite":true,"duration_ms":35332,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["67.85.-d","03.75.Hh"],"model":"deepseek-v4-flash","headline":"A tilted polarization can replace one critical point with three critical lines in a dipolar supersolid phase diagram.","keywords":["dipolar Bose-Einstein condensates","supersolids","tilted polarization","phase diagram","critical lines","stripes phase","hexagonal and honeycomb phases","anisotropic superfluidity"],"falsifier":"Perform an unconstrained ground-state search on the same mean-field functional, allowing arbitrary lattice symmetries, oblique wavevectors, and stripe orientations not aligned with the in-plane polarization. If any pattern outside the three ansätze has lower energy in the regions the paper assigns to compressed hexagonal, stripes, or stretched honeycomb—especially near the two triple points—then the three-critical-lines topology is not the complete ground-state diagram.","tokens_in":12675,"feed_emoji":"🧲","tokens_out":7113,"duration_ms":71372,"temperature":0.7,"pith_summary":"This paper studies a planar (quasi-2D) dipolar Bose-Einstein condensate and asks what happens when the polarization direction is tilted relative to the plane's normal. For strictly perpendicular polarization, the known phase diagram shows hexagonal, stripes, and honeycomb supersolid phases meeting at a single critical point. The paper claims that any in-plane polarization component destroys this point and replaces it with three critical lines, each separating two phases, with first- and second-order segments and two triple points. The tilt also compresses the hexagonal lattice and stretches the honeycomb lattice along the in-plane polarization axis, making the superfluid fraction direction-dependent: superfluidity is enhanced along the tilt direction and suppressed perpendicular to it. If correct, this gives a practical handle for tuning both the order of supersolid transitions and the anisotropy of supersolid transport in ultracold dipolar gases.","feed_headline":"Tilting polarization splits a critical point into three lines","feed_subtitle":"An in-plane polarization deforms the supersolid lattices and makes superflow direction-dependent.","key_machinery":"The central object is an effective 2D mean-field energy functional for the in-plane condensate wave function, obtained by integrating out a Thomas–Fermi transverse profile and expressing the dipole–dipole interaction through an effective potential with Fourier kernel F(qx, qy, α) = ∫ dq/2π [3(−q cos q + sin q)/$q^{3}$]^2 ($q_y^{2}$ + (q_x − q tan α)^2)/($q_y^{2}$ + $q^{2}$ + (q_x − q tan α)^2), together with the Lee–Huang–Yang quantum-fluctuation term. The variational search compares three ansätze: the homogeneous state, stripes oriented along the in-plane polarization (the x-axis), and hexagonal lattices that can be compressed or stretched along the same axis via the angle θ. Minimizing over Fourier amplitudes, the modulation wavevector k0, and θ locates the phases. The superfluid fraction tensor is computed by solving ∇·(ρ(r)∇Ki) = ∂ρ(r)/∂xi through minimization of a convex action, which needs only the Fourier harmonics present in ∂ρ(r)/∂xi.","core_discovery":"The authors establish that the ground-state phase diagram of a planar dipolar Bose gas with polarization tilted by angle α from the plane normal is organized by three critical lines instead of the single critical point present at α = 0. In the density–scattering-length plane (ρ, as/add), the homogeneous superfluid meets the stripes phase along one continuous critical line, while a compressed-hexagonal phase and a stretched-honeycomb phase each meet the stripes phase along their own critical lines; first-order transitions continue these boundaries away from criticality, and triple points mark where homogeneous, stripes, and each lattice phase coexist. The same effective anisotropy along the in-plane polarization deforms the density patterns, giving the lattice spacing ratio Dx/Dd = 2 sin(θ), with the hexagonal phase more strongly deformed than the honeycomb phase. The authors compute the resulting superfluid fraction tensor and find that f_xx increases and f_yy decreases with α in the lattice phases, while the stripes phase keeps f_xx = 1. A converged many-mode calculation for α = 30° confirms the single-mode topology, with the stripes–homogeneous critical line unaffected by higher-order harmonics.","pith_inferences":["The splitting of the critical point can be understood as the in-plane tilt breaking the rotational degeneracy of the soft mode: a free-energy expansion around the multicritical point should reproduce the fan of three critical lines and the two triple points.","Because the contrast of the density modulation changes little while f_xx and f_yy diverge with α, superfluid anisotropy is not simply a proxy for modulation depth; measuring dipole-oscillation frequencies or moment of inertia along x and y could test the predicted directional response.","The paper's restriction to three ansätze leaves room for oblique or square modulated states; an unconstrained ground-state search near the triple points could either confirm the topology or reveal additional phases that the present variational family misses.","The mechanism is generic for long-range interacting planar systems: any in-plane easy axis in a dipolar supersolid should generically induce anisotropic superfluidity and split degenerate multicritical points."],"forward_implications":["For any nonzero in-plane polarization component, the single α = 0 critical point splits into three critical lines and two triple points, changing the topology of the supersolid phase diagram.","The homogeneous-to-stripes transition becomes second order over a wide intermediate-density region, resembling phase diagrams of quasi-one-dimensional dipolar systems.","The hexagonal and honeycomb phases are axially deformed along the in-plane polarization direction, with the spacing ratio set by Dx/Dd = 2 sin(θ), producing anisotropic long-distance properties.","In the deformed lattice phases, the superfluid fraction increases along the in-plane polarization direction and decreases perpendicular to it; stripes remain fully superfluid along the stripe direction.","The parameter regime is compatible with current experiments using 162Dy, with lattice spacings around 4.6 µm, making the predicted anisotropic superfluidity experimentally accessible."],"supporting_citations":[{"why":"Establishes the perpendicular-polarization phase diagram including the stripes phase and the single critical point that this work generalizes to tilted polarization.","marker":"[49]"},{"why":"Reports the original mean-field critical point and hexagonal/honeycomb/stripes phase diagrams for the untilted case, providing the baseline that the tilt modifies.","marker":"[48]"},{"why":"Predicts the modulated phases in planar dipolar BECs whose competition is central to the phase diagram studied here.","marker":"[43]"},{"why":"Supplies the Lee–Huang–Yang quantum-fluctuation correction term used in the energy functional to stabilize the modulated phases.","marker":"[33]"},{"why":"Provides the moving-walls formulation of the superfluid fraction tensor that the paper adapts to compute anisotropic superfluidity.","marker":"[60]"},{"why":"Defines the collapse criterion for quasi-2D dipolar systems used to mark the collapse region in the phase diagrams.","marker":"[65]"}],"fun_headline_variants":["Tilted polarization splits critical point into three lines","Three critical lines emerge when dipolar BEC polarization tilts","In-plane tilt deforms supersolid lattices and direction-dependent superflow","Polarization tilt radically alters dipolar phase diagram","Anisotropic superfluid from tilted polarization in planar dipolar gas"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation compares only three ansätze—homogeneous, stripes fixed along the in-plane polarization, and hexagonal lattices deformed only along that same axis—so if some oblique, rotated, or square lattice has lower energy anywhere in the (ρ, as/add) plane, the reported phase boundaries and even the critical-line topology could change.","fun_headline_variants_meta":{"raw":{"variants":["Tilted polarization splits critical point into three lines","Three critical lines emerge when dipolar BEC polarization tilts","In-plane tilt deforms supersolid lattices and direction-dependent superflow","Polarization tilt radically alters dipolar phase diagram","Anisotropic superfluid from tilted polarization in planar dipolar gas"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000206,"raw_usage":{"total_tokens":1393,"prompt_tokens":937,"completion_tokens":456,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":553,"completion_tokens_details":{"reasoning_tokens":371}},"tokens_in":553,"tokens_out":456,"duration_ms":4776,"temperature":1.0,"reasoning_tokens":371,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T19:48:33.114783+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform an unconstrained ground-state search on the same mean-field functional, allowing arbitrary lattice symmetries, oblique wavevectors, and stripe orientations not aligned with the in-plane polarization. If any pattern outside the three ansätze has lower energy in the regions the paper assigns to compressed hexagonal, stripes, or stretched honeycomb—especially near the two triple points—then the three-critical-lines topology is not the complete ground-state diagram.","supporting_citations":[{"cited_title":"Raghunandan, C","cited_arxiv_id":null,"evidence_quote":"Defines the collapse criterion for quasi-2D dipolar systems used to mark the collapse region in the phase diagrams."}],"review_version":1}