Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interatomic potentials”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An Ab Initio -Derived Force Field for Amorphous Silica Interfaces for Use in Molecular Dynamics Simulations

Here, we present a classical interatomic force field, silica-DDEC, to describe the interactions of amorphous and crystalline silica surfaces, parametrized using density functional theory-based charges. Charge schemes for silica surfaces were developed using the density-derived electrostatic and chemical (DDEC) method, which reproduces atomic charges of the periodic models as well as the electrostatic potential away from the atom sites. Lennard–Jones parameters were determined by requiring the correct description of (i) the amorphous silica density, coordination defects, and local coordination geometry, relative to experimental measurements, and (ii) water-silica interatomic distances compared with ab initio results. Deprotonated surface silanol sites are also described within the model based on DDEC charges. The result is a general electronic structure-derived model for describing fully flexible amorphous and crystalline silica surfaces and interactions of liquids with silica surfaces of varying structure and protonation state.

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

Phase Boundaries, Isotope Effect and Superconductivity of Lithium Under Hydrostatic Conditions

We present theoretical and experimental studies of superconductivity and low temperature structural phase boundaries in lithium. We mapped the structural phase diagram of 6Li and 7Li under hydrostatic conditions between 5 – 55 GPa and within the temperature range of 15 – 75 K, observing the FCC → hR1 → cI16 phase transitions. 6Li and 7Li show some differences at the structural boundaries, with a potential shift of the phase boundaries of 6Li to lower pressures. Density functional theory calculations and topological analysis of the electron density elucidates the superconducting properties and interatomic interactions within these phases of lithium.

Deemyad, Shanti↗
DOE OSTI2025

Experimental confirmation of first-principles thermal conductivity in Zirconium-doped ThO 2

The degradation of thermal conductivity in advanced nuclear fuels due to the accumulation of fission products and irradiation-induced defects is inevitable, and must be considered as part of safety and efficiency analyses of nuclear reactors. Here, this study examines the thermal conductivity of a zirconium-doped ThO 2 crystal, synthesized via the hydrothermal method using a spatial domain thermoreflectance technique. Zirconium is one of the soluble fission products in oxide fuels that can effectively scatter heat-carrying phonons in the crystalline lattice of fuel. Thus, thermal property measurements of zirconium-doped ThO 2 single crystals provide insights into the effects of substitutional zirconium doping, isolated from extrinsic factors such as grain boundary scattering. The experimental results are compared with first-principles calculations of the lattice thermal conductivity of ThO 2 , employing an iterative solution of the Peierls-Boltzmann transport equation. Additionally, the non-perturbative Green's function methodology is utilized to compute phonon-point defect scattering rates, accounting for local distortions around point defects, including mass difference changes, interatomic force constants, and structural relaxation. The congruence between the predicted results from first-principles calculations and the measured temperature-dependent thermal conductivity validates the computational methodology. Furthermore, the methodologies employed in this study enable systematic investigations of thermal conductivity reduction by fission products, potentially leading to the development of more accurate fuel performance codes.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Effect of Redox-Inactive Metal Ion–Nickel(III) Interactions on the Redox Properties and Proton-Coupled Electron Transfer Reactivity

We report Mononuclear nickel(II) and nickel(III) complexes of a bisamidate-bisalkoxide ligand, (NMe 4 ) 2 [Ni- II (HMPAB)] (1) and (NMe 4 )[Ni III (HMPAB)] (2), respectively, have been synthesized and characterized by various spectroscopic techniques including X-ray crystallography. The reaction of redox-inactive metal ions (M n+ = Ca 2+ , Mg 2+ , Zn 2+ , Y 3+ , and Sc 3+ ) with 2 resulted in 2-M n+ adducts, which was assessed by an array of spectroscopic techniques including X-ray absorption spectroscopy (XAS), electron paramagnetic resonance (EPR), and reactivity studies. The X-ray structure of Ca 2+ coordinated to Ni(III) complexes, 2Ca 2+ T, was determined and exhibited an average Ni-Ca distance of 3.1253 angstrom, close to the metal ions' covalent radius. XAS analysis of 2-Ca 2+ and 2-Y 3+ in solution further revealed an additional coordination to Ca and Y in the 2-M n+ adducts with shortened Ni-M distances of 2.15 and 2.11 Å, respectively, implying direct bonding interactions between Ni and Lewis acids (LAs). Such a short interatomic distance between Ni(III) and M is unprecedented and was not observed before. EPR analysis of 2 and 2-M n+ species, moreover, displayed rhombic signals with g av > 2.12 for all complexes, supporting the +III oxidation state of Ni. The Ni III /Ni II redox potential of 2 and 2-M n+ species was determined, and a plot of E 1/ 2 of 2-M n+ versus pK a of [M(H 2 O) n ] m+ exhibited a linear relationship, implying that the Ni III /Ni II potential of 2 can be tuned with different redox-inactive metal ions. Reactivity studies of 2 and 2-M n+ with different 4-X-2,6-ditert-butylphenol (4-XDTBP) and other phenol derivatives were performed, and based on kinetic studies, we propose the involvement of a proton-coupled electron transfer (PCET) pathway. Analysis of the reaction products after the reaction of 2 with 4-OMe-DTBP showed the formation of a Ni(II) complex (1a) where one of the alkoxide arms of the ligand is protonated. A pK a value of 24.2 was estimated for 1a. The reaction of 2-M n+ species was examined with 4-OMe-DTBP, and it was observed that the k 2 values of 2-M n+ species increase by increasing the Lewis acidity of redox-inactive metal ions. However, the obtained k 2 values for 2-M n+ species are much lower compared to the k 2 value for 2. Such a variation of PCET reactivity between 2 and 2-M n+ species may be attributed to the interactions between Ni(III) and LAs. Our findings show the significance of the secondary coordination sphere effect on the PCET reactivity of Ni(III) complexes and furnish important insights into the reaction mechanism involving high-valent nickel species, which are frequently invoked as key intermediates in Ni-mediated enzymatic reactions, solar-fuel catalysis, and biomimetic/synthetic transformation reactions.

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

Periodic GFN1-xTB Tight Binding: A Generalized Ewald Partitioning Scheme for the Klopman–Ohno Function

A novel formulation is presented for the treatment of electrostatics in the periodic GFN1-xTB tight-binding model. Periodic GFN1-xTB is hindered by the functional form of the second-order electrostatics, which only recovers Coulombic behavior at large interatomic distances and lacks a closed-form solution for its Fourier transform. We address this by introducing a binomial expansion of the Klopman–Ohno function to partition short- and long-range interactions, enabling the use of a generalized Ewald summation for the solution of the electrostatic energy. This approach is general and is applicable to any damped potential of the form |R n + c| –m . Benchmarks on the X23 molecular crystal dataset and a range of prototypical bulk semiconductors demonstrate that this systematic treatment of the electrostatics eliminates unphysical behavior in the equation of state curves. In the bulk systems studied, we observe a mean absolute error in total energy of 35 meV/atom, comparable to the machine-learned universal force field, M3GNet, and sufficiently precise for structure relaxation. These results highlight the promising potential of GFN1-xTB as a universal tight-binding parametrization.

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

Local Atomic Configuration in Pristine and A-Site Doped Silver Niobate Perovskite Antiferroelectrics

Antiferroelectrics have attracted increasing research interests in recent years due to both their great potential in energy storage applications and intriguing structural characteristics. However, the links between the electrical properties and structural characteristics of distorted perovskite antiferroelectrics are yet to be fully deciphered. Here, we adopt local-structure methods to elucidate the nanoscale atomic structure of AgNbO3-based antiferroelectrics and their structural evolution upon La doping. The local structural features including interatomic distance distributions and atomic displacements have been analyzed using neutron small-box pair distribution function (PDF) refinement in conjunction with large-box Reverse Monte Carlo modelling. Our results highlight the correlation of cation displacements in AgNbO 3 and its disruption by the incorporation of La, apparently in corroboration with the observed anomalous dielectric properties. Spatial ordering of cation vacancies is observed in La-doped AgNbO 3 samples, which coordinates with oxygen octahedral tilting to relieve lattice strain. These results provide renewed insights into the atomic structure and antiferroelectric phase instabilities of AgNbO 3 and relevant perovskite materials, further lending versatile opportunities for enhancing their functionalities.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Complexity of many-body interactions in transition metals via machine-learned force fields from the TM23 data set

Abstract This work examines challenges associated with the accuracy of machine-learned force fields (MLFFs) for bulk solid and liquid phases ofd-block elements. In exhaustive detail, we contrast the performance of force, energy, and stress predictions across the transition metals for two leading MLFF models: a kernel-based atomic cluster expansion method implemented using sparse Gaussian processes (FLARE), and an equivariant message-passing neural network (NequIP). Early transition metals present higher relative errors and are more difficult to learn relative to late platinum- and coinage-group elements, and this trend persists across model architectures. Trends in complexity of interatomic interactions for different metals are revealed via comparison of the performance of representations with different many-body order and angular resolution. Using arguments based on perturbation theory on the occupied and unoccupieddstates near the Fermi level, we determine that the large, sharpddensity of states both above and below the Fermi level in early transition metals leads to a more complex, harder-to-learn potential energy surface for these metals. Increasing the fictitious electronic temperature (smearing) modifies the angular sensitivity of forces and makes the early transition metal forces easier to learn. This work illustrates challenges in capturing intricate properties of metallic bonding with current leading MLFFs and provides a reference data set for transition metals, aimed at benchmarking the accuracy and improving the development of emerging machine-learned approximations.

Chemistry↗
DOE OSTI2023

Reconsidering short-range order in complex concentrated alloys

Abstract The seemingly contradictory state of research on short-range order in many-component alloys is addressed through a critical review of the characterization of face-centered-cubic 3 d systems. Despite the paucity of direct observations, the ordering of many widely studied alloys is argued to be primarily interesting for its potential ubiquity. To clarify this situation, future research directions are proposed with reference to historical results, including a review of the fundamental principles of ordering and clustering. Graphical abstract

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ambient and High Pressure CuNiSb 2 : Metal-Ordered and Metal-Disordered NiAs-Type Derivative Pnictides

The mineral Zlatogorite, CuNiSb 2 , was synthesized in the laboratory for the first time by annealing elements at ambient pressure (CuNiSb 2 -AP). Rietveld refinement of synchrotron powder X-ray diffraction data indicates that CuNiSb 2 -AP crystallizes in the NiAs-derived structure ( P 3 m 1, #164) with Cu and Ni ordering. The structure consists of alternate NiSb 6 and CuSb 6 octahedral layers via face-sharing. The formation of such structure instead of metal disordered NiAs-type structure ( P 6 3 / mm c, #194) is validated by the lower energy of the ordered phase by first-principle calculations. Interatomic crystal orbital Hamilton population, electron localization function, and charge density analysis reveal strong Ni-Sb, Cu-Sb, and Cu-Ni bonding and long weak Sb-Sb interactions in CuNiSb 2 -AP. The magnetic measurement indicates that CuNiSb 2 -AP is Pauli paramagnetic. First-principle calculations and experimental electrical resistivity measurements reveal that CuNiSb 2 -AP is a metal. The low Seebeck coefficient and large thermal conductivity suggest that CuNiSb 2 is not a potential thermoelectric material. Single crystals were grown by chemical vapor transport. The high pressure sample (CuNiSb 2 -8 GPa) was prepared by pressing CuNiSb 2 -AP at 700 °C and 8 GPa. However, the structures of single crystal and CuNiSb 2 -8 GPa are best fit with a disordered metal structure in the P 3 m 1 space group, corroborated by transmission electron microscopy.

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

Minimization of Disorder as a Key Design Principle for Natural Sizes of Light Harvesting 2 Complexes

The light harvesting 2 (LH2) complex of purple bacteria has excellent energy conversion efficiency. Clarifying the design principle behind such efficiency at the atomistic level is crucial for understanding its structure–function relationship and can be utilized for the design of artificial light harvesting systems. To this end, we conducted comprehensive computational investigation of the dynamical and statistical nature of electronic excited states of pigment molecules in a natural LH2 complex with 9-fold symmetry and its two non-natural in silico analogues with 6- and 12-fold symmetries. To ensure reliable and efficient all-atomistic molecular dynamics simulations, we combined a well established interpolation approach for the construction of the potential energy surface with a neural network machine learning approach. Outcomes of these calculations clarify that non-natural forms of LH2-type complexes have significantly larger quasistatic disorder than those for the natural one. In addition, non-natural systems have more disruptions of the hydrogen bonding, underscoring its crucial role for reducing the disorder. On the other hand, local environmental dynamics are relatively insensitive to the structural changes although there is moderate enhancement in the anharmonic or interatomic components for the synthetic ones. These findings based on all-atomistic simulations provide direct computational evidence that the structure and sizes of natural LH2 complexes are designed to minimize the energetic disorder. We analyze quantitative implications of these for the energy transferring capability of the LH2 complex.

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

Multitarget Rydberg gates via spatial blockade engineering

Multi-target gates offer the potential to reduce gate depth in syndrome extraction for quantum error correction. Although neutral-atom quantum computers have demonstrated native multi-qubit gates, existing approaches that avoid additional control or multiple atomic species have been limited to single-target gates. We propose single-control-multi-target CZ^n gates on a single-species neutral-atom platform that require no extra control and have gate durations comparable to standard CZ gates. Our approach leverages tailored interatomic distances to create an asymmetric blockade between the control and target atoms. Using a GPU-accelerated pulse synthesis protocol, we design smooth control pulses for CZZ and CZZZ gates, achieving fidelities of up to 99.55% and $99.24\%$, respectively, even in the presence of simulated atom placement errors and Rydberg-state decay. Our approach is most effective for N=2 (CZZ) and N=3 targets (CZZZ); for larger N, increasing spatial crowding of the targets introduces significant challenges for maintaining the required blockade asymmetry. This work presents a practical path to implementing low-overhead multi-target gates in single-species neutral-atom systems, significantly reducing the resource overhead for syndrome extraction. To motivate the impact of these gates, we apply a greedy scheduling algorithm and we demonstrate that our proposed gates can reduce the number of atom reconfiguration costs by up to 50% for color code syndrome extraction of code distances greater than 5.

Stein, Samuel A.↗
DOE OSTI2023

Architector for high-throughput cross-periodic table 3D complex building

Abstract Rare-earth and actinide complexes are critical for a wealth of clean-energy applications. Three-dimensional (3D) structural generation and prediction for these organometallic systems remains a challenge, limiting opportunities for computational chemical discovery. Here, we introduce Architector , a high-throughput in-silico synthesis code for s-, p-, d-, and f-block mononuclear organometallic complexes capable of capturing nearly the full diversity of the known experimental chemical space. Beyond known chemical space, Architector performs in-silico design of new complexes including any chemically accessible metal-ligand combinations. Architector leverages metal-center symmetry, interatomic force fields, and tight binding methods to build many possible 3D conformers from minimal 2D inputs including metal oxidation and spin state. Over a set of more than 6,000 x-ray diffraction (XRD)-determined complexes spanning the periodic table, we demonstrate quantitative agreement between Architector-predicted and experimentally observed structures. Further, we demonstrate out-of-the box conformer generation and energetic rankings of non-minimum energy conformers produced from Architector , which are critical for exploring potential energy surfaces and training force fields. Overall, Architector represents a transformative step towards cross-periodic table computational design of metal complex chemistry.

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

Crystal Prediction and Design of Tunable Light Emission in BTB–Based Metal–Organic Frameworks

Metal-organic frameworks (MOFs) have recently been shown to exhibit unique mechanisms of luminescence based on charge transfer between structural units in the framework. These MOFs have the potential to be structural tuned for targeted emission with little or no metal participation. Here a computationally led, material design and synthesis methodology is presented that elucidates the mechanisms of light emission in interpenetrated structures comprised of metal centers (M = In, Ga, InGa, InEu) and BTB (1,3,5-Tris(4-carboxyphenyl)benzene) linkers, forming unique luminescent M-BTB MOF frameworks. Gas phase and periodic electronic structure calculations indicate that the intensity of the emission and the wavelength are overwhelmingly controlled by a combination of the number of interacting stacked linkers and their interatomic spacings, respectively. In the MOF, the ionic radii of the metal centers primarily control the expansion or shrinkage of the linker stacking distances. Experimentally, multiple M-BTB-based MOFs are synthesized and their photoluminescence was tested. Experiments validated the modeling by confirming that shifts in the crystal structure result in variations in light emission. Through this material design method, the mechanisms of tuning luminescence properties in interpenetrated M-BTB MOFs have been identified and applied to the design of MOFs with specific wavelength emission based on their structure.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Density functional modeling of the binding energies between aluminosilicate oligomers and different metal cations

Interactions between negatively charged aluminosilicate species and positively charged metal cations are critical to many important engineering processes and applications, including sustainable cements and aluminosilicate glasses. In an effort to probe these interactions, here we have calculated the pair-wise interaction energies (i.e., binding energies) between aluminosilicate dimer/trimer and 17 different metal cations M n+ (M n+ = Li + , Na + , K + , Cu + , Cu 2+ , Co 2+ , Zn 2+ , Ni 2+ , Mg 2+ , Ca 2+ , Ti 2+ , Fe 2+ , Fe 3+ , Co 3+ , Cr 3+ , Ti 4+ and Cr 6+ ) using a density functional theory (DFT) approach. Analysis of the DFT-optimized structural representations for the clusters (dimer/trimer + M n+ ) shows that their structural attributes (e.g., interatomic distances) are generally consistent with literature observations on aluminosilicate glasses. The DFT-derived binding energies are seen to vary considerably depending on the type of cations (i.e., charge and ionic radii) and aluminosilicate species (i.e., dimer or trimer). A survey of the literature reveals that the difference in the calculated binding energies between different M n+ can be used to explain many literature observations associated with the impact of metal cations on materials properties (e.g., glass corrosion, mineral dissolution, and ionic transport). Analysis of all the DFT-derived binding energies reveals that the correlation between these energy values and the ionic potential and field strength of the metal cations are well captured by 2nd order polynomial functions ( R 2 values of 0.99–1.00 are achieved for regressions). Given that the ionic potential and field strength of a given metal cation can be readily estimated using well-tabulated ionic radii available in the literature, these simple polynomial functions would enable rapid estimation of the binding energies of a much wider range of cations with the aluminosilicate dimer/trimer, providing guidance on the design and optimization of sustainable cements and aluminosilicate glasses and their associated applications. Finally, the limitations associated with using these simple model systems to model complex interactions are also discussed.

Gong, Kai↗
DOE OSTI2025

Hund’s flat band in a frustrated spinel oxide

Electronic flat bands associated with quenched kinetic energy and heavy electron mass have attracted great interest for promoting strong electronic correlations and emergent phenomena such as high-temperature charge fractionalization and superconductivity. Intense experimental and theoretical research has been devoted to establishing the rich nontrivial metallic and heavy fermion phases intertwined with such localized electronic states. Here, in this study, we investigate the transition metal oxide spinel LiV 2 O 4 , an enigmatic heavy fermion compound lacking localized f orbital states. We use angle-resolved photoemission spectroscopy and dynamical mean-field theory to reveal a kind of correlation-induced flat band with suppressed interatomic electron hopping arising from intra-atomic Hund’s coupling. The appearance of heavy quasiparticles is ascribed to a proximate orbital-selective Mott state characterized by fluctuating local moments as evidenced by complementary magnetotransport measurements. The spectroscopic fingerprints of long-lived quasiparticles and their disappearance with increasing temperature further support the emergence of a high-temperature “bad” metal state observed in transport data. This work resolves a long-standing puzzle on the origin of heavy fermion behavior and unconventional transport in LiV 2 O 4 . Simultaneously, it opens a path to achieving flat bands through electronic interactions in d-orbital systems with geometrical frustration, potentially enabling the realization of exotic phases of matter such as the fractionalized Fermi liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-Principles Calculations of the Structural, Electronic, Optical, and Mechanical Properties of 21 Pyrophosphate Crystals

Pyrophosphate crystals have a wide array of applications in industrial and biomedical fields. However, fundamental understanding of their electronic structure, optical, and mechanical properties is still scattered and incomplete. In the present research, we report a comprehensive theoretical investigation of 21 pyrophosphates A 2 M (H 2 P 2 O 7 ) 2 •2H 2 O with either triclinic or orthorhombic crystal structure. The molecule H 2 P 2 O 7 is the dominant molecular unit, whereas A = (K, Rb, NH 4 , Tl), M = (Zn, Cu, Mg, Ni, Co, Mn), and H 2 O stand for the cation elements, transition metals, and the water molecules, respectively. The electronic structure, interatomic bonding, partial charge distribution, optical properties, and mechanical properties are investigated by first-principles calculations based on density functional theory (DFT). Most of these 21 crystals are theoretically investigated for the first time. The calculated results show a complex interplay between A, M, H 2 P 2 O 7 , and H 2 O, resulting in either metallic, half-metallic, or semi-conducting characteristics. The novel concept of total bond order density (TBOD) is used as a single quantum mechanical metric to characterize the internal cohesion of these crystals to correlate with the calculated properties, especially the mechanical properties. This work provides a large database for pyrophosphate crystals and a road map for potential applications of a wider variety of phosphates.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗
DOE OSTI2025

Influence of isoelectronic mass modulation on phonon decay and transport in layered pnictogen-carbon systems

This study systematically investigates the structural, electronic, and thermal transport properties of layered pnictogen–carbon (Pn 2 C 2 ; Pn = P, As, Sb, Bi) monolayers, revealing a profound influence of isoelectronic mass modulation on phonon dynamics and thermal conductivity. Through first-principles calculations and lattice dynamics analysis, we demonstrate that increasing the atomic mass of pnictogens leads to elongated Pn-C bonds, weakened interatomic interactions, and enhanced phonon-phonon scattering, resulting in a dramatic reduction in lattice thermal conductivity (κ). Specifically, κ decreases exponentially from 65.6 W/m-K in P 2 C 2 to an ultralow 0.37 W/m-K in Bi 2 C 2 , driven by suppressed acoustic-optical phonon gaps Δ$^{Γ}_{A-0}$, increased anharmonicity, and reduced phonon lifetimes and group velocities. The transition from semiconducting to metallic behavior in heavier pnictogen-based systems further enhances phonon-electron scattering, contributing to thermal conductivity suppression. Structural analysis highlights the strengthening of out-of-plane C—C bonds and the contraction of C—Pn—C bond angles, which disrupt in-plane phonon transport. Further, these findings establish a design framework for engineering ultralow-κ materials through chemical substitution and structural tuning, offering significant potential for thermoelectric and thermal management applications. This work provides fundamental insights into phonon decay mechanisms and thermal transport in 2D materials, paving the way for advanced phononic and energy-efficient devices.

2D materials↗