{"id":"64ffd392-31c6-4b26-a0c6-0ea93ec09007","arxiv_id":"1908.00596","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Simulations map the bistability boundary between oblate and prolate vesicle shapes in wall-bounded shear flow and classify the transition as a saddle-node bifurcation with a diverging transition time.","lead":"Vesicles are soft fluid bags whose shape under flow depends on how hard the flow pulls. Using numerical simulations, this paper maps the flow strength at which a flattened vesicle snaps from an oblate (disc-like) shape to a prolate (cigar-like) shape near a bottom wall, and shows how this transition changes the vesicle's lift from the wall.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Near-wall Ca_c is extracted from a non-autonomous lift trajectory; the reported saddle-node boundary may be an artifact of h changing during the long transition.","rationale":"The reader's weakest assumption, that the exponent -1/2 is assumed by the fit and not independently verified, is valid and concerning. I do not dispute it. However, the more fundamental issue is that the near-wall data used to extract Ca_c are generated in a non-autonomous system: the vesicle's height, and hence the effective wall correction to Ca_c, changes during the long transition. The paper's own lift data (Section 4) show substantial vertical migration over the same dimensionless times as ttr. If h(ttr) is large near criticality, the fitted divergence is not evidence for a saddle-node of a fixed-confinement system, and the reported boundary in Fig. 5b could be a trajectory effect. The free-space comparison does not resolve this, because it removes the wall entirely and does not test the near-wall boundary. The proposed check, reporting h(ttr) and re-fitting on runs with small height change or clamping h, directly tests whether the reported near-wall Ca_c is a true bifurcation boundary or an artifact of lift. This is load-bearing because Fig. 5b and the claimed wall-induced 15% reduction in Ca_c are central results. The verdict remains CONDITIONAL rather than moving to REJECT, since the numerical method is validated and a saddle-node may well exist in free space; the condition is to supply the h-sensitive analysis or a fixed-height control.","tokens_in":10515,"tokens_out":13561,"duration_ms":144809,"concrete_test":"Report h(ttr) for every data point in Fig. 4a and re-fit the transition-time divergence using only runs in which the centroid height has changed by less than 0.1 R0 before β crosses zero. If the re-fitted Ca_c or the exponent differs from the reported values, the near-wall boundary is confounded by lift; the decisive control is to repeat the study with the centroid height held fixed at h = 1.1 R0 (e.g., via a vertical constraint force) and compare the resulting Ca_c and divergence exponent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim treats the near-wall transition-time divergence as a saddle-node bifurcation at fixed flow conditions, summarized by the normal form ẋ = r + x² with constant r (Eq. 9). This requires the control parameter Ca − Ca_c to be constant over the transition. In the simulations, the vesicle is released at h = 1.1 R0 and lifts freely while the shape evolves (Section 4). Since the wall shifts Ca_c by about 15% (v = 0.59: Ca_c ≈ 0.88 near wall vs 1.05 free-space; Fig. 5b), Ca_c is a function of h. For Ca close to the fitted Ca_c, ttr is large, so h increases substantially before β = 0, making Ca − Ca_c(h(t)) time-dependent. The measured 1/√(Ca−Ca_c) divergence and the near-wall boundary in Fig. 5b may therefore be a property of the lift trajectory rather than a genuine saddle-node of the fixed-confinement system. The paper reports no h(ttr) values and no test controlling h; the free-space fits are cleaner but do not validate the near-wall boundary.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses boundary-element simulations of a vesicle in wall-bounded shear flow to study the transition from oblate to prolate shapes. It reports a divergence of the transition time as ttr ~ 1/sqrt(Ca - Ca_c), interprets this as a saddle-node bifurcation, and constructs a phase diagram in the (v, Ca) plane separating bistable from prolate-only regions. It also studies lift dynamics, finding that near the wall an oblate vesicle lifts faster than a prolate one despite having a smaller far-field stresslet, and reports transient asymmetric and tethered shapes at high Ca.","tokens_in":10692,"tokens_out":3266,"duration_ms":33121,"significance":"If the reported phase diagram and bifurcation classification are correct, they provide a useful characterization of shape bistability for vesicles under shear, complementing earlier studies by Spann et al. and Noguchi and Gompper. The paper's strengths include the use of a validated numerical method (reference [34]), external comparison with Spann et al.'s capillary-number range, and the collapse of shape-parameter curves in Fig. 5a. The lift-velocity comparison between oblate and prolate initial shapes is a potentially interesting result. However, the central quantitative outputs—the critical capillary numbers and the saddle-node exponent—are obtained from fits with no error bars, no convergence study, and an assumed functional form, and the near-wall measurements are complicated by the fact that the vesicle lifts freely during the transition. These issues must be addressed before the phase diagram can be considered firmly established.","major_comments":[{"comment":"The critical capillary numbers Ca_c(v) are extracted by fitting the transition time to ttr = A/sqrt(Ca - Ca_c) + t_infinity, with Ca_c, A, and t_infinity as free parameters. The exponent -1/2 is assumed a priori and is then used to classify the transition as a saddle-node bifurcation. The paper reports no confidence intervals, no goodness-of-fit measures, and no test of whether the data could be described equally well by other divergences (e.g., a logarithmic divergence or an exponent -1). Because the phase boundary in Fig. 5b and the bifurcation classification both rest on this fit, the authors should provide error estimates for Ca_c and, if possible, a discriminating test of the exponent using the same simulation data.","section":"Section 3, Fig. 4a"},{"comment":"The near-wall simulations release the vesicle at h = 1.1 R0 and allow it to lift freely while the shape evolves. Because the wall shifts Ca_c by about 15% (v = 0.59: Ca_c ≈ 0.88 near wall versus 1.05 in free space), Ca_c is a function of height. For Ca close to the fitted Ca_c, the transition time is long, so h increases substantially before β crosses zero, making Ca - Ca_c(h(t)) time-dependent. The measured 1/sqrt(Ca - Ca_c) divergence and the near-wall boundary in Fig. 5b may therefore reflect the non-autonomous lift trajectory rather than a genuine saddle-node of a fixed-confinement system. The paper reports no values of h at the transition time and no test with the height constrained. Please quantify the height change during the transition and provide at least one controlled-height simulation or an alternative check to separate the lift effect from the intrinsic bifurcation.","section":"Section 4, Fig. 4a (near-wall data)"},{"comment":"The phase boundary in Fig. 5b is drawn from the fits discussed above, while the symbols indicate simulation points labeled 'bistable' or 'prolate'. The criteria for this labeling are not specified (e.g., how long a run must remain oblate to count as bistable, and how the final β sign is determined). The agreement between the fitted boundary and the simulation points is presented as good, but with no error bars on the fitted Ca_c values and no statement of run durations, this agreement is only qualitative. The manuscript should define the labeling protocol, report the number of runs per state point, and show error bars or confidence intervals on the fitted boundary.","section":"Fig. 5b, Section 3"}],"minor_comments":[{"comment":"There is a typo in the sentence 'Examples of faccurate codes for three-dimensional simulations...' — 'faccurate' should be 'accurate'.","section":"Section 2"},{"comment":"The manuscript uses nonstandard spellings such as 'bi-dimensionnal' and 'aforementionned'; these should be corrected to standard English.","section":"Abstract and Section 1"},{"comment":"The fit of ttr versus v at fixed Ca is described as also a saddle-node transition, but no critical vc value or fit parameters are reported. Please provide the fitted expression or a table of values.","section":"Fig. 4b and surrounding text"},{"comment":"The collapse in Fig. 5a uses Ca_c from the fits in Fig. 4a and β/βmax, so the collapse is partly ensured by the rescaling itself. Please clarify whether the collapse is a test of the saddle-node form or merely a restatement of the fit.","section":"Fig. 5a, Section 3"},{"comment":"The sensitivity check uses only two perturbation amplitudes (1% and 4%) for a single reduced volume (v = 0.63). Please report how many independent runs were performed for each amplitude and whether the perturbations were varied systematically in direction as well as amplitude.","section":"Section 3, sensitivity check"},{"comment":"The claim that the oblate vesicle has a roughly 50% higher lift velocity near the wall is based on single trajectories. Please state whether this difference is robust to numerical noise and, if possible, provide error bars or repeated runs.","section":"Section 4, Fig. 7a"},{"comment":"Reference [29] is incomplete: 'J. Comp. Phys. , (2011)' lacks volume and page numbers. Please update all references to complete bibliographic information.","section":"References"},{"comment":"The critical capillary number for v = 0.635 is quoted as approximately 0.52 in the Fig. 7b caption, while Fig. 5b suggests a value near 0.515. Please ensure numerical values are consistent between the text, captions, and figures.","section":"Section 3, Fig. 5b and Fig. 7b caption"}],"recommendation":"major_revision","confidential_remarks":"The core idea is interesting and the numerical method is well established, but the manuscript's main quantitative claim—the saddle-node bifurcation and the near-wall phase boundary—rests on fits without error bars and on a non-autonomous measurement protocol. The issues are fixable with additional analysis (error estimates, exponent discrimination, controlled-height tests), so I recommend major revision rather than rejection. The authors should also be encouraged to make the simulation data for Fig. 4 and Fig. 5 available or at least tabulate the fitted parameters with uncertainties."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know about arXiv:1908.00596. It does something genuinely new: it maps the oblate-prolate bistability boundary across reduced volumes v in [0.59, 0.75] for a vesicle in wall-bounded shear flow, and it identifies the transition as a saddle-node bifurcation through the divergence of the transition time. The second thing is a caution: the quantitative boundary from the near-wall simulations is not as solid as the qualitative picture, because the vesicle is allowed to lift during the long transition and the paper never checks whether that biases the fits.\n\nWhat the paper does well. The numerical method is validated in a separate publication [34], and the boundary element scheme is appropriate. The phase diagram in Fig. 5b is a useful new map, and the authors are honest about separating directly simulated points from the fitted boundary. Comparing wall-bounded and free-space results is a nice check: for v = 0.59, Ca_c drops from about 1.05 to 0.88 near the wall, a 15% shift. The near-wall lift-velocity inversion (oblate faster than prolate near the wall) is an interesting qualitative result backed by a simple lubrication-layer argument.\n\nSoft spots, in descending order. (1) The central Ca_c values come from nonlinear least-squares fits of ttr to A/sqrt(Ca - Ca_c) + t_infinity. There are no error bars, no convergence study, and no test of whether the exponent is genuinely -1/2 or just fitted as such. The data points in Fig. 4a are sparse; the fits look plausible, but that is not the same as a measured critical exponent. (2) The stress-test concern is real. The vesicle starts at h = 1.1 R0 and lifts while the shape evolves. For Ca near Ca_c, ttr is large, so h moves by a non-negligible amount before beta crosses zero. Since Ca_c itself depends on h (the free-space value is ~15% higher), the near-wall fits describe a system with a time-dependent control parameter. The paper does not report h at ttr or any test holding h fixed. The free-space fits are cleaner; the near-wall boundary in Fig. 5b may partially absorb the lift dynamics. (3) The bistability sensitivity check uses only two perturbation amplitudes at one reduced volume. Thin, but the simulation points in Fig. 5b independently label bistable versus prolate regions, so the qualitative claim does not rest on the perturbation test alone. (4) No code or data release, so the fits are not checkable from the paper.\n\nBottom line: the qualitative picture is likely correct, and the paper is a solid contribution to vesicle dynamics. The near-wall quantitative boundary needs an error analysis, a convergence check, and an explicit treatment of the h-dependence. That is exactly what peer review is for: send it out, ask for the fit data and a check of the divergence exponent, and the paper will be stronger for it.","headline":"Genuinely new (v, Ca) bistability phase diagram and saddle-node classification for a vesicle in wall-bounded shear flow, built on a validated numerical method, but the near-wall Ca_c values come from fits that may be biased by the vesicle lifting during the transition.","tokens_in":11264,"tokens_out":6143,"would_cite":true,"duration_ms":53966,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["87.16.Dg","47.15.G-"],"model":"deepseek-v4-flash","headline":"A deflated vesicle in shear flow leaves its flattened disk shape for a stretched cigar through a saddle-node bifurcation.","keywords":["vesicles","shear flow","oblate-prolate transition","saddle-node bifurcation","bistability","capillary number","reduced volume","lift dynamics"],"falsifier":"Run the same initial condition at $v=0.59$ for capillary numbers a few percent above the fitted $Ca_c\\approx0.88$ and perturb the oblate shape by a small bump; if the vesicle returns to an oblate steady state rather than switching to prolate, the assumed saddle-node picture is wrong, and a free-exponent fit of $t_{tr}$ versus $Ca-Ca_c$ would settle whether the divergence is really $-1/2$.","tokens_in":10264,"feed_emoji":"🫧","tokens_out":11369,"duration_ms":99482,"temperature":0.7,"pith_summary":"The paper establishes, by numerical simulation, that a deflated vesicle placed close to a wall can keep either an oblate (disk-like) or a prolate (cigar-like) shape under weak shear flow, and that the two shapes coexist over a bounded region of the reduced-volume/capillary-number plane. It claims that the irreversible switch from oblate to prolate is a saddle-node bifurcation: the time needed for the switch diverges at a critical capillary number as $t_{tr}\\sim 1/\\sqrt{Ca-Ca_c(v)}$. The result maps, in these coordinates, where the bistable zone lies for reduced volumes between about 0.59 and 0.75, and it shows what a near-wall vesicle does during the switch. It also reports that, close to the wall, the flattened oblate shape lifts away faster than the stretched prolate one, opposite to the far-field rule based on stresslet strength.","feed_headline":"Disk-shaped vesicle in shear flow snaps to cigar at a critical stress","feed_subtitle":"Near the critical capillary number, the switch time diverges; the paper maps where disk and cigar shapes can coexist.","key_machinery":"The central object is the shape parameter $\\beta$, a signed measure of triaxial deformation that converts a shape transition into a zero-crossing event. The mechanism carrying the argument is the normal form of a saddle-node bifurcation, the collision in which a stable and an unstable equilibrium meet and annihilate: $\\dot{x}=r+x^2$ with $r\\propto Ca-Ca_c(v)$; for $r<0$ there are two fixed points, for $r>0$ none, and integrating the rescaled form $\\dot{X}=1+X^2$ between fixed states gives the $1/\\sqrt{r}$ divergence of the transition time. Numerically, the argument is carried by a boundary-element solution of the Stokes equations that uses the wall Green's function to include the bottom wall, together with a finite-element discretization of the Helfrich bending energy under a surface incompressibility constraint.","core_discovery":"Starting each run from a relaxed oblate shape, the authors follow the shape parameter $\\beta=(L_1/R_0-1)(L_2/R_0-1)(L_3/R_0-1)$, whose sign distinguishes an oblate ($\\beta<0$) from a prolate ($\\beta>0$) vesicle. For weak flow the vesicle stays oblate; above a critical capillary number $Ca_c(v)$ it stretches, thins in the vorticity direction, and switches sign of $\\beta$. The dimensionless time $t_{tr}$ at which $\\beta=0$ diverges as $Ca\\to Ca_c^+$ according to $t_{tr}\\sim 1/\\sqrt{Ca-Ca_c}$ plus a finite reorientation time $t_\\infty$; the same $1/\\sqrt{}$ law is reported when the reduced volume approaches its critical value at fixed $Ca$. The paper reads this as the normal form $\\dot{x}=r+x^2$ of a saddle-node bifurcation, with $r\\propto Ca-Ca_c(v)$, so the oblate solution disappears discontinuously beyond $Ca_c$ while the prolate branch persists for all shear rates. A near-wall $v=0.59$ vesicle has $Ca_c\\approx0.88$ in wall-bounded flow versus about $1.05$ in unbounded flow, a roughly 15% reduction attributed to the lubrication layer underneath the vesicle.","pith_inferences":["If the saddle-node normal form is exact, the basin of attraction of the oblate state should shrink as $\\sqrt{Ca_c-Ca}$ near the transition, so the size of a perturbation needed to trigger the switch should vanish like $\\sqrt{Ca_c-Ca}$; the paper's two perturbed runs are consistent with that but do not measure the scaling.","A microfluidic application could exploit the boundary $Ca_c(v)$: choosing a shear rate slightly below $Ca_c$ would make identical-looking vesicles switch or not switch depending on small shape perturbations, effectively amplifying initial-condition differences.","The reported 15% wall-induced drop in $Ca_c$ suggests confinement can be used as a control parameter; varying the wall distance should smoothly shift the bistable boundary, a prediction that could be tested in the same numerical setup.","If the exponent $-1/2$ is confirmed experimentally, transition-time measurements could serve as a non-invasive way to infer the bending modulus $\\kappa$ from the measured location of $Ca_c$."],"forward_implications":["The bistable region in the $(v,Ca)$ plane is bounded from above by $Ca_c(v)$; within it, history matters and a vesicle can settle into either shape depending on the initial condition.","For reduced volumes above about $0.74$, even weak hydrodynamic stress triggers the oblate-to-prolate transition, while for $v\\simeq0.59$ the oblate shape survives up to $Ca_c\\simeq0.9$.","The oblate branch of solutions ceases to exist discontinuously past $Ca_c$, so the transition is an abrupt switch rather than a continuous drift of the shape.","Near a wall at the same $Ca$ and $v$, the oblate vesicle has a lift velocity about 50% higher than the prolate one close to the wall, even though its stresslet is smaller; far from the wall the ordering reverses.","At high capillary numbers, wall-induced asymmetry can produce transient tethered shapes independently of whether the vesicle started oblate or prolate."],"supporting_citations":[{"why":"Establishes the equilibrium bending-energy shapes of vesicles, providing the oblate and prolate branches whose coexistence the flow modifies.","marker":"[23]"},{"why":"Gives the earlier free-space phase diagram for biconcave-to-prolate transitions that the new critical capillary numbers are compared against after rescaling.","marker":"[25]"},{"why":"Supplies the Green's function for Stokes flow near a no-slip wall, the basis for the boundary-element implementation of confinement.","marker":"[28]"},{"why":"Describes and validates the numerical method (boundary element plus loop-subdivision finite elements) that generates the simulation data.","marker":"[34]"},{"why":"Anchors the bistable region by identifying the reduced-volume threshold below which oblate shapes are local minima of bending energy.","marker":"[35]"},{"why":"Provides the near-wall lift-velocity results and far-field stresslet scaling that the oblate-versus-prolate lift comparison extends and reverses close to the wall.","marker":"[13]"}],"fun_headline_variants":["Critical shear flips vesicle shape, transition time diverges","Vesicle's oblate-prolate switch: critical point, divergent time","Disk to cigar at critical shear: vesicle shape transition","Near critical shear, vesicle shape change time diverges","Vesicle shape transition in shear: a saddle-node bifurcation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the transition time diverges as $(Ca-Ca_c)^{-1/2}$, because every reported critical capillary number is obtained by fitting the simulated transition times to that assumed law, and the paper gives no confidence interval for the fitted values.","fun_headline_variants_meta":{"raw":{"variants":["Critical shear flips vesicle shape, transition time diverges","Vesicle's oblate-prolate switch: critical point, divergent time","Disk to cigar at critical shear: vesicle shape transition","Near critical shear, vesicle shape change time diverges","Vesicle shape transition in shear: a saddle-node bifurcation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000319,"raw_usage":{"total_tokens":1817,"prompt_tokens":980,"completion_tokens":837,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":596,"completion_tokens_details":{"reasoning_tokens":749}},"tokens_in":596,"tokens_out":837,"duration_ms":8584,"temperature":1.0,"reasoning_tokens":749,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:44:02.092804+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same initial condition at $v=0.59$ for capillary numbers a few percent above the fitted $Ca_c\\approx0.88$ and perturb the oblate shape by a small bump; if the vesicle returns to an oblate steady state rather than switching to prolate, the assumed saddle-node picture is wrong, and a free-exponent fit of $t_{tr}$ versus $Ca-Ca_c$ would settle whether the divergence is really $-1/2$.","supporting_citations":[{"cited_title":"Kraus, U","cited_arxiv_id":null,"evidence_quote":"Establishes the equilibrium bending-energy shapes of vesicles, providing the oblate and prolate branches whose coexistence the flow modifies."},{"cited_title":"Noguchi, G","cited_arxiv_id":null,"evidence_quote":"Supplies the Green's function for Stokes flow near a no-slip wall, the basis for the boundary-element implementation of confinement."},{"cited_title":"Gounley, G","cited_arxiv_id":null,"evidence_quote":"Describes and validates the numerical method (boundary element plus loop-subdivision finite elements) that generates the simulation data."},{"cited_title":"Boedec, M","cited_arxiv_id":null,"evidence_quote":"Anchors the bistable region by identifying the reduced-volume threshold below which oblate shapes are local minima of bending energy."},{"cited_title":"farutin, C","cited_arxiv_id":null,"evidence_quote":"Provides the near-wall lift-velocity results and far-field stresslet scaling that the oblate-versus-prolate lift comparison extends and reverses close to the wall."}],"review_version":1}