{"id":"7f86a298-7436-4657-be02-6d7d9a373cb6","arxiv_id":"2608.04920","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Equilibrium molecular fluctuations in a confined charged fluid are converted into a space-time-resolved response matrix that reveals distinct, nonlocal pathways for momentum, charge, solute, and heat transport.","lead":"Using equilibrium molecular dynamics simulations, the authors extract a space-time-resolved matrix of coupled transport responses for a charged fluid confined between two walls. The framework shows that momentum and charge flow through distinct nonlocal pathways, which matters for designing nanofluidic devices such as ion channels, batteries, and water filters.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Thermostat distortion of long-time correlations is the key unresolved risk; the unsteady-Stokes fit is supportive but cannot rescue the central claim if the hydrodynamic tail is an artifact of the Nosé–Hoover bath.","rationale":"The paper's mathematical core—the Green–Kubo relation Eq. (4)—is standard and sound for Hamiltonian dynamics in equilibrium. The novelty is the spatiotemporal resolution of coupled response kernels, obtained by not integrating the correlations. Formally, the argument is: (i) use linear response theory to legitimize the equilibrium correlation as a response kernel; (ii) compute it from thermostat-regulated MD; (iii) interpret the slowly decaying diagonal particle kernel via an unsteady Stokes fit. The weakest step is the identification of the simulated stationary distribution/dynamics with the equilibrium ensemble of the unperturbed Hamiltonian. A Nosé–Hoover thermostat with damping time 0.1 is aggressive: it acts on the dynamical time scale of molecular motion and can equilibrate the system too strongly, altering the very long-time correlations (hydrodynamic tails, slowly relaxing modes) that the paper claims to resolve. The authors do not perform a thermostat-dependence check, nor do they report energy drift in an NVE validation, and the code/data are not yet public. The unsteady-Stokes fit is interesting but does not independently confirm the intrinsic nature of the tail: it only shows that a two-parameter fit can describe the extracted kernel. Given the framing as a new method for extracting spatiotemporal transport laws, reproducing the headline long-time hydrodynamic mode under different thermostats or in NVE is essential before the claim can be accepted as robust. The reader's concern on this point is exactly the most load-bearing one, so my verdict stays CONDITIONAL: the paper should be accepted only if these checks are added and the code/data released. If the thermostat test were to show the long-time mode persists, the conditional concerns would be fully resolved and the verdict could become ACCEPT.","tokens_in":18982,"tokens_out":1747,"duration_ms":19778,"concrete_test":"Re-run the core diagonal and off-diagonal correlation functions with at least two alternatives: (a) velocity-rescaling thermostat at the same temperature with a much longer damping time (e.g., tau_T = 5–10) and (b) microcanonical (NVE) production after equilibration, using identical starting configurations and binning. If the long-time tail of C_nn(t), the shape of K_nn(z,t), and the charge-sector correlations change materially (e.g., the apparent hydrodynamic relaxation time shifts by more than the reported statistical error), the central claim needs to be rephrased as thermostat-dependent. A second check: compute the same kernels from NVE segments shorter than the typical Lyapunov time to control energy drift.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that momentum transport appears as a long-lived, nonlocal hydrodynamic mode resolved from equilibrium fluctuations. This requires the measured long-time current correlations to be genuine intrinsic dynamics of the confined fluid. The simulations use a Nosé–Hoover thermostat with a very short damping time of 0.1 (Materials and Methods, S4), which strongly couples the fluid to a heat bath on the same scale as the molecular relaxation. Standard linear-response theory (SI Eqs. S34–S41) assumes Hamiltonian dynamics with a fixed equilibrium distribution; a deterministic thermostat modifies the Liouville dynamics and can renormalize or suppress hydrodynamic tails and alter off-diagonal couplings. The paper gives no direct test of thermostat dependence, so the 'long-lived mode'—the empirical basis of the strongest claim—may be an artifact of the thermostat rather than intrinsic confined-fluid memory. The unsteady-Stokes fit (SI Eqs. S20–S28) shows the data are consistent with a chosen hydrodynamic model after fitting H_eff and nu, but consistency with a fitted model does not establish that the correlations are unperturbed by the integrator/thermostat; a distorted slow tail could still be fitted. The reader flagged this as a weak assumption, and I agree it is the single most load-bearing concern.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a space–time-resolved generalization of the Onsager response matrix for confined fluids, in which transport is described by kernels R_ab(r,r',t) extracted from equilibrium current correlations via Eq. (4). The authors simulate a charged Lennard-Jones fluid in a slit pore, compute local–global and local–local response projections, and report a striking separation: the particle channel exhibits a long-lived, nonlocal momentum response that they identify with unsteady Stokes hydrodynamics, while the charge channel relaxes rapidly and locally. Off-diagonal particle–charge kernels are used to examine electro-osmotic flow and streaming current pathways, and a full 4×4 matrix for particle, solute, heat, and charge transport is presented, together with an excess-flux transformation in the Supplementary Material. The SI contains the linear-response derivation, detailed simulation parameters, uncertainty statements, analytical reference models, and a candid Limitations section.","tokens_in":19281,"tokens_out":5985,"duration_ms":76054,"significance":"If the central results are reliable, the framework is a valuable step beyond the usual reduction of nanoscale transport to effective coefficients: it makes the full spatiotemporal structure of coupled transport accessible from equilibrium MD and clarifies where Onsager reciprocity lives after coarse-graining. The paper is commendably explicit in deriving Eq. (4) from Liouville dynamics, in describing the 256-replica sampling protocol, in separating measured kernels from fitted reference models, and in acknowledging representation dependence and statistical cost. The commitment to public deposition of simulation data and analysis scripts, if honored, is an additional strength. The significance is currently tempered by two load-bearing concerns: the possible distortion of long-time correlations by the Nosé–Hoover thermostat, and the large statistical uncertainty in the off-diagonal charge sector that nonetheless carries several of the paper's interpretive claims.","major_comments":[{"comment":"The identification of measured correlations with response kernels in Eq. (4) is derived under Hamiltonian dynamics (SI Eqs. (S34)–(S41)), but the production runs use a Nosé–Hoover thermostat with damping time 0.1, which can couple to molecular and possibly hydrodynamic time scales. The paper provides no test of thermostat dependence, so the 'long-lived hydrodynamic mode' that anchors the main claim could be an artifact of the thermostat rather than intrinsic confined-fluid memory. Please add an NVE comparison or a thermostat-damping variation (for example, damping times from 0.1 to 10) and show that the slow kernels and the integrated L_nn(t) are unchanged within uncertainty.","section":"SI Materials and Methods, MD simulations; SI Eqs. (S34)–(S41)"},{"comment":"The off-diagonal charge entries in Table S2 are statistically indistinguishable from zero: L_cn/V = (-1.8 ± 2.8)×10^-2 and L_ch/V = (-2.3 ± 2.6)×10^-1. Since the off-diagonal electrokinetic and thermo-electric 'distinct pathways' are central results, these uncertainties need to be propagated into the kernels and profiles; Figs. 2 and 3 currently show no error bars for the off-diagonal panels. At minimum, report uncertainties for those panels and restrict the pairwise 'distinct pathway' claims to channels where the signal exceeds the noise.","section":"Table S2; Figs. 2 and 3"},{"comment":"The interpretation of the particle response as a 'hydrodynamic mode' relies on a two-parameter no-slip unsteady Stokes fit (H_eff and nu). A good fit to a flexible reference model does not by itself establish that the observed long-time tail is hydrodynamic, because the fitted parameters could absorb errors from the thermostat or from the Irving–Kirkwood current convention. Please provide an out-of-sample test, such as predicting the local–local G_nn(z,z') from the same fitted parameters, and report fit residuals as a function of z and t.","section":"SI Analytical reference models, Eqs. (S20)–(S28); Fig. 1C"}],"minor_comments":[{"comment":"In Eq. (5), the chain R_ab(r,r',t)=R†_ba(r',r,-t)=R†_ba(r',r,t) conflates stationarity with time-reversal symmetry; for cross-correlations the two relations are not redundant, and the chain as written implies R_ba(r',r,t)=R_ba(r',r,-t), which is not generally true. Please state the symmetry more carefully.","section":"Eq. (5)"},{"comment":"The color scales in the matrix plots are not defined: the main text says plotted quantities are rescaled for visualization, but the reader cannot determine zero levels, sign conventions, or amplitude ratios without consulting many separate captions. Please add explicit color bars and state the normalization for each panel.","section":"Fig. 3 and SI Fig. S2"},{"comment":"The uncertainty statement says errors are ±2 SE over the ensemble samples, but it is not clear whether the 95% confidence intervals include both replica variance and the two lateral directions as independent samples, nor whether block averaging over time is used. Please clarify the estimator.","section":"SI Materials and Methods, Extraction of response kernels"},{"comment":"The caption gives units as [a][b]/(k_B T σ τ), but for charge entries involving e^2 the reader must infer e=1; please state explicitly that e=1 in reduced units.","section":"Table S2 caption"},{"comment":"The excess-flux transformation uses bulk densities ρ_bulk_s and ρ_bulk_h; please state explicitly whether these are evaluated at the pore center and discuss the sensitivity of the excess-basis conclusions to this choice.","section":"SI Eqs. (S46)–(S49)"}],"recommendation":"major_revision","confidential_remarks":"The thermostat sensitivity is the key blocking issue: if the authors can demonstrate NVE or thermostat-invariance of the slow kernels, and if the off-diagonal figures are accompanied by propagated uncertainties, the paper would be a strong contribution. The manuscript does not show circularity, because the kernels are measured directly and the hydrodynamic model is used only as a reference. One editorial caution: the visual strength of Fig. 3 may exceed the statistical significance of the off-diagonal channels, so the final version should make the uncertainty status of each panel explicit."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe thing to know: this paper takes a textbook result—the Kubo relation between equilibrium current correlations and response functions—and actually carries it out in field-resolved form for a confined charged fluid. That sounds incremental, but the resulting object, the full 4x4 space-time response kernel, is not something I've seen before, and the physical contrast between the long-lived nonlocal momentum channel and the short-lived localized charge channel is a clean, useful demonstration. The paper is honest and careful; the limitations section is unusually candid.\n\nWhat's genuinely good: the derivation is clean, the MD implementation is described in detail (256 replicas, nonuniform bins, Irving-Kirkwood currents), and the excess-flux transformation in the SI shows they understand the representation-dependence of their object. The fact that the kernels are direct equilibrium correlations, not fitted to the hydrodynamic model, deals with the circularity concern—the Stokes fit is only an interpretive overlay.\n\nSoft spots, in order of seriousness. First, the Nosé–Hoover thermostat with damping 0.1 is aggressive, and the paper's central claim about a long-lived hydrodynamic mode rests on correlation tails out to t > 1000. A thermostat that strongly couples to the fluid can renormalize or suppress those tails, and the paper gives no test of thermostat dependence. This is a real gap, not a nitpick. It doesn't sink the paper—the charge-localization result is less sensitive to this—but it means the headline 'hydrodynamic mode' needs a control. Second, the off-diagonal charge channel has enormous error bars (Table S2: L_cn/V = -1.8e-2 ± 2.8e-2), and Fig. 3 is shown without error bars. That's a presentation issue as much as a statistical one. Third, the heat current is convention-dependent; they acknowledge it, but it limits the quantitative force of the heat sector. Finally, data and code are promised but not yet public; for a paper this simulation-heavy, that should be a condition of acceptance.\n\nThe thermodynamic interpretation (the backflow asymmetry between EOF and SC) is suggestive but not fully nailed down; I'd trust the kernels more than the narrative.\n\nWho is this for? Nanofluidics and MD practitioners who want a more complete microscopic description than integrated Onsager coefficients. It's a proof of principle, not a breakthrough, but it's a useful one.\n\nMy recommendation: send it out for peer review, but require the thermostat test and error bars on the full matrix before acceptance. The framework deserves to be evaluated by the community.\n\nBest.","headline":"Textbook Kubo theory, but the first full 4x4 space-time response matrix for a confined charged fluid — worth refereeing, with a thermostat control and error bars required.","tokens_in":19739,"tokens_out":3235,"would_cite":true,"duration_ms":38695,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["05.60.-k","47.61.-k"],"model":"deepseek-v4-flash","headline":"Coupled transport in nanoconfined fluids is a space-time-resolved Onsager response kernel extracted from equilibrium molecular fluctuations, with momentum long-lived and nonlocal, charge fast and local, and the two linked through distinct…","keywords":["nanoconfined fluids","Onsager response matrix","space-time response kernels","molecular dynamics","electrokinetic coupling","Green-Kubo relations","nonlocal transport","electrical double layer"],"falsifier":"Run the same confined system under a weak, localized, time-dependent force (for instance a moving barrier or a small oscillating electric field localized near one wall) and compare the resulting flux at position z and time t with the convolution of the equilibrium kernel against the driving gradient; if the long-time tail predicted by the kernel does not match the nonequilibrium measurement, the identification of correlations with response fails. A simpler check is to re-extract $K_{nn}(z,t)$ in the microcanonical ensemble: the claimed long-lived hydrodynamic mode should persist without any thermostat.","tokens_in":18787,"feed_emoji":"🔬","tokens_out":6847,"duration_ms":78319,"temperature":0.7,"pith_summary":"Nanoconfined flow, solute, heat, and charge transport are usually compressed into bulk coefficients and effective interfacial parameters such as slip lengths, zeta potentials, and resistances. The paper argues this compression discards the essential physics: under nanoconfinement the response of a fluid is spatially heterogeneous and has memory, so the right object is a space–time-resolved Onsager response matrix $R_{ab}(\\mathbf r,\\mathbf r',t)$ extracted from equilibrium current correlations. Using molecular dynamics of a charged Lennard–Jones fluid in a slit pore, it shows that momentum transport is a long-lived, nonlocal hydrodynamic mode across the pore, charge transport relaxes rapidly through localized ionic friction, and the off-diagonal electrokinetic kernels connect these two behaviors. If correct, every familiar transport coefficient is a fully integrated limit of this kernel, and nonlocal, history-dependent constitutive laws become the natural description at the nanoscale. The payoff would be a direct microscopic route to effective interfacial parameters and a new target for coarse-grained and data-driven transport models.","feed_headline":"Nanopore transport is a space-time kernel, not a coefficient","feed_subtitle":"Equilibrium fluctuations reveal how momentum, heat, and charge couple nonlocally under confinement.","key_machinery":"The central object is the space–time-resolved Onsager response kernel $R_{ab}(\\mathbf r,\\mathbf r',t)$, the equilibrium correlation of local flux densities divided by $k_B T$, defined through Eq. (4). It carries the argument because all familiar transport quantities are projections of it: local–global kernels $K_{ab}(\\mathbf z,t)$, stationary profiles $M_{ab}(\\mathbf z)$, two-point responses $G_{ab}(\\mathbf z,\\mathbf z')$ (via Einstein–Helfand), global correlations $C_{ab}(t)$, and the Onsager matrix $L_{ab}$ as the full large-time integral. The paper deliberately stops before coarse graining and treats these correlations as the central observables, using the unsteady Stokes Green function only as a reference to interpret the hydrodynamic sector.","core_discovery":"The central claim, stated on its own terms, is that coupled transport in a nanoconfined fluid is fully described by the equilibrium correlation kernel $R_{ab}(\\mathbf r,\\mathbf r',t) = (k_B T)^{-1}\\langle \\mathbf j_a(\\mathbf r,t)\\,\\mathbf j_b^\\dagger(\\mathbf r',0)\\rangle_{\\rm eq}$, and that this kernel, rather than its space/time integrals, is the physically meaningful response function under confinement. The paper shows in a charged Lennard-Jones slit that particle ($n$), solute ($s$), heat ($h$), and charge ($c$) channels organize in a block structure: $n$, $s$, and $h$ share long-lived, pore-spanning hydrodynamic relaxation while $c$ relaxes fast and locally; the reciprocal electro-osmotic and streaming-current kernels have the same integrated Onsager coefficient but distinct spatial and temporal pathways, with a backflow-like negative region in the electro-osmotic profile. Because the kernel is the pre-integrated object, Onsager reciprocity holds nonlocally, $\\mathbf G_{ab}(\\mathbf z,\\mathbf z') = \\mathbf G_{ba}(\\mathbf z',\\mathbf z)$, while one-sided projections need not be symmetric. Conventional coefficients, including the full Onsager matrix, are recovered by the successive integrations of Eq. (3), and effective parameters such as slip length or zeta potential are reinterpreted as coarse-grained summaries of this underlying field-level response.","pith_inferences":["An extension the authors only gesture at: because Eq. (4) is a field-level Green–Kubo relation, the same protocol can yield frequency- and wavevector-dependent transport functions, turning the kernel into a microscopic impedance or memory function for nanofluidic devices.","The fitted unsteady-Stokes comparison suggests a testable prediction: the effective kinematic viscosity and confinement extracted from the time-dependent kernel should coincide with values from independent nonequilibrium or microcanonical measurements; disagreement would pinpoint thermostat artifacts.","The formalism implies that apparent violations of Onsager symmetry in projected or coarse-grained coefficients are projection artifacts; a direct check would be to construct coarse-grained coefficients from one-sided projections and verify that the full kernel still satisfies $G_{ab}(\\mathbf z,\\mathbf z')=G_{ba}(\\mathbf z',\\mathbf z)$.","A concrete extension to real materials would replace the minimal Lennard-Jones electrolyte with water-like models and check whether rotational and hydration modes add spectrally separated features to the heat and charge kernels, as the authors themselves note."],"forward_implications":["Conventional transport coefficients follow as fully integrated limits of the response kernels: the Onsager matrix $L_{ab}$, conductivity, permeability, and mobility profiles all emerge from successive spatial and temporal integration.","Nanoscale momentum transport is genuinely nonlocal and long-lived (hydrodynamic modes crossing the pore), so single-point local closures and steady-state effective parameters miss the physics that matters.","Charge transport is fast, localized ionic friction; its stationary profile is near-bulk in the pore center, so charge relaxation can be modeled locally once the interfacial structure is retained.","Electrokinetic coupling is interfacial but nonlocal: EOF and streaming current have different dynamical pathways, and the sign of the integrated electrokinetic response depends on EDL spatial organization and hydrodynamic weighting, not just net charge.","Reciprocity lives at the kernel level; projected one-sided quantities need not be symmetric, giving a microscopic basis for nonlocal, history-dependent transport laws."],"supporting_citations":[{"why":"Supplies the linear-response derivation (Eqs. S38–S41) and the molecular dynamics protocol behind the kernel extraction.","marker":"(30)"},{"why":"Provides the Kubo linear-response relation on which the kernel expression Eq. (4) and its time-reversal symmetry Eq. (5) rest.","marker":"(32)"},{"why":"Defines the Onsager matrix and reciprocal relations that the space-time kernel generalizes and recovers upon full integration.","marker":"(27)"},{"why":"Basis for the Green–Kubo and Kirkwood–Buff constructions that the paper shows are projections of the same field-level correlations.","marker":"(18,31)"},{"why":"Irving–Kirkwood construction used to define local density and current fields (Eqs. S1–S8).","marker":"(47)"},{"why":"Provides the molecular dynamics simulation engine used to generate the equilibrium trajectories.","marker":"(49)"},{"why":"Einstein–Helfand relation used to extract stationary local-local responses $G_{ab}(z,z')$ from the long-time growth rate of Helfand covariances.","marker":"(50)"},{"why":"Reference for diffuse-charge dynamics and ionic friction used to interpret the rapid localized charge response.","marker":"(34)"}],"fun_headline_variants":["Confined fluids: transport as a space-time kernel, not a coefficient","Molecular fluctuations map nonlocal transport beneath the nanopore","Space-time kernel replaces coefficients for nanoconfined transport","Coupled transport in confinement is a space-time memory kernel","Nonlocal, history-dependent transport from equilibrium noise"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole derivation assumes the equilibrium current correlations measured in the thermostatted simulations are the true response kernels, which requires that the Nosé–Hoover damping does not alter the long-time dynamics behind the claimed hydrodynamic mode and that the chosen Irving–Kirkwood definitions, especially for the heat current, capture the physically relevant fluxes.","fun_headline_variants_meta":{"raw":{"variants":["Confined fluids: transport as a space-time kernel, not a coefficient","Molecular fluctuations map nonlocal transport beneath the nanopore","Space-time kernel replaces coefficients for nanoconfined transport","Coupled transport in confinement is a space-time memory kernel","Nonlocal, history-dependent transport from equilibrium noise"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000379,"raw_usage":{"total_tokens":2035,"prompt_tokens":983,"completion_tokens":1052,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":599,"completion_tokens_details":{"reasoning_tokens":971}},"tokens_in":599,"tokens_out":1052,"duration_ms":11244,"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-06T13:41:32.481783+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same confined system under a weak, localized, time-dependent force (for instance a moving barrier or a small oscillating electric field localized near one wall) and compare the resulting flux at position z and time t with the convolution of the equilibrium kernel against the driving gradient; if the long-time tail predicted by the kernel does not match the nonequilibrium measurement, the identification of correlations with response fails. A simpler check is to re-extract $K_{nn}(z,t)$ in the microcanonical ensemble: the claimed long-lived hydrodynamic mode should persist without any thermostat.","supporting_citations":[],"review_version":1}