{"id":"37b6cedf-6339-4e71-a912-f00139f71051","arxiv_id":"2502.04578","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"ViPS imaging combines vibrational photothermal heating with optical diffraction tomography to visualize label-free thermophoretic transport inside living cells and measure intracellular diffusion and Soret coefficients.","lead":"This paper demonstrates a microscope method that uses infrared light to heat a small spot inside living cells and watches how biomolecules migrate along the resulting temperature gradient. It is the first label-free way to image intracellular thermophoresis and to measure how fast those molecules diffuse.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation 7 excludes the pre-existing dry-mass diffusion term ∇²σ; if that force-balance assumption fails, the fitted D and S_T absorb baseline intracellular inhomogeneity and the central quantitative claim is not robust.","rationale":"The reader identified exactly the same load-bearing assumption in Eq. 7, and I agree it is the most serious weakness. The stated goal of the paper is to quantify intracellular thermophoresis, and the numerical model is the only link between the measured Δσ images and the reported D and S_T values. The explicit exclusion of ∇²σ is not a mathematical error, but it is an unvalidated physical assumption about the intracellular force balance. Since the manuscript shows substantial dry-mass heterogeneity and attributes thermophoresis to mobile molecules, the condition that pre-existing gradients are balanced by non-thermal forces is not safe without a direct check. I considered other potential concerns: the 1064-nm control addresses the optical-tweezer alternative, the 5-ms thermal steady-state claim is supported by the authors' prior work, and the use of a single representative cell for the headline fit is a limitation but not an internal inconsistency. The ∇²σ omission is the one place where the central numerical result could be systematically biased rather than merely noisy. The reader's conditional verdict already reflects this weak spot, so my recommendation is UNCHANGED rather than a further downgrade: the paper should be published only if the model assumption is tested or the quantitative claim is softened.","tokens_in":11401,"tokens_out":3953,"duration_ms":46216,"concrete_test":"Re-run the fitting procedure for the data in Fig. 2e with the simulation equation modified to include the omitted diffusion source: ∂Δσ/∂t = D∇²(Δσ) + D∇²σ(x) + D S_T ∇·[∇T_s(x)(σ(x) + Δσ)], where σ(x) is the measured pre-heating dry-mass map from ODT. Refit D and S_T; if the fitted values move by more than the reported uncertainty (0.10 μm²/s for D, 0.10 × 10⁻³ K⁻¹ for S_T), the reported coefficients are not robust. As a decisive synthetic control, generate time-series data with known D and S_T from a realistic σ(x) with nonzero ∇²σ, then fit using both the original Eq. 7 and the corrected equation; quantify the bias introduced by the omission.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the quantitative isolation of thermophoretic transport from the measured time-resolved refractive-index change, yielding D ≈ 4.5 μm²/s and S_T ≈ 9.3 × 10⁻³ K⁻¹. The linchpin of that extraction is the simulation model, Eq. 7, whose diffusion term acts only on the perturbation Δσ: ∂Δσ/∂t = D∇²(Δσ) + D S_T ∇·[∇T_s(σ + Δσ)]. The manuscript explicitly excludes the diffusion of the pre-existing steady-state dry-mass distribution, i.e., D∇²σ(x), on the assumption that diffusive forces from that distribution are balanced by forces other than the temperature gradient. That assumption is load-bearing because cells are strongly heterogeneous: the same paper shows RI maps with substantial intranuclear contrast, and Fig. 4 shows that heterogeneity grows over hours. If σ(x) has non-negligible second derivatives, then even with no heating the mobile population would undergo diffusive redistribution; the difference between the heated and unheated states would then include a term of order D∇²σ t. Over 2.4 s, D∇²σ t ≈ 11 μm² · ∇²σ; using the measured RI-derived σ map, this can easily be comparable to the observed Δσ signal, especially at the edge of the heating spot where σ gradients are steep. If the omitted term is comparable to the fitted signal, the reported D and S_T are not measured values but effective parameters that absorb the cell's baseline inhomogeneity. The confidence intervals (±0.10 μm²/s, ±0.10 × 10⁻³ K⁻¹) are residual-based fit uncertainties and do not include this model-error contribution. The paper provides no independent evidence that the baseline σ gradients are force-balanced, nor does it show the magnitude of ∇²σ in the analyzed nucleus. This is a correctness risk for the central numbers, distinct from the well-controlled 1064-nm tweezing control and the plausible thermal-expansion separation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript introduces ViPS (vibrational photothermal Soret) microscopy, which combines continuous mid-IR heating with time-resolved optical diffraction tomography (ODT) to separate the fast thermal-expansion refractive-index change (Δn_ρ) from the slower dry-mass-concentration change (Δn_σ) driven by thermophoresis. The authors report intracellular diffusion and Soret coefficients in the nucleus of a living COS7 cell, D ≈ 4.6 μm² s⁻¹ and S_T ≈ 0.01 K⁻¹, obtained by fitting a 2D advection-diffusion model to the measured Δσ(t) curve. They also observe a reversed transport direction in the cytoplasm, which they attribute to diffusiophoresis, and a time-lapse reduction of thermophoretic activity in dying cells, which they interpret as glass formation. The central quantitative claims rest on the model in Eq. 7 and on a single-cell fit.","tokens_in":11716,"tokens_out":7446,"duration_ms":79554,"significance":"If the quantitative extraction is valid, ViPS imaging would be a significant methodological advance: it is label-free, provides simultaneous temperature and concentration maps, and extends thermophoresis measurements from in vitro systems to living cells. The time-scale separation is physically motivated, the 1064-nm off-resonance control addresses an important artifact, and the RNA-depletion experiment is a reasonable first step toward molecular attribution. However, the two headline coefficients come from a single representative cell, and the governing equation explicitly excludes a term that may be comparable to the signal; these issues must be resolved before the quantitative claims can be accepted.","major_comments":[{"comment":"The model excludes the diffusion term D∇²σ(x) of the pre-existing dry-mass distribution, with the stated assumption that the forces maintaining this distribution are balanced by factors other than the temperature gradient. This assumption is load-bearing: the measured Δσ signal is at the percent level, and the RI maps show substantial intranuclear inhomogeneity, so D∇²σ(x) over 2.4 s could easily contribute a signal comparable to the observed Δσ at the edge of the heating spot. The manuscript provides no direct evidence that σ is stationary over the measurement window in the absence of heating, nor a robustness test of the fitted D and S_T when the term is included or when a residual force term is added. Please provide such a control (e.g., a no-heating time series of σ from the same cells) or a sensitivity analysis of the fit; otherwise the reported D and S_T are effective parameters that may absorb baseline cell inhomogeneity.","section":"Methods, Eq. 7"},{"comment":"The reported D = 4.49 ± 0.10 μm²/s and S_T = 9.28 ± 0.10 × 10⁻³ K⁻¹ are obtained from a single representative cell's Δσ(t) curve, and the quoted uncertainties are residual-based confidence intervals, not cell-to-cell variability. The abstract presents these values as the method's quantitative output and uses them to infer a molecular weight around 100 kDa. The later 12-cell measurement (Fig. 3c) uses a different, more approximate estimator (−Δσ/σ/ΔT_s) and does not provide a complementary D distribution. Please report fitted D and S_T values across a population of cells, or explicitly state that the values are a single-cell demonstration and provide population-based estimates using a consistent analysis pipeline.","section":"Results, Fig. 2e"},{"comment":"The simulation is two-dimensional, but the measured quantity is a depth-integrated projection of a three-dimensional transport process, and the driving temperature field is the depth-averaged ΔT_s from a 4.2-μm-thick slab. The no-flux boundary condition (Eq. 8) constrains only the lateral nuclear boundary; axial thermophoretic and diffusive fluxes through the top and bottom of the integration volume are not modeled or bounded. This collapse from 3D to 2D may introduce systematic errors in D and S_T that are not reflected in the fit residuals. The authors should quantify this error, for example by performing a 3D simulation on the same cell geometry or by showing that the axial temperature gradient is negligible over the integration thickness.","section":"Eq. 7 and Fig. 2b"}],"minor_comments":[{"comment":"The statement that S_T ≈ 0.01 K⁻¹ means '~1% of the molecules migrate after 2.4 seconds of heating with a 1 K temperature difference' is imprecise; S_T relates the steady-state concentration gradient to the temperature gradient, not a finite-time fractional migration. Consider rephrasing to 'a steady-state relative concentration change of ~1% per kelvin'.","section":"Results"},{"comment":"The diffusion coefficient is compared with FRAP measurements that are typically performed at 37°C, while the present measurements were performed at 24°C. A brief comment on the expected temperature dependence of D would help the reader assess the molecular-weight inference.","section":"Results"},{"comment":"The cooling duration of 10.2 s is stated to ensure that the cell returns to its initial state, but no return-to-baseline data are shown. A supplementary figure of the temporal recovery of σ after heating would support the repeatability of the measurement cycle.","section":"Supplementary Note 1"},{"comment":"The interpretation of the time-lapse data as 'glass formation during the dying process' is presented as a conclusion from four cells and indirect RI evidence. Consider softening the claim or adding a quantitative metric of mobility change (e.g., a time-dependent apparent diffusion coefficient) to strengthen the connection.","section":"Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The core idea is promising and the experimental platform is impressive, but the two quantitative headline values rest on a single-cell fit and on a model with an explicitly excluded term that could be of the same order as the signal. These issues are addressable within the manuscript's scope through additional controls and population statistics, which is why I recommend major revision rather than rejection. I also note that reference 12 is the authors' own unpublished preprint; the editor may wish to verify that this prior work is publicly available and that the current paper's novelty is clearly delineated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, it is a real technical advance: combining vibrational photothermal microscopy with optical diffraction tomography to image label-free Soret transport inside living cells, while deriving the temperature map from the fast thermal-refractive-index component. That temperature map is independent of the transport model, and the 1064-nm control cleanly rules out optical tweezing as the cause of the slow concentration changes. Second, the headline numbers (D ≈ 4.6 μm²/s, S_T ≈ 9.3×10⁻³ K⁻¹) are less solid than they look. They come from fitting a 2D advection–diffusion model (Eq. 7) to what appears to be a single representative cell's Δσ curve, and the model explicitly excludes ∇²σ, the diffusion of the pre-existing dry-mass distribution, on a force-balance assumption. The authors state this assumption openly, but they give no evidence that the baseline dry-mass gradients are actually balanced by non-thermal forces. If the assumption fails, the fitted D and S_T absorb the cell's baseline inhomogeneity rather than pure thermophoresis. The reader's stress-test note quantifies this risk: over 2.4 s, D∇²σ t can be comparable to the observed signal near steep RI gradients. That is a legitimate model-error concern, not a refutation—it is testable by showing the magnitude of ∇²σ in the analyzed region or by fitting multiple cells and checking stability.\n\nWhat else is good: the RNA-reduction experiment (Actinomycin D) is a nice perturbation, giving a chemically specific handle on what is migrating. The observed opposite direction in cytoplasm versus nucleus is qualitatively compelling, and the time-lapse decline in thermophoresis over hours is suggestive, though only four cells and no significance testing. The nucleus–cytoplasm Soret values in Fig. 3c are single numbers without error bars, and the fitted D and S_T come without cell-to-cell statistics. No code or data are shipped; \"available upon request\" is weaker than it should be.\n\nOverall, the core physical separation of timescales is sound, the method is new, and the phenomena are interesting. The quantitative fragility is real but addressable, and the paper's own discussion shows the authors know some limits. This deserves peer review; a good referee should push for validation of the force-balance assumption, multi-cell statistics, and a look at the baseline ∇²σ magnitude. I would cite it if I worked on thermophoresis or label-free intracellular imaging, and it is worth a reading-group slot.","headline":"A genuinely new label-free way to image intracellular thermophoresis with simultaneous temperature mapping; the headline D and S_T values rest on a model assumption the authors state but do not validate.","tokens_in":12307,"tokens_out":2826,"would_cite":true,"duration_ms":31032,"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":"Time-resolved vibrational photothermal optical diffraction tomography separates thermal-expansion refractive index changes from concentration-driven changes, enabling label-free measurement of thermophoretic transport of biomolecules…","keywords":["thermophoresis","Soret effect","label-free imaging","optical diffraction tomography","vibrational photothermal microscopy","living cells","diffusiophoresis","dry mass transport"],"falsifier":"Compare ViPS fits on the same nucleus before and after disrupting the cytoskeleton or otherwise altering the steady-state dry-mass distribution while keeping the temperature map unchanged; if the recovered D and S_T shift systematically, the excluded ∇²σ term in Eq. 7 is not actually balanced. Alternatively, run the Eq. 7 simulation with and without the steady-state dry-mass gradient term on measured σ images and test whether the parameter search still returns the same D and S_T.","tokens_in":11152,"feed_emoji":"🔬","tokens_out":4692,"duration_ms":46262,"temperature":0.7,"pith_summary":"The paper introduces ViPS imaging, which uses a focused infrared beam to heat water inside a living cell and an optical diffraction tomography microscope to watch how the refractive index changes over time. By separating the fast refractive index change from thermal expansion from the slower change caused by molecules physically moving along the temperature gradient, the method extracts the thermophoretic (Soret) transport directly. In living COS7 cells the authors measure a nuclear diffusion coefficient of about 4.6 μm² s⁻¹ and a Soret coefficient of about 0.01 K⁻¹, implying a 1% dry-mass redistribution per kelvin. They also report reversed transport in the cytoplasm, which they attribute to thermophoresis-induced diffusiophoresis, and a loss of thermophoretic activity in dying cells under CO₂ depletion, consistent with glass-like aggregation.","feed_headline":"Laser heating exposes how molecules drift inside living cells","feed_subtitle":"One label-free scan yields the diffusion and Soret coefficients of nuclear biomolecules in living COS7 cells.","key_machinery":"The central object is the decomposition of the photothermal refractive-index change Δn into a fast thermal-expansion part Δn_ρ and a slow dry-mass-concentration part Δn_σ (Eq. 1), enabled by time-resolved ODT with 20 ms resolution. The steady-state temperature map ΔT_s is recovered from the instantaneous Δn_ρ image via the thermo-optic coefficient (Eq. 3), and the Δn_σ images are integrated over depth to give dry-mass concentration change Δσ (Eq. 2). These data are then fitted with a two-dimensional diffusion-drift simulation (Eq. 7) in which the Soret coefficient S_T and diffusion coefficient D are the free parameters.","core_discovery":"This paper claims that a single time-resolved refractive-index dataset from a vibrational photothermal optical diffraction tomography microscope can separate the thermal-expansion contribution from the concentration contribution to the photothermal signal, because the two evolve on very different time scales (microseconds vs. milliseconds to seconds). Exploiting this separation allows the temperature-rise map to be recovered from the fast component and the dry-mass transport to be followed in the slow component, so that both the temperature gradient driving thermophoresis and the molecular response to it are measured in the same label-free experiment. Using this approach the authors determine diffusion and Soret coefficients in the nucleus of living COS7 cells (D ≈ 4.6 μm² s⁻¹, S_T ≈ 0.01 K⁻¹), observe a negative Soret effect in the cytoplasm consistent with diffusiophoresis, and document a marked decline of thermophoretic mobility in dying cells.","pith_inferences":["Because the method records all molecular species simultaneously rather than a single fluorescent label, it could in principle be extended to survey thermophoretic mobility across the whole proteome, provided the diffusion-and-Soret fit is generalized to multi-component mixtures.","If the approximation in Eq. 7 is relaxed to include the pre-existing dry-mass gradient, the same dataset could yield a map of intracellular forces that balance baseline diffusion, turning the current assumption into a measurable quantity.","The reported disappearance of thermophoretic mobility in dying cells could be developed into a label-free viability assay, with the time constant of mobility loss potentially predicting the onset of irreversible aggregation.","Molecular-species-specific ViPS, using bond-selective mid-infrared excitation instead of the water overtone band, could separate the thermophoretic contributions of proteins, lipids, and nucleic acids and directly test the diffusiophoresis interpretation."],"forward_implications":["Intracellular thermophoresis can be quantified label-free in a single measurement, giving both the driving temperature field and the molecular response.","A 1 K intracellular temperature gradient redistributes roughly 1% of nuclear dry mass within seconds for molecules near 100 kDa.","RNA contributes an appreciable fraction of the observed nuclear thermophoretic signal, based on the Actinomycin D depletion experiment.","The reversal of transport direction in the cytoplasm indicates that diffusiophoresis can dominate thermophoresis in crowded intracellular environments.","Time-lapse ViPS imaging can track the loss of intracellular fluidity during cell death, offering a direct assay of glass-like transitions."],"supporting_citations":[{"why":"Supplies the ViP-ODT system used here and establishes the microsecond time scale of thermal diffusion that lets the fast and slow refractive-index components be separated.","marker":"[12]"},{"why":"Provides the theoretical basis for why molecules move along a temperature gradient, grounding the interpretation of the observed transport as thermophoresis.","marker":"[13]"},{"why":"Gives the refractive-index increment α used to convert Δn_σ images into depth-integrated dry-mass concentration changes Δσ.","marker":"[18]"},{"why":"Establishes the baseline direction of biomolecular thermophoresis toward lower temperature in aqueous environments, which the nuclear measurements reproduce.","marker":"[20]"},{"why":"Provides intracellular diffusion coefficients for GFP and IgG that are used to estimate the molecular weight of the thermophoretic species from the measured D.","marker":"[21]"},{"why":"Describes the thermophoresis-induced diffusiophoresis mechanism that the authors invoke to explain the negative Soret effect in the cytoplasm.","marker":"[26]"},{"why":"Supports the interpretation that the loss of thermophoretic mobility in dying cells reflects glass-forming behavior of the cytoplasm due to polymer aggregation.","marker":"[27]"},{"why":"The only prior study of intracellular thermophoresis, using fluorescence, which the present label-free approach extends by adding temperature-gradient measurement and species-blind mass transport.","marker":"[28]"}],"fun_headline_variants":["Soret coefficients mapped in live cells with photothermal microscopy","Label-free Soret imaging tracks molecular drift in live cells","Optical heating reveals intracellular thermophoresis without labels","Photothermal Soret microscopy reads molecular transport in living cells"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The simulation assumes that the cell's pre-existing uneven dry-mass distribution is already balanced by forces unrelated to the temperature gradient, so any baseline concentration gradient does not relax on its own during the measurement; if that balance does not hold, the fitted diffusion and Soret coefficients would absorb the cell's baseline inhomogeneity instead of pure thermophoresis.","fun_headline_variants_meta":{"raw":{"variants":["Soret coefficients mapped in live cells with photothermal microscopy","Label-free Soret imaging tracks molecular drift in live cells","Optical heating reveals intracellular thermophoresis without labels","Photothermal Soret microscopy reads molecular transport in living cells"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000225,"raw_usage":{"total_tokens":1472,"prompt_tokens":962,"completion_tokens":510,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":578,"completion_tokens_details":{"reasoning_tokens":444}},"tokens_in":578,"tokens_out":510,"duration_ms":6054,"temperature":1.0,"reasoning_tokens":444,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T22:14:18.081794+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare ViPS fits on the same nucleus before and after disrupting the cytoskeleton or otherwise altering the steady-state dry-mass distribution while keeping the temperature map unchanged; if the recovered D and S_T shift systematically, the excluded ∇²σ term in Eq. 7 is not actually balanced. Alternatively, run the Eq. 7 simulation with and without the steady-state dry-mass gradient term on measured σ images and test whether the parameter search still returns the same D and S_T.","supporting_citations":[{"cited_title":"Label-free mid-infrared photothermal microscopy revisits intracellular thermal dynamics: what do fluorescent nanothermometers measure?","cited_arxiv_id":"2406.16265","evidence_quote":"Supplies the ViP-ODT system used here and establishes the microsecond time scale of thermal diffusion that lets the fast and slow refractive-index components be separated."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the theoretical basis for why molecules move along a temperature gradient, grounding the interpretation of the observed transport as thermophoresis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the refractive-index increment α used to convert Δn_σ images into depth-integrated dry-mass concentration changes Δσ."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the baseline direction of biomolecular thermophoresis toward lower temperature in aqueous environments, which the nuclear measurements reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Describes the thermophoresis-induced diffusiophoresis mechanism that the authors invoke to explain the negative Soret effect in the cytoplasm."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supports the interpretation that the loss of thermophoretic mobility in dying cells reflects glass-forming behavior of the cytoplasm due to polymer aggregation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The only prior study of intracellular thermophoresis, using fluorescence, which the present label-free approach extends by adding temperature-gradient measurement and species-blind mass transport."}],"review_version":1}