Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio data”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Performance analysis of distributed symmetric sparse matrix vector multiplication algorithm for multi-core architectures

In this article, sparse matrix vector multiply (SpMVM) is an important kernel that frequently arises in high performance computing applications. Due to its low arithmetic intensity, several approaches have been proposed in literature to improve its scalability and efficiency in large scale computations. In this paper, our target systems are high end multi-core architectures and we use messaging passing interface + open multiprocessing hybrid programming model for parallelism. We analyze the performance of recently proposed implementation of the distributed symmetric SpMVM, originally developed for large sparse symmetric matrices arising in ab initio nuclear structure calculations. We also study important features of this implementation and compare with previously reported implementations that do not exploit underlying symmetry. Our SpMVM implementations leverage the hybrid paradigm to efficiently overlap expensive communications with computations. Our main comparison criterion is the "CPU core hours" metric, which is the main measure of resource usage on supercomputers. We analyze the effects of topology-aware mapping heuristic using simplified network load model. Furthermore, we have tested the different SpMVM implementations on two large clusters with 3D Torus and Dragonfly topology. Our results show that the distributed SpMVM implementation that exploits matrix symmetry and hides communication yields the best value for the "CPU core hours" metric and significantly reduces data movement overheads.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Coherence Spectroscopy in the Condensed Phase: Insights into Molecular Structure, Environment, and Interactions

The role of coherences, or coherently excited superposition states, in complex condensed-phase systems has been the topic of intense interest and debate for a number of years. In many cases, coherences have been utilized as spectators of ultrafast dynamics or for identifying couplings between electronic states. In rare cases, they have been found to drive excited state dynamics directly. Interestingly though, the utilization of coherences as a tool for high-detail vibronic spectroscopy has largely been overlooked until recently, despite their encoding of key information regarding molecular structure, electronically sensitive vibrational modes, and intermolecular interactions. Furthermore, their detection in the time domain makes for a highly comprehensive spectroscopic technique wherein the phase and dephasing times are extracted in addition to amplitude and intensity, an element not afforded in analogous frequency domain “steady-state” measurements. However, practical limitations arise in disentangling the large number of coherent signals typically accessed in broadband nonlinear spectroscopic experiments, often complicating assignment of the origin and type of coherences generated. Two-dimensional electronic spectroscopy (2DES) affords an avenue by which to disperse and decompose the large number of coherent signals generated in nonlinear experiments, facilitating the assignment of various types of quantum coherences. 2DES takes advantage of the broad bandwidth necessary for achieving the high time resolution desired for ultrafast dynamics and coherence generation by resolving the excitation axis to detect all excitation channels independently. This feature is beneficial for following population dynamics such as electronic energy transfer, and 2DES has become the choice method for such studies. Simultaneously, coherences arise as oscillations at well-defined coordinates across the 2D map often atop those evolving population signals. By isolating the coherent contribution to the 2DES data and Fourier transforming along the population time, a 3D spectral representation of the coherent 2D data is generated, and coherences are then ordered by their oscillation frequency, ν2. Individual coherences can then be selected by their frequency and evaluated via their distinct “2D coherence” spectra, yielding a significantly more distinctive set of spectroscopic signatures over other 1D methodologies and single-point 2DES analysis. Given that coherences of different origin result in unique 2D coherence spectra, these characteristics can be catalogued and compared directly against experiment for prompt assignment, a strategy not afforded by traditional 2DES analysis. In this Account, a structure-driven time-independent spectral model is discussed and employed to compare the 2D fingerprints of various coherences to experimental 2D coherence spectra. The frequency-domain approach can easily integrate ab initio derived vibronic parameters, and its correspondence with experimental coherence spectra of a model compound is demonstrated. Several examples and applications are discussed herein, from 2D Franck–Condon analysis of a model compound, to identifying the signatures of interpigment vibronic coupling in a photosynthetic light-harvesting complex. The 3D spectral approach to 2DES provides remarkable spectroscopic detail, in turn leading to new insights in molecular structure and interactions, which complement the time-resolved dynamics simultaneously recorded. The approach presented herein has the potential to distill down the convoluted set of nonlinear signals appearing in 2D coherent spectra, making the technique more amenable to high-detail vibronic spectroscopy in inherently complex condensed phase systems.

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

High Resolution Photoelectron Spectroscopy of the Acetyl Anion

High-resolution photoelectron spectra of cryogenically cooled acetyl anions (CH 3 CO - ) obtained using slow photoelectron velocity-map imaging are reported. The high resolution of the photoelectron spectrum yields a refined electron affinity of 0.4352 ± 0.0012 eV for the acetyl radical as well as the observation of a new vibronic structure that is assigned based on ab initio calculations. Three vibrational frequencies of the neutral radical are measured to be 1047 ± 3 cm -1 (ν 6 ), 834 ± 2 cm -1 (ν 7 ), and 471 ± 1 cm -1 (ν 8 ). Here this work represents the first experimental measurement of the ν 6 frequency of the neutral. The measured electron affinity is used to calculate a refined value of 1641.35 ± 0.42 kJ mol -1 for the gas-phase acidity of acetaldehyde. Analysis of the photoelectron angular distributions provides insight into the character of the highest occupied molecular orbital of the anion, revealing a molecular orbital with strong d-character. Additionally, details of a new centroiding algorithm based on finite differences, which has the potential to decrease data acquisition times by an order of magnitude at no cost to accuracy, are provided.

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

Low-temperature paddlewheel effect in glassy solid electrolytes

Glasses are promising electrolytes for use in solid-state batteries. Nevertheless, due to their amorphous structure, the mechanisms that underlie their ionic conductivity remain poorly understood. Here, ab initio molecular dynamics is used to characterize migration processes in the prototype glass, 75Li2S–25P2S5. Lithium migration occurs via a mechanism that combines concerted motion of lithium ions with large, quasi-permanent reorientations of PS4 3- anions. This latter effect, known as the ‘paddlewheel’ mechanism, is typically observed in high temperature crystalline polymorphs. In contrast to the behavior of crystalline materials, in the glass paddlewheel dynamics contribute to Lithium-ion mobility at room temperature. Paddlewheel contributions are confirmed by characterizing spatial, temporal, vibrational, and energetic correlations with Lithium motion. Furthermore, the dynamics in the glass differ from those in the stable crystalline analogue, γ-Li3PS4, where anion reorientations are negligible and ion mobility is reduced. These data imply that glasses containing complex anions, and in which covalent network formation is minimized, may exhibit paddlewheel dynamics at low temperature. Consequently, these systems may be fertile ground in the search for new solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Quantum mechanical calculation of the collision-induced absorption spectra of N{sub 2}–N{sub 2} with anisotropic interactions

We present quantum mechanical calculations of the collision-induced absorption spectra of nitrogen molecules, using ab initio dipole moment and potential energy surfaces. Collision-induced spectra are first calculated using the isotropic interaction approximation. Then, we improve upon these results by considering the full anisotropic interaction potential. We also develop the computationally less expensive coupled-states approximation for calculating collision-induced spectra and validate this approximation by comparing the results to numerically exact close-coupling calculations for low energies. Angular localization of the scattering wave functions due to anisotropic interactions affects the line strength at low energies by two orders of magnitude. The effect of anisotropy decreases at higher energy, which validates the isotropic interaction approximation as a high-temperature approximation for calculating collision-induced spectra. Agreement with experimental data is reasonable in the isotropic interaction approximation, and improves when the full anisotropic potential is considered. Calculated absorption coefficients are tabulated for application in atmospheric modeling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A transferable classical force field to describe glyme based lithium solvate ionic liquids

A non-polarizable force field for lithium (Li + ) and bis(trifluoromethanesulfonyl)imide (TFSI – ) ions solvated in diglyme at around 0.2 mol fraction salt concentration was developed based on ab initio molecular dynamics (AIMD) simulations and a modified polymer consistent force field model. A force–torque matching based scheme, in conjunction with a genetic algorithm, was used to determine the Lennard-Jones (LJ) parameters of the ion–ion and ion–solvent interactions. This force field includes a partial charge scaling factor and a scaling factor for the 1–4 interactions. The resulting force field successfully reproduces the radial distribution function of the AIMD simulations and shows better agreement compared to the unmodified force field. The new force field was then used to simulate salt solutions with glymes of increasing chain lengths and different salt concentrations. The comparison of the MD simulations, using the new force field, with experimental data at different salt concentrations and AIMD simulations on equimolar concentrations of the triglyme system demonstrates the transferability of the force field parameters to longer glymes and higher salt concentrations. Furthermore, the force field appears to reproduce the features of the experimental x-ray structure factors, suggesting accuracy beyond the first solvation shell, for equimolar salt solutions using both triglyme and tetraglyme as the solvent. Altogether, the new force field was found to accurately reproduce the molecular descriptions of LiTFSI-glyme systems not only at various salt concentrations but also with glymes of different chain lengths. Thus, the new force field provides a useful and accurate tool to perform in silico studies of this family of systems at the atomistic level.

25 ENERGY STORAGE↗
DOE OSTI2017

Magnetic properties of the CrMnFeCoNi high-entropy alloy

In this paper, we present experimental data showing that the equiatomic CrMnFeCoNi high-entropy alloy undergoes two magnetic transformations at temperatures below 100 K while maintaining its fcc structure down to 3 K. The first transition, paramagnetic to spin glass, was detected at 93 K and the second transition of the ferromagnetic type occurred at 38 K. Field-assisted cooling below 38 K resulted in a systematic vertical shift of the hysteresis curves. Strength and direction of the associated magnetization bias was proportional to the strength and direction of the cooling field and shows a linear dependence with a slope of 0.006 ± 0.001 emu T. The local magnetic moments of individual atoms in the CrMnFeCoNi quinary fcc random solid solution were investigated by ab initio (electronic density functional theory) calculations. Results of the numerical analysis suggest that, irrespective of the initial configuration of local magnetic moments, the magnetic moments associated with Cr atoms align antiferromagnetically with respect to a cumulative magnetic moment of their first coordination shell. The ab initio calculations further showed that the magnetic moments of Fe and Mn atoms remain strong (between 1.5 and 2 μ B ), while the local moments of Ni atoms effectively vanish. Finally, these results indicate that interactions of Mn- and/or Fe-located moments with the surrounding magnetic structure account for the observed macroscopic magnetization bias.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development and transferability of neural-network models for plasma-surface interactions

Plasma-surface interactions are increasingly critical to modern technologies; yet, accurate molecular dynamics simulations remain limited by the capabilities of interatomic potentials. Deep Potentials (DPs) promise to revolutionize the field by providing a systematic method for producing accurate interatomic potentials. The primary challenge of DP development is selecting a dataset, which efficiently spans the set of atomic environments one expects to encounter in the subsequent molecular dynamics simulations. The computational cost of density functional theory calculations, which are the typical basis for DP development, makes it impossible to directly verify the quality of a given DP. To address this challenge, we explore the development of a deep-learned interatomic potential, “DeepREBO,” trained to reproduce the behavior of the REBO2 empirical potential, enabling direct validation of training methodology and transferability. Using an active learning framework, we begin with a minimal dataset and iteratively expand it to train a Deep Potential-Smooth Edition model that faithfully reproduces REBO2 results for 25 eV hydrogen bombardment of diamond (001), a particularly challenging case. We show that small, carefully curated datasets can outperform large, unguided ones, with effective models requiring fewer than 15 000 snapshots. Subsequent transferability tests demonstrate that while DeepREBO generalizes well to diamond (111) surfaces, performance degrades for amorphous carbon or higher-energy impacts, highlighting the need for use-case-specific training data. We also evaluate methods to improve short-range repulsion. This study outlines best practices for training robust deep potentials and underscores the importance of dataset design for predictive plasma simulations.

Ab-initio molecular dynamics↗
DOE OSTI2022

Ab Initio Phase Diagram of Chromium to 2.5 TPa

Chromium possesses remarkable physical properties such as hardness and corrosion resistance. Chromium is also a very important geophysical material as it is assumed that lighter Cr isotopes were dissolved in the Earth’s molten core during the planet’s formation, which makes Cr one of the main constituents of the Earth’s core. Unfortunately, Cr has remained one of the least studied 3d transition metals. In a very recent combined experimental and theoretical study (Anzellini et al., Scientific Reports, 2022), the equation of state and melting curve of chromium were studied to 150 GPa, and it was determined that the ambient body-centered cubic (bcc) phase of crystalline Cr remains stable in the whole pressure range considered. However, the importance of the knowledge of the physical properties of Cr, specifically its phase diagram, necessitates further study of Cr to higher pressure. In this work, using a suite of ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries, we obtain the theoretical phase diagram of Cr to 2.5 TPa. We calculate the melting curves of the two solid phases that are present on its phase diagram, namely, the lower-pressure bcc and the higher-pressure hexagonal close-packed (hcp) ones, and obtain the equation for the bcc-hcp solid–solid phase transition boundary. We also obtain the thermal equations of state of both bcc-Cr and hcp-Cr, which are in excellent agreement with both experimental data and QMD simulations. We argue that 2180 K as the value of the ambient melting point of Cr which is offered by several public web resources (“Wikipedia,” “WebElements,” “It’s Elemental,” etc.) is most likely incorrect and should be replaced with 2135 K, found in most experimental studies as well as in the present theoretical work.

equation of state↗
DOE OSTI2016

Ab initio generation of thermal neutron scattering law for uranium dioxide

Neutronic analysis of moderated nuclear fission systems, e.g. nuclear reactors, relies upon an accurate prediction of the thermal neutron spectrum. Thermal neutrons have wavelengths on the order of the interatomic spacing between atoms in matter. Consequently, the scattering cross-sections of these neutrons are affected by the scattering law of the material (i.e. S(α,β), where α and β are dimensionless momentum and energy transfer variables), which describes the atomic structure and thermodynamic motion of the atomic nuclei. Uranium-dioxide (UO{sub 2}) is the primary fuel source for nuclear power reactors moderated by light water. In such systems (e.g. nuclear power reactors, spent fuel pools, etc.) the nuclear fuel comprises a large fraction of the system volume, approximately 1/3 in light-water reactors. Therefore, the scattering law of the fuel material may impact the thermal neutron spectrum and criticality of these systems. In this work previously established ab initio lattice dynamics (AILD) methods were utilized to calculate the phonon density of states (DOS), which acts as a fundamental input to the scattering law. Subsequently, the S(α,β) of UO{sub 2} was generated for publication as an ENDF thermal scattering library (TSL), which will be available for use in standard neutronics computer applications. (authors)

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

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

Study on impact of Cr and Mo on diffusion of H in 2.25Cr1Mo steel using first-principle calculations

Chrome-molybdenum steel (2.25Cr1Mo steel) is one of the main materials of a steam generator (SG) in high-temperature reactor-pebblebed modules (HTR-PM). It is essential to analyze the source term of tritium in this material, because the behavior of tritium in a SG is important for performing source term analysis in normal and accident conditions. In this article, the diffusion behavior of H atom in 2.25Cr1Mo steel was calculated to estimate the diffusivity of tritium using first-principles density functional theory. To develop and simplify the model of hydrogen diffusion in 2.25Cr1Mo steel, the impact of Cr and Mo on the diffusion of hydrogen in bcc-Fe were first calculated, all the possible diffusion paths were considered, and the minimum energy path was obtained. The diffusion activation energy and pre-exponential factor of the diffusion coefficient were obtained from Vienna Ab initio Simulation Package combined with the climbing image-nudged elastic band method. Here the results indicate that the minimum energy path for the impurity H atom is from one tetrahedral interstitial site to an adjacent tetrahedral interstitial site. The function of the diffusion coefficient of H in 2.25Cr1Mo steel with temperature T can be expressed as D = 1.486 X 10 -7 x ($\frac{-16.350kJ/mol}{RT}$) (M2/s). The diffusion coefficient of our calculation and some of the previous experiments have an excellent quantitative agreement, which indicates the reliability of our crystalline model and the practicability of the present theoretical approach. More importantly, the computational results in this work can be treated as a good screening method to collect reasonable experimental data, which will provide a good reference for tritium source term evaluation in the SG of the HTR-PM.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Rovibrational Spectroscopic Constants and Fundamental Vibrational Frequencies for Isotopologues of Cyclic and Bent Singlet HC2N isomers

Through established, highly-accurate ab initio quartic force fields (QFFs), a complete set of fundamental vibrational frequencies, rotational constants, and rovibrational coupling and centrifugal distortion constants have been determined for both the cyclic 1(sup 1) 1A' and bent 2(sup 1)A' DCCN, H(C13)CCN, HC(C-13)N, and HCC(N-15) isotopologues of HCCN. Spectroscopic constants are computed for all isotopologues using second-order vibrational perturbation theory (VPT2), and the fundamental vibrational frequencies are computed with VPT2 and vibrational configuration interaction (VCI) theory. Agreement between VPT2 and VCI results is quite good with the fundamental vibrational frequencies of the bent isomer isotopologues in accord to within a 0.1 to 3.2 / cm range. Similar accuracies are present for the cyclic isomer isotopologues. The data generated here serve as a reference for astronomical observations of these closed-shell, highly-dipolar molecules using new, high-resolution telescopes and as reference for laboratory studies where isotopic labeling may lead to elucidation of the formation mechanism for the known interstellar molecule: X 3A0 HCCN.

molecular data↗
DOE OSTI2017

Ab initio molecular dynamics simulations of liquid water using high quality meta-GGA functionals

We have used ab initio molecular dynamics (AIMD) to characterize water properties using two meta-generalized gradient approximation (meta-GGA) functionals, M06-L-D3 and B97M-rV, and compared their performance against a standard GGA corrected for dispersion, revPBE-D3, at ambient conditions (298 K, and 1 g cm –3 or 1 atm). Simulations of the equilibrium density, radial distribution functions, self-diffusivity, the infrared spectrum, liquid dipole moments, and characterizations of the hydrogen bond network show that all three functionals have overcome the problem of the early AIMD simulations that erroneously found ambient water to be highly structured, but they differ substantially among themselves in agreement with experiment on this range of water properties. We show directly using water cluster data up through the pentamer that revPBE-D3 benefits from a cancellation of its intrinsic functional error by running classical trajectories, whereas the meta-GGA functionals are demonstrably more accurate and would require the simulation of nuclear quantum effects to realize better agreement with all cluster and condensed phase properties.

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

Origins of Singlet Fission in Solid Pentacene from an ab initio Green’s Function Approach

We develop a new first-principles approach to predict and understand rates of singlet fission with an ab initio Green's-function formalism based on many-body perturbation theory. Starting with singlet and triplet excitons computed from a GW plus Bethe-Salpeter equation approach, we calculate the exciton--bi-exciton coupling to lowest order in the Coulomb interaction, assuming a final state consisting of two non-interacting spin-correlated triplets with finite center-of-mass momentum. For crystalline pentacene, symmetries dictate that the only purely Coulombic fission decay from a bright singlet state requires a final state consisting of two inequivalent nearly degenerate triplets of nonzero, equal and opposite, center-of-mass momenta. For such a process, we predict a singlet lifetime of 40 to 150 fs, in very good agreement with experimental data, indicating that this process can dominate singlet fission in crystalline pentacene. As a result, our approach is general and provides a framework for predicting and understanding multiexciton interactions in solids.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Rovibrational spectroscopic constants and fundamental vibrational frequencies for isotopologues of cyclic and bent singlet HC{sub 2}N isomers

Through established, highly accurate ab initio quartic force fields, a complete set of fundamental vibrational frequencies, rotational constants, and rovibrational coupling and centrifugal distortion constants have been determined for both the cyclic 1 {sup 1} A' and bent 2 {sup 1} A' DCCN, H{sup 13}CCN, HC{sup 13}CN, and HCC{sup 15}N isotopologues of HCCN. Spectroscopic constants are computed for all isotopologues using second-order vibrational perturbation theory (VPT2), and the fundamental vibrational frequencies are computed with VPT2 and vibrational configuration interaction (VCI) theory. Agreement between VPT2 and VCI results is quite good, with the fundamental vibrational frequencies of the bent isomer isotopologues in accord to within a 0.1-3.2 cm{sup –1} range. Similar accuracies are present for the cyclic isomer isotopologues. The data generated here serve as a reference for astronomical observations of these closed-shell, highly dipolar molecules using new, high-resolution telescopes and as reference for laboratory studies where isotopic labeling may lead to elucidation of the formation mechanism for the known interstellar molecule: X {sup 3} A' HCCN.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Towards the extraction of the crystal cell parameters from pair distribution function profiles

The approach based on atomic pair distribution function (PDF) has revolutionized structural investigations by X-ray/electron diffraction of nano or quasi-amorphous materials, opening up the possibility of exploring short-range order. However, the ab initio crystal structural solution by the PDF is far from being achieved due to the difficulty in determining the crystallographic properties of the unit cell. A method for estimating the crystal cell parameters directly from a PDF profile is presented, which is composed of two steps: first, the type of crystal cell is inferred using machine-learning approaches applied to the PDF profile; second, the crystal cell parameters are extracted by means of multivariate analysis combined with vector superposition techniques. The procedure has been validated on a large number of PDF profiles calculated from known crystal structures and on a small number of measured PDF profiles. The lattice determination step has been benchmarked by a comprehensive exploration of different classifiers and different input data. The highest performance is obtained using the k -nearest neighbours classifier applied to whole PDF profiles. Descriptors calculated from the PDF profiles by recurrence quantitative analysis produce results that can be interpreted in terms of PDF properties, and the significance of each descriptor in determining the prediction is evaluated. The cell parameter extraction step depends on the cell metric rather than its type. Monometric, dimetric and trimetric cells have top-1 estimates that are correct 40, 20 and 5% of the time, respectively. Promising results were obtained when analysing real nanocrystals, where unit cells close to the true ones are found within the top-1 ranked solution in the case of monometric cells and within the top-6 ranked solutions in the case of dimetric cells, even in the presence of a crystalline impurity with a weight fraction up to 40%.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structural and electronic characteristics of amorphous Ge 8 Sb 2 Te 11

GeTe-rich phase-change materials can be utilized in rewriteable optical memory due to the large contrast in reflectivity between amorphous and crystalline phases. Here we explored the structure and electronic properties of amorphous Ge 8 Sb 2 Te 11 using ab initio molecular dynamics simulations. The results indicate that amorphous Ge 8 Sb 2 Te 11 is dominantly composed of Ge-, Sb- and Te-centered octahedrons with distortions, while 30.4 % of Ge-centered clusters are in the form of tetrahedrons which are randomly distributed. The 5-fold rings possess a large proportion, and the Gesingle bondTe and Sbsingle bondTe bonds present larger formation energies than other bonds, leading to the ABAB bonding sequence (A: Ge and Sb, B: Te). The lone pair electrons locating at the opposite direction of bonds possess a large fraction of 14.8 %, which may enhance the distortions of local clusters. Importantly, these structural properties lead to the high stability of amorphous Ge 8 Sb 2 Te 11 and thus long data retention in the information storage.

36 MATERIALS SCIENCE↗