Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “average atom”

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 145 records · Page 8

Development of a wide bandwidth heterodyne dispersion interferometer for electron density measurement of atmospheric pressure plasmas

One of the challenges of electron density measurements of an atmospheric pressure plasma (APP) with a laser interferometer is the significant and unwanted phase shift caused by changes in the neutral gas density. These unwanted phase shifts can be mitigated and plasma density measured using an interferometer architecture called a dispersion interferometer (DI). A DI is composed of two nonlinear orientation patterned GaAs crystals for frequency doubling and measures the phase shift induced by a plasma in the interference signal between two second-harmonic beams. Measurement of plasma dynamics or a short plasma pulse in less than a millisecond is enabled with heterodyne detection in a DI using an acousto-optic cell with a frequency of 40 MHz. This heterodyne DI (HDI) is targeted to measure APPs in an electron density range of 10 20 –10 24 m −3 . Finally, the HDI achieves a line-integrated density resolution of 2 × 10 15 m −2 (a phase resolution of 0.005°) with a time constant of 1 μs using ensemble averaging techniques.

atmospheric pressure plasma↗

Quantitative three-dimensional imaging of chemical short-range order via machine learning enhanced atom probe tomography

Abstract Chemical short-range order (CSRO) refers to atoms of specific elements self-organising within a disordered crystalline matrix to form particular atomic neighbourhoods. CSRO is typically characterized indirectly, using volume-averaged or through projection microscopy techniques that fail to capture the three-dimensional atomistic architectures. Here, we present a machine-learning enhanced approach to break the inherent resolution limits of atom probe tomography enabling three-dimensional imaging of multiple CSROs. We showcase our approach by addressing a long-standing question encountered in body-centred-cubic Fe-Al alloys that see anomalous property changes upon heat treatment. We use it to evidence non-statistical B 2 -CSRO instead of the generally-expected D0 3 -CSRO. We introduce quantitative correlations among annealing temperature, CSRO, and nano-hardness and electrical resistivity. Our approach is further validated on modified D0 3 -CSRO detected in Fe-Ga. The proposed strategy can be generally employed to investigate short/medium/long-range ordering phenomena in different materials and help design future high-performance materials.

36 MATERIALS SCIENCE↗

Medium-independent hydrogen atom binding isotherms of nickel oxide electrodes

Adsorption and transfer of hydrogen atoms at solid/solution interfaces are fundamental to heterogeneous catalysis for chemical energy transformations and other processes. Reported here are electrochemical and spectroelectrochemical measurements of the thermodynamics of H-atom binding to nickel oxide electrodes, both the average and the distribution of NiO–H bond dissociation free energies (BDFEs). These are perhaps the first measurements of binding isotherms at non-metal electrodes. Remarkably, both the BDFEs and the non-Langmuirian isotherms are the same in water, acetonitrile, and dimethylformamide, and with different buffers and proton activities. Such medium independence of the BDFEs and isotherms has not been previously reported for any binary material. The medium independence supports the common use of computed hydrogen binding energies as intrinsic descriptors of surface reactivity, while the broadened isotherms add a level of complexity to such analyses. Furthermore, this work demonstrates the capability to derive key thermodynamic parameters at chemically reactive solid-liquid interfaces.

13 HYDRO ENERGY↗

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗

Accurate Prediction of Voltage of Battery Electrode Materials Using Attention-Based Graph Neural Networks

Performing first-principles calculations to discover electrodes’ properties in the large chemical space is a challenging task. While machine learning (ML) has been applied to effectively accelerate those discoveries, most of the applied methods ignore the materials’ spatial information and only use predefined features: based only on chemical compositions. Here, we propose two attention-based graph convolutional neural network techniques to learn the average voltage of electrodes. Our proposed methods, which combine both atomic composition and atomic coordinates in 3D-space, improve the accuracy in voltage prediction significantly when compared to composition-based ML models. The first model directly learns the chemical reaction of electrodes and metal ions to predict their average voltage, whereas the second model combines electrodes’ ML predicted formation energy (E form ) to compute their average voltage. Our E form -based model demonstrates improved accuracy in transferability from our subset of learned Li ions to Na ions. Moreover, we predicted the theoretical voltage of 10 Na x MPO 4 F (M = Ti, Cr, Fe, Cu, Mn, Co, and Ni) fluorophosphate battery frameworks, which are unavailable in the Material Project database. It could be shown that we can expect average voltages higher than 3.1 V from those Na battery frameworks except from the NaTiPO 4 F and TiPO 4 F pair of electrodes, which offer an average voltage of 1.32 V.

25 ENERGY STORAGE↗

Attosecond Time Delay Trends across the Isoelectronic Noble Gas Sequence

The analysis and measurement of Wigner time delays can provide detailed information about the electronic environment within and around atomic and molecular systems, with one the key differences being the lack of a long-range potential after a halogen ion undergoes photoionization. In this work, we use relativistic random-phase approximation to calculate the average Wigner delay from the highest occupied subshells of the atomic pairings (2p, 2s in Fluorine, Neon), (3p, 3s in Chlorine, Argon), (4p, 4s, 3d, in Bromine, Krypton), and (5p, 5s, 4d in Iodine, Xenon). The qualitative behaviors of the Wigner delays between the isoelectronic pairings were found to be similar in nature, with the only large differences occurring at photoelectron energies less than 20 eV and around Cooper minima. Interestingly, the relative shift in Wigner time delays between negatively charged halogens and noble gases decreases as atomic mass increases. All atomic pairings show large differences at low energies, with noble gas atoms showing large positive Wigner delays, while negatively charged halogen ions show negative delays. The implications for photoionization studies in halide-containing molecules is also discussed.

Grafstrom, Brock↗

Observational constraints on the multiphase nature of outflows using large spectroscopic surveys at z ∼ 0

ABSTRACT Mass outflow rates and loading factors are typically used to infer the quenching potential of galactic-scale outflows. However, these generally rely on observations of a single gas phase that can severely underestimate the total ejected gas mass. To address this, we use observations of high mass (≥1010 M⊙), normal star-forming galaxies at z ∼ 0 from the MaNGA, xCOLD GASS, xGASS, and ALFALFA surveys and a stacking of Na d, Hα, CO(1–0), and H i 21 cm tracers with the aim of placing constraints on an average, total mass outflow rate, and loading factor. We find detections of outflows in both neutral and ionized gas tracers, with no detections in stacks of molecular or atomic gas emission. Modelling of the outflow components reveals velocities of |vNa d| = 131 km s−1 and |vHα| = 439 km s−1 and outflow rates of $\dot{M}_{\rm {Na\,\small{D}}}$ = 7.55 M⊙ yr−1 and $\dot{M}_{\text{H}\alpha }$ = 0.10 M⊙ yr−1 for neutral and ionized gas, respectively. Assuming a molecular/atomic outflow velocity of 200 km s−1, we derive upper limits of $\dot{M}_{\text{CO}}\lt 19.43$ M⊙ yr−1 and $\dot{M}_{\rm {H\,\small {I}}}\lt $ 26.72 M⊙ yr−1 for the molecular and atomic gas, respectively. Combining the detections and upper limits, we find average total outflow rates of $\dot{M}_{\text{tot}}\lesssim$27 M⊙ yr−1 and a loading factor of ηtot ≲ 6.39, with molecular gas likely contributing ≲72 per cent of the total mass outflow rate, and neutral and ionized gas contributing ∼28 and <1 per cent, respectively. Our results suggest that, to first order, a degree of quenching via ejective feedback could occur in normal galaxies when considering all gas phases, even in the absence of an active galactic nucleus.

Roberts-Borsani, G. W.↗

Synthesis and Characterization of Pt‐Ag Icosahedral Nanocages with Enhanced Catalytic Activity toward Oxygen Reduction

Abstract There is an urgent need to develop cost‐effective electrocatalysts based on Pt for a broad spectrum of applications, including those vital to the operation of fuel cells. Hollowing out the interior of Pt nanocrystals offers a simple and viable strategy for maximizing the utilization efficiency of this precious metal while enhancing the electrocatalytic performance. Herein, we report the synthesis and electrocatalytic evaluation of Pt−Ag icosahedral nanocages with an average wall thickness of 1.6 nm. The Pt atoms are coated on the surface of Ag icosahedral seeds, leading to the formation of Ag@Pt nL core‐shell icosahedral nanocrystals with tunable shell thicknesses. The core‐shell nanocrystals are then converted to icosahedral nanocages by selectively etching away the Ag in the core. The as‐obtained nanocages with a composition of Pt 4.5 Ag exhibit an almost 3‐fold enhancement in specific activity toward oxygen reduction relative to the commercial Pt/C in acid media.

Wang, Wenxia↗

Codesigning Alloy Compositions of CdSe y Te 1− y Absorbers and Mg x Zn 1− x O Contacts to Increase Solar Cell Efficiency

Thin‐film solar cells such as CdTe are a major commercial photovoltaic technology, with more than 25 GW installed worldwide and levelized costs of electricity competitive with fossil fuels. Further progress may result from integrating CdSe y Te 1− y absorbers with Mg x Zn 1− x O contacts, but the device efficiency is difficult to maximize due to coupled dependence on chemical composition of both alloys. Herein, a high‐throughput approach is demonstrated to codesign chemical compositions in alloyed Mg x Zn 1− x O/CdSe y Te 1− y thin‐film solar cells, using combinatorial libraries of PV devices with orthogonal composition gradients in CdSe y Te 1− y absorbers and Mg x Zn 1− x O contacts. It is found that the solar cell performance is a strong and coupled function of both elemental compositions, with efficiency up to 17.7% ( V OC = 836 mV, fill factor = 69%, J SC = 30.6 mA cm −2 ) at atomic compositions of Mg/(Mg + Zn) ≈18% and average Se/(Se + Te) ≈4%. These performance trends among >100 devices are explained by >100 ns lifetime of photoexcited charge carriers at the Mg x Zn 1− x O/CdSe y Te 1− y interface where strong Se accumulation is also observed. This study reports the optimal compositions of the commercially relevant Mg x Zn 1− x O/CdSe y Te 1− y solar cells and demonstrates a general approach to codesigning performance of alloyed thin‐film solar cells and other optoelectronic devices.

14 SOLAR ENERGY↗

Machine-learning-guided descriptor selection for predicting corrosion resistance in multi-principal element alloys

More than $270 billion is spent on combatting corrosion annually in the USA alone. As such, we present a machine-learning (ML) approach to down select corrosion-resistant alloys. Our focus is on a non-traditional class of alloys called multi-principal element alloys (MPEAs). Given the vast search space due to the variety of compositions and descriptors to be considered, and based upon existing corrosion data for MPEAs, we demonstrate descriptor optimization to predict corrosion resistance of any given MPEA. Our ML model with descriptor optimization predicts the corrosion resistance of a given MPEA in the presence of an aqueous environment by down selecting two environmental descriptors (pH of the medium and halide concentration), one chemical composition descriptor (atomic % of element with minimum reduction potential), and two atomic descriptors (difference in lattice constant (Δa) and average reduction potential). Our findings show that, while it is possible to down select corrosion-resistant MPEAs by using ML from a large search space, a larger dataset and higher quality data are needed to accurately predict the corrosion rate of MPEAs. This study shows both the promise and the perils of ML when applied to a complex chemical phenomenon like corrosion of alloys.

36 MATERIALS SCIENCE↗

Ultrafast lattice disordering can be accelerated by electronic collisional forces

In the prevalent picture of ultrafast structural phase transitions, the atomic motion occurs in a slowly varying potential energy surface determined adiabatically by the fast electrons. However, this ignores non-conservative forces caused by electron-lattice collisions, which can significantly influence atomic motion. Most ultrafast techniques only probe the average structure and are less sensitive to random displacements, and therefore do not detect the role played by non-conservative forces in phase transitions. Here we show that the lattice dynamics of the prototypical insulator-to-metal transition of VO 2 cannot be described by a potential energy alone. We use the sample temperature to control the preexisting lattice disorder before ultrafast photoexcitation across the phase transition and our ultrafast diffuse scattering experiments show that the fluctuations characteristic of the rutile metal develop equally fast (120 fs) at initial temperatures of 100 K and 300 K. This indicates that additional non-conservative forces are responsible for the increased lattice disorder. In conclusion, these results highlight the need for more sophisticated descriptions of ultrafast phenomena beyond the Born-Oppenheimer approximation as well as ultrafast probes of spatial fluctuations beyond the average unit cell measured by diffraction.

36 MATERIALS SCIENCE↗

Physics of Ultra-High Energy Density Relativistic Plasmas from Ordered Nanostructures

This award funded a four-year program at Colorado State University on the physics of ultra-high energy density plasmas generated by focusing femtosecond, >1021 Wcm-2 laser pulses onto ordered nanowire arrays. We used the ALEPH petawatt-class Ti:Sa laser with its high contrast 400 nm second harmonic beamline to irradiate arrays of Ni and TiO2 nanowires made by template assisted electrodeposition and atomic layer deposition over a wide range of average densities. We met two of the three original project goals: radial ion acceleration via doppler shifted K shell emission, and time and spectrally resolved K shell measurements. The main experimental advance was integrating the Lawrence Livermore National Laboratory (LLNL) TREX ultrafast x-ray streak camera with a suite of spherically bent quartz Bragg crystal spectrometers built using Princeton Plasma Physics Laboratory (PPPL) methodology. This is the first TREX-on-bent-quartz setup fielded at a femtosecond, >1021 Wcm-2 facility. It is now a permanent diagnostic at ALEPH and is offered to the LaserNetUS user community. The program produced three peer reviewed journal articles that acknowledge this award. The platform will continue at the upgraded CSU ATLAS facility, where the ALEPH laser is being upgraded to 2 PW.

Hollinger, Reed [Colorado State University]↗

Combining machine-learned and empirical force fields with the parareal algorithm: application to the diffusion of atomistic defects

We numerically investigate an adaptive version of the parareal algorithm in the context of molecular dynamics. This adaptive variant has been originally introduced in [1]. We focus here on test cases of physical interest where the dynamics of the system is modelled by the Langevin equation and is simulated using the molecular dynamics software LAMMPS. In this work, the parareal algorithm uses a family of machine-learning spectral neighbor analysis potentials (SNAP) as fine, reference, potentials and embedded-atom method potentials (EAM) as coarse potentials. We consider a self-interstitial atom in a tungsten lattice and compute the average residence time of the system in metastable states. Our numerical results demonstrate significant computational gains using the adaptive parareal algorithm in comparison to a sequential integration of the Langevin dynamics. We also identify a large regime of numerical parameters for which statistical accuracy is reached without being a consequence of trajectorial accuracy.

36 MATERIALS SCIENCE↗

Electron–Ion Covariance Reveals Geometry-Resolved Inner-Shell Spectra in CS 2 Photodissociation

The chemical shifts of inner-shell atomic orbitals are highly sensitive to chemical bonding and molecular structure. In evolving systems, however, the spectra of distinct molecular species and geometries overlap, obscuring the underlying chemical dynamics. Here we demonstrate the use of electron–ion covariance analysis to combine the structural sensitivity of Coulomb explosion imaging with inner-shell spectroscopy, yielding geometry-specific spectra of transient and product species. We apply this approach to the excited state dynamics of CS 2 probed by ionization above the S 2p edge. Electron–ion covariance with time- and momentum-selected S + and S 2+ ions isolates distinct S 2p photoelectron spectra for ground-state CS 2 , bent photoexcited CS 2 , the CS photoproduct, and bare atomic sulfur ─ species whose spectra overlap strongly in the channel-averaged measurement. Clear chemical shifts are observed in the covariance photoelectron spectrum for each of these species, all of which are consistent with high-level calculations. Here, by extracting the atomic S contribution to the photoelectron spectrum in a finely time-resolved manner we can disentangle this contribution to the overall time-resolved photoelectron spectrum as the photodissociation proceeds. These results demonstrate the promise of electron–ion covariance as a general approach to geometry-resolved inner-shell spectroscopy, opening a route to tracking structural evolution through chemical shifts in complex photoexcited molecules.

Ionization↗

The relation between lignin sequence and its 3D structure

Background: Lignin, the second most abundant biopolymer on earth, plays a major structural role in plants, conferring mechanical strength and regulating water conduction. Understanding the three-dimensional structure of lignin is important for fundamental reasons as well as engineering plants towards lignin valorization. Lignin lacks a specific primary sequence, making its average chemical composition the focus of most recent studies. However, it remains unclear whether the 3D structure of lignin molecules depends on their sequence. Methods: In this work, we performed all-atom molecular dynamics simulation of three S/G-lignin molecules with the same average composition but different sequence. Results: A detailed statistical analysis of the radius of gyration and relative shape anisotropy reveals that the lignin sequence has no statistically significant effect on the global three-dimensional structure. We found however, that homopolymers of C-lignin with the same molecular weight have smaller radii of gyration than S/G-lignin. We attribute this to lower hydroxyl content of C-lignin, which makes it more compact and rigid. Conclusions: The 3D structure of lignin is influenced by the overall content of monomeric units and interunit linkages and not by its precise primary sequence. General Significance: Lignin is assumed to not have a well-defined primary structure. The results presented here demonstrate there are no significant differences in the global 3D structure of lignin molecules with the same average composition but different primary sequence.

59 BASIC BIOLOGICAL SCIENCES↗

Atomic-Scale Investigation of Oxidation at the Black Phosphorus Surface

Black phosphorus (BP) exhibits extraordinary electronic properties that are desirable for a wide variety of electronic and optoelectronic applications. However, applications of BP are hindered by its rapid degradation in ambient conditions. Despite significant advances that have been made in understanding the degradation mechanism, no consensus has yet been reached on how BP oxidation occurs at the atomic scale as experimental studies have been mostly restricted to averaged effects of degradation over a micron- to millimeter-sized region. Here, BP oxidation is investigated using scanning tunneling microscopy/spectroscopy (STM/S). Introducing O 2 gas to the BP surface in ultrahigh vacuum at a pressure of 10 –5 mbar for 1 min creates two new types of defects on the surface. We identify these defects as dangling atomic oxygen and phosphorus multivacancies using density functional theory simulations. In addition to the structural changes to the surface, the electronic structure is also drastically altered by the introduction of oxygen. The 300 meV band gap of BP is lifted due to dosing. This change in the electronic structure is reversible through STM tip manipulation. As a result, these are the first experimental results showing the atomic-scale oxidation of BP, an important step toward understanding the degradation process

2D materials↗

Synthesis of platinum nanoparticles on strontium titanate nanocuboids via surface organometallic grafting for the catalytic hydrogenolysis of plastic waste

Pt/SrTiO 3 (Pt/STO), prepared on small scale by atomic layer deposition (ALD), is a capable heterogeneous catalyst for the selective hydrogenolysis of polyolefins to hydrocarbon oils, providing a promising approach for upcycling plastic waste. However, because deposition by ALD is costly and resource-intensive, a new synthesis of Pt/STO is needed to effectively scale catalyst production and pursue the commercialization of upcycling processes. To that effect, this work details a scalable deposition method for Pt/STO made by surface organometallic chemistry (SOMC) techniques using Pt(II) acetylacetonate or and trimethyl(methylcyclopentadienyl)platinum. The STO support was calcined (550 °C), treated with ozone (200 °C), and finally steamed (200 °C) to afford a clean STO surface populated with only hydroxyl groups. Pt precursors were dissolved in toluene and deposited onto STO. After reduction at 300 °C, the STO support was decorated with 1.0–1.5 nm Pt nanoparticles. The size and loading of these nanoparticles were varied by employing a multi-cycle deposition and oxidation and/or reduction process designed to ALD techniques. In conclusion, these Pt/STO catalysts hydrogenolyzed isotactic polypropylene into liquid products (>95% yield) with average molecular weights of 200–300 Da (~25 carbon atoms) and narrow size distributions at 300 °C and 180 psi H 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic crystallography reveals spontaneous anisotropy in cubic GeTe

Cubic energy materials such as thermoelectrics or hybrid perovskite materials are often understood to be highly disordered. In GeTe and related IV–VI compounds, this is thought to provide the low thermal conductivities needed for thermoelectric applications. Since conventional crystallography cannot distinguish between static disorder and atomic motions, we develop the energy-resolved variable-shutter pair distribution function technique. This collects structural snapshots with varying exposure times, on timescales relevant for atomic motions. In disagreement with previous interpretations, we find the time-averaged structure of GeTe to be crystalline at all temperatures, but with anisotropic anharmonic dynamics at higher temperatures that resemble static disorder at fast shutter speeds, with correlated ferroelectric fluctuations along the <100> c direction. We show that this anisotropy naturally emerges from a Ginzburg–Landau model that couples polarization fluctuations through long-range elastic interactions. By accessing time-dependent atomic correlations in energy materials, we resolve the long-standing disagreement between local and average structure probes and show that spontaneous anisotropy is ubiquitous in cubic IV–VI materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗