{"id":"8407635c-8e66-4cb4-b5f7-5cb02a331e06","arxiv_id":"2506.22695","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A hybrid particle-continuum simulation method couples protein drift-diffusion to fluctuating concentration and temperature fields in membranes, and is demonstrated on protein localization, thermal gradient sensing, and hot Brownian motion.","lead":"The authors build a computer simulation framework that tracks individual membrane proteins while also modeling how surrounding chemical concentrations and temperatures fluctuate over time. The goal is to study biological processes like thermal sensing, protein positioning, and hot-particle escape in membranes out of equilibrium.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Direct substitution shows Appendix A K^(2) is inconsistent with Eq. 3: the cross-coupling term has the opposite sign and the mechanical-heating term differs by a factor c0/q.","rationale":"The reader's weakest assumption was that the Appendix A operators are asserted rather than derived. I checked the mapping directly and found a concrete mismatch, not merely a missing derivation: the theta_C equation obtained from Eq. 31 differs from the printed Eq. 3 in the sign of the concentration-temperature coupling and in the factor multiplying the mechanical heating term. This is internal inconsistency, not a disagreement with external consensus, and it is load-bearing because Eq. 7 converts the same K operators into fluctuation amplitudes. The Appendix C covariance test validates the discretization and integrator against 2kB K, but it does not validate that K is the dissipative operator of the stated model; a covariance test of a K that mismatches Eq. 3 would pass while the physics written in Eq. 3 is wrong. I therefore keep the reader's conditional verdict, with the explicit added condition that Eqs. 3 and 31 be reconciled and affected simulations rerun. If the authors demonstrate that Eq. 3 is only a typographical error and the code implements the correct sign and factor, the central self-consistency claim could stand.","tokens_in":19688,"tokens_out":17163,"duration_ms":170464,"concrete_test":"Perform the term-by-term expansion of K^(2)DS (Eq. 31 with DS from Eq. 15) and compare with the deterministic parts of Eqs. 2-3; also check the GENERIC degeneracy DE^T K^(2) = 0. If, as direct substitution shows, the theta_C update contains +c0 grad Phi : kappa_bar grad q / cC and c0 q kappa_bar |grad Phi|^2/(cC theta_C), then Eq. 3 must be corrected and the affected simulations (thermal-gradient sensing and hot Brownian motion) rerun, reporting the change in sensed signals and escape times.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Expanding the deterministic part of Eq. 6 for the theta_C component using K^(2) from Eq. 31 and DS from Eq. 15 gives dtheta_C/dt = div(kappa0 grad theta_C) + c0 grad Phi : kappa_bar grad q / cC + c0 q kappa_bar |grad Phi|^2 / (cC theta_C), with kappa_bar = theta_C/gamma. Equation 3 instead has -c0 grad Phi : kappa_bar grad q / cC and +(c0 grad Phi):(1/gamma)(c0 grad Phi)/cC, i.e. the sign of the first coupling term is opposite and the second term is c0^2/gamma |grad Phi|^2/cC rather than c0 q/gamma |grad Phi|^2/cC. Since the stochastic driving fields in Eq. 7 are defined through B^(j)B^(j),T = 2kB K^(j), a mismatch between the stated model and K^(2) breaks the claimed fluctuation-dissipation balance. If Eq. 3 is the intended model, K^(2) is not its dissipative operator; if K^(2) is correct, Eq. 3 misstates the deterministic drift. Either way, the central self-consistency claim is not supported by the printed equations, and the deferred derivation cannot be patched without changing one of the two.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a GENERIC-based stochastic hybrid model that couples a discrete protein position X to a fluctuating concentration field q and three temperature variables (protein temperature θP, membrane temperature field θC, and interfacial temperature θI). Fluctuations are generated from dissipative operators K^(j) via the fluctuation-dissipation relation B^(j)B^(j),T = 2kB K^(j), with analytic factorizations for efficient sampling. The numerical discretization is tested with a heat-equation convergence study (second-order) and a covariance test (error 6.8e-9). Three applications are presented: protein localization in concentration gradients, thermal-gradient sensing with a temporal filter, and escape kinetics from energy wells under local heating that yields hot Brownian motion.","tokens_in":20016,"tokens_out":8504,"duration_ms":86382,"significance":"If the construction were internally consistent, the framework would be a valuable contribution to non-equilibrium membrane simulation: it couples particle and field fluctuations while preserving a GENERIC structure, and the numerical validation is strong (second-order convergence; covariance error 6.8e-9), with the methodology described in sufficient detail to be reproduced. The applications are plausible demonstrations of the intended physics, but they are qualitative rather than quantitative biological predictions. The paper's central claim of self-consistency is, however, not supported by the printed equations because the dissipative operator K^(2) does not match the deterministic dynamics of Eq. (3); this must be resolved before the manuscript can be assessed.","major_comments":[{"comment":"Direct substitution of K^(2) from Eq. (31) and DS from Eq. (15) into Eq. (6) gives, for the θC component, +c0 ∇Φ : ¯κ ∇q / cC and c0 q ¯κ |∇Φ|^2 / (cC θC), with ¯κ = θC/γ. Equation (3) instead has −c0 ∇Φ : ¯κ ∇q / cC and c0^2 |∇Φ|^2 / (γ cC). The concentration-coupling term has the opposite sign and the mechanical-heating term differs by a factor c0/q. Because the stochastic fields are defined through B^(j)B^(j),T = 2kB K^(j) in Eq. (7), the simulated fluctuations are consistent with a deterministic model different from Eqs. (1)–(4). This contradicts the statement that Appendix A provides the operators for the model in equations 1–3 and undermines the central self-consistency claim.","section":"Appendix A, Eq. (31) vs Eq. (3)"},{"comment":"The numerical validation does not test the mapping from K^(j) to Eqs. (1)–(4). Figure 9 verifies second-order convergence for the scalar heat equation only, and Figure 10 verifies that the sampled increments reproduce the covariance 2kB K^(j) of the discretized operators. Neither test can detect a sign or factor error in K^(2) itself. Once the inconsistency in Major Comment 1 is resolved, the authors should add a direct comparison between Σ_j K^(j) DS and the deterministic right-hand sides of Eqs. (1)–(4).","section":"Appendix C, Figs. 9–10"}],"minor_comments":[{"comment":"Figures 4 and 8 do not report error bars or the number of independent realizations; for stochastic simulations these are needed to judge whether the reported differences (e.g., in escape times) are statistically significant.","section":"Results, Figs. 4 and 8"},{"comment":"Equation (39) reads h(2) = R1ξ1, which appears to be a typo for h(2) = R^(2)ξ^(2); please correct this to avoid confusion with h(1) in Eq. (37).","section":"Appendix B, Eq. (39)"},{"comment":"The block notation in Eq. (31) is ambiguous: the θC,θC entry appears to combine the operator −∇·(κ0θC^2∇)/cC with a multiplication by θC, and the □ convention for operator action is not defined precisely; please specify exactly how each block acts on the entropy-gradient input fields.","section":"Appendix A, Eq. (31)"},{"comment":"There are minor grammatical errors (e.g., 'phenomena arises' in the Abstract) that should be corrected during copyediting.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The sign and factor inconsistency between Eq. (3) and K^(2) is straightforward to verify by direct substitution and is blocking the paper's main claim. The paper should be revised rather than rejected: the authors can correct either Eq. (3) or K^(2), but they must then check whether the simulation code and the reported results (Figs. 4, 6, 8) use the corrected deterministic dynamics, re-running simulations as needed. I would also encourage the authors to state explicitly whether the simulation code is available, as this would aid verification, but the scientific content can be assessed without it."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things you should know. First, this is a genuine extension of the authors' earlier SELM work: it couples a single protein's drift-diffusion to a fluctuating concentration field, a membrane temperature field, and an interfacial lipid temperature, with fluctuations derived from a GENERIC dissipative structure. The numerical validation is careful – second-order convergence for the heat discretization and a covariance check hitting its target to 6.8e-9. That part is solid.\n\nSecond, the central self-consistency claim does not survive contact with the appendix. I ran the substitution in the stress-test note and it holds. Using DS from Eq. 15 in K^(2) from Eq. 31 gives a θC drift where the c0∇Φ:κ̄∇q term has the opposite sign from Eq. 3, and the mechanical heating term is c0 q|∇Φ|²/(γ cC) instead of c0²|∇Φ|²/(γ cC). Since the stochastic fields are defined through B^(j)B^(j)^T = 2k_B K^(j), this breaks the fluctuation-dissipation balance. Either Eq. 3 or Eq. 31 is wrong; the paper does not bridge the gap. The covariance test in Fig. 10 only checks that the integrator samples the covariance defined by K^(j); it does not validate against the physical drift in Eqs. 1–4. The convergence test validates a heat equation, not the coupling terms.\n\nThe application studies are demonstrations rather than tests: thermal sensing depends on an assumed exponential filter, and the escape-time plots have no error bars. No code or data are provided.\n\nSo the paper is novel and the numerics are careful, but as printed the load-bearing mapping is inconsistent. I would not cite it in its current form. It does deserve a serious referee – the fix may be straightforward and the framework is useful – but the referee should demand a corrected derivation and a check that the GENERIC drift matches Eqs. 1–4.","headline":"A novel SELM extension to thermal and concentration gradients, but the key GENERIC operator in Appendix A contradicts the stated membrane-temperature equation, so the fluctuation-dissipation claim does not hold as printed.","tokens_in":20516,"tokens_out":10862,"would_cite":false,"duration_ms":100114,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82C31","65C30","92C05"],"pacs":["87.16.D-","05.40.Jc"],"model":"deepseek-v4-flash","headline":"A hybrid stochastic model whose noises are derived from the dissipative operators themselves gives self-consistent protein drift-diffusion coupled to fluctuating concentration and temperature fields.","keywords":["non-equilibrium statistical mechanics","membrane proteins","drift-diffusion","temperature gradients","concentration gradients","fluctuation-dissipation balance","hot Brownian motion","stochastic hybrid methods"],"falsifier":"Use the paper's own operator recipe: combine the three dissipative operators $K^{(1)}$, $K^{(2)}$, $K^{(3)}$ to compute the deterministic dynamics they imply, and check term by term whether those dynamics match the governing equations; if any term of the concentration or membrane-temperature equation fails to reappear, the fluctuation amplitudes are inconsistent with the modeled dynamics. A complementary experimental falsifier would be to place a well-characterized Brownian particle in a steady, known temperature gradient, measure its position and local temperature increments, and compare their covariance with the prediction forced by $B^{(j)}B^{(j),T} = 2k_B K^{(j)}$.","tokens_in":19496,"feed_emoji":"🌡️","tokens_out":16021,"duration_ms":155565,"temperature":0.7,"pith_summary":"The paper sets out to build simulation methods that track a single membrane protein while the membrane around it changes: the concentration of a signaling species diffuses and fluctuates, the temperature field heats, cools, and fluctuates, and the protein exchanges heat with both. The central claim is that the coupled equations are self-consistent in a precise statistical-mechanical sense: the random forces, fluxes, and heat fluctuations are not appended by hand but derived from the same dissipative operators that produce the deterministic drift, diffusion, and heat exchange. That matters for biology because protein positioning, thermal-gradient sensing, and escape from membrane energy wells are governed by how a protein's Brownian motion couples to fluctuations in its environment. The paper shows the model reproduces qualitative expectations for these three phenomena, such as faster escape from heated wells, and argues the same framework can be used for other biological systems and soft materials with mechanical-thermal coupling.","feed_headline":"One hybrid stochastic model tracks protein, concentration, and heat","feed_subtitle":"Self-consistent equations keep fluctuation-dissipation balance and reproduce heat-driven escape and gradient sensing.","key_machinery":"The load-bearing object is a set of three dissipative operators, $K^{(1)}$, $K^{(2)}$, $K^{(3)}$, together with the fluctuation-dissipation relation $B^{(j)}B^{(j),T} = 2k_B K^{(j)}$ that derives the stochastic driving fields from them. $K^{(1)}$ encodes the protein's overdamped motion and its thermal exchange with the interface; $K^{(2)}$ encodes diffusion of the concentration field, thermal conduction in the membrane, and their coupling; $K^{(3)}$ encodes heat exchange between the protein, the interface, and the membrane. Because each noise term is built from the same operator that gives the corresponding deterministic term, the fluctuations and the dissipation stay in balance even when the fields are spatially varying. On the numerical side, the machinery is completed by a finite-volume discretization with the discrete gradient being the negative adjoint of the discrete divergence, and by a two-stage stochastic integrator that reuses the same Wiener increments in both stages so the noise-induced drift is handled consistently.","core_discovery":"On the paper's own terms, the discovery is a way of writing the protein-membrane system in which the statistics of every fluctuation is forced to match the physics of the corresponding dissipation. The model has five state components: the protein position $X$, the concentration field $q$, the protein temperature $\\theta_P$, the interfacial temperature $\\theta_I$, and the membrane temperature field $\\theta_C$. The protein moves by overdamped Langevin dynamics with a mobility $M_{XX}$; the concentration field obeys a fluctuating advection-diffusion equation driven by the protein's chemical potential; and the temperature variables satisfy heat-exchange equations whose deterministic terms include mechanical work done by the protein and chemical fluxes. The noise amplitudes in all of these are obtained from the relation $B^{(j)}B^{(j),T} = 2k_B K^{(j)}$, where the $K^{(j)}$ are the operators describing the irreversible exchanges, which is what makes the model self-consistent rather than a set of equations with ad hoc noises. The paper then demonstrates the scheme on concentration-gradient-driven positioning, thermal-gradient sensing under fluctuations, and hot Brownian motion in energy wells, reporting that stronger external heating reduces well-escape times and that the escape time becomes dominated by how quickly the particle heats up once heating is strong.","pith_inferences":["The covariance check reported for one parameter set could be turned into a quantitative test of the central claim by comparing the joint position-temperature covariance of a hot Brownian particle against experiments with independently measured mobilities and thermal conductivities.","The framework extends naturally to reaction-driven heating: making the heat source in the temperature equations depend on a local reaction rate would cover catalytic enzymes and self-thermophoretic particles without changing the operator-based derivation of the noises.","Because the paper notes that mean-square displacement is an ambiguous descriptor out of equilibrium, a useful follow-up would be to convert the first-passage escape times into a position-dependent renormalized diffusivity and test whether that diffusivity carries memory from the particle's heating history.","If the correspondence between the dissipative operators and the deterministic equations were checked automatically during model construction, the same recipe could be applied to other soft-matter systems, such as active fluids or curved membranes, with noises obtained from the relation $B^{(j)}B^{(j),T} = 2k_B K^{(j)}$ each time."],"forward_implications":["Protein positioning in a concentration gradient is controlled by the ratio of the signaling molecule's diffusion time scale to the protein's motion time scale: fast-diffusing signals migrate to the protein, slow ones let the protein move to the signal, and intermediate cases meet between the two starting points.","An array of thermal-sensing proteins can recover a spatial temperature gradient that is hidden by fluctuations, provided downstream reaction chemistry time-averages the local temperature signal; too much filtering suppresses the resolvable gradient while too little leaves noise, so sensing involves a real trade-off.","Localized heating in an energy well reduces protein escape times, and when heating is strong the escape kinetics become controlled by how quickly the particle heats to the point where $k_B T$ is comparable to the barrier height, with only weak dependence on well depth.","Because the protein, concentration, and temperature fields evolve together, the method can capture effects that reduced single-particle descriptions miss, such as a hot protein locally heating membrane regions it has visited and thereby altering its own later escape kinetics.","The same construction is intended to transfer to other biological systems and soft materials where particle drift-diffusion couples to energy and mass exchange, not just to the membrane-protein examples used for demonstration."],"supporting_citations":[{"why":"Supplies the non-equilibrium statistical mechanics formalism whose dissipative operators and fluctuation relation are used to derive all noise terms.","marker":"[61, 62]"},{"why":"Provides the hybrid stochastic continuum-particle methodology for coupling discrete particles to fluctuating continuum fields.","marker":"[48, 49]"},{"why":"Gives the stochastic immersed boundary and finite-volume discretization conventions used for the spatial fields and particle-field forces.","marker":"[43, 44]"},{"why":"Grounds the overdamped mobility description of the protein and the fluctuation-dissipation link for its thermal noise.","marker":"[39, 40, 58]"},{"why":"Defines hot Brownian motion theory and experiments whose thermal-gradient phenomena the model extends by evolving particle and field temperatures.","marker":"[51, 52]"},{"why":"Supplies the reaction-diffusion reductions and stochastic finite-volume methods used for concentration coupling and the signaling filter.","marker":"[60]"},{"why":"Provides the Euler-Heun stochastic integrator whose two-stage update with shared Wiener increments is used to time-step the system.","marker":"[65]"}],"fun_headline_variants":["Self-consistent model merges protein, heat, and concentration noise","Hybrid model keeps fluctuation-dissipation balance across protein, heat, and density","New stochastic framework couples protein drift, thermal noise, and gradients","Model ties protein motion to heat and concentration fluctuations rigorously","Unified stochastic scheme for protein drift, heat flow, and gradient sensing"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire derivation rests on the claim that the three dissipative operators written in the appendix are exactly the irreversible parts of the governing equations; the paper asserts this correspondence rather than demonstrating it step by step, so if any operator is mis-specified the noise it generates would not match the dynamics it is meant to model.","fun_headline_variants_meta":{"raw":{"variants":["Self-consistent model merges protein, heat, and concentration noise","Hybrid model keeps fluctuation-dissipation balance across protein, heat, and density","New stochastic framework couples protein drift, thermal noise, and gradients","Model ties protein motion to heat and concentration fluctuations rigorously","Unified stochastic scheme for protein drift, heat flow, and gradient sensing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000544,"raw_usage":{"total_tokens":2601,"prompt_tokens":941,"completion_tokens":1660,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":1569}},"tokens_in":557,"tokens_out":1660,"duration_ms":13362,"temperature":1.0,"reasoning_tokens":1569,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:59:59.403888+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Use the paper's own operator recipe: combine the three dissipative operators $K^{(1)}$, $K^{(2)}$, $K^{(3)}$ to compute the deterministic dynamics they imply, and check term by term whether those dynamics match the governing equations; if any term of the concentration or membrane-temperature equation fails to reappear, the fluctuation amplitudes are inconsistent with the modeled dynamics. A complementary experimental falsifier would be to place a well-characterized Brownian particle in a steady, known temperature gradient, measure its position and local temperature increments, and compare their covariance with the prediction forced by $B^{(j)}B^{(j),T} = 2k_B K^{(j)}$.","supporting_citations":[{"cited_title":"Spatially Adaptive Stochastic Numerical Methods for Intrinsic Fluctuations in Reaction-Diffusion Systems","cited_arxiv_id":null,"evidence_quote":"Supplies the reaction-diffusion reductions and stochastic finite-volume methods used for concentration coupling and the signaling filter."},{"cited_title":"Numerical solution of stochastic differential equations","cited_arxiv_id":null,"evidence_quote":"Provides the Euler-Heun stochastic integrator whose two-stage update with shared Wiener increments is used to time-step the system."}],"review_version":1}