Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Material Point Method”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Leveraging Electrons for Electrochemical CO 2 Capture Using a Hemi‐Labile Iron Complex

Climate change, driven by anthropogenic carbon emissions, demands urgent action to prevent a 2050 tipping point. With CO 2 levels at 427 ppm (50% above pre-industrial levels), deploying energy-efficient carbon capture technologies is crucial. Electrochemical carbon capture processes that have been touted to have the potential to meet these needs rely on the applied cell voltage, and electron utilization (CO 2 molecules separated per electron), which has generally been asserted to have a theoretical limit of one. Here, we introduce an electron-leveraging strategy to enhance electron utilization beyond this limit to 1.43 by employing Fe-EDDHA, a redox-active coordination complex having a ligand with multiple hemi-labile coordination sites. The reversibility and robustness of the system were enabled by the efficient prevention of CO 2 reduction upon the introduction of nicotinamide as a guardian of the iron(2+) center. The proof-of-concept cyclic system exhibits a minimum operational energy of 22.6 kJ e mol −1 and an average of 63.7 kJ e mol −1 over 29 cycles, using a simulated flue gas (15% CO 2 ). Our electron-leveraging strategy holds promise for advancing energy-efficient electrochemical carbon capture technologies, and offers an alternative to prevalent redox potential shifting methods proposed to mitigate undesired electron transfer reactions in redox-active materials across diverse operational conditions.

carbon storage↗
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fundamental Interactions of Bimetallic Cu x Pd y ( x + y = 4) Clusters Supported on the α-WC(0001) Surface and Their Performance for CO 2 Adsorption and Dissociation

The tungsten carbide α-WC(0001) surface, an active system for the activation of H 2 and important hydrogenation processes involving unsaturated hydrocarbons, can serve as a support of bimetallic clusters to produce materials with unique catalytic properties, opening routes for a wide range of technical applications. In particular, Cu x Pd y clusters are of particular interest because they combine metals with different properties. A stochastic method was applied to obtain the geometry of Cu x Pd y (x + y = 4) bare clusters, evaluating thousands of possibilities to obtain stable structures, yielding one isomer for Cu 4 , Cu 2 Pd 2 , Cu 1 Pd 3 , and Pd 4 and two isomers for Cu 3 Pd 1 . These clusters were supported on C and W terminations of the tungsten carbide (0001) surface, exploring all of the binding possibilities. The adsorption energies on the C and W terminations are in the ranges from −2.51 to −3.02 eV and from −2.26 to −3.30 eV, respectively. The strongest and weakest binding was seen for monometallic Cu 4 and Pd 4 clusters on both C and W terminations, while the Cu-Pd bimetallics have intermediate adsorption energies but lack a clear trend in terms of composition. The location of Cu x Pd y clusters over the (0001) surface induces a decrease in the work function relative to the pristine surface, while the cluster-surface Bader charge transfer and variations in the partial density of states point to changes in the electronic structure of the carbide atoms upon binding of the metallic clusters. The d-band center of the Cu x Pd y deposited on WC(0001) indicates an intermediate reactivity among Cu(111) and Pd(111) surfaces, modulating the reactivity with small numbers of Cu and Pd atoms, i.e., atom economy in catalyst design. The likelihood of existence of the most stable Cu x Pd y (x + y = 4) clusters in the temperature range of 298-400 K is 100%. The composite Cu x Pd y /α-WC(0001) (x + y = 4), is a nontrivial system since 22 isomers are needed to completely describe its structural properties. Among the isomers, seven structures are necessary to represent Cu 3 Pd 1 /α-WC(0001), five for Pd 4 /α-WC(0001), two for Cu 4 /α-WC(0001), and four for Cu 2 Pd 2 /α-WC(0001) and Cu 1 Pd 3 /α-WC(0001). The large number of cluster isomers supported on the tungsten carbide surface opens the door for several applications in the heterogeneous catalysis of the Cu x Pd y /α-WC(0001) composite, with the possibility of modulating the geometric, electronic, and chemical properties according to a desired application. Test studies for the adsorption of CO 2 indicate that the Cu x Pd y /α-WC(0001) composites are highly active for the adsorption and decomposition of the molecule, with bimetallic and admetal-carbide interactions playing a key role in the binding performance. In conclusion, this high activity indicates that these systems should be useful as catalysts for the conversion of CO 2 to oxygenates or light alkanes.

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

Heat Transfer Experiments of Ribbed, Serpentine Cooling Passages with Supercritical CO2

The results of heat transfer experiments with supercritical CO2 (sCO2) to capture the performance of ribbed, serpentine cooling passages for Reynolds numbers in the range of 100,000 - 400,000 are presented. The oxy-combustion turbine operating in an Allam-Fetvedt cycle has a turbine inlet temperature exceeding 1100°C, necessitating internally cooled turbine blades due to material creep strength limits. Reynolds numbers for sCO2, with its combination of high fluid density and low viscosity, can be up to an order of magnitude higher than comparable internal cooling conditions for an air-breathing engine where the majority of experimental data exists. To improve the prediction of cooling performance through turbulated passages, a test rig was designed and constructed with interchangeable insert capability for operation in sCO2 with heat transfer from a heater outlet flow stream to a lower temperature recuperator outlet flow stream. The inserts included serpentine passage geometry with five passes, one of the inserts including chevron ribs and the other a plain wall. The multi-pass serpentine geometry allowed for the evaluation of the effects 180 deg. tip turns have on the return passage flow characteristics, as they can depart significantly from fully-developed flow characteristics. The test procedure included the use of a modified Wilson plot method, with hot flow stream conditions kept effectively constant so that monitored changes in the overall thermal resistance could be attributed to manipulations in cooling flow stream conditions. Conditions during testing consisted of test section pressures of approximately 180 bar, cooling flow stream temperature of 175°C and hot flow stream temperature of 410°C. Test points at discrete cool flow Reynolds numbers up to 400,000 demonstrated an expected trend of decreased fluid thermal resistance with increasing Reynolds number, while at comparable Reynolds number there was an exhibited lower thermal resistance and higher passage pressure drop for the ribbed insert compared to the plain-walled insert. Processing of test data included establishment of Nusselt number enhancement ratios, so that for different internal cooling passage dimensions in the actual turbine blade at comparable Reynolds number the ratio could be applied to predict cooling performance.

Marshall, Michael↗
DOE OSTI2022

Hot-Casting-Assisted Liquid Additive Engineering for Efficient and Stable Perovskite Solar Cells

High-performance inorganic–organic lead halide perovskite solar cells (PSCs) are often fabricated with a liquid additive such as dimethyl sulfoxide (DMSO), which retards crystallization and reduces roughness and pinholes in the perovskite layers. However, DMSO can be trapped during perovskite film formation and induce voids and undesired reaction byproducts upon later processing steps. Here, it is shown that the amount of residual DMSO can be reduced in as-spin-coated films significantly through use of preheated substrates, or a so-called hot-casting method. Hot casting increases the perovskite film thickness given the same concentration of solutions, which allows for reducing the perovskite solution concentration. By reducing the amount of DMSO in proportion to the concentration of perovskite precursors and using hot casting, it is possible to fabricate perovskite layers with improved perovskite–substrate interfaces by suppressing the formation of byproducts, which increase trap density and accelerate degradation of the perovskite layers. The best-performing PSCs exhibit a power conversion efficiency (PCE) of 23.4% (23.0% stabilized efficiency) under simulated solar illumination. Furthermore, encapsulated devices show considerably reduced post-burn-in decay, retaining 75% and 90% of their initial and post-burn-in efficiencies after 3000 h of operation with maximum power point tracking (MPPT) under high power of ultraviolet (UV)-containing continuous light exposure.

perovskite-substrate interface↗
DOE OSTI2025

Strong Kitaev Interaction in BaCo 2⁢ (AsO 4 ) 2

The inelastic neutron scattering results and their analysis unequivocally point to a dominant Kitaev interaction in the honeycomb-lattice cobaltate BaCo 2 ⁢(AsO 4 ) 2 . Our anisotropic-exchange model closely describes all available neutron scattering data in the material’s field-polarized phase. Furthermore, the density-matrix renormalization group results for our model are in close accord with the unusual double-zigzag magnetic order and the low in-plane saturation field of BaCo 2 ⁢(AsO 4 ) 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Analysis of phase stability and chemical segregation in the Mo-V alloys using a generalized embedded atom method potential

A new interatomic potential for the Mo-V system is introduced to facilitate the study of phase stability and mechanical properties at lower temperatures. This potential is based on a generalization of the embedded atom method and includes contributions from embedding energy, explicit two- and three-body interactions and nonlocal many-body interaction terms. The parameters of the potential are optimized by using data from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, such as elastic constants, equation of states, phonon dispersion curves, point defect properties and melting temperatures for different compositions. Even though our potential is trained on a small dataset, its accuracy is comparable to available machine learning potentials for Mo and V. Furthermore, our results show that an ordered B2 phase is stable at low temperatures in alloys containing 50% V, but the solid solution phase is stable above 800 K. However, such long-range ordering is not observed in V-rich or Mo-rich alloys. In addition, our results show that V segregates to dislocation cores and grain boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optimally tuned starting point for single-shot GW calculations of solids

The dependence of ab initio many-body perturbation theory within the $\textit{GW}$ approximation on the eigensystem used in calculating quasiparticle corrections limits this method's predictive power. Here, we investigate the accuracy of the recently developed Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional as a generalized Kohn-Sham starting point for single-shot $\textit{GW}$ ($G_0W_0$) calculations for a range of semiconductors and insulators. Comparison to calculations based on well-established functionals, namely, PBE, PBE0, and HSE, as well as to self-consistent GW schemes and to experiment, shows that band gaps computed via $G_0W_0$@WOT-SRSH have a level of precision and accuracy that is comparable to that of more advanced methods such as quasiparticle self-consistent $\textit{GW}$ and eigenvalue self-consistent $\textit{GW}$. In this work, we also find that $G_0W_0$@WOT-SRSH improves the description of states deeper in the valence band manifold. Finally, we show that $G_0W_0$@WOT-SRSH significantly reduces the sensitivity of computed band gaps to ambiguities in the underlying WOT-SRSH tuning procedure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tailoring optical response of MXene thin films

Due to their attractive optical properties, 2D MXenes have garnered interest in nanophotonic and optoelectronic applications. However, tuning their properties typically requires the iterative synthesis of MXenes with a specific set of properties, such as the absorption band position, electronic conductivity, and dielectric constant. We demonstrate how to tailor the optical properties of MXene thin films over a broad 1500-nm wavelength range by mixing different ratios of highly conductive Ti 3 C 2 T x with poorly conductive Nb 2 CT x . By changing the MXene film composition, the epsilon-near-zero (ENZ) point, where the optical properties transit from dielectric to metallic, was varied in the spectral range from 1.1 to 2.6 µm. Additionally, we observed a reduction in absorption in some compositions compared to the absorption of the pure MXene films. Compared to other methods, this approach enables simple and continual tuning of MXene optical properties without requiring multiple time-consuming synthesis steps.

2D materials↗
DOE OSTI2024

High pressure melt line of nickel using a generalized embedded atomic method potential

As the second most abundant metal in the Earth's core, nickel plays an important role in determining the structure and temperature of the Earth's core. Yet, the melt line of Ni at pressures corresponding to the Earth's core has not been explored in the literature. Many previous experimental and simulation efforts have reported the melting point of Ni at pressures below 100 GPa, but there exist large discrepancies, most of which have persisted due to various experimental and simulation bottlenecks in handling extreme pressure and temperature conditions. We adopted the generalized embedded atom method, which overcomes the limitations of existing interatomic potentials, to probe phase stability and phase boundaries of Ni at pressures between 50 and 500 GPa. Further, the potential was validated by comparing the cold curves, phonon dispersion curves, and enthalpies of fusion with ab initio density functional theory calculations. Our analysis shows that face centered cubic (FCC) is stable, and the hexagonal close packed (HCP) and body centered cubic (BCC) phases are metastable close to the melt line. Melting temperatures at different pressures were obtained from two-phase co-existence simulations and take the following functional form: $T$ m = $1969.23+19.15P-0.012P$ 2 . In contrast to iron, differences between the melting points of the stable and metastable phases of Ni are less than 250 K at 300 GPa, and the difference in melting points of the metastable BCC and HCP phases changes sign at 500 GPa, which implies that the phase transition mechanisms during solidification can be very complex.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantum spin liquids bootstrapped from Ising criticality in Rydberg arrays

Arrays of Rydberg atoms constitute a highly tunable, strongly interacting venue for the pursuit of exotic states of matter. We develop a strategy for accessing a family of fractionalized phases known as quantum spin liquids in two-dimensional Rydberg arrays. We specifically use effective field theory methods to study arrays assembled from Rydberg chains tuned to an Ising phase transition that famously hosts emergent fermions propagating within each chain. This highly entangled starting point allows us to naturally access spin liquids familiar from Kitaev's honeycomb model—albeit from an entirely different framework. In particular, we argue that finite-range repulsive Rydberg interactions, which frustrate nearby symmetry-breaking orders, can enable coherent propagation of emergent fermions between the chains in which they were born. Delocalization of emergent fermions across the full two-dimensional Rydberg array yields a gapless Z 2 spin liquid with a single massless Dirac cone. Here, the Rydberg occupation numbers exhibit universal power-law correlations that provide a straightforward experimental diagnostic of this phase. We further show that explicitly breaking symmetries perturbs the gapless spin liquid into gapped, topologically ordered descendants: Breaking lattice symmetries generates toric-code topological order, whereas introducing Floquet-mediated chirality generates non-Abelian Ising topological order. In the toric-code phase, we analytically construct microscopic incarnations of non-Abelian defects, which can be created and transported by dynamically controlling the atom positions in the array. Our work suggests that appropriately tuned Rydberg arrays provide a cold-atoms counterpart of solid-state “Kitaev materials” and, more generally, it spotlights a different angle for pursuing experimental platforms for Abelian and non-Abelian fractionalization.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Development of an MC&A toolbox for liquid-fueled molten salt reactors with online reprocessing (Final Report)

A critical barrier to the deployment of MSRs is the absence of a well-defined nuclear material control and accounting (MC&A) approach, a vital prerequisite to meet NRC licensing requirements as well as facilitating future international exports. Liquid-fueled MSR variants (especially those incorporating online refueling or reprocessing) present a unique set of challenges to traditional nuclear material control and accountancy. Unlike solid-fueled light-water reactors, traditional item-counting methods cannot be applied as an accountancy strategy. Rather, MC&A approaches to MSR variants (including both uranium and thorium fueled designs) are more analogous to bulk material handling facilities (e.g., enrichment and reprocessing); yet further complicating matters is the fact that fuel medium is also highly radioactive. Meanwhile, the space of MSRs covers a broad range of design parameters, including thermal and fast spectra designs; operation in actinide breeder or burner modes, choice of actinide fuel used (e.g., 235 U, 232 Th / 233 U, denatured 233 U, etc.), pool or loop-type configuration, and even different salt chemistry. Each of these design choices introduce significant challenges to MC&A approaches within MSR facilities. We propose to bridge this gap for liquid-fueled MSRs by developing a modular, component-based test framework for evaluating viable process monitoring and MC&A techniques specifically suited to liquid-fueled MSR system variants employing online refueling or reprocessing. This test platform will consist of a toolbox of independent process modules representing discrete physical units (such as the reactor core, off-gas processing, decay tanks, and actinide separation units), each with its own self-contained physics responsive to the input mass flow, along with appropriate measurement models that can be coupled to key flow points. These dynamic physical signatures thus afford the ability to test the viability and efficacy of potential accountancy techniques under the full range of reactor operating conditions. As process modules are connected via mass flows, the result is a reconfigurable, generic MSR mass flow model capable of serving as an MC&A test platform for a broad spectrum of possible MSR configurations. The resulting MSR MC&A toolbox will enable robust assessment of accountancy strategies for this unique facility type, including analysis of physical feedbacks arising both from depletion of the fuel over time as well as from potential off-normal events, including those introduced by equipment failures (e.g., a pump failure) as well as by malicious action (i.e., attempts to divert material). The proposed toolbox both addresses a critical needs area for the MPACT analysis toolkit while leveraging existing MPACT-sponsored tools, especially with respect to simulation of measurement and accountancy techniques for advanced fuel cycle facilities. Beyond enabling new analysis capabilities for MSR systems, the design of this toolbox will be to such to maximize compatibility with existing MPACT tools, such to enhance existing facility MC&A analysis capabilities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Atomistic simulations and machine learning of solute grain boundary segregation in Mg alloys at finite temperatures

Understanding solute segregation thermodynamics is the first step in investigating grain boundary (GB) properties, such as strong yttrium (Y) effects on grain growth and texture evolution in micro-scale polycrystalline magnesium (Mg) alloys. To estimate the average GB segregation behavior in low-solute-concentration Mg alloys (e.g., 2 at.% Y), a state-of-the-art spectral approach is applied based on a per-site segregation energy spectrum for Y solute atoms at zero K obtained from molecular statistics (MS) simulations of ~10 4 GB sites in Mg symmetric tilt GBs (STGBs). Although selected MS simulation results are consistent with verification by density functional theory (DFT) calculations, estimates of average segregation tendency based on the zero-K energy spectrum deviate from experimental observations. To resolve this problem, thermodynamic integration (TI) methods based on molecular dynamics (MD) simulations are used to determine the per-site segregation free energies of Y at representative GB sites, which show contributions beyond harmonic approximations can be important for certain GB sites at high temperatures. A surrogate model of per-site segregation free energy is constructed from a small set of TI data points using stacking cross-validation regressors and physics-informed descriptors. This model is applied to predict the Y segregation free energy spectra for thousands of GB sites in Mg STGBs with uncertainty quantification. Finally, the average segregation tendency predicted by the spectral approach based on the free energy spectra agrees well (within the uncertainty range) with experimental observations for micro-scale polycrystalline Mg alloys at typical thermomechanical processing temperatures (500 ~ 800 K), where thermodynamic equilibrium states are likely to be achieved due to fast diffusion.

Atomistic simulations↗
DOE OSTI2020

Geometry-independent microstructure optimization for electron beam powder bed fusion additive manufacturing

The ability to manipulate the formation of microstructure is one of the potential advantages of additive manufacturing. Yet, the additive manufacturing process is riddled with complex interactions between processing parameters, geometry, and spot melting sequence which makes the task of controlling microstructure challenging. This paper proposes a method to solve the microstructure texture optimization problem by leveraging established optimization techniques and increasingly fast and accurate analytical heat transfer models. Respectively, we apply two different optimization algorithms, genetic algorithm and concurrent greedy line search, in order to induce arbitrarily shaped microstructure regions and we use a heat transfer model to evaluate the fitness of solutions in an iterative approach that is far more efficient than pure trial and error. The objective of the optimization algorithm is to control local thermal conditions – hence microstructure formation – during the melting process by manipulating the order and the dwell time of all the points of the infill pattern sequence. We then analyze the optimized sequence to determine important path characteristics that affect the microstructure formation. Through examples, we show that relative order and frequency, as measured by weighted average and variance of spot order, are important factors that influence the final characteristics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Erosion-Corrosion Detection in Carbon Steel Pipe Loops using SRNLs Thickness and Mass Loss Measurement Coupons - 20464

Structural health of waste transfer infrastructure is of utmost importance to DOE and communities surrounding their active sites. A Fitness-for-Service (FFS) program for the Waste Transfer System has been implemented by the US Department of Energy (DOE) for the purposes of evaluating and quantifying the degradation of critical components. Currently, the method of approach is to section and evaluate pipelines that have been removed at end service. In order to augment data collection methods applied to making lifespan estimates of waste transfer system components, Florida International University (FIU) has been investigating new-technology based strategies for quantifying the rate of degradation. Present research considers the option of a coupon system installed on a section of pipe at points of interest that will provide easy access to high resolution measurements of erosion/corrosion rates during operation. Researchers at Savannah River National Laboratory (SRNL) have developed this novel solution for in situ erosion and corrosion measurements in many applicable locations such as tanks, pipes, ducts, and injection points. The coupons are approximately 12.7 mm (1/2 inch) in diameter and 25.4 mm (1 inch) in length and can be placed in the pipe sections using Swagelok compression fittings. The coupon design allows for monitoring of the total degradation resulting from corrosion and erosion. High-resolution Ultrasonic Thickness (UT) testing allows in-situ monitoring of both degradation mechanisms. When removed from the system mass loss resulting from corrosion and erosion may also be evaluated. Experimental testing at FIU implemented the SRNL coupon system in a 6.10 meter (20 foot) pipe loop system comprised of 50.8 mm (2 inch) and 76.2 mm (3 inch) pipe sections in order to replicate expected use cases at the sites. Operational conditions were emulated by pumping abrasive materials through the pipeline at full capacity and collecting measurements from the SRNL coupons. Initial testing was conducted by circulating water followed by a mixture of sand and water for uniform wear rates. Thickness changes and mass loss were measured, and the surface exposed to the process fluid was observed for visual evidence of characteristic wear patterns resulting from an erosive environment. The test matrix consisted of several hours of operation and using fine, medium and coarse grain sands. Test results were analyzed to check for all three parameters (mass loss/gain, thickness change and surface patters) which could provide estimates of the remaining useful life of the pipe sections. Thus, the research presents the findings of this evaluation of the SRNL coupons as a means for conducting high accuracy, cost effective in-situ analysis of pipeline and waste transfer systems for erosion and corrosion susceptibility. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Landau level of fragile topology

Here, we study the Hofstadter butterfly and Landau levels of the twisted bilayer graphene (TBG). We show that the nontrivial fragile topology of the lowest two bands near the charge neutral point makes their Hofstadter butterfly generically connected with higher bands, closing the gap between the first and second conduction (valence) bands at a certain magnetic flux per unit cell. We also develop a momentum space method for calculating the TBG Hofstadter butterfly, from which we identify three phases where the Hofstadter butterflies of the lowest two bands and the higher bands are connected in different ways. We show this leads to a crossing between the $\textit{ν}$ = 4 Landau fan from the charge neutral point and the zero-field band gap at one flux per moiré unit cell, which corresponds to a magnetic field 25$θ^2$ T (twist angle $\textit{θ}$ in degrees). This provides an experimentally testable feature of the fragile topology. In general, we expect it to be a generic feature that the Hofstadter butterfly of the topological bands are connected with the Hofstadter spectra of other bands. We further show the TBG band theory with Zeeman splitting being the most sizable splitting could result in Landau fans at the charge neutral point and half fillings near the magic angle, and we predict their variations under an in-plane magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

From Anderson localization on random regular graphs to many-body localization

Highlights: • Anderson transition from ergodicity to localization on random regular graphs (RRG). • Analytical, pool method, and exact-diagonalization study of correlations on RRG. • Many-body localization (MBL) to ergodicity transition: quantum dots and spin chains. • Dynamical eigenstate correlation functions in RRG and MBL problems. • Anderson localization on RRG as a toy model for MBL. The article reviews the physics of Anderson localization on random regular graphs (RRG) and its connections to many-body localization (MBL) in disordered interacting systems. Properties of eigenstate and energy level correlations in delocalized and localized phases, as well at criticality, are discussed. In the many-body part, models with short-range and power-law interactions are considered, as well as the quantum-dot model representing the limit of the “most long-range” interaction. Central themes – which are common to the RRG and MBL problems – include ergodicity of the delocalized phase, localized character of the critical point, strong finite-size effects, and fractal scaling of eigenstate correlations in the localized phase.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Determination of single-crystal elastic moduli of Li RE F 4 ( RE =Y, Gd, and Tb) by resonant ultrasound spectroscopy

The tetragonal fluoro-scheelite Li RE F 4 compounds (RE = rare earth) have been shown to exhibit a variety of useful optical and magnetic properties. While LiYF 4 has been widely studied, many of the fundamental thermodynamic properties of other members of this family remain unknown. Here, we report the complete elastic tensors (C ij ) of single-crystalline LiYF 4 , LiTbF 4 , and LiGdF 4 using resonant ultrasound spectroscopy and density functional theory (DFT). Here we compare the results for LiYF 4 with prior experimental results using time-of-flight ultrasound methods. This is the first report, however, of the experimental elastic tensors of LiTbF 4 and LiGdF 4 . The present results point to a softening of the elastic moduli of the Li RE F 4 system when Y is replaced by the larger ionic radius of Tb or Gd. Furthermore, we find that just 0.3% doping with Nd on the Y site also leads to a slight softening of the moduli. The variation of the elastic moduli as a function of temperature up to 216°C was also measured. A nearly linear softening of all seven independent elastic moduli was observed with increasing temperature. Phonon dispersions and phonon density of states obtained by DFT support the experimental finding of a significantly higher sound velocity due to lighter Y atoms in LiYF 4 , as compared to heavier LiTbF 4 and LiGdF 4 .

36 MATERIALS SCIENCE↗