{"id":"469bb2bc-fd78-415e-acdd-492cc696b313","arxiv_id":"2506.00803","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"An approximate 3D channel model for molecular communication with a wall-attached absorbing ring receiver in a tube, derived under flow-dominated conditions and checked against particle simulations.","lead":"This paper derives approximate analytical formulas for how molecules released in a tube move with the flow and diffuse toward a ring-shaped absorbing receiver attached to the inner wall. The formulas give communication engineers a closed-form way to compute molecular signals in pipes and blood-vessel-like environments.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The flow-dominated one-crossing heuristic and the deterministic-time switching in Section III-B are not validated against any test where axial diffusion is non-negligible; the paper only reports high-Pe regimes, so the central approximation's error bound is untested.","rationale":"The reader's weakest_assumption correctly identified the deterministic-time switching and the conditional-density reinterpretation as the load-bearing heuristic. My stress-test agrees and sharpens it: the paper's own validation only spans Pe from 14.3 to 600, all satisfying the flow-dominated assumption, so the approximation error is not probed where the heuristic could fail. The concrete test I propose directly targets that gap by lowering Pe and by computing the multiple-crossing fraction. I also note the paper includes an explicit limitation statement in the conclusion ('extending the framework to diffusion-dominated regimes remains an open challenge'), which I weigh as honest scope-setting but which does not obviate the need to quantify the error within the claimed flow-dominated regime. No machine-checked proof or code is provided, and the derivations are heuristic with omitted details; independent numerical verification of the key formulas (e.g., the exact coefficients in (16) and (18)) would also be valuable, but the single most load-bearing concern remains the uncontrolled deterministic-switching approximation. Therefore I recommend keeping the verdict CONDITIONAL, consistent with the reader, with the condition being the demonstration that the approximation error remains small as Pe decreases toward the boundary of the flow-dominated regime.","tokens_in":9714,"tokens_out":1892,"duration_ms":15821,"concrete_test":"Run particle-based simulations for a set of parameter combinations with much lower Pe, e.g., v=200 um/s, D=700 um^2/s, rho=10 um (Pe=2.86), and also v=500 um/s, D=700 um^2/s (Pe=7.14), keeping d1, d2, N, M as in the paper. Compute the NRMSE between the theoretical R(t) from (25) and the simulated absorption ratio. If the NRMSE degrades significantly below the 0.95 mean reported in Table III, the flow-dominated deterministic-switching approximation is not robust and the paper should state a quantitative Pe validity range.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that the approximate concentration (20), arrival probability (25), and arrival rate (26) describe the channel response in the flow-dominated regime. The load-bearing heuristic in Section III-B is that each molecule crosses z=d1 and z=d2 exactly once at deterministic times t1≈d1/v and t2≈d2/v, so the 3D problem with z-dependent boundary conditions can be replaced by a 2D radial diffusion with BCs switching at fixed times. This is an uncontrolled approximation: even for high Pe, the axial Gaussian (19) assigns nonzero probability to molecules that cross d1 or d2 at times far from d1/v or d2/v, or that cross multiple times due to axial diffusion. The paper provides no quantitative bound on the error introduced by this deterministic-switching step. The validation reported in Table III (mean NRMSE 0.9504) covers only Pe from 14.3 to 600, all with flow strongly dominating diffusion along the tube axis in the sense that axial displacement by drift over the Rx length (20-1000 um) far exceeds the diffusive spread. Table III additionally shows NMSE 0.9966 and NRMSE mean 0.9504, but the least accurate case has absorption ratio as low as 0.0634, suggesting the model's match is weakest precisely where the survival analysis matters most. The paper does not report any test in a regime where axial diffusion is comparable to flow advection, nor does it provide an analytic error bound as a function of Pe, D, d1, d2, or v. Without such a check, the claimed validity of the approximate solution in the flow-dominated regime is not established beyond the narrow parameter range tested.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a three-dimensional channel model for molecular communication in a cylindrical tube with an absorbing ring-shaped receiver attached to the inner wall. The authors formulate an advection-diffusion PDE with heterogeneous boundary conditions, then derive approximate analytical expressions for the molecular concentration, the arrival probability, and the arrival rate under an assumed flow-dominated regime. The key approximation replaces the z-dependent absorbing boundary condition with a time-dependent boundary condition that switches at deterministic crossing times t1 ≈ d1/v and t2 ≈ d2/v, factorizes the solution into radial and axial parts, and then applies a survival analysis that integrates over inverse-Gaussian crossing times. The model is validated against particle-based simulations for six parameter sets, with mean NRMSE 0.9504 and mean NMSE 0.9966 reported in Table III.","tokens_in":10089,"tokens_out":17350,"duration_ms":168433,"significance":"If the approximations can be properly justified, the paper would be a useful contribution: it claims to be the first analytic treatment of an absorbing ring-shaped receiver in a cylindrical MC channel, and the final formulas (20), (25), and (26) are explicit and involve no fitted parameters, with the coefficients determined solely by the initial condition, boundary conditions, and geometry. The particle-based validation is external to the derivation, and the comparison metrics are reported in detail. However, the significance is currently conditional because the derivation contains a load-bearing gap in the survival analysis and an unquantified approximation in the flow-dominated regime, so the contribution is not yet established at the level required for publication.","major_comments":[{"comment":"The Laplacian operator is misprinted in the governing PDE: the radial part should contain (1/r)∂/∂r, but instead it contains (1/r)∂/∂θ and omits the former term. As written, Eq. (2) is not the cylindrical advection-diffusion equation, and the Bessel eigenfunctions used in Eqs. (13)-(18) solve a different operator. The subsequent mathematics appears to be consistent with the correct operator, so this is likely a typo, but it must be corrected because the stated PDE is the foundation of the formulation.","section":"Section III-A, Eqs. (2) and (9)"},{"comment":"The reinterpretation of the unconditioned concentration (20) as the conditional probability density p_{R,Θ,Z|T1,T2}(r,θ,z|t;t1,t2) given the first passage times T1=t1 and T2=t2 is not justified and is in fact mathematically incorrect: the axial factor (19) is the marginal density of Z(t), whereas the density conditional on T1=t1 has, for t>t1, mean d1+v(t-t1) and variance 2D(t-t1), and conditioning on T2=t2 imposes further bridge-type constraints. Substituting the marginal density for the conditional density changes the survival function in (21) and therefore the averaged formulas (25) and (26). The authors must either derive the correct conditional axial density or provide a rigorous asymptotic argument, with an explicit error bound, showing that the unconditional Gaussian can be used in the flow-dominated limit.","section":"Section III-C, Eqs. (19)-(22)"},{"comment":"The deterministic-switching and one-crossing heuristic is an uncontrolled approximation. The paper provides no quantitative characterization of the flow-dominated regime and no error bound as a function of v, D, d1, and d2; the validation in Table II covers only six scenarios with Pe values from 14.3 to 600, and the least accurate case has an absorption ratio as low as 0.0634, where the survival analysis is most consequential. Please add an analysis of the approximation error (e.g., via an axial Péclet number v(d2-d1)/D or a bound on the multiple-crossing probability) or extend the simulation study to regimes where axial diffusion is non-negligible, so that the claimed validity domain of the model is actually established.","section":"Section III-B and Section IV-B"}],"minor_comments":[{"comment":"The phrases 'details omitted' appear at several key points, notably after Eq. (18) and before Eqs. (23) and (26); please move these derivations to an appendix or supplementary material so that the results are independently verifiable.","section":"Section III-B and III-C"},{"comment":"The expression ∫_0^t I_{[t1,t2]}(τ)dτ denotes the length of the intersection of [0,t] and [t1,t2]; it would be helpful to state explicitly that it equals 0 for t<t1, t-t1 for t1≤t≤t2, and t2-t1 for t>t2.","section":"Eq. (23)"},{"comment":"The heading 'Tx-Rx Distance' is misleading because d1 is the distance to the upstream edge of the ring, not to the center of the receiver; consider renaming it to 'distance to the start of the ring'.","section":"Table I"},{"comment":"The same symbol t is used both for the time argument of the PDE solution and for the survival-time variable in Section III-C; using a different symbol for the integration variable would improve readability.","section":"Notation"},{"comment":"The figure captions could state that solid lines are theory and crosses are simulation, and may also mention that the theoretical curves use N=M=10 terms; this is currently clear from the text but would be useful in the captions themselves.","section":"Fig. 2"}],"recommendation":"major_revision","confidential_remarks":"The paper's novelty claim is plausible and the simulation results are encouraging, but the survival-analysis derivation in Section III-C is built on an unjustified identification of the unconditional concentration with the conditional density given the first passage times. This is a mathematical error rather than a mere missing bound, so the revision needs to be substantial. If the authors can correct that step or recast the result as an asymptotic approximation with a stated error order, the paper could be viable for this journal; otherwise, the central formulas for arrival probability and rate are not supported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a genuine first: an approximate analytical channel model for an absorbing ring-shaped receiver in a cylindrical tube with heterogeneous wall conditions, and it backs the model with particle-based simulation. The central trick is a time-segmented separation of variables: the radial/angular diffusion sees reflecting BCs before and after the ring and an absorbing BC on the ring, with the switch times fixed by the drift times d1/v and d2/v; the axial factor is the free-space Gaussian. The geometry is new relative to earlier analytical tube models with homogeneous BCs, and the final series expressions for concentration, arrival probability, and arrival rate are explicit and checkable. The simulation agreement is respectable: mean NRMSE around 0.95 across six parameter sets, and the paper honestly points out that low-absorption cases need larger truncation orders.\n\nThe soft spots are real but not disqualifying. The switching heuristic is uncontrolled: there is no error bound in Pe, D, d1, or d2, and all tests are at high Pe. The stress-test note is right that no test addresses a regime where axial diffusion is comparable to drift; I would have liked at least a moderate-Pe data point to probe where the approximation degrades. The paper does restrict itself to the flow-dominated regime, so this is a boundary of the stated claim rather than a contradiction of it.\n\nSecond, the PDE written in (2) and (9) is wrong as displayed: it omits the (1/r) partial/r term and includes a spurious (1/r) partial/theta. The subsequent Bessel eigenfunctions are consistent with the correct operator, so the error is a typo, but it will trip readers and needs fixing.\n\nThird, the reinterpretation of the concentration as a conditional density given T1 and T2 is stated without rigorous justification. It is a modeling assumption inside an already heuristic construction, so I can live with it, but the paper should flag it as an approximation rather than present it as a formal identity. The omissions of derivation details (labeled 'details omitted' in III-B and III-C) also make independent checking slower, though the final formulas look plausible.\n\nWho is this for? Researchers modeling MC in blood-vessel-like or pipe-like tubes, especially those who want a closed-form channel response for a wall-attached ring receiver. It is an engineering model with clear scope, not a rigorous boundary-value solution.\n\nI would send it to peer review. A serious referee can request the corrected PDE, an error analysis or additional Pe values, and a more careful statement about the conditional-density step. The core idea is solid and the simulation support is decent, so it deserves referee time.","headline":"A useful, honest engineering model for a new receiver geometry in MC tubes; the heuristic time-switching step is uncontrolled and the PDE has a typo, but the simulation support is decent and the paper deserves peer review.","tokens_in":10555,"tokens_out":4311,"would_cite":true,"duration_ms":40542,"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 derives the first analytical channel model for molecular communication in a tube whose receiver is an absorbing ring on the wall, using a flow-dominated approximation that splits the problem into three time periods and yields…","keywords":["molecular communication","channel modeling","tubular environment","absorbing ring receiver","heterogeneous boundary conditions","advection-diffusion equation","survival analysis","flow-dominated regime"],"falsifier":"Run particle-based simulations for parameter sets with decreasing Péclet number (for instance, holding tube radius and flow fixed while increasing the diffusion coefficient from 100 to 400 to 700 µm²/s, and then beyond) and compare the simulated arrival probability with formula (25); if the normalized RMSE worsens sharply as Pe drops below the values in Table II, that locates the point where the deterministic-crossing heuristic fails.","tokens_in":9482,"feed_emoji":"🧪","tokens_out":9587,"duration_ms":84362,"temperature":0.7,"pith_summary":"Molecular communication inside fluid-filled tubes needs channel models that say how many emitted molecules arrive at a receiver over time. This paper proposes the first analytical model for a receiver shaped as an absorbing ring on the tube wall, where the wall reflects molecules except on the ring. Because the tube wall behaves differently on different axial stretches, the advection–diffusion equation carries heterogeneous boundary conditions, which are normally hard to solve; the paper's key move is to split time into three periods using the assumption that flow dominates diffusion. That yields approximate formulas for molecular concentration, arrival probability, and arrival rate, which the paper reports agree with particle-based simulations across six parameter sets. If correct, the formulas give an inexpensive way to design and evaluate tube-based molecular communication systems, such as sensors in blood vessels or fluid pipelines.","feed_headline":"First channel formula for a wall-mounted ring receiver","feed_subtitle":"Closed-form formulas predict how many molecules the ring captures over time, matching simulations in six flow settings.","key_machinery":"The load-bearing device is a piecewise separation of variables in the flow-dominated regime. The heuristic is that each molecule, carried by the axial flow, crosses the planes $z=d_1$ and $z=d_2$ exactly once, at approximately $t_1=d_1/v$ and $t_2=d_2/v$; therefore the radial-angular part of the concentration evolves under a reflecting wall before $t_1$, an absorbing wall between $t_1$ and $t_2$, and a reflecting wall after $t_2$, with continuity patched at the switching times. The radial solutions are Bessel-function eigenfunction expansions (13), (15), (17), and the axial factor is the free drifting Gaussian (19). The second piece is survival analysis: the crossing times $T_1$ and $T_2$ are modeled as inverse Gaussian random variables, and the concentration is treated as the conditional density of a molecule's position given $T_1$ and $T_2$, which converts the PDE solution into arrival probability (25) and arrival rate (26).","core_discovery":"The paper's central claim is that the channel response of an absorbing ring-shaped receiver in a cylindrical tube can be approximated by separating the advection–diffusion problem radially with boundary conditions that switch at the deterministic times $t_1 \\approx d_1/v$ and $t_2 \\approx d_2/v$, while the axial factor is the Gaussian (19). Reinterpreting the resulting concentration as a conditional density given these crossing times, and averaging over the inverse-Gaussian crossing-time distributions, produces expressions (20), (25), and (26) for concentration, arrival probability, and arrival rate. The paper states this is the first theoretical study of heterogeneous boundary conditions in cylindrical molecular-communication environments and validates the expressions against particle-based simulations over six parameter sets, reporting a mean NRMSE of 0.9504. On the paper's terms, the result makes the channel impulse response of a wall-attached ring receiver analytically computable in the flow-dominated regime.","pith_inferences":["A natural testable extension is to apply the same switching-time heuristic in channels with multiple absorbing patches or non-circular cross-sections, since the paper's framework suggests the switch times would still be $d_1/v$ and $d_2/v$ but does not test that case.","Because the paper explicitly leaves diffusion-dominated regimes open, a direct follow-up is to scan the Péclet number downward and map where the reported fit degrades, quantifying the regime boundary.","If the model carries over to blood-vessel-scale parameters, wall-attached ring receivers would give implantable molecular communication a physically fixed alternative to floating point receivers, though the paper does not itself make this engineering claim."],"forward_implications":["For the six parameter combinations tested in Table II, the paper reports that the approximate arrival probability tracks particle-based simulations with a mean NRMSE of 0.9504.","The arrival probability at the ring receiver saturates below one, quantifying the fraction of molecules that are swept past the ring before being absorbed.","The derivation framework is claimed to extend to tubes with more than three heterogeneous wall sections and to Robin boundary conditions.","The arrival-rate expression requires only a one-dimensional numerical integration once inverse-Gaussian CDFs are available, making the channel impulse response practical to evaluate."],"supporting_citations":[{"why":"Supplies the standard advection–diffusion PDE formulation (its Eq. (17)) and the flow-dominance/Péclet-number criterion that defines the regime of validity.","marker":"[3]"},{"why":"Provides the cylindrical-tube separation-of-variables approach with Bessel eigenfunctions that the three-period radial solution adapts.","marker":"[8]"},{"why":"Provides the tube channel model with Robin boundary conditions whose separation-of-variables technique the piecewise radial solution builds on.","marker":"[10]"},{"why":"Supplies the exact cylindrical channel model for homogeneous wall conditions that this paper generalizes to heterogeneous conditions.","marker":"[12]"},{"why":"Is the earlier theoretical treatment of heterogeneous boundaries (in spherical geometry) that this work extends to cylindrical tubes.","marker":"[18]"},{"why":"Provides the experimental pipe-flow scenario and parameter ranges (initial distribution, laminar flow, Reynolds number) used for the simulations and validation.","marker":"[19]"},{"why":"Supplies the inverse Gaussian first-passage-time distribution and survival-analysis framework used for $T_1$ and $T_2$ in the arrival probability and rate.","marker":"[21]"},{"why":"Provides the strong Markov property of Brownian motion, letting the paper treat $T_2 - T_1$ as independent of $T_1$.","marker":"[22]"}],"fun_headline_variants":["Ring receiver channel response now computable in tubes","First 3D model for ring receivers in molecular tubes","Tube molecular channel formula with ring receiver","New analytical model for ring receivers in tubular MC","First ring-receiver channel model for tube molecular comms"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole derivation depends on the assumption that axial flow is so strong relative to diffusion that every molecule crosses the two planes bounding the ring exactly once, at the deterministic times $d_1/v$ and $d_2/v$, so that the wall condition can be treated as switching at fixed times; if diffusion is not negligible, the piecewise separation breaks down.","fun_headline_variants_meta":{"raw":{"variants":["Ring receiver channel response now computable in tubes","First 3D model for ring receivers in molecular tubes","Tube molecular channel formula with ring receiver","New analytical model for ring receivers in tubular MC","First ring-receiver channel model for tube molecular comms"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000222,"raw_usage":{"total_tokens":1431,"prompt_tokens":899,"completion_tokens":532,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":515,"completion_tokens_details":{"reasoning_tokens":459}},"tokens_in":515,"tokens_out":532,"duration_ms":5747,"temperature":1.0,"reasoning_tokens":459,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:58:19.107518+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run particle-based simulations for parameter sets with decreasing Péclet number (for instance, holding tube radius and flow fixed while increasing the diffusion coefficient from 100 to 400 to 700 µm²/s, and then beyond) and compare the simulated arrival probability with formula (25); if the normalized RMSE worsens sharply as Pe drops below the values in Table II, that locates the point where the deterministic-crossing heuristic fails.","supporting_citations":[{"cited_title":"Channel Modeling for Diffusive Molecular Communication-A Tutorial Review,","cited_arxiv_id":null,"evidence_quote":"Supplies the standard advection–diffusion PDE formulation (its Eq. (17)) and the flow-dominance/Péclet-number criterion that defines the regime of validity."},{"cited_title":"Diffusive Molecular Communica- tion in Biological Cylindrical Environment,","cited_arxiv_id":null,"evidence_quote":"Provides the cylindrical-tube separation-of-variables approach with Bessel eigenfunctions that the three-period radial solution adapts."},{"cited_title":"Modeling Molecular Communications in Tubes with Poiseuille Flow and Robin Boundary Condition,","cited_arxiv_id":null,"evidence_quote":"Provides the tube channel model with Robin boundary conditions whose separation-of-variables technique the piecewise radial solution builds on."},{"cited_title":"Transfer Function Models for Cylindrical MC Channels with Diffusion and Laminar Flow,","cited_arxiv_id":null,"evidence_quote":"Supplies the exact cylindrical channel model for homogeneous wall conditions that this paper generalizes to heterogeneous conditions."},{"cited_title":"Channel Modeling of Diffusive Molecular Communication with Heterogeneous Spherical Boundaries,","cited_arxiv_id":null,"evidence_quote":"Is the earlier theoretical treatment of heterogeneous boundaries (in spherical geometry) that this work extends to cylindrical tubes."},{"cited_title":"Ex- perimental System for Molecular Communication in Pipe Flow with Magnetic Nanoparticles,","cited_arxiv_id":null,"evidence_quote":"Provides the experimental pipe-flow scenario and parameter ranges (initial distribution, laminar flow, Reynolds number) used for the simulations and validation."},{"cited_title":"Molecular Communi- cation in Fluid Media: The Additive Inverse Gaussian Noise Channel,","cited_arxiv_id":null,"evidence_quote":"Supplies the inverse Gaussian first-passage-time distribution and survival-analysis framework used for $T_1$ and $T_2$ in the arrival probability and rate."},{"cited_title":"M ¨orters and Y","cited_arxiv_id":null,"evidence_quote":"Provides the strong Markov property of Brownian motion, letting the paper treat $T_2 - T_1$ as independent of $T_1$."}],"review_version":1}