{"id":"c749aaa9-867f-4696-a288-afa8125b49ca","arxiv_id":"2501.07816","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A hybrid particle-continuum simulation of microswimmers in a nematic liquid crystal reproduces 2D living-nematic experiments and predicts 3D defect-loop dynamics controlled by winding profile, swimmer density, and reversal time.","lead":"This paper reports a computer model of bacteria swimming in a liquid crystal, mixing particle-based swimmers with a fluid dynamics model of the anisotropic liquid. It reproduces known 2D experiments and then predicts how 3D loop-shaped defects in the liquid crystal move and change when swimmers are present.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 3D loop-defect conclusions may be artifacts of hand-imposed director profiles: Fig. 4(b) shows even the swimmer-free loop relaxes quickly, so the winding-controlled dynamics need a passive-relaxation control.","rationale":"The reader's weakest assumption is essentially the concern I would stress-test: the 3D loop predictions depend on hand-constructed initial director fields that may not be representative or dynamically maintained. The 2D validation against spiral and C-pattern experiments is credible and independently grounds the hybrid method, so I do not object to the method itself. The 3D section is where the central design-rule claim lives, and there the evidence depends wholly on the Eq. (13) ansatz and on a small number of deterministic runs. The fact that a swimmer-free loop in Fig. 4(b) shrinks rapidly shows the initial field is not force-balanced, making it impossible to separate intrinsic winding-controlled active-nematic dynamics from relaxation of the imposed ansatz without a control. One might argue that all active-defect simulations initialize loops and let them evolve; however, the paper's claim is stronger, namely that the winding profile is a usable design rule for living-nematic experiments. Because swimmers align to the initial director and create active stress from their local distribution, the ansatz also preselects swimmer accumulation sites, which makes the winding-control interpretation more fragile rather than less. This is not a fatal flaw; it is a missing control. The reader's CONDITIONAL verdict is therefore appropriate: the 3D predictions should be treated as preliminary until either a passive-relaxation control or three-dimensional experimental confirmation shows that the prescribed β/α profiles are dynamically maintained and that the observed destinies are reproducible.","tokens_in":19534,"tokens_out":6464,"duration_ms":78738,"concrete_test":"Run a passive-relaxation control: initialize the exact wedge-twist and pure-twist loops of Figs. 4–6 with ρ0=0, record β(s,t) along the loop and the loop shape over the timescale used in the active runs, and determine how much the prescribed β/α profile drifts before swimmers would begin acting. If β changes substantially (e.g., by more than ~10% of π/2) before t≈50τ, the imposed winding profile is not dynamically maintained, and the measured loop dynamics are not controlled by the stated initial profile. As a complementary check, prepare the same loops by relaxing an initially straight disclination with the appropriate boundary anchoring and add swimmers only after the loop reaches a slowly evolving quasi-static state; if the expansion, shrinkage, or buckling in Fig.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central 3D claim is that loop destiny is controlled by the local winding profile encoded in Eq. (13). The evidence for this rests on initializing the loop directors by hand with prescribed β and α profiles, then letting the loop evolve with swimmers. These initial configurations are not force-balanced metastable states: Fig. 4(b) shows that at ρ0=0 the imposed wedge-twist loop shrinks rapidly under elastic forces. If the elastic relaxation timescale is comparable to or shorter than the timescale on which swimmers redistribute and generate active stress, then the β-profile at the time swimmers actually act differs from the prescribed profile. In that case, the observed expansion, shrinkage, and buckling would reflect relaxation from the hand-built initialization rather than an intrinsic property of a loop in a living nematic. The concern is compounded by the coupling: swimmers align to the local director field (Eq. 5) and the active stress is proportional to their local distribution, so an artificially planted splay/bend pattern directly seeds the swimmer accumulation that then drives the loop. Thus the conclusion that the winding profile dictates loop morphodynamics risks inheriting the initial condition rather than establishing a robust property of the system. A secondary but related weakness is that each parameter set is represented by a single deterministic trajectory without ensemble statistics, so run-to-run variability of these driven-relaxation processes is unknown.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a hybrid simulation method in which active Brownian particles (microswimmers) are coupled to a lattice-Boltzmann nematic hydrodynamics solver. The microswimmers are modeled as extensile point-force dipoles that align with the local director, while their active stress enters the nematic Navier-Stokes equation through a Gaussian kernel. The method is first validated against published quasi-2D experiments on spiral-patterned and C-patterned cells, reproducing dilute following, stable swirling, undulated jets, and a wavelength scaling law. The authors then apply the method to 3D disclination loops initialized with prescribed wedge-twist and pure-twist director profiles described by Eq. (13), and report that loop expansion, shrinkage, buckling, and tumbling depend on the loop's winding profile (β and α) as well as swimmer activity, size, and initial distribution. A final section adds stochastic swimming-direction reversals and predicts density pumping in a hybrid cell, including accumulation near homeotropic walls for long reversal times.","tokens_in":19983,"tokens_out":5715,"duration_ms":59750,"significance":"If the 3D predictions hold, the paper would provide a validated simulation platform for living nematics with particle-level resolution, extending prior continuum and particle-based approaches. The explicit benchmarks against independent experiments (spiral and C-pattern) are a strength, as is the demonstration of stabilized undulations and the wavelength scaling. The 3D loop-defect study goes beyond uniform-activity active nematics by letting the swimmer distribution be modulated by the nematic structure, which is a promising direction. However, the quantitative 2D agreement is only semi-quantitative (factor of about two in critical density; free critical concentration in the fit), and the 3D conclusions rest on hand-imposed initial director fields and single trajectories, so the predictive claims require additional controls.","major_comments":[{"comment":"The 3D loop-defect dynamics are initialized with hand-imposed director profiles that are not equilibrium states. The manuscript itself states that this configuration is inherently unstable and that at ρ0=0 the loop rapidly shrinks due to elastic forces (Fig. 4(b)). Because swimmers align to the local director (Eq. 5) and the active stress is proportional to their local distribution, an artificially planted splay/bend pattern can seed swimmer accumulation before the loop has relaxed. The central claim that winding profile dictates expansion, shrinkage, and buckling therefore requires a passive-relaxation control: either relax the bare loop for a time comparable to the swimmer redistribution time before introducing swimmers, or measure the β and α profiles at the time active stress becomes significant and show that they match the prescribed profiles. Without such a control, the observed dynamics may be a relaxation artifact of the initialization rather than an intrinsic property of the living nematic.","section":"III.B, Eq. (13), Fig. 4(b)"},{"comment":"Each parameter set is represented by a single deterministic trajectory, with no independent realizations or ensemble statistics. This is particularly concerning near threshold densities, e.g., the transition between shrinkage and expansion around ρ0=2.3×10^-3 in Fig. 4(b), where stochastic fluctuations (the Dr term in Eq. 5 and the random reversal process in III.C) could move a single realization across the boundary. The 'destiny' classifications (expansion, shrinkage, buckling, tumbling) and the density-pumping prediction should be supported by at least a few independent realizations, or the authors should explicitly state that the results are single-trajectory observations whose run-to-run variability has not been assessed.","section":"III.B, Figs. 4-7"},{"comment":"The quantitative validation is weaker than the text suggests. The simulated critical concentration for the C-pattern undulation is reported as 4.45×10^10 m^-2, compared with 2.36×10^10 m^-2 from experiment, a factor of about 1.9 discrepancy. In addition, the wavelength scaling λ ∝ 1/sqrt(ρ_s^j - ρ_sc^j) uses ρ_sc^j = 0.34 as a free fitting parameter (Fig. 2(h)). Since the abstract and conclusion describe the method as a 'faithful tool' validated against experiments, the authors should either qualify the 2D validation as semi-quantitative or discuss the origins of the factor-of-two threshold mismatch and the role of the free critical concentration in the scaling fit.","section":"III.A.2, Fig. 2(h)"}],"minor_comments":[{"comment":"There are typos in the text, for example 'protypical' in the first sentence of the Introduction and 'T o pological Defects' in the Section III.B header; please proofread the manuscript.","section":"Introduction and III.B headings"},{"comment":"The caption of Fig. 6(e) refers to the 'left side of the simulation box (y < Ny/2)' containing 'triple as many microswimmers' as the right side, while the text says the '+y half' contains 3Np/4 and the '−y half' contains Np/4; please reconcile the notation so the nonuniform initial distribution is unambiguous.","section":"III.B, Fig. 6(e) caption"},{"comment":"Reference [62] is cited as 'in preparation' to support the traveling-wave interpretation of the stabilized undulation; such an unpublished reference cannot be verified and should be updated or removed before publication.","section":"III.A.2, Ref. [62]"},{"comment":"The parameter l is used in figure captions and in Eq. (12) but is never explicitly defined; please state that g⊥ = l and give its value in the parameter list.","section":"II.B, Eq. (12)"},{"comment":"There is a garbled sentence in Section III.B: 'with trev ∼ 2 − 12 sin a nematic DSCG solution' appears to be missing a word; in addition, the expression for Πp contains tensor-index notation that is difficult to parse and should be checked for missing parentheses and consistent index ordering.","section":"II.B, passive stress tensor"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern is on target: the manuscript explicitly admits that the ρ0=0 loop is unstable and shrinks under elastic forces, yet no control separates the effect of the hand-imposed winding profile from the effect of relaxation of that profile. I would ask for a passive-relaxation control and at least a small ensemble of trajectories before treating the 3D loop-destiny predictions as design rules. The 2D benchmarks are a solid foundation and make the paper suitable for Soft Matter after revision, but the quantitative 'faithful tool' claim should be softened or supported with more discussion."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look, and worth a referee. The hybrid method is genuinely new: coupling active Brownian particles to a lattice Boltzmann nematic solver with a spatially distributed active stress is not something I've seen done this way. The 2D validation against the spiral and C-pattern experiments is credible, and the simulations actually reproduce the qualitative state transitions plus details like velocity profiles and the stabilized undulation. The factor-of-two mismatch in the critical density for the C-pattern is the kind of semi-quantitative agreement you'd expect from a coarse-grained model, and the wavelength fit with a free rho_c is a bit soft, but those are not fatal. Credit where due: the paper is honest about the threshold mismatch and includes a comment that 3D experimental validation is still needed.\n\nThe 3D loop-defect work is the interesting part. The claim that winding profile controls expansion, shrinkage, and buckling is plausible and the schematics are clear. But the stress-test concern about hand-imposed director fields has real substance. The initial loops in Eq. (13) are not metastable; the rho0=0 curves in Figs. 4(b), 5(b), and 6(b) show that the loops relax under elastic forces on the simulation timescale. So the passive-relaxation control is actually present--the paper does include swimmer-free runs. What is missing is not the control, but the statistical and dynamical evidence to separate relaxation from winding-driven effects. Each parameter set is a single deterministic trajectory, so we don't know run-to-run variability, and the loop's winding profile at the time swimmers act is not shown to match the prescribed profile. The swimmers align to the local director and generate active stress proportional to their local density, so the initial splay/bend structure directly seeds the swimmer accumulation; the claim that winding profile dictates destiny could be partly inherited from the initial condition. That's the core soft spot.\n\nThe citation pattern is fine. The self-citations are to the authors' own LBM method, which is appropriate, and the unpublished companion reference [62] is used for extra detail, not as load-bearing evidence. No code or data are released, which would help but is not a requirement for a methods paper of this type.\n\nBottom line: this deserves peer review. A good referee should push for ensemble statistics, a demonstration that the imposed loops are either metastable or that the active dynamics are robust to the initialization, and ideally a comparison with a uniformly-active control. The 2D validation alone makes it publishable somewhere reasonable; the 3D results could make it a more important paper if those concerns are addressed.","headline":"A useful hybrid ABP-LBM method for living nematics with credible 2D validation, but the 3D winding-profile claims rest on single trajectories from hand-built director fields and need a stronger control.","tokens_in":20403,"tokens_out":1981,"would_cite":true,"duration_ms":23830,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76A15"],"pacs":["83.80.Xz","47.63.mf"],"model":"deepseek-v4-flash","headline":"A hybrid simulation coupling active Brownian particles to a lattice-Boltzmann nematic model shows that the winding profile of a 3D loop defect dictates whether confined microswimmers expand, shrink, or buckle the loop.","keywords":["living nematic","active Brownian particles","lattice Boltzmann method","nematic liquid crystal hydrodynamics","3D topological defects","loop defect dynamics","microswimmer guidance","bacterial reversal time"],"falsifier":"Prepare a living nematic with a single, well-characterized 3D loop defect of known winding profile and track its lateral size over time across a range of bacterial densities; for example, a homeotropic-cell wedge-twist loop that shrinks at a density where Fig. 4(b) predicts expansion would falsify the claim. A second check is to measure bacterial density near the hybrid wall as a function of reversal time, since the claim predicts that for large $t_{\\mathrm{rev}}$ bacteria accumulate at bend-dominated regions, so observing splay-only accumulation regardless of reversal time would disprove it.","tokens_in":19346,"feed_emoji":"🦠","tokens_out":5994,"duration_ms":52510,"temperature":0.7,"pith_summary":"The paper develops a simulation method for living nematics, where bacteria or synthetic swimmers are dispersed in a liquid crystal, and claims that the method is faithful enough to explain existing 2D experiments and to predict new 3D behavior. The method couples active Brownian particle dynamics for the swimmers to a lattice-Boltzmann hydrodynamic model of the nematic director and flow. It reproduces stable swirling, polar jets, and their bend-induced undulation in patterned cells, with quantitative wavelengths and thresholds matching experiments. When applied to 3D loop defects, the simulations show that a loop's fate depends on its winding profile: +1/2 wedge windings and radial twist windings control expansion, shrinkage, and buckling, and the swimmer distribution feeds back on these dynamics. The paper also predicts that with stochastic swimming-direction reversals, bacteria do not necessarily accumulate in splay regions, challenging a common expectation.","feed_headline":"3D loop windings decide whether bacteria expand or shrink them","feed_subtitle":"A hybrid nematic simulation reproduces bacterial swirls and jets, then predicts loop fate from wedge and twist windings.","key_machinery":"The carrying object is the hybrid method itself: point-like active Brownian particles (overdamped translation, orientation coupling to the director and to the flow via a Jeffery/Bretherton torque, and orientational noise) whose active stress feeds back through an anisotropic-Gaussian kernel into a lattice-Boltzmann solver of the Beris-Edwards nematodynamic equations, including the $Q$-tensor, Landau-de Gennes free energy, and Navier-Stokes flow. For the 3D loop predictions, the load-bearing construction is the local director profile of a disclination winding, Eq. (13), with tangent $\\mathbf{t}$, rotation axis $\\boldsymbol{\\Omega}$, phase offset $\\alpha$, and twist angle $\\beta$; $\\beta=0$ or $\\pi$ gives wedge windings, $\\beta=\\pi/2$ pure-twist windings, $\\alpha=0$ radial twist, and $\\alpha=\\pi/2$ tangential twist. The argument proceeds by imposing these profiles as initial conditions and letting the coupled particle-nematic dynamics evolve the loop, then attributing expansion, shrinkage, and buckling to the active flow generated by the +1/2 wedge or radial twist windings.","core_discovery":"The central claim is that the morphodynamics and destiny of a 3D disclination loop in a living nematic are governed by the loop's local winding profile as well as by swimmer activity, effective size, and initial distribution. Using the director-field construction $\\mathbf{n} = \\cos(\\phi/2)\\mathbf{m} + \\sin(\\phi/2)(\\cos\\beta\\,\\mathbf{m}\\times\\mathbf{t} + \\sin\\beta\\,\\mathbf{t})$ (Eq. 13) to prescribe wedge and twist windings, the authors find that in a wedge-twist loop the +1/2 wedge drives expansion or shrinkage depending on its orientation, while in a pure-twist loop radial twist windings act like the wedge and other windings generate out-of-plane flows that buckle the loop. They also find that swimmers migrate to splay regions, which can accelerate loop motion and, when the +1/2 wedge forms a dense cluster, can eject child loops; with large swimmer spacing the reduced active stress fails to stabilize the loop. When swimmers reverse direction stochastically with characteristic time $t_{\\mathrm{rev}}$, the accumulation pattern changes, and for long $t_{\\mathrm{rev}}$ swimmers may accumulate where bend deformation dominates, so splay accumulation is not universal.","pith_inferences":["Because the method retains particle-level resolution, it could serve as a testbed for designing surface-anchor patterns that steer bacterial jets or stabilize defect loops in 3D, for example by combining patterned anchoring with time-varying reversal times.","The finding that nonuniform swimmer distributions can convert a pure-twist loop into a wedge-twist loop suggests that local activity gradients, not just global density, can rewire defect topology, a design handle the paper does not fully develop.","The predicted reversal-time-dependent accumulation could be tested with bacteria whose reversal frequency is controlled chemically or genetically, giving a direct experimental route to validate the model beyond reproduction of existing patterns.","The point-force treatment of swimmers likely underestimates excluded-volume and anchoring effects at high local densities; near dense clusters at +1/2 wedges, a squirmer-based model might change quantitative expansion rates while preserving the qualitative winding dependence."],"forward_implications":["In thin patterned cells, the hybrid method reproduces the experimentally observed transition from dilute director-following to dense stable circulation, and then to undulated swirling, with quantitatively captured dominant undulation wavelength.","On periodic C-patterns, the method captures polar jet formation in splay regions and the density-dependent bend undulation, including the $\\lambda \\propto 1/\\sqrt{\\rho_j^s - \\rho_j^{sc}}$ wavelength scaling and wave rupture at high density.","For 3D wedge-twist loops, the +1/2 wedge direction determines whether the loop expands (outward wedge, homeotropic cell) or shrinks (inward wedge, planar cell); swimmer clustering near the wedge can create child loops, while large swimmer spacing weakens active stress and accelerates shrinkage.","For pure-twist loops, radial twist windings control expansion and shrinkage, non-radial windings generate out-of-plane flows that buckle the loop, and an initially nonuniform swimmer distribution can transform a pure-twist loop into a wedge-twist loop through torsional buckling.","Varying the reversal time of swimming direction changes where bacteria accumulate, so the common picture that bacteria concentrate in splay regions is not universal; for large $t_{\\mathrm{rev}}$ they can pile up at bend-dominated walls."],"supporting_citations":[{"why":"Supplies the experimental spiral-pattern system with undulated swirling around a +1 defect that the 2D validation targets.","marker":"[22]"},{"why":"Provides measured bacterial trajectories and swirling/undulation states used to validate the spiral-pattern simulations.","marker":"[26]"},{"why":"Provides the experimental C-pattern polar-jet system and its undulation threshold and wavelength that the simulation reproduces.","marker":"[29]"},{"why":"Establishes living liquid crystals and the bend instability of extensile swimmers, the background physics for the undulation mechanisms.","marker":"[24]"},{"why":"Documents topological defects ensnaring bacteria in 2D and supplies reversal-time data used in the density-pumping section.","marker":"[27]"},{"why":"Gives the director-profile parametrization (Eq. 13) and predicts active-flow directions for 3D defect loops that the paper extends to living nematics.","marker":"[65]"},{"why":"Supplies the lattice-Boltzmann scheme for liquid-crystal hydrodynamics underlying the continuum solver.","marker":"[48]"},{"why":"Provides the specific lattice-Boltzmann implementation with finite anchoring used for the nematic flow field.","marker":"[53]"}],"fun_headline_variants":["Winding profile dictates loop dynamics in living nematic","Bacteria expand or shrink 3D defects depending on winding","Hybrid simulation reveals loop fate from wedge and twist windings","Splay accumulation not universal: bacteria seek bend","Hybrid model guides 3D living nematic loop control"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The 3D loop-defect predictions rest on hand-constructed initial director fields (Eq. 13); if real living-nematic loops are not well represented by those prescribed wedge or twist profiles, the finding that winding profile controls loop dynamics could be an artifact of the initialization rather than a property of the system.","fun_headline_variants_meta":{"raw":{"variants":["Winding profile dictates loop dynamics in living nematic","Bacteria expand or shrink 3D defects depending on winding","Hybrid simulation reveals loop fate from wedge and twist windings","Splay accumulation not universal: bacteria seek bend","Hybrid model guides 3D living nematic loop control"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000695,"raw_usage":{"total_tokens":3242,"prompt_tokens":1145,"completion_tokens":2097,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":761,"completion_tokens_details":{"reasoning_tokens":2016}},"tokens_in":761,"tokens_out":2097,"duration_ms":16162,"temperature":1.0,"reasoning_tokens":2016,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:35:22.448383+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Prepare a living nematic with a single, well-characterized 3D loop defect of known winding profile and track its lateral size over time across a range of bacterial densities; for example, a homeotropic-cell wedge-twist loop that shrinks at a density where Fig. 4(b) predicts expansion would falsify the claim. A second check is to measure bacterial density near the hybrid wall as a function of reversal time, since the claim predicts that for large $t_{\\mathrm{rev}}$ bacteria accumulate at bend-dominated regions, so observing splay-only accumulation regardless of reversal time would disprove it.","supporting_citations":[{"cited_title":"Ishimoto and E","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental spiral-pattern system with undulated swirling around a +1 defect that the 2D validation targets."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides measured bacterial trajectories and swirling/undulation states used to validate the spiral-pattern simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental C-pattern polar-jet system and its undulation threshold and wavelength that the simulation reproduces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes living liquid crystals and the bend instability of extensile swimmers, the background physics for the undulation mechanisms."},{"cited_title":"Sokolov, S","cited_arxiv_id":null,"evidence_quote":"Documents topological defects ensnaring bacteria in 2D and supplies reversal-time data used in the density-pumping section."},{"cited_title":"Ramaswamy,The mechanics and statistics of active matter, Annu","cited_arxiv_id":null,"evidence_quote":"Gives the director-profile parametrization (Eq. 13) and predicts active-flow directions for 3D defect loops that the paper extends to living nematics."},{"cited_title":"Kozhukhov, B","cited_arxiv_id":null,"evidence_quote":"Supplies the lattice-Boltzmann scheme for liquid-crystal hydrodynamics underlying the continuum solver."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the specific lattice-Boltzmann implementation with finite anchoring used for the nematic flow field."}],"review_version":1}