{"id":"4dd9623a-296b-46b5-9629-2f289a312288","arxiv_id":"1908.09147","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Diffusiophoresis modifies the colloid Batchelor scale, and an effective Péclet number built on it collapses chaotic-flow mixing times onto the no-salt curve.","lead":"This paper shows that when colloids and salt are mixed together in a flow, the colloid mixing time can be collapsed onto the same curve as without salt by using an effective Péclet number that absorbs the diffusiophoretic drift. This gives a practical way to predict how salt gradients accelerate or delay colloidal mixing in microfluidic and chaotic flows.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The effective-Péclet derivation assumes a time-independent diffusiophoretic drift Vdp~Ddp/ℓ_s, but the salt field that drives it decays on the salt mixing time, which is much shorter than the colloid mixing time in these runs; the paper never checks that the steady-drift approximation holds.","rationale":"Both the reader and I locate the weakest point in the order-of-magnitude treatment of §3.3–3.4; my specific concern is the unstated assumption that the salt-gradient magnitude (and hence Vdp) is quasi-steady over the colloid mixing time. Because the salt is initialized as a single Fourier mode and has no source, |∇ ln S| decays on the salt mixing time ≈5–9, while the diffusiophoretic colloid mixing times extend to ≈30–40; treating Vdp as a constant Ddp/ℓ_s is thus not self-evidently justified. This is a genuine gap in the derivation of Eqs. (3.16) and (3.18). However, the central claim is the empirical collapse shown in Fig. 2(right) and the Dc-independence plateau in Fig. 3, which are direct numerical results that do not depend on the derivation being rigorous. The collapse over five orders of magnitude in Pe_eff and the plateau in the salt-repelling case are non-trivial and support the effective-Péclet description even if the quoted derivation is heuristic. I therefore do not think this gap warrants changing the reader's ACCEPT verdict; it does suggest that a follow-up check of the drift time-dependence is worthwhile. My proposed test would settle whether the constant-drift approximation is the reason the collapse works or whether the collapse is more robust.","tokens_in":11764,"tokens_out":24459,"duration_ms":235963,"concrete_test":"In the existing simulations, compute the Lagrangian time average of |∇ ln S| over the colloid mixing interval 0 ≤ t ≤ Tmix,c for every run, and form ⟨Vdp⟩/Ddp = ⟨|∇ ln S|⟩. Compare this to the assumed 1/ℓ_s = √Pe_s/L used in (3.8). If the ratio ⟨|∇ ln S|⟩ / (√Pe_s/L) is significantly below 1 for runs with large Pe_eff, or if the ratio correlates with the residuals of the collapsed data in Fig. 2(right), the constant-drift assumption fails. Otherwise the steady-drift approximation is adequate and the derivation is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.2 estimates the diffusiophoretic drift as Vdp ≃ Ddp/ℓ_s and, in the quasi-static balances of §3.3 and §3.4, treats this as a constant when equating terms (a) and (c) (Ddp>0) or (b) and (c) (Ddp<0). This is the step that produces equations (3.16) and (3.18). However, the salt field S is released as an initial 1+sin x disturbance and decays by the same chaotic mixing; there is no source to sustain the salt gradients. In the numerical runs Pe_s ∈ [600,2500], so the salt mixing time Tmix,s = 3.2 ln(Pe_s/120) ≈ 5–9, while the diffusiophoretic colloid mixing times reach ≈ 30–40 (Pe_eff up to 10^8). Thus Vdp(t) ≈ Ddp |∇ ln S| decays substantially before the colloids reach their modified Batchelor scale, so ℓ_c,diff is not set by a quasi-steady drift. If Vdp is not approximately constant, the effective Péclet numbers (3.16) and (3.18) are not the correct scales for the whole mixing history. The numerical collapse in Fig. 2(right) is the only evidence, and it could be dominated by an early-time compression followed by diffusion at a different scale. This missing check on the time dependence of Vdp is the most load-bearing gap between the derivation and the central claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the joint mixing of colloids and salt in a linear strain flow and in a globally chaotic flow, in the presence of diffusiophoresis. The authors show that for a Gaussian patch in pure strain, the mixing time is rescaled by an effective Péclet number Pe_eff = Pe_c(1 + Ddp/Ds) built from the modified Batchelor scale. For chaotic advection, they derive effective Péclet numbers by balancing terms in the equation for concentration gradients: Pe_eff ~ Pe_c Ddp^2/(Dc Ds) for salt-attracting (Ddp > 0) and Pe_eff ~ Pe_c Dc Ds/Ddp^2 for salt-repelling (Ddp < 0), valid when Ddp^2/(Dc Ds) >> 1. Numerical simulations of a sine flow with random phase show that the colloid mixing time collapses onto the no-salt curve when plotted against these effective Péclet numbers over five orders of magnitude.","tokens_in":12116,"tokens_out":11989,"duration_ms":115923,"significance":"If the proposed rescaling is correct, it provides a simple predictive tool for colloidal mixing times in microfluidic and environmental flows, and it extends earlier work by Deseigne et al. (2014) to the salt-attracting configuration. The pure-strain part is supported by the analytical solution of Raynal et al. (2018) and shows an impressive collapse over eight decades. The chaotic-advection part is supported by a numerical parameter sweep with Ddp, Dc, and Ds varied independently, and the collapse over five orders of magnitude is a strong empirical indication that the effective Péclet number is the right control parameter. The gradient-equation derivation is elegant and connects the result to the classical Batchelor-scale theory.","major_comments":[{"comment":"In §3.2–§3.4, the derivation of the effective Péclet numbers assumes that the diffusiophoretic drift Vdp ∼ Ddp/ℓs is present in the quasi-static regime when the colloids reach their modified Batchelor scale. This requires the salt gradients to survive at least as long as the colloid compression. In the simulations of §3.5, Pe_s ∈ [600, 2500], so the salt mixing time estimated from Eq. (3.19) is τ_s ≈ 3.2 ln(Pe_s/120) ≈ 5–9, while the reported colloid mixing times for large Pe_eff reach ≈ 30–40. The salt field that drives Vdp therefore decays substantially before the colloids attain the scales described by Pe_eff, and Eqs. (3.16) and (3.18) are not derived for the parameter regime simulated. The paper does not test the time dependence of Vdp (for instance, by measuring |∇ ln S| or the actual drift velocity during the runs). Without this check, the collapse in Fig. 2(right) is consistent with, but does not uniquely support, the quasi-static mechanism; an early-time compression followed by ordinary diffusion of the resulting fine structures could produce a similar rescaling. The authors should either demonstrate that Vdp remains of order Ddp/ℓs throughout the relevant phase or restate Eqs. (3.16) and (3.18) as empirical scalings.","section":"§3.2–§3.4"},{"comment":"In §3.3, the condition Ddp^2/(Dc Ds) >> 1 is stated as the validity condition for the attracting case, but this is not sufficient. The quasi-static balance (a)∼(c) in Eq. (3.7) can only hold if the salt is still present when the colloids approach their modified Batchelor scale, i.e., if τ_s ≳ τ_c,eff, which in terms of Péclet numbers requires Pe_s ≳ Pe_eff. In the simulations, Pe_eff^attract can exceed Pe_s by orders of magnitude (Pe_eff up to 10^8 versus Pe_s ≈ 600–2500), so the derivation's premise is violated exactly for the points that most strongly test the collapse. Please state the full validity condition, verify it against the numerical parameters, or modify the theoretical claim accordingly.","section":"§3.3"}],"minor_comments":[{"comment":"Equations (3.16) and (3.18) write 'Pe' without a subscript; the symbol should be Pe_c to distinguish the colloid Péclet number from Pe_s and from the effective Péclet number.","section":"Eqs. (3.16), (3.18)"},{"comment":"Page 822 contains the typo 'height decades'; it should read 'eight decades'.","section":"p. 822"},{"comment":"The reference to 'Frish (1995)' should be 'Frisch (1995)' as in the standard turbulence textbook.","section":"References"},{"comment":"The condition Ddp^2/DsDc should be written with parentheses, e.g., Ddp^2/(Dc Ds), to avoid ambiguity.","section":"Abstract, §4"},{"comment":"The numerical method section states that the same code as in Volk et al. (2014) is used, but it does not report the grid resolution or any convergence test; a sentence on the resolution and validation would strengthen reproducibility.","section":"§3.5"}],"recommendation":"major_revision","confidential_remarks":"The timescale mismatch between salt and colloid mixing is the key technical issue; I would ask the authors to address it directly. The paper is within the journal's scope and the numerical evidence is good, so the work is likely acceptable after a careful revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The main thing to know: this is a solid subfield paper with one genuinely new result—the salt-attracting effective Péclet number in chaotic advection—and a convincing numerical collapse of mixing times onto the no-salt curve. The salt-repelling expression was already in Deseigne et al. (2014), so the novelty is narrower than the abstract implies. But the unified treatment of both regimes on one curve is useful.\n\nThe derivation follows Raynal & Gence's gradient equation, and the order-of-magnitude estimates leading to (3.16) and (3.18) are internally consistent. The numerics are wide-ranging (Pe_c from 7.8e3 to 3e5, Ddp^2/DcDs up to large values) and the collapse over five decades in Pe_eff is impressive. The self-citation is present but appropriate: the analytical solution from Raynal et al. (2018) and the code from Volk et al. (2014) are central to the study, and the paper cites the earlier salt-repelling result from Deseigne et al. properly.\n\nThe soft spot is the quasi-steady assumption. In the chaotic flow, the salt field is an initial 1+sin x disturbance and decays on a time scale Tmix,s ≈ 3.2 ln(Pe_s/120) ≈ 5–9 for their parameters, while the colloid mixing times reach 30–40. The drift Vdp ~ Ddp/l_s is therefore not constant over the colloid mixing process; it decays substantially before the colloids reach the modified Batchelor scale. The derivation of (3.16) and (3.18) treats the drift as quasi-steady, and the paper never checks the time-dependent balance. The collapse may still be real—the initial compression phase could set the scale—but the current argument does not explain why a decaying drift should produce a single effective Péclet number. This is the main thing I would want a referee to investigate.\n\nMinor points: the baseline curve (3.19) is empirical, so the collapse is only as good as that fit; and only one chaotic flow (random-phase sine flow) is tested, so the generality of the scaling is unproven. Neither is fatal.\n\nWho is this for: people doing microfluidic mixing with salt gradients, and to a lesser extent the scalar-mixing theory crowd. If the quasi-steady issue is resolved, it is a citable result. I would send it to peer review; a good referee can test the timing argument with the existing data or a simple modification.","headline":"The salt-attracting effective Péclet scaling is new and the numerical collapse is impressive, but the quasi-steady drift assumption deserves scrutiny because the salt gradients decay long before the colloids finish mixing.","tokens_in":12657,"tokens_out":6826,"would_cite":true,"duration_ms":72199,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76F25","76R50"],"pacs":[],"model":"deepseek-v4-flash","headline":"Effective Péclet numbers collapse diffusiophoretic mixing times onto the salt-free law.","keywords":["diffusiophoresis","chaotic advection","Batchelor scale","effective Péclet number","colloid mixing time","salt-attracting configuration","salt-repelling configuration","scalar gradient equation"],"falsifier":"Run the same random-phase sine-flow simulations and evaluate each of the six terms in the gradient equation (3.7) at the moment the colloid patch first reaches its Batchelor scale, for $D_{\\mathrm{dp}}>0$ and $D_{\\mathrm{dp}}<0$ with $D_{\\mathrm{dp}}^2/(D_c D_s)$ just above 1; if terms (d), (e), or (f) are comparable to term (c), the dominant-balance assignment behind the effective Péclet numbers is wrong and the collapse should break.","tokens_in":11571,"feed_emoji":"🧂","tokens_out":16274,"duration_ms":132082,"temperature":0.7,"pith_summary":"This paper attempts to establish that, in a laminar flow with global chaos, diffusiophoresis—the drift of colloids along salt gradients—changes colloid mixing only through the Batchelor scale, not through the functional form of the mixing-time law. The authors define effective Péclet numbers from the diffusiophoretically modified Batchelor scale, one for salt-attracting colloids ($D_{\\mathrm{dp}} > 0$, delayed mixing) and one for salt-repelling colloids ($D_{\\mathrm{dp}} < 0$, accelerated mixing), and show numerically that the mixing time then follows the same logarithmic curve as in the absence of salt. If correct, the practical prediction of mixing in microfluidic and environmental suspensions reduces to computing a modified small scale and reusing existing mixing correlations.","feed_headline":"Effective Péclet numbers collapse diffusiophoretic mixing","feed_subtitle":"A modified Batchelor scale makes delayed or accelerated colloid mixing obey the salt-free law.","key_machinery":"The key object is the effective Péclet number $Pe_{\\mathrm{eff}}$, defined by replacing the Batchelor scale in the standard relation $\\ell_{c}/L \\sim 1/\\sqrt{Pe_c}$ with the diffusiophoretically modified scale, $\\ell_{c,\\mathrm{diff}}/L \\sim 1/\\sqrt{Pe_{\\mathrm{eff}}}$. In the chaotic case the derivation runs through the equation for the gradient $\\mathbf{G} = \\nabla C$ of the colloid concentration, whose terms describe diffusion, stretching, and the drift $\\mathbf{v}_{\\mathrm{dp}} = D_{\\mathrm{dp}} \\nabla \\ln S$. In the quasi-static regime, where production and dissipation balance, the dominant balance determines the scaling: for $D_{\\mathrm{dp}}>0$ the drift term competes with diffusion, while for $D_{\\mathrm{dp}}<0$ it competes with stretching. The order-of-magnitude estimate $V_{\\mathrm{dp}} \\sim D_{\\mathrm{dp}}/\\ell_s$ is what converts these balances into the two effective Péclet numbers.","core_discovery":"In a globally chaotic flow, the paper's central claim is that the colloid mixing time is the same function of an effective Péclet number as in the absence of diffusiophoresis, $T_{\\mathrm{mix}} = f(Pe_{\\mathrm{eff}})$, with $Pe_{\\mathrm{eff}}$ built from the modified Batchelor scale. For salt-attracting colloids the effective number is $Pe_{\\mathrm{eff}}^{\\mathrm{attract}} = Pe_c \\, D_{\\mathrm{dp}}^2/(D_c D_s)$, which is much larger than $Pe_c$ and so delays mixing; for salt-repelling colloids it is $Pe_{\\mathrm{eff}}^{\\mathrm{repell}} = Pe_c \\, D_c D_s/D_{\\mathrm{dp}}^2$, much smaller than $Pe_c$ and so accelerates mixing. Both scalings follow from balancing the diffusive dissipation term against the diffusiophoretic gradient-production term (salt-attracting) or the stretching production term against the same drift term (salt-repelling), under the condition $D_{\\mathrm{dp}}^2/(D_c D_s) \\gg 1$. In the numerical sine flow with random phases the calculated mixing times collapse onto the no-salt curve $T_{\\mathrm{mix}} = 3.2 \\log(Pe_{\\mathrm{eff}}/120)$ over five decades of effective Péclet number. For the pure-strain stagnation point the same idea works analytically, with $Pe_{\\mathrm{eff}} = Pe_c(1 + D_{\\mathrm{dp}}/D_s)$ and an exact Gaussian solution.","pith_inferences":["Because the gradient-balance argument depends on the flow only through a characteristic stretching rate, the same two effective Péclet numbers should collapse mixing-time data in other globally chaotic laminar flows, such as deterministic blinking-vortex flows, provided the stretching rate is taken from the relevant Lyapunov exponent.","Near the crossover $D_{\\mathrm{dp}}^2 \\approx D_c D_s$ neither balance is dominant; an interpolation formula combining the three leading terms of the gradient equation is a natural testable extension.","The authors' closing suggestion can be made concrete: for patches much larger than the flow scale, homogenization with imposed mean salt and colloid gradients should produce effective diffusivities whose dependence on $D_c$ matches the reciprocal factors $D_c D_s/D_{\\mathrm{dp}}^2$ and $D_{\\mathrm{dp}}^2/(D_c D_s)$ seen here."],"forward_implications":["With $D_{\\mathrm{dp}}>0$, the effective Péclet number is $Pe_c \\, D_{\\mathrm{dp}}^2/(D_c D_s)$, so salt-attracting colloids are predicted to mix much more slowly than the same colloids without salt.","With $D_{\\mathrm{dp}}<0$, the effective Péclet number is $Pe_c \\, D_c D_s/D_{\\mathrm{dp}}^2$, so salt-repelling colloids mix faster, and the theory is equivalent to an effective diffusivity $D_{\\mathrm{eff}} \\sim D_{\\mathrm{dp}}^2/D_s$ independent of $D_c$.","In a globally chaotic flow the mixing time is the same logarithmic function of $Pe_{\\mathrm{eff}}$ as without diffusiophoresis, so existing no-salt correlations or fits can be reused once $Pe_{\\mathrm{eff}}$ is known.","Large scales in the dilating direction remain essentially unaffected by diffusiophoresis; the action of the salt gradient is confined to the Batchelor scale, which is what the effective Péclet number encodes.","The paper's regime of validity is $D_{\\mathrm{dp}}^2/(D_c D_s) \\gg 1$, and numerical data with $D_{\\mathrm{dp}}^2/(D_c D_s) \\gtrsim 10$ already show the predicted plateau in the salt-repelling case."],"supporting_citations":[{"why":"Defines the diffusiophoretic drift that appears in the colloid transport equation, a velocity proportional to the gradient of the logarithm of salt concentration.","marker":"Anderson (1989)"},{"why":"Supplies the scalar-gradient equation and the quasi-static production–dissipation balance used to derive the Batchelor scale and connect it to mixing time in chaotic flows.","marker":"Raynal & Gence (1997)"},{"why":"Provides the exact Gaussian-patch solution in a pure strain flow and the modified Batchelor scale used to define the stagnation-point effective Péclet number.","marker":"Raynal et al. (2018)"},{"why":"Proposed the salt-repelling effective Péclet number through the Ranz stretch model and reported experiments with the square of the diffusiophoretic coefficient about 30 times the product of the colloid and salt diffusivities, which the present theory recovers.","marker":"Deseigne et al. (2014)"},{"why":"Supplies the random-phase sine-flow numerical configuration and the baseline chaotic-mixing data for effective compressible flows.","marker":"Volk et al. (2014)"},{"why":"Provides the stretch model of mixing that motivated the earlier effective-Péclet rescaling and is the comparison point for the gradient-balance approach.","marker":"Ranz (1979)"},{"why":"Establishes the globally chaotic random-phase sine flow whose logarithmic mixing-time scaling defines the reference curve for the collapse.","marker":"Pierrehumbert (1994, 2000)"}],"fun_headline_variants":["Colloid mixing obeys one law with effective Péclet number","Diffusiophoresis rescales Péclet number for mixing","Salt-driven drift shifts mixing time via modified Péclet","Effective Péclet predicts colloid mixing under salt","Diffusiophoretic mixing collapses to universal curve"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the order-of-magnitude balance in the gradient equation: the diffusiophoretic drift is estimated as $V_{\\mathrm{dp}} \\sim D_{\\mathrm{dp}}/\\ell_s$, and among the six gradient terms only the diffusion term and one drift or stretching term survive in the quasi-static regime, with terms (d), (e), and (f) negligible; if this ordering fails, equations (3.16) and (3.18) do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Colloid mixing obeys one law with effective Péclet number","Diffusiophoresis rescales Péclet number for mixing","Salt-driven drift shifts mixing time via modified Péclet","Effective Péclet predicts colloid mixing under salt","Diffusiophoretic mixing collapses to universal curve"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000945,"raw_usage":{"total_tokens":4157,"prompt_tokens":1185,"completion_tokens":2972,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":801,"completion_tokens_details":{"reasoning_tokens":2891}},"tokens_in":801,"tokens_out":2972,"duration_ms":20165,"temperature":1.0,"reasoning_tokens":2891,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:20:19.913449+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same random-phase sine-flow simulations and evaluate each of the six terms in the gradient equation (3.7) at the moment the colloid patch first reaches its Batchelor scale, for $D_{\\mathrm{dp}}>0$ and $D_{\\mathrm{dp}}<0$ with $D_{\\mathrm{dp}}^2/(D_c D_s)$ just above 1; if terms (d), (e), or (f) are comparable to term (c), the dominant-balance assignment behind the effective Péclet numbers is wrong and the collapse should break.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the diffusiophoretic drift that appears in the colloid transport equation, a velocity proportional to the gradient of the logarithm of salt concentration."},{"cited_title":"& Gence, J.-N","cited_arxiv_id":null,"evidence_quote":"Supplies the scalar-gradient equation and the quasi-static production–dissipation balance used to derive the Batchelor scale and connect it to mixing time in chaotic flows."},{"cited_title":"Journal of Fluid Mechanics 847 , 228–243","cited_arxiv_id":null,"evidence_quote":"Provides the exact Gaussian-patch solution in a pure strain flow and the modified Batchelor scale used to define the stagnation-point effective Péclet number."},{"cited_title":"pinch of salt","cited_arxiv_id":null,"evidence_quote":"Proposed the salt-repelling effective Péclet number through the Ranz stretch model and reported experiments with the square of the diffusiophoretic coefficient about 30 times the product of the colloid and salt diffusivities, which the present theory recovers."},{"cited_title":", Mauger, C","cited_arxiv_id":null,"evidence_quote":"Supplies the random-phase sine-flow numerical configuration and the baseline chaotic-mixing data for effective compressible flows."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the stretch model of mixing that motivated the earlier effective-Péclet rescaling and is the comparison point for the gradient-balance approach."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the globally chaotic random-phase sine flow whose logarithmic mixing-time scaling defines the reference curve for the collapse."}],"review_version":1}