{"id":"352eff60-07a6-4383-878d-e75df588e691","arxiv_id":"2501.04918","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A linear stability analysis of a passive droplet in a confined active nematic yields explicit growth rates showing stable and unstable translation and deformation modes in both extensile and contractile systems.","lead":"This paper computes how a round drop of ordinary fluid behaves when it sits inside a ring of active fluid whose microscopic swimmers push on it, and finds that small bumps on the drop surface can grow or shrink depending on how the swimmers push. The math is a step toward explaining why passive droplets inside living cells, such as nucleoli, move and deform on their own.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The dispersion relation is internally consistent in the stated Pe→0 limit; the load-bearing risk is the strong-decoupling assumption itself, which is acknowledged but not validated at finite Pe.","rationale":"The reader identified the Pe→0 decoupling as the weakest assumption, and my stress-test concurs. This is the most load-bearing premise: all subsequent results depend on the director being slaved to the instantaneous geometry through Eq. (15). The derivation appears internally consistent, so the correctness risk lies in the applicability of the limit rather than in a detectable mathematical error. The paper is transparent about the assumption and its possible consequences, which supports a conditional acceptance pending a finite-Pe check or explicit statement of scope. No new concern that would require rejecting the paper or changing the verdict was found.","tokens_in":18313,"tokens_out":18493,"duration_ms":183782,"concrete_test":"Perform a finite-Pe linear stability analysis in the same annular geometry by solving the coupled active-nematic/Stokes eigenvalue problem with a spectral method at Pe = 0.01, 0.1, 1, and 10, extracting the leading growth rates for m = 1 and m = 2, and comparing them to Λ_m from Eq. (54). If Λ_m(Pe) does not approach the Pe→0 prediction as Pe→0, or if finite Pe introduces new unstable modes, then the central claim is confined to an asymptotically small-Pe regime that is not verified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the dispersion relation (53)-(54), whose coefficients Ξc_m, Ξa_m, Ξe_m are derived under the assumption Pe→0 (strong nematic relaxation), introduced at Eq. (14) and used at Eq. (15). In this limit the director is enslaved to the instantaneous geometry via ∇²β=0, and all flow-alignment terms in the Ericksen–Leslie equation and stress become subdominant (they are O(Pe) relative to the elastic-relaxation terms, given ξ=O(1)). The derivation is long but self-consistent: the base state, linearized stresses, reciprocal-theorem reduction, and final growth rates all follow the declared conventions, and I found no algebraic or sign inconsistency. However, the Pe→0 limit is structurally load-bearing: if Pe is not small, the director can be advected and distorted by the very flows generated by active and elastic stresses, and the decoupling used throughout no longer holds. The authors acknowledge this in Sec. V, noting that additional instabilities such as spontaneous swirling flows could arise. Because the paper provides no finite-Pe check, numerical validation, or estimate of the Pe range over which the predictions remain valid, the stability conclusions are conditional on a limit that may not be realized in the biological and experimental systems cited as motivation. This is a legitimate scope limitation, not an internal inconsistency.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a two-dimensional linear stability analysis of a passive viscous Newtonian droplet surrounded by an active nematic liquid crystal in circular confinement, and of the inverted configuration in which an active nematic droplet is surrounded by a passive viscous layer. In the sharply aligned limit with strong anchoring and Pe→0, the director field is slaved to the instantaneous geometry through Laplace's equation. Using the Lorentz reciprocal theorem for Stokes flow, the authors reduce the linearized interfacial problem to a first-order amplitude equation for azimuthal Fourier modes and obtain closed-form growth rates Λ_m = (1/Ca)Ξc_m + SΞa_m + ξ²Ξe_m as functions of the radius ratio η, viscosity ratio λ, mode number m, and elastic parameter ξ. They identify stable and unstable regions for the translational mode (m=1) and deformation modes (m≥2) for both extensile (S=+1) and contractile (S=-1) systems, and they visualize the bulk force densities and surface tractions to interpret the mechanisms.","tokens_in":18562,"tokens_out":14621,"duration_ms":144689,"significance":"If correct, this is a useful analytical reference calculation for active-passive interfaces in confinement. Its strengths are the self-contained derivation of the auxiliary Stokes problem, the absence of any fitted parameters, and the explicit closed-form growth rates covering both geometries and both signs of activity. The reciprocal-theorem formulation is elegant and could be reused in related problems. The significance is partly conditional: because the entire calculation is performed in the Pe→0 (strong elastic relaxation) limit and with an isotropic viscous response, the stability phase diagrams apply only in that regime, and the paper does not demonstrate how robust the predictions are at finite Pe.","major_comments":[{"comment":"The central dispersion relation (53)-(57) is derived in the strict limit Pe→0, in which the director is slaved to the instantaneous shape. The authors acknowledge in Sec. V that finite-Pe effects can introduce spontaneous-flow instabilities, but they do not provide any estimate of Pe for the biological systems cited in the abstract and introduction, nor any small-Pe correction. Because the active and elastic stresses generate flows that can advect the director when Pe is not small, the stability phase diagrams (Figs. 2, 5, 7) are conditional on this limit. I request that the authors quantify the regime of validity (parameter estimates for bacterial suspensions and chromatin systems) and either provide a finite-Pe analysis or explicitly restrict the biological claims to Pe≪1.","section":"Sec. II.D (Eqs. 14-15) and Sec. V"},{"comment":"The closed-form growth rates are long and are not validated against any known limit or direct numerical simulation. In particular, setting S=0 and ξ=0 should recover the capillary stability of a confined passive droplet, and the m=1 mode should be neutrally stable for a passive system; the paper does not verify these limits. A few numerical spot checks of the dispersion relation against a direct solution of the linearized Stokes problem would substantially increase confidence in the algebra and in the phase diagrams.","section":"Sec. III.F (Eqs. 53-57) and Sec. IV.C"}],"minor_comments":[{"comment":"The caption refers to 'translational modes' for m≥2 and labels the capillary contribution as Ξa_m; both should be corrected to 'deformation modes' and Ξc_m, respectively.","section":"Fig. 4 caption"},{"comment":"The text states 'we assume that ξ is of order unity', but several figures (e.g., Figs. 5 and 7) use ξ²=0.005 and 0.1; please reconcile this statement with the parameter range shown.","section":"Sec. II.D"},{"comment":"The pressure inside the drop is denoted p0 without the overbar convention established for inner variables; use \\bar{p}_0 for consistency.","section":"Eq. (23)"},{"comment":"There are two typos: 'in unclear' should be 'is unclear' and 'an useful tool' should be 'a useful tool'.","section":"Sec. V"},{"comment":"In the text preceding Eq. (63), the notation Ξa_m is used for both the active contribution and, in the figure caption, for the capillary contribution; please clarify the notation consistently.","section":"Sec. III.F.2 and Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is a clean analytical contribution with no signs of circularity or fitted parameters. I found the derivation self-consistent, but the structural Pe→0 assumption and the absence of numerical validation are significant enough that the stability claims should be more carefully qualified or supplemented before publication. The paper fits the journal's scope and does not need rejection, but it needs revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, it delivers a genuinely new analytic result: the full linear stability problem for a passive drop in a sharply aligned active nematic annulus, and the inverse geometry, with growth rates decomposed explicitly into capillary, active, and elastic contributions. The Lorentz reciprocal theorem route is clean, the appendix gives the auxiliary Stokes solutions in detail, and there are no fitted parameters. I checked the main steps and found no algebraic or sign inconsistencies. The phase diagrams (Figs. 2, 5, 7, 8) are concrete enough to guide future simulations or experiments.\n\nSecond, the load-bearing assumption is Pe→0, introduced at Eq. (14) and used to reduce the director equation to Laplace's equation. In that limit the director is enslaved to the instantaneous geometry, and all flow-alignment terms drop out at leading order. The authors acknowledge in Sec. V that finite Pe could introduce new instabilities such as swirling flows, but they give no estimate of the Pe range where their predictions survive, and there is no numerical or experimental validation. That is a real limitation, though it is stated. What is less explicitly stated is the fate of the flow-alignment parameter ϖ: it appears in the original equations and in the full Ericksen–Leslie stress, but never reappears in the reduced stress or the dispersion relations. The paper should say plainly that ϖ is set to zero, or explain why it is subdominant, because a reader cannot infer that from the scaling as written.\n\nSeveral smaller points are worth noting. The elastic pressure diverges at η→0, making small drops always unstable to translation; that may be an artifact of the sharp-alignment singular geometry rather than a physical prediction. The m=1 and m≥2 modes are treated separately, and the paper notes that the fastest-growing mode is what matters, but the coupling between modes is not discussed—reasonable for a linear stability analysis. The viscosity dependence is presented as color maps, but the asymptotic limits in Eq. (60) are useful and could have been used more to build intuition.\n\nWho is this for? Active-matter theorists and soft-matter fluid dynamicists who want an exact linear benchmark for co-annular active–passive systems. It is also relevant to the nuclear-organization community, though the connection there is motivational rather than quantitative. The paper deserves a serious referee. It is honest, self-contained, and likely correct within its declared limit. I would accept it for peer review, with the main referee asks being: clarify the ϖ reduction and add a section on the expected validity of Pe→0, either by scaling arguments or a simple finite-Pe perturbation. I would not cite it in my own work within the next year unless I were actively working on this exact geometry.","headline":"A careful and transparent analytic linear stability theory for a confined active–passive interface, internally consistent but resting on an untested Pe→0 limit and a silently dropped flow-alignment parameter.","tokens_in":19072,"tokens_out":2594,"would_cite":false,"duration_ms":29270,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76A15","76D07","35B35"],"pacs":[],"model":"deepseek-v4-flash","headline":"A circular droplet inside a confined active nematic can be destabilized by activity in both extensile and contractile systems, with each Fourier mode's growth rate fixed by a compact dispersion relation.","keywords":["active nematic","passive droplet","linear stability","Lorentz reciprocal theorem","Stokes flow","extensile vs contractile","interfacial instability","confinement"],"falsifier":"Run the same two-annulus system in a finite-$\\mathrm{Pe}$ simulation or a microtubule/motor experiment with an extensile nematic, $\\eta\\approx0.2$, small $\\xi^2$, and large $\\mathrm{Ca}$: the paper predicts the $m=2$ deformation mode should grow with positive $\\Lambda_2$. If the droplet remains circular or the fastest-growing mode differs, the dispersion relation is wrong.","tokens_in":54,"feed_emoji":"🫧","tokens_out":8383,"duration_ms":121481,"temperature":0.7,"pith_summary":"The paper asks what happens to a circular passive liquid droplet sitting inside a confined active nematic, a fluid made of self-propelled elongated units, when the outer fluid is active and the whole system is squeezed into a circular domain. Using a linear stability analysis in the sharply aligned, strong-relaxation limit, it derives a dispersion relation for the growth of each interfacial Fourier mode. The central result is that activity destabilizes the droplet in both extensile and contractile systems, but in opposite ways: the translational mode ($m=1$) is stabilized by extensile activity and destabilized by contractile activity, while shape-deforming modes ($m\\ge2$) are destabilized by extensile activity for small to intermediate droplets. Capillary and elastic stresses typically resist deformation, and the theory also covers the inverse geometry of an active nematic droplet surrounded by a passive viscous layer, where extensile activity is always destabilizing. These predictions matter because the same stress balance is thought to govern passive compartments in cells, from droplets in bacterial suspensions to nucleoli and heterochromatin inside the nucleus.","feed_headline":"Passive droplets can destabilize inside active liquid crystals","feed_subtitle":"A linear model predicts when a confined droplet starts drifting or deforming, with hints for cell nuclei.","key_machinery":"The load-bearing object is the auxiliary Stokes flow used with the Lorentz reciprocal theorem. Instead of solving the actual perturbed velocity field, the paper defines a passive Newtonian problem in the same annulus with a prescribed normal traction jump $\\Sigma_0\\cos m\\theta$ at the interface, whose streamfunction is $G_m(r)$ and solves the biharmonic equation with no-slip at $r=1$ and matching at $r=\\eta$. The reciprocal theorem re-expresses the unknown real flow in terms of contour integrals around the circular interface, which involve the already-known linearized active and elastic stresses and $G_m(r)$; the arbitrary $\\Sigma_0$ cancels, leaving the ordinary differential equation and dispersion relation for $\\tilde\\zeta_1$. The same auxiliary problem serves both geometries, with the active stresses moved from the outer annulus to the droplet interior in the inverse case. Also load-bearing is the limit $\\mathrm{Pe}\\to0$, which makes the director angle $\\beta$ satisfy Laplace's equation and turns the nematic orientation into a pure boundary-value problem driven by the instantaneous droplet shape.","core_discovery":"Under the strong relaxation limit $\\mathrm{Pe}\\to0$, the nematic director is determined purely by the geometry: it is the solution of Laplace's equation with strong tangential anchoring at the droplet surface and at the outer boundary. The paper linearizes about the circular concentric base state and represents the interface as $\\zeta(\\theta,t)=\\epsilon\\tilde\\zeta_1(t)\\cos m\\theta$. Its central result is the first-order evolution equation $d\\tilde\\zeta_1/dt=\\Lambda_m\\tilde\\zeta_1$ with $\\Lambda_m=\\frac{1}{\\mathrm{Ca}}\\Xi_m^c+S\\Xi_m^a+\\xi^2\\Xi_m^e$, where $\\mathrm{Ca}$ is the active capillary number, $S=+1$ for extensile and $S=-1$ for contractile activity, and $\\xi$ is the ratio of the active length scale to the confinement radius. The functions $\\Xi_m^c$, $\\Xi_m^a$, and $\\Xi_m^e$ are explicit integrals of the streamfunction of an auxiliary Newtonian flow and depend on the radius ratio $\\eta$ and viscosity ratio $\\lambda$. The paper computes their signs across the parameter space and concludes that the translational mode is stable under extensile activity but unstable under contractile activity, whereas deformation modes are destabilized by extensile activity in small to intermediate drops; in the inverse geometry an extensile active droplet is always destabilized, and elasticity is predominantly stabilizing. The instability picture is therefore set by a competition of capillary, elastic, and active tractions with bulk active forces.","pith_inferences":["Extending beyond the paper: at finite $\\mathrm{Pe}$ the director is no longer enslaved to the geometry, so the same dispersion relation would acquire corrections from flow alignment and director advection; recomputing $\\Lambda_m$ at small finite $\\mathrm{Pe}$ would show whether the $m=1$ stability reversal survives.","Extending beyond the paper: the strong-anchoring assumption fixes the director orientation at the interfaces; allowing weak anchoring would introduce a surface energy term that could soften the elastic traction and modify the small-drop limit where the elastic growth rate diverges.","Extending beyond the paper: the linear analysis treats all modes independently, but in the inverted problem a $+1$ topological defect sits at the droplet center; defect motion at nonlinear order could break circular symmetry before the predicted linear instability takes over.","Extending beyond the paper: the result that the active contribution to translation depends only weakly on viscosity ratio suggests an experimentally robust test: the sign of $\\Lambda_1$ should be nearly independent of drop viscosity, which could be checked directly in a microfluidic device."],"forward_implications":["A passive droplet in an extensile nematic annulus should not drift on its own at linear order, but should begin to deform through $m\\ge2$ modes once the radius ratio falls in the predicted unstable window.","In a contractile environment the same droplet should instead be set into translation first, since the active contribution to the translational mode changes sign with $S=-1$.","The stability thresholds shift with the viscosity ratio and the relative elastic stress $\\xi^2$, so confinement size and drop rheology could be used to tune whether a droplet moves or changes shape.","Because all modes grow independently in the linear theory, the mode with the largest $\\Lambda_m$ is expected to dominate the observed dynamics in an experiment.","The inverse configuration suggests that an extensile active droplet inside a passive layer will generically deform or translate, which is relevant to nucleoli and other active inclusions in the cell nucleus."],"supporting_citations":[{"why":"Supplies the Lorentz reciprocal theorem that converts the unknown Stokes flow into contour integrals around the interface.","marker":"[30]"},{"why":"Provides the continuum active-nematic equations for the director and stress used throughout the model.","marker":"[31–37]"},{"why":"Supplies the adopted form of the elastic stress tensor and its derivation.","marker":"[40]"},{"why":"Provides the expression for the elastic surface traction at the anchored interface.","marker":"[41]"}],"fun_headline_variants":["Droplet instabilities in active liquid crystals","Linear stability analysis of drops in active nematics","Active stresses can destabilize confined droplets","Extensile and contractile activity deform droplets","Predicting droplet drift and deformation in active fluids"],"cache_read_input_tokens":21248,"weakest_assumption_plain":"The load-bearing premise is that the nematic director relaxes infinitely fast ($\\mathrm{Pe}\\to0$), so the flow never advects or bends the director; the director field at every instant is just the shape-determined solution of Laplace's equation.","fun_headline_variants_meta":{"raw":{"variants":["Droplet instabilities in active liquid crystals","Linear stability analysis of drops in active nematics","Active stresses can destabilize confined droplets","Extensile and contractile activity deform droplets","Predicting droplet drift and deformation in active fluids"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000263,"raw_usage":{"total_tokens":1654,"prompt_tokens":1055,"completion_tokens":599,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":671,"completion_tokens_details":{"reasoning_tokens":530}},"tokens_in":671,"tokens_out":599,"duration_ms":6394,"temperature":1.0,"reasoning_tokens":530,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:22:07.907966+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same two-annulus system in a finite-$\\mathrm{Pe}$ simulation or a microtubule/motor experiment with an extensile nematic, $\\eta\\approx0.2$, small $\\xi^2$, and large $\\mathrm{Ca}$: the paper predicts the $m=2$ deformation mode should grow with positive $\\Lambda_2$. If the droplet remains circular or the fastest-growing mode differs, the dispersion relation is wrong.","supporting_citations":[{"cited_title":"A drop of active matter,","cited_arxiv_id":null,"evidence_quote":"Supplies the Lorentz reciprocal theorem that converts the unknown Stokes flow into contour integrals around the interface."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the adopted form of the elastic stress tensor and its derivation."},{"cited_title":"Free object in a confined active contractile nematic fluid: Fixed-point and limit-cycle behaviors,","cited_arxiv_id":null,"evidence_quote":"Provides the expression for the elastic surface traction at the anchored interface."}],"review_version":1}