{"id":"84b362c2-f722-4e0a-b077-e722b2013da4","arxiv_id":"2411.14363","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Measurements show that the hindered settling of discs and rods collapses onto the sphere curve when rescaled by particle volume relative to a sphere.","lead":"The paper measures how fast discs and rods settle in a viscous fluid when packed together, and finds they slow down less than spheres do. The authors show all their shape data collapse onto the sphere curve if shifted by a factor set by particle volume, which could let engineers predict settling of many natural particles.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The single-particle Stokes normalization may not match the suspension orientation state; if so, the fitted shift factors S(A) are biased, and the claimed A^-1/A^2 volume scaling is not yet established.","rationale":"Agree with the reader that the most fragile step is the normalization by US. The paper's Table I and the methods paragraph show an explicit asymmetry: rods are normalized by their vertical settling speed, flat particles by an average over initial orientations. In the suspension, the mean orientation can differ and can depend on phi and A, so the measured H is not necessarily the true H for the suspension's orientation state. This is not a stylistic concern; it feeds directly into the fitted S(A) because a vertical rescaling of H can be traded against a horizontal shift when matching the R-Z curve. The disc branch is especially sensitive: only three aspect ratios are used, the S range is only a factor of about 3, and the flat particles are exactly the ones for which US is least controlled. The End Matter even states that a systematic error in Us affects the scaling quality. A concrete orientation measurement and renormalization would settle the issue. Other concerns, such as missing error bars on S and the data not being available in the arXiv version, are real but secondary; they affect the strength of the evidence rather than the correctness of the normalization. The reader's CONDITIONAL verdict remains appropriate, so no change is recommended.","tokens_in":7474,"tokens_out":27651,"duration_ms":280396,"concrete_test":"Measure the orientation distribution of discs and rods during steady settling at each phi, for example by imaging with a second camera or using index-matched 3D particle tracking; compute the orientation-averaged Stokes velocity <US> appropriate to that distribution from single-particle drag coefficients; and re-normalize H = U/<US>. Then re-fit S(A) to the Richardson-Zaki n=4.5 curve. If the corrected S(A) still scales as A^-1 for discs and A^2 for rods within about 20% per point, the volume-scaling claim survives; if the slope changes by more than the spread of the three points, the claimed collapse is an artifact of the normalization.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (H collapses onto spheres with a horizontal shift S(A) equal to the sphere/shape volume ratio) is only as clean as the single-particle Stokes velocities used to normalize H. In Table I and the paragraph starting \"We determined the stokes velocity from measurements of individual particles,\" US is measured for vertically oriented rods, but for flat particles it is an average over initial orientations. In the settling suspension, both rods and discs are in a semi-dilute regime where free rotation is suppressed; the orientation distribution is set by the flow and steric constraints, and it can vary with aspect ratio and volume fraction. If the mean settling velocity U(phi) corresponds to a different orientation state than the US used, then H is multiplied by a factor f(phi,A) that is not constant. Because S(A) is obtained by minimizing the deviation between the shifted H(a*) and the R-Z sphere curve, a vertical rescaling of each shape's data can be partly absorbed by a different horizontal shift. The disc branch is the vulnerable one: only three aspect ratios (A=0.05, 0.07, 0.15) are used, and S(A) proportional to A^-1 spans only a factor of about 3; a modest systematic error in US for the largest or smallest discs could change the fitted slope. The paper acknowledges this in the End Matter (\"an error in Us is a systematic error in H for a given particle type\") but does not quantify it. Unless the orientation-state correction is shown to be negligible, the A^-1/A^2 volume scaling is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript presents sedimentation experiments on monodisperse suspensions of three flat particle types (hexagons with aspect ratio A = 0.15, 0.07 and a disc with A = 0.05) and three rod types (A = 3, 10, 19), all in the Stokes regime. The authors extract the mean settling velocity from interface tracking, with PIV used at low volume fraction, and normalize by the measured single-particle Stokes velocity to form the hindered settling function H = U/U_S. They plot H against the normalized interparticle separation a* = a/L and compare with the literature sphere Richardson–Zaki curve. They find that the shaped-particle data collapse onto the sphere curve after a horizontal shift S(A), and that S(A) is proportional to A^-1 for flat particles and A^2 for rods, which they identify with the ratio of the sphere volume to the particle volume at fixed largest dimension. From this they argue that backflow proportional to particle volume is the dominant mechanism of hindered settling for nonpolar axisymmetric particles.","tokens_in":7772,"tokens_out":17316,"duration_ms":168864,"significance":"If the proposed volume-ratio collapse is correct, the paper provides a simple, useful rule for estimating hindered settling of rods and platelets and extends the classical sphere-based empirical framework to anisotropic shapes. The experimental work has clear strengths: the particles are well characterized and monodisperse, Reynolds numbers are very small, the settling interface and PIV measurements provide complementary determinations of U, finite-size checks are reported, and the comparison target is an independently compiled sphere Richardson–Zaki curve, so the collapse is not internally circular. The principal risk is that the fitted shift factors S(A), which carry the entire volume-scaling claim, are sensitive to the normalization of H by a single-particle Stokes velocity measured in a specific orientation state and to the conversion from volume fraction to interparticle separation. These systematic uncertainties are acknowledged in the End Matter but not quantified. Because the multi-shape collapse is the central result, the scaling law is not yet established at the level claimed.","major_comments":[{"comment":"The central claim depends on the normalization H = U/U_S, but U_S is measured for vertically oriented rods and as an average over initial orientations for the flat particles. In a settling semi-dilute suspension, steric and hydrodynamic interactions constrain particle rotations, so the orientation distribution may differ from the one used in the single-particle measurement. Since the settling velocity of an anisotropic particle depends on orientation, this multiplies H by a factor f(phi,A) that is not necessarily constant. The shift factor S(A) is then obtained by minimizing the deviation of the shifted H(a*) from the sphere R-Z curve, and a vertical rescaling of a data set can be partly absorbed as a different horizontal shift. The End Matter notes that an error in U_S is a systematic error in H for a given particle type, but it does not quantify the size, sign, or phi-dependence of the effect. The authors should measure or tightly bound the orientation distribution in the settling suspension (the particles are visible in the images) and report the range of single-particle settling velocities over orientations. Without this, the fitted S(A), and hence the A^-1 and A^2 scaling, may be biased.","section":"Normalization paragraph after Fig. 2; Table I; End Matter"},{"comment":"The horizontal shift S(A) is not a directly measured quantity: it is the one-parameter horizontal translation that best fits each data set to the sphere R-Z curve. This procedure can absorb systematic errors in the conversion from volume fraction to interparticle separation, a* = k(A) phi^{-1/3}. The factor k(A) is taken from the literature, and an error in k(A) or in phi (both acknowledged as systematic in the End Matter) is largely indistinguishable from a shift in S. The paper reports no residuals, no confidence intervals for S(A), and no propagation of the stated uncertainties in U, U_S, phi, and k(A). It also includes literature data [22,29] without stating whether the same definitions of L, a*, and U_S were used. To make the collapse convincing, the authors should tabulate S(A) with uncertainties, show residual plots of the shifted data against the R-Z curve, and demonstrate that plausible variations in k(A) and phi do not change the fitted power-law exponents.","section":"Fig. 3 and End Matter"},{"comment":"The flat-particle branch of the volume-ratio law rests on only three aspect ratios (A = 0.15, 0.07, 0.05), spanning a factor of about three in A, and two of these three points are hexagons rather than discs. For a regular hexagonal prism with circumdiameter L, the geometric volume is (3*sqrt(3)/8) A L^3, whereas the disc volume used in the text is (pi/4) A L^3; at fixed A the sphere-to-particle volume ratio differs by about 20%. If the scaling is truly volume-only, the hexagon points should be placed using their actual volume, not the disc formula, and the statement that the flat shapes are 'effectively axisymmetric' should be checked against the single-particle Stokes velocity. Moreover, Fig. 3(c) shows no error bars on S(A), so the A^-1 slope for the flat branch is not yet supported with stated precision. The authors should report fitted prefactors and confidence intervals, and ideally add at least one more disc aspect ratio to separate a volume effect from a shape effect.","section":"Fig. 3(c); Table I"}],"minor_comments":[{"comment":"Table I contains typographical spacing errors in the Reynolds number entries (e.g., '6 .33 × 10−5' and '4 .96 × 10−4'); these should be corrected.","section":"Table I"},{"comment":"The power-law fits should be stated with numerical prefactors and uncertainties; the text says the prefactors are 'of order unity' without giving values, and the abstract's phrase 'exactly the ratio' is stronger than the data support.","section":"Fig. 3(c)"},{"comment":"The text says the collapse is achieved for 'all six of our particle shapes' but also includes two literature data sets; clarify that the comparison in Fig. 3(b) involves eight sets and state the aspect ratios and normalization used for Refs. [22,29].","section":"Fig. 3(b)"},{"comment":"The interpretive statement that backflow is the 'dominant contribution' and that the effect 'emerges from terms simply proportional to volume' goes beyond what an empirical horizontal-shift collapse can demonstrate; it would be better framed as a consequence of the volume-ratio scaling rather than as an independently established mechanism.","section":"Abstract and conclusion"},{"comment":"The paper should clarify the definition of 'effectively axisymmetric' for hexagons and report how the hexagon's single-particle Stokes velocity compares with a disc of the same circumdiameter and volume.","section":"Fig. 1 and Table I"}],"recommendation":"major_revision","confidential_remarks":"To the editor: This is a carefully executed experimental study with a genuine external check against sphere data. My reservation is not about the experimental technique but about the strength of the inference from fitted horizontal shifts to a universal volume-ratio law. The orientation-normalization issue can likely be addressed with the existing images, and the uncertainty analysis can be done with the existing data. I would therefore treat this as a major revision rather than a rejection, provided the authors either supply the orientation and error analysis or substantially temper the central claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nHere's my read of arXiv:2411.14363.\n\nThe genuinely new piece is the first controlled set of hindered-settling measurements for flat particles (discs and hexagons) at finite volume fraction, plus new rod data. The experimental work looks careful: monodisperse particles, low Reynolds numbers, interface tracking and PIV cross-checks, and a check of finite-size effects. The central observation—that all data collapse onto the sphere Richardson-Zaki curve with a horizontal shift S(A) tracking the particle volume relative to a sphere—is simple and would be very useful if it holds. Using the independent sphere curve as the benchmark is a real external test, and including the earlier rod data from Herzhaft & Guazzelli and Turney et al. strengthens the case.\n\nThe soft spots are about how much weight the quantitative scaling can carry. First, S(A) is extracted by fitting each shape's data to the sphere curve, and the exponents A^{-1} and A^2 are read off from only three aspect ratios per branch, with no reported uncertainties on S(A). A modest systematic error in one shape's data could change the fitted slope. Second, the normalization by the single-particle Stokes velocity is a real concern: US is measured for vertically oriented rods, but for the flat particles it is an average over initial orientations. In the settling suspension, the orientation distribution is set by flow and steric constraints. The paper itself acknowledges in the End Matter that an error in US is a systematic error in H for that particle type. If the effective US in the suspension differs from the measured one, the shift S(A) could absorb a spurious vertical rescaling. The disc branch is the vulnerable one: the A^{-1} law spans only a factor of about 3 in S(A) over the measured range.\n\nNone of this is fatal. The collapse across several independent data sets is evidence that the volume-based picture is on the right track, and the orientation issue is addressable with a check or a bound. But the quantitative claim—the exact A^{-1}/A^2 scaling—is not yet nailed down. The authors should be asked to provide error bars on S(A), make the data available, and discuss the orientation-state issue.\n\nBottom line: this deserves a serious referee. It is a solid experimental contribution with a potentially important unifying observation, but the exponents need more support before I'd treat it as a law. I would bring it to a reading group, and I'd cite it for the new flat-particle data even while cautioning on the collapse.","headline":"New flat-particle settling data and a plausible volume-based collapse, but the exponents rest on three points and a possibly biased Stokes normalization.","tokens_in":8293,"tokens_out":3373,"would_cite":true,"duration_ms":29331,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that the hindered settling of discs and rods collapses onto the sphere curve after a horizontal shift equal to the particle volume relative to a sphere, so mean settling speed is set by volume, not orientation.","keywords":["hindered settling","Stokes flow","anisotropic particles","discs","rods","sedimentation","Richardson-Zaki","aspect ratio"],"falsifier":"Measure the orientation distribution of rods and discs during settling and recompute H using a single-particle Stokes velocity matched to that distribution; if the required shift S(A) changes with the orientation state, or if a cylinder with aspect ratio one does not land at the predicted S=2/3, the volume-only collapse is an artifact of the normalization.","tokens_in":7253,"feed_emoji":"⬇️","tokens_out":7889,"duration_ms":71923,"temperature":0.7,"pith_summary":"This paper measures how suspensions of discs and rods settle in viscous fluid at low Reynolds number. It claims that all six particle shapes follow the same hindered settling function as spheres once each data set is shifted horizontally by a factor S(A), and that this factor is just the volume of a sphere of the same largest dimension divided by the particle volume. The result would matter because it turns a many-body hydrodynamic problem into a volume-ratio correction: for any nonpolar axisymmetric particle, the mean settling curve could be predicted from sphere data without accounting for orientation dynamics. The paper also argues that this supports backflow, proportional to particle volume, as the dominant hindering mechanism.","feed_headline":"Discs and rods settle like spheres after a volume shift","feed_subtitle":"Six disc and rod suspensions collapse onto the sphere settling curve using only the particle volume ratio.","key_machinery":"The argument is carried by three objects: the hindered settling function H=U/U_S, the normalized interparticle separation a* = a/L = k(A) $φ^{{-1/3}}$, and the horizontal shift factor S(A) found by total least squares collapse of each shape's H(a*) onto the sphere Richardson-Zaki curve (1-φ)^4.5. S(A) absorbs all shape dependence, and its measured $A^{{-1}}$ and $A^{2}$ scalings match the ratio of sphere volume to particle volume. Physically, the paper invokes Batchelor's dilute-limit calculation, which attributes -5.5 of the -6.55 φ slope to upward backflow proportional to particle volume; the successful volume-only collapse suggests this backflow term dominates even at finite φ.","core_discovery":"The central discovery is that, over aspect ratios from thin flat particles to long rods, the hindered settling function H(φ)=U(φ)/U_S is the same universal curve as for spheres once plotted against normalized interparticle separation and shifted by a shape-dependent factor S(A). For flat particles S(A) follows $A^{{-1}}$ and for rods $A^{2}$, with prefactor 2/3, which is exactly the ratio of the volume of a sphere of diameter L to the volume of the particle. The collapse includes previously published fiber data and holds up to the semi-dilute regime where the mean separation is comparable to the particle size. The authors conclude that the dominant contribution to hindered settling is volume-proportional backflow, and that orientational degrees of freedom do not change the mean settling velocity.","pith_inferences":["If the volume-only rule holds more generally, it provides a cheap predictive protocol: calibrate once with spheres and compute particle volume to estimate settling in any new suspension.","A natural theoretical target would be a Batchelor-style dilute-limit calculation for anisotropic particles; the data suggest the backflow coefficient is shape-independent, while only the volume normalization changes.","The rule may extend to non-axisymmetric shapes such as rectangular platelets or ellipsoidal grains, but the present data only directly support effectively axisymmetric particles.","Orientation may still matter for quantities the paper did not measure, such as the sharpness of the settling interface and velocity fluctuations, even if it drops out of the mean settling speed."],"forward_implications":["The mean hindered settling of any nonpolar axisymmetric particle, such as a clay platelet, a paper fiber, or a red blood cell, can be predicted from the sphere master curve using only the ratio of the sphere volume to the particle volume.","Existing empirical sphere fits, such as Richardson-Zaki with n≈4.5, can serve as the universal reference curve for shaped particles after the S(A) shift.","At equal normalized interparticle separation, flat particles and rods hinder settling less than spheres, with hindering beginning only when the typical separation falls below one largest dimension.","The dilute-limit backflow contribution identified by Batchelor for spheres remains the dominant hindering mechanism in the semi-dilute regime for anisotropic shapes."],"supporting_citations":[{"why":"Provides the compiled sphere sedimentation data and the Richardson-Zaki fits that define the reference master curve for the collapse.","marker":"[8]"},{"why":"Supplies the empirical (1-φ)^n form used as the sphere reference curve for the total least squares fit.","marker":"[9]"},{"why":"Gives the dilute-limit -6.55φ result and the backflow decomposition used to interpret the volume scaling.","marker":"[10]"},{"why":"Contributes the fiber suspension settling data included in the collapse to test the volume-only shift.","marker":"[22]"},{"why":"Contributes the magnetic-resonance rod settling data included in the collapse.","marker":"[29]"},{"why":"Supplies the total least squares method used to compute the optimal horizontal shift S(A).","marker":"[27, 28]"}],"fun_headline_variants":["Volume shift makes discs and rods settle like spheres","Discs and rods hindered settling matches spheres after volume scaling","Volume alone predicts hindered settling of discs and rods","Non-spherical settling: one volume scaling fitted all","Hindered settling of discs and rods collapses to sphere curve"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the single-particle Stokes velocity used to normalize each suspension is the correct reference for the particles' orientation state while settling, since rods were normalized by the vertically oriented velocity and flat particles by an average over initial orientations.","fun_headline_variants_meta":{"raw":{"variants":["Volume shift makes discs and rods settle like spheres","Discs and rods hindered settling matches spheres after volume scaling","Volume alone predicts hindered settling of discs and rods","Non-spherical settling: one volume scaling fitted all","Hindered settling of discs and rods collapses to sphere curve"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000486,"raw_usage":{"total_tokens":2342,"prompt_tokens":837,"completion_tokens":1505,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":453,"completion_tokens_details":{"reasoning_tokens":1428}},"tokens_in":453,"tokens_out":1505,"duration_ms":11070,"temperature":1.0,"reasoning_tokens":1428,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:16:35.757489+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the orientation distribution of rods and discs during settling and recompute H using a single-particle Stokes velocity matched to that distribution; if the required shift S(A) changes with the orientation state, or if a cylinder with aspect ratio one does not land at the predicted S=2/3, the volume-only collapse is an artifact of the normalization.","supporting_citations":[{"cited_title":"Brzinski III and D","cited_arxiv_id":null,"evidence_quote":"Provides the compiled sphere sedimentation data and the Richardson-Zaki fits that define the reference master curve for the collapse."},{"cited_title":"Richardson and W","cited_arxiv_id":null,"evidence_quote":"Supplies the empirical (1-φ)^n form used as the sphere reference curve for the total least squares fit."},{"cited_title":"Batchelor, Sedimentation in a dilute dispersion of spheres, Journal of fluid mechanics 52, 245 (1972)","cited_arxiv_id":null,"evidence_quote":"Gives the dilute-limit -6.55φ result and the backflow decomposition used to interpret the volume scaling."},{"cited_title":"Herzhaft and ´E","cited_arxiv_id":null,"evidence_quote":"Contributes the fiber suspension settling data included in the collapse to test the volume-only shift."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Contributes the magnetic-resonance rod settling data included in the collapse."}],"review_version":1}