{"id":"d95d4eb0-281f-43fd-b036-ff13fcf74600","arxiv_id":"2607.25558","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"QMC equations of state for 39K Bose-Bose mixtures in the 3D-to-2D crossover show the 2D LHY functional is only reliable for very tight confinement; the paper's QMC-based fitted 2D functionals cover the crossover.","lead":"Quantum Monte Carlo simulations of ultracold 39K Bose-Bose mixtures squeezed into quasi-2D sheets show that the widely used 2D Lee-Huang-Yang (LHY) functional only matches the simulations at very strong confinement. The paper builds QMC-based density functionals for the 3D-to-2D crossover and finds that modeling droplets or vortices with the 2D LHY functional outside a narrow window can give qualitatively wrong answers.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Droplet/vortex predictions depend on an untested polynomial EoS at low density; bulk departure from 2D LHY is less exposed.","rationale":"The reader's weakest assumption identifies exactly the same bridge between the DMC bulk data and the paper's droplet/vortex predictions. I agree with it. The manuscript's own evidence supports the narrow-LHY claim through raw bulk DMC (Fig. 3), but the abstract's droplet/vortex statements are computed through the polynomial fit. The text provides no error bars on DMC quantities and no sensitivity check against alternative fit forms; the large-N validation is circular, as the functional was fit to the same EoS. A fit-form sensitivity check is the cheapest, most direct way to determine whether the concern lands. A direct DMC simulation of a trapped droplet would be stronger but much more expensive. Therefore the existing CONDITIONAL verdict is appropriate and needs no adjustment.","tokens_in":10674,"tokens_out":7793,"duration_ms":82991,"concrete_test":"Refit the DMC EoS for f=0.15 and f=0.05 with an alternative form that includes the expected 2D low-density log term—e.g., E/N = αρ + βργ + δ ρ² ln(ρ/ρ0)—or with a monotone spline through the raw DMC points, and rerun the DFT of Figs. 2, 4 and 5. If the central density, E/N versus N, or the vortex-state ordering changes by more than the DMC statistical uncertainty, the polynomial ansatz is load-bearing; if profiles/energies are unchanged, the concern is retired.","verdict_should_be":"UNCHANGED","load_bearing_attack":"All droplet and vortex results (Figs. 2, 4, 5) are generated by inserting the per-confinement fit E/N = αρ + βργ (Table VI) into the LDA DFT, Eq. (2). The fit is anchored to DMC energies near the bulk equilibrium, while a self-bound droplet edge and a vortex core reach arbitrarily small ρ. There the ansatz gives E/N → αρ (negative linear), with no ρ² ln ρ log term that the 2D LHY EoS (Eq. 11) and the true 2D EoS possess. Since γ ≈ 1.02–1.26 and α < 0, the low-density extrapolation is uncontrolled and could materially alter the outermost shells of the profile, the droplet energy, and core energetics. The large-N 'validation' (saturation density → bulk DMC value) is circular: the functional is fit to that same bulk EoS. Thus the quantitative magnitude of the claimed departure and the Fig. 5 vortex-state ordering rest on an ansatz that is not tested against alternative forms or direct trapped-droplet DMC. The bulk EoS comparison in Fig. 3 is less exposed, so the 'narrow LHY' effect itself is not overturned, but the droplet/vortex consequences are conditional.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports diffusion Monte Carlo (DMC) equations of state for ultradilute 39K Bose-Bose mixtures in a one-dimensional harmonic confinement, spanning the crossover from three to two dimensions. Using two different interaction potentials with the same s-wave scattering length and effective range, the authors establish that the EoS is universal in these two parameters down to a_z ≈ 9–10 a_11 (f = 0.05), and that below this confinement two-parameter universality breaks down. The DMC EoS is fitted to E/N = αρ + βργ and used to construct a 2D QMC density functional within the local density approximation; this functional is then applied to compute droplet density profiles, energies, and vortex states. The central conclusions are that the 2D LHY functional is applicable only in a very narrow range of very tight confinements, and that the 2D QMC functional provides more reliable droplet and vortex predictions than 2D LHY in the quasi-2D regime.","tokens_in":10839,"tokens_out":6268,"duration_ms":68086,"significance":"If the results hold, the paper provides a practical QMC-based density functional for quasi-2D Bose-Bose droplets and directly challenges the common practice of using the 2D LHY functional in crossover and quasi-2D settings. The claim that the range of validity of 2D LHY is very narrow has important consequences for predictions of droplet profiles, energies, and vortex stability, and is supported by the paper's two-potential universality check and by independent PIMC results. The manuscript's strengths include fully tabulated model potentials, a systematic test of universality with two potential families, and a direct comparison with prior PIMC and quasi-2D LHY studies. The main vulnerability is the untested low-density extrapolation of the fitted EoS used to build the density functional, which directly affects all droplet and vortex predictions.","major_comments":[{"comment":"All droplet and vortex results (Figs. 2, 4, 5) are obtained from the fitted EoS E/N = αρ + βργ, with parameters in Table VI. The fits are anchored to DMC data near the bulk equilibrium density, but a self-bound droplet edge and a vortex core sample arbitrarily small ρ. There the ansatz yields E/N ≈ αρ (linear and negative), whereas the 2D LHY EoS (Eq. 11) and the expected 2D equation of state contain a nonanalytic ρ² ln ρ term. Since γ > 1 and α < 0 in all cases, the low-density behavior of the functional is an uncontrolled extrapolation. No sensitivity test against alternative fit forms (e.g., including a logarithmic term) or against direct trapped-droplet DMC is provided. The Fig. 5 vortex-state ordering and the quantitative droplet energies in Fig. 4 therefore rest on an ansatz whose low-density shape is not validated. This is the main load-bearing gap in the paper.","section":"Results / Eq. (2) / Table VI"},{"comment":"The statement that saturation densities obtained from the 2D QMC functional 'reach the quasi-2D bulk DMC value, demonstrating functional’s reliability for large droplets' is circular. The functional was constructed by fitting to that same bulk DMC EoS; in the large-N limit the DFT necessarily recovers the minimum of the input EoS. This check confirms internal consistency, but it does not independently validate the functional for inhomogeneous configurations, where the low-density region matters. An independent validation, such as a DMC calculation of a trapped droplet or a direct comparison with PIMC for the same mixture, would be needed to support the droplet-level conclusions.","section":"Results / Fig. 3"},{"comment":"No statistical error bars are reported for the DMC energies per particle or for the fitted parameters α, β, γ in Table VI. The claim that below f = 0.05 two scattering parameters are 'not sufficient' is based on the difference between POT1 and POT2 results in Fig. 1. Without error bars or fit residuals, the reader cannot judge whether the apparent POT1/POT2 discrepancy at f = 0.01 is statistically significant. Given that this is a load-bearing part of the universality claim, the authors should provide at least representative error bars on the EoS data points and on the fitted parameters.","section":"Fig. 1 / Table VI / DMC implementation"}],"minor_comments":[{"comment":"Typographical errors: 'Paramatres' should be 'Parameters', 'unists' should be 'units'. Please proofread the End Matter.","section":"End Matter / Tables I–V"},{"comment":"The sentence 'The liquid distribution in the direction of squeezing follows has the single-particle gaussian shape' is grammatically awkward. Also, the inset would benefit from a definition of the plotted quantity (probability density versus z/a_z, presumably).","section":"Fig. 1 caption"},{"comment":"The confinement is given only as a_z in µm. Since the text uses the dimensionless factor f = a_z/0.639 µm, it would help to list f alongside a_z in the table, or at least define the corresponding f values in the caption.","section":"Table VI"},{"comment":"The top and bottom panels show different squeezings, but the atom number N is only stated in the text for one case. Please label or state N explicitly for each curve so the comparison is unambiguous.","section":"Fig. 2"},{"comment":"Ref. [33] is cited as Phys. Rev. A 109, 013313 (2014); PR A volume 109 corresponds to 2024, so the year appears to be a typo. Also, Ref. [35] and [39] would benefit from article titles, which are currently omitted.","section":"References"},{"comment":"The comparison with the PIMC results of Ref. [40] is indirect: that work uses a symmetric mixture, different scattering parameters, and no finite-range effects. The agreement with 3D functionals is suggestive but not a quantitative test of the proposed 2D QMC functional. This should be stated more cautiously.","section":"Comparison with Ref. [40]"}],"recommendation":"major_revision","confidential_remarks":"The bulk EoS part of the paper is solid and the central universality claim is credible; the two-potential check is a good methodological choice. However, the paper's headline droplet and vortex conclusions are built on an EoS fit whose low-density behavior is not tested, and the large-N 'validation' is circular. This is fixable with additional sensitivity tests (e.g., alternate fit forms, direct trapped-droplet DMC, or a stronger PIMC comparison), so major revision rather than rejection is appropriate. The absence of error bars on the DMC data is a separate concern that should be addressed during revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing you should know: this paper makes a solid case that the 2D LHY functional, used in a lot of quasi-2D droplet work, only works at very tight confinement, and it offers a cheap QMC-based alternative. The bulk equation-of-state evidence is the strong part. The weak part is the bridge from that bulk EoS to the droplet and vortex predictions.\n\nWhat is new: a systematic DMC study of the 3D-to-2D crossover for 39K mixtures, using two potential models with the same a and r_eff, showing universality down to a_z ≈ 9–10 a_11 and building per-confinement 2D functionals. The two-potential check is a sound methodological choice, and the qualitative departure from 2D LHY is corroborated by independent work (Refs. 31, 32, 40). The paper is honest about its limits: it explicitly says it cannot compare directly to Ref. [40], notes that below f = 0.05 two scattering parameters are insufficient, and shows 3D functionals fail under the strongest squeezing.\n\nThe soft spot is the DFT step. All droplet and vortex results come from inserting a three-parameter fit E/N = αρ + βργ into the local-density DFT. The fit is anchored around bulk equilibrium densities, but droplet edges and vortex cores sample arbitrarily low densities, where the true 2D EoS has a logarithmic term that this ansatz lacks. The paper gives no sensitivity test against alternative fit forms and no error bars on any DMC quantity, so the factor-4 and factor-10 discrepancy claims carry unknown uncertainty. The large-N validation—saturation density reaching the bulk value—is circular, since the functional was fit to that same bulk EoS. These are real limitations, but they do not overturn the central effect: the bulk comparison in Fig. 3 is less exposed, and the narrow window of 2D LHY validity looks robust across several cross-checks.\n\nThis paper is for people working on quasi-2D droplets, vortices, and LHY functionals, and for anyone using 2D LHY in excitation or turbulence studies. It deserves a serious referee. The referees should push for error bars, an alternative low-density fit form, and a check on the vortex ordering, but the work is worth engaging with now.","headline":"A credible QMC-based case that 2D LHY functionals fail outside very tight confinement, but the droplet/vortex predictions lean on an untested low-density EoS ansatz.","tokens_in":11552,"tokens_out":1683,"would_cite":true,"duration_ms":16955,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["67.85.-d","03.75.Hh","05.30.Jp"],"model":"deepseek-v4-flash","headline":"2D LHY functional is accurate only for the tightest quasi-2D confinements","keywords":["Bose-Bose mixtures","quantum droplets","Lee-Huang-Yang correction","quasi-2D confinement","diffusion Monte Carlo","density functional theory","effective range","universality"],"falsifier":"A direct diffusion Monte Carlo simulation of a self-bound droplet at f = 0.15 with N ≈ 20,000 atoms, avoiding the density-functional bridge entirely, would settle whether the 2D QMC functional's predicted central density, energy, and vortex stability are correct; alternatively, measuring the low-density tail of the droplet's equation of state and comparing it with the fit E/N = αρ + βρ^γ would test the ansatz that carries the argument.","tokens_in":10405,"feed_emoji":"⚛️","tokens_out":5640,"duration_ms":53399,"temperature":0.7,"pith_summary":"The paper shows that the standard Lee-Huang-Yang density functional for 2D Bose-Bose droplets, the 2D LHY functional, is accurate only for very strongly squeezed quasi-2D systems, a much narrower range than commonly assumed. Using diffusion Monte Carlo at zero temperature with two different interaction potentials, the authors establish that equations of state are universal in terms of the s-wave scattering length and effective range down to a harmonic oscillator length of about 9–10 times the intraspecies scattering length, but diverge for tighter confinement. From these equations of state they construct a new 2D quantum-Monte-Carlo-based density functional for each confinement strength, and with it compute droplet density profiles, energies, and vortex stability. Their functional agrees with 3D functionals at moderate squeezing and approaches the bulk Monte Carlo results for large droplets, while the 2D LHY functional overestimates central densities and self-binding energies by large factors in that regime. If correct, the work provides a computationally cheap and systematically more reliable route to quasi-2D droplet physics and changes expectations for vortex studies.","feed_headline":"2D LHY functional applies only at very tight confinement","feed_subtitle":"Quantum Monte Carlo data provide a better density functional for droplets across the 3D-to-2D crossover.","key_machinery":"The key machinery is a confinement-specific 2D density functional built from zero-temperature diffusion Monte Carlo (DMC) equations of state. For each squeezing strength, the bulk energy per particle from DMC is fitted to the three-parameter form E/N = αρ + βρ^γ, with γ between about 1.02 and 1.26 for the tabulated cases, and this expression enters the DFT equation through E_int = ρ E/N. The universal range is established by comparing two interaction potential models (POT1 and POT2) matched to the same scattering length and effective range. The functional is then used, together with the local density approximation, to compute droplet profiles, liquid-drop energies, and vortex-core structure.","core_discovery":"The central quantitative finding is a universality boundary: for 39K Bose-Bose mixtures squeezed by a harmonic trap, equations of state computed with two different potential models that share the same scattering length and effective range agree down to a confinement length a_z ≈ 9–10 a_11 (f = 0.05), but below that the two-parameter description breaks down and details of the interaction potential matter. In the universal regime the paper builds a 2D QMC functional by fitting the DMC energy per particle to E/N = αρ + βρ^γ for each confinement. Using this functional inside a density-functional equation with the local density approximation, the paper shows that the widely used 2D LHY functional","pith_inferences":["The paper's fit form E/N = αρ + βρ^γ contains no logarithmic term that the exact 2D equation of state possesses; a natural extension would be to test whether including a log-density term changes droplet-edge properties or the vortex conclusions.","Because the functional was validated on large droplets by showing that saturation density approaches the bulk DMC value, the same quantity it was fitted to, an independent check such as a direct DMC simulation of a finite self-bound droplet would strengthen the case that the LDA-plus-fit bridge is faithful in low-density regions.","If the universality breakdown below a_z ≈ 10 a_11 is a general feature, analogous functional-construction efforts for other mixtures would need to map out their own universality windows before applying the 2D LHY functional.","The vortex-stability reversal between 2D LHY and QMC functionals suggests that experimental searches for vortices in droplets should target intermediate confinements, where the competing functionals give distinct density profiles and where the QMC functionals predict vortex states to be energetically favored."],"forward_implications":["If the paper is right, vortex and excitation studies of quasi-2D droplets relying on the 2D LHY functional need to be revisited in the moderate-confinement regime, where that functional can flip the sign of the vortex-energy comparison.","The 2D QMC functional provides a computationally efficient alternative that remains accurate across the crossover from 3D to 2D, matching 3D functional results at moderate confinement and the bulk DMC limit for large droplets.","For the tightest confinements (a_z below about 10 a_11), two scattering parameters are no longer sufficient; any theoretical or experimental analysis in that regime must account for details of the interatomic potentials.","The construction used here, computing the bulk equation of state with QMC, fitting a simple empirical form, and using it as a density functional, can be applied to other mixtures and to dipolar molecular droplets, where beyond-LHY effects are stronger."],"fun_headline_variants":["2D LHY valid only in tight squeeze","Quantum Monte Carlo functional beats 2D LHY in crossover","39K droplets: QMC functional spans 3D-2D crossover","Universality breaks below a_z≈9a_11 in 2D limit","2D LHY functional fails beyond tightest confinement"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"All droplet and vortex predictions are obtained by feeding the fitted form E/N = αρ + βρ^γ into a density functional and applying the local density approximation across the density profile; if that simple fit misrepresents the true equation of state at low densities or near equilibrium, the computed profiles, energies, and vortex-stability conclusions inherit the error.","fun_headline_variants_meta":{"raw":{"variants":["2D LHY valid only in tight squeeze","Quantum Monte Carlo functional beats 2D LHY in crossover","39K droplets: QMC functional spans 3D-2D crossover","Universality breaks below a_z≈9a_11 in 2D limit","2D LHY functional fails beyond tightest confinement"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000173,"raw_usage":{"total_tokens":1147,"prompt_tokens":807,"completion_tokens":340,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":551,"completion_tokens_details":{"reasoning_tokens":252}},"tokens_in":551,"tokens_out":340,"duration_ms":3797,"temperature":1.0,"reasoning_tokens":252,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T02:05:18.379406+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct diffusion Monte Carlo simulation of a self-bound droplet at f = 0.15 with N ≈ 20,000 atoms, avoiding the density-functional bridge entirely, would settle whether the 2D QMC functional's predicted central density, energy, and vortex stability are correct; alternatively, measuring the low-density tail of the droplet's equation of state and comparing it with the fit E/N = αρ + βρ^γ would test the ansatz that carries the argument.","supporting_citations":[],"review_version":1}