Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic simulation”

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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Visualizing Anisotropic Oxygen Diffusion in Ceria under Activated Conditions

Oxygen reactivity plays a key role in the performance of ceria-based catalysts. Aberration-corrected transmission electron microscopy and molecular dynamics simulations were used to study the oxygen atom diffusion in ceria under activated conditions. Reactive oxygen atom and its real-time diffusion were visualized. Additionally, the interplay between cerium and oxygen atoms originating from a Coulomb interaction was revealed by the out-of-plane buckling of cerium atoms associated with oxygen transport. Anisotropic oxygen atom diffusion that depends on crystal orientations was discovered, demonstrating a preferential [001] crystallographic diffusion pathway. These findings reveal prospects for applications of anisotropic orientation-relevant fluorite-structured oxides.

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

Bulk viscosity of the rigid rotor one-component plasma

Bulk viscosity of a plasma consisting of strongly coupled diatomic ions is computed using molecular dynamics simulations. The simulations are based on the rigid rotor one-component plasma, which is introduced as a model system that adds two degrees of molecular rotation to the traditional one-component plasma. It is characterized by two parameters: the Coulomb coupling parameter, Γ, and the bond length parameter, Ω. Results show that the long-range nature of the Coulomb potential can lead to long rotational relaxation times, which in turn yield large values for bulk viscosity. Here, the bulk-to-shear viscosity ratio is found to span from small to large values depending on the values of Γ and Ω. Although bulk viscosity is often neglected in plasma modeling, these results motivate that it can be large in molecular plasmas with rotational degrees of freedom.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the role of methyl groups in the molecular architectures of mesophase pitches

The role of methyl groups on the liquid–crystal structure of mesophase pitches was investigated by combining experimental characterizations and atomic-scale computational modeling, using three pitches synthesized from different precursors. Of the three pitches, C-9 alkyl benzene and naphthalene-based pitches have 13 and 7 methyl groups per average polyaromatic hydrocarbon, respectively. By contrast, mesophase produced from a coal-tar pitch has about one methyl group. The coal tar–based mesophase pitch is hydrogen deficient or more aromatic compared with C-9 alkyl benzene- and naphthalene-based pitches. Additionally, X-ray diffraction data showed that average coherent domain sizes of C-9 alkyl benzene (3.7 nm) and naphthalene-based (3.6 nm) pitches with more methyl groups are larger than that of coal tar–based mesophase (2.4 nm). Based on the identified features, the influence of the methyl group on the layering structures was investigated via molecular dynamics simulations. The results revealed that methyl groups are critical in mesophase layering in C-9 alkyl benzene- and naphthalene-based pitches, by reducing CH-π interaction. However, similar alignment could be achieved without the same degree of methyl substitutions for the coal tar-based pitch because of stronger π-π interaction than the other precursors. The insights from this study contribute to our understanding of the formation of conventional mesophase pitch and have implications for the processing of coal-derived materials. In conclusion, this knowledge is vital to produce valuable products like carbon fiber and graphite from pitches.

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

Exploring the relation between transonic dislocation glide and stacking fault width in FCC metals

Theory predicts limiting gliding velocities that dislocations cannot overcome. Computational and recent experiments have shown that these limiting velocities are soft barriers and dislocations can reach transonic speeds in high rate plastic deformation scenarios. In this paper we systematically examine the mobility of edge and screw dislocations in several face centered cubic (FCC) metals (Al, Au, Pt, and Ni) in the extreme large-applied-stress regime using molecular dynamics simulations. Our results show that edge dislocations are more likely to move at transonic velocities due to their high mobility and lower limiting velocity than screw dislocations. Importantly, among the considered FCC metals, the dislocation core structure determines the dislocation’s ability to reach transonic velocities. This is likely due to the variation in stacking fault width due to relativistic effects near the limiting velocities.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Defect-Free Growth of Decagonal Quasicrystals around Obstacles

Quasicrystals are aperiodic solids with exotic physical properties. The mechanisms driving their growth are far from understood, particularly in the presence of rigid obstacles. Here, using in situ synchrotron x-ray tomography and molecular dynamics simulations, we investigate the interaction between decagonal quasicrystals and obstacles (shrinkage pores) during solidification of an Al 79 ⁢Co 6 ⁢Ni 15 alloy. The results reveal defect-free quasicrystal growth around the pores, independent of pore size and geometry. Our findings point to a universal feature of aperiodic solids: their ability to accommodate structural disruptions due to additional (phasonic) degrees of freedom. This Letter opens new possibilities for synthesizing large-scale, single quasicrystals from a liquid metal precursor.

alloys↗
DOE OSTI2023

Polarity-dependence of the nonlinear dielectric response in interfacial water

Molecular dynamics simulations are used to study the nonlinear dielectric responses of a confined aqueous film in a planar nanopore under perpendicular electric fields at varied voltages between confining graphene sheets. Dielectric saturation reminiscent of the bulk phase behavior is prevalent at very strong fields, whereas we observe a nonmonotonic permittivity dependence on the electric field at intermediate strengths where field-alignment and spontaneous polarization of interfacial water are of comparable magnitude. The coupling between the two effects results in distinct dielectric responses at opposite confinement walls. The normal component of both the differential dielectric constant and dielectric difference constant tensors averaged over the region closer to the wall under an incoming electric field (field pointing from the liquid to the solid phase) initially increases with the strength of the imposed field. The differential permittivity peaks at a field strength previously shown to offset the surface-induced orientation bias of hydration molecules at this wall. Further strengthening of the field results in a conventional saturation behavior. At the opposite wall (subject to outgoing field) and in the central region of the water slab, the nonlinear dielectric response resembles bulklike saturation. The conditions at the permittivity extremum coincide with the window of accelerated reorientation rates of interfacial water molecules under an incoming field we uncovered in earlier molecular dynamics analyses.

Chemistry↗
DOE OSTI2019

Elucidating Interfacial Phenomena between Solid-State Electrolytes and the Sulfur-Cathode of Lithium–Sulfur Batteries

Sulfide solid-state electrolytes are promising materials to be used in Li–S batteries because they have the potential to mitigate some of the challenges currently faced by this technology such as the dissolution of polysulfide species and the flammability of conventional organic-based electrolytes. Nevertheless, the electrode/electrolyte interfaces are not well-understood yet. In this study, we use density functional theory and ab initio molecular dynamics simulations to investigate the formation of the cathode–electrolyte interface at two limiting states of charge: S 8 and Li 2 S. β-Li 3 PS 4 and two lithium halide-doped sulfides (Li 6 PS 5 Cl and Li 7 P 2 S8I) were selected because of their Li-ion superionic conductivity and compatibility with sulfur electrodes. The work of adhesion and interfacial energy of the interfaces were calculated as descriptors of the interfacial properties, and then the mechanisms of the interfacial reactions taking place were characterized. Lastly, a charge transfer analysis is presented, suggesting that, overall, the S 8 electrode will oxidize the electrolyte, whereas the fully discharged Li 2 S cathode will reduce it at different extents depending on the reactivity of the exposed facet.

25 ENERGY STORAGE↗
DOE OSTI2020

Mode-specific pressure effects on the relaxation of an excited nitromethane molecule in an argon bath

The vibrational and rotational mode-specific relaxations of CH 3 NO 2 with 50 kcal/mol of initial internal energy in an argon bath is computed at 300 K at pressures of 10-400 atm. This work uses archived information from our previously published [J. Chem. Phys. 142, 014303 (2015)] molecular dynamics simulations and employs our previous published [J. Chem. Phys. 151, 034303 (2019)] method for projecting time-dependent Cartesian velocities onto normal mode eigenvectors. The computed relaxations cover three types of energies: vibrational, rotational, and Coriolis. In general, rotational and Coriolis relaxations in all modes are initially fast followed by an orders of magnitude slower relaxation. For all modes, that slower relaxation rate is approximately comparable to the vibrational relaxation rate. For all three types of energies, there are small-scale mode-to-mode variations. Of particular prominence is the exceptionally fast relaxation shared in common by the external rotation about the C-N axis, the internal hindered rotation of the CH 3 group relative to the NO 2 group, and the symmetric stretch of the CH 3 group.

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

Elucidating the water–anatase TiO 2 (101) interface structure using infrared signatures and molecular dynamics

The structure and dynamics of water on solid surfaces critically affect the chemistry of materials in ambient and aqueous environments. Here, we investigate the hydrogen bonding network of water adsorbed on the majority (101) surface of anatase TiO 2 , a widely used photocatalyst, using polarization- and azimuth-resolved infrared spectroscopy combined with neural network potential molecular dynamics simulations. Our results show that one monolayer of water saturates the undercoordinated titanium (Ti 5c ) sites, forming one-dimensional chains of molecule hydrogen bonded to surface undercoordinated bridging oxygen (O 2c ) atoms. As the coverage increases, water adsorption on O 2c sites leads to significant restructuring of the water monolayer and the formation of a two-dimensional hydrogen bond network characterized by tightly bound pairs of water molecules on adjacent Ti 5c and O 2c sites. This structural motif likely persists at ambient conditions, influencing the reactions occurring there. In conclusion, the results reported here provide critical details of the structure of the water–anatase (101) interface that were previously hypothesized but unconfirmed experimentally.

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

Two-dimensional mapping of absolute OH densities in an atmospheric pressure plasma effluent via planar laser-induced fluorescence: effects of He/H 2 O and He/O 2 mixtures in N 2 and air, with and without solid targets

Planar laser-induced fluorescence (LIF) was employed to measure the absolute density of hydroxyl radicals (OH) in the effluent of the COST Reference Microplasma Jet for two feed gas mixtures: He/H 2 O and He/O 2 . Experiments were conducted with the effluent propagating into air and N 2 environments. For the He/H 2 O case, measurements were also performed with the effluent impinging on a solid target at varying distances from the jet nozzle. Calibration of the OH-LIF signal from the COST-Jet was achieved by comparing it to a reference signal generated by the photofragmentation of H 2 O 2 . Results demonstrated that OH densities were sustained longer when the effluent propagates in a nitrogen environment compared to air, particularly with water added to the feed gas. The broader OH distribution in N 2 suggests slower consumption due to the absence of oxygen, which accelerates OH depletion in air via reactions involving O 2 and HO 2 . Even when water was not added to the feed, as in the He/O 2 case, appreciable OH densities were observed, due to gas impurities and reactive species interactions with atmospheric humidity, forming reaction fronts that delineate the gas flow. Two-dimensional fluid dynamics simulations elucidated the influence of atmospheric gas entrainment and solid targets on the OH distribution. Experimental trends were further compared with a zero-dimensional chemistry model to explore OH production and consumption mechanisms in air and nitrogen environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effects of Ligand Chemistry on Ion Transport in 2D Hybrid Organic–Inorganic Perovskites

2D hybrid organic–inorganic perovskites are potentially promising materials as passivation layers that can enhance the efficiency and stability of perovskite photovoltaics. The ability to suppress ion transport is proposed as a stabilization mechanism, yet an effective characterization of relevant modes of halide diffusion in 2D perovskites is nascent. In light of this knowledge gap, molecular dynamics simulations with enhanced sampling and experimental validation to systematically characterize how ligand chemistry in seven (R-NH 3 ) 2 PbI 4 systems impacts halide diffusion, particularly in the out-of-plane direction is combined. It is found that increasing stiffness and length of ligands generally inhibits ion transport, while increasing ligand polarization generally enhances it. Structural and energetic analyses of the migration pathways provide quantitative explanations for these trends, which reflect aspects of the disorder of the organic layer. Overall, this mechanistic analysis greatly enhances the current understanding of halide migration in 2D hybrid organic–inorganic perovskites and yields insights that can inform the design of future passivation materials.

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

Tuning the Functionalities of Porous Liquids for Emergent Gas-Capture Properties

Type 3 Porous Liquids (PLs) are a class of materials with the potential to revolutionize gas capture, storage, and utilization. These PLs are formed by suspending sorbent nanoparticles (e.g., metal−organic frameworks) in sterically excluded solvents, creating permanent porosity for gas capture in a processable, low-viscosity phase. Herein, a computational study revealed sorbent surface functionalization strategies to enhance CO 2 sorption, and the molecular structural signatures underpinning the enhancements in gas uptake. PLs composed of a ZIF-8 surface functionalized with 3-amino-1,2,4-triazole (Atz) in glyceryl triacetate were targeted for emergent CO 2 capture, exceeding that of the unfunctionalized ZIF-8 PL. ZIF-8 was surface functionalized with Atz at various surface coverage fractions (f), and classical molecular dynamics simulations predicted an increase in CO 2 sorption capacity with increasing f, up to f = 0.75. Additionally, detailed structural analyses revealed that solvent orientational order, derived from the solvent triplet-angle distribution, can identify the gas-capture potential of a PL without requiring computationally expensive direct modeling of the CO 2 sorption. Combined with experimental validation, initial computational screening of PL compositions promises to accelerate the discovery of PL compositions for novel gas separation materials platforms.

gas capture↗
DOE OSTI2025

Carbon-carbon interactions in warm dense titanium carbide

Encouraged by experimental reports of diamond precipitation in carbon-containing materials under warm dense matter conditions and in shock compression experiments, we examine the behavior of carbon in TiC using density functional theory-molecular dynamics simulations. Two polymorphs of TiC are considered, the ambient-pressure B1 ($Fm$$\overline{3}$$m$) structure in which C-C interactions are prevented by geometric constraints, and a high-pressure Cmcm structure in which C atoms condense into zigzag chains. Chemistry-inspired bonding analyses confirm the covalently bound nature of these chains and further illuminate important interatomic interactions in both structures. Upon melting of B1 TiC, new short-range C-C interactions develop in the liquid, while the pre-existing C-C interactions in Cmcm TiC persist in the liquid. Here, the resulting carbon networks in the melt provide a promising environment for eventual diamond nucleation.

Ab initio molecular dynamics↗
DOE OSTI2020

North Dakota Integrated Carbon Storage Complex Feasibility Study. Final report

In spring 2017, the Energy & Environmental Research Center (EERC) initiated an effort to determine the feasibility of developing a commercial-scale CO 2 geologic storage complex able to store 50+ million tonnes (Mt) of CO 2 in central North Dakota safely, permanently, and economically. The objective was to fulfill the goals of the U.S. Department of Energy (DOE) Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative and address technical and nontechnical challenges specific to commercial-scale deployment of a CO2 storage project. The findings clearly show that the concept of capturing CO 2 from a lignite-fired electrical generation facility in central North Dakota and safely and permanently storing the CO 2 in the deep subsurface is indeed technically, economically, and socially feasible. This project evaluated two study areas and their respective geologic storage complexes located adjacent to separate coal-fired facilities in North Dakota: the Basin Electric Power Cooperative (BEPC)-owned Great Plains Synfuels Plant (GPSP) and the Minnkota Power Cooperative (Minnkota)-owned Milton R. Young Station (MRYS). These locations, one with CO 2 capture in place and an existing CO 2 pipeline, are bolstered by progressive North Dakota pore space ownership and long-term liability laws. These elements and a motivated team created an ideal synergistic scenario for ensuring success of the CarbonSAFE Initiative and promoting North Dakota’s statewide vision for carbon management. The project included drilling two new geologic characterization wells, integrating an existing 3-D seismic survey, creating a geologic model subsequently used for injection simulation, a risk assessment, public outreach, and generating a site development plan based on results. In addition, the performance of select National Risk Assessment Partnership tools was evaluated. The geologic characterization wells were drilled ~5600 feet deep to the Broom Creek Formation; ~350 feet of core was retrieved from each well. The core included the Broom Creek (targeted injection zone) and a portion of the overlying Opeche Shale (seal). The Flemmer-1 well, west of Beulah, North Dakota, yielded 169 feet of sandstone. The BNI-1 well located south of Center, North Dakota, yielded 124 feet of sandstone. In each case, laboratory analysis of the sandstone showed permeability in the 300–1000-mD range, with porosity of 20%–30%. The Flemmer-1 well was sited within the boundaries of an existing 3-D seismic survey. Colocating the well with the seismic survey maximized the relationship between new and legacy data and developed a first-of-its-kind interpretation of the geologic fabric of the Broom Creek. Geologic characterization data were integrated into a 5544-mi 2 geocellular model that encompassed both new wells and stratigraphy from the surface to the Amsden Formation (underlying the Broom Creek). The model was later expanded vertically to include the deeper Black Island–Deadwood interval to examine its potential viability as a storage target. The geocellular model provided the foundation for dynamic simulation of CO 2 into the Broom Creek. Results of the simulation suggest that the Broom Creek could accept the DOE target rate of 2 Mt/yr of CO 2 into as few as two wells. To bracket the expected capture from MRYS, simulations were also investigated for a 4-Mt/yr rate near MRYS. Although more wells are needed (two additional), the Broom Creek still has the storage resource to accept the CO 2 at the increased rate. A risk assessment exercise was conducted to identify and assess technical and nontechnical risks that could prevent potential candidate storage complexes within the study area from serving as commercial storage sites. The assessment identified and evaluated six technical risk categories: 1) CO 2 injectivity, 2) storage capacity, 3) lateral migration of CO 2 , 4) lateral pressure propagation, 5) vertical migration of CO 2 or formation brine, and 6) induced seismicity. Following two rounds of analysis and scoring, no risks were determined to preclude continued efforts to develop carbon capture, utilization, and storage (CCUS) in central North Dakota. The risk assessment results will be used to guide future site characterization, modeling and simulation, and monitoring activities. A specific nontechnical strategic risk based on challenges that may be realized in amalgamation of pore space resulted in vertically expanding the geologic model to incorporate the potential for stacked storage in multiple saline reservoirs. By including the Black Island–Deadwood interval (the basal sedimentary reservoir in this region), the amalgamated areal extent could be reduced by as much as 45%. Working with a smaller geographic area reduces risks and costs associated with monitoring and pore space leasing. An economic evaluation incorporating capture; transport (<5 mi); Class VI wells; permitting; and monitoring, verification, and accounting suggests implementing commercial-scale CCUS is economically attractive if the federal tax benefits of 45Q are included. This is validated by Minnkota’s continued pursuit of CO 2 capture and geologic storage at MRYS through its Project Tundra initiative, indicating that there is a business case for CCUS in central North Dakota. Currently, North Dakota is the only state with underground injection control (UIC) Class VI primacy. Built into the North Dakota Century Code is a series of regulatory requirements that guide the process to obtain a Class VI CO 2 storage facility permit. As part of this project, a site development plan was compiled to assure compliance with North Dakota’s requirements to permit a commercial-scale CO 2 storage operation and includes a prospective time line encompassing a general breakdown of activities. In total, an estimated 30 months would be needed to execute the necessary steps to attain a North Dakota CO 2 storage facility permit. Outreach was an integral part of the project and encompassed any project-related activity that had contact or exposure beyond the project team. The goals of outreach were to foster an environment from which stakeholders could make informed decisions about the project and gauge community receptiveness to a CCUS project. A consistent set of messages and outreach products were developed in conjunction with an outreach advisory board that integrated project partners and team members. To gauge public acceptability of geologically storing CO 2 , 5611 households in the project area were invited to participate in an online survey. The survey results indicate that the public attitude regarding CCUS is neutral to positive, with strong sentiment that CO 2 capture and storage may be an approach to maintain the economic vitality of the region. To achieve project objectives, critical support in the form of financial backing, engineering evaluations, site access, outreach collaboration, operations data, risk assessment/evaluation, and software access was provided by BEPC, the North Dakota Industrial Commission Lignite Research Council, ALLETE Clean Energy, BNI Energy, North American Coal Corporation, Minnkota, Prairie Public Broadcasting, Computer Modelling Group Ltd., and Schlumberger.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Dynamic Copper Site Redispersion through Atom Trapping in Zeolite Defects

Single-site copper-based catalysts have shown remarkable activity and selectivity for a variety of reactions. However, deactivation by sintering in high-temperature reducing environments remains a challenge and often limits their use due to irreversible structural changes to the catalyst. Here, we report zeolite-based copper catalysts in which copper oxide agglomerates formed after reaction can be repeatedly redispersed back to single sites using an oxidative treatment in air at 550 °C. Under different environments, single-site copper in Cu–Zn–Y/deAlBeta undergoes dynamic changes in structure and oxidation state that can be tuned to promote the formation of key active sites while minimizing deactivation through Cu sintering. For example, single-site Cu 2+ reduces to Cu 1+ after catalyst pretreatment (270 °C, 101 kPa H 2 ) and further to Cu0 nanoparticles under reaction conditions (270–350 °C, 7 kPa EtOH, 94 kPa H 2 ) or accelerated aging (400–450 °C, 101 kPa H 2 ). After regeneration at 550 °C in air, agglomerated CuO was dispersed back to single sites in the presence and absence of Zn and Y, which was verified by imaging, in situ spectroscopy, and catalytic rate measurements. Ab initio molecular dynamics simulations show that solvation of CuO monomers by water facilitates their transport through the zeolite pore, and condensation of the CuO monomer with a fully protonated silanol nest entraps copper and reforms the single-site structure. Importantly, the capability of silanol nests to trap and stabilize copper single sites under oxidizing conditions could extend the use of single-site copper catalysts to a wider variety of reactions and allows for a simple regeneration strategy for copper single-site catalysts.

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

Slip-free multiplication and complexity of dislocation networks in FCC metals

Abstract During plastic deformation of crystalline solids, intricate networks of dislocation lines form and evolve. To capture dislocation density evolution, prominent theories of crystal plasticity assume that 1) multiplication is driven by slip in active slip systems and 2) pair-wise slip system interactions dominate network evolution. In this work, we analyze a massive database of over 100 discrete dislocation dynamics simulations (with cross-slip suppressed), and our findings bring both of these assumptions into question. We demonstrate that dislocation multiplication is commonly observed on slip systems with no applied stress and no plastic strain rate, a phenomenon we refer to as slip-free multiplication. We show that while the formation of glissile junctions provides one mechanism for slip-free multiplication, additional mechanisms which account for the influence of coplanar interactions are needed to fully explain the observations. Unlike glissile junction formation which results from a binary reaction between a pair of slip systems, these new multiplication mechanisms require higher order reactions that lead to complex network configurations. While these complex configurations have not been given much attention previously, they account for about 50% of the line intersections in our database.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Impact of large-scale effects on mass transfer and concentration polarization in Reverse Osmosis membrane systems

Here, we present well-resolved computational fluid dynamics simulations of a large-scale reverse osmosis membrane-spacer configuration (1 m). Our computational model solves the flow and transport equations with variable solute-dependent properties. We utilize a high resolution computational mesh to resolve all relevant length scales associated with spacer-induced mixing and thin concentration boundary layers. An important contribution of this work is the development of a modified mass-transfer correlation that accounts for the development of the concentration boundary layer along the channel. A set of 2D axisymmetric simulations were performed for a spiral wound module layer with varying cross-flow conditions and spacer diameters which indicate a significant entrance length effect for concentration profile development at lower flow rates while mixing effects dominate at higher flow rates. The mass-transfer correlations at higher flow rates compare well with published correlations while a surrogate model for Sherwood number was obtained that depends on an additional similarity variable that accounted for entrance length effects at lower flow rates. Finally, a large-scale membrane-spacer design relevant to high-pressure reverse osmosis is studied with a non-uniform arrangement of spacers, which indicate a substantial saving in pressure drop (~40%) compared to traditional uniformly spaced pattern with minor variations (~2%) in concentration polarization, product water quality (~1%) and water recovery (~7%) compared to a uniform spacer pattern.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Structure, bonding and electronic characteristics of amorphous Se

Ovonic threshold switching (OTS) selector can effectively improve the storage density and suppress the leakage current of advanced phase-change memory. As a prototypical OTS material, amorphous GeSe is widely investigated. But the attentions paid to amorphous Se (i.e., the functional constituent in amorphous GeSe) is very limited up to now. Here, in this work, we have explored the structure, bonding and electronic characteristics of amorphous Se using ab initio molecular dynamics simulations. The results reveal that the Se atoms in amorphous Se tend to form the 2-coordinated configurations, and they connect with each other to form the long chains. The fraction of vibrational density of state located in high frequency range is relatively large, and the formation energy of Se-Se bond is as large as 4.44 eV, hinting that the Se-Se bonds in chains possess a high stability. In addition, the mid-gap state related to the OTS behavior is also found in the amorphous Se despite of the small proportion. Our findings enrich the knowledge of the amorphous Se, which contributes to the applications of Se-based OTS selectors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗