Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Clemson”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Irradiation-induced formation of G-phase precipitates and M 2 X carbides in self-ion irradiated HT-9

Ferritic-martensitic steels with high chromium content are a promising material group for advanced nuclear systems due to their high temperature strength and good irradiation tolerance. HT-9 is an optimized and often-studied alloy in this group, but additional studies are required on its radiation response under extreme conditions to be experienced in various types of nuclear reactors, especially with respect to phase stability under irradiation. Self-ion irradiation of HT-9 by 5 MeV Fe ions was used to simulate neutron-induced behavior reaching peak doses of 100 and 300 dpa at temperatures ranging from 450 to 550 °C. M 23 C 6 carbides that existed prior to irradiation were found to remain stable under all examined irradiation conditions. As irradiation progressed at 450 and 500 °C, however, formation of spherical-like G-phase precipitates and needle-like M 2 X carbides was observed. G-phase precipitates were found to be enriched in Ni, Si, and Mn, and show no interface segregation, whereas needle-like M 2 X carbides were rich in Cr and Mo and clearly displayed interface segregation of Ni and Si. M 2 X carbide formation is believed to be assisted by vacancies, while G-phase precipitation is thought to be assisted by interstitials. Finally, this difference in defect-mediated formation leads to a difference in distribution with depth. M 2 X carbides are distributed over shallower depths than that of G-phase precipitates, consistent with defect imbalance predictions that consider the influence of the injected interstitial effect.

atom probe tomography↗
DOE OSTI2025

X-ray Scintillating Lanthanide-Based Coordination Polymers

Reported are the syntheses and characterization of 25 lanthanide-containing coordination polymers exhibiting four structure types: type I′: [Ln(TFTP) 1.5 (bpy)] (Ln = La 3+ –Nd 3+ ), type I: [Ln(TFTP) 1.5 (bpy)(H 2 O)] (Ln = Sm 3+ –Lu 3+ ), type II′: [Pr(TFTP) 1.5 (phen)], and type II: [Ln(TFTP) 1.5 (phen)(H 2 O)] (Ln = Nd 3+ –Yb 3+ , excluding Pm). These compounds feature lanthanide metal cations bound to tetrafluoroterephthalate linkers (TFTP) and either, 2,2′-bipyridine (bpy) or 1,10-phenanthroline (phen) capping ligands. Initial calculations using crystallographic data revealed large void spaces, ranging from 86 Å 3 in type I to 333 Å 3 in type II. Yet, adsorption isotherm data indicate very little gas uptake of both materials with a BET surface area of 9 m 2 /g for I-Tb (8) and 19 m 2 /g for II-Tb (20). X-ray irradiation of I-Eu (6), II-Eu (18), I-Tb (8), and II-Tb (20) produced visible red and green emission, respectively. Furthermore, radioluminescence studies of both type I and II with Sm 3+ , Eu 3+ , and Tb 3+ metal centers indicate that both structure types featuring Tb 3+ are brighter than those with Eu 3+ or Sm 3+ . Notably, structure type I with Tb 3+ is the brightest scintillator with a relative luminosity of 1.7× that of bismuth germanate (BGO).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Understanding and Tuning the Electronic Properties of Prussian Blue Analogues

Prussian blue analogues (PBAs) have attracted increasing interest owing to their potential applications in various fields such as energy storage and conversion, neuromorphic computing, and magnetic switching. With a general formula of A x M N [M C (CN) 6 ], they feature an open framework that provides abundant channels for diffusion of alkali metal ions A and allows flexible compositional control of transition metal ions M N and M C . The oxidation states of transition metal ions can be tuned by adjusting the amount (x) of alkali ions A. Here, we carried out density functional theory calculations combined with experimental measurements to investigate the effects of transition metal ions, alkali ions, and oxidation states on the electronic properties of PBAs. Our calculations found that the band gaps of PBAs can be tuned from close to 0 eV to more than 4 eV. Experimentally, we introduced the synthesis/characterization of five previously unreported PBAs (M N = Ru, Os; M C = Fe, Ru, and Os) to complete the nine stable M N :M C transition metal combinations in group VIII of the periodic table. The optically measured intervalence charge transfer excitation energies of group VIII PBAs are consistent with calculated band gaps. They demonstrate wide band gap tunability by adjusting transition metals and oxidation states, enabling semiconductor-to-metal transitions for memristor applications and enhancing electronic conductivity for battery applications. In conclusion, this work provides a computational/experimental database of electronic properties versus structural compositions for PBAs.

Electrical conductivity↗
DOE OSTI2025

Multiple Sources of Riparian Wetland Suspended Solids during Episodic Rain Events: Influence on Uranium Transport

Suspended solids can be the primary vector for transporting contaminants in streams. The objective of this study was to determine whether changes in the properties of suspended solids during rain events impacted contaminant transport. Stream water was collected during five episodic events downstream from a U-contaminated wetland located in South Carolina, USA. The suspended particles were initially composed of Fe-flocs (particles formed in situ prior to the rain event) that had significantly greater Fe, Mn, organic-C, and U content than particles collected later during a sampling rain event. XANES and EXAFS revealed that U in the Fe-flocs was U(VI) and that it was not incorporated in a mineral structure but existed as inner- or outer-sphere adsorbed uranyl species associated with organic matter and Fe-oxides. The uranyl had an extraordinarily high affinity for the suspended solids, with solid to liquid U ratios of >72,000 (μg/kg)/(μg/L). After the initial flush of Fe-flocs, a greater fraction of the suspended solids had lower organic-C, Fe, Mn, and amorphous phases and were composed of more quartz, kaolinite, and gibbsite, resulting in lower U concentrations than those in the solids collected earlier in the rain event. This study highlights the importance of understanding suspended solids as transport vectors and their potential dynamic nature during rain events.

computer simulations↗
DOE OSTI2026

An Algorithm for Atom-Centered Lossy Compression of the Atomic Orbital Basis in Density Functional Theory Calculations

Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-ζ (QZ) size are required to adequately converge Kohn–Sham density functional theory (DFT) calculations toward the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. Here, in this work, we develop and test a threshold-based natural atomic orbital (NAO) scheme in which ϵ-NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. This enables compression of the AO basis that is optimal for a given threshold, 10 –ϵ , by discarding NAOs with occupation numbers below that threshold. Extensive pilot test calculations using the Hartree–Fock functional and taking the converged density matrix as input suggest that a threshold of 10 –5 can yield a compression factor (ratio of AO to compressed ϵ-NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold 10 –7 , while the compression factor is typically between 2 and 2.5.

basis sets↗
DOE OSTI2025

Crystal Nucleation Kinetics and Mechanism: Influence of Interaction Potential

Modulating liquid-to-solid transitions and the resulting crystalline structure for tailored properties is much desired. Colloidal systems are exemplary to this end, but the fundamental knowledge gaps in relating the influence of intermolecular interactions to crystallization behavior continue to hinder progress. In this study, we address this knowledge gap by studying nucleation and growth in systems with modified Lennard-Jones potential. Specifically, we study the commonly used 12-6 potential and a softer 7-6 potential. The thermodynamic state point for the study is chosen such that both systems are investigated at the same level of supercooling and pressure. Under these conditions, we find that the nucleation rate for both systems is comparable. Interestingly, the nucleation pathways and resulting crystal structures are different. In the 12-6 system, nucleation and growth occur predominantly through the FCC structure. Softening the potential alters the critical nucleus composition and introduces two distinct nucleation pathways. One pathway predominantly leads to the nucleus with a body-centered cubic (BCC) structure, while the other favors the face-centered cubic (FCC) arrangement. Our study illustrates that polymorph selection can be achieved through modifications to intermolecular interactions without impacting nucleation kinetics. The results have significant implications in designing approaches for polymorph selection and modulating self-assembly mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Revolutionizing Methane Transformation with the Dual Production of Aromatics and Electricity in a Protonic Ceramic Electrocatalytic Membrane Reactor

Reducing the energy and carbon intensity of the conventional chemical processing industry can be achieved by electrochemically transforming natural gases into higher-value chemicals with higher efficiency and near-zero emissions. In this work, the direct conversion of methane to aromatics and electricity has been achieved in a protonic ceramic electrocatalytic membrane reactor through the integration of a proton-conducting membrane assembly and a trimetallic Pt–Cu/Mo/ZSM-5 catalyst for the nonoxidative methane dehydro-aromatization reaction. In this integrated system, a remarkable 15.6% single-pass methane conversion with an 11.4% benzene yield has been demonstrated, while a peak power density of 276 mW cm –2 is obtained at 700 °C. The enhanced 15.7% increase in conversion and 16.0% improvement in the yield are observed when compared with the thermochemical process, which is attributed to the shift of reaction equilibrium by the removal of hydrogen through the protonic membrane. Concurrently, the faster H2 removal at a higher electrical current gave rise to a higher methane conversion and benzene yield. Furthermore, the catalyst can be efficiently regenerated by eliminating carbon deposition. A stable cell potential is maintained for 45 h under a constant current load of 0.13 A cm –2 . Lastly, the dual production of aromatics and electricity in the electrocatalytic membrane reactor has been demonstrated to be an attractive approach for decarbonizing chemical processing.

aromatic compounds↗
DOE OSTI2026

Nanohybrid of Cu2O-Ti3C2Tx as a Silver-Free MXene Sorbent for Iodine Gas Capture from Nuclear Waste

The capture of volatile radioiodine from nuclear fuel reprocessing off-gas streams remains a critical challenge due to the high volatility, long half-life of 129I, and biological uptake of iodide from the environment. Although silver-based sorbents provide strong iodine chemisorption, their high cost and regulatory classification as mixed radioactive-hazardous waste motivate the development of alternative materials. Here, we report a silver-free Cu2O-Ti3C2Tx MXene hybrid for iodine gas capture at 150 °C. Structural and compositional analyses confirm the formation of Cu2O nanoparticles on Ti3C2Tx nanosheets and their subsequent conversion to thermodynamically stable CuI upon static iodine gas exposure, achieving an iodine mass loading of up to 1115 mg/g. These results demonstrate the potential of Cu2O-Ti3C2Tx MXene as a copper-based alternative to silver sorbents for elevated-temperature iodine gas capture.

iodine gas capture↗
DOE OSTI2025

Investigation of Americium-Containing Phosphates, Silicates, Borates, Molybdates, and Fluorides Synthesized via High-Temperature Flux Crystal Growth

The crystal chemistry of americium-containing extended structures was investigated, and several classes of americium-containing solid-state oxide materials were obtained in single-crystal form via high-temperature flux crystal growth. This enabled the structural characterization of rare examples of ternary, quaternary, and penternary americium-containing silicates K 3 Am (Si 2 O 7 ) and Cs 6 Am 2 Si 21 O 48 , phosphates Na 3 Am (PO 4 ) 2 and K 3 Am (PO 4 ) 2 , borates Ba 3 Am 2 (BO 3 ) 4 and AmBO 3 , borate halides Ca 5 Am(BO 3 ) 4 Cl, molybdates Li 0.5 Am 0.5 MoO 4 , and fluorides CsAm 2 F 7 . Using these crystallographic data, the ionic radii of Am 3+ with coordination numbers of six (0.975 Å), seven (1.052 Å), and nine (1.162 Å) were established. A maximum entropy method (MEM) analysis was performed on the single-crystal X-ray diffraction data that were collected for K 3 Nd(PO 4 ) 2 /K 3 Am(PO 4 ) 2 , K 3 NdSi 2 O 7 /K 3 AmSi 2 O 7 , and NdBO 3 /AmBO 3 , to qualitatively compare the ionicities of the Nd–O and Am–O bonds. In conclusion, Raman spectroscopy data were collected on single crystals of K 3 Am(PO 4 ) 2 and compared to the calculated Raman spectrum of K 3 Am(PO 4 ) 2 obtained from DFT calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Strong electron–phonon coupling in magic-angle twisted bilayer graphene

The unusual properties of superconductivity in magic-angle twisted bilayer graphene (MATBG) have sparked considerable research interest. However, despite the dedication of intensive experimental efforts and the proposal of several possible pairing mechanisms, the origin of its superconductivity remains elusive. Here, by utilizing angle-resolved photoemission spectroscopy with micrometre spatial resolution, we reveal flat-band replicas in superconducting MATBG, where MATBG is unaligned with its hexagonal boron nitride substrate. These replicas show uniform energy spacing, approximately 150 ± 15 meV apart, indicative of strong electron–boson coupling. Strikingly, these replicas are absent in non-superconducting twisted bilayer graphene (TBG) systems, either when MATBG is aligned to hexagonal boron nitride or when TBG deviates from the magic angle. Calculations suggest that the formation of these flat-band replicas in superconducting MATBG are attributed to the strong coupling between flat-band electrons and an optical phonon mode at the graphene K point, facilitated by intervalley scattering. These findings, although they do not necessarily put electron–phonon coupling as the main driving force for the superconductivity in MATBG, unravel the electronic structure inherent in superconducting MATBG, thereby providing crucial information for understanding the unusual electronic landscape from which its superconductivity is derived.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

International database of reference gamma spectra for nuclear safeguards applications

Nuclear safeguards missions use gamma spectroscopy as a non-destructive measurement technique for examining nuclear materials. Despite advances in the development of detection equipment as well as software codes, one of the concerns is the lack of well-documented spectra needed to test and validate isotopic analysis codes for their applicability. To address this need, this work introduces IDB, an international database of the reference gamma spectra of uranium, plutonium and mixed oxide nuclear material samples. IDB provides access to well-characterized sets of gamma spectra described by rich metadata, including information on the sample, measurement configuration and detector specifications. These spectra are accessible in different formats, also compatible with analysis code standards, thus promoting their sustainability and maintenance.

None, Dipti [International Atomic Energy Agency (I↗
DOE OSTI2026

Nanoscale control over single vortex motion in an unconventional superconductor

Precise control of superconducting vortices is crucial for studying vortex dynamics and vortex braiding. We propose a new method to pull vortex lines in the layered superconductor FeSe using a scanning tunneling microscope (STM) tip. Weak contact of the STM tip with the FeSe surface locally reduces the superconducting gap, thereby creating a tunable vortex pinning potential on the nanometer scale. This enables controlled vortex line deformation even in dense vortex lattices. Analytical modeling reveals that the deformation strength scales logarithmically with conductance and depends on tip geometry. Our findings point to local strain-induced gap suppression as the mechanism for STM-mediated vortex manipulation, providing a fundamental insight relevant to vortex behavior in the context of quantum information studies.

Song, Sang Yong [Oak Ridge National Laboratory (OR↗
DOE OSTI2024

The putative center in NGC 1052

Context. Many active galaxies harbor powerful relativistic jets, however, the detailed mechanisms of their formation and acceleration remain poorly understood.Aims. To investigate the area of jet acceleration and collimation with the highest available angular resolution, we study the innermost region of the bipolar jet in the nearby low-ionization nuclear emission-line region (LINER) galaxy NGC 1052.Methods. We combined observations of NGC 1052 taken with VLBA, GMVA, and EHT over one week in the spring of 2017. Our study is focused on the size and continuum spectrum of the innermost region containing the central engine and the footpoints of both jets. We employed a synchrotron-self absorption model to fit the continuum radio spectrum and we combined the size measurements from close to the central engine out to ∼1 pc to study the jet collimation.Results. For the first time, NGC 1052 was detected with the EHT, providing a size of the central region in-between both jet bases of 43 μas perpendicular to the jet axes, corresponding to just around 250 RS (Schwarzschild radii). This size estimate supports previous studies of the jets expansion profile which suggest two breaks of the profile at around 3 × 103 RS and 1 × 104 RS distances to the core. Furthermore, we estimated the magnetic field to be 1.25 Gauss at a distance of 22 μas from the central engine by fitting a synchrotron-self absorption spectrum to the innermost emission feature, which shows a spectral turn-over at ∼130 GHz. Assuming a purely poloidal magnetic field, this implies an upper limit on the magnetic field strength at the event horizon of 2.6 × 104 Gauss, which is consistent with previous measurements.Conclusions. The complex, low-brightness, double-sided jet structure in NGC 1052 makes it a challenge to detect the source at millimeter (mm) wavelengths. However, our first EHT observations have demonstrated that detection is possible up to at least 230 GHz. This study offers a glimpse through the dense surrounding torus and into the innermost central region, where the jets are formed. This has enabled us to finally resolve this region and provide improved constraints on its expansion and magnetic field strength.Key words: methods: observational / techniques: high angular resolution / techniques: interferometric / galaxies: active / galaxies: jets / galaxies: Seyfert⋆⋆ NASA Hubble Fellowship Program, Einstein Fellow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The persistent shadow of the supermassive black hole of M87. II. Model comparisons and theoretical interpretations

The Event Horizon Telescope (EHT) observation of M87∗ in 2018 has revealed a ring with a diameter that is consistent with the 2017 observation. The brightest part of the ring is shifted to the southwest from the southeast. In this paper, we provide theoretical interpretations for the multi-epoch EHT observations for M87∗ by comparing a new general relativistic magnetohydrodynamics model image library with the EHT observations for M87∗ in both 2017 and 2018. The model images include aligned and tilted accretion with parameterized thermal and nonthermal synchrotron emission properties. The 2018 observation again shows that the spin vector of the M87∗ supermassive black hole is pointed away from Earth. A shift of the brightest part of the ring during the multi-epoch observations can naturally be explained by the turbulent nature of black hole accretion, which is supported by the fact that the more turbulent retrograde models can explain the multi-epoch observations better than the prograde models. The EHT data are inconsistent with the tilted models in our model image library. Assuming that the black hole spin axis and its large-scale jet direction are roughly aligned, we expect the brightest part of the ring to be most commonly observed 90 deg clockwise from the forward jet. This prediction can be statistically tested through future observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A multifrequency study of sub-parsec jets with the Event Horizon Telescope

Context. The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed, including the main science targets, Sgr A* and M 87, and various calibrators. Sixteen sources were AGN.Aims. We investigated the morphology of the sixteen AGN in the EHT 2017 data set, focusing on the properties of the VLBI cores: size, flux density, and brightness temperature. We studied their dependence on the observing frequency in order to compare it with the Blandford-Königl (BK) jet model. In particular, we aimed to study the signatures of jet acceleration and magnetic energy conversion.Methods. We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components (1749+096, 1055+018, BL Lac, J0132–1654, J0006–0623, CTA 102, and 3C 454.3) and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2–86 GHz range. Combining these data into a multifrequency EHT+ data set, we studied the dependences of the VLBI core component flux density, size, and brightness temperature on the frequency measured in the AGN host frame (and hence on the distance from the central black hole), characterizing them with power law fits. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole.Results. Our observations spanning event horizon to parsec scales indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Only some of the discrepancies may be alleviated by tweaking the model parameters or the jet collimation profile. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations. For our sample, we estimate a general radial dependence of the Doppler factor δ ∝ r≤0.5. This interpretation is consistent with a magnetically accelerated sub-parsec jet. We also estimate a steep decrease of the magnetic field strength with radius B ∝ r−3, hinting at jet acceleration or efficient magnetic energy dissipation.Key words: techniques: interferometric / galaxies: active / galaxies: jets / galaxies: nuclei / quasars: general / quasars: supermassive black holes⋆⋆ Deceased.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

The stellar origins of 96 Zr excesses in presolar graphites from the Murchison meteorite

Context. Zirconium-96 is a stable isotope that can be synthesized under different neutron-rich nucleosynthetic conditions. Astrophysical models predict its production to occur in various stellar environments: from low-to-intermediate-mass asymptotic giant branch (AGB) stars to massive stars and core-collapse supernovae. Aims. Detections of 96Zr excesses, in combination with other isotopic measurements from presolar grains can provide unique constraints on its stellar origin. Presolar grains are microscopic particles found in primitive Solar System materials, which formed in stellar winds and supernova ejecta. The isotopic composition of each grain can provide us a snapshot of the nucleosynthetic processes that took place during the parent star’s lifetime. Methods. In this study, we measured the stable isotopes of C, N, O, Mo, Zr, and Ru in high-density presolar graphite grains from the Murchison meteorite and found four grains that contain positive isotopic anomalies in 96 Zr carried by their internal subgrains. We analyzed multi-element isotopic datasets from each grain to explore the source of the observed 96 Zr excesses. Results. Comparisons with stellar models indicate that two grains likely condensed in an intermediate-mass AGB star with initial metallicity of Z ≤ 0.014. Their 96 Zr/ 94 Zr ratios also match those predicted for born-again AGB stars undergoing a very late thermal pulse and rapidly accreting white dwarfs. After comparing the relative populations of the aforementioned dust-producing stars, we propose rapidly accreting white dwarfs as a new, and more likely, stellar source for one of the presolar grains. The remaining two grains could have originated in the supernova ejecta of massive stars, due to correlated excesses in the p-nuclides, 92, 94 Mo. Thus, grains with 96 Zr anomalies can have a variety of stellar origins, in agreement with theoretical studies. Conclusions. Our study highlights the importance of multi-element analysis in constraining the types of stars where presolar grains have condensed. These data will help improve our understanding of various nucleosynthesis processes in different stellar phases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Advancing microelectronics through nanoscale science: A perspective on needs and opportunities from the nanoscale science research centers

Microelectronics are the cornerstone of the modern world, enhancing our daily lives by providing services such as communications and datacenters. These resources are accessible thanks to the continual pursuit of a deeper understanding of the chemical and physical phenomena underlying the materials synthesis approaches and fabrication processes used to create microelectronic components and subsequently the components' responses to electrical, optical, and other stimuli that are utilized within microelectronic systems. Today, further development of microelectronics requires multidisciplinary expertise across scientific disciplines and fields of study—synthesis, materials characterization, nanoscale fabrication, and performance characterization—with focus placed on comprehending the nanoscale forms and features of microelectronic components. The Nanoscale Science Research Centers (NSRCs) are Department of Energy, Office of Science user facilities that support the international scientific community in advancing nanoscale science and technology. As a key component of the U.S. Government's National Nanotechnology Initiative, the NSRCs enable transformative discoveries by providing world-class facilities, expertise, and collaborative opportunities. Here, in this perspective, we showcase a non-exhaustive cross-section of the capabilities housed at and developed by the NSRCs and their user communities to address fundamental synthesis, metrology, fabrication, and performance considerations toward advancing the development of new microelectronics. Finally, we provide a timely outlook on the next major areas of necessary development in nanoscale sciences to continue the innovation of microelectronics into the next generation.

2D materials↗
DOE OSTI2025

Phonon Olympics: Phonon property and lattice thermal conductivity benchmarking from open-source packages

Three widely used open-source packages for determining phonon properties and lattice thermal conductivities (ALAMODE, phono3py, and ShengBTE) are benchmarked by teams of expert users and the package developers. The phonons for Ge, RbBr, monolayer MoSe 2 , and AlN are modeled at zero temperature, and they scatter through three-phonon and phonon-isotope processes, with thermal conductivities obtained from the linearized Peierls–Boltzmann transport equation with input from density functional theory calculations. Over a wide range of temperatures, the thermal conductivities calculated by the teams fall within at most ±15% of their mean values for each of the four materials. The phonon frequencies, obtained from the harmonic force constants, do not show large differences between the calculations, indicating that the modal heat capacities and group velocities are not responsible for the thermal conductivity variations. It is the lifetimes associated with three-phonon scattering, obtained from the cubic force constants, that drive the variations. The many decisions required to calculate the cubic force constants (e.g., supercell size, atomic displacement, neighbor cutoff, and application of symmetries) make identification of the precise origin of the thermal conductivity variations challenging. The calculated thermal conductivities do not generally show agreement with experimental measurements, which is attributed to the limitations of the density functional theory calculations. Guidance for the development of best practices is provided, which will help to standardize protocols needed for building thermal conductivity databases. The results provide a baseline for future benchmarking of other packages and more advanced calculations.

McGaughey, Alan J. H. [Carnegie Mellon Univ., Pitt↗