Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first-principles model”

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 361 records · Page 20

Non-Electricity Based Renewable Fuels: Theory and Computation for Solar Thermochemical Hydrogen

Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox active defects start to interact with each other, thereby impeding the formation of additional defects. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.

first-principles calculations↗

Be₂C as a neutron moderator for molten salt reactors: A first-principles study of structural, electronic, and defect properties

Beryllium carbide (Be 2 C), valued for its high neutron moderation efficiency and low absorption cross section, is a promising high-temperature neutron moderator for molten salt reactors. Its practical adoption, however, demands significant technological maturation, requiring comprehensive theoretical and experimental studies of its response to different conditions, including high temperature and irradiation. Here, we report initial results on the fundamental properties and radiation-induced defects of Be 2 C, focusing on antisites, vacancies, interstitial atoms, and Frenkel pairs in the Be and C sublattices. Using density functional theory (DFT) and ab initio molecular dynamics (AIMD), we calculate the defects formation and binding energies, evaluating their dependence on the supercell size, charge states, and chemical environment. In general, carbon defects exhibit higher formation energies, greater sensitivity to cell size, and stronger impacts on the density of states compared to beryllium defects, with charged state the effects being more pronounced. Static DFT reveals multiple metastable interstitial configurations, while AIMD identifies ground states as C-C <100> dumbbells and octahedral Be interstitials. In conclusion, the diversity of metastable configurations and defect states complicates the diffusion mechanisms, requiring further molecular dynamics analysis to elucidate the mechanisms and rates of radiation-induced atomic transport, as well as the structural stability of Be2C.

Ab initio modeling↗

Momentum-Resolved Direct Observation of Chiral Phonons in Elemental Tellurium

Chiral phonons carry finite phonon angular momentum, yet their momentum-resolved behavior in acoustic branches remains largely unexplored. In this work, we report the direct detection of chiral acoustic phonons in elemental Te. Through a combination of high-energy-resolution inelastic x-ray scattering experiments and atomistic modeling, we confirm that longitudinal acoustic modes develop pronounced circular atomic motion by mixing with transverse modes. Our results establish a direct link between unexpected scattering intensity in nominally forbidden geometries and circular phonons, offering a framework for characterizing momentum-resolved phonon chirality and highlighting the emergence of intrinsic phonon angular momentum in acoustic branches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effect of XC functionals and dispersion corrections on the DFT‐computed structural and vibrational properties of SrCl 2 –NaCl and ZrF 4 –LiF

Density functional theory (DFT) calculations were performed to examine the impact of exchange–correlation (XC) functionals and van der Waals corrections (specifically the D3 method) on the structural and vibrational properties of the SrCl 2 –NaCl and ZrF 4 –LiF salt systems. Multiple XC functionals, including the local density approximation (LDA), the generalized gradient approximation using the Perdew–Burke–Ernzerhof (PBE) model, and its modified form suitable for solids (PBEsol), the dispersion-corrected PBE-D3 and PBEsol-D3, were considered. Of these functionals, LDA was found to exhibit the highest degree of error, while PBEsol and PBE-D3 displayed the least error. Underestimated lattice parameters compared with experimental values were observed to result in higher force constants, leading to an overprediction of vibrational frequencies. Conversely, an overestimation of lattice parameters was associated with lower vibrational frequencies. The methodology presented in this study yielded results that are in good agreement with experiment, irrespective of the method (finite differences vs. density functional perturbation theory) employed for calculating infrared and Raman spectra. It was further demonstrated that for alkali halides with weak Raman scattering, utilizing a supercell constructed from primitive cells better predicts Raman features than does the use of conventional cells.

Raman↗

Origin of light-induced metastability in ZrTe 5

Here we study the nonequilibrium electronic structure of a model Dirac semimetal ZrTe 5 by using time-and-angle resolved photoemission spectroscopy and density functional theory–based electron and phonon calculations. By measuring the electronic dispersion near the Γ point at time delays up to 10 picoseconds, we discovered that the band spectral weight does not recover during the measured temporal window, revealing the existence of a light-induced metastable state in the electronic structure of this material. Our calculations find that the photoexcited A 1⁢g phonon mode leads to a band renormalization that both supports our experimental observations at the zone center and predicts changes to the band structure outside of our experimental window, ultimately showing the evolution from a direct to an indirect gap semimetal; such band renormalization dramatically reduces the electron-hole recombination rate giving rise to the metastability in this system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

First-Principles Simulation of Beam-Induced Processes Underlying Atomic Manipulation in Electron Microscopes

The development of experimental methods and apparatuses capable of promoting atomically precise material manipulations holds great promise for realizing the ultimate limit of feature miniaturization in materials and devices. The ability to modify materials atom by atom is anticipated to usher in new technologies in areas as diverse as separation science, medicine, and quantum information science. Historically, scanning probe-based techniques have been the most prominent approaches in this space. However, these methods are best suited for the manipulation of surface-exposed regions of materials, as the strong perturbations required for bond scission are delivered most effectively to atoms in the near-proximity to the scanning probe. In contrast, convergent electron beams with energies tuned slightly below the threshold for inducing irreversible knock-on damage have recently been employed (within scanning transmission electron microscopy) to promote atomic-scale bond rearrangements in various beam-stable solids. Currently, however, the efficiency and selectivity of beam-induced atomic manipulation processes with focused electron beams are such that long irradiation times are required to induce a desired atomic rearrangement. With a better understanding of the underlying physics dictating the outcome of a given irradiation event, methods can be devised to improve the efficiency of these techniques so that their promise can be fully realized through widespread adoption.To this end, this Account details our recent efforts to develop and apply tractable first-principles simulation approaches for studying the response of materials to electric beam-like external electric potentials applied in real space. We briefly review the concepts and capabilities in the area of atomically precise materials manipulation and review the early demonstrations of accomplishments in this area, focusing on studies using scanned convergent electron beam probes in particular. We expound upon the depth of the challenge and identify critical shortcomings of theoretical methods that have previously been employed in the simulation of beam-induced processes. We then describe the computational methods that we have generalized from the concepts and tools most commonly applied to the study of molecular photochemistry and how our adaptations of these methods can be employed to capture the relevant dynamical phenomena for beam-induced processes ranging from the initial electron scattering to the ensuing multistate reactions. Here, we contextualize these methods within the current state of the art in this area, which has historically focused primarily on the simulation of inelastic image formation in the electron microscope for the purpose of interpreting the results of quantitative electron microscopy experiments. We demonstrate that the spatial distribution of state-specific excitation rates due to the presence of an external (probe) electric charge is inhomogeneous, such that irradiation at particular locations in materials can favor specific electronic transitions (and disallow others). In addition to the potential for excited-state reaction pathways to be accessed through the initial inelastic scattering of the tightly focused electron beam from the targeted atoms, we also identify favorable conditions for the electronically nonadiabatic evolution of the highly vibrationally excited system to open complex multistate reaction pathways. Implications of the early results for understanding the mechanisms and potential routes to improved efficiency and selectivity in beam-induced reactions are discussed. We conclude with a summary of the current state of theory and modeling capabilities in this area and provide our perspective on future directions for theoretical and experimental developments that we view as crucial to advancing the use of convergent electron beams in mode-specific, atomically precise platforms for direct-write materials modifications.

36 MATERIALS SCIENCE↗

Effective many-body interactions in reduced-dimensionality spaces through neural network models

Accurately describing properties of challenging problems in physical sciences often requires complex mathematical models that are unmanageable to tackle head on. Therefore, developing reduced-dimensionality representations that encapsulate complex correlation effects in many-body systems is crucial to advance the understanding of these complicated problems. However, a numerical evaluation of these predictive models can still be associated with a significant computational overhead. To address this challenge, in this paper we discuss a combined framework that integrates recent advances in the development of active-space representations of coupled cluster (CC) downfolded Hamiltonians with neural network approaches. The primary objective of this effort is to train neural networks to eliminate the computationally expensive steps required for evaluating hundreds or thousands of Hugenholtz diagrams, which correspond to multidimensional tensor contractions necessary for evaluating a many-body form of downfolded effective Hamiltonians. Using small molecular systems (the H 2 O and HF molecules) as examples, we demonstrate that training neural networks employing effective Hamiltonians for a few nuclear geometries of molecules can accurately interpolate or extrapolate their forms to other geometrical configurations characterized by different intensities of correlation effects. We also discuss differences between effective interactions that define CC downfolded Hamiltonians with those of bare Hamiltonians defined by Coulomb interactions in the active spaces. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗

Multimodal Approach Reveals the Symmetry-Breaking Pathway to the Broken Helix in EuIn 2 ⁢As 2

Understanding and manipulating emergent phases, which are themes at the forefront of quantum-materials research, rely on identifying their underlying symmetries. This general principle has been particularly prominent in materials with coupled electronic and magnetic degrees of freedom, in which magnetic order influences the electronic band structure and can lead to exotic topological effects. However, identifying symmetry of a magnetically ordered phase can pose a challenge, particularly in the presence of small domains. Here we introduce a multimodal approach for determining magnetic structures, which combines symmetry-sensitive optical probes, scattering, and group-theoretical analysis. We apply it to EuIn 2 ⁢As 2 , a material that has received attention as a candidate axion insulator. While first-principles calculations predict this state on the assumption of a simple collinear antiferromagnetic structure, subsequent neutron-scattering measurements reveal a much more intricate magnetic ground state characterized by two coexisting magnetic wave vectors reached by successive thermal phase transitions. The proposed high- and low-temperature phases are a spin helix and a state with interpenetrating helical and Néel antiferromagnetic order termed a “broken helix,” respectively. Employing a multimodal approach, we identify the magnetic structure associated with these two phases of EuIn 2 ⁢As 2 . We find that the higher-temperature phase is characterized by a variation of the magnetic moment amplitude from layer to layer, with the moment vanishing entirely in every third Eu layer. The lower-temperature structure is similar to the broken helix, with one important difference: Because of local strain, the relative orientation of the magnetic structure and the lattice is not fixed. Consequently, the symmetry required to protect the axion phase is not generically protected in EuIn 2 ⁢As 2 , but we show that it can be restored if the magnetic structure is tuned with uniaxial strain. Finally, we present a spin Hamiltonian that identifies the spin interactions that account for the complex magnetic order in EuIn 2 ⁢As 2 . Our work highlights the importance of a multimodal approach in determining the symmetry of complex order parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning

Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Panoramic View of MXenes via an Atomic Coordination‐Based Design Strategy

Two‐dimensional (2D) transition metal carbides and nitrides, known as MXenes, possess unique physical and chemical properties, enabling diverse applications in fields ranging from energy storage to communication, catalysis, sensing, healthcare, and beyond. Despite extensive research and notable advancements, a fundamental understanding of MXenes’ phase diversity and its connection to their hierarchical precursors, including the intermediate MAX phases and the ancestral bulk phases, remains limited. Here, in this study, it is hypothesized that the atomic coordination environments adopted by transition metal and nonmetallic atoms in their three‐dimensional (3D) bulk precursors may persist in 2D MXenes to govern their phase diversity. Using high‐throughput modeling based on first‐principles density functional theory, a wide range of MXene phases is unveiled and comprehensively evaluate their relative stabilities across a large chemical space. The key to the approach lies in considering various atomic coordination environments drawn from four types of ancestral bulk phases. Through this comprehensive structural library of MXenes, general guiding principles are uncovered, such as a close alignment between the phase stability of MXenes and that of their 3D precursors. These findings introduce a new design strategy in which the atomic coordination environments in bulk phases can serve as reliable predictors for accessing the diverse structural landscape of MXenes.

MXenes↗

First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport

The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

First-principles study of the structures and redox mechanisms of Ni-rich lithium nickel manganese cobalt oxides

To reduce the cobalt (Co) content in lithium-ion batteries, Ni-rich (high-Ni) lithium nickel manganese cobalt oxides (NMC) are pursued as one of the next-generation cathode materials. However, there is still debate on the crystal and electronic structures of the baseline, LiNiO 2 . Density Functional Theory (DFT) calculations were performed to provide a theoretical understanding of Ni-rich NMC. First, it was found that the commonly used $R\bar{3}m$ structure for LiNiO 2 is metallic, contrary to the experimentally reported mix-conducting behavior. Among the four different space groups, $R\bar{3}m$, C2/m, P2 1 /c, and P2/c, P2/c with charge disproportionation of Ni 2+ and Ni 4+ is the most energetically stable and semiconducting structure of LiNiO 2 . Therefore, the atomic structures of representative Ni-rich NMC were built by partially replacing Ni with Co or Mn in the P2/c LiNiO 2 to form Li x Ni y Mn z Co 1-y-z O 2 . In the fully lithiated (x=1.0) high Ni content NMC (y>0.5), the oxidation state of all Mn ions becomes 4+, while Co ions still maintain 3+, and part of the Ni ions become 3+ to compensate for the charge. Upon delithiation, the local environment shows more variation of the charge states on the transition metal (TM) ions. The average oxidation on each TM follows a sequence of losing electrons that starts from Ni 2+ to Ni 3+ , then oxidizing Ni 3+ and Co 3+ , while Mn 4+ remains electrochemically inactive till x=0. Finally, a general relationship for the oxidation state change in each TM as a function of x is derived and shows agreement with both modeling and experimental data.

25 ENERGY STORAGE↗

Python Library for Monte Carlo Simulations with Ab Initio and Machine-Learned Interatomic Potentials

There is a growing need in the simulation community for software that provides a transparent, reproducible, usable, and extensible (TRUE) Monte Carlo (MC) simulation framework employing energies from ab initio methods and machine-learning interatomic potentials (MLIPs). We introduce a Python library (ASE-MC) that adds Monte Carlo functionality to the Atomic Simulation Environment (ASE) package. Now, we can combine the powerful tools used to build systems and perform ab initio and MLIP in ASE with MC simulation algorithms to sample the configurational space with a concise Python script. After presenting the design philosophy, we demonstrate the flexibility of our approach using selected examples. These example simulations include liquid water described with a message-passing MLIP in the canonical and isothermal–isobaric ensembles, sampling the characteristic dihedral angle of biphenyl and comparing an MLIP to first-principles calculations, and a grand canonical Monte Carlo simulation of ammonia adsorption on Pt(111). These examples showcase the main features of the software, which include flexibility in the choice of ab initio or MLIP engine, ab initio or MLIP grand canonical MC with cavity bias insertions and deletions, the ability to add custom MC moves to the move set, and how users can condense complex MC workflows into a single Python script. Finally, this library serves as a framework for reproducible Monte Carlo simulations, facilitating easy reproduction of the work and application to new systems.

97 MATHEMATICS AND COMPUTING↗

Designing multicomponent hydrides with potential high T c superconductivity

While hydrogen-rich materials have been demonstrated to exhibit high T c superconductivity at high pressures, there is an ongoing search for ternary, quaternary, and more chemically complex hydrides that achieve such high critical temperatures at much lower pressures. First-principles searches are impeded by the computational complexity of solving the Eliashberg equations for large, complex crystal structures. Here, we adopt a simplified approach using electronic indicators previously established to be correlated with superconductivity in hydrides. This is used to study complex hydride structures, which are predicted to exhibit promisingly high critical temperatures for superconductivity. In particular, we propose three classes of hydrides inspired by the Fm$\overline{3}$ m RH 3 structures that exhibit strong hydrogen network connectivity, as defined through the electron localization function. The first class [RH 11 X 3 Y] is based on a Pm$\overline{3}$m structure showing moderately high T c , where the T c estimate from electronic properties is compared with direct Eliashberg calculations and found to be surprisingly accurate. The second class of structures [(RH 11 ) 2 X 6 YZ] improves on this with promisingly high density of states with dominant hydrogen character at the Fermi energy, typically enhancing T c . The third class [(R 1 H 11 )(R 2 H 11 )X 6 YZ] improves the strong hydrogen network connectivity by introducing anisotropy in the hydrogen network through a specific doping pattern. These design principles and associated model structures provide flexibility to optimize both T c and the structural stability of complex hydrides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Simulations of edge and SOL turbulence in diverted negative and positive triangularity plasmas

Optimizing the performance of magnetic confinement fusion devices is critical to achieving an attractive fusion reactor design. Negative triangularity (NT) scenarios have been shown to achieve excellent levels of energy confinement, while avoiding edge localized modes. Modeling turbulent transport in the edge and SOL is key in understanding the impact of NT on turbulence and extrapolating the results to future devices and regimes. Previous gyrokinetic turbulence studies have reported beneficial effects of NT across a broad range of parameters. However, most simulations have focused on the inner plasma region, neglecting the impact of NT on the outermost edge. In this work, we investigate the effect of NT in edge and scrape-off layer simulations, including the magnetic X-point and separatrix. For the first time, we employ a multi-fidelity approach, combining global, non-linear gyrokinetic simulations with drift-reduced fluid simulations, to gain a deeper understanding of the underlying physics at play. First-principles simulations using the GENE-X code demonstrate that in comparable NT and PT geometries, similar profiles are achieved, while the turbulent heat flux is reduced by more than 50% in NT. Comparisons with results from the drift-reduced fluid turbulence code GRILLIX suggest that the turbulence is driven by trapped electron modes. The parallel heat flux width on the divertor targets is reduced in NT, primarily due to a lower spreading factor S.

GENE-X↗

The search for high-entropy fuel-cell catalysts using disorder descriptors

The transition to a hydrogen economy depends on efficient, affordable catalysts for fuel cells. Platinum—the industry standard for fuel-cell electrodes—is costly and scarce, highlighting the need for practical alternatives. High-entropy alloys offer vast compositional diversity and tunable properties that can mitigate these issues, yet their chemical complexity and configurational disorder have hindered rational discovery. Here, we introduce a data-driven framework that couples machine learning with first-principles disorder descriptors—including the entropy forming ability, disordered enthalpy-entropy descriptor, and electronic-structure similarity metrics to platinum—to predict alloy synthesizability and catalytic performance. These descriptors are applied for the first time in the context of fuel-cell catalyst discovery. The workflow rapidly screens more than 20 000 compositions and identifies several platinum-free candidates that are economically viable, readily scalable, and exhibit promising predicted activity. These results demonstrate that disorder descriptors are reliably predicted by machine learning models and can be effectively integrated into materials-discovery pipelines, accelerating innovation across complex compositional spaces.

fuel-cell catalysts↗

Group Conductivity and Nonadiabatic Born Effective Charges of Disordered Metals, Warm Dense Matter, and Hot Dense Plasma

The average ionization state is a critical parameter in plasma models for charged particle transport, equations of state, and optical response. The dynamical or nonadiabatic Born effective charge (NBEC), calculated via first principles time-dependent density-functional theory, provides exact ionic partitioning of bulk electron response for both metallic and insulating materials. The NBEC can be transformed into a “group conductivity,” i.e., the electron conductivity ascribed to a subset of ions. We show that for disordered metallic systems, such as warm dense matter (WDM) and hot dense plasma, the static limit of the NBEC is different from the average ionization states, but that the ionization state can be extracted from the group conductivity even in mixed systems. Here we demonstrate this approach using a set of archetypical examples, including cold and warm aluminium, low- and high-density WDM carbon, and a WDM carbon-beryllium-hydrogen mixture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Synthesis and characterization of the novel breathing pyrochlore compound Ba 3 ⁢Tm 2⁢ Zn 5⁢ O 11

In this study, a novel material from the rare-earth based breathing pyrochlore family, Ba 3 ⁢Tm 2 ⁢Zn 5 ⁢O 11 , was successfully synthesized. Powder x-ray diffraction and high-resolution powder neutron diffraction confirmed phase purity and the $F\bar4$3⁢𝑚 breathing pyrochlore crystal structure, while thermogravimetric analysis revealed incongruent melting behavior compared to its counterpart, Ba 3 ⁢Yb 2⁢ Zn 5⁢ O 11 . High-quality single crystals of Ba 3 ⁢Tm 2 ⁢Zn 5 ⁢O 11 were grown using the traveling solvent floating zone technique and assessed using Laue x-ray diffraction and single crystal x-ray diffraction. Further, thermodynamic characterization indicated paramagnetic behavior down to 0.05 K, and inelastic neutron scattering measurements identified distinct dispersionless crystal electric field energy bands, with the fitted crystal electric field model predicting a single-ion singlet ground state and an energy gap of ∼9 meV separating it from the first excited (singlet) state. Additional low-energy excitation studies on single crystals revealed dispersionless bands at 0.8 and 1 meV. Computed phonon dispersions from first-principles calculations ruled out phonons as the origin of these modes, further illustrating the puzzling and unique properties of Ba 3 ⁢Tm⁢ 2 Zn 5 ⁢O 11 .

36 MATERIALS SCIENCE↗