{"id":"3cf1744e-4772-4d43-933d-44541c1e294e","arxiv_id":"2507.02505","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"In a 2D two-component polariton model, increasing spin-orbit coupling turns quantum droplets into vortices, and high interactions re-form ordered droplet lattices.","lead":"A numerical study of two-component polariton condensates shows that increasing spin-orbit coupling can switch the condensate from smooth self-bound droplet states to vortex states, while stronger interactions favor ordered droplet lattices. The result maps a potential control knob for topological and self-bound states in photonic quantum fluids.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The droplet-to-vortex transition is inferred from early-time transient snapshots that the text itself says are not stationary solutions, so the central phase-control claim lacks stability evidence.","rationale":"The reader's weakest assumption concerns the conservative reduction from the driven-dissipative polariton model. That concern is important for experimental relevance, but it does not threaten the claim within the reduced model itself. The more load-bearing issue is internal: the paper's own text says the reported states are transient and not stationary solutions, while the conclusion and abstract claim stable vortex/droplet phases. A transition between phases cannot be established by early-time snapshots without long-time evolution or stationary-state computation. This reinforces the reader's overall CONDITIONAL verdict rather than changing it, since the missing evidence is exactly the kind of verification needed before the phase diagram can be accepted. I did not find grounds to reject the conservative-model results outright; the numerical method is standard and the model is plausible, but the central claim is not yet supported at the claimed level. The reader's rationale already noted missing convergence and long-time stability, so there is partial agreement, though the reader's formal weakest-assumption statement points elsewhere.","tokens_in":12490,"tokens_out":4053,"duration_ms":55529,"concrete_test":"Re-run the Fig. 5 setup (g = 100, delta_g = 500; sigma = 0.05, 2, 11) using the same split-step method, but evolve Eq. (7) to t ~ 1000 while monitoring the conserved norm N, the angular momentum L_z = integral psi* (x d_y - y d_x) psi dx dy, and the phase winding around every density minimum at regular intervals. Separately, compute stationary states by imaginary-time propagation or Newton relaxation seeded with the sigma = 11 snapshot. If the vortex core and winding survive past t = 1000 and match a stationary solution, the transition claim is supported; if the core decays to a droplet or the winding number changes, the central claim is a transient artifact.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section II explicitly states that all presented numerical results are 'early-time dynamical states' and that the displayed vortex and droplet structures 'are not stationary solutions, but transient, self-organized configurations.' Yet the conclusion repeatedly calls these 'stable vortex states' and 'stable droplet configurations,' and the central claim treats the spin-orbit strength sigma as a control parameter for a droplet-to-vortex phase transition. A phase transition requires the states to persist or correspond to attractors/stationary solutions of Eq. (7); single-time snapshots from a Gaussian initial condition do not establish this. In particular, the vortex rings reported for sigma = 11 in Fig. 5(c) could be intermediate structures that later shed angular momentum or relax into the droplet state, which would invalidate the claimed transition. No convergence tests, evolution timescales, norm-vs-time traces, angular-momentum diagnostics, or winding-number measurements are provided, so the reader cannot distinguish a genuine topological phase from a transient density dip. This missing stability evidence is the weakest load-bearing link in the paper's central argument.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a two-component (spin-1/2) polariton condensate in 2D, modeled by a Gross-Pitaevskii equation with photonic spin-orbit coupling and Lee-Huang-Yang-type quantum fluctuations. Starting from a driven-dissipative reservoir model, the authors impose a critical-balance condition and a strong-interaction inequality to reduce the system to a closed conservative GPE, Eq. (7). They then perform real-time split-step Fourier simulations from real Gaussian initial states and report density and phase profiles for varying intra- and inter-component interaction strengths g and δg and spin-orbit strength σ. The central claimed result is a droplet-to-vortex transition driven by increasing σ, with high interaction strengths producing ordered droplet lattices.","tokens_in":12742,"tokens_out":6904,"duration_ms":77787,"significance":"If substantiated, the claimed control of droplet versus vortex phases via the photonic spin-orbit strength σ would be a useful extension of quantum-droplet physics to spin-orbit-coupled polariton systems, and the numerical exploration of large parameter ranges is a step toward that goal. The paper is honest in one important respect: it explicitly labels the displayed structures as early-time, non-stationary configurations, and it does not fit parameters to pre-selected output patterns. However, the manuscript's central physical claim is currently supported only by transient snapshots, and it is weakened by internal sign inconsistencies and an unvalidated reduction from the open reservoir model. The significance will only be established after the stability and consistency issues are addressed.","major_comments":[{"comment":"The central claim that increasing σ drives a transition from a self-bound droplet to stable vortex states is not established by the evidence shown. The text after Eq. (11) states that the presented results are 'early-time dynamical states' and 'not stationary solutions, but transient, self-organized configurations,' yet the Conclusion repeatedly calls the σ = 11 configurations 'stable vortex states' and describes a 'topological transition.' No long-time evolution traces, norm-versus-time plots, angular-momentum diagnostics, winding-number measurements, or convergence tests in time step and grid resolution are provided. The vortex rings in Fig. 5(c) could be transient structures that later relax into the droplet state. Please provide long-time simulations (e.g., peak density and norm versus time, phase singularities tracked over time) and resolution/convergence checks for at least the representative parameters of Figs. 5 and 7, and either use 'stable' only where the dynamics justify it or reformulate the claim as applying to transient self-organized structures.","section":"Section II, Figs. 5-7 and Conclusion"},{"comment":"There is a sign inconsistency between the working GPE and the chemical-potential expression used for Fig. 1. Substituting the plane-wave ansatz ψ_l = A_l e^{-iμt} into Eq. (7) gives μ = (δg/2)(A_1^2 + A_2^2) + (-1)^l g(A_l^2 - A_{3-l}^2) - sqrt(σ_l/σ_{3-l}) p log p, whereas Eq. (10) has a minus sign on the g-term and a plus sign on the LHY term. In addition, the SOC derivative operator changes from -σ(∂/∂x - i(-1)^l ∂/∂y)^2 in Eq. (1) to +σ(∂/∂x + i(-1)^l ∂/∂y)^2 in Eq. (7). These signs cannot all be correct. Please reconcile the model equations, the chemical-potential formula, and the code; specify which sign convention was actually implemented in the simulations, since Fig. 1 and the phase/density dynamics depend on it.","section":"Section I, Eq. (7) vs. Eq. (10)"},{"comment":"The reduction from the open reservoir model to the conservative GPE is a load-bearing step but is not validated. The critical-balance condition (4) and the inequality (6) are stated, but no parameter values are given to show that the simulated regime satisfies g(|ψ_1|^2 + |ψ_2|^2) >> g_R n_l, g_R n_{3-l}, nor is any comparison made with the full Eqs. (1)-(2). In actual polariton condensates, pumping and decay typically dominate, so the observed droplet and vortex structures may be properties of a different closed system rather than of photonic polariton condensates. Please either perform simulations of the full driven-dissipative model (or of a controlled approximation with small gain/loss) and show that the reported patterns are robust, or explicitly restrict the conclusions to the conservative GPE without claiming direct polariton experimental relevance.","section":"Section I, Eqs. (1)-(7)"},{"comment":"The manuscript does not report the numerical setup needed to assess the reliability of the results: domain size, grid points, time step, integration time, boundary conditions, or any convergence study. The only statement is that parameters were 'adjusted based on the interaction parameters.' Given that the central results are complex self-organized patterns, it is essential to demonstrate that the observed droplet lattices and vortex rings are not numerical artifacts. Please provide a table of representative numerical parameters and at least one convergence test (e.g., doubling spatial resolution and halving time step for the parameters of Figs. 5 and 7).","section":"Section II, numerical methods and Fig. 7"}],"minor_comments":[{"comment":"The text repeatedly refers to Eq. (3) when discussing the conservative GPE, but Eq. (3) is the reservoir steady-state condition; the reduced GPE is Eq. (7). Similarly, the Gaussian initial state is called 'Eq. (7)' in Section II, but it is Eq. (11). Please correct all equation cross-references.","section":"Section I and II, cross-references"},{"comment":"The caption contains typos: 'for gδg = 500 and g = 100' should read 'for δg = 500 and g = 100', and 'σ − 1.5' should read 'σ = 1.5'.","section":"Fig. 8 caption"},{"comment":"The sentence 'The density distributions (|ψ1|2, |ψ2|2) remain uniform, suggesting a near-ground-state configuration with minimal influence from spin-orbit coupling' is repeated verbatim twice; please delete the duplicate.","section":"Section II, after Fig. 2"},{"comment":"References [60] and [61] are identical (Madimabe, Tabi, Tiofack, Kofané, Phys. Rev. B 107, 184502 (2023)); please merge or renumber and ensure each cited item is distinct.","section":"References"},{"comment":"The phase profiles are described qualitatively as 'spiral phase windings,' but no color bar, phase range, or winding-number label is provided, making it difficult to verify the sign and magnitude of the topological charge from the figures.","section":"Figs. 3 and 6"}],"recommendation":"major_revision","confidential_remarks":"The core idea is testable and the parameter exploration is potentially useful, but the present evidence does not yet support the claimed phase transition. I recommend major revision rather than rejection because the missing long-time data, convergence checks, and sign reconciliation are all feasible within the manuscript's scope. If the authors cannot reproduce the transition in longer-time simulations or cannot resolve the Eq. (7)/Eq. (10) inconsistency, the manuscript should not be accepted. I also note that the paper leans heavily on the authors' own prior works (e.g., Refs. [5], [15]-[17], [26], [65]); the novelty relative to Ref. [65] should be clarified."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know first: this is a parameter-sweep paper on a known two-component GPE with photonic spin-orbit coupling and LHY correction. The genuinely new part is the transient dynamics scan over sigma, g, and delta-g, and the droplet-to-vortex trend visible in the figures looks plausible. But the paper's own methods section says the displayed states are early-time transients, not stationary solutions, while the conclusion repeatedly calls them stable vortex states and stable droplet configurations. That mismatch is the load-bearing weakness.\n\nThe paper does a few things right. It is transparent about the transient nature of the snapshots in Section II, it explicitly cites the companion stationary-state paper [65], and it does not fit parameters to data—the reported patterns emerge from the simulations. There is no circularity problem here.\n\nThe stress-test note is on target. A phase-transition claim needs persistence or attractor evidence; none is supplied. No convergence tests, no norm-vs-time traces, no angular-momentum or winding-number diagnostics. The vortex rings at sigma = 11 could be intermediate structures that later relax into a droplet. The conclusion overstates what the evidence supports. There is also a sign mismatch between Eq. (7) and Eq. (10): the LHY term flips sign, and the g-dependent term in the chemical potential has the opposite sign relative to the equation of motion. These are mechanical errors, but they need fixing before the numerics can be trusted. The reduction to a conservative closed system via the critical balance condition is a standard move, but it is not validated against any dissipative simulation, so the polariton relevance is conditional rather than demonstrated.\n\nWho is this for: people actively working on SOC quantum droplets in polariton or BEC systems. They will find the parameter sweep a useful map, but the central claim of a robust droplet-vortex phase transition is not yet supported. A serious referee should engage, but the paper needs revision: stability diagnostics, sign corrections, and conclusions that match the transience disclaimer. I would not cite the stability claim as established, but the paper deserves a proper review rather than a desk rejection.","headline":"A numerically plausible but overclaimed droplet-to-vortex trend in a known polariton model; the stability evidence is missing and the text contradicts its own transience disclaimer.","tokens_in":13267,"tokens_out":3099,"would_cite":false,"duration_ms":35660,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Raising the photonic spin-orbit coupling strength switches a binary polariton condensate from a self-bound droplet into vortex states, while strong atomic interactions lock it into ordered droplet lattices.","keywords":["spin-orbit coupling","quantum droplets","vortex states","exciton-polariton condensates","Lee-Huang-Yang correction","two-component Gross-Pitaevskii equation","droplet lattices","two-dimensional quantum fluids"],"falsifier":"Re-solve the full open system, Eqs. (1)-(2), with finite reservoir density, nonzero polariton decay $\\gamma_c$, and pump rate $P_m$ at the same parameter values as Figs. 2-8, and check whether the droplet-to-vortex transition and the ordered droplet lattices survive; if these structures appear only when $R n_m/2 = \\gamma_c$ holds exactly, the central claim collapses. A complementary experiment would image the polariton density and phase in a planar microcavity with measurable TE-TM splitting, varying the effective spin-orbit strength at fixed interactions to look for the predicted switch from a vortex-free droplet to spiral phase windings.","tokens_in":12307,"feed_emoji":"🌀","tokens_out":13704,"duration_ms":125955,"temperature":0.7,"pith_summary":"This paper sets out to show that the strength of photonic spin-orbit coupling acts as a switch between two kinds of self-organization in a two-component (2+1)-dimensional exciton-polariton condensate: a self-bound quantum droplet at low coupling, and quantized vortex states once the coupling is raised. It also claims that strong atomic interactions convert the vortices into ordered lattices of quantum droplets, with the inter-component interaction playing a secondary role in nucleating and then deforming the vortices. These results come from numerically solving the coupled Gross-Pitaevskii equation with a Lee-Huang-Yang quantum-fluctuation correction, in the conservative limit where pumping and decay are taken to cancel exactly. If the claims hold, spin-orbit coupling becomes a control parameter for selecting the topological phase of a low-dimensional quantum fluid in polariton and cold-atom settings.","feed_headline":"Spin-orbit strength flips polariton droplets into vortices","feed_subtitle":"One tunable knob selects between droplets, vortex rings, and droplet lattices in a 2D condensate.","key_machinery":"The load-bearing object is the reduced conservative two-component Gross-Pitaevskii equation (Eq. (7) in the text) for spin-up and spin-down polariton fields $\\psi_1$, $\\psi_2$. Three terms carry the physics: the mean-field interaction terms $(\\delta g/2)(|\\psi_l|^2 + |\\psi_{3-l}|^2) + (-1)^l g(|\\psi_l|^2 - |\\psi_{3-l}|^2)$; the Lee-Huang-Yang quantum-fluctuation correction $-p \\log p$ weighted by $\\sqrt{\\sigma_l/\\sigma_{3-l}}$, with $p = (\\sigma_1 |\\psi_1|^2 + \\sigma_2 |\\psi_2|^2)/(2\\sqrt{\\sigma_1 \\sigma_2})$, whose repulsive beyond-mean-field pressure balances the mean-field attraction and stabilizes self-bound droplets; and the photonic spin-orbit coupling operator $\\sigma(\\partial_x - i(-1)^l \\partial_y)^2 \\psi_{3-l}$, which couples the two components and generates the phase winding that turns droplets into vortices. The argument also depends on the reduction step that removes gain and loss: at the steady reservoir density with $R n_m/2 = \\gamma_c$ the gain-loss term vanishes and the system becomes conservative. The split-step Fourier method propagates the Gaussian inputs in real time, and the emergent structures are read from density profiles and phase plots.","core_discovery":"The central claim is that the competition among photonic spin-orbit coupling, mean-field interactions, and Lee-Huang-Yang quantum fluctuations decides which phase a binary polariton condensate occupies: a uniform condensate, a self-bound droplet, a vortex-bearing state, or a lattice of droplets. Starting from a dissipative two-component Gross-Pitaevskii model with an excitonic reservoir, the authors impose the critical balance condition $R n_m/2 = \\gamma_c$ and neglect reservoir-mediated interactions, leaving a closed conservative equation for the two spin components with the LHY term $-\\sqrt{\\sigma_l/\\sigma_{3-l}}\\, p \\log p$. In real-time simulations from Gaussian initial states with no imposed phase winding, a stable droplet without phase singularities forms at low spin-orbit strength ($\\sigma = 0.05$); increasing $\\sigma$ to 2 and then 11 splits the droplet, develops spiral phase windings, and produces vortex rings, which the authors interpret as the spin-orbit term injecting quantized angular momentum. Increasing the inter-component interaction $\\delta g$ at fixed $\\sigma$ first nucleates vortices and then deforms them until the spiral fringes dissolve, while at high spin strength and fixed $\\delta g$, raising the intra-component interaction $g$ from 20000 to 50000 yields progressively more ordered droplet lattices with fewer phase singularities.","pith_inferences":["A natural consequence the paper does not spell out: because the TE-TM splitting in a microcavity is set largely by the cavity structure, the effective spin-orbit strength could be engineered to dial between droplet and vortex phases in a single device without retuning interactions.","Since the simulations start from phase-winding-free Gaussians, the vortices must be nucleated by the spin-orbit term itself; a direct test would be to track the total winding number as $\\sigma$ grows and see whether it increases continuously or in integer jumps.","The paper's own closing remarks concede that reservoir dynamics and gain-loss imbalance are beyond its scope; the sharpest extension would be to re-add pump and decay and check whether the droplet and vortex structures survive, since real polariton condensates are open, driven-dissipative systems."],"forward_implications":["The spin-orbit strength $\\sigma$ acts as a genuine control parameter: the same interaction parameters that give a vortex-free self-bound droplet at $\\sigma = 0.05$ produce vortex rings at $\\sigma = 11$, so the condensate topology can be switched without retuning the atomic interactions.","The reported structures are transient, self-organized configurations that emerge dynamically from Gaussian inputs carrying no initial phase winding, not stationary solutions; the paper explicitly leaves stationary-state dynamics to a companion study.","Raising the inter-component interaction $\\delta g$ at fixed $\\sigma$ first nucleates vortices and then dissolves their spiral fringes, so $\\delta g$ controls the deformation as well as the creation of topological defects.","Increasing the intra-component interaction $g$ at high $\\sigma$ makes droplet lattices more ordered and suppresses phase singularities, which the authors read as stronger interactions stabilizing coherence in the self-trapped regime.","If the claims hold, the tunable droplet-vortex-droplet-lattice sequence offers a route to vortex-based precision sensing and to studying superfluid-to-self-trapped transitions in low-dimensional quantum fluids."],"supporting_citations":[{"why":"Supplies the starting coupled Gross-Pitaevskii model with photonic spin-orbit coupling and LHY correction that the paper reduces and solves.","marker":"[61]"},{"why":"Source of the parameter $p$ expressing the Lee-Huang-Yang quantum-fluctuation correction used in the equation of motion.","marker":"[62]"},{"why":"Provides the steady-state reservoir condition and critical balance treatment used to drop the gain-loss terms.","marker":"[63]"},{"why":"Earlier conservative modeling of polariton condensates that the reduction procedure is said to align with.","marker":"[31]"},{"why":"The split-step Fourier method used for the real-time numerical propagation of the Gross-Pitaevskii equation.","marker":"[64]"},{"why":"Companion study of stationary-state dynamics that the paper distinguishes from its transient self-organized states.","marker":"[65]"},{"why":"Establishes that quantum fluctuations stabilize self-bound droplets in two-component condensates, the mechanism the simulations invoke.","marker":"[66]"},{"why":"Dipolar-condensate context for LHY-stabilized droplets against which the reported droplet behavior is compared.","marker":"[67]"}],"fun_headline_variants":["SOC knob dials condensate from droplets to vortex rings","Interplay of SOC and quantum fluctuations selects polariton phases","One parameter sweeps binary condensate from droplet to vortex lattice","Spin-orbit coupling and LHY interactions settle condensate phase","Vortex droplets emerge from photonic spin-orbit competition"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that pumping, decay, and reservoir interactions can be discarded completely: the reduction requires the gain to cancel the loss exactly, $R n_m/2 = \\gamma_c$, and the polariton interactions to dominate the reservoir interactions, $g(|\\psi_1|^2 + |\\psi_2|^2) \\gg g_R n_l, \\bar{g}_R n_{3-l}$, conditions the paper states and then relies on, even though in real polariton microcavities gain and loss usually dominate.","fun_headline_variants_meta":{"raw":{"variants":["SOC knob dials condensate from droplets to vortex rings","Interplay of SOC and quantum fluctuations selects polariton phases","One parameter sweeps binary condensate from droplet to vortex lattice","Spin-orbit coupling and LHY interactions settle condensate phase","Vortex droplets emerge from photonic spin-orbit competition"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000232,"raw_usage":{"total_tokens":1529,"prompt_tokens":1022,"completion_tokens":507,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":424}},"tokens_in":638,"tokens_out":507,"duration_ms":5958,"temperature":1.0,"reasoning_tokens":424,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:27:40.595366+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-solve the full open system, Eqs. (1)-(2), with finite reservoir density, nonzero polariton decay $\\gamma_c$, and pump rate $P_m$ at the same parameter values as Figs. 2-8, and check whether the droplet-to-vortex transition and the ordered droplet lattices survive; if these structures appear only when $R n_m/2 = \\gamma_c$ holds exactly, the central claim collapses. A complementary experiment would image the polariton density and phase in a planar microcavity with measurable TE-TM splitting, varying the effective spin-orbit strength at fixed interactions to look for the predicted switch from a vortex-free droplet to spiral phase windings.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the starting coupled Gross-Pitaevskii model with photonic spin-orbit coupling and LHY correction that the paper reduces and solves."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the parameter $p$ expressing the Lee-Huang-Yang quantum-fluctuation correction used in the equation of motion."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the steady-state reservoir condition and critical balance treatment used to drop the gain-loss terms."},{"cited_title":"Sakaguchi, B","cited_arxiv_id":null,"evidence_quote":"Earlier conservative modeling of polariton condensates that the reduction procedure is said to align with."},{"cited_title":"Abdolabadi, A","cited_arxiv_id":null,"evidence_quote":"The split-step Fourier method used for the real-time numerical propagation of the Gross-Pitaevskii equation."},{"cited_title":"Sanjay, S","cited_arxiv_id":null,"evidence_quote":"Companion study of stationary-state dynamics that the paper distinguishes from its transient self-organized states."},{"cited_title":"Chomaz, S","cited_arxiv_id":null,"evidence_quote":"Dipolar-condensate context for LHY-stabilized droplets against which the reported droplet behavior is compared."}],"review_version":1}