Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulated annealing”

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 91 records · Page 5

Computational Fermi Level Engineering and Doping-Type Conversion of Ga2O3 Via Three-Step Processing

Ga2O3 is being actively explored for power electronics, and other applications due to its ultra- wide bandgap and low projected fabrication cost of high-quality crystals. N-type doping of Ga2O3 can be achieved and tuned, but p-type doping faces fundamental obstacles. Successful engineering of Ga2O3 based devices requires critical control of doping density, Fermi level position, providing opportunities for predictive process simulation. We use first-principles defect equilibrium calculations to simulate a 3-step growth- annealing-quench protocol for hydrogen assisted Mg doping in Ga2O3, taking into account the hydrogen- oxygen-water gas phase equilibrium. We predict type conversion to a net p-type regime following O-rich annealing after growth under reducing conditions in the presence of H2. We show that there is an optimal temperature that maximizes the net acceptor density during the equilibrium annealing step for a given Mg doping level. Quenching of non-equilibrium annealed samples then results in a Fermi level EF below mid-gap down to about EV +1.5 eV, creating a significant number of uncompensated neutral MgGa0 acceptors. This type converted Ga2O3 can create significant built-in field in a p-n junction with an adjoining n-type material and could enable vertical MOSFETs.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗

Computational Fermi Level Engineering and Doping-Type Conversion of Ga2O3 Via Three-Step Processing

Ga2O3 is being actively explored for power electronics, deep-ultraviolet optoelectronics, and other applications due to its ultra-wide bandgap and low projected fabrication cost of large-size and high-quality crystals. N-type doping of Ga2O3 can be achieved and tuned, but p-type doping faces fundamental obstacles due to deep character of acceptor levels and polaron transport of resulting holes. However, successful engineering of Ga2O3 based devices requires critical control of doping density, Fermi level position, and free carrier concentration, providing opportunities for predictive process simulation. We use first-principles defect theory and defect equilibrium calculations to simulate a 3-step growth-annealing-quench protocol for hydrogen assisted Mg doping in Ga2O3, taking into account the hydrogen-oxygen-water gas phase equilibrium. We predict type conversion to a net p-type regime following O-rich annealing after growth under reducing conditions in the presence of H2. This process is similar to the Mg acceptor activation by H removal in GaN. We show that there is an optimal temperature that maximizes the net acceptor density during the equilibrium annealing step for a given Mg doping level. Quenching of non-equilibrium annealed samples then results in a Fermi level EF below mid-gap down to about EV +1.5 eV, creating a significant number of uncompensated neutral MgGa0 acceptors. The resulting free hole concentration in quenched samples is very low (109 – 1012 cm-3) due to deep energy level of these Mg acceptors, but this type converted Ga2O3 material can create a significant built-in field in a p-n junction with an adjoining n-type material. Additionally, the electron concentration is greatly suppressed in such acceptor-doped Ga2O3, which could enable the use as a current blocking layer to fabricate normally-off (enhancement-mode) vertical Ga2O3 based metal-oxide-semiconductor field effect transistors (MOSFETs), and as a guard ring for edge termination in MOSFETs and Schottky barrier diodes (SBDs) with increased breakdown voltage.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗

3-regular three-XORSAT planted solutions benchmark of classical and quantum heuristic optimizers

With current semiconductor technology reaching its physical limits, special-purpose hardware has emerged as an option to tackle specific computing-intensive challenges. Optimization in the form of solving quadratic unconstrained binary optimization problems, or equivalently Ising spin glasses, has been the focus of several new dedicated hardware platforms. These platforms come in many different flavors, from highly-efficient hardware implementations on digital-logic of established algorithms to proposals of analog hardware implementing new algorithms. In this work, we use a mapping of a specific class of linear equations whose solutions can be found efficiently, to a hard constraint satisfaction problem (three-regular three-XORSAT, or an Ising spin glass) with a 'golf-course' shaped energy landscape, to benchmark several of these different approaches. We perform a scaling and prefactor analysis of the performance of Fujitsu's digital annealer unit (DAU), the D-Wave advantage quantum annealer, a virtual MemComputing machine, Toshiba's simulated bifurcation machine (SBM), the SATonGPU algorithm from Bernaschi et al, and our implementation of parallel tempering. We identify the SATonGPU and DAU as currently having the smallest scaling exponent for this benchmark, with SATonGPU having a small scaling advantage and in addition having by far the smallest prefactor thanks to its use of massive parallelism. Furthermore, our work provides an objective assessment and a snapshot of the promise and limitations of dedicated optimization hardware relative to a particular class of optimization problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evidence for ionization damage in mid-wave infrared nBn detectors

The effects of 63 MeV proton and 60Co gamma irradiation on the operation of 4.3 μm cutoff nBn photodetectors are demonstrated separately, and both are shown to yield a total ionizing dose (TID) effect. The effect is shown here in an InAsSbBi nBn detector and has been observed in other bulk alloy nBn detectors, and is unusual as it is notably absent in superlattice nBn detectors. The non-antireflection coated detectors exhibit a pre-radiation quantum efficiency of 17% at 3.3 μm wavelength and a dark current density of 50 μA/cm2, or roughly 300× the Rule 07 expectation, at their ideal operating voltage of −0.4 V bias at 150 K. Step-wise proton irradiation and in situ measurement indicate that the dark current increases to about 400× Rule 07 at the highest proton dose level of 150 krad(Si) (9.10 × 1011 p+/cm2), while the quantum efficiency is degraded at a relatively faster rate than the majority of analogous detectors characterized by our lab. Both the photocurrent and dark current are also shown to exhibit a turn-on voltage magnitude reduction of 100 mV following either gamma or proton irradiation, a trend which is attributable to negative trapped charge at the barrier interface (TID effect). This theory is further supported by an observed capacitance density magnitude reduction with dose and affirmed with Silvaco TCAD simulations. Following both proton exposure and subsequent anneal and gamma exposure and subsequent anneal; dark current, photocurrent, and CV all approach their pre-radiation baseline values.

Physics↗

Role of mechanical stress localizations on the radiation hardness of AlGaN/GaN high electron mobility transistors

Multi-material, multi-layered systems such as AlGaN/GaN high electron mobility transistors (HEMTs) contain residual mechanical stresses that arise from sharp contrasts in device geometry and materials parameters. These stresses, which can be either tensile or compressive, are difficult to detect and eliminate because of their highly localized nature. We propose that their high-stored internal energy makes potential sites for defect nucleation sites under radiation, particularly if their locations coincide with the electrically sensitive regions of a transistor. In this study, we validate this hypothesis with molecular dynamic simulation and experiments exposing both pristine and annealed HEMTS to 2.8 MeV Au +3 irradiation. Our unique annealing process uses mechanical momentum of electrons, also known as the electron wind force (EWF) to mitigate the residual stress at room temperature. High-resolution transmission electron microscopy and cathodoluminescence spectra reveal the reduction of point defects and dislocations near the two-dimensional electron gas region of EWF-treated devices compared to pristine devices. The EWF-treated HEMTs showed relatively higher resilience with approximately 10% less degradation of drain saturation current and ON-resistance and 5% less degradation of peak transconductance. Both mobility and carrier concentration of the EWF-treated devices were less impacted compared to the pristine devices. Our results suggest that the lower density of nanoscale stress localization contributed to the improved radiation tolerance of the EWF-treated devices. Intriguingly, the EWF is found to modulate the defect distribution by moving the defects to electrically less sensitive regions in the form of dislocation networks, which act as sinks for the radiation induced defects and this assisted faster dynamic annealing.

AlGaN/GaN HEMTs↗

Intergranular Oxidation of Unstressed Alloy 600 SA and Alloy 600TT in PWR Primary Water

Unstressed coupons of Alloy 600 SA (solution-annealed) and Alloy 600TT were exposed to simulated pressurized water reactor (PWR) primary water at 360oC. The depth of intergranular oxide penetration, as well as the microstructure and microchemistry associated with intergranular (IG) oxidation of selected samples, were characterized by SEM, serial FIB sectioning, TEM and APT. There is no clear trend to indicate that IG oxidation depths in unstressed coupons are different between Alloy 600 SA and TT. At the same time, all SA+TT materials examined evidenced significant interaction between oxidation and Cr-rich GB precipitates. In addition, the spatial extent of GB compositional changes ahead of the oxidation front (e.g. Cr depletion and Ni enrichment) appears different for Alloy 600 SA vs. SA+TT. In SA heats exposed for only 1000 h (42 days), the composition of some GBs changes over distances of 500–1000 nm. SA+TT materials appear not to develop such long-range GB compositional changes (10s nm alteration distances only).

Alloy 600, IGA, PWR, Oxidation, FIB, TEM, APT↗

Two-step annealing induced structural rejuvenation: A cause for memory effect in metallic glasses

The structural origin of the memory effect is still elusive in glassy materials. Here, in this letter, the memory effect in a series of metallic glasses (MGs) is observed by two step annealing using differential scanning calorimetry (DSC) and molecular dynamics simulations. It demonstrates that the Kovacs peak can be reflected from not only the system potential energy but also the atomic dynamics, showing the memory effect could originate from the formation of more loose packing configurations around mobile atoms by structural rejuvenation in the second step of high temperature annealing. In contrast, the high-low temperature annealing protocol can only induce the monotonic relaxation. Therefore, our results pave a pathway to better understand the origin of memory effect in MGs and other glassy materials from their dynamics and local atomic structure.

Memory effect↗

Evolution of dislocation loops and voids in post-irradiation annealed ThO 2 : A combined in-situ TEM and cluster dynamics investigation

The effect of isochronal annealing on the evolution of dislocation loop and void population in proton irradiated ThO 2 has been investigated. Post-irradiation annealing in other actinide oxides like UO 2 shows significant loop coarsening. ThO 2 samples were irradiated with 2 MeV protons up to a dose of 0.1 dpa at 600°C. Post-irradiation isochronal annealing was performed at 600, 800, 1000 and 1100°C for 1 h at each temperature using in-situ TEM. Only faulted 1/3<111> type dislocation loops were observed, and their sizes and distribution were characterized. Further, the population of self-interstitial atom (SIA) dislocation loops did not show any significant growth and coarsening. Additionally, nanometric voids were observed at annealing temperatures of 1000 and 1100°C. Using cluster dynamics (CD), we have studied the nucleation and growth of point defects and defect clusters, i.e., SIA prismatic dislocation loops and nanometric and sub-nanometric voids in proton irradiated ThO 2 . The CD model was further utilized to predict the growth and coarsening of loops and voids during isochronal annealing at the experimental and higher temperatures. The model did not predict significant SIA loop growth which closely corresponds to the TEM observations. CD predicted SIA loop coarsening is insignificant even at high annealing temperature of 1500°C because the model only considers the growth of defect clusters by absorption of like point defects, i.e., SIA loops absorb interstitials and voids absorb vacancies, and cannot account for their migration and coalescence due to elastic interaction. The CD model also predicts the evolution of nanometric voids having mean size within the error bounds of TEM observations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Revealing a Pathway for Low-Temperature Recrystallization in Germanium

Thermally activated annealing in semiconductors faces inherent limitations, such as dopant diffusion. Here, a nonthermal pathway is demonstrated for a complete structural restoration in predamaged germanium via ionization-induced recovery. By combining experiments and modeling, this study reveals that the energy transfer of only 2.4 keV nm −1 from incident ions to target electrons can effectively annihilate pre-existing defects and restore the original crystalline structure at room temperature. Moreover, it is revealed that the irradiation-induced crystalline-to-amorphous (c/a) transformation in Ge is reversible, a phenomenon previously considered unattainable without additional thermal energy imposed during irradiation. For partially damaged Ge, the overall damage fraction decreases exponentially with increasing fluence. Surprisingly, the recovery process in preamorphized Ge starts with defect recovery outside the amorphous layer and a shrinkage of the amorphous thickness. After this initial stage, the remaining damage decreases slowly with increasing fluence, but full restoration of the pristine state is not achieved. These differences in recovery are interpreted in the framework of structural differences in the initial defective layers that affect recovery kinetics. This study provides new insights on reversing the c/a transformation in Ge using highly-ionizing irradiation and has broad implications across materials science, radiation damage mitigation, and fabrication of Ge-based devices.

athermal recovery↗

Additive‐Free, Low‐Temperature Crystallization of Stable α‐FAPbI 3 Perovskite

Abstract Formamidinium lead triiodide (FAPbI 3 ) is attractive for photovoltaic devices due to its optimal bandgap at around 1.45 eV and improved thermal stability compared with methylammonium‐based perovskites. Crystallization of phase‐pure α‐FAPbI 3 conventionally requires high‐temperature thermal annealing at 150 °C whilst the obtained α‐FAPbI 3 is metastable at room temperature. Here, aerosol‐assisted crystallization (AAC) is reported, which converts yellow δ‐FAPbI 3 into black α‐FAPbI 3 at only 100 °C using precursor solutions containing only lead iodide and formamidinium iodide with no chemical additives. The obtained α‐FAPbI 3 exhibits remarkably enhanced stability compared to the 150 °C annealed counterparts, in combination with improvements in film crystallinity and photoluminescence yield. Using X‐ray diffraction, X‐ray scattering, and density functional theory simulation, it is identified that relaxation of residual tensile strains, achieved through the lower annealing temperature and post‐crystallization crystal growth during AAC, is the key factor that facilitates the formation of phase‐stable α‐FAPbI 3 . This overcomes the strain‐induced lattice expansion that is known to cause the metastability of α‐FAPbI 3 . Accordingly, pure FAPbI 3 p–i–n solar cells are reported, facilitated by the low‐temperature (≤100 °C) AAC processing, which demonstrates increases of both power conversion efficiency and operational stability compared to devices fabricated using 150 °C annealed films.

36 MATERIALS SCIENCE↗

Athermal annealing of pre-existing defects in crystalline silicon

Systematic investigations of electronic energy loss (S e ) effects on pre-existing defects in crystalline silicon (Si) are crucial to provide reliance on the use of ionizing irradiation to anneal pre-existing defects, leading to successful implementation of this technology in the fabrication of Si-based devices. In this regard, the S e effects on nonequilibrium defect evolution in pre-damaged Si single crystals at 300 K has been investigated using intermediate-energy ions (12 MeV O and Si ions) that interact with the pre-damaged surface layers of Si mainly by ionization, except at the end of their range where the nuclear energy loss (S n ) is no longer negligible. Furthermore, under these irradiation conditions, experimental results and molecular dynamics simulations have revealed that pre-existing disorder in Si can be almost fully annealed by subsequent irradiation with intermediate-energy incident ions with S e values as low as 1.5–3.0 keV/nm. Selective annealing of pre-existing defect levels in Si at room temperature can be considered as an effective strategy to mediate the transient enhanced diffusion of dopants in Si. This approach is more desirable than the regular thermal annealing, which is not compatible with the processing requirements that fall below the typical thermal budget.

36 MATERIALS SCIENCE↗

Flow annealed importance sampling bootstrap meets differentiable particle physics

High-energy physics requires the generation of large numbers of simulated data samples from complex but analytically tractable distributions called matrix elements. Surrogate models, such as normalizing flows, are gaining popularity for this task due to their computational efficiency. We adopt an approach based on flow annealed importance sampling bootstrap (FAB) that evaluates the differentiable target density during training and helps avoid the costly generation of training data in advance. We show that FAB reaches higher sampling efficiency with fewer target evaluations in high dimensions in comparison to other methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Mechanical Response of High Burnup 17x17 PWR Fuel Rods Under Bending

As a part of the US Department of Energy Office of Nuclear Energy High Burnup Spent Fuel Data Project, Oak Ridge National Laboratory has performed 4pt-bend tests on a series of 17x17 PWR fuel rods that were discharged with rod average burnups between 50 and 60 GWd/MTU. The rods include four different kinds of fuel cladding: standard Zircaloy-4 (Zirc-4), low-tin (LT) Zirc-4, ZIRLO, and M5. Six rods of each type (combining Zirc-4 and LT Zirc-4 as one type) were tested to failure at room temperature and 200 oC. Three of each type experienced a full-length heat treatment cycle from RT to 400 oC to RT with an 8 hr soak at 400 oC to simulate the expected limiting vacuum drying condition, which potentially induces hydride reorientation and causes some annealing of irradiation damage in the cladding. Overall, the fuel rods show significant capacity to deflect in bending and exhibit non-linear deformation before failing. The calculated failure strains range from 1.5% to 3.8%. The simulated vacuum drying heat treatment had the largest effect on fuel rods with ZIRLO and M5 cladding. In both cases the heat treatment caused increases in the plastic strain at failure, suggesting some annealing of irradiation damage is occurring. The finite element modeling of the 4pt bending test assume the pellet is a collection of pellet wedge pieces. The results for the Zircaloy-4 rods suggest that the pellets are bonded to the cladding and friction controls the relative motion of the pellet wedge pieces.

Cantonwine, Paul↗

Aggregation of Size-Selected Oxide Clusters Deposited onto Au(111)

Kinetic Monte Carlo (kMC) simulations along with density functional theory (DFT) calculations were used to investigate the aggregation of size-selected Nb 3 O y (y = 5, 6, 7) clusters deposited onto the Au(111) surface. Recent STM experiments showed that the cluster binding sites and sizes of the cluster assemblies on the Nb 3 O y /Au(111) surfaces strongly depend on the stoichiometry of the clusters, i.e., the oxygen-to-niobium ratio. To better understand the origins of these differences, kMC simulations of the nucleation and growth of cluster assemblies were performed using energy barriers for diffusion and intercluster interactions estimated from DFT calculations of cluster binding and dimerization energies, respectively. Comparisons of the kMC simulations with STM images of the as-deposited Nb 3 O y /Au(111) surfaces at RT and after high temperature annealing were used to further optimize the energetics and gauge the importance of nearest neighbor interactions. The kMC simulations demonstrate that the assembly of Nb 3 O y clusters on Au(111) are largely controlled by the magnitude of the barriers for diffusion and interparticle-bond formation, while changes at higher temperatures are sensitive to the binding energies between nearest neighbors. Simulations for the Nb 3 O 5 and Nb 3 O 6 clusters, which exhibit smaller cluster assembly sizes in STM, required larger diffusion barriers as well as different barriers for interparticle binding, which reflected differences in DFT calculated dimerization energies. Furthermore, the results demonstrate the effectiveness of combined DFT and kMC calculations for understanding how the stoichiometry affects the aggregation of small oxide clusters on a metal surface.

36 MATERIALS SCIENCE↗

Computing molecular excited states on a D-Wave quantum annealer

Abstract The possibility of using quantum computers for electronic structure calculations has opened up a promising avenue for computational chemistry. Towards this direction, numerous algorithmic advances have been made in the last five years. The potential of quantum annealers, which are the prototypes of adiabatic quantum computers, is yet to be fully explored. In this work, we demonstrate the use of a D-Wave quantum annealer for the calculation of excited electronic states of molecular systems. These simulations play an important role in a number of areas, such as photovoltaics, semiconductor technology and nanoscience. The excited states are treated using two methods, time-dependent Hartree–Fock (TDHF) and time-dependent density-functional theory (TDDFT), both within a commonly used Tamm–Dancoff approximation (TDA). The resulting TDA eigenvalue equations are solved on a D-Wave quantum annealer using the Quantum Annealer Eigensolver (QAE), developed previously. The method is shown to reproduce a typical basis set convergence on the example $$\hbox {H}_2$$ H 2 molecule and is also applied to several other molecular species. Characteristic properties such as transition dipole moments and oscillator strengths are computed as well. Three potential energy profiles for excited states are computed for $$\hbox {NH}_3$$ NH 3 as a function of the molecular geometry. Similar to previous studies, the accuracy of the method is dependent on the accuracy of the intermediate meta-heuristic software called qbsolv.

74 ATOMIC AND MOLECULAR PHYSICS↗

Prediction of softening kinetics and recrystallization texture in non-isothermally annealed bulged tubes using CPFEM and CA models

A hierarchically coupled cellular automata (CA) model, crystal plasticity finite element method (CPFEM), and thermal finite element (FE) model is developed to predict the softening kinetics of the bulged steel tube during non-isothermal annealing. Through the developed model, the kinetics of softening mechanisms including static recovery (SRV) and static recrystallization (SRX), as well as the recrystallization texture are predicted. Later, the Johnson-Mehl-Avrami-Kohnogorov (JMAK) model based on the predicted SRX data is developed to interpret the recrystallization behavior of the material. To perform this study, diverse experimental tests including tube hydroforming (THF), annealing, uniaxial tensile test, hardness test, as well as microstructure observations through optical microscopy and Electron Backscatter Diffraction (EBSD) tests on steel tube are performed. The obtained experimental data are utilized to calibrate and verify the implemented CPFEM model for simulation of THF process, thermal FE model for prediction of the local temperature over annealing time, and CA algorithm for modeling of the softening kinetics and texture evolution throughout the annealing process. The study shows that the predicted deformation characteristics, softening kinetics, recrystallization texture and temperature profile during non-isothermal annealing are in good agreement with experimental data. During the annealing process, a total of four stages for the kinetics of softening mechanisms is observed: No softening; SRV only; SRV dominant; and SRX dominant. During the progress of SRX, the behavior of recrystallization is abruptly changed, confirming that two different mechanisms are controlling the kinetics of transformation.

36 MATERIALS SCIENCE↗

Electronic structure theory with molecular point group symmetries on quantum annealers

Quantum computation has the potential to revolutionize quantum chemistry through major speedups in computation times and an exponential reduction in computational resources. Here, we combine the symmetry-adapted Jordan–Wigner encoding based on the full Boolean symmetry group $\mathbb{Z}$$^{k}_{2}$ with our new implementation of the Xia–Bian–Kais (XBK) method for improving the efficiency of electronic structure theory calculations on quantum annealers, particularly by reducing the number of qubits needed to achieve the same accuracy. By providing a more extensive symmetry-adapted encoding (SAE) than previous work, we are able to simulate molecules larger than those previously reported that have been studied using methods developed for quantum annealers and without using an active space. We calculated the potential energy surfaces of H 2 , LiH, He 2 , H 2 O, O 2 , N 2 , Li 2 , F 2 , CO, BH 3 , NH 3 , and CH 4 , with the largest molecule in the STO-6G basis set requiring 16 qubits with our SAE, and compared them with full configuration interaction results. The application of SAE to the XBK method provides an exponential reduction in the size of the Hilbert space and scales well with the size of the problem. It does not introduce significant additional errors for even or large values of a key variational parameter that determines the number of ancilla qubits used in the XBK method’s Hamiltonian embedding, or for certain molecules such as He 2 and H 2 O. Here, we provide an explanation for this behavior and a recommendation on the usage of our method. In addition, we briefly discuss the potential of extracting electronic excited states from our method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Active doping controls the mode of failure in dense colloidal gels

Mechanical properties of disordered materials are governed by their underlying free energy landscape. In contrast to external fields, embedding a small fraction of active particles within a disordered material generates nonequilibrium internal fields, which can help to circumvent kinetic barriers and modulate the free energy landscape. In this work, we investigate through computer simulations how the activity of active particles alters the mechanical response of deeply annealed polydisperse colloidal gels. We show that the “swim force” generated by the embedded active particles is responsible for determining the mode of mechanical failure, i.e., brittle vs. ductile. We find, and theoretically justify, that at a critical swim force the mechanical properties of the gel decrease abruptly, signaling a change in the mode of mechanical failure. The weakening of the elastic modulus above the critical swim force results from the change in gel porosity and distribution of attractive forces among gel particles, while below the critical swim force, the ductility enhancement is caused by an increase of gel structural disorder. Above the critical swim force, the gel develops a pronounced heterogeneous structure characterized by multiple pore spaces, and the mechanical response is controlled by dynamical heterogeneities. We contrast these results with those of a simulated monodisperse gel that exhibits a nonmonotonic trend of ductility modulation with increasing swim force, revealing a complex interplay between the gel energy landscape and embedded activity.

Zhou, Tingtao (ORCID:000000021766719X)↗