Pub Date : 2025-04-22DOI: 10.1007/s10955-025-03438-5
Elena Rufeil Fiori, Christian Maes, Robbe Vidts
Applying a time-periodic magnetic field to the standard ferromagnetic Curie–Weiss model brings the spin system in a steady out-of-equilibrium condition. We recall how the hysteresis gets influenced by the amplitude and the frequency of that field, and how an amplitude- and frequency-dependent (dynamical) critical temperature can be discerned. The dissipated power measures the area of the hysteresis loop and changes with temperature. The excess heat determines a nonequilibrium specific heat giving the quasistatic thermal response. We compute that specific heat, which appears to diverge at the critical temperature, quite different from the equilibrium case.
{"title":"Specific Heat of the Driven Curie–Weiss Model","authors":"Elena Rufeil Fiori, Christian Maes, Robbe Vidts","doi":"10.1007/s10955-025-03438-5","DOIUrl":"10.1007/s10955-025-03438-5","url":null,"abstract":"<div><p>Applying a time-periodic magnetic field to the standard ferromagnetic Curie–Weiss model brings the spin system in a steady out-of-equilibrium condition. We recall how the hysteresis gets influenced by the amplitude and the frequency of that field, and how an amplitude- and frequency-dependent (dynamical) critical temperature can be discerned. The dissipated power measures the area of the hysteresis loop and changes with temperature. The excess heat determines a nonequilibrium specific heat giving the quasistatic thermal response. We compute that specific heat, which appears to diverge at the critical temperature, quite different from the equilibrium case.</p></div>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":"192 5","pages":""},"PeriodicalIF":1.3,"publicationDate":"2025-04-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143861374","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-04-21DOI: 10.1007/s10955-025-03442-9
Vebjørn H. Bakkestuen, Mihály A. Csirik, Andre Laestadius, Markus Penz
A detailed analysis of density-functional theory for quantum-electrodynamical model systems is provided. In particular, the quantum Rabi model, the Dicke model, and a generalization of the latter to multiple modes are considered. We prove a Hohenberg–Kohn theorem that manifests the magnetization and displacement as internal variables, along with several representability results. The constrained-search functionals for pure states and ensembles are introduced and analyzed. We find the optimizers for the pure-state constrained-search functional to be low-lying eigenstates of the Hamiltonian and, based on the properties of the optimizers, we formulate an adiabatic-connection formula. In the reduced case of the Rabi model we can even show differentiability of the universal density functional, which amounts to unique pure-state v-representability.
{"title":"Density-Functional Theory for the Dicke Hamiltonian","authors":"Vebjørn H. Bakkestuen, Mihály A. Csirik, Andre Laestadius, Markus Penz","doi":"10.1007/s10955-025-03442-9","DOIUrl":"10.1007/s10955-025-03442-9","url":null,"abstract":"<div><p>A detailed analysis of density-functional theory for quantum-electrodynamical model systems is provided. In particular, the quantum Rabi model, the Dicke model, and a generalization of the latter to multiple modes are considered. We prove a Hohenberg–Kohn theorem that manifests the magnetization and displacement as internal variables, along with several representability results. The constrained-search functionals for pure states and ensembles are introduced and analyzed. We find the optimizers for the pure-state constrained-search functional to be low-lying eigenstates of the Hamiltonian and, based on the properties of the optimizers, we formulate an adiabatic-connection formula. In the reduced case of the Rabi model we can even show differentiability of the universal density functional, which amounts to unique pure-state <i>v</i>-representability.</p></div>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":"192 5","pages":""},"PeriodicalIF":1.3,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10955-025-03442-9.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143856621","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-04-21DOI: 10.1007/s10955-025-03443-8
Weifeng Jiang, Zhengxing Lian, Yujun Zhu
In this paper, the rotational entropy (h_r(varphi )) is investigated for a random dynamical system (varphi ) on the torus. The formula of (h_r(varphi )) is obtained for (varphi ) which satisfies certain assumptions, and the lower and upper bounds of (h_r(varphi )) are given for more general (varphi ). Several examples are presented to show that these results may not hold without the assumptions. This work can be seen as a random version of the previous work (Jiang et al. in J Differ Equ 379:862–883, 2024), in which the rotational entropy was introduced and investigated as a homotopy invariant for any torus map.
{"title":"Rotational Entropy for Random Torus Maps","authors":"Weifeng Jiang, Zhengxing Lian, Yujun Zhu","doi":"10.1007/s10955-025-03443-8","DOIUrl":"10.1007/s10955-025-03443-8","url":null,"abstract":"<div><p>In this paper, the rotational entropy <span>(h_r(varphi ))</span> is investigated for a random dynamical system <span>(varphi )</span> on the torus. The formula of <span>(h_r(varphi ))</span> is obtained for <span>(varphi )</span> which satisfies certain assumptions, and the lower and upper bounds of <span>(h_r(varphi ))</span> are given for more general <span>(varphi )</span>. Several examples are presented to show that these results may not hold without the assumptions. This work can be seen as a random version of the previous work (Jiang et al. in J Differ Equ 379:862–883, 2024), in which the rotational entropy was introduced and investigated as a homotopy invariant for any torus map.</p></div>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":"192 5","pages":""},"PeriodicalIF":1.3,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143852594","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-04-21DOI: 10.1007/s10955-025-03445-6
Xiaoxiao Nie, Yu Huang
As an extension of weighted entropy, the weighted topological sequence entropy and the weighted measure-theoretic sequence entropy are defined. A variational principle of relating the two weighted sequence entropies is established. The weighted maximal pattern entropy is also defined. It is shown that for homeomorphism dynamical systems the weighted maximal pattern entropy is equal to the supremum of the weighted sequence entropies over all strictly increasing sequences in integers both in topological and measure-theoretic settings.
{"title":"Weighted Sequence Entropy and Maximal Pattern Entropy","authors":"Xiaoxiao Nie, Yu Huang","doi":"10.1007/s10955-025-03445-6","DOIUrl":"10.1007/s10955-025-03445-6","url":null,"abstract":"<div><p>As an extension of weighted entropy, the weighted topological sequence entropy and the weighted measure-theoretic sequence entropy are defined. A variational principle of relating the two weighted sequence entropies is established. The weighted maximal pattern entropy is also defined. It is shown that for homeomorphism dynamical systems the weighted maximal pattern entropy is equal to the supremum of the weighted sequence entropies over all strictly increasing sequences in integers both in topological and measure-theoretic settings.</p></div>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":"192 5","pages":""},"PeriodicalIF":1.3,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143856671","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-04-17DOI: 10.1007/s10955-025-03436-7
Evan Habbershaw, Cory D. Hauck, Jingwei Hu, Jeffrey R. Haack
A multi-species Fokker–Planck model for simulating particle collisions in a plasma is presented. The model includes various parameters that must be tuned. Under reasonable assumptions on these parameters, the model satisfies appropriate conservation laws, dissipates an entropy, and satisfies an (mathcal {H})-Theorem. In addition, the model parameters provide the additional flexibility that is used to match simultaneously momentum and temperature relaxation formulas derived from the Boltzmann collision operator for a binary mixture with Coulomb potential. A numerical method for solving the resulting space-homogeneous kinetic equation is presented and two examples are provided to demonstrate the relaxation of species bulk velocities and temperatures to their equilibrium values.
{"title":"A Nonlinear, Conservative, Entropic Fokker–Planck Model for Multi-species Collisions","authors":"Evan Habbershaw, Cory D. Hauck, Jingwei Hu, Jeffrey R. Haack","doi":"10.1007/s10955-025-03436-7","DOIUrl":"10.1007/s10955-025-03436-7","url":null,"abstract":"<div><p>A multi-species Fokker–Planck model for simulating particle collisions in a plasma is presented. The model includes various parameters that must be tuned. Under reasonable assumptions on these parameters, the model satisfies appropriate conservation laws, dissipates an entropy, and satisfies an <span>(mathcal {H})</span>-Theorem. In addition, the model parameters provide the additional flexibility that is used to match simultaneously momentum and temperature relaxation formulas derived from the Boltzmann collision operator for a binary mixture with Coulomb potential. A numerical method for solving the resulting space-homogeneous kinetic equation is presented and two examples are provided to demonstrate the relaxation of species bulk velocities and temperatures to their equilibrium values.</p></div>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":"192 5","pages":""},"PeriodicalIF":1.3,"publicationDate":"2025-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143840293","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-04-16DOI: 10.1007/s10955-025-03441-w
Alexander I. Bufetov, Yosuke Kawamoto
We investigate the intertwining of Laguerre processes of parameter (alpha ) in different dimensions. We introduce a Feller kernel that depends on (alpha ) and intertwines the (alpha )-Laguerre process in (N+1) dimensions and that in N dimensions. When (alpha ) is a non-negative integer, the new kernel is interpreted in terms of the conditional distribution of the squared singular values: if the singular values of a unitarily invariant random matrix of order ((N+alpha +1) times (N+1)) are fixed, then the those of its ((N+alpha ) times N ) truncation matrix are given by the new kernel.
{"title":"The Intertwining Property for Laguerre Processes with a Fixed Parameter","authors":"Alexander I. Bufetov, Yosuke Kawamoto","doi":"10.1007/s10955-025-03441-w","DOIUrl":"10.1007/s10955-025-03441-w","url":null,"abstract":"<div><p>We investigate the intertwining of Laguerre processes of parameter <span>(alpha )</span> in different dimensions. We introduce a Feller kernel that depends on <span>(alpha )</span> and intertwines the <span>(alpha )</span>-Laguerre process in <span>(N+1)</span> dimensions and that in <i>N</i> dimensions. When <span>(alpha )</span> is a non-negative integer, the new kernel is interpreted in terms of the conditional distribution of the squared singular values: if the singular values of a unitarily invariant random matrix of order <span>((N+alpha +1) times (N+1))</span> are fixed, then the those of its <span>((N+alpha ) times N )</span> truncation matrix are given by the new kernel.</p></div>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":"192 5","pages":""},"PeriodicalIF":1.3,"publicationDate":"2025-04-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10955-025-03441-w.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143840359","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-04-11DOI: 10.1007/s10955-025-03422-z
Gregory Wimsatt, Alexander B. Boyd, James P. Crutchfield
The Trajectory Class Fluctuation Theorem (TCFT) presents equalities between thermodynamic quantities, such as work costs and free energy changes, and the probabilities of classes of system-state trajectories in equilibrium-steady-state nonequilibrium processes. Conceptually, the TCFT unifies a host of previously-established fluctuation theorems, interpolating from Crooks’ Detailed Fluctuation Theorem (single trajectories) to Jarzynski’s Equality (full trajectory ensembles). Leveraging coarse-grained information about how systems evolve, the TCFT provides a substantial strengthening of the Second Law of Thermodynamics—that, in point of fact, can be a rather weak bound between requisite work and free energy change. It also can be used to improve empirical estimates of free energies, a task known to be statistically challenging, by diverting attention from rare, work-dominant trajectories in convenient but highly nonequilibrium processes. The TCFT also reveals new forms of free energy useful for bounding work costs when computing with systems whose microscopic details are difficult to ascertain—forms that can be solved analytically and practically estimated. For engineered systems more generally, it connects the role of system state trajectories in system functionality to the particular work costs required to evolve those trajectories. Previously, the TCFT was used to connect the microscopic dynamics of experimentally-implemented Josephson-junction information engines with the mesoscopic descriptions of how information was processed. The development here justifies that empirical analysis, explicating its mathematical foundations.
{"title":"Trajectory Class Fluctuation Theorem","authors":"Gregory Wimsatt, Alexander B. Boyd, James P. Crutchfield","doi":"10.1007/s10955-025-03422-z","DOIUrl":"10.1007/s10955-025-03422-z","url":null,"abstract":"<div><p>The Trajectory Class Fluctuation Theorem (TCFT) presents equalities between thermodynamic quantities, such as work costs and free energy changes, and the probabilities of classes of system-state trajectories in equilibrium-steady-state nonequilibrium processes. Conceptually, the TCFT unifies a host of previously-established fluctuation theorems, interpolating from Crooks’ Detailed Fluctuation Theorem (single trajectories) to Jarzynski’s Equality (full trajectory ensembles). Leveraging coarse-grained information about how systems evolve, the TCFT provides a substantial strengthening of the Second Law of Thermodynamics—that, in point of fact, can be a rather weak bound between requisite work and free energy change. It also can be used to improve empirical estimates of free energies, a task known to be statistically challenging, by diverting attention from rare, work-dominant trajectories in convenient but highly nonequilibrium processes. The TCFT also reveals new forms of free energy useful for bounding work costs when computing with systems whose microscopic details are difficult to ascertain—forms that can be solved analytically and practically estimated. For engineered systems more generally, it connects the role of system state trajectories in system functionality to the particular work costs required to evolve those trajectories. Previously, the TCFT was used to connect the microscopic dynamics of experimentally-implemented Josephson-junction information engines with the mesoscopic descriptions of how information was processed. The development here justifies that empirical analysis, explicating its mathematical foundations.</p></div>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":"192 4","pages":""},"PeriodicalIF":1.3,"publicationDate":"2025-04-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10955-025-03422-z.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143822014","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-04-10DOI: 10.1007/s10955-025-03421-0
Jing Ouyang, Changguo Xiao
We are concerned with the boundary layer problem for steady relativistic Boltzmann equation in half-space. By introducing some special test functions and choosing suitable estimate orders, we derive a crucial bound on the macroscopic part (textbf{P}f). Under certain admissible conditions and assuming exponential decay for the source term, we demonstrate the existence, continuity and the exponential spacial decay of a boundary layer solution for both linear and nonlinear steady relativistic Boltzmann equations with hard potential collision kernel. The uniqueness is also established under specific constraint conditions.
{"title":"The Steady Relativistic Boltzmann Equation in Half-Space","authors":"Jing Ouyang, Changguo Xiao","doi":"10.1007/s10955-025-03421-0","DOIUrl":"10.1007/s10955-025-03421-0","url":null,"abstract":"<div><p>We are concerned with the boundary layer problem for steady relativistic Boltzmann equation in half-space. By introducing some special test functions and choosing suitable estimate orders, we derive a crucial bound on the macroscopic part <span>(textbf{P}f)</span>. Under certain admissible conditions and assuming exponential decay for the source term, we demonstrate the existence, continuity and the exponential spacial decay of a boundary layer solution for both linear and nonlinear steady relativistic Boltzmann equations with hard potential collision kernel. The uniqueness is also established under specific constraint conditions.</p></div>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":"192 4","pages":""},"PeriodicalIF":1.3,"publicationDate":"2025-04-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143818311","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-04-05DOI: 10.1007/s10955-025-03424-x
Renato Spacek, Pierre Monmarché, Gabriel Stoltz
In molecular dynamics, transport coefficients measure the sensitivity of the invariant probability measure of the stochastic dynamics at hand with respect to some perturbation. They are typically computed using either the linear response of nonequilibrium dynamics, or the Green–Kubo formula. The estimators for both approaches have large variances, which motivates the study of variance reduction techniques for computing transport coefficients. We present an alternative approach, called the transient subtraction technique (inspired by early work by Ciccotti and Jaccucci in Phys Rev Lett 35(12):789–792, 1975, https://doi.org/10.1103/PhysRevLett.35.789), which amounts to simulating a transient dynamics started off equilibrium and relaxing towards the equilibrium state, from which we subtract a sensibly coupled equilibrium trajectory, resulting in an estimator with smaller variance. We present the mathematical formulation of the transient subtraction technique, give error estimates on the bias and variance of the associated estimator, and demonstrate the relevance of the method through numerical illustrations for various systems.
{"title":"Transient Subtraction: A Control Variate Method for Computing Transport Coefficients","authors":"Renato Spacek, Pierre Monmarché, Gabriel Stoltz","doi":"10.1007/s10955-025-03424-x","DOIUrl":"10.1007/s10955-025-03424-x","url":null,"abstract":"<div><p>In molecular dynamics, transport coefficients measure the sensitivity of the invariant probability measure of the stochastic dynamics at hand with respect to some perturbation. They are typically computed using either the linear response of nonequilibrium dynamics, or the Green–Kubo formula. The estimators for both approaches have large variances, which motivates the study of variance reduction techniques for computing transport coefficients. We present an alternative approach, called the <i>transient subtraction technique</i> (inspired by early work by Ciccotti and Jaccucci in Phys Rev Lett 35(12):789–792, 1975, https://doi.org/10.1103/PhysRevLett.35.789), which amounts to simulating a transient dynamics started off equilibrium and relaxing towards the equilibrium state, from which we subtract a sensibly coupled equilibrium trajectory, resulting in an estimator with smaller variance. We present the mathematical formulation of the transient subtraction technique, give error estimates on the bias and variance of the associated estimator, and demonstrate the relevance of the method through numerical illustrations for various systems.</p></div>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":"192 4","pages":""},"PeriodicalIF":1.3,"publicationDate":"2025-04-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143778035","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-04-05DOI: 10.1007/s10955-025-03417-w
Nitay Ben-Shachar, Joseph T. Johnson, Douglas R. Brumley, Jason Nassios, John E. Sader
Asymptotic analyses of the Boltzmann equation for near-continuum low-Mach-number gas flows predominantly assume diffuse scattering from solid surfaces, i.e., complete surface accommodation, despite gas scattering often deviating from this idealized behavior in practice. While some results for arbitrary surface accommodation exist to second order in small Knudsen number, the full theory to this order is yet to be reported. Here, we present a matched asymptotic expansion of the linearized Boltzmann–BGK equation that generalizes existing theories to Maxwell-type boundary conditions with arbitrary accommodation at solid surfaces. This is performed to second order in small Knudsen number for smooth solid surfaces, and holds for steady and unsteady flow at oscillatory frequencies far smaller than the molecular collision frequency. In contrast to diffuse scattering, we find that the second-order Knudsen layer functions vary as (eta log ^2eta ) for incomplete but arbitrary accommodation at a curved surface, where (eta ) is the dimensionless normal coordinate. A modified refined moment method is developed to numerically handle this spatial dependency. Analytical formulas for all velocity slip and temperature jump coefficients for the Hilbert region are reported that exhibit accuracies greater than 99.9%. This resolves conflicting literature reports on the second-order velocity slip and temperature jump coefficients.
{"title":"Near-Continuum Gas Flows to Second Order in Knudsen Number with Arbitrary Surface Accommodation","authors":"Nitay Ben-Shachar, Joseph T. Johnson, Douglas R. Brumley, Jason Nassios, John E. Sader","doi":"10.1007/s10955-025-03417-w","DOIUrl":"10.1007/s10955-025-03417-w","url":null,"abstract":"<div><p>Asymptotic analyses of the Boltzmann equation for near-continuum low-Mach-number gas flows predominantly assume diffuse scattering from solid surfaces, i.e., complete surface accommodation, despite gas scattering often deviating from this idealized behavior in practice. While some results for arbitrary surface accommodation exist to second order in small Knudsen number, the full theory to this order is yet to be reported. Here, we present a matched asymptotic expansion of the linearized Boltzmann–BGK equation that generalizes existing theories to Maxwell-type boundary conditions with arbitrary accommodation at solid surfaces. This is performed to second order in small Knudsen number for smooth solid surfaces, and holds for steady and unsteady flow at oscillatory frequencies far smaller than the molecular collision frequency. In contrast to diffuse scattering, we find that the second-order Knudsen layer functions vary as <span>(eta log ^2eta )</span> for incomplete but arbitrary accommodation at a curved surface, where <span>(eta )</span> is the dimensionless normal coordinate. A modified refined moment method is developed to numerically handle this spatial dependency. Analytical formulas for all velocity slip and temperature jump coefficients for the Hilbert region are reported that exhibit accuracies greater than 99.9%. This resolves conflicting literature reports on the second-order velocity slip and temperature jump coefficients.</p></div>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":"192 4","pages":""},"PeriodicalIF":1.3,"publicationDate":"2025-04-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143778028","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}