Engineering Papers⌕ Search

Engineering topics

Shang, Shun-Li

Publications and source records attributed to Shang, Shun-Li.

Synthesis of epitaxial TaO$_2$ thin films on Al$_2$O$_3$ by suboxide molecular-beam epitaxy and thermal laser epitaxy

Tantalum dioxide (TaO2) is a metastable tantalum compound. Here, we report the epitaxial stabilization of TaO2 on Al2O3 (1-102) (r-plane sapphire) substrates using suboxide molecular-beam epitaxy (MBE) and thermal laser epitaxy (TLE), demonstrating single-oriented, monodomain growth of anisotropically strained thin films. Microstructural investigation is performed using synchrotron X-ray diffraction and scanning transmission electron microscopy. The tetravalent oxidation state of tantalum is confirmed using X-ray absorption and photoemission spectroscopy as well as electron energy-loss spectroscopy. Optical properties are investigated via spectroscopic ellipsometry and reveal a 0.3 eV Mott gap of the tantalum 5d electrons. Density-functional theory and group theoretical arguments are used to evaluate the limited stability of the rutile phase and reveal the potential to unlock a hidden metal-insulator transition concomitant with a structural phase transition to a distorted rutile phase, akin to NbO2. Our work expands the understanding of tantalum oxides and paves the way for their integration into next-generation electronic and photonic devices.

FOS: Physical sciences↗

Predicting phase transitions in PbTi⁢O 3 using zentropy through quasiharmonic phonon calculations

According to x-ray diffraction (XRD) measurements, PbTi⁢O 3 undergoes a phase transition from a tetragonal ferroelectric (FE) phase to a cubic paraelectric phase at 763 K. However, x-ray absorption fine-structure (XAFS) measurements indicate that PbTi⁢O 3 is locally tetragonal even after the phase transition. The difference in these results is because XAFS measurements can probe local features of a structure, while XRD averages over such local features. For both measurements to be consistent, PbTi⁢O 3 is macroscopically cubic but locally tetragonal after the phase transition. Despite this, most models, such as the Landau-Ginsburg-Devonshire theory and effective Hamiltonians, are still unable to explain this phenomenon. Moreover, these methods involve model parameters fitted to experimental or theoretical data and do not consider other tetragonal configurations, such as domain walls, to predict the phase transition. In our previous study, we used our zentropy approach to predict the phase transition by considering the tetragonal FE ground-state configuration and the tetragonal 90° and 180° domain wall configurations with their total energies at 0 K. Here, in this paper, the Helmholtz energies of the three configurations are obtained from density functional theory calculations through energy-volume curves and phonon calculations. The predicted phase transition temperature using the meta-GGA 𝑟 2⁢ SCAN and revised multiplicities of configurations is 716 K, showing good agreement with the experimental value of 763 K.

36 MATERIALS SCIENCE↗

MaterialsMap: A CALPHAD-based tool to design composition pathways through feasibility map for desired dissimilar materials, demonstrated with resistance spot welding joining of Ag-Al-Cu

Assembly of dissimilar metals can be achieved by different methods, for example, casting, welding, and additive manufacturing (AM). However, undesired phases formed in liquid-phase assembling processes due to solute segregation during solidification diminish mechanical and other properties of the processed parts. In the present work, an open-source software named MaterialsMap, has been developed based on the CALculation of Phase Diagrams (CALPHAD) approach. The primary objective of MaterialsMap is to facilitate the design of an optimal composition pathway for assembling dissimilar alloys with liquid-phases based on the formation of desired and undesired phases along the pathway. In MaterialsMap, equilibrium thermodynamic calculations are used to predict equilibrium phases formed at slow cooling rate, while Scheil-Gulliver simulations are employed to predict non-equilibrium phases formed during rapid cooling. Additionally, by combining these two simulations, MaterialsMap offers a thorough guide for understanding phase formation in various manufacturing processes, assisting users in making informed decisions during material selection and production. As a demonstration of this approach, a compositional pathway was designed from pure Al to pure Cu through Ag using MaterialsMap. The design was experimentally verified using resistance spot welding (RSW).

36 MATERIALS SCIENCE↗

Electronic-grade epitaxial (111) KTaO 3 heterostructures

KTaO 3 heterostructures have recently attracted attention as model systems to study the interplay of quantum paraelectricity, spin-orbit coupling, and superconductivity. However, the high and low vapor pressures of potassium and tantalum present processing challenges to creating heterostructure interfaces clean enough to reveal the intrinsic quantum properties. Here, we report superconducting heterostructures based on high-quality epitaxial (111) KTaO 3 thin films using an adsorption-controlled hybrid PLD to overcome the vapor pressure mismatch. Electrical and structural characterizations reveal that the higher-quality heterostructure interface between amorphous LaAlO 3 and KTaO 3 thin films supports a two-dimensional electron gas with substantially higher electron mobility, superconducting transition temperature, and critical current density than that in bulk single-crystal KTaO 3 -based heterostructures. Our hybrid approach may enable epitaxial growth of other alkali metal–based oxides that lie beyond the capabilities of conventional methods.

36 MATERIALS SCIENCE↗

Additively manufactured Ni-20Cr to V functionally graded material: Computational predictions and experimental verification of phase formations

Throughout this work, a database for the Cr-Ni-V system was constructed by modeling the binary Cr-V and ternary Cr-Ni-V systems using the CALPHAD approach aided by density functional theory (DFT)-based first-principles calculations and ab initio molecular dynamics (AIMD) simulations. To validate this new database, a functionally graded material (FGM) using Ni-20Cr and elemental V was fabricated using directed energy deposition additive manufacturing (DED AM) and experimentally characterized. The deposited Ni-20Cr was pure fcc phase, while increasing the amount of V across the gradient resulted in the formation of sigma phase, followed by the bcc phase. The experimentally measured phase data was compared with computational predictions made using a Cr-Ni-V thermodynamic database from the literature as well as the database developed in the present work. The newly developed database was shown to better predict the experimentally observed phases due to its accurate modeling of binary systems within the database and the ternary liquid phase, which is critical for accurate Scheil calculations.

36 MATERIALS SCIENCE↗

Zentropy theory for accurate prediction of free energy, volume, and thermal expansion without fitting parameters

Based on statistical mechanics, a macroscopically homogeneous system, i.e., a single phase in the present context, is composed of many independent configurations that the system embraces. The macroscopical properties of the system are determined by the properties and statistical probabilities of those configurations with respect to external conditions. The volume of a single phase is thus the weighted sum of the volumes of all configurations. Consequently, the derivative of the volume to temperature of a single phase depends on both the derivatives of the volumes of every configuration to temperature and the derivatives of their statistical probabilities to temperature, with the latter introducing nonlinear emergent behaviors. It is shown that the derivative of the volume to the temperature of the single phase can be negative, i.e., negative thermal expansion, due to the symmetry-breaking non-ground-state configurations with smaller volumes than that of the ground-state configuration and the rapid increase of the statistical probabilities of the former, and negative thermal expansion can be predicted without fitting parameters from the zentropy theory that combines quantum mechanics and statistical mechanics with the free energy of each configuration predicted from quantum mechanics and the partition function of each configuration calculated from its free energy.

36 MATERIALS SCIENCE↗

Thermodynamic modeling with uncertainty quantification using the modified quasichemical model in quadruplet approximation: Implementation into PyCalphad and ESPEI

The modified quasichemical model in the quadruplet approximation (MQMQA) considers the first- and the second-nearest-neighbor coordination and interactions, particularly useful in describing short-range ordering (SRO) in complex liquids such as molten salts, slag in metal processing, and electrolytic solutions. Here, the present work implements the MQMQA into the Python based open-source software PyCalphad for thermodynamic calculations. This endeavor facilitates the development of MQMQA-based thermodynamic database with uncertainty quantification (UQ) and propagation (UP) using the open-source software ESPEI. A new database structure based on Extensible Markup Language (XML) is proposed for ESPEI evaluation of MQMQA model parameters. Using the KF-NiF 2 , KCl-NaCl-MgCl 2 , and CaCl 2 -CaF 2 -LiCl-LiF salt systems as examples, we demonstrate the successful implementation of MQMQA in PyCalphad through thermodynamic calculations of Gibbs energy, equilibrium quadruplet fractions, and phase diagram, as well as database development with UQ and UP using ESPEI. Furthermore, as an application of the present implementation, both the LiF–TbF 3 and LiF-HoF 3 systems have been modeled by MQMQA for the first time, which are in good agreement with experiments. The present implementation hence offers an open-source capability for performing CALPHAD modeling for complex liquids with SRO using MQMQA plus a new XML database structure.

36 MATERIALS SCIENCE↗

Determination of Kinetic Properties of Ni(II) Ions in Molten LiF-NaF-KF via Voltammetry

Kinetic properties of Ni(II) in eutectic LiF-NaF-KF (FLiNaK) molten salt were determined at T = 748–823 K using cyclic voltammetry (CV), square wave voltammetry (SWV), and chronoamperometry (CA) measurements using a glassy C working electrode, Ni(II)/Ni reference electrode, and Ni counter electrode. Reduction of Ni(II) to Ni(s) was determined to be a single step, two-electron transfer process. Diffusivity values were calculated using the Berzins and Delahay equation and semi-integral electroanalysis from the CV measurements as well as using the Cottrell equation from the CA measurements. Diffusivity of Ni(II) in molten FLiNaK at T = 748–823 K was determined to be 3.27 × 10 –7 –3.04 × 10 –6 cm 2 s –1 with an activation energy of 62–104 kJ mol –1 . The estimated kinetic properties varied appreciably among methodologies possibly due to inherent assumptions in theory regarding reversibility and unit activity during Ni metal deposition on the working electrode.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamic modeling of the Nb-Ni system with uncertainty quantification using PyCalphad and ESPEI

Here, the Nb–Ni system is remodeled with uncertainty quantification (UQ) using software tools of PyCalphad and ESPEI (the Extensible, Self-optimizing Phase Equilibria Infrastructure) with the presently implemented capability of modeling site fraction based on Wyckoff positions. The five- and three-sublattice models are used to model the topologically close pack (TCP) μ-Nb 7 Ni 6 and δ-NbNi 3 phases according to their Wyckoff positions. The inputs for CALPHAD-based thermodynamic modeling include the thermochemical data as a function of temperature predicted by first-principles and phonon calculations based on density functional theory (DFT), ab initio molecular dynamics (AIMD) simulations, together with phase equilibrium and site fraction data in the literature. In addition to phase diagram and thermodynamic properties, the CALPHAD-based predictions of site fractions of Nb in μ-Nb 7 Ni 6 agree well with experimental data. Furthermore, the UQ estimation using the Markov Chain Monte Carlo (MCMC) method as implemented in ESPEI is applied to study the uncertainty of site fraction in μ-Nb 7 Ni 6 and enthalpy of mixing (ΔH mix ) in liquid.

36 MATERIALS SCIENCE↗

Parameter-free prediction of phase transition in PbTiO 3 through combination of quantum mechanics and statistical mechanics

Thermodynamics of ferroelectric materials and their ferroelectric to paraelectric (FE-PE) transitions is commonly described by the phenomenological Landau theory and more recently by effective Hamiltonian and various potentials, all with model parameters fitted to experimental or theoretical data. Here, in this work, we show that the zentropy theory, which considers the total entropy of a system as a weighted sum of entropies of configurations that the system may experience and the statistical entropy among the configurations, can predict the FE-PE transition without fitting parameters. For PbTiO 3 , the configurations are identified as the FE configurations with 90- or 180° domain walls in addition to the ground state FE configuration without domain wall. With the domain wall energies predicted from first-principles based on density functional theory in the literature as the only inputs, the FE-PE transition for PbTiO 3 is predicted showing remarkable agreement with experiments, unveiling the microscopic fundamentals of the transition.

36 MATERIALS SCIENCE↗