{"id":"629d5788-d94c-4664-8f04-2cd88b420f70","arxiv_id":"1908.05193","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"low","formal_verification":"none","parameter_count":5,"one_line_summary":"Paramagnetic sessile drops elongate linearly with the square of an applied magnetic field, an effect that scales with drop volume and magnetic susceptibility.","lead":"This paper shows that sessile drops of paramagnetic salt solutions elongate measurably in uniform magnetic fields, with the shape change growing with the square of the field strength. The authors derive a modified Young-Laplace equation from the electromagnetic stress tensor and fit it to drop outlines, pointing toward contact-free shaping of drops for microfluidics and liquid optics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (3) cannot follow from the stated EMST boundary conditions: it gives a stress proportional to -B_t^2 that vanishes at the apex, while the text and Fig.","rationale":"The paper's experimental scaling results, especially the linear decrease of D and D/D0 with B^2 and the correlation with magnetic Bond number, are directly measured and are not undermined by the theoretical concern. My concern targets claim (i): that the deformation is due to the magnetic stress difference across the interface as derived from the EMST. The printed Eq. (3) is inconsistent with a direct substitution of Eq. (1) and the stated boundary conditions; the correct small-chi magnetic stress jump contains both B_n^2 and B_t^2 contributions, whereas Eq. (3) contains only B_t^2. This is not a subtle sign convention issue because flipping the sign of the difference still leaves the normal-field term absent. The text and Fig. 2, which describe a stress largest at the apex and proportional to B_n^2, confirm that the paper itself is not using Eq. (3) literally at that point. Since the numerical fitting code is not public, the most reasonable reading is that either the implementation differs from the printed equation or the equation is wrong; in either case the central model-validation claim cannot be verified from the manuscript. The reader's weakest assumption identified the need for B_n specification, which is related but different; for chi << 1 the field-distortion part of that concern is mild. I therefore recommend UNVERDICTED pending the authors' clarification of the exact stress expression used and a corrected derivation.","tokens_in":8332,"tokens_out":19097,"duration_ms":187343,"concrete_test":"Re-derive the magnetic stress jump from Eq. (1) with the stated boundary conditions; then re-run the fitting procedure on the published 60 microliter GdCl3 drop at B = 0.6 T using (a) the printed Eq. (3) and (b) the corrected EMST expression, keeping all other fitting choices identical. Report the residuals, optimized gamma and optimized xi0*rho for both versions. If the two fits are indistinguishable after fitting xi0*rho, the shape fit cannot discriminate the stress model; if they differ, determine which expression reproduces the experimental outline and the arrows in Fig. 2.","verdict_should_be":"UNVERDICTED","load_bearing_attack":"Applying Eq. (1) to the liquid (permeability mu = mu0(1+chi)) and to air (mu0), with H_t continuous and B_n continuous across the interface, the magnetic part of the normal stress difference (inside minus outside) is -mu0*chi/2 * [(1+chi) H_n^2 + H_t^2]. In terms of the air-side flux density this is -chi/(2 mu0) * [B_n^2/(1+chi) + B_t^2]. Equation (3) as printed instead gives chi/mu0 * (B_n^2 - B^2) = -chi/mu0 * B_t^2, which vanishes exactly at the apex, where B_t = 0, and has its largest magnitude at the equator. The text and Fig. 2 state the opposite: the magnetic stress is proportional to the normal component B_n^2 and is largest at the apex. This is an internal inconsistency, not a small-chi effect: ignoring the (1+chi) factor and the missing factor 1/2 changes the functional form of the stress, removing the normal-field contribution altogether. Because the paper does not state which expression was actually solved in the numerical fits, and code and data are not public, the reader cannot tell whether the fitted outlines validate the EMST or merely an ad hoc B_t^2 term. The reader's demagnetization concern is less load-bearing here, since chi << 1 makes field distortion O(chi), but the algebraic form of Eq. (3) is a first-order issue.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a theory and experiments for the elongation of sessile drops of paramagnetic salt solutions in a uniform magnetic field. The authors derive an augmented Young–Laplace equation by adding an electromagnetic stress difference, obtained from the electromagnetic stress tensor (EMST) of Stierstadt and Liu, to the usual surface-tension and gravitational terms. They solve this equation numerically and fit the drop outline to side-view photographs. They report that the elongation is linear in B^2 and increases with drop volume and magnetic susceptibility, and they combine these data in a master curve against a magnetic Bond number. The paper suggests applications in drop actuation and liquid optics.","tokens_in":8627,"tokens_out":23795,"duration_ms":237495,"significance":"The paper's strength is the direct experimental characterization of a weakly magnetic, particle-free system: the linear B^2 scaling in Fig. 3 and the monotonic dependence on volume and susceptibility are clean, falsifiable observations, and the authors are careful to note the reversibility of the deformation and the role of contact-line friction. A correct theoretical description of this system would be of practical interest to the lab-on-a-chip community. However, the central theoretical claim is currently not supported because the magnetic-stress expression used in the shape equation is inconsistent with the EMST boundary conditions stated in the paper, and the numerical fitting does not disclose which expression was actually implemented. The manuscript therefore is not acceptable in its present form.","major_comments":[{"comment":"The derivation of the magnetic stress is internally inconsistent. Starting from Eq. (1) and imposing the stated boundary conditions (H_t continuous, B_n continuous) for a linear paramagnet (χ << 1), the magnetic part of the normal stress difference (inside minus outside) is -χ/(2μ0)[B_n^2/(1+χ)+B_t^2] (with B_n and B_t evaluated on the air side) plus the thermodynamic terms, not χ/μ0 (B_n^2 - B^2). The printed Eq. (3) reduces to -χ/μ0 B_t^2, which vanishes at the apex (B_t = 0) and has its maximum at the equator. This directly contradicts the text and Fig. 2, which state that the magnetic stress is proportional to B_n^2 and is largest at the apex. The algebra connecting Eq. (2) and Eq. (3) is correct, but Eq. (2) itself is not the EMST result for the stated boundary conditions. Because the manuscript does not state whether the numerical fits actually solved Eq. (3) as printed or another expression, and because no code or data are provided, the experimental validation of Eq. (7) cannot be assessed. Moreover, if B in Eq. (3) is taken as the undisturbed uniform applied field, the correct leading-order stress is isotropic (proportional to -χ B0^2/(2μ0) with corrections of order χ^2), so the shape change would be second order in χ; the paper's O(χ) elongation from a B_t^2 term is therefore not the correct EMST prediction for a uniform applied field. This is a load-bearing issue that the revised manuscript must address by rederiving the stress and specifying exactly which field is used in the numerical solution.","section":"§2, Eq. (2)–(3)"},{"comment":"The quantitative validation of Eq. (7) is weakened by the fitting strategy. For the zero-field outline, γ is floated; for the field-on outlines, ξ0ρ is floated; and the triple-contact-line diameter is also optimized. Since ξ0ρ is a spatially constant term, it can absorb any constant error in the pressure balance, so the successful fit of Eq. (7) does not independently determine the absolute magnitude of the magnetic stress; it mainly tests the functional form of its dependence on the surface normal. The paper reports no goodness-of-fit statistics (e.g., residual errors, R^2, or confidence intervals on the fitted shapes), and the optimized triple-line diameter changes by 0.3 mm while the measured change is 0.05 mm, which indicates a systematic inconsistency in the boundary condition at the contact line. The independent checks in Figs. 3 and 4 (linear B^2 scaling, volume and susceptibility trends) are valuable and are not vitiated by this concern, but they cannot by themselves validate the specific stress expression in Eq. (3).","section":"§3, Fitting methodology"},{"comment":"The manuscript never states how B_n (or H_n) in Eq. (3) is computed for a deformed drop. If B is the known uniform applied field, then B_n^2 + B_t^2 = B0^2 and Eq. (3) is purely tangential; if B is instead the self-consistent field including the drop's magnetization, the magnetostatic problem must be solved together with Eq. (7), but no such calculation is described. This ambiguity matters because the sign and spatial distribution of the magnetic stress depend on this choice. Please specify the field model used in the Runge-Kutta solve and, ideally, release the code or provide the numerical method in an appendix.","section":"§2–3, Field evaluation"}],"minor_comments":[{"comment":"The abbreviation 'ppw' is not defined at first use; it presumably means parts per weight, but this should be stated explicitly.","section":"Experimental section"},{"comment":"The labels V1,2,3 and χ1,2,3 are introduced only in the caption; please define them in the main text and include units on the axes.","section":"Fig. 3"},{"comment":"The text says the magnetic stress is 'directed outwards from the drop along the magnetic field lines'; since the stress in Eq. (3) is a normal traction, the wording should be revised to describe a normal stress difference rather than a force along field lines.","section":"Fig. 2 discussion"},{"comment":"Reference [31] is to the authors' previous diamagnetic study; a brief comparison of the two results would help readers assess the generality of the EMST description.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper's experimental observations may be of interest if reanalyzed with a correct stress expression, but the central equation as printed is inconsistent with the stated boundary conditions. I am not convinced that Eq. (3) is what was solved in the numerical fits; the discrepancy between the printed equation and the described stress distribution (largest at the apex) suggests that a different expression may have been implemented. I recommend that the revised manuscript be checked by a referee with specific expertise in magnetic-fluid statics, and that the authors be required to provide the code or a detailed numerical algorithm."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe one thing to know: the experimental data is probably fine, but the theory section doesn't hold together as written. The paper reports that paramagnetic sessile drops elongate linearly with B^2 and that the effect scales with volume and susceptibility. Those trends are direct measurements and look credible. The images support it.\n\nWhat's actually new: this is the first demonstration of shape change for paramagnetic salt solution drops in a uniform field, with systematic variation of volume and susceptibility. The scaling with B^2 and the magnetic Bond number collapse are useful. The authors credit Rowghanian et al. for the augmented Young-Laplace derivation and their own prior diamagnetic work, so they aren't overclaiming novelty.\n\nThe soft spot is Eq. (3). As printed, it gives the magnetic stress as χ/µ0 (B_n^2 - B^2) = -χ/µ0 B_t^2. That vanishes at the apex, where the tangential field is zero, and is largest at the equator. But the text and Fig. 2 state the opposite: stress proportional to B_n^2 and largest at the apex. You can't have both. If Eq. (3) was actually the expression solved in the numerical fits, the predicted shapes would be oblate, not prolate—the opposite of what is observed. The paper doesn't state which expression was implemented, and no code or data are available. So the model validation is not reproducible from the text.\n\nThere are also lesser issues: the fits let surface tension and the field-independent chemical potential term float against the same outlines being compared, so the good agreement is partly circular. The linear B^2 scaling and the volume/χ trends are not circular—those are measured directly. The demagnetization concern is minor because χ ~ 1e-4, but the paper should still say how the interface field is obtained.\n\nNet: the experimental scaling result is worth keeping, but the theoretical framework needs a serious correction before the paper can claim a validated model. A referee should ask for the derivation of Eq. (3) to be fixed or the equation removed, and for the fitting code and data to be posted.\n\nSend to peer review, but with the expectation that revision should be substantial.","headline":"The experimental scaling result is real and useful, but Eq. (3) contradicts the paper's own interpretation of the magnetic stress, and the model fitting is partly circular.","tokens_in":9151,"tokens_out":29103,"would_cite":false,"duration_ms":258367,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Paramagnetic sessile drops elongate linearly with the square of an applied magnetic field, and a modified Young–Laplace equation predicts their outlines.","keywords":["paramagnetic drops","electromagnetic stress tensor","Young–Laplace equation","magnetic Bond number","drop shaping","sessile drops","superhydrophobic surfaces","particle-free magnetic actuation"],"falsifier":"Map the magnetic flux density in and around a deformed paramagnetic drop—for example with a Hall-probe scan or field imaging—and compare the interface value of $B_n$ with the uniform applied field assumed in the fits. If the measured $B_n$ deviates enough to change the magnetic stress difference by more than the reported fit error, the predicted outlines and the linear $D$ versus $B^2$ relation should fail at higher fields, higher susceptibilities, or larger volumes.","tokens_in":8111,"feed_emoji":"🧲","tokens_out":15492,"duration_ms":126345,"temperature":0.7,"pith_summary":"This paper establishes that a sessile drop of paramagnetic salt solution, sitting on a superhydrophobic surface, elongates in a uniform magnetic field because of a mismatch in magnetic stress across the liquid–vapour interface. The shape change follows a modified Young–Laplace equation derived from the full electromagnetic stress tensor, and the paper confirms experimentally that the deformation grows linearly with the square of the applied field and scales with drop volume and magnetic susceptibility. The practical payoff is that magnetic fields could shape and move droplets without electrodes, contact lines, or embedded magnetic particles, complementing electrowetting and dielectrowetting for liquid lenses, liquid displays, and microfluidic mixing.","feed_headline":"Drop stretching grows with magnetic field strength squared","feed_subtitle":"A stress-balance equation predicts it, enabling particle-free magnetic shaping for liquid optics and microfluidics.","key_machinery":"The central object is the full electromagnetic stress tensor for quasi-static, non-dissipative media. Setting $D=E=0$ and using the Clausius–Mossotti approximation $\\rho(\\delta\\chi/\\delta\\rho)\\approx\\chi$ for $\\chi\\ll1$ reduces its normal stress difference across the liquid–vapour interface to a term proportional to $B_n^2 - B^2$. This term is inserted into the augmented Young–Laplace equation $0 = \\Delta\\sigma^{\\mathrm{surf}}_{nn} + \\Delta\\sigma^{\\mathrm{grav}}_{nn} + \\Delta\\sigma^{\\mathrm{EM}}_{nn}$, whose numerical solution for the axisymmetric outline $r(z)$ is fitted to drop photographs; the apex radius of curvature $b$ fixes the field-independent potential difference through $2\\gamma b^{-1} = a^l_0 - a^v_0$. The machinery converts a field measurement plus surface tension, density, volume, and susceptibility into a predicted drop outline.","core_discovery":"The paper's central claim is that the equilibrium shape of a paramagnetic sessile drop in a static uniform magnetic field is determined by the stress balance $0 = \\Delta\\sigma^{\\mathrm{surf}}_{nn} + \\Delta\\sigma^{\\mathrm{grav}}_{nn} + \\Delta\\sigma^{\\mathrm{EM}}_{nn}$, where the magnetic term, derived from the full electromagnetic stress tensor, reduces for $\\chi\\ll1$ to $\\Delta\\sigma^{\\mathrm{EM}}_{nn} = a^l_0 - a^v_0 - \\xi_0\\rho + \\chi\\mu_0^{-1}(B_n^2 - B^2)$. Numerically integrating this augmented Young–Laplace equation for the axisymmetric outline $r(z)$ and fitting it to photographed side profiles of drops made from manganese chloride and gadolinium chloride salt solutions, the authors find that the dimensionless shape parameter $D=(w-h)/(w+h)$ decreases linearly with $B^2$, that the drop elongates more when volume or magnetic susceptibility is larger, and that the normalized deformation collapses onto a single linear trend against the magnetic Bond number, $D_m/D_0 = -0.44B_m + 0.97$, with a mean fit error of $0.01$. The deformation is reversible as long as the drop volume is constant. The authors further note that the same stress-balance statement holds for ferro- and diamagnetic drops, and they propose a unified electromagnetic treatment as a route to combining electric and magnetic drop actuation.","pith_inferences":["A consequence the authors do not pursue is that the same stress-balance argument predicts a field tilted relative to the drop's symmetry axis should produce both elongation and a lateral force, so the setup could be extended from shaping to field-controlled drop transport.","The symmetry between electric and magnetic terms in the stress tensor suggests a dielectric analogue: a dielectric liquid drop in a uniform electric field should show the same linear collapse of normalized shape against an electric Bond number, and comparing the slopes would test how universal the fitted slope is.","The fits treat the triple contact line as mobile, while the paper notes its motion is inhibited by surface friction; including contact-line pinning in the model would predict a field threshold below which reversible elongation is suppressed, which is testable on rougher versus smoother superhydrophobic coatings."],"forward_implications":["For any paramagnetic sessile drop with $\\chi\\ll1$, the field-induced change in shape is proportional to $B^2$, so once surface tension, density, volume, and susceptibility are known, the outline at any field strength is fixed by the same fitted equation.","Because the same stress tensor yields the modified Young–Laplace equation for para-, dia-, and ferro-magnetic liquids, the formalism unifies magnetic and electric drop actuation, and the two can in principle be combined for shaping modes neither alone provides.","On superhydrophobic substrates, paramagnetic salt solutions give particle-free magnetic actuation: no ferrofluid particles are needed, and the deformation is reversible provided the drop volume is constant.","The linear collapse of normalized shape change against the magnetic Bond number means measurements at one volume, susceptibility, and field can be rescaled to predict another drop's deformation."],"supporting_citations":[{"why":"It supplies the full electromagnetic stress tensor and the boundary conditions from which the magnetic stress difference in the augmented Young–Laplace equation is derived.","marker":"20"},{"why":"It provides the dimensionless derivation of the modified Young–Laplace equation from the electromagnetic stress tensor for ferrofluid drops, which the paper extends to paramagnetic sessile drops.","marker":"30"},{"why":"It is the authors' prior experimental validation of the modified Young–Laplace equation for diamagnetic sessile drops, and the present work generalizes that validation to paramagnetic drops.","marker":"31"},{"why":"It defines the magnetic Bond number as the ratio of magnetic to surface energy, which the paper uses for its scaling collapse.","marker":"39"},{"why":"It gives the literature molar magnetic susceptibilities of the salts and water used to compute the drop susceptibilities.","marker":"40"},{"why":"It supplies the Axisymmetric Drop Shape Analysis methodology in electric fields on which the authors' fitting algorithm is based.","marker":"37"},{"why":"It outlines the Axisymmetric Drop Shape Analysis procedure that the fitting algorithm adapts to magnetic fields.","marker":"38"}],"fun_headline_variants":["Magnetic stretch: drop shape follows B²","Drop elongation scales with field squared","B² law predicts magnetic drop shaping","Field-driven drop reshaping: the B² effect"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction rests on treating the applied magnetic field as uniform and unaffected by the drop's own magnetization; if the drop distorts the field significantly, the magnetic stress and the linear $B^2$ scaling would change.","fun_headline_variants_meta":{"raw":{"variants":["Magnetic stretch: drop shape follows B²","Drop elongation scales with field squared","B² law predicts magnetic drop shaping","Field-driven drop reshaping: the B² effect"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000198,"raw_usage":{"total_tokens":1352,"prompt_tokens":911,"completion_tokens":441,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":527,"completion_tokens_details":{"reasoning_tokens":387}},"tokens_in":527,"tokens_out":441,"duration_ms":5363,"temperature":1.0,"reasoning_tokens":387,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:20:28.374462+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Map the magnetic flux density in and around a deformed paramagnetic drop—for example with a Hall-probe scan or field imaging—and compare the interface value of $B_n$ with the uniform applied field assumed in the fits. If the measured $B_n$ deviates enough to change the magnetic stress difference by more than the reported fit error, the predicted outlines and the linear $D$ versus $B^2$ relation should fail at higher fields, higher susceptibilities, or larger volumes.","supporting_citations":[{"cited_title":"Stierstadt \\ and\\ author M","cited_arxiv_id":null,"evidence_quote":"It supplies the full electromagnetic stress tensor and the boundary conditions from which the magnetic stress difference in the augmented Young–Laplace equation is derived."},{"cited_title":"Rowghanian , author C","cited_arxiv_id":null,"evidence_quote":"It provides the dimensionless derivation of the modified Young–Laplace equation from the electromagnetic stress tensor for ferrofluid drops, which the paper extends to paramagnetic sessile drops."},{"cited_title":"Dodoo \\ and\\ author A","cited_arxiv_id":null,"evidence_quote":"It is the authors' prior experimental validation of the modified Young–Laplace equation for diamagnetic sessile drops, and the present work generalizes that validation to paramagnetic drops."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It defines the magnetic Bond number as the ratio of magnetic to surface energy, which the paper uses for its scaling collapse."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It gives the literature molar magnetic susceptibilities of the salts and water used to compute the drop susceptibilities."},{"cited_title":"Bateni , author S","cited_arxiv_id":null,"evidence_quote":"It supplies the Axisymmetric Drop Shape Analysis methodology in electric fields on which the authors' fitting algorithm is based."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It outlines the Axisymmetric Drop Shape Analysis procedure that the fitting algorithm adapts to magnetic fields."}],"review_version":1}