{"id":"d5b9a101-32c1-484b-b100-6a1c9e6b3c86","arxiv_id":"2411.13241","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Plaque size and the systole-diastole phase of particle release measurably change the molecular communication channel impulse response in a simulated stenosed carotid artery.","lead":"This paper simulates and analyzes how a plaque narrowing a carotid artery changes the arrival times of nanoparticles released into the blood, treating the artery as a molecular communication channel. It derives analytical models for non-Newtonian blood flow and shows that plaque size and the timing of release within the heart cycle alter the received signal, a step toward future nano-device-based atherosclerosis detection.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Laminar-flow assumption likely breaks down at the 75% stenosis under peak systole, so the reported speedups and CIR detection signatures for the most severe case may not represent realistic carotid hemodynamics.","rationale":"I read the paper in good faith. The analytical CIR derivations (Eqs. 8 and 14) are internally consistent under the stated assumptions, and the simulation work is a reasonable extension of the authors' prior study. However, the most consequential vulnerability is the laminar-flow premise, because it underpins every quantitative result presented as a plaque-detection signature. The Reynolds-number estimate at the stenosis throat indicates the premise is likely violated exactly in the regime (large plaque, peak systole) that produces the paper's headline speedups. This is not a fabrication or an internal contradiction—the authors disclose the assumption—but it means the results do not yet establish the claimed sensing capability for realistic carotid hemodynamics. The proposed LES/RANS comparison would directly test whether the laminar assumption changes the conclusions. I therefore agree with the reader's weakest-assumption identification and see no reason to alter the CONDITIONAL verdict; the paper should be accepted only with the understanding that the quantitative claims are contingent on laminar validity and need a turbulence-aware follow-up.","tokens_in":64,"tokens_out":8826,"duration_ms":197330,"concrete_test":"Run the same 75%-stenosis geometry and pulsatile inlet waveform with a scale-resolving simulation (e.g., LES or a validated low-Re RANS transition model) for release at peak systole. Compare the time at which the first 1% of particles reach the RX and the full CIR against the laminar result. If the first-arrival time or CIR shape changes by more than ~20% (or if the simulation develops non-laminar flow structures), the laminar assumption is invalid for this regime and the quantitative detection claims in Section IV-B require revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that plaque size and release phase within the cardiac cycle significantly alter the MC channel impulse response, yielding detectable speedups (5.5%, 33%, 43% in Section IV-B) and CIR shape changes—rests on the assumption that blood flow remains laminar in all simulated scenarios, including a 75% stenosis under pulsatile inflow. This assumption is stated explicitly in Section II: 'We assume laminar flow is present in all simulation scenarios...' and turbulence modeling is deferred to future work (Section V). For the most severe case, rp=0.75×rc, the throat radius is 0.75 mm (diameter 1.5 mm). With a mean inlet velocity of 34.2 cm/s, the throat velocity is ~5.5 m/s; at peak systole (Fig. 2 shows a factor of roughly 2–3 above the mean) it can reach ~11–16 m/s. Using blood density 1050 kg/m³ and dynamic viscosity ~4×10⁻³ Pa·s, the throat Reynolds number is Re ≈ 5000–8000, far above the ~2000 threshold for pipe-flow transition. Even without fully developed turbulence, a severe stenosis produces flow separation, recirculation, and unsteady vortices that a laminar solver cannot capture. These features would directly affect particle transport: recirculation zones can trap particles, turbulent diffusion enhances radial mixing, and the flattened turbulent velocity profile changes the convective delay distribution. Thus the simulated first-arrival times, the speedups quoted above, and the claim that PS release gives fastest arrivals are all potentially artifacts of the laminar model. Because the abstract and conclusion generalize to 'crucial indicators for early non-invasive plaque detection,' this is a load-bearing gap between the model and the intended physical scenario.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper models a 50 mm carotid-artery segment as a molecular communication channel with a circular transmitter, a passive cross-sectional receiver, and an idealized annular plaque. It derives closed-form channel impulse responses (CIRs) for Newtonian, power-law, and Herschel-Bulkley fluids in the flow-dominated regime, adds a Venturi-based travel-time reduction model, and compares the analytical velocity profiles with OpenFOAM simulations. It then imposes a human carotid pulsatile inlet velocity profile and simulates particle release at peak systole, early diastole, and late diastole for four plaque sizes. The main reported results are that the plaque accelerates traversal by 5.5%, 33%, and 43% for rp = 0.25, 0.5, and 0.75 times the vessel radius, that peak-systole release gives the earliest arrivals, and that these CIR changes could serve as a non-invasive plaque-detection metric.","tokens_in":10707,"tokens_out":9270,"duration_ms":95086,"significance":"If the results hold, the paper offers a physically plausible and novel sensing signature for molecular-communication-based plaque detection: plaque-induced changes in flow profile and CIR, including Venturi-style speedups and sensitivity to the cardiac phase at release. The analytical CIR derivations in Section III-A are transparent and reduce correctly to the Newtonian limit, the Venturi model gives a simple falsifiable prediction, and the authors explicitly compare analytical and simulated profiles rather than hiding the mismatch. The release of the simulation dataset [7] is a genuine reproducibility strength. However, the most severe stenosis case is very likely outside the laminar regime on which all simulations rest, and the simulation results are presented without uncertainty quantification, so the quantitative claims currently outrun the evidence.","major_comments":[{"comment":"The load-bearing assumption that flow remains laminar in all scenarios is almost certainly violated for the most severe stenosis. With rc = 3.0 mm and rp = 0.75×rc, the lumen radius at the throat is 0.75 mm; conservation of mass gives an average throat velocity of about 5.5 m/s at the stated mean inlet speed of 34.2 cm/s, and the pulsatile peak in Fig. 2 is several times higher. With blood density 1050 kg/m³ and viscosity near 4×10⁻³ Pa·s, the throat Reynolds number is of order 10³–10⁴, far above the transitional threshold for pipe flow. A laminar solver cannot capture flow separation, recirculation, or turbulent mixing in this regime, all of which directly alter particle residence times and CIR shape. Since the 43% speedup and the peak-systole fastest-arrival claim for rp = 0.75×rc rest on these simulations, the central quantitative claim for severe stenosis is not supported. The authors should either restrict the conclusions to stenosis levels where the flow remains laminar, or add a transitional/turbulence-resolving validation for the worst case.","section":"§II, §IV-B, Fig. 6"},{"comment":"The speedups of 5.5%, 33%, and 43% are quoted to two significant figures, yet Fig. 6 shows single curves with no error bars, confidence intervals, or repetition count, and the text does not define the \"traversal speed\" metric (first-arrival time, median arrival, or shift in the full CIR). Without this information, a reader cannot distinguish a genuine plaque-induced channel change from sampling noise or from the choice of metric. Please state whether the runs are deterministic, report the number of repetitions and the dispersion across them, and define exactly how the traversal speed was computed.","section":"§IV-B, Fig. 6"},{"comment":"The piecewise radius function in Eq. (15) is inconsistent with the definition of rp in §II. rp is defined as the radial expansion of the plaque, so the lumen radius in Region 3 should be rc − rp, not rp; this is also the only interpretation consistent with the factor r_c²/r(x)² = 1/0.75² = 1.78 quoted for rp = 0.25×rc in Fig. 5. As written, Eq. (15) makes Eq. (16) and the derived 15%/30%/43% Venturi predictions irreproducible. Please correct the notation (for example, introduce r_lumen = rc − rp) and specify the integration domain for each region.","section":"§III-B, Eq. (15)"}],"minor_comments":[{"comment":"The condition \"for t ≥ d/u0t\" is dimensionally wrong; it should be t ≥ d/u0. In addition, the symbol a in the formula is never defined and should be set to a = rc.","section":"§III-A3, Eq. (14)"},{"comment":"The text says that particles \"arrive later at the TX\" for release time tED; the receiver is at the RX, not the TX, so this should read \"arrive later at the RX.\"","section":"§IV-B"},{"comment":"The caption says the power-law and Herschel-Bulkley models have been fitted to the simulated profile at x = 0.01 m, but the fitted parameter values and the fitting procedure are not reported; please provide them.","section":"Fig. 5 caption"},{"comment":"The domain intervals for the four regions in the piecewise radius function are not stated; please define them explicitly in terms of lc, lp,outer, and lp,inner.","section":"§III-B, Eq. (15)"},{"comment":"The simulations are said to use the Casson model while the analytical comparison uses power-law and Herschel-Bulkley models; please clarify whether the same rheological parameters were used in the pulsatile simulations and how the Casson parameters relate to the analytical model parameters.","section":"§II and §III"},{"comment":"The text mentions an \"irregularity in the upper third of the PS plot curve\" but does not explain its physical origin; please annotate or discuss this feature.","section":"Fig. 6"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for a molecular-communication venue. The main risk is overclaiming the severe-stenosis results under an explicitly laminar framework; the authors should either rescope the claims or provide a transitional-flow validation. The Venturi notation error is easily fixed. I do not see citation-pattern problems: the self-citations are to the authors' own baseline, solver, and dataset."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things you should know. First, the closed-form CIRs for power-law and Herschel-Bulkley blood models are genuinely new in this MC setting, and I checked the derivations: Eqs. (8) and (14) follow from the stated flow-dominated, uniform-release assumptions. Second, the pulsatile-flow simulations are a real step beyond their earlier constant-flow paper. The qualitative result — plaque size and release phase within the cardiac cycle change arrival statistics — is plausible and visually supported.\n\nThe paper earns credit for shipping reproducible artifacts: the dataset is public, the derivations are parameter-free, and the Fig. 5 comparison honestly shows that analytical profiles deviate from simulation, concluding that simulation is needed for the full scenario. Self-citations to [6], [7], and [25] are used as baselines and prior work, not to inflate novelty. That is fine.\n\nThe soft spots are real but not fatal. The stress-test concern lands: at rp = 0.75rc the throat diameter is 1.5 mm, continuity gives about 5.5 m/s at mean flow, and the peak-systole factor of 2–3 puts the throat Reynolds number around 4000–6000, well into transitional/turbulent territory. Flow separation and recirculation at that stenosis would change particle trapping and arrival dispersion. The laminar assumption is stated explicitly and deferred to future work, so the paper is honest, but the 43% speedup and the claim that peak-systole release gives fastest arrivals are laminar-model predictions, not robust hemodynamic signatures. Also, Fig. 6 shows single runs with no error bars, and the analytical-model comparison in Fig. 5 fits the models to the simulation before comparing. These issues weaken, but do not destroy, the central qualitative claim. The abstract's phrase \"crucial indicators\" outruns the evidence; \"potential indicators\" would match what is actually shown.\n\nWho is this for? MC/IoBNT modelers, especially people building simulation pipelines for in-body communication channels. Not clinicians, and not yet a detection method. The paper deserves a serious referee: the derivations and dataset warrant evaluation, and the laminar/turbulent question is exactly what a good reviewer should push on. Revisions should either add a transitional-flow comparison or restrict the detection claims to milder stenoses, and they should add uncertainty quantification to the simulation curves.","headline":"Solid analytical extension of the authors' prior MC plaque work, with correct new CIR derivations; the quantitative claims for severe stenosis are weakened by an acknowledged but likely violated laminar-flow assumption.","tokens_in":11417,"tokens_out":2186,"would_cite":false,"duration_ms":26366,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper argues that plaque-induced changes in blood flow and nanoparticle arrival times form a measurable molecular-communication signature for detecting atherosclerosis non-invasively.","keywords":["atherosclerosis","molecular communication","channel impulse response","Internet of Bio-Nano Things","non-Newtonian blood flow","pulsatile flow","plaque detection","OpenFOAM simulation"],"falsifier":"Run a turbulence-resolving or experimentally validated simulation of 50-nm SPIONs through a 75% stenosis at peak systole with the same carotid waveform and compare the first-arrival traversal speedup and full CIR against the unobstructed case; if the roughly 43% speedup or the clean separation of CIRs by plaque radius is not reproduced, the proposed detection signature does not survive realistic high-systole flow.","tokens_in":10197,"feed_emoji":"🫀","tokens_out":5560,"duration_ms":53621,"temperature":0.7,"pith_summary":"This paper tries to establish that atherosclerosis can be sensed through molecular communication itself: as plaque narrows a carotid artery, the channel that flowing nanoparticles experience changes in measurable ways, and a receiver counting particle arrivals can pick up the difference. The authors simulate a plaque-obstructed blood vessel with pulsatile blood flow and find that traversal speed rises by 5.5%, 33%, and 43% as the plaque radius grows from 25% to 50% to 75% of the vessel radius. They also find that timing matters: particles released at peak systole arrive fastest, while early- and late-diastole releases arrive slower or approach the constant-flow baseline. If these channel changes hold in real arteries, an internet-of-bio-nano-things network could detect and perhaps locate plaque non-invasively by watching how injected nanoparticles move.","feed_headline":"Plaque changes blood's molecular channel by up to 43%","feed_subtitle":"Simulations show stenosis and release timing alter particle arrival, a signature future bio-nano sensors could detect.","key_machinery":"The central object is the channel impulse response $h(t)$, the fraction of released nanoparticles that have crossed the receiver plane at time $t$. The argument is carried by three interlocking pieces: the flow-dominated-regime CIR derivation adapted from pipe-flow molecular communication, with closed-form expressions for Newtonian, power-law, and Herschel-Bulkley fluids; the Venturi-effect model that treats the plaque as a piecewise radial constriction and predicts traversal-time reduction $T=\\int_0^{l_c} dx/u(x)$; and the OpenFOAM MPPIC particle simulation driven by a realistic carotid pulsatile flow profile with particle release at peak systole, early diastole, and late diastole. The simulation supplies what the analytical models cannot: the lasting asymmetric flow profile downstream of the plaque.","core_discovery":"In the paper's own terms, the discovery is that plaque-induced stenosis leaves a detectable fingerprint in the molecular communication channel impulse response. Using a CFD simulation with a non-Newtonian Casson blood model, a human-carotid pulsatile inlet waveform, and superparamagnetic iron-oxide nanoparticles as information carriers, the authors show that the fraction of received particles over time separates by plaque size: first-arrival traversal speed increases by 5.5% for $r_p=0.25\\,r_c$, 33% for $r_p=0.5\\,r_c$, and 43% for $r_p=0.75\\,r_c$, close to the 15%, 30%, and 43% predicted by their Venturi model. They further derive closed-form channel impulse responses for power-law and Herschel-Bulkley non-Newtonian models, show that the two models nearly coincide in this parameter regime, and demonstrate that the simulated flow profile stays skewed and asymmetric downstream of the plaque, in line with in-vivo observations. The intended upshot is that flow-profile shape and CIR changes are usable indicators for early, non-invasive plaque detection.","pith_inferences":["An implicit design consequence is that an IoBNT detector must either control or measure the release phase: at small plaque sizes, the shift between systole and diastole release can be as large as the plaque-induced shift, so phase-unaware measurements could confound plaque size estimation.","The Venturi traversal-time formula could be inverted: given a measured speedup and known channel length, a device could estimate the effective plaque radius without imaging.","A natural extension would test wall compliance via fluid-solid interaction; the published dataset already includes a first FSI trial, so the signature's robustness to moving walls is checkable.","If turbulence modeling is added, the same simulation pipeline could quantify whether the peak-systole detection window remains the best or shifts to another phase."],"forward_implications":["A plaque that narrows the vessel to 75% of its normal radius can speed up nanoparticle traversal by about 43%, making traversal time a direct physical indicator of stenosis severity.","Particle release should be synchronized with peak systole for the fastest transport through an obstructed vessel, while early-diastole release produces the slowest arrivals.","The largest separation between plaque and no-plaque channel impulse responses occurs at peak systole, making that phase the most useful detection window for future IoBNT sensors.","Analytical non-Newtonian CIRs and the Venturi traversal-time estimate give quick approximations that match simulation well for larger plaques, even though full CFD is needed to capture the skewed downstream flow profile."],"supporting_citations":[{"why":"Provides the previous simulation scenario and baseline constant-flow results that this paper extends to pulsatile flow.","marker":"[6]"},{"why":"Supplies the closed-form Hagen-Poiseuille solutions for power-law and Herschel-Bulkley fluids used to derive analytical CIRs.","marker":"[15]"},{"why":"Gives the flow-dominated CIR derivation for pipe flow that the paper adapts to the plaque channel.","marker":"[20]"},{"why":"Supplies in-vivo evidence of skewed asymmetric carotid flow profiles with plaque, supporting the simulated metric.","marker":"[24]"},{"why":"Supplies the OpenFOAM MPPIC simulation methodology and its experimental validation.","marker":"[25]"},{"why":"Supplies the measured human carotid blood-flow waveform used for the pulsatile inlet boundary condition.","marker":"[12]"},{"why":"Releases the simulation dataset including pulsatile flow modeling for reproducibility.","marker":"[7]"}],"fun_headline_variants":["Blood's molecular channel reveals plaque size","Plaque leaves a 43% signature in molecular blood flow","Bio-nano sensors could spot plaque via channel changes","Plaque skews blood's molecular signals for early detection","Molecular communication exposes plaque in arteries"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result rests on the assumption that flow remains laminar even at peak systole through a 75% stenosis; turbulence would change the arrival statistics and could erase the speedup signature.","fun_headline_variants_meta":{"raw":{"variants":["Blood's molecular channel reveals plaque size","Plaque leaves a 43% signature in molecular blood flow","Bio-nano sensors could spot plaque via channel changes","Plaque skews blood's molecular signals for early detection","Molecular communication exposes plaque in arteries"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000379,"raw_usage":{"total_tokens":2039,"prompt_tokens":994,"completion_tokens":1045,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":610,"completion_tokens_details":{"reasoning_tokens":971}},"tokens_in":610,"tokens_out":1045,"duration_ms":10802,"temperature":1.0,"reasoning_tokens":971,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:39:52.976867+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a turbulence-resolving or experimentally validated simulation of 50-nm SPIONs through a 75% stenosis at peak systole with the same carotid waveform and compare the first-arrival traversal speedup and full CIR against the unobstructed case; if the roughly 43% speedup or the clean separation of CIRs by plaque radius is not reproduced, the proposed detection signature does not survive realistic high-systole flow.","supporting_citations":[{"cited_title":"A Molecular Communication Perspective on Detecting Arterial Plaque Formation,","cited_arxiv_id":null,"evidence_quote":"Provides the previous simulation scenario and baseline constant-flow results that this paper extends to pulsatile flow."},{"cited_title":"Predicting laminar–turbulent transition in Poiseuille pipe flow for non-Newtonian fluids,","cited_arxiv_id":null,"evidence_quote":"Supplies the closed-form Hagen-Poiseuille solutions for power-law and Herschel-Bulkley fluids used to derive analytical CIRs."},{"cited_title":"Ex- perimental System for Molecular Communication in Pipe Flow With Magnetic Nanoparticles,","cited_arxiv_id":null,"evidence_quote":"Gives the flow-dominated CIR derivation for pipe flow that the paper adapts to the plaque channel."},{"cited_title":"In vivo three-dimensional blood velocity profile shapes in the human common, internal, and external carotid arteries,","cited_arxiv_id":null,"evidence_quote":"Supplies in-vivo evidence of skewed asymmetric carotid flow profiles with plaque, supporting the simulated metric."},{"cited_title":"OpenFOAM Simulation of Microfluidic Molecular Communications: Method and Experimental Validation,","cited_arxiv_id":null,"evidence_quote":"Supplies the OpenFOAM MPPIC simulation methodology and its experimental validation."},{"cited_title":"Characterization of Common Carotid Artery Blood-Flow Waveforms in Normal Human Subjects,","cited_arxiv_id":null,"evidence_quote":"Supplies the measured human carotid blood-flow waveform used for the pulsatile inlet boundary condition."},{"cited_title":"Dataset for Advanced Plaque Modeling for Atherosclerosis Detection using Molecular Communication","cited_arxiv_id":null,"evidence_quote":"Releases the simulation dataset including pulsatile flow modeling for reproducibility."}],"review_version":1}