{"id":"2e2cbea6-91dd-4db5-b045-2e0a72266c5a","arxiv_id":"1908.08202","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For a solid sphere in a concentric drop, Stokes-flow calculations give drag, torque, drop deformation dynamics, and a pulsatile-flow protocol that suppresses droplet breakup.","lead":"A theoretical study derives exact low-Reynolds-number flow fields for a solid particle suspended inside a droplet, and uses them to predict when the droplet breaks apart in external flows. It then proposes turning the flow on and off in pulses as a practical way to transport such compound particles without rupturing the confining drop.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Single-mode, first-order Ca interface equation is used up to D=1 (drop-inclusion contact), where the deformation is O(1); Eq. (18) even predicts a negative breakup time, so the breakup and pulsatile-flow claims are not controlled.","rationale":"The reader's conditional verdict is appropriate. The weakest assumption identified by the reader (single-quadrupole truncation) is a manifestation of a deeper issue: the perturbative expansion in Ca is used all the way to the breakup point D=1, where the interface displacement is O(1) relative to the film thickness and higher-order modes are no longer small corrections. My concern sharpens the reader's point by also flagging concrete internal inconsistencies in the breakup-time predictions: Eq. (18) has a sign error that contradicts the stated qualitative behavior, and Table I's formulas for D and tb appear inconsistent with Eq. (16). These issues directly affect the pulsatile-flow claim, which is the central novel recommendation. However, the paper still contains valid exact solutions for rotation and translation, and the small-deformation relaxation dynamics may be correct for Ca well below the breakup threshold. The appropriate response is to require the authors to fix the breakup-time formulas and to validate the breakup criterion against fully nonlinear simulations or a higher-order analytical treatment before the design rules are used. This is exactly the sense of the reader's conditional verdict, so no change in the verdict category is needed.","tokens_in":14941,"tokens_out":19211,"duration_ms":195769,"concrete_test":"Perform a boundary-integral simulation of the compound particle (α=1.5, λ=1) in simple shear at Ca=0.2, starting from a sphere, and measure (i) the time t_contact at which the interface first touches the solid inclusion, and (ii) the spherical-harmonic amplitudes of the interface shape when rmin/a = 1. Compare t_contact to the single-mode prediction from Eq. (14) and Table I. If t_contact differs by more than 20%, or if the fourth-order harmonic (P4) amplitude is more than 10% of the quadrupole amplitude at contact, the single-mode breakup criterion is invalid for the pulsatile-flow parameters. Repeat for the pulsatile square wave with T=5 and T=10 to check whether contact actually occurs.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central dynamical claims—exponential relaxation, breakup time tb, Cacrit, and the pulsatile-flow transport protocol—all come from Eq. (14), which is derived by linearizing the interface shape about a sphere as S = r − αa(1 + bCa x·E·x). The kinematic condition is then reduced to a single scalar ODE for b. This is standard at O(Ca), where only the quadrupole mode appears. The problem is that the breakup criterion D=1 is applied far outside the perturbative regime. Contact requires bCa ≈ 2(α−1)/α (from Table I's rmin formula). For α = 1.5 and Ca = 0.2 (the values used in Fig. 7), bCa ≈ 0.66, and for α = 2, bCa = 1. Thus the interface displacement is O(1) relative to the film thickness, not small compared to the drop radius. The expansion in Ca has no small parameter at the breakup point, so higher-order harmonics and O(Ca^2) corrections to the flow are uncontrolled. The paper's own clarification in Section II.C notes that the thin-film limit requires (α−1) ~ O(Ca) for the expressions to hold, but the same equations are used for α = 1.5 and 2 at Ca = 0.2 to predict Cacrit, tb, and the pulsatile-flow threshold. This is an internal inconsistency in the domain of validity. Additionally, the breakup-time formulas show concrete signs of error: Eq. (18) gives tb/G^-1 = −4/[15λ(α−1)] + ..., which is negative for α→1+, while the text states that the thin-film breakup time is very large. Table I also appears to list D and tb expressions that are inconsistent with Eq. (16) and with the values plotted in Fig. 6; for simple shear the tabulated D = bCa/(2−2α) is negative for α>1. If the breakup time is wrongly signed or defined, the pulsatile-flow criterion T < tb in Section II.E is not reliable. This is load-bearing because the novelty of the paper relative to earlier compound-drop work is precisely the transient breakup dynamics and the pulsatile transport recommendation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the low-Reynolds-number hydrodynamics of a concentric compound particle consisting of a rigid spherical inclusion inside a viscous drop. It reports analytical Stokes-flow solutions for rotation, translation, and general linear ambient flows, gives closed-form expressions for the viscous torque and drag, and derives an ordinary differential equation for the time evolution of the confining drop shape in the small-capillary-number limit. On the basis of this evolution equation, the paper predicts a critical capillary number and a breakup time at which the deformed interface touches the solid inclusion, and it proposes pulsatile flow as a method to transport compound particles without breakup. The manuscript also compares simple shear, uniaxial, and biaxial flows and discusses thin-film limits for all computed quantities.","tokens_in":15299,"tokens_out":7645,"duration_ms":73826,"significance":"If the results are valid, the paper provides useful closed-form reference solutions for a class of compound-particle problems that are usually treated numerically or in bispherical series. The derivation is self-contained, has no fitted parameters, and correctly recovers known compound-droplet limits in several places, which are genuine strengths. The pulsatile-flow idea is original and experimentally testable. However, the significance and the practical claims rest on extrapolating a first-order small-deformation calculation to the point of drop-inclusion contact, and that extrapolation is not controlled by the analysis as presented. The correctable nature of the identified problems makes a major revision the appropriate outcome.","major_comments":[{"comment":"The central evolution equation is obtained from the O(Ca) single-quadrupole shape ansatz S(x,t) = r - alpha a (1 + b Ca x.E.x). The paper then integrates this equation up to D = 1, the point at which the interface touches the solid inclusion. At contact, b Ca = 2(alpha - 1)/alpha from the rmin formula in Table I. For the parameters used in Fig. 7 (alpha = 1.5, Ca = 0.2) this gives b Ca = 0.67, and for alpha = 2 it gives b Ca = 1. These values are O(1), so the interface displacement is not small compared with the film thickness, and neither the restriction to a single quadrupole harmonic nor the neglect of O(Ca^2) corrections to the flow is justified. The clarification after Eq. (19) explicitly restricts the thin-film limit to (alpha - 1) ~ O(Ca), yet the same equations are used for alpha = 1.5 and alpha = 2 at Ca = 0.2 to predict Cacrit, tb, and the pulsatile-flow threshold. This is an internal inconsistency in the domain of validity and it undermines the breakup and transport claims.","section":"Section II.C, Eq. (14) and Fig. 7"},{"comment":"The thin-film breakup time is reported as tb/G^-1 = -4/[15 lambda (alpha - 1)] - 4(15 lambda - 7)/(45 lambda) + O(alpha - 1). For alpha -> 1+ and lambda > 0, this expression tends to -infinity, which contradicts the statement in the same paragraph that the breakup time is very large. The sign is also inconsistent with the logarithmic expression for tb in Table I, whose argument gives a positive, large value in the thin-film limit. This equation needs to be re-derived and corrected before the breakup-time predictions can be trusted.","section":"Section II.C, Eq. (18)"},{"comment":"The paper equates D = 1, i.e., contact between the deformed interface and the solid inclusion, with breakup of the confining drop. Contact is a geometric condition, not a dynamical rupture condition: whether the film then drains, dewets, or ruptures depends on additional physics that is not part of the model. Since the pulsatile-flow protocol is designed specifically to prevent D = 1, this identification is load-bearing and should be either justified with a separate argument or softened to a statement about onset of contact rather than breakup.","section":"Section II.C, Eq. (16) and definition of breakup"}],"minor_comments":[{"comment":"Equation (11) appears to be inconsistent with the expression for V_D in Appendix A: the numerator contains 9 alpha^5 - 5 alpha^5 - 4, while the appendix gives 9 alpha^5 - 5 alpha^3 - 4. This is presumably a typographical error, but it should be corrected and the inequality V_D/V_P < 1 should be verified with the corrected expression.","section":"Section II.B, Eq. (11)"},{"comment":"Several entries in Table I are typeset in a way that is hard to parse; for example, the simple-shear expression for D reads as b Ca/(2 - 2 alpha), which is negative for alpha > 1, and the uniaxial and biaxial D entries are ambiguous. Please reformat the table so that each fraction is unambiguous.","section":"Table I"},{"comment":"The design criterion at the end of the pulsatile-flow section states that stability is maintained when T > tau, but the preceding discussion and Fig. 7 show that large T leads to breakup and small T prevents it. The inequality should read T < tau (or the definition of tau should be clarified), and the statement should be checked against the values used in Fig. 7.","section":"Section II.E"}],"recommendation":"major_revision","confidential_remarks":"The domain-of-validity issue is the main reason for major revision. The paper's linearized evolution equation is a legitimate leading-order result, but the breakup and pulsatile-flow conclusions require either a restricted parameter regime, a higher-order calculation, or numerical verification. I would not recommend rejection because the underlying derivation is structured and the identified errors appear correctable within the scope of the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. The solid half is the exact Stokes-flow part: rotation of a concentric solid-in-drop particle, the modified torque, the drag formulas (including the connection to Rushton–Davies and Sadhal–Oguz), and the comparison of steady deformations in shear, uniaxial, and biaxial flows. Those are real, reproducible additions, and the authors are honest about which pieces reduce to earlier compound-droplet results. The transient interface equation at O(Ca) is standard single-mode perturbation theory, and within its stated small-deformation regime it is fine.\n\nThe soft spots are where the paper tries to extract breakup and a pulsatile-flow protocol from that first-order expansion. Contact between interface and inclusion means bCa is O(1) relative to the film thickness; at α=1.5, Ca=0.2 the deformation is not small, so higher harmonics and O(Ca^2) flow corrections are uncontrolled. The paper's own clarification that the thin-film limit requires (α−1)~O(Ca) does not rescue the α=1.5 and 2 plots. There are also concrete internal errors: Eq. (11) does not match the Appendix A expression for VD/VP; Table I gives simple-shear D = bCa/(2−2α), which is negative for α>1, and the breakup-time column has the same sign problem; Eq. (18) gives negative tb for α→1+, opposite to the text's claim that the breakup time is very large. These are not cosmetic. The breakup time and Cacrit are load-bearing for the pulsatile-flow recommendation, and right now they are not reliable. The pulsatile idea is plausible—switching off the flow before the interface reaches the inclusion should work if the linear relaxation is trusted—but as a published design rule it is premature without a controlled nonlinear analysis or numerical simulation.\n\nWho should read this: people working on compound droplets, core-shell particles, or microfluidic encapsulation who want closed-form drag and torque and an orienting comparison of flow types. The exact rotation and translation parts will get cited. The breakup and pulsatile part needs repair before it is used.\n\nRecommendation: send to peer review, not desk-reject. A good referee should ask for corrected Eq. (11), corrected Table I, a sign-checked re-derivation of tb, and either a matched-asymptotics treatment of the near-contact regime or numerical evidence that the single-mode equation remains predictive at D=1.","headline":"The exact rotation/translation results are solid and citable, but the breakup-time and pulsatile-flow claims rest on a first-order shape expansion pushed past its validity, plus some concrete sign/transcription errors.","tokens_in":15868,"tokens_out":3637,"would_cite":true,"duration_ms":38264,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76D07","76T99"],"pacs":["47.55.dr","47.63.mc"],"model":"deepseek-v4-flash","headline":"This paper derives exact low-Reynolds solutions for a rigid sphere inside a liquid drop, and shows that a pulsatile flow can shuttle the compound particle without rupturing its drop.","keywords":["compound particle","Stokes flow","drop deformation","concentric inclusion","breakup","pulsatile flow","capillary number","microfluidics"],"falsifier":"Measure the time-dependent shape of a compound particle in a steady shear flow at moderate capillary number. If the relaxation from a deformed state is not a single exponential with the predicted time constant $\\tau$, and if the steady shape exceeds the quadrupole prediction with noticeable higher-order modes before the critical capillary number is reached, then the central truncation fails.","tokens_in":1596,"feed_emoji":"💧","tokens_out":2105,"duration_ms":62202,"temperature":0.7,"pith_summary":"This paper asks how a solid sphere confined inside a liquid drop, a compound particle, behaves when pushed, spun, or placed in an imposed flow. It solves the low-Reynolds Stokes equations analytically for a concentric inclusion and derives closed-form expressions for the flow fields, the drag and torque on the particle, and the time-dependent shape of the confining interface. The central results are that rotation preserves the concentric configuration, translation does not unless an extra body force is applied to the drop, and the deformed interface always relaxes exponentially toward a steady shape with a relaxation time set by the viscosity ratio and size ratio. On that basis it claims that a pulsatile flow whose period is shorter than the breakup time can transport the compound particle while keeping the drop intact. The practical interest is that these formulas give material-selection and operating-time guidelines for preparing and moving encapsulated cells, organisms, or colloids in microfluidic devices.","feed_headline":"Pulsatile flow can carry a particle in a drop without breaking it","feed_subtitle":"New exact solutions show when a compound particle stays intact in shear, extension, and biaxial flow.","key_machinery":"The load-bearing object is the assumed interface shape truncated to a single quadrupole mode, $S(x,t)=r-\\alpha a\\,(1+b\\,Ca\\,x\\cdot E\\cdot x)$, where $E$ is the symmetric part of the imposed velocity-gradient tensor and $b$ is the single time-dependent amplitude. The kinematic boundary condition at the interface turns this shape ansatz into one scalar ordinary differential equation, $Ca\\,\\mathrm{d}b/\\mathrm{d}t=-b\\,g(\\alpha,\\lambda)-f(\\alpha,\\lambda)/h(\\alpha,\\lambda)$, whose coefficients are rational functions of the size ratio $\\alpha$ and viscosity ratio $\\lambda$. All subsequent statements, the exponential relaxation, the critical capillary number, the breakup time, the comparison among shear, uniaxial, and biaxial flows, and the pulsatile-flow protocol, follow algebraically once this scalar equation is accepted.","core_discovery":"On the paper's own terms, the discovery is that the entire low-Reynolds dynamics of a concentric compound particle reduces to a small set of explicit formulas. For rotation, the spherical interface is an exact solution: the particle and drop rotate at different rates, and the torque on the particle is $L_P=-8\\pi\\Omega_P\\mu a^3 \\lambda\\alpha^3/[1+\\lambda(\\alpha^3-1)]$. For translation, the concentric state is not steady because the drop lags behind the solid; the paper calculates the extra force $F_{\\mathrm{ext}}$ needed to make $V_P=V_D$ and gives closed-form drag laws for the thin-film and interface-confinement limits. In a general linear flow, the interface shape $S(x,t)=r-\\alpha a\\,(1+b\\,Ca\\,x\\cdot E\\cdot x)$ evolves by $Ca\\,\\mathrm{d}b/\\mathrm{d}t=-b\\,g(\\alpha,\\lambda)-f(\\alpha,\\lambda)/h(\\alpha,\\lambda)$, an exponential relaxation with time scale $\\tau=Ca\\,h(\\alpha,\\lambda)/g(\\alpha,\\lambda)$. Setting the deformation parameter $D=1$ marks breakup when the deformed interface touches the solid inclusion, and the resulting critical capillary number and breakup time are tabulated for simple shear, uniaxial, and biaxial flows. Finally, because the interface needs a time $\\tau$ to move, a square-wave pulsatile flow with period $T