Engineering Papers⌕ Search

Engineering topics

Qi, Liang

Publications and source records attributed to Qi, Liang.

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Mechanism and Kinetics of Ethanol–Acetaldehyde Conversion to 1,3-Butadiene over Isolated Lewis Acid La Sites in Silanol Nests in Dealuminated Beta Zeolite

Biomass-derived ethanol (EtOH) and acetaldehyde (AcH) conversion to 1,3-butadiene (1,3-BD) is an alternative process for 1,3-BD production. The present investigation reports the preparation and characterization of isolated La sites introduced into the silanol nests in DeAlBEA as well as detailed studies of the mechanism and kinetics for the conversion of an EtOH-AcH mixture to 1,3-BD. La sites supported on DeAlBEA are found to be present as (≡SiO) 2 La-OH groups that are H-bonded with adjacent Si-OH groups, possessing high C-C coupling activity and stability, superior to state-of-the-art Y-DeAlBEA. La sites supported on silica (La-SiO 2 ) with a similar chemical structure but no H-bonding interaction with Si-OH groups were prepared for comparison. Lewis acid La sites promote AcH aldol condensation, and the activity of such sites is nearly identical for both La-DeAlBEA and La-SiO 2 . Further, the rate of C 4 product formation increases by a factor of 4.8 upon addition of EtOH to the feed of AcH over La-DeAlBEA, whereas that over La/SiO 2 remains unchanged. Investigation of the mechanism and kinetics of AcH aldol condensation and EtOH-AcH conversion to 1,3-BD revealed two C-C bond forming pathways-AcH aldol condensation by Lewis acid La sites and direct coupling of EtOH-AcH over H-bonded (≡SiO) 2 La-OH···HO-Si≡ sites. This study provides important information about the role of the local environment of isolated Lewis acid sites and their effects on the direct coupling of EtOH and AcH to form 1,3-BD.

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

Native Defects Hybridization and Electronic Transport in Cu 2 Se–CuGaSe 2 Hierarchical Composites

Engineering the energy distribution of electronic defects in semiconductors can afford the coexistence of degenerate and nondegenerate transport behavior. By leveraging native defects in Cu 2 Se and CuGaSe 2 phases, the tunability of the concentration and energy distribution of active electronic defects in hierarchical (1-x)Cu 2 Se-(x)CuGaSe 2 composites is demonstrated. This study found that the electrical conductivity of various composites decreases with rising temperature indicating degenerate semiconducting behavior, whereas the carrier density increases with temperature, which is consistent with non-degenerate semiconductivity. Remarkably, this study found that while the carrier density and electrical conductivity drop with the increasing CuGaSe 2 content for samples with x ≤ 0.45 is consistent with an increase in the activation energy of “active” electronic defects, the sudden increase by 140% in the electrical conductivity and by 109% in the carrier density for near equimolar composition suggests the formation of a large density of electronic defects with lower activation energy. This unusual behavior is understood within the context of the hybridization of coexisting native electronic defects to form degenerate hybrid acceptor states with lower activation energy. The reported unique electronic transport can be leveraged for the development of more versatile electronic and optoelectronic devices with superior performance.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomistic simulations and machine learning of solute grain boundary segregation in Mg alloys at finite temperatures

Understanding solute segregation thermodynamics is the first step in investigating grain boundary (GB) properties, such as strong yttrium (Y) effects on grain growth and texture evolution in micro-scale polycrystalline magnesium (Mg) alloys. To estimate the average GB segregation behavior in low-solute-concentration Mg alloys (e.g., 2 at.% Y), a state-of-the-art spectral approach is applied based on a per-site segregation energy spectrum for Y solute atoms at zero K obtained from molecular statistics (MS) simulations of ~10 4 GB sites in Mg symmetric tilt GBs (STGBs). Although selected MS simulation results are consistent with verification by density functional theory (DFT) calculations, estimates of average segregation tendency based on the zero-K energy spectrum deviate from experimental observations. To resolve this problem, thermodynamic integration (TI) methods based on molecular dynamics (MD) simulations are used to determine the per-site segregation free energies of Y at representative GB sites, which show contributions beyond harmonic approximations can be important for certain GB sites at high temperatures. A surrogate model of per-site segregation free energy is constructed from a small set of TI data points using stacking cross-validation regressors and physics-informed descriptors. This model is applied to predict the Y segregation free energy spectra for thousands of GB sites in Mg STGBs with uncertainty quantification. Finally, the average segregation tendency predicted by the spectral approach based on the free energy spectra agrees well (within the uncertainty range) with experimental observations for micro-scale polycrystalline Mg alloys at typical thermomechanical processing temperatures (500 ~ 800 K), where thermodynamic equilibrium states are likely to be achieved due to fast diffusion.

Atomistic simulations↗
DOE OSTI2023

Multi-scale investigation of short-range order and dislocation glide in MoNbTi and TaNbTi multi-principal element alloys

Refractory multi-principal element alloys (RMPEAs) are promising materials for high-temperature structural applications. Here, we investigate the role of short-range ordering (SRO) on dislocation glide in the MoNbTi and TaNbTi RMPEAs using a multi-scale modeling approach. Monte carlo/molecular dynamics simulations with a moment tensor potential show that MoNbTi exhibits a much greater degree of SRO than TaNbTi and the local composition has a direct effect on the unstable stacking fault energies (USFEs). From mesoscale phase-field dislocation dynamics simulations, we find that increasing SRO leads to higher mean USFEs and stress required for dislocation glide. The gliding dislocations experience significant hardening due to pinning and depinning caused by random compositional fluctuations, with higher SRO decreasing the degree of USFE dispersion and hence, amount of hardening. Finally, we show how the morphology of an expanding dislocation loop is affected by the applied stress.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dehydrogenation of Propane and n -Butane Catalyzed by Isolated PtZn 4 Sites Supported on Self-Pillared Zeolite Pentasil Nanosheets

Propene and 1,3-butadiene are important building-block chemicals that can be produced by dehydrogenation of propane and butane over Pt catalysts. The challenge is to develop highly active and selective catalysts that are resistant to deactivation by Pt sintering and coke formation. We have recently shown that these objectives can be met for propane dehydrogenation using atomically dispersed Pt atoms anchored to neighboring ≡SiOZn-OH groups bonded to the framework dealuminated zeolite BEA. In the present study, we demonstrate that significantly superior performance can be achieved by using self-pillared pentasil (SPP) zeolite nanosheets. Following catalyst reduction in H 2 , atomic-resolution, scanning transmission electron microscopy (STEM) and X-ray absorption spectroscopy (XAS) indicate that Pt is stabilized in structures well approximated as (≡Si-O-Zn) 4-5 Pt. These species are highly active, selective, and stable for propane dehydrogenation (PDH) to give propene and for n-butane dehydrogenation (BDH) to give 1,3-butadiene. No catalyst deactivation was observed after 12 days of time on stream, and the selectivity remained at nearly 100% for PDH conducted at 823 K and a weight hourly space velocity (WHSV) of 1,350 h -1 . The apparent rate coefficient for PDH over this catalyst is significantly higher than that reported previously for Pt-containing catalysts. For BDH at 823 K and a WHSV of 3560 h -1 , the selectivity to butene isomers and 1,3-butadiene is 98.9%, and the selectivity to 1,3-butadiene is 45%. Here, we propose that the high catalyst stability observed during PDH and BDH is a consequence of a large fraction of the Pt-containing centers being located on the external surface of the zeolite nanosheets, where nascent coke precursors can desorb before condensing to form coke.

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

High temperature stability and transport characteristics of hydrogen in alumina via multiscale computation

Here, the impact of hydrogen charge states on the stability and transport characteristics of hydrogen interstitials in alumina polymorphs is evaluated by multiscale computational methods including density functional theory (DFT), ab initio molecular dynamics (AIMD) and machine learned force fields. Thermodynamic calculations show that the protonic H i +1 interstitial is the most stable defect species for most values of the electronic bandgap in both and amorphous alumina (Al 2 O 3 ). Further, active learned Gaussian approximation potentials (GAP) were developed using AIMD data to study temperature dependent long time proton diffusion in alumina. Diffusivity calculations from GAP-MD simulations are found to be comparable with of the AIMD data, while being ~340 times faster and scalable to larger systems. Comparisons with diffusivity values for other interstitial charge states (H i 0 and H i -1 ) and published experimental literature indicate that H i +1 diffusion is the likely mechanism of hydrogen transport. A good agreement is obtained between H i +1 diffusivity calculated in α-Al 2 O 3 from DFT: 5.05 10 -3 exp(-0.81 eV/k B /T) cm 2 /s and reported experiment: 9.7X10 -4 exp(-0.83 eV/k B /T) cm 2 /s. Computationally and experimentally calculated energy barriers (0.81 and 0.83 eV respectively) only differ by 2.5%. Similarly, the pre-exponential diffusion coefficients only differ by 0.5 orders of magnitude. Moreover, the diffusivity of H i +1 in amorphous Al 2 O 3 in the 1000–2000 K range is calculated to be 2.53X10 -2 exp(-0.89 eV/k B /T), just one order of magnitude higher than the corresponding value in α-Al 2 O 3 . This suggests that local structural disorder does not significantly affect the energy landscape and diffusion behavior of H i +1 in Al 2 O 3 . Overall, these results show promise for the application of alumina polymorphs as hydrogen permeation barriers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Formation of I 1 stacking fault by deformation defect evolution from grain boundaries in Mg

I 1 stacking faults (SFs) in Mg alloys are regarded as the nucleation sites of $\langle \text{c+a} \rangle$ dislocations that are critical for these alloys to achieve high ductility. Previously it was proposed that the formation of I 1 SFs requires the accumulations of a large number of vacancies, which are difficult to achieve at low temperatures. In this study, molecular dynamics (MD) and molecular statics (MS) simulations based on empirical interatomic potentials were applied to investigate the deformation defect evolutions from the symmetric tilt grain boundaries (GBs) in Mg and Mg-Y alloys under external loading along $\langle c \rangle$-axis. The results show the planar faults (PFs) on Pyramidal I planes first appear due to the nucleation and glide of $\langle {\frac {1}{2}}c + p\rangle$ partial dislocations from GBs, where $\langle p\rangle$ = ${\frac {1}{3}} \langle10\bar{1}0\rangle$. These partial dislocations with pyramidal PFs interact with other defects, including pyramidal PFs themselves, GBs, and $\langle p\rangle$ partial dislocations, generating a large amount of I 1 SFs. Detailed analyses show the nucleation and growth of I 1 SFs are achieved by atomic shuffle events and deformation defect reactions without the requirements of vacancy diffusion. Our simulations also suggest the Y clusters at GBs can reduce the critical stress for the formation of pyramidal PFs and I 1 SFs, which provide a possible reason for the experimental observations that Y promotes the $\langle{\text{c+a}}\rangle$ dislocation activities.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enhanced mobility of cations and anions in the redox state: The polaronium mechanism

Hugely enhanced slow-ion diffusivity has been widely observed under extreme redox conditions and for unclear reasons. Aided by first-principles calculations on model systems of ZrO 2 , CeO 2 , BaTiO 3 and Li 4/3 Mn 2/3 O 2 , here we successfully explained the intriguing phenomenon by a polaronium mechanism. We found a polaronium, defined as a transitory complex of a polaronic electron or hole and a migrating counterion, becomes highly mobile when the counterion comes from a d 0 or f 0 cation (e.g., Zr 4+ , Ti 4+ and Ce 4 ) or a p 6 anion (e.g., O 2- ) in the host compound. Upon a redox reaction, the complex attains a d 1 /f 1 or p 5 configuration, which spontaneously forms because it is favored by an electron-phonon interaction (manifest as the Jahn-Teller effect in high symmetry systems) that enables local relaxation and lowers the system energy. Here our calculations found such interaction reaching its peak at the saddle point where the local environment is softest, so soft that it allows a reorientation of the anisotropic d/f/p orbital to minimize the electron repulsion locally. Since the complex may dissolve after a successful ion-migration event, the redox electron/hole can be recycled to form another free-radical-like polaronium elsewhere, thereby enhancing ion migration repeatedly. The proposed polaronium mechanism, which also operates in ceramics doped with mixed-valence cations, is most relevant under dynamic and extreme thermal/field/irradiation conditions where extra electrons/holes are abundantly generated by non-equilibrium redox reactions. For such operations, some with emerging applications, our diffusion-enhancing mechanism may provide new theoretical insight to help understand their material/microstructure stability and performance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Few‐Unit‐Cell MFI Zeolite Synthesized using a Simple Di‐quaternary Ammonium Structure‐Directing Agent

Abstract Synthesis of a pentasil‐type zeolite with ultra‐small few‐unit‐cell crystalline domains, which we call FDP (few‐unit‐cell crystalline domain pentasil), is reported. FDP is made using bis‐1,5(tributyl ammonium) pentamethylene cations as structure directing agent (SDA). This di‐quaternary ammonium SDA combines butyl ammonium, in place of the one commonly used for MFI synthesis, propyl ammonium, and a five‐carbon nitrogen‐connecting chain, in place of the six‐carbon connecting chain SDAs that are known to fit well within the MFI pores. X‐ray diffraction analysis and electron microscopy imaging of FDP indicate ca. 10 nm crystalline domains organized in hierarchical micro‐/meso‐porous aggregates exhibiting mesoscopic order with an aggregate particle size up to ca. 5 μm. Al and Sn can be incorporated into the FDP zeolite framework to produce active and selective methanol‐to‐hydrocarbon and glucose isomerization catalysts, respectively.

Lu, Peng↗
DOE OSTI2021

Few‐Unit‐Cell MFI Zeolite Synthesized using a Simple Di‐quaternary Ammonium Structure‐Directing Agent

Abstract Synthesis of a pentasil‐type zeolite with ultra‐small few‐unit‐cell crystalline domains, which we call FDP (few‐unit‐cell crystalline domain pentasil), is reported. FDP is made using bis‐1,5(tributyl ammonium) pentamethylene cations as structure directing agent (SDA). This di‐quaternary ammonium SDA combines butyl ammonium, in place of the one commonly used for MFI synthesis, propyl ammonium, and a five‐carbon nitrogen‐connecting chain, in place of the six‐carbon connecting chain SDAs that are known to fit well within the MFI pores. X‐ray diffraction analysis and electron microscopy imaging of FDP indicate ca. 10 nm crystalline domains organized in hierarchical micro‐/meso‐porous aggregates exhibiting mesoscopic order with an aggregate particle size up to ca. 5 μm. Al and Sn can be incorporated into the FDP zeolite framework to produce active and selective methanol‐to‐hydrocarbon and glucose isomerization catalysts, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dislocation cross-slip in precipitation hardened Mg–Nd alloys

An integrated experimental characterization and molecular dynamics (MD) simulation approach was used to explore the dislocation-precipitate interactions in a dilute Magnesium–Neodymium (Mg–Nd) precipitation hardened alloy. In situ indentation in a transmission electron microscope (TEM) and postmortem TEM characterization of deformed samples and MD simulations revealed that basal < a > type dislocations interacted differently with β 1 (Mg 3 Nd) and β"(Mg 3-7 Nd) precipitates. For β" precipitates, the basal dislocations directly shear the precipitates. For β 1 precipitates, such shearing becomes much more difficult because it requires the creation of antiphase boundaries in the ordered lattice of β 1 precipitates. Screw dislocations were observed to cross-slip from basal to the prismatic plane, which could be parallel to the broad facet of β 1 precipitates. It is postulated that double cross-slip (basal to prismatic to basal) via the Hirsch mechanism may enable screw dislocations to overcome the β 1 precipitates. MD simulations also revealed that an edge dislocation is unable to bypass the precipitate via the generation of screw dislocation segments plus the double cross-slip mechanism. The edge dislocation can bow around the precipitate, and the segment that adopts a screw character can cross-slip to the prismatic plane but with increasing applied stress, cross-slips back to the original basal plane to continue glide via precipitate shearing. Finally, the implication of glide dislocation - β 1 -precipitate interaction mechanisms on the strength and ductility of β 1 precipitate dominant microstructures is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Irradiation-Induced Extremes Create Hierarchical Face-/Body-Centered-Cubic Phases in Nanostructured High Entropy Alloys

A nanoscale hierarchical dual-phase structure is reported to form in a nanocrystalline NiFeCoCrCu high-entropy-alloy (HEA) film via ion irradiation. Under the extreme energy deposition and consequent thermal energy dissipation induced by energetic particles, a fundamentally new phenomenon is revealed, in which the original single-phase face-centered-cubic (FCC) structure partially transforms into alternating nanometer layers of a body-centered-cubic (BCC) structure. The orientation relationship follows the Nishiyama–Wasser-man relationship, that is, (011) BCC || (1¯1¯1) FCC and [100] BCC || [11¯0] FCC . Simulation results indicate that Cr, as a BCC stabilizing element, exhibits a tendency to segregate to the stacking faults (SFs). Furthermore, the high densities of SFs and twin boundaries in each nanocrystalline grain serve to accelerate the nucleation and growth of the BCC phase during irradiation. By adjusting the irradiation parameters, desired thicknesses of the FCC and BCC phases in the laminates can be achieved. Finally, this work demonstrates the controlled formation of an attractive dual-phase nanolaminate structure under ion irradiation and provides a strategy for designing new derivate structures of HEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Grain boundary structure search by using an evolutionary algorithm with effective mutation methods

Grain boundaries (GBs) accommodate the misorientation between adjacent grains in a polycrystalline material. GBs are geometrically described by the macroscopic and microscopic degrees of freedom. Besides, at the atomistic level, GBs exhibit complicated behaviors under varying thermodynamic conditions. The complexity of atomistic GB structures demands stochastic searching for possible states. The effectiveness of stochastic search methods relies on techniques to recreate and select atomistic structures. In this work, we developed a new mutation operator that can induce direct and collective atomistic structure changes to boost the search efficiency of exploring GB structures with evolutionary algorithms (EA). We implemented the mutation methods along with innovative selection, crossover, boundary condition preprocessing methods to form an EA-based package to explore GB structures in grand canonical ensembles with atomistic simulations. We used this package to study the [001] symmetric tilt grain boundaries (STGBs) in FCC copper (Cu), the [110] STGBs in BCC tungsten (W), and the $[1\bar{2}10]$ STGBs in HCP magnesium (Mg). The results show that our design and implementation based on new mutation procedures, selection, and boundary conditions provide a high-quality search of atomistic GB structures in the grand canonical ensemble for different crystal lattices.

36 MATERIALS SCIENCE↗