{"id":"99dd9450-eb5b-49a5-9961-db565d54270c","arxiv_id":"2605.25167","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"Develops an energetic variational formulation for electrohydrodynamics of surfactant-laden droplets that derives coupled Stokes equations, interfacial conditions with Marangoni and Maxwell stresses, electrostatics, surfactant transport, and contact-line dynamics, plus reduced models and a numerical ","lead":"The paper develops an energetic variational framework for surfactant-laden droplet dynamics coupled with electric fields, deriving the governing equations from Onsager's principle. This provides a unified way to model fluid flow, surfactant transport, and electrostatic effects relevant to microfluidic technologies.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Whether variation of the chosen Rayleighian w.r.t. interface position recovers the exact Marangoni + Maxwell jump condition without omitted or sign-flipped terms.","rationale":"The reader's weakest_assumption correctly isolates the variational construction as the single point on which all listed equations depend. This is an internal-correctness issue rather than an external-consensus disagreement. The numerical illustrations in the graph-reduced model provide indirect support but cannot substitute for verifying the variation step itself.","tokens_in":1715,"tokens_out":357,"duration_ms":25377,"concrete_test":"From the derivation section, extract the expression obtained by varying the Rayleighian w.r.t. interface velocity and compare term-by-term to the standard physical condition [[-pI + 2μD(u) - ε(EE - ½|E|²I)]]·n = σκn + ∇_sσ + contact-line contributions. If the Maxwell or Marangoni term is absent or carries an incorrect coefficient, the central claim fails.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The strongest claim requires that the specific free-energy (surfactant-dependent surface tension plus electrostatic) and dissipation (bulk viscous, possibly interface/contact-line) functionals, when minimized subject to incompressibility, produce the full set of equations including the interfacial stress balance. This step is load-bearing because any algebraic error in the functional derivatives with respect to normal interface velocity would omit the Marangoni term (surface gradient of tension) or fail to produce the correct Maxwell stress jump [[T^M]].n. The paper performs this construction once; the claim stands or falls on that single calculation being free of hidden assumptions about how the electric potential couples to the interface motion.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript develops an energetic variational formulation for the electrohydrodynamics of surfactant-laden droplets based on Onsager's principle. The governing equations are obtained by minimizing the Rayleighian (rate of free-energy change plus dissipation) subject to incompressibility; this is claimed to simultaneously produce the bulk Stokes equations in each phase, the interfacial stress balance incorporating Marangoni and Maxwell stresses, the electrostatic equation, the surface transport equation for insoluble surfactant, and moving contact-line dynamics. The paper also derives a reduced mean-curvature-flow model by replacing viscous dissipation with Rayleigh dissipation, reduces the sessile-droplet problem to a 1D graph formulation, and presents a first-order IMEX scheme with illustrative numerical results.","tokens_in":1868,"tokens_out":588,"duration_ms":21898,"significance":"If the central variational derivation is free of algebraic error, the work supplies a thermodynamically consistent framework that unifies bulk flow, interface conditions, surfactant transport, and electrostatics. Such energetic formulations are valuable for constructing structure-preserving discretizations and for extending to more complex multiphysics droplet problems. The reduction steps to mean-curvature flow and the 1D graph model, together with the numerical illustration, add practical utility.","major_comments":[{"comment":"The load-bearing step is the functional derivative of the Rayleighian with respect to the normal interface velocity (the step that produces the interfacial stress balance). The manuscript must explicitly display this calculation, including the contributions from the surfactant-dependent surface tension in the free energy and from the electrostatic energy, to confirm that the Marangoni term (surface gradient of tension) and the Maxwell-stress jump [[T^M]]·n appear with correct signs and without omitted terms. Any hidden assumption about how the electric potential couples to interface motion would undermine the strongest claim.","section":"Variational derivation (Rayleighian minimization)"},{"comment":"The dissipation functional must be stated with sufficient precision (bulk viscous dissipation plus any interfacial or contact-line contributions) so that the reader can verify that the resulting bulk Stokes equations and the normal-stress jump are recovered exactly when the variation is performed subject to the incompressibility constraint.","section":"Definition of dissipation functional"}],"minor_comments":[{"comment":"The abstract and introduction should cite the specific prior energetic-variational works on two-phase flow or electrohydrodynamics that the present construction extends.","section":"Introduction"},{"comment":"In the 1D graph reduction, the notation for the reduced variables (height h, surfactant concentration, potential) should be introduced with a clear diagram or coordinate definition to avoid ambiguity when the curvature and normal vectors are written in 1D form.","section":"Graph formulation"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the careful reading and constructive comments. The suggestions will improve the clarity of the variational derivation in the manuscript. We address each major comment below.","responses":[{"response":"We agree that explicitly displaying the functional derivative calculation is essential for rigor. In the revised manuscript, we will add a dedicated subsection or appendix detailing the variation of the Rayleighian with respect to the normal interface velocity. This will include the explicit contributions from the surfactant-dependent surface tension term in the free energy, which produces the Marangoni stress as the surface gradient of tension, and from the electrostatic energy, which produces the Maxwell stress jump [[T^M]]·n with the correct signs. The electric potential is determined by solving the electrostatic equations in the bulk domains with appropriate interface conditions, and its coupling to the interface motion enters solely through the energy variation without additional assumptions.","revision_made":"yes","referee_comment":"[Variational derivation (Rayleighian minimization)] The load-bearing step is the functional derivative of the Rayleighian with respect to the normal interface velocity (the step that produces the interfacial stress balance). The manuscript must explicitly display this calculation, including the contributions from the surfactant-dependent surface tension in the free energy and from the electrostatic energy, to confirm that the Marangoni term (surface gradient of tension) and the Maxwell-stress jump [[T^M]]·n appear with correct signs and without omitted terms. Any hidden assumption about how the electric potential couples to interface motion would undermine the strongest claim."},{"response":"We will revise the manuscript to state the dissipation functional with greater precision, explicitly writing the bulk viscous dissipation as the integral over each phase of 2μ|D(u)|^2 and specifying any interfacial dissipation or contact-line dissipation terms. This precise definition will enable readers to directly verify the recovery of the Stokes equations in the bulk and the normal stress jump at the interface upon performing the variation subject to incompressibility.","revision_made":"yes","referee_comment":"[Definition of dissipation functional] The dissipation functional must be stated with sufficient precision (bulk viscous dissipation plus any interfacial or contact-line contributions) so that the reader can verify that the resulting bulk Stokes equations and the normal-stress jump are recovered exactly when the variation is performed subject to the incompressibility constraint."}],"tokens_in":1467,"tokens_out":494,"duration_ms":25380,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main takeaway is a unified energetic variational model for two-phase Stokes flow with insoluble surfactants and electrostatic effects on moving interfaces. They define a free energy combining surfactant-dependent surface tension with electrostatic energy, pair it with a viscous dissipation functional, and apply Onsager's principle by minimizing the Rayleighian under incompressibility. This produces the bulk Stokes equations, the interfacial stress balance that includes both Marangoni and Maxwell stresses, the electrostatic equation, surfactant transport on the interface, and contact-line motion.\n\nThe same framework also yields a reduced mean-curvature model by swapping the dissipation and a 1D graph formulation for sessile droplets. They close with a first-order implicit-explicit scheme and a few numerical examples showing coupled surfactant and electric-field effects.\n\nThe strength is the systematic construction: everything comes from the same energy and dissipation choices without separate ad-hoc interface conditions. The reductions show the approach is flexible rather than one-off.\n\nThe load-bearing step is the variation with respect to interface position that must recover the exact jump condition with the correct Marangoni term and Maxwell stress. The abstract states it works, but the algebra has to be free of omitted gradients or sign errors; that calculation is where any mistake would show up. The numerics are illustrative only, so they do not test the derivation.\n\nThis is for readers already working with variational methods in interfacial flows or with electrohydrodynamic droplet models. Someone in that area would find the unified treatment and the reduced models useful.\n\nIt should go to referees. The modeling is coherent and the target application is concrete.","headline":"This paper gives a variational derivation that pulls the full set of electrohydrodynamic equations for surfactant droplets out of a single Rayleighian minimization.","tokens_in":2344,"tokens_out":394,"would_cite":false,"duration_ms":30280,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Minimizing the Rayleighian derives the full set of equations for surfactant-laden electrohydrodynamic droplets.","keywords":["energetic variational formulation","electrohydrodynamics","surfactant-laden droplets","Onsager principle","Stokes flow","moving interface","Maxwell stresses","droplet dynamics"],"falsifier":"A direct check showing that the derived interfacial stress balance fails to recover the standard combination of viscous, Marangoni, and Maxwell stresses would falsify the claim.","tokens_in":2614,"feed_emoji":"💧","tokens_out":694,"duration_ms":20142,"temperature":0.7,"pith_summary":"The paper develops an energetic variational framework for two-phase Stokes flow coupled with surfactant transport on a moving interface and electrostatic effects. It applies Onsager's principle by minimizing the Rayleighian, which sums the rate of change of free energy and a dissipation functional, subject to the incompressibility constraint. This single minimization produces the bulk Stokes equations, the interfacial stress balance that includes both Marangoni and Maxwell stresses, the electrostatic equation, the surface transport equation for insoluble surfactant, and the moving contact-line dynamics. The same structure also yields a reduced mean-curvature model and a one-dimensional graph formulation for sessile droplets, together with a numerical scheme. The approach matters because it builds thermodynamic consistency directly into models used for digital microfluidics and related liquid-handling technologies.","feed_headline":"Energy minimization yields full equations for surfactant droplets in fields","feed_subtitle":"Rayleighian minimization under incompressibility produces bulk flow, interface stresses, electrostatics, surfactant transport and contact-li","key_machinery":"The Rayleighian functional, minimized subject to the incompressibility constraint following Onsager's principle.","core_discovery":"By minimizing the Rayleighian, defined as the sum of the rate of change of the free energy and the dissipation functional, subject to the incompressibility constraint, the formulation simultaneously yields the Stokes equations in each bulk phase, the interfacial stress-balance condition incorporating Marangoni and Maxwell stresses, the electrostatic equation, the surface transport equation for insoluble surfactant concentration, and the moving contact-line dynamics.","pith_inferences":["The same Rayleighian construction may extend to three-dimensional non-graph geometries while preserving the same stress and transport equations.","Energy-stable discretizations could be designed by mimicking the continuous variational structure at the discrete level.","The framework might accommodate additional physics, such as soluble surfactants, by adjusting only the free-energy functional."],"forward_implications":["Replacing the viscous dissipation functional with Rayleigh dissipation produces a reduced model in which surfactant-laden droplets evolve by motion by mean curvature.","Representing sessile droplets as graphs reduces the full system to a one-dimensional coupled model for liquid height, surfactant concentration, and electric potential.","A first-order implicit-explicit scheme can be applied to the reduced graph system to compute the coupled effects of surfactant transport and electric fields.","The variational structure automatically incorporates energy dissipation into the dynamics of the interface and the surfactant."],"fun_headline_variants":["Rayleighian minimization derives surfactant droplet equations in fields","Energy minimization derives equations for surfactant-laden droplets in fields","Variational Rayleighian yields coupled equations for droplet electrohydrodynamics","Minimizing Rayleighian gives electrohydrodynamic equations for surfactant droplets"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The system dynamics are obtained by minimizing the Rayleighian subject to the incompressibility constraint, following Onsager's principle, with specific choices for the free energy and dissipation functionals.","fun_headline_variants_meta":{"raw":{"variants":["Rayleighian minimization derives surfactant droplet equations in fields","Energy minimization derives equations for surfactant-laden droplets in fields","Variational Rayleighian yields coupled equations for droplet electrohydrodynamics","Minimizing Rayleighian gives electrohydrodynamic equations for surfactant droplets"]},"model":"grok-4.3","cost_usd":0.011965,"raw_usage":{"total_tokens":5221,"prompt_tokens":655,"num_sources_used":0,"completion_tokens":65,"cost_in_usd_ticks":119649500,"prompt_tokens_details":{"text_tokens":655,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":4501,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":655,"tokens_out":65,"duration_ms":37149,"temperature":1.0,"reasoning_tokens":4501,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-29T23:32:31.442980+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A direct check showing that the derived interfacial stress balance fails to recover the standard combination of viscous, Marangoni, and Maxwell stresses would falsify the claim.","supporting_citations":[],"review_version":1}