Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic bonding”

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 73 records · Page 4

Plutonium oxide melt structure and covalency

Advances in nuclear power reactors include the use of mixed oxide fuel, containing uranium and plutonium oxides. The high-temperature behaviour and structure of PuO 2-x above 1,800 K remain largely unexplored, and these conditions must be considered for reactor design and planning for the mitigation of severe accidents. Here, in this study, we measure the atomic structure of PuO 2-x through the melting transition up to 3,000 +/- 50 K using X-ray scattering of aerodynamically levitated and laser-beam-heated samples, with O/Pu ranging from 1.57 to 1.76. Liquid structural models consistent with the X-ray data are developed using machine-learned interatomic potentials and density functional theory. Molten PuO 1.76 contains some degree of covalent Pu-O bonding, signalled by the degeneracy of Pu 5f and O 2p orbitals. The liquid is isomorphous with molten CeO 1.75 , demonstrating the latter as a non-radioactive, non-toxic, structural surrogate when differences in the oxidation potentials of Pu and Ce are accounted for. These characterizations provide essential constraints for modelling pertinent to reactor safety design. The molten structure of plutonium oxide-a component of mixed oxide nuclear fuels-is measured, showing some degree of covalent bonding. Its atomic structure is similar to that of cerium oxide, which could be a non-radioactive structural surrogate.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Growth of Mesoscale Ordered Two-Dimensional Hydrogen-Bond Organic Framework with the Observation of Flat Band

Flat bands (FBs), presenting a strongly interacting quantum system, have drawn increasing interest recently. However, experimental growth and synthesis of FB materials have been challenging and have remained elusive for the ideal form of monolayer materials where the FB arises from destructive quantum interference as predicted in 2D lattice models. Here, we report surface growth of a self-assembled monolayer of 2D hydrogen-bond (H-bond) organic frameworks (HOFs) of 1,3,5-tris(4-hydroxyphenyl)benzene (THPB) on Au(111) substrate and the observation of FB. High-resolution scanning tunneling microscopy or spectroscopy shows mesoscale, highly ordered, and uniform THPB HOF domains, while angle-resolved photoemission spectroscopy highlights a FB over the whole Brillouin zone. Density-functional-theory calculations and analyses reveal that the observed topological FB arises from a hidden electronic breathing-kagome lattice without atomically breathing bonds. In conclusion, our findings demonstrate that self-assembly of HOFs provides a viable approach for synthesis of 2D organic topological materials, paving the way to explore many-body quantum states of topological FBs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ultrafast Optomechanical Strain in Layered GeS

Strong coupling between light and mechanical strain forms the foundation for next-generation optical micro- and nano-electromechanical systems. Such optomechanical responses in two-dimensional materials present novel types of functionalities arising from the weak van der Waals bond between atomic layers. In this paper, by using structure-sensitive megaelectronvolt ultrafast electron diffraction, we report the experimental observation of optically driven ultrafast in-plane strain in the layered group IV monochalcogenide germanium sulfide (GeS). Surprisingly, the photoinduced structural deformation exhibits strain amplitudes of order 0.1% with a 10 ps fast response time and a significant in-plane anisotropy between zigzag and armchair crystallographic directions. Rather than arising due to heating, experimental and theoretical investigations suggest deformation potentials caused by electronic density redistribution and converse piezoelectric effects generated by photoinduced electric fields are the dominant contributors to the observed dynamic anisotropic strains. Our observations define new avenues for ultrafast optomechanical control and strain engineering within functional devices.

2D materials↗

Hydrogen separation with a graphenylene monolayer: Diffusion Monte Carlo study

Here, we performed fixed-node diffusion Monte Carlo (DMC) calculations to investigate structural and energetic properties of graphenylene (GPNL), a two-dimensional network of sp 2 -bonded carbon atoms with large near-circular pores, and its H 2 separation performance for gas mixtures. We have found that the energetic stability of a GPNL monolayer is comparable to that of γ-graphyne, as evidenced by its large cohesive energy of 6.755(3) eV/atom. Diffusion barriers of several gas molecules, including hydrogen, through a GPNL membrane were determined from the analysis of their adsorption energies depending on the adsorption distance, which led to our estimation for hydrogen selectivity with respect to other target molecules. DMC hydrogen selectivity of a GPNL monolayer was found to be exceptionally high at 300 K, as high as 10 10 –10 11 against CO and N 2 gases. This, along with high hydrogen permeance due to its generic pore structure, leads us to conclude that GPNL is a promising membrane to be used as a high-performance hydrogen separator from gas mixtures. We find that when compared to our DMC results, DFT calculations tend to overestimate H 2 selectivity, which is mostly due to their inaccurate description of short-range repulsive interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Large Angle Rocking Beam Electron Diffraction Utilizing Electron Direct Detector

Electron diffraction of a crystal is fundamentally a function of the potential of that crystal. The intensity of electron diffraction patterns is, as a result, sensitive to the charge density of atoms and bonding inside crystals. Experimentally, the traditional method to probe this information is quantitative Convergent Beam Electron Diffraction (CBED). Quantitative CBED or QCEBD is a method which uses dynamic diffraction theory to quantify CBED intensities and to extract information about crystal structure and bonding. The sensitivity of dynamical scattering is leveraged to measure crystal symmetry and crystal structure factors. At its limit, electron structure factors are measured at high accuracy allowing access to chemical bonding information, which corroborate theoretical calculations through multipole model refinements of the experimental charge density. The primary limit to QCBED is the requirement of a non-overlapping convergent beam. This subsequently limits QCBED to crystals with small unit cells which are stable under the focused probe.

Busch, Robert↗

Artificial intelligence “sees” split electrons

Chemical bonds between atoms are stabilized by the exchange-correlation (xc) energy, a quantum-mechanical effect in which “social distancing” by electrons lowers their electrostatic repulsion energy. Kohn-Sham density functional theory (DFT) states that the electron density determines this xc energy, but the density functional must be approximated. Furthermore, this is usually done by satisfying exact constraints of the exact functional (making the approximation predictive), by fitting to data (making it interpolative), or both. Two exact constraints—the ensemble-based piecewise linear variation of the total energy with respect to fractional electron number and fractional electron z-component of spin —require hard-to-control nonlocality. On page 1385 of this issue, Kirkpatrick et al. have taken a big step toward more accurate predictions for chemistry through the machine learning of molecular data plus the fractional charge and spin constraints, expressed as data that a machine can learn.

74 ATOMIC AND MOLECULAR PHYSICS↗

graphenv: a Python library for reinforcement learning on graph search spaces

Many important and challenging problems in combinatorial optimization (CO) can be expressed as graph search problems, in which graph vertices represent full or partial solutions and edges represent decisions that connect them. Graph structure not only introduces strong relational inductive biases for learning (Battaglia et al., 2018) - in this context, by providing a way to explicitly model the value of transitioning (along edges) between one search state (vertex) and the next - but lends itself to problems both with and without clearly defined algebraic structure. For example, classic CO problems on graphs such as the Traveling Salesman Problem (TSP) can be expressed as either pure graph search or integer programs. Other problems, however, such as molecular optimization, do no have concise algebraic formulations and yet are readily implemented as a graph search (V. et al., 2022; Zhou et al., 2019). Such "model-free" problems constitute a large fraction of modern reinforcement learning (RL) research owing to the fact that it is often much easier to write a forward simulation that expresses all of the state transitions and rewards, than to write down the precise mathematical expression of the full optimization problem. In the case of molecular optimization, for example, one can use domain knowledge alongside existing software libraries to model the effect of adding a single bond or atom to an existing but incomplete molecule, and let the RL algorithm build a model of how good a given decision is by "experiencing" the simulated environment many times through. In contrast, a model-based mathematical formulation that fully expresses all the chemical and physical constraints is intractable. In recent years, RL has emerged as an effective paradigm for optimizing searches over graphs and led to state-of-the-art heuristics for games like Go and chess, as well as for classical CO problems such as the TSP. This combination of graph search and RL, while powerful, requires non-trivial software to execute, especially when combining advanced state representations such as Graph Neural Networks (GNN) with scalable RL algorithms.

97 MATHEMATICS AND COMPUTING↗

Atomistic insight into the ferroelastic post-stishovite transition by high-pressure single-crystal X-ray diffraction

Abstract The post-stishovite transition is a classic pseudo-proper typed ferroelastic transition with a symmetry-breaking spontaneous strain. This transition has been studied using high-pressure spontaneous strains, optic modes, and elastic moduli (Cij) based on the Landau modeling, but its atomistic information and structural distortion remain poorly understood. Here we have conducted synchrotron single-crystal X-ray diffraction measurements on stishovite crystals up to 75.3 GPa in a diamond-anvil cell. Analysis of the data reveals atomic positions, bond lengths, bond angles, and variations of SiO6 octahedra across the transition at high pressure. Our results show that the O coordinates split at ~51.4 GPa, where the apical and equatorial Si-O bond lengths cross over, the SiO6 octahedral distortion vanishes, and the SiO6 octahedra start to rotate about the c axis. Moreover, distortion mode analysis shows that an in-plane stretching distortion (GM1+ mode) occurs in the stishovite structure at high pressure while a rotational distortion (GM2+ mode) becomes dominant in the post-stishovite structure. These results are used to correlate with elastic moduli and Landau parameters (symmetry-breaking strain e1–e2 and order parameter Q) to provide atomistic insight into the ferroelastic transition. When the bond lengths of two Si-O bonds are equal due to the contribution from the GM1+ stretching mode, C11 converges with C12, and the shear wave VS1[110] polarizing along [110] and propagating along [110] vanishes. Values of e1–e2 and Q are proportional to the SiO6 rotation angle from the occurrence of the GM1+ rotational mode in the post-stishovite structure. Our results on the pseudo-proper type transition are also compared with that for the proper type in albite and improper type in CaSiO3 perovskite. The symmetry-breaking strain, in all these types of transitions, arises as the primary effect from the structural angle (such as SiO6 rotation or lattice constant angle) and its relevant distortion mode in the low-symmetry ferroelastic phase.

Geochemistry & Geophysics↗

Fostering a Guiding Multiscale Model for the Development of Advanced MgB 2 Hydrogen Storage Materials (Final Technical Report)

Project Goal and Objective. The demand for energy and for an upgraded energy infrastructure has steadily grown, as have the needs for energy independence and alternatives to our reliance on petroleum. Hydrogen is considered the most viable fuels for wide-scale implementation in the near future as it is less-polluting, non-toxic, and has more stored energy than petroleum. It is envisioned that hydrogen can eventually become the prime energy carrier, integrating the transportation, grid, and chemical sectors in a way that improves resiliency, diversifies feedstocks, and affords new economic opportunities. A key remaining challenge is the development materials with enhanced gravimetric and volumetric hydrogen storage capacities that offer a higher performance than compressed gas. These materials would eliminate the need for large-scale compression, thereby dramatically reducing the footprint and cost of gas storage. The high gravimetric and volumetric hydrogen capacities of complex hydrides has prompted an intensive investigation of the potential of this class of materials as hydrogen storage media over the past 25 years. Among the many complex hydrides that have been explored, magnesium borohydride, Mg(BH 4 ) 2 , has been found to possess the best combination of practical thermodynamic properties. These include a gravimetric H 2 density of 14.9 wt% H 2 and thermodynamics for the dehydrogenation of Mg(BH 4 ) 2 to MgB 2 (equation 1) (ΔH° = 39 kJ/mol H 2 , ΔS = 112 J/K mol H 2 ) which lie in the narrow window required Mg(BH 4 ) 2 $\Leftrightarrow$ MgB 2 + 4 H 2 (1) for reversibility under moderate pressure and temperature. However, overcoming the extremely slow kinetics of the reversible release of hydrogen by this material in the solid state is a daunting challenge. At temperatures greater than 400 °C, the borohydride releases up to 14 wt% hydrogen giving MgB 2 . We discovered that the direct re-hydrogenation of MgB 2 to Mg(BH 4 ) 2 can be accomplished under 950 bar H 2 at 400 °C. While this demonstrated that complete reversibility can be achieved, the conditions employed are far too extreme for commercial hydrogen storage applications. More recently, we found through US DOE funded research projects (EERE HyMARC and HySCOR), that hydrogen cycling, can be accomplish at much milder conditions upon modification of the borohydride or boride. Guided by these discoveries these discoveries, the objective of this research project was to obtain key information that will enable the development of a model of reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 . The ultimate goal of our efforts is to attain a model of this transformation that can be utilized to accelerate development further advanced materials. This project directly follows on discoveries that were made over the course of a US DOE, EERE HyMARC project that was focused on improvement of the hydrogen cycling kinetics of modified MgB 2 . We found that that mechanical milling with graphene results the desired, pronounced kinetic enhancement. The dramatic lowering of the conditions required for the hydrogenation of MgB 2 is a significant step towards overcoming its chemical inertness allowing its development as a practical onboard hydrogen storage material. However, the exact nature of the modification(s) of MgB 2 that is responsible for its activation towards hydrogenation is completely unknown. This situation is not unique, as efforts to develop hydrogen storage materials typically have a narrow focus rather than a comprehensive approach that takes atomic level bonding and structure; molecular dynamics; long range, nano- and mesoscale-structure and their interconnection all into account. The goal of this project was the development of a comprehensive, multi-scale computational model of reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 that can be utilized for development of higher performance versions of the modified material. Development of the model requires determination of: 1) the bulk, nano-scale, and meso-scale structural changes occurring at elevated pressure following mechano-chemical modification of MgB 2 ; 2) the reaction pathway of the reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 ; 3) the effect of elevated pressure and mechano-chemical modification on the chemical reaction pathways; 4) the interactions at solid-gas interfaces; and particle surfaces; and 5) the kinetics and thermodynamic parameters associated with each step of the hydrogenation reaction pathway. This investigation required advanced techniques as preliminary, standard XRD, 11 B NMR, and FTIR analysis showed no signs of material modification. In order to gain this level of understanding of modified MgB 2 , required the teaming of a diverse group of experts and state-of-the art experimental capabilities at the University of Hawaii at Manoa (UHM) and collaborating National Laboratories: Craig Jensen , Department of Chemistry (PI and Project Director), solid state, solution, and high pressure NMR spectroscopy; solid-state synthesis; and high pressure hydrogenation (collaboration with SNL); Godwin Severa , Hawaii Natural Energy Institute (co-PI) calorimetry; infrared and Raman spectroscopy (collaboration with NREL); Dera , high pressure X-ray diffraction including in situ experiments (collaboration with ANL); Hope Ishii , Hawaii Institute of Geophysics electron microscopy investigations (collaboration with LBNL); and Joe Brown , Mechanical Engineering , material electronic structure and electric field effects.

08 HYDROGEN↗

Diamond-like carbon coating for passive and active electronics

Systems and methods for building passive and active electronics with diamond-like carbon (DLC) coatings are provided herein. DLC may be layered upon substrates to form various components of electronic devices. Passive components such as resistors, capacitors, and inductors may be built using DLC as a dielectric or as an insulating layer. Active components such as diodes and transistors may be built with the DLC acting substantially like a semiconductor. The amount of sp 2 and sp 3 bonded carbon atoms may be varied to modify the properties of the DLC for various electronic components.

Timpson, Erik Joseph↗

Binuclear Acetoniminato Derivatives of Iron and Chromium Carbonyls: A Theoretical Study

The structures and energetics of binuclear acetoniminato metal carbonyl derivatives of the types (Me 2 C═N) 2 Fe 2 (CO) n and (Me 2 C═N) 2 Cr 2 (CO) n have been investigated by density functional theory. The (Me 2 C═N) 2 Fe 2 (CO) 6 structure, with bridging Me 2 C═N ligands using only their nitrogen atoms to bridge an Fe═Fe bond related to experimentally known species, as well as a related (Me 2 C═N) 2 Cr 2 (CO) 8 structure, appear to be energetically preferred. Among carbonyl-richer systems, the iron system energetically prefers a {Me 2 C═NC(O)} 2 Fe 2 (CO) 6 structure with bridging Me 2 C═NC(O) ligands, whereas the chromium system energetically prefers a structure with a 2,3-diazabutadiene ligand having each nitrogen atom bonded to a Cr(CO) 5 moiety. Preferred structures for systems with fewer CO groups include the tetracarbonyl (Me 2 C═N) 2 Fe 2 (CO) 4 in which one of the bridging Me 2 C═N uses its C═N double bond to bond to an iron atom, as well as (Me 2 C═N)2Cr 2 (CO) n (n = 7, 6) structures with chromium–chromium multiple bonds bridged by Me2C═N groups.

Carbonyls↗

Structural and thermal properties and the origin of the ultralow thermal conductivity in the defect stannite CuIn 2 Se 4

Phase-pure CuIn 2 Se 4 , a ternary metal chalcogenide that forms in a disordered stannite crystal structure, was synthesized to investigate the structure thermal property relationships as well as reveal the origin of the ultralow thermal conductivity this material possesses over a large temperature range. Modeling of the temperature-dependent heat capacity and thermal conductivity revealed distinctive thermal properties and large lattice anharmonicity. Electron localization function calculations highlight the asymmetric bonding inherent to CuIn 2 Se 4 , which together with lattice anharmonicity directly impacts the thermal properties. Our findings reveal the specific atomic arrangement and bonding governing the thermal properties of this ternary metal chalcogenide. Our findings underscore the specific atomic arrangement and bonding governing the thermal properties in this ternary chalcogenide. This study advances the fundamental understanding of stannites, and our findings can be applied to these and other multinary metal chalcogenides of interest for applications where low thermal conductivity is desirable.

36 MATERIALS SCIENCE↗

Atomic Structure, Dynamics, Changes in Chemical Bonding and Semiconductor-Metal Transition in Sb 2 Se 3 : A Remarkable Material for Quantum Networks and Energy Applications

Antimony sesquiselenide has become an outstanding functional material for photovoltaics, energy storage and transformation, memory and photonic applications. Sb 2 Se 3 is one of the most successful emerging solar light absorbers and has also been identified as a highly promising ultralow-loss phase-change material (PCM) for next-generation coherent nanophotonic processors, photonic tensor cores, quantum and neuromorphic networks. Unlike benchmark telluride PCMs, Sb 2 Se 3 features a quasi-one-dimensional (1D) crystalline structure consisting of (Sb 4 Se 6 ) ∞ ribbons, lacks the typical PCM chemical bonding, and undergoes an extended semiconductor-metal transition above the melting point. Consequently, the origin of high optical contrast between crystalline (SET) and amorphous (RESET) logic states remains elusive and presents a significant challenge. Using high-energy X-ray diffraction and Raman spectroscopy over a wide temperature range, supported by first-principles simulations and complemented by thermal, optical and electrical measurements, as well as by 121 Sb-Mossbauer spectroscopy, the quasi-1D network of orthorhombic antimony sesquiselenide was found to undergo significant evolution in amorphous and supercooled Sb 2 Se 3 , leading to lower coordination, shorter interatomic distances and a higher p-electron density on antimony, indicating changes in chemical bonding. The observed novel Sb 2 Se 3 nanocrystalline polymorph, characterized by trigonal antimony coordination and more isolated Sb-Se ribbons, could help reduce multiple trapping defect states in the bandgap, which are typical of orthorhombic Sb 2 Se 3 , thereby enhancing the power-conversion efficiency of photovoltaic devices. Semimetallic and metallic liquid Sb 2 Se 3 exhibit a gradual transformation into a denser 2D and/or 3D network with higher antimony coordination. Localized electron states in the pseudogap are becoming extended, leading to an increase in electronic conductivity σ following the relationship σ ∝ N(E F ) 2 . Liquid Sb 2 Se 3 also appears to be strongly fragile, with a nonmonotonic change in viscosity and higher atomic mobility in the metallic liquid. Furthermore, these results explain extraordinary functionalities of Sb 2 Se 3 for photonic and energy applications.

antimony↗

Ultrafast Correlation Energy Estimator

A virtually no-cost method is proposed that can compute the correlation energies of general, covalently bonded, organic, and inorganic molecules (including conjugated π-electron systems) with a well-defined dominant Lewis structure at the accuracy of 99.5% of the near-exact values determined by the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] in the complete-basis-set (CBS) limit. This Correlation Energy Per Bond (CEPB) method assigns a partial correlation energy to each bond type (characterized by the identities of the two atoms forming the bond and its integer bond order) and to a lone pair, regardless of the bond length, bond angle, sp-hybridization, π-electron conjugation, ionicity, noncovalent interactions, etc. At its current stage, the method is mainly suitable for near-equilibrium geometries. The correlation energies per bond are determined by a fit to the CCSD(T)/CBS benchmarks. It can neither improve the equilibrium structures nor discern conformers or positional isomers, yet its accuracy for reaction energies rivals that of the second-order Møller–Plesset perturbation theory, which is far more expensive. Its promising performance underscores the possibility that surprisingly compact, chemically intuitive molecular fragments exist into which correlation energies can be partitioned, leading to various ultrafast correlation-energy estimators tailored to different purposes.

Binding energy↗

Weakening the N–H Bonds of NH 3 Ligands: Triple Hydrogen-Atom Abstraction to Form a Chromium(V) Nitride

Weakening and cleaving N-H bonds is crucial for improving molecular ammonia oxidation catalysts. In this work, we report the synthesis and hydrogen atom abstraction chemistry of bis(ammonia) Cr porphyrin complexes Cr(TPP)(NH 3 ) 2 and Cr(TMP)(NH 3 ) 2 (TPP = 5,10,15,20-tetraphenyl-meso-porphyrin and TMP = 5,10,15,20-tetramesityl-meso-porphyrin) using bulky aryloxyl radicals. The triple hydrogen atom abstraction reaction results in the formation of Cr V (por)(≡N), with the nitride deriving from ammonia, as indicated by UV-Vis and FT-IR spectroscopy, and a single crystal structural determination of Cr(TPP)(≡N). Subsequent oxidation of this chromium(V)-nitrido complex results in the formation of Cr III (por), with scission of the Cr≡N bond. A computational analysis illustrates the progression from Cr II to Cr V and evaluates the energetics of abstracting hydrogen atoms from Cr II -NH 3 to generate Cr V ≡N. The formation and isolation of Cr V (por)(≡N) illustrates the stability of these species and the need to chemically activate the nitride moiety for atom transfer or N-N coupling reactivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bond Dissociation Energies of the Actinide Halides AnX, An = Ac–Lr and X = F–I, Utilizing Relativistic Composite Coupled Cluster Approaches

Bond dissociation energies (BDEs) have been calculated for the set of actinide halides AnX with An=Ac, Pa, and Np-Lr and X=F-I. Two composite thermochemistry methods based on the Feller-Peterson-Dixon (FPD) approach have been utilized, one involving spinor-based relativistic CCSD(T) calculations where spin-orbit (SO) was included at the orbital level and another using scalar relativistic CCSD(T) with a posteriori SO contributions based on 2-component multireference configuration interaction calculations. The method that was chosen for a given actinide halide was based on which representation yielded the best single determinant reference determinant for the coupled cluster calculation. The spinor-based method was chosen for all cases except for AmX, CmX, and BkX. Both composite approaches included contributions accounting for basis set truncation, outer-core correlation, the Gaunt interaction, and QED. The scalar FPD results, as well as the spinor-based calculations for AcF, also included higher order electron correlation up through CCSDT(Q). In addition to BDEs, CCSD(T) equilibrium bond lengths, harmonic frequencies, and vibrational anharmonicity constants are reported for all species. Last, the FPD BDEs for the fluorides were used to confirm the trend across the actinide series previously predicted by Gibson using bonding models based atomic promotion energies that provide a single 6d electron for bonding. In particular the local minimum in the BDEs at AmF is confirmed in the present calculations. Furthermore, the BDEs for LrX are predicted to be slightly larger than those of AcX, making them the largest in the actinide halide series.

Actinides↗

How arsenic makes amorphous GeSe a robust chalcogenide glass for advanced memory integration

The 3D integration technology in semiconductor fabrication requires a key component, the ovonic threshold switching (OTS) selector, to suppress the current leakage. The As doped amorphous (a-) GeSe glass is a commercialized OTS material in 3D phase-change memory, but the understanding of such a doping mechanism is still inadequate. Here we systematically explore the effect of As doping on the structural, bonding, and dynamics properties of a-GeAsSe using ab initio molecular dynamics simulations. The results reveal that As atoms form strong bonds with both Ge and Se atoms. The distorted octahedral structures and the 5-fold rings linked by atoms are increased. All of these structural features lead to a more disordered configuration. Moreover, as atoms have notably slowed down the atomic mobility, rendering a-GeAsSe a high stability. Overall, our studies offer insightful understanding of As-doping in OTS materials, paving the way for the design and application of advanced selector devices.

36 MATERIALS SCIENCE↗