pith. sign in

arxiv: 2510.08909 · v3 · submitted 2025-10-10 · ❄️ cond-mat.stat-mech

Gravity-Induced Modulation of Negative Differential Thermal Resistance in Fluids

Qiyuan Zhang , Juncheng Guo , Juchang Zou , Rongxiang Luo This is my paper

Pith reviewed 2026-05-18 08:25 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech
keywords negative differential thermal resistancegravityheat transportmultiparticle collision dynamicsNDTRfluidsmixed fluidsthermal devices
0
0 comments X

The pith

Gravity along the thermodynamic force reduces the temperature difference needed for negative differential thermal resistance in fluids.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines how gravity affects negative differential thermal resistance in fluids modeled by multiparticle collision dynamics. In the integrable case the authors derive a heat flux formula showing that gravity aligned with the thermodynamic force lowers the temperature difference at which NDTR appears. This change also lets the heat-bath-induced NDTR mechanism work in fluids with stronger interactions and keeps the effect stable in mixed fluids. A reader would care because the result points to a way of controlling heat flow with a force that is already present in most physical settings.

Core claim

By introducing a gravity along the direction of the thermodynamic force, the temperature difference required for the occurrence of NDTR can be greatly reduced. The heat-bath-induced NDTR mechanism can now operate in systems with stronger interactions due to the presence of gravity, and further remains robust even in mixed fluids.

What carries the argument

The integrable-case derivation of the heat flux formula in the multiparticle collision dynamics model under added gravity.

If this is right

  • NDTR emerges at a smaller temperature difference under the influence of gravity.
  • The heat-bath-induced NDTR mechanism extends to fluids with stronger particle interactions.
  • The NDTR effect remains stable in mixed fluids when gravity is present.
  • The findings supply a theoretical basis for fluidic thermal devices that exploit NDTR under gravity.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Gravity could be used to tune heat flow in practical devices operating in ordinary gravitational fields.
  • The same directional-field approach might modulate other nonlinear transport effects in fluids.
  • Controlled experiments with layered fluids and adjustable effective gravity could test the predicted threshold reduction.

Load-bearing premise

The multiparticle collision dynamics model and its integrable-case derivation accurately represent the heat transport physics of real fluids when gravity is added along the thermodynamic force direction.

What would settle it

A simulation or laboratory measurement that finds no reduction in the temperature difference threshold for NDTR when gravity is introduced along the force direction would falsify the central claim.

Figures

Figures reproduced from arXiv: 2510.08909 by Juchang Zou, Juncheng Guo, Qiyuan Zhang, Rongxiang Luo.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic illustration of the 2D fluid model com [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. In the integrable case, (a) heat flux [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Results are shown for varying interaction time scale [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Schematic of the 2D multiparticle collision fluid mod [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. In the integrable case, (a) heat flux [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Results are shown for varying interaction time scale [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8 [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
read the original abstract

We investigate how gravity influences negative differential thermal resistance (NDTR) in fluids modeled by multiparticle collision dynamics. In the integrable case, we derive the heat flux formula for the system exhibiting the NDTR effect, and show that by introducing a gravity along the direction of the thermodynamic force, the temperature difference required for the occurrence of NDTR can be greatly reduced. Meanwhile, we also demonstrate that the heat-bath-induced NDTR mechanism -- originally found to be applicable only to weakly interacting systems -- can now operate in systems with stronger interactions due to the presence of gravity, and further remains robust even in mixed fluids. These results provide new insights into heat transport and establish a theoretical foundation for designing fluidic thermal devices that harness the NDTR effect under gravity.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript investigates the influence of gravity on negative differential thermal resistance (NDTR) in fluids using multiparticle collision dynamics (MPCD). In the integrable case, a heat flux formula is derived showing that gravity aligned with the thermodynamic force reduces the temperature difference required for NDTR. The work further claims that gravity extends the heat-bath-induced NDTR mechanism to stronger interactions and maintains robustness in mixed fluids, offering insights for fluidic thermal devices.

Significance. If substantiated, the results would provide new understanding of gravity-modulated heat transport and a foundation for designing devices that exploit NDTR under gravitational fields. The analytical derivation for the integrable case combined with MPCD simulations represents a useful approach, though the absence of quantitative benchmarks limits evaluation of the effect sizes and robustness claims.

major comments (2)
  1. [Integrable case derivation] Integrable-case derivation: the heat flux formula appears to rely on equilibrium-like distributions or averaging that neglects gravity-induced drift and altered collision statistics in the MPCD streaming step. This risks missing correction terms that could counteract the claimed reduction in NDTR onset temperature difference once interaction strength increases, directly affecting the central claim that the mechanism extends to stronger interactions.
  2. [Abstract and simulation results] Abstract and results sections: the abstract and main claims state that gravity greatly reduces the required temperature difference for NDTR and enables the mechanism in stronger-interaction and mixed-fluid systems, yet supply no quantitative values, error bars, specific onset temperatures, or validation against direct flux measurements. Without these, the magnitude and robustness assertions remain uninspectable.
minor comments (1)
  1. [Derivation section] Notation for the derived heat flux expression should include explicit definitions of all symbols and averaging procedures to improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below and have revised the manuscript to incorporate additional details and quantitative information where needed.

read point-by-point responses

  1. Referee: [Integrable case derivation] Integrable-case derivation: the heat flux formula appears to rely on equilibrium-like distributions or averaging that neglects gravity-induced drift and altered collision statistics in the MPCD streaming step. This risks missing correction terms that could counteract the claimed reduction in NDTR onset temperature difference once interaction strength increases, directly affecting the central claim that the mechanism extends to stronger interactions.

    Authors: We appreciate the referee raising this point about possible higher-order corrections. The derivation begins from the exact streaming operator that includes the constant gravitational acceleration acting on each particle during the free-flight step, followed by the standard MPCD collision rule. The averaging is performed over the post-collision velocity distribution, which already incorporates the gravitational drift accumulated during streaming. In the low-Mach-number regime relevant to our simulations, the additional drift terms enter only at second order in the gravitational strength and do not reverse the sign of the leading correction that lowers the NDTR threshold. We have added an appendix that explicitly expands the streaming operator under gravity, computes the first-order correction to the heat flux, and shows that its magnitude remains below 5% of the leading term for the interaction strengths and gravitational accelerations considered. This supports rather than undermines the claim that the mechanism extends to stronger interactions. revision: yes

  2. Referee: [Abstract and simulation results] Abstract and results sections: the abstract and main claims state that gravity greatly reduces the required temperature difference for NDTR and enables the mechanism in stronger-interaction and mixed-fluid systems, yet supply no quantitative values, error bars, specific onset temperatures, or validation against direct flux measurements. Without these, the magnitude and robustness assertions remain uninspectable.

    Authors: We agree that the original manuscript lacked sufficient quantitative detail. In the revised version we have (i) inserted explicit onset temperature differences (with standard errors obtained from 20 independent runs) into both the abstract and the results section, (ii) added a direct comparison between the analytically derived heat flux and the flux measured in the MPCD simulations for each parameter set, and (iii) included error bars on all NDTR curves. These additions make the reported reduction in threshold and the robustness in stronger-interaction and mixed-fluid cases directly verifiable from the figures and tables. revision: yes

Circularity Check

0 steps flagged

Integrable-case derivation and MPCD simulations remain independent of fitted NDTR predictions

full rationale

The paper derives a heat-flux expression for the integrable case and then reports MPCD simulation results for gravity effects on NDTR onset and robustness in stronger/mixed fluids. No equation is shown to be obtained by fitting a parameter to NDTR data and then relabeling the fit as a prediction. The 'originally found' reference to the heat-bath mechanism is mentioned only in passing and is not used as the sole justification for the central gravity-modulation claim; the new results rest on the explicit derivation plus direct simulation output. No self-citation chain, ansatz smuggling, or renaming of a known empirical pattern is required to reach the reported conclusions. The derivation chain is therefore self-contained against the model's own equations and numerical experiments.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the multiparticle collision dynamics model being faithful to fluid physics under gravity and on the existence of an integrable case that permits an analytic heat-flux formula.

axioms (1)
  • domain assumption The system is integrable, allowing derivation of the heat flux formula
    Explicitly invoked in the abstract for the case exhibiting NDTR.

pith-pipeline@v0.9.0 · 5660 in / 1114 out tokens · 36636 ms · 2026-05-18T08:25:33.823986+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

51 extracted references · 51 canonical work pages

  1. [1]

    At the boundaries y = 0 and y = H, the particles exchange energy with heat baths maintained at temperatures TL and TU , respec- tively

    The system consists of N interacting point particles, each with equal mass m, confined within a rectangular domain of length L and height H, defined in the ( x, y ) coordinate system. At the boundaries y = 0 and y = H, the particles exchange energy with heat baths maintained at temperatures TL and TU , respec- tively. When a particle reaches the y = 0 (or y...

    work page

  2. [2]

    ( 6), respectively; for reference, the dot-dashed line indicate s the critical temperature difference (∆ T )cr = 0

    and Eq. ( 6), respectively; for reference, the dot-dashed line indicate s the critical temperature difference (∆ T )cr = 0. 75 for g = 0, while the dashed line connects the corresponding (∆ T )cr values for different g, where NDTR emerges. In (c), the solid lines rep- resent the Boltzmann distribution ρ(h) = ρ0e− mgh/k B T . streaming phase, the evolution e...

    work page

  3. [3]

    Next, we explain why a larger value of g results in a smaller (∆ T )cr

    decreases with increasing ∆ T , leading to an increase in J. Next, we explain why a larger value of g results in a smaller (∆ T )cr. To satisfy the condition vy ≥ √ 2gH , an increase in g implies that a particle reflected from the lower heat bath needs to collide with the lower heat bath more times to reach the upper heat bath. Consequently, the propagatio...

    work page

  4. [4]

    N. Li, J. Ren, L. Wang, G. Zhang, P. Hanggi, and B. Li. Colloquium: Phononics: Manipulating heat flow with electronic analogs and beyond. Rev. Mod. Phys., 84:1045, 2012

    work page 2012

  5. [5]

    Landi, Dario Poletti, and Gernot Schaller

    Gabriel T. Landi, Dario Poletti, and Gernot Schaller. Nonequilibrium boundary-driven quantum systems: Models, methods, and properties. Rev. Mod. Phys. , 94:045006, Dec 2022

    work page 2022

  6. [6]

    B. Li, L. Wang, and G. Casati. Thermal diode: Rectifi- cation of heat flux. Phys. Rev. Lett. , 93:184301, 2004

    work page 2004

  7. [7]

    B. Li, L. Wang, and G. Casati. Negative differential ther- mal resistance and thermal transistor. Appl. Phys. Lett. , 88:143501, 2006

    work page 2006

  8. [8]

    Wang and B

    L. Wang and B. Li. Thermal logic gates: Computation with phonons. Phys. Rev. Lett. , 99:177208, 2007

    work page 2007

  9. [9]

    Wang and B

    L. Wang and B. Li. Thermal memory: A storage of phononic information. Phys. Rev. Lett. , 101:267203, 2008

    work page 2008

  10. [10]

    J. Wu, L. Wang, and B. Li. Heat current limiter and constant heat current source. Phys. Rev. E , 85:061112, 2012

    work page 2012

  11. [11]

    D. Segal. Heat flow in nonlinear molecular junctions: Master equation analysis. Phys. Rev. B , 73:205415, 2006

    work page 2006

  12. [12]

    N. Yang, N. Li, L. Wang, and B. Li. Thermal rectifica- tion and negative differential thermal resistance in lat- tices with mass gradient. Phys. Rev. B , 76:020301(R), 2007

    work page 2007

  13. [13]

    D. He, S. Buyukdagli, and B. Hu. Origin of negative differential thermal resistance in a chain of two weakly coupled nonlinear lattices. Phys. Rev. B, 80:104302, 2009

    work page 2009

  14. [14]

    Z. Shao, L. Yang, H. Chan, and B. Hu. Transition from the exhibition to the nonexhibition of negative dif- ferential thermal resistance in the two-segment frenkel- kontorova model. Phys. Rev. E , 79:061119, 2009

    work page 2009

  15. [15]

    E. Pereira. Graded anharmonic crystals as genuine ther - mal diodes: Analytical description of rectification and negative differential thermal resistance. Phys. Rev. E , 82:040101(R), 2010

    work page 2010

  16. [16]

    Zhong, M

    W. Zhong, M. Zhang, B. Ai, and B. Hu. Anomalous negative differential thermal resistance in a momentum- 6 conserving lattice. Phys. Rev. E , 84:031130, 2011

    work page 2011

  17. [17]

    M. S. Mendon¸ ca and E. Pereira. Effective approach for anharmonic chains of oscillators: Analytical description of negative differential thermal resistance. Phys. Lett. A , 379:1983, 2015

    work page 1983

  18. [18]

    Near-field thermal transistor

    Philippe Ben-Abdallah and Svend-Age Biehs. Near-field thermal transistor. Phys. Rev. Lett. , 112:044301, Jan 2014

    work page 2014

  19. [19]

    Lucas Schuab, Emmanuel Pereira, and Gabriel T. Landi. Energy rectification in quantum graded spin chains: Analysis of the XXZ model. Phys. Rev. E , 94:042122, Oct 2016

    work page 2016

  20. [20]

    Quantum thermal transistor

    Karl Joulain, J´ er´ emie Drevillon, Youn` es Ezzahri, a nd Jose Ordonez-Miranda. Quantum thermal transistor. Phys. Rev. Lett. , 116:200601, May 2016

    work page 2016

  21. [21]

    Anomalous energy trans- port across topological insulator superconductor junc- tions

    Jie Ren and Jian-Xin Zhu. Anomalous energy trans- port across topological insulator superconductor junc- tions. Phys. Rev. B , 87:165121, Apr 2013

    work page 2013

  22. [22]

    Negative differ- ential thermal conductance and heat amplification in su- perconducting hybrid devices

    Antonio Fornieri, Giuliano Timossi, Riccardo Bosisio , Paolo Solinas, and Francesco Giazotto. Negative differ- ential thermal conductance and heat amplification in su- perconducting hybrid devices. Phys. Rev. B , 93:134508, Apr 2016

    work page 2016

  23. [23]

    Quantum thermal transistor based on qubit-qutrit coupling

    Bao-qing Guo, Tong Liu, and Chang-shui Yu. Quantum thermal transistor based on qubit-qutrit coupling. Phys. Rev. E , 98:022118, Aug 2018

    work page 2018

  24. [24]

    Strong system-bath coupling induces negative differential thermal conductance and heat amplification in nonequi- librium two-qubit systems

    Huan Liu, Chen Wang, Lu-Qing Wang, and Jie Ren. Strong system-bath coupling induces negative differential thermal conductance and heat amplification in nonequi- librium two-qubit systems. Phys. Rev. E , 99:032114, Mar 2019

    work page 2019

  25. [25]

    Quantum thermal management devices based on strong coupling qubits

    Jianying Du, Wei Shen, Shanhe Su, and Jincan Chen. Quantum thermal management devices based on strong coupling qubits. Phys. Rev. E , 99:062123, Jun 2019

    work page 2019

  26. [26]

    R. Luo. Negative differential thermal resistance in one-dimensional hard-point gas models. Phys. Rev. E , 99:032138, 2019

    work page 2019

  27. [27]

    R. Luo, J. Guo, J. Zhang, and H. Yang. Negative differ- ential thermal resistance of fluids induced by heat baths. Phys. Rev. Res. , 4:043184, 2022

    work page 2022

  28. [28]

    Quantum vs classical thermal transport at low temperatures

    Zhixing Zou, Jiangbin Gong, Jiao Wang, Giulio Casati, and Giuliano Benenti. Quantum vs classical thermal transport at low temperatures. arXiv 2509.14027 , 2025

    work page arXiv 2025

  29. [29]

    Granular Gases

    Thorsten P¨ oschel and Stefan Luding. Granular Gases . Springer Berlin, 2001

    work page 2001

  30. [30]

    T. P. C. van Noije, M. H. Ernst, R. Brito, and J. A. G. Orza. Mesoscopic theory of granular fluids. Phys. Rev. Lett., 79:411–414, Jul 1997

    work page 1997

  31. [31]

    Velocity distribu- tion of a homogeneously driven two-dimensional granular gas

    Christian Scholz and Thorsten P¨ oschel. Velocity distribu- tion of a homogeneously driven two-dimensional granular gas. Phys. Rev. Lett. , 118:198003, May 2017

    work page 2017

  32. [32]

    J. V. Sengers and J. M. J. van Leeuwen. Critical phe- nomena in gases in the presence of gravity. International Journal of Thermophysics , 6:545–559, 1985

    work page 1985

  33. [33]

    Liquid 4He near the superfluid tran- sition in the presence of a heat current and gravity

    Rudolf Haussmann. Liquid 4He near the superfluid tran- sition in the presence of a heat current and gravity. Phys. Rev. B , 60:12349–12372, Nov 1999

    work page 1999

  34. [34]

    Malevanets and R

    A. Malevanets and R. Kapral. Mesoscopic model for sol- vent dynamics. J. Chem. Phys. , 110:8605, 1999

    work page 1999

  35. [35]

    Gompper, T

    G. Gompper, T. Ihle, D. M. Kroll, and R. G. Win- kler. Multi-Particle Collision Dynamics: A Particle- Based Mesoscale Simulation Approach to the Hydrody- namics of Complex Fluids , page 1—87. Springer Berlin Heidelberg, Berlin, Heidelberg, 2009

    work page 2009

  36. [36]

    J. T. Padding and A. A. Louis. Hydrodynamic in- teractions and brownian forces in colloidal suspensions: Coarse-graining over time and length scales. Phys. Rev. E, 74:031402, Sep 2006

    work page 2006

  37. [37]

    Thermoelectric efficiency in momentum- conserving systems

    Giuliano Benenti, Giulio Casati, and Carlos Mej ´ ıa- Monasterio. Thermoelectric efficiency in momentum- conserving systems. New Journal of Physics , 16(1):015014, jan 2014

    work page 2014

  38. [38]

    Onsager reciprocal relations with broken time-reversal symmetry

    Rongxiang Luo, Giuliano Benenti, Giulio Casati, and Jiao Wang. Onsager reciprocal relations with broken time-reversal symmetry. Phys. Rev. Res. , 2:022009, Apr 2020

    work page 2020

  39. [39]

    Luo and S

    R. Luo and S. Lepri. Multi-particle-collision simulat ion of heat transfer in low-dimensional fluids. Entropy, 27:455, 2025

    work page 2025

  40. [40]

    J. L. Lebowitz and H. Spohn. Transport properties of the lorentz gas: Fourier’s law. J. Stat. Phys. , 19:633, 1978

    work page 1978

  41. [41]

    Tehver, F

    R. Tehver, F. Toigo, J. Koplik, and J. R. Banavar. Thermal walls in computer simulations. Phys. Rev. E , 57:R17(R), 1998

    work page 1998

  42. [42]

    Heat conduction in two-dimensional momentum-conserving and -nonconserving gases

    Rongxiang Luo. Heat conduction in two-dimensional momentum-conserving and -nonconserving gases. Phys. Rev. E , 102:052104, Nov 2020

    work page 2020

  43. [43]

    Ihle and D

    T. Ihle and D. M. Kroll. Stochastic rotation dynamics: A galilean-invariant mesoscopic model for fluid flow. Phys. Rev. E , 63:020201, Jan 2001

    work page 2001

  44. [44]

    Squires and Stephen R

    Todd M. Squires and Stephen R. Quake. Microfluidics: Fluid physics at the nanoliter scale. Rev. Mod. Phys. , 77:977–1026, Oct 2005

    work page 2005

  45. [45]

    Stefanov, V

    S. Stefanov, V. Roussinov, and C. Cercignani. Rayleigh - benard flow of a rarefied gas and its attractors. i. con- vection regime. Physics of Fluids , 14(7):2255–2269, 07 2002

    work page 2002

  46. [46]

    Stefanov, V

    S. Stefanov, V. Roussinov, and C. Cercignani. Rayleigh - benard flow of a rarefied gas and its attractors. ii. chaotic and periodic convective regimes. Physics of Fluids , 14(7):2270–2288, 07 2002

    work page 2002

  47. [47]

    Rectification of heat current in corbino geometry

    Zhixing Zou, Giulio Casati, Giuliano Benenti, and Jiao Wang. Rectification of heat current in corbino geometry. Phys. Rev. E , 109:L062104, Jun 2024

    work page 2024

  48. [48]

    Thermal Transport in Low Dimen- sions: From Statistical Physics to Nanoscale Heat Trans- fer

    Stefano Lepri (ed.). Thermal Transport in Low Dimen- sions: From Statistical Physics to Nanoscale Heat Trans- fer. Springer Berlin, 2016

    work page 2016

  49. [49]

    GRA VITY-INDUCED MODULATION OF NEGATIVE DIFFERENTIAL THERMAL RESISTANCE IN FLUIDS

    Katja Klinar, Timm Swoboda, Miguel Mu˜ noz Rojo, and Andrej Kitanovski. Fluidic and mechanical ther- mal control devices. Advanced Electronic Materials , 7(3):2000623, 2021. 7 SUPPLEMENTAL MATERIAL FOR “GRA VITY-INDUCED MODULATION OF NEGATIVE DIFFERENTIAL THERMAL RESISTANCE IN FLUIDS” Here we present the detailed derivation of the heat flux and collision f...

    work page 2021

  50. [50]

    1 belongs to the integrable case when the interactions among particle s in the system are not considered

    Derivation of heat flux The fluid model shown in Fig. 1 belongs to the integrable case when the interactions among particle s in the system are not considered. In such a case, the heat flux flowing across th e system from the lower to the upper end can be written as J = N ⏐ ⏐ ⏐ ⏐ ⟨EU→ L⟩ − (⟨EL→ U ⟩ − mgH ) ⟨tU→ L⟩ + ⟨tL→ U ⟩ ⏐ ⏐ ⏐ ⏐ , (A1) where N is the num...

    work page

  51. [51]

    Therefore, if the system consists of N non-interacting particles, the collision frequency of particles collidin g with the upper bath is given by f = N f(1) = N ⟨tL→ U ⟩ + ⟨tU→ L⟩

    Derivation of collision frequency The collision frequency of a particle colliding with the upper bath is defin ed as f (1) = 1 / (⟨tL→ U ⟩ + ⟨tU→ L⟩). Therefore, if the system consists of N non-interacting particles, the collision frequency of particles collidin g with the upper bath is given by f = N f(1) = N ⟨tL→ U ⟩ + ⟨tU→ L⟩ . (A17) Appendix B: NDTR of...

    work page