4.7 Article

Phonon Transport Simulator (PhonTS)

Journal

COMPUTER PHYSICS COMMUNICATIONS
Volume 192, Issue -, Pages 196-204

Publisher

ELSEVIER SCIENCE BV
DOI: 10.1016/j.cpc.2015.01.008

Keywords

Thermal conductivity; Boltzmann transport equation

Funding

  1. US Government under DOE [DE-AC07-05ID14517]
  2. US Government under the Energy Frontier Research Center (Office of Science, Office of Basic Energy Science) [FWP 1356]

Ask authors/readers for more resources

Thermal, conductivity prediction remains an important subject in many scientific and engineering areas. Only recently has such prediction become possible on the basis of the first principles calculations, thus ensuring high quality results. Implementation of the methodology, however, is technically challenging and requires a lengthy development process. We thus introduce the Phonon Transport Simulator (PhonTS), a Fortran90, fully parallel code to perform such calculations. PhonTS possesses a large array of options and returns the thermal conductivity tensor together with related quantities, such as spectral thermal conductivity, phonon lifetimes, mean free paths and Gruneisen parameters. First principles calculations are implemented via convenient interfaces to widely-used third-party codes, while many classical potentials are included in PhonTS itself. The code is carefully validated against data published in the literature from various thermal conductivity computational techniques and against experimental data. Program summary Program title: PhonTS Catalogue identifier: AEVO_v1_0 Program-summary-URL: http://cpc.cs.qub.ac.uk/summaries/AEVO_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License, version 3 No. of lines in distributed program, including test data, etc.: 7060788 No. of bytes in distributed program, including test data, etc.: 58119383 Distribution format: tar.gz Programming language: Fortran90. Computer: Linux cluster. Operating system: Linux. Has the code been vectorised or parallelized?: Yes, via the MPI-libraries RAM: 200 MB Classification: 7.8, 7.9. External routines: LAPACK, MPI-libraries Nature of problem: Computes thermal conductivity in crystal solids from the level of the interatomic interactions. Solution method: Iterative or variational solution of the Boltzmann Transport Equation for phonons, by treating anharmonicity as a perturbation to the harmonic problem. Additional comments: Due to the large file size, PhonTS is not delivered directly when download or Email is requested. Instead an html file giving details of how the program can be obtained is sent. Running time: Problem dependent. For solid argon (4 atoms in the conventional cell, k-space mesh of 9 x 9 x 9, classical potential, fully convergent results) simulations take similar to 10(4) seconds on a single processor. (C) 2015 Elsevier B.V. All rights reserved.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

Article Physics, Applied

Thermal conductivity of the n=1-5 and 10 members of the (SrTiO3)nSrO Ruddlesden-Popper superlattices

Natalie M. Dawley, Ella K. Pek, Che-Hui Lee, Eugene J. Ragasa, Xue Xiong, Kiyoung Lee, Simon R. Phillpot, Aleksandr V. Chernatynskiy, David G. Cahill, Darrell G. Schlom

Summary: In this study, the thermal conductivity of (SrTiO3)(n)SrO Ruddlesden-Popper superlattice thin films with different numbers of layers was measured, revealing a decrease in thermal conductivity with increasing interface density. Characterization by X-ray diffraction and scanning transmission electron microscopy confirmed the Ruddlesden-Popper superlattice structure of the samples.

APPLIED PHYSICS LETTERS (2021)

Article Chemistry, Physical

Generation and characterization of an improved carbon fiber model by molecular dynamics

Linyuan Shi, Marina Sessim, Michael R. Tonks, Simon R. Phillpot

Summary: A high-fidelity model of carbon fibers was developed using kinetic Monte Carlo combined with large-scale Molecular Dynamics, with the generated fiber core and thin fiber models showing similar densities and structural characteristics to experimental structures. Introducing artificial defects in the models was shown to heal during structural equilibration, resulting in Young's moduli within the experimental range.

CARBON (2021)

Article Physics, Applied

Lattice stability and point defect energetics of TiSi2 and TiGe2 allotropes from first-principles calculations

David L. Brown, Kevin S. Jones, Simon R. Phillpot

Summary: This study used density functional theory to investigate the phase stabilities and point defect energetics of TiSi2 and TiGe2 allotropes, focusing on the C49 and C54 phases. The results revealed that the ground state structure of TiGe2 is the C54 phase with low sheet resistance, while TiSi2 predominantly forms the less desirable C49 phase with higher resistance. Ge atoms were found to introduce additional covalent bond stability for the C54 phase, impacting its enthalpy of formation.

JOURNAL OF APPLIED PHYSICS (2021)

Article Physics, Applied

Effect of the initial temperature on the shock response of Cu50Zr50 bulk metallic glass by molecular dynamics simulation

Peng Wen, Brian Demaske, Douglas E. Spearot, Simon R. Phillpot, Gang Tao

Summary: The effect of initial temperature on the shock response of Cu50Zr50 bulk metallic glass was investigated through molecular dynamics simulations. It was found that the shock Hugoniot relationship shows temperature dependence, with flow stress decreasing at high temperatures, indicating lower shear resistance. Additionally, a shock-induced melting regime was identified at high temperatures and pressures.

JOURNAL OF APPLIED PHYSICS (2021)

Article Chemistry, Physical

High-temperature oxidation of carbon fiber and char by molecular dynamics simulation

Linyuan Shi, Marina Sessim, Michael R. Tonks, Simon R. Phillpot

Summary: Reactive MD simulations were used to study the initial stage of carbon fiber and amorphous carbon char oxidation with atomic oxygen. Carbon monoxide was found to be the primary product in both systems. Oxygen was adsorbed on the surfaces of both materials, with significantly higher amounts on the surface of amorphous carbon char compared to carbon fiber.

CARBON (2021)

Article Physics, Applied

Molecular dynamics simulation of the shock response of materials: A tutorial

Peng Wen, Gang Tao, Douglas E. Spearot, Simon R. Phillpot

Summary: The aim of this Tutorial is to help new researchers understand how to perform molecular dynamics simulations of the shock response of materials and to highlight the importance of MD in studying shock physics. Two approaches, non-equilibrium and equilibrium MD simulations, for generating shock waves are reviewed. By using analysis techniques, MD simulations can provide a powerful mechanistic understanding of shock that complements shock experiments.

JOURNAL OF APPLIED PHYSICS (2022)

Article Materials Science, Multidisciplinary

Phase-field modeling of carbon fiber oxidation coupled with heat conduction

Marina Sessim, Linyuan Shi, Simon R. Phillpot, Michael R. Tonks

Summary: In this paper, a multiphysics mesoscale model was developed to study the oxidation of carbon fibers in a phenolic impregnated carbon ablator thermal protection system. The presented model used the phase-field method to capture the reduction of carbon fibers due to the oxidation reaction and was fully coupled with heat transport in the system. Sensitivity analysis revealed that the reaction rate parameters have the most impact on the oxidation time.

COMPUTATIONAL MATERIALS SCIENCE (2022)

Article Chemistry, Physical

Ultralow Lattice Thermal Conductivity in the Aikinite Structure Family, CuXPbXBi2-XS3, and Thermoelectric Properties of Cu0.14Pb0.14Bi1.86S3

Srikanth Balijapelly, Ashlee Hauble, Santhoshkumar Sundaramoorthy, Jeremy Lee Watts, Susan M. Kauzlarich, Aleksandr Chernatynskiy, Amitava Choudhury

Summary: In this article, the synthesis, characterization, and ultralow thermal conductivity of three complex quaternary chalcogenide compounds were reported. These compounds possess complex structures with mixed occupied sites, leading to ultralow thermal conductivities. Spectroscopy and charge transport properties reveal narrow band gaps, low electrical conductivities, and high Seebeck coefficients. One of the compounds exhibits a promising figure of merit at high temperatures.

ACS APPLIED ENERGY MATERIALS (2022)

Article Materials Science, Multidisciplinary

A phase-field study of stainless-steel oxidation from high-temperature carbon dioxide

Xueyang Wu, Iman Abdallah, Wen Jiang, Robert S. Ullberg, Simon R. Phillpot, Adrien Couet, John H. Perepezko, Michael R. Tonks

Summary: An electrochemical phase-field model is used to investigate the oxidation mechanisms of the 21-2N valve stainless steel alloy exposed to carbon dioxide at 973 K. The model includes three observed oxide phases: Mn3O4, Cr2O3, and MnCr2O4. The sensitivity of oxidation processes to diffusion mobilities is examined and the oxidation rate is calibrated against experimental data. It is found that both inward oxygen and outward metal diffusion are important for oxidation, and the order of initial oxide layers impacts the diffusion of Mn.

COMPUTATIONAL MATERIALS SCIENCE (2023)

Article Chemistry, Inorganic & Nuclear

NaGaSe2: A Water-Loving Multifunctional Non-van der Waals Layered Selenogallate

Srikanth Balijapelly, Santhoshkumar Sundaramoorthy, Dibya Jyoti Mondal, Sanjit Konar, Nikolay Gerasimchuk, Aleksandr Chernatynskiy, Amitava Choudhury

Summary: A missing member of ternary chalcometallates, NaGaSe2, has been synthesized and its crystal structure has been analyzed. The compound has the ability to absorb water molecules to form hydrated phases, which can undergo reversible structural transformations. The hydrated phase of NaGaSe2 exhibits significantly increased Na ionic conductivity and selective water absorption. The compound also shows the potential for solid-state ion exchange and has a band gap of about 3 eV. The sorption studies confirm its selective water absorption over other solvents.

INORGANIC CHEMISTRY (2023)

Article Materials Science, Multidisciplinary

Dislocation formation in the heteroepitaxial growth of PbSe/PbTe systems

Yang Li, Boyang Gu, Adrian Diaz, Simon R. Phillpot, David L. Mcdowell, Youping Chen

Summary: The paper presents a multiscale study of the kinetic processes of the heteroepitaxial growth of the PbSe/PbTe (111) and PbTe/PbSe(001) systems, using the Concurrent Atomistic-Continuum (CAC) method as the simulation tool. The simulations have reproduced the growth mode and layer morphology observed in experiments, and visualized the formation and evolution of dislocations.

ACTA MATERIALIA (2023)

Article Materials Science, Multidisciplinary

Binding of uranyl cations to a Zr-based metal-organic framework by density functional theory

Yuan Liu, An T. Ta, Shubham Pandey, Kyoung Chul Park, Shenyang Hu, Natalia B. Shustova, Simon R. Phillpot

Summary: The bonding nature of uranyl cations with Zr-MOF was investigated using density functional theory. The results revealed that the binding energy between uranyl cations and Zr-MOF depended on the specific bonding site and the degree of deprotonation in Zr-MOF. These findings are significant for the design of Zr-MOFs with efficient capture of uranyl cations.

COMPUTATIONAL MATERIALS SCIENCE (2023)

Article Materials Science, Multidisciplinary

Migration velocities of intergranular He gas bubbles under thermal gradients in Fe by phase-field modeling

Yixi Shen, Peng Wen, An T. Ta, Simon R. Phillpot, Douglas E. Spearot

Summary: In this study, we use phase-field modeling to calculate the migration velocity of intergranular He bubbles in Fe under thermal gradients. Grain boundaries with energies ranging from 0.88 to 1.6 J/m2 are considered. We address the influences of the phase-field parameterization strategy, thermal conductivity in grain boundaries, and GB diffusion coefficients on the geometry of intergranular He bubbles. Our simulations show good agreement with a theoretical model for the contact angle at the bubble/GB junction. We find that intergranular He bubbles migrate slower than He bubbles in the Fe matrix due to smaller temperature gradients at the bubble/GB junction. We observe a linear relationship between bubble migration velocity and temperature gradient for intergranular He bubbles, consistent with theoretical predictions.

JOURNAL OF NUCLEAR MATERIALS (2023)

Article Chemistry, Applied

Binding of radionuclides and surrogate to 18-crown-6 ether by density functional theory

Yuan Liu, An T. Ta, Kyoung Chul Park, Shenyang Hu, Natalia B. Shustova, Simon R. Phillpot

Summary: We use density functional theory to investigate the interactions of cerium, americium, and curium cations with crown ethers. Our results demonstrate that crown ethers can capture cerium, americium, and curium ions, and modifying the crown ether structure by substituting nitrogen atoms for oxygen atoms significantly increases their binding energies with radionuclides.

MICROPOROUS AND MESOPOROUS MATERIALS (2024)

Article Chemistry, Inorganic & Nuclear

Lithium selenometallates of triel elements, Li5MSe4 (M = Al and Ga), aliovalent doping and their ionic conductivity

Santhoshkumar Sundaramoorthy, Aleksandr V. Chernatynskiy, Nikolay Gerasimchuk, Amitava Choudhury

Summary: Ternary selenometallates, Li5MSe4 (M = Al(I) and Ga(II)), have been synthesized for the first time, and their potential for ionic conductivity and optical properties have been investigated. Aliovalent substitution of Sn has been shown to significantly enhance the ionic conductivity of the materials.

DALTON TRANSACTIONS (2022)

Article Computer Science, Interdisciplinary Applications

Modeling and meshing for tokamak edge plasma simulations

Usman Riaz, E. Seegyoung Seol, Robert Hager, Mark S. Shephard

Summary: The accurate representation and effective discretization of a problem domain into a mesh are crucial for achieving high-quality simulation results and computational efficiency. This work presents recent developments in extending an automated tokamak modeling and meshing infrastructure to better support the near flux field following meshing requirements of the XGC Gyro-kinetic Code.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

Electron-phonon coupling from GW perturbation theory: Practical workflow combining BerkeleyGW, ABINIT, and EPW

Zhenglu Li, Gabriel Antonius, Yang-Hao Chan, Steven G. Louie

Summary: This article presents a workflow for practical calculations of electron-phonon coupling and includes the effect of many-electron correlations using GW perturbation theory. The workflow combines different software packages to enable accurate calculations at the level of quasiparticle band structures.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

MASAP: A package for atomic scattering amplitude in solids

Akihiro Koide, Sara Rabouli, Pierre Le Meur, Sylvain Tricot, Philippe Schieffer, Didier Sebilleau, Calogero R. Natoli

Summary: We present the MsSpec Atomic Scattering Amplitude Package (MASAP), which includes a computation program and a graphical interface for generating atomic scattering amplitude (ASA). The study investigates the applicability of plane wave (PW) and curved spherical wave (SW) scattering in describing electron propagation. The results show that the imaginary part of the optical potential enhances the elastic scattering in the forward direction but causes damping effects in other directions.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

A Bi-directional method for evaluating integrals involving higher transcendental functions. HyperRAF: A Julia package for new hyper-radial functions

A. Bagci, Gustavo A. Aucar

Summary: The electron repulsion integrals over Slater-type orbitals with non-integer principal quantum numbers are investigated in this study. These integrals are important in calculations of many-electron systems. New relationships free from hyper-geometric functions are derived to simplify the calculations. With the use of auxiliary functions and straightforward recurrence relationships, these integrals can be efficiently computed, providing initial conditions for the evaluation of expectation values and potentials.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

First-principles calculations of specular reflection of high-energy electrons during the two-dimensional crystal growth

Andrzej Daniluk

Summary: RHEED_DIFF_2D is an open-source software for qualitative numerical simulations of RHEED oscillation intensity changes with layer deposition, used for interpreting heteroepitaxial structures under different scattering crystal potential models.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

An incremental singular value decomposition approach for large-scale spatially parallel & distributed but temporally serial data - applied to technical flows ☆

Niklas Kuehl, Hendrik Fischer, Michael Hinze, Thomas Rung

Summary: The article presents a strategy and algorithm for simulation-accompanying, incremental Singular Value Decomposition (SVD) for time-evolving, spatially parallel discrete data sets. The proposed method improves computational efficiency by introducing a bunch matrix, resulting in higher accuracy and practical applicability.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

TRAVOLTA: GPU acceleration and algorithmic improvements for constructing quantum optimal control fields in photo-excited systems

Jose M. Rodriguez-Borbon, Xian Wang, Adrian P. Dieguez, Khaled Z. Ibrahim, Bryan M. Wong

Summary: This paper presents an open-source software package called TRAVOLTA for massively parallelized quantum optimal control calculations on GPUs. The TRAVOLTA package is an improvement on the previous NIC-CAGE algorithm and incorporates algorithmic improvements for faster convergence. Three different variants of GPU parallelization are examined to evaluate their performance in constructing optimal control fields in various quantum systems. The benchmarks show that the GPU-enhanced TRAVOLTA code produces the same results as previous CPU-based algorithms but with a speedup of more than ten times. The GPU enhancements and algorithmic improvements allow large quantum optimal control calculations to be efficiently executed on modern multi-core computational hardware.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

MCNOX: A code for computing and interpreting ultrafast nonlinear X-ray spectra of molecules at the multiconfigurational level

Weijie Hua

Summary: This work introduces a program called MCNOX for computing and analyzing ultrafast nonlinear X-ray spectra. It is designed for cutting-edge applications in photochemistry/photophysics enabled by X-ray free-electron lasers and high harmonic generation light sources. The program can calculate steady-state X-ray absorption spectroscopy and three types of ultrafast nonlinear X-ray spectra, and it is capable of identifying major electronic transitions and providing physical and chemical insights from complex signals.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

PLQ-sim: A computational tool for simulating photoluminescence quenching dynamics in organic donor/acceptor blends

Leandro Benatto, Omar Mesquita, Lucimara S. Roman, Rodrigo B. Capaz, Graziani Candiotto, Marlus Koehler

Summary: Photoluminescence Quenching Simulator (PLQ-Sim) is a user-friendly software for studying the dynamics of photoexcited states at the interface between organic semiconductors. It provides important information on organic photovoltaic and photothermal devices and calculates transfer rates and quenching efficiency.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

A method of calculating bandstructure in real-space with application to all-electron and full potential

Dongming Li, James Kestyn, Eric Polizzi

Summary: This study introduces a practical and efficient approach to calculate the all-electron full potential band structure in real space using a finite element basis. Instead of the k-space method, this method solves the Kohn-Sham equation self-consistently within a larger finite system enclosing the unit-cell. Non-self-consistent calculations are then performed in the Brillouin zone to obtain the band structure results, which are found to be in excellent agreement with the pseudopotential k-space method. Furthermore, the study successfully observes the band bending of core electrons.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

EUTERPE: A global gyrokinetic code for stellarator geometry

R. Kleiber, M. Borchardt, R. Hatzky, A. Koenies, H. Leyh, A. Mishchenko, J. Riemann, C. Slaby, J. M. Garcia-Regana, E. Sanchez, M. Cole

Summary: This paper describes the current state of the EUTERPE code, focusing on the implemented models and their numerical implementation. The code is capable of solving the multi-species electromagnetic gyrokinetic equations in a three-dimensional domain. It utilizes noise reduction techniques and grid resolution transformation for efficient computation. Additionally, various hybrid models are implemented for comparison and the study of plasma-particle interactions. The code is parallelized for high scalability on multiple CPUs.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

SMIwiz: An integrated toolbox for multidimensional seismic modelling and imaging

Pengliang Yang

Summary: This paper presents an open source software called SMIwiz, which combines seismic modelling, reverse time migration, and full waveform inversion into a unified computer implementation. SMIwiz supports both 2D and 3D simulations and provides various computational recipes for efficient calculation. Its independent processing and batchwise job scheduling ensure scalability, and its viability is demonstrated through applications on benchmark models.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

Generating and grading 34 optimised norm-conserving Vanderbilt pseudopotentials for actinides and super-heavy elements in the PseudoDojo

Christian Tantardini, Miroslav Ilias, Matteo Giantomassi, Alexander G. Kvashnin, Valeria Pershina, Xavier Gonze

Summary: Material discovery has been an active research field, and this study focuses on developing pseudopotentials for actinides and super-heavy elements. These pseudopotentials are crucial for accurate first-principles calculations and simulations.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

Generalisation of splitting methods based on modified potentials to nonlinear evolution equations of parabolic and Schrödinger type

S. Blanes, F. Casas, C. Gonzalez, M. Thalhammer

Summary: This paper explores the extension of modified potential operator splitting methods to specific classes of nonlinear evolution equations. Numerical experiments confirm the advantages of the proposed fourth-order modified operator splitting method over traditional splitting methods in dealing with Gross-Pitaevskii systems.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

Pole-fitting for complex functions: Enhancing standard techniques by artificial-neural-network classifiers and regressors *

Siegfried Kaidisch, Thomas U. Hilger, Andreas Krassnigg, Wolfgang Lucha

Summary: Motivated by a use case in theoretical hadron physics, this paper revisits an application of a pole-sum fit to dressing functions of a confined quark propagator. Specifically, it investigates approaches to determine the number and positions of singularities closest to the origin for a function known numerically on a specific grid on the positive real axis. Comparing the efficiency of standard techniques to a pure artificial-neural-network approach and a combination of both, it finds that the combined approach is more efficient. This approach can be applied to similar situations where the positions of poles need to be estimated quickly and reliably from real-axis information alone.

COMPUTER PHYSICS COMMUNICATIONS (2024)