Fluctuation-based evaluation of thermal conductivity using deviational phonon Monte Carlo simulations
Abstract
We develop a fluctuation based Monte Carlo method for evaluating thermal conductivity within the phonon Boltzmann transport equation (BTE) framework. Starting from the energy based variance reduced formulation of the BTE, we derive a relation between the spatial spreading of deviational phonon energy and the macroscopic heat conduction equation. The resulting formulation expresses thermal conductivity in terms of the mean square displacement of deviational phonon particles, without imposing thermal reservoirs. Consequently, thermal conductivity along multiple directions can be evaluated from a single simulation. The method is demonstrated for crystalline silicon at room temperature in bulk, square nanowires, and anisotropic polycrystalline nanostructures. The calculated thermal conductivities agree with theoretical predictions, including the size dependence induced by phonon boundary scattering and the mode-resolved contributions. For anisotropic polycrystalline structures, the method simultaneously resolves the directional thermal conductivities and captures their dependence on grain morphology. These results demonstrate the validity of the proposed fluctuation-based approach for evaluating thermal conductivity across different transport regimes and structural configurations.
Related
- Non-destructive evaluation of thermal conductivity of wood and insulation materials using permittivity measurements with ground penetrating radar
- Synergistic Suppression of Ionic Conduction and Lattice Thermal Conductivity Through Trigonally Coordinated [CuI 3 ] 2− Units and Zero‐Dimensionality in Cs 3 Cu 2 I 5
- High-resolution measurements of thermal conductivity matrix and search for thermal Hall effect in La2CuO4
- Predictive Structure to Thermal Conductivity Modeling Framework for BEOL Interconnect Stacks in Advanced Technology Nodes Enabled by Extensive Layer Resolved Thermal Measurements
- Lattice dynamics of disordered phases of silicon and spatially resolved computations of thermal conductivity
BibTeX
@misc{hori2026fluctuationbased,
title = {{Fluctuation-based evaluation of thermal conductivity using deviational phonon Monte Carlo simulations}},
author = {Takuma Hori and Yuzuho Oki},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.14433},
eprint = {2610.14433},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.14433},
}
RIS
TY - UNPB
TI - Fluctuation-based evaluation of thermal conductivity using deviational phonon Monte Carlo simulations
AU - Takuma Hori
AU - Yuzuho Oki
PY - 2026
DA - 2026/10/10
PB - arXiv
DO - 10.48550/arxiv.2610.14433
UR - https://arxiv.org/abs/2610.14433
AB - We develop a fluctuation based Monte Carlo method for evaluating thermal conductivity within the phonon Boltzmann transport equation (BTE) framework. Starting from the energy based variance reduced formulation of the BTE, we derive a relation between the spatial spreading of deviational phonon energy and the macroscopic heat conduction equation. The resulting formulation expresses thermal conductivity in terms of the mean square displacement of deviational phonon particles, without imposing thermal reservoirs. Consequently, thermal conductivity along multiple directions can be evaluated from a single simulation. The method is demonstrated for crystalline silicon at room temperature in bulk, square nanowires, and anisotropic polycrystalline nanostructures. The calculated thermal conductivities agree with theoretical predictions, including the size dependence induced by phonon boundary scattering and the mode-resolved contributions. For anisotropic polycrystalline structures, the method simultaneously resolves the directional thermal conductivities and captures their dependence on grain morphology. These results demonstrate the validity of the proposed fluctuation-based approach for evaluating thermal conductivity across different transport regimes and structural configurations.
ER -