{"id":"629de4fc-2877-4941-9b19-94c6dc470641","arxiv_id":"2411.14636","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Surface-attached active nematic drops, simulated without the thin-film approximation, adopt a wide range of stable shapes and flows whose symmetry is set by boundary anchoring.","lead":"A model of a droplet of living or synthetic material that generates its own internal flow is simulated on a flat surface. The drop settles into many different stable shapes depending on how its internal units are anchored at the boundaries, offering a route to control soft active materials.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (4) is not the free-energy minimizer for a unit director; the simulations solve an unconstrained vector Laplace problem, so the computed state diagram may not describe the claimed active nematic drops.","rationale":"The reader's verdict is sound in direction and confidence; my disagreement is only about which assumption is the weakest. The paper's central claim — that the two-dimensional state diagram and the boundary-alignment symmetry-breaking rule are generic properties of surface-attached active nematic drops — rests on the director being computed correctly. As written, Eq. (4) is not the Euler-Lagrange equation for F in Eq. (3) when p is a unit director: the constraint yields ∇²p + |∇p|²p = 0, not ∇²p = 0, and the weak form in Eq. (S5a) confirms that the code solves the unconstrained vector problem. This is an internal-consistency issue rather than a modeling choice; if p is instead an unconstrained order parameter, then the active stress should involve a Q-tensor-like normalization, and the interpretation of winding numbers as physical angles is not justified. The paper's own caveat about flow coupling is real but secondary: the director should first be the correct minimizer under the unit constraint. The proposed check — rerunning the same grid with p=(cos θ, sin θ) and ∇²θ=0, and reporting all panels of Fig. 2 — will settle whether the state diagram is robust. Until then, the appropriate verdict remains the reader's CONDITIONAL. I give credit where due: the numerical method is standard and clearly described, and the contrast with the thin-film approximation is a legitimate motivation.","tokens_in":12967,"tokens_out":12085,"duration_ms":135665,"concrete_test":"Independently re-derive the Euler-Lagrange equation for F in Eq. (3) with |p|=1, then recompute the director field as p=(cos θ, sin θ) with ∇²θ=0 inside the drop and θ fixed by the boundary winding numbers (θ=w_i π/2 on the interface, θ=w_s π/2 on the substrate). Rerun the full Fig. 2 parameter grid at |Caα| = 2.5, 4.5, and 7.5 for all six (ws, wi) pairs, keeping Eqs. (1)-(2) and all other boundary conditions unchanged, and compare steady shapes, symmetry-breaking status, and the claimed (ws, wi) -> (ws+1, wi+1) with Caα -> −Caα equivalence against the published diagrams. If any qualitative transition or symmetry rule changes, the published state diagram is an artifact of solving the unconstrained p equation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Even taking the strong-elastic limit at face value, Eq. (4) is not the first variation of F in Eq. (3) for a nematic director, because p is a unit vector (|p|=1). Minimizing F under this constraint gives ∇²p + |∇p|²p = 0, equivalently ∇²θ=0 for p=(cos θ, sin θ); the linear vector Laplace equation ∇²p=0 is the minimizer only for an unconstrained order-parameter vector. The weak form in Eq. (S5a) confirms that the code solves the unconstrained problem, so the magnitude of p is free to vary inside the drop even though the boundary fixes |p|=1. The active stress in Eq. (2), −αpp, then depends on this uncontrolled magnitude, and the drawn 'orientation' is not the director field used in the equations. Since every transition line in Fig. 2, the 'symmetry breaking is set by boundary alignment only' rule, and the period-2 mapping (ws, wi) -> (ws+1, wi+1) with Caα -> −Caα are computed with this director field, the central claim is not yet supported by the model as stated. The caveat the authors do state — that the period-2 equivalence 'is unlikely to hold if the strong elastic limit is relaxed' — is a separate but related worry: the director is slaved and never reacts to the flow it drives. The most efficient check is to rerun the same grid with the correct unit-director equation; if the state diagram changes, the main conclusions are artifacts of the omitted constraint.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies a two-dimensional Stokes-flow model of a surface-attached drop of an active nematic fluid. The model combines incompressibility and Stokes equations with an active stress -alpha p p, kinematic and stress boundary conditions at the free surface, and a director field that is assumed to relax instantaneously to the minimum of a one-constant Frank energy. The authors simulate the full equations without a thin-film approximation, sweep over extensile/contractile activity and discrete quarter-turn anchoring conditions at the substrate and interface, and report steady-state shapes and flows, a state diagram, scaling laws for maximum velocity and dissipation, and a demonstration of reversible switching by changing the interfacial anchoring. The central claims are that stable steady states exist for all tested parameters, that symmetry breaking is controlled by the boundary-alignment mismatch, and that (ws, wi) states are equivalent under (ws, wi) -> (ws+1, wi+1) with Ca_alpha -> -Ca_alpha.","tokens_in":13253,"tokens_out":7411,"duration_ms":74275,"significance":"If the results are correct, they substantially extend the thin-film theory of active drops by showing multiple stable morphologies and quantitative scaling laws in a simple parameter-free model. The design-principle claim, especially the reversible anchoring switching, is potentially important for experiments on active droplets and colonies. The paper has genuine strengths: the model is minimal, the parameters are few, the simulations are time-dependent with interface tracking, and many results are presented as falsifiable predictions, including the state diagram, the scaling exponents, and the symmetry-breaking rule. However, the significance is conditional on correcting the director-field equation described below; as written, the computation may solve a different model from the one stated, and the absence of a convergence study and of deposited code or data weakens the quantitative claims.","major_comments":[{"comment":"The director field is treated as an unconstrained vector. Eq. (4) is presented as the first variation of F in Eq. (3), but for a nematic director with |p| = 1 the first variation under the unit-vector constraint gives p x (curl p) = 0, equivalently Laplace's equation for the angle theta, not the vector Laplace equation del^2 p = 0. The weak form in Eq. (S5a) confirms that the code solves the unconstrained vector Laplace problem, with unit magnitude imposed only on the Dirichlet boundaries. Because the active stress in Eq. (2) is -alpha p p, an uncontrolled interior magnitude of p directly changes the flow forcing, so the shapes, flows, state diagram in Fig. 2, the \"symmetry breaking is set by boundary alignment only\" rule, and the period-2 equivalence under Ca_alpha -> -Ca_alpha are computed for a different order-parameter model than the nematic director model described. Please rerun the parameter sweep with the constrained director equation, or explicitly state that p is an unconstrained order-parameter vector and justify that choice, and report whether the state diagram, the equivalence mapping, and the scaling laws survive.","section":"Model section, Eq. (4) and Supplementary Eq. (S5a)"},{"comment":"The claims that the drop adopts equilibrium shapes for all values of Ca_alpha and w, and that the state diagram in Fig. 2 uncovers all possible states, are stronger than the actual sampling. The simulations use only ws in {0,1}, wi in {0,1,2,3}, and a limited range of |Ca_alpha| (up to 4.5 in Fig. 2 and up to 10^2 in Fig. 3). No grid-convergence study, temporal-refinement check, or uncertainty estimates are reported for the phase boundaries or the scaling exponents. Since the scaling laws in Fig. 3 are a central quantitative result, please provide a resolution study, report the parameter ranges actually simulated, and state the numerical uncertainty behind each claimed \"all possible states\" statement.","section":"Fig. 3 and the section on scaling laws"},{"comment":"The paper correctly notes that the equivalence (ws, wi) -> (ws+1, wi+1) with Ca_alpha -> -Ca_alpha is unlikely to hold if the strong elastic limit is relaxed. This caveat is not limited to that mapping: the entire state diagram, the symmetry-breaking rule, and the scaling laws are computed with a director field that is slaved to the free-energy minimum and never coupled back to the flow. Please either add a test of flow-alignment coupling for a subset of parameters or explicitly restrict the \"quantitative principles\" and \"all possible states\" claims to the strong-elastic, flow-uncoupled model, including in the abstract and conclusions.","section":"Homeotropic substrate anchoring paragraph"}],"minor_comments":[{"comment":"The text refers to \"Fig. 1b\" through \"Fig. 1e\" when describing the time-dependent switching sequence; these references should be to Fig. 4b-e.","section":"Reversible shape and flow control section"},{"comment":"The quantity rho T s_dot_tot is called the total entropy production rate, but it can be negative in the contractile case (Fig. 3b), so it is a balance of viscous dissipation and active work rather than a total entropy production rate; please define it more precisely or use a less misleading name.","section":"Eq. (7) and Fig. 3b"},{"comment":"The phrase \"all possible states\" overstates the discrete anchoring set that was sampled; please rephrase as \"all states in the sampled anchoring set\" or otherwise qualify the claim.","section":"Abstract and Fig. 2 introduction"},{"comment":"The boundary term is written as an integral of sigma dot n dot phi_u and then stated to be expressed as an integral of the surface gradient; the equivalence should be shown explicitly in the supplement.","section":"Supplementary Eq. (S5c)"},{"comment":"There are small typographical issues, such as \"4th-order\" instead of \"fourth-order\" and an incompletely parenthesized thin-film equation in the section comparing with previous work; please correct these.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The main risk is the unconstrained director calculation: the state diagram and scaling laws may be artifacts of solving a vector Laplace problem for a non-unit order parameter. If constrained-director runs reproduce the qualitative state diagram, the paper would be a solid contribution to soft active matter. I would also encourage the editor to require either a convergence study or code/data deposition before acceptance, because the quantitative scaling claims currently lack reproducibility checks."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First: the paper has a load-bearing error in the director model. The authors write a one-constant Frank energy (Eq. 3) for a nematic director and claim the first variation gives ∇²p = 0 (Eq. 4). For a unit director, that is not the Euler–Lagrange equation; the constraint |p|=1 changes it to ∇²p + |∇p|² p = 0, equivalently ∇²θ=0 for p=(cosθ, sinθ). The weak form in the SI (S5a) confirms the code actually solves the unconstrained vector Laplace problem, so |p| varies inside the drop and the active stress −αpp is not the intended nematic stress. This isn't a cosmetic detail: the state diagram, symmetry-breaking rule, and the (ws,wi)→(ws+1,wi+1) mapping are all computed with this director field. The central claim that these are results for surface-attached active nematic drops is not yet supported.\n\nWhat is genuinely useful: the paper takes a well-studied thin-film problem and does full 2D Stokes-flow simulations with a deformable interface, showing stable steady-state shapes, symmetry breaking set by boundary anchoring, and reversible switching via time-dependent anchoring. That is a real step beyond Loisy et al. and the presentation is clear. The qualitative contrast with thin-film predictions is likely to survive a corrected director equation, even if the details change.\n\nOther soft spots: 'all possible states' is based on two substrate winding numbers, three interfacial winding numbers, and a limited |Caα| range. There's no convergence study, no error bars, no code or data. The strong-elastic-limit assumption is acknowledged by the authors and is a separate concern.\n\nWho is this for: people working on active wetting and biofilm/colony models will want to know the thin-film results may be incomplete. It deserves a serious referee, but the director equation has to be fixed—or the model must be explicitly redefined as an unconstrained vector order parameter with the free energy extended accordingly—before the results can be taken at face value. I'd send it to review, with a strong request for a corrected model and a convergence check. Worth a discussion in our group as a lesson in variational constraints.","headline":"Useful beyond-thin-film step, but the director equation is wrong for a nematic, so the state diagram as presented doesn't yet describe the claimed system.","tokens_in":13797,"tokens_out":8520,"would_cite":false,"duration_ms":80295,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A full 2D model of an active nematic drop on a rigid substrate reveals stable steady-state shapes and flows whose symmetry is set by boundary anchoring alone.","keywords":["active matter","active nematic drops","thin-film approximation","active Capillary number","winding number","boundary anchoring","symmetry breaking","surface-attached drop"],"falsifier":"Run a high-resolution simulation of the same model at $\\mathrm{Ca}_\\alpha=100$ with $w_s=w_i=0$: the paper's state diagram predicts a mirror-symmetric steady state, so observing spontaneous left-right symmetry breaking or persistent oscillations instead would contradict the claim that boundary mismatch is necessary for symmetry breaking.","tokens_in":12733,"feed_emoji":"🫧","tokens_out":11939,"duration_ms":94878,"temperature":0.7,"pith_summary":"Most theoretical work on surface-attached active drops has relied on the thin-film approximation, which assumes the drop is shallow and slowly varying. This paper drops that assumption and simulates the full two-dimensional equations for a nematic active drop on a rigid substrate. It finds that such drops settle into a wide variety of stable steady-state shapes and internal flows, selected by the angle (winding number) the active units make with the substrate and the liquid-air interface and by the active Capillary number. The paper claims that whether the drop breaks left-right symmetry is determined solely by a mismatch between substrate and interface anchoring, and that drops can be switched reversibly between distinct shapes by changing the interfacial anchoring. These behaviors are absent from thin-film predictions, which give a unique shape set by $|\\mathrm{Ca}_\\alpha|$ alone.","feed_headline":"Active drops on surfaces reach stable shapes, defying thin-film theory","feed_subtitle":"Boundary anchoring alone decides symmetry, and switching it reversibly reshapes the drop.","key_machinery":"The central object is the full 2D continuum model of an active nematic drop: the incompressible Stokes equations with an active stress $\\sigma_a = -\\alpha \\mathbf{p}\\mathbf{p}$, coupled to a Laplacian orientation field $\\nabla^2\\mathbf{p}=0$ (the strong elastic limit) and to a deformable interface with capillary stress. The behavior is governed by three dimensionless parameters: the active Capillary number $\\mathrm{Ca}_\\alpha = \\alpha/(\\gamma/R)$ comparing active stress to capillary pressure, and two winding numbers $w_s$ and $w_i$ that specify, in quarter turns, the angle the director makes with the solid substrate and the liquid-air interface. These winding numbers act as boundary conditions that determine whether the director field has defects at the contact points or in the bulk, and thus control whether the resulting flow is symmetric or not. The key methodological step is relaxing the thin-film (lubrication) approximation and solving the full equations numerically with a finite-element method on a deforming mesh, which is what reveals the multiplicity of stable states.","core_discovery":"The paper's central claim is that a surface-attached two-dimensional drop of a nematic active fluid—described by Stokes flow with an active stress $-\\alpha \\mathbf{p}\\mathbf{p}$ and an orientation field that relaxes instantly to its free-energy minimum, so $\\nabla^2 \\mathbf{p}=0$—reaches a stable steady-state shape and flow for every tested combination of the active Capillary number $\\mathrm{Ca}_\\alpha=\\alpha/(\\gamma/R)$ and boundary winding numbers $(w_s,w_i)$. The steady states are far richer than thin-film theory predicts: they include symmetric mushroom-like contractile shapes, flattened extensile lobes connected by a thin film, and asymmetric shapes with single vortices or spiral defects. Symmetry breaking occurs only when the orientation imposed at the substrate differs from that at the interface, and the handedness of the resulting asymmetry is set by the rotation sense of the units and the sign of activity. The paper also establishes an equivalence between planar and homeotropic substrate anchoring, $(w_s,w_i)\\to(w_s+1,w_i+1)$ with $\\mathrm{Ca}_\\alpha\\to-\\mathrm{Ca}_\\alpha$, and demonstrates that changing the interfacial winding number over time reversibly toggles the drop between its stable states.","pith_inferences":["If the strong elastic limit is relaxed to include flow-alignment coupling, the clean equivalence between anchoring configurations and the criterion that symmetry breaking requires a boundary mismatch could break down; the paper itself flags this as unlikely to hold, and it is the first assumption worth testing.","The predicted reversible switching suggests an experimental protocol: pattern or photocontrol the anchoring at the liquid-air interface of a surface-attached active nematic drop (for example, with light-sensitive surfactants) and observe whether the drop toggles between the corresponding steady states with the reported time scales.","The scaling laws found here could be used to infer the active stress magnitude in biological systems: measuring the maximum internal flow speed and drop radius in a microbial colony or cell aggregate would give $\\alpha$ from $|\\mathbf{u}|_{\\max}\\sim\\mathrm{Ca}_\\alpha$ in the low-activity regime.","Because the model is two-dimensional, the corresponding 3D surface-attached drop might show additional instabilities (e.g., azimuthal symmetry breaking) that the 2D state diagram cannot capture; extending these simulations to 3D would test whether the boundary-mismatch criterion survives."],"forward_implications":["Thin-film-based predictions for surface-attached active drops are incomplete: they miss the stable symmetric states, the symmetry-breaking criterion, and the multiplicity of shapes that the full 2D model produces.","Boundary anchoring alone—not the magnitude or sign of the active stress—determines whether a surface-attached active drop breaks symmetry, so surface chemistry can be used as a control knob.","The equivalence $(w_s,w_i)\\to(w_s+1,w_i+1)$ with $\\mathrm{Ca}_\\alpha\\to-\\mathrm{Ca}_\\alpha$ means that a contractile drop with homeotropic substrate anchoring behaves like an extensile drop with planar anchoring, halving the number of state-diagram quadrants that need to be computed or measured.","Reversible switching of the interfacial anchoring toggles the drop between stable shapes even after large deformations, suggesting a route to reconfigurable active droplets that behave as soft actuators.","For $|\\mathrm{Ca}_\\alpha|\\lesssim1$, internal flow speed and energy dissipation follow simple scaling laws ($|\\mathbf{u}|_{\\max}\\sim|\\mathrm{Ca}_\\alpha|$, $\\dot s_v\\sim\\mathrm{Ca}_\\alpha^2$, $\\dot s_a\\sim-\\mathrm{Ca}_\\alpha^2$), so the cost of driving a given shape can be estimated from activity alone."],"supporting_citations":[{"why":"Provides the continuum active-matter framework and the definition of extensile and contractile active stresses that the model uses.","marker":"[1]"},{"why":"An earlier thin-film model of active drops that the paper positions itself against as incomplete.","marker":"[32]"},{"why":"The thin-film model of surface-attached active droplets whose prediction of unconditional symmetry breaking is the paper's central contrast.","marker":"[34]"},{"why":"Companion thin-film analysis that supplies the quarter-turn winding-number boundary conditions the authors adopt.","marker":"[36]"},{"why":"Introduces the active Capillary number concept used to define $\\mathrm{Ca}_\\alpha$.","marker":"[47]"},{"why":"Shows that free active nematic drops undergo fingering instabilities, used to emphasize that surface-attached drops instead reach stable states.","marker":"[48]"}],"fun_headline_variants":["Active drops morph into rich shapes, guided by boundary anchoring","Symmetry breaking in active drops: anchoring controls shape","Reversible reshaping of active drops via boundary anchoring","Active drops show intricate steady states beyond thin-film limit","Boundary anchoring flips symmetry of active drops"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's results assume that the orientation field of the active units relaxes instantly to its lowest-energy configuration and is never dragged or rotated by the flow, so it obeys Laplace's equation rather than evolving with the fluid.","fun_headline_variants_meta":{"raw":{"variants":["Active drops morph into rich shapes, guided by boundary anchoring","Symmetry breaking in active drops: anchoring controls shape","Reversible reshaping of active drops via boundary anchoring","Active drops show intricate steady states beyond thin-film limit","Boundary anchoring flips symmetry of active drops"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000448,"raw_usage":{"total_tokens":2247,"prompt_tokens":919,"completion_tokens":1328,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":535,"completion_tokens_details":{"reasoning_tokens":1253}},"tokens_in":535,"tokens_out":1328,"duration_ms":9621,"temperature":1.0,"reasoning_tokens":1253,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:04:43.180383+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a high-resolution simulation of the same model at $\\mathrm{Ca}_\\alpha=100$ with $w_s=w_i=0$: the paper's state diagram predicts a mirror-symmetric steady state, so observing spontaneous left-right symmetry breaking or persistent oscillations instead would contradict the claim that boundary mismatch is necessary for symmetry breaking.","supporting_citations":[{"cited_title":"Ben Amar and L","cited_arxiv_id":null,"evidence_quote":"An earlier thin-film model of active drops that the paper positions itself against as incomplete."},{"cited_title":"Loisy, J","cited_arxiv_id":null,"evidence_quote":"The thin-film model of surface-attached active droplets whose prediction of unconditional symmetry breaking is the paper's central contrast."},{"cited_title":"Loisy, J","cited_arxiv_id":null,"evidence_quote":"Companion thin-film analysis that supplies the quarter-turn winding-number boundary conditions the authors adopt."},{"cited_title":"Alert, J","cited_arxiv_id":null,"evidence_quote":"Shows that free active nematic drops undergo fingering instabilities, used to emphasize that surface-attached drops instead reach stable states."}],"review_version":1}