Official Resources
- Homepage: https://phonopy.github.io/phono3py/
- Documentation: https://phonopy.github.io/phono3py/
- Source Repository: https://github.com/phonopy/phono3py
- License: BSD 3-Clause License
Overview
phono3py is a code for computing lattice thermal conductivity and related properties from first principles using three-phonon interactions. It extends Phonopy to include anharmonic effects via third-order force constants, enabling calculations of phonon lifetimes, thermal conductivity, and other anharmonic properties essential for thermoelectric materials and thermal management applications.
Scientific domain: Anharmonic lattice dynamics, thermal conductivity, phonon-phonon interactions
Target user community: Researchers studying thermal transport, thermoelectric materials, and anharmonic phonon properties
Theoretical Methods
- Third-order force constants (3rd IFCs)
- Phonon-phonon interaction strength
- Phonon Boltzmann transport equation (BTE)
- Relaxation time approximation (RTA)
- Iterative solution of BTE
- Phonon linewidth and lifetime
- Lattice thermal conductivity tensor
- Mode-dependent thermal conductivity
- Cumulative thermal conductivity
- Group velocities from phonon dispersion
Capabilities (CRITICAL)
- Lattice thermal conductivity calculation
- Phonon lifetimes and linewidths
- Phonon-phonon scattering rates
- Temperature-dependent thermal conductivity
- Directional thermal conductivity (tensor)
- Mode-resolved contributions
- Cumulative thermal conductivity
- Spectral thermal conductivity
- Phonon-phonon interaction strengths
- Phonon linewidths and lifetimes
- Lattice thermal conductivity tensor (RTA and full BTE solution)
- Cumulative thermal conductivity analysis
- Mode-resolved thermal conductivity contributions
- Thermal conductivity vs temperature
- Isotope scattering effects
- Boundary scattering models
- Group velocity and mean free path
- Phase space volume for three-phonon processes
- Mode Gruneisen parameters from 3rd order IFCs
- Interface to multiple DFT codes (same as phonopy)
- HDF5 output for large datasets
- Parallelization support
Sources: Official phono3py documentation, GitHub, cited in 7/7 source lists
Key Strengths
- Production quality: Standard tool for thermal conductivity
- phonopy integration: Seamless workflow with phonopy
- Comprehensive: Full BTE solution with mode-resolved analysis
- Well-documented: Extensive documentation and examples
- Multi-code support: Interfaces to 12+ DFT codes
Inputs & Outputs
-
Input formats:
- POSCAR (VASP structure format)
- FORCES_FC3 (forces on displaced atoms for 3rd order)
- FORCES_FC2 (forces for 2nd order force constants)
- fc3.hdf5 (precomputed third-order force constants)
- fc2.hdf5 (precomputed second-order force constants)
- DFT code outputs via interfaces
-
Output data types:
- kappa-*.hdf5 (thermal conductivity data)
- fc3.hdf5 (third-order force constants)
- gamma-*.hdf5 (phonon linewidths)
- Thermal conductivity vs temperature
- Mode-resolved thermal conductivity
- Cumulative thermal conductivity plots
Interfaces & Ecosystem
Workflow and Usage
Typical phono3py Workflow:
# 1. Generate displaced supercells
phono3py -d --dim="2 2 2" -c POSCAR
# 2. Run DFT on all displaced structures
# (many calculations required)
# 3. Create FORCES_FC3
phono3py --cf3 disp-{00001..NNNNN}/vasprun.xml
# 4. Calculate thermal conductivity
phono3py --mesh="11 11 11" --fc3 --fc2 --br
Advanced Features
- Iterative BTE: Full solution beyond RTA
- Spectral analysis: Frequency-resolved thermal conductivity
- Cumulative functions: Mean free path accumulation
- Isotope scattering: Natural isotope disorder effects
- Boundary scattering: Size-dependent thermal conductivity
- Python API: Programmatic access for automation
Performance Characteristics
- Computational cost: Very high for 3rd order IFCs
- Scalability: Handles standard systems; large systems challenging
- Parallelization: OpenMP support
- Typical runtime: Days to weeks for complete workflow
Computational Cost
- 3rd order force constants: Extremely expensive (N³ scaling)
- BTE solution: Hours to days
- Dense q-grids increase cost significantly
- Supercell size critical for accuracy
Limitations & Known Constraints
- Computational cost: Extremely expensive; requires forces for numerous displaced supercells (scales as N³)
- Supercell size: Large supercells needed for convergence; 2×2×2 or 3×3×3 minimum for many systems
- Cutoff distance: Third-order cutoff must be carefully converged; long-range interactions may be important
- Three-phonon processes only: Neglects four-phonon and higher-order processes (important at high T)
- Harmonic phonon requirement: Requires stable harmonic phonons; unstable modes cause failures
- Classical treatment: Uses classical Bose-Einstein statistics; quantum corrections not included
- Isotope scattering: Simplified model; detailed isotope configurations not considered
- Boundary scattering: Phenomenological models; not fully first-principles
- Memory: HDF5 files can become very large for fine q-point meshes
- Convergence: Requires extensive convergence testing (supercell size, cutoffs, q-mesh)
Verification & Sources
Primary sources:
- Official documentation: https://phonopy.github.io/phono3py/
- GitHub repository: https://github.com/phonopy/phono3py
- A. Togo et al., Phys. Rev. B 91, 094306 (2015) - phono3py methodology
- L. Chaput, Phys. Rev. Lett. 110, 265506 (2013) - Direct BTE solution method
Secondary sources:
- phono3py examples and tutorials
- Comparison studies with ShengBTE and ALAMODE
- High-throughput thermal conductivity databases using phono3py
- Confirmed in 7/7 source lists (claude, g, gr, k, m, q, z)
Confidence: CONFIRMED - Appears in all 7 independent source lists
Verification status: ✅ VERIFIED
- Official homepage: ACCESSIBLE
- Documentation: COMPREHENSIVE and ACCESSIBLE
- Source code: OPEN (GitHub)
- Community support: Active (GitHub issues, shared with phonopy)
- Academic citations: >400 (Google Scholar)
- DFT interfaces: Verified for 10+ codes
Best Practices
- Converge supercell size carefully
- Systematic q-point grid convergence
- Test RTA vs iterative BTE
- Validate against experimental data
- Appropriate cutoff distances for 3rd order IFCs
- Use HDF5 for large datasets
Community and Support
- Open-source (BSD 3-Clause)
- Active GitHub development
- Comprehensive documentation
- Large user community (shared with phonopy)
- Mailing list and GitHub issues
- Workshop materials available
Development
- Atsushi Togo (lead developer)
- Active development
- Regular updates
- Standard tool in thermal transport community
Research Impact
phono3py is the standard tool for first-principles lattice thermal conductivity calculations, widely cited in thermoelectric and thermal transport literature with >400 citations.