Researcher profile

Zi-Kui Liu

Zi-Kui Liu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2025arXiv

ZENN: A Thermodynamics-Inspired Computational Framework for Heterogeneous Data-Driven Modeling

Traditional entropy-based methods - such as cross-entropy loss in classification problems - have long been essential tools for representing the information uncertainty and physical disorder in data and for developing artificial intelligence algorithms. However, the rapid growth of data across various domains has introduced new challenges, particularly the integration of heterogeneous datasets with intrinsic disparities. To address this, we introduce a zentropy-enhanced neural network (ZENN), extending zentropy theory into the data science domain via intrinsic entropy, enabling more effective learning from heterogeneous data sources. ZENN simultaneously learns both energy and intrinsic entropy components, capturing the underlying structure of multi-source data. To support this, we redesign the neural network architecture to better reflect the intrinsic properties and variability inherent in diverse datasets. We demonstrate the effectiveness of ZENN on classification tasks and energy landscape reconstructions, showing its superior generalization capabilities and robustness-particularly in predicting high-order derivatives. ZENN demonstrates superior generalization by introducing a learnable temperature variable that models latent multi-source heterogeneity, allowing it to surpass state-of-the-art models on CIFAR-10/100, BBCNews, and AGNews. As a practical application in materials science, we employ ZENN to reconstruct the Helmholtz energy landscape of Fe$_3$Pt using data generated from density functional theory (DFT) and capture key material behaviors, including negative thermal expansion and the critical point in the temperature-pressure space. Overall, this work presents a zentropy-grounded framework for data-driven machine learning, positioning ZENN as a versatile and robust approach for scientific problems involving complex, heterogeneous datasets.

preprint2022arXiv

A Solution to the Temperature Evolution of Multi-well Free-energy Landscape

It has been a grand challenge to resolve the temperature evolution of multi-well free-energy landscape which is fundamentally relevant to phase transitions and associated critical phenomena as listed by Ginzburg [Rev. Mod. Phys. 76 (2004) 981]. To address this challenge, here we provide a simple solution based on a priori concept of Boltzmann thermal mixing among multiple parabolic potentials. The success and the impact of the present approach have been extensively demonstrated using a variety of materials, including Nb, YBa2Cu3O(7-x), Ca3Ru2O7, Ni, BiFeO3, and PbTiO3.

preprint2022arXiv

Insight into Ideal Shear Strength of Ni-based Dilute Alloys using First-Principles Calculations and Correlational Analysis

The present work examines the effect of alloying elements (denoted X) on the ideal shear strength for 26 dilute Ni-based alloys, Ni$_{11}$X, as determined by first-principles calculations of pure alias shear deformations. The variations in ideal shear strength are quantitatively explored with correlational analysis techniques, showing the importance of atomic properties such as size and electronegativity. The shear moduli of the alloys are affirmed to show a strong linear relationship with their ideal shear strengths, while the shear moduli of the individual alloying elements were not indicative of alloy shear strength. Through combination with available ideal shear strength data on Mg alloys, a potential application of the Ni alloy data is demonstrated in the search for a set of atomic features suitable for machine learning applications to mechanical properties. As another illustration, the predicted Ni ideal shear strengths play a key role in a predictive multiscale framework for deformation behavior of single crystal alloys at large strains, as shown by simulated stress-strain curves.

preprint2022arXiv

Predictive Crystal Plasticity Modeling of Single Crystal Nickel Based on First-Principles Calculations

To reduce reliance on experimental fitting data within the crystal plasticity finite element method (CPFEM), an approached is proposed that integrates first-principles calculations based on density functional theory (DFT) to predict the strain hardening behavior of pure Ni single crystals. Flow resistance was evaluated through the Peierls-Nabarro equation using the ideal shear strength and elastic properties calculated by DFT-based methods, with hardening behavior modeled by imposing strains on supercells in first-principles calculations. Considered alone, elastic interactions of pure edge dislocations capture hardening behavior for small strains on single slip systems. For larger strains, hardening is captured through a strain-weighted linear combination of edge and screw flow resistance components. The rate of combination is not predicted in the present framework, but agreement with experiments through large strains (~0.4) for multiple loading orientations demonstrates a possible route for more predictive crystal plasticity modeling through incorporation of analytical models of mesoscale physics.

preprint2022arXiv

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

The Nb-Ni system has been remodeled with uncertainty quantification (UQ) by using the presently upgraded software tools of PyCalphad and ESPEI that contain the new capability to model site occupancy of Wyckoff position for the phases of interest. Specifically, the five- and three-sublattice models are used to model the topologically close pack (TCP) phases of μ-Nb7Ni6 and δ-NbNi3, respectively, according to exactly their Wyckoff positions; where the inputs for CALPHAD-based modeling include the presently predicted thermochemical data as a function of temperature by density functional theory (DFT) based first-principles and phonon calculations together with both phase equilibrium and site occupancy data in the literature. Besides phase diagram and thermodynamic properties, the present CALPHAD predictions of site occupancies are also agreed well with experimental data such as the measured Nb sites in μ-Nb7Ni6. In addition, the predicted UQ values using the Markov Chain Monte Carlo (MCMC) method as implemented in ESPEI make it possible to quantify uncertainties in the Nb-Ni system, such as site occupancies in μ-Nb7Ni6 and enthalpy of mixing in liquid.

preprint2022arXiv

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 in complex liquids such as molten salts, slag in metal processing, and electrolytic solutions. 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) 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-NiF2 system as an example, 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 using ESPEI. The present implementation offers an open-source capability for performing CALPHAD modeling for complex liquids with short-range ordering using MQMQA.

preprint2021arXiv

Generative deep learning as a tool for inverse design of high-entropy refractory alloys

Generative deep learning is powering a wave of new innovations in materials design. In this article, we discuss the basic operating principles of these methods and their advantages over rational design through the lens of a case study on refractory high-entropy alloys for ultra-high-temperature applications. We present our computational infrastructure and workflow for the inverse design of new alloys powered by these methods. Our preliminary results show that generative models can learn complex relationships in order to generate novelty on demand, making them a valuable tool for materials informatics.

preprint2020arXiv

An ab initio molecular dynamics exploration of associates in Ba-Bi liquid with strong ordering trends

Fictive associates are widely used to describe and model liquid phases with strong ordering trends. However, little evidence is known about the assumed associates in most cases. In the present work, an ab initio molecular dynamics (AIMD) study is employed to investigate the characters of the Ba-Bi liquid, in which associates have been assumed in existing thermodynamic modeling. It is found that in the Ba rich melt, the Bi atoms are almost completely surrounded by Ba atoms. The Bi-centered coordination polyhedrons are strongly associated to crystalline structures of Ba5Bi3 and Ba4Bi3 with a longer lifetime than other polyhedrons during the AIMD simulations. In addition, these Bi-centered polyhedrons in Ba rich melt connect with each other through vertex, edge, face, and/or bipyramid sharing to form medium range orders (MRO). In the Bi rich melt, the Ba-centered polyhedrons also form MROs, but they are both structurally and compositionally diverse with a shorter lifetime. These findings from AIMD study provide evidences that there exist a strongly ordering Ba4Bi3 associate and a weakly ordering BaBi3 associate in the Ba-Bi liquid. The predicted enthalpy of mixing in the liquid agrees well with the results by the CALPHAD modeling in the literature.

preprint2020arXiv

Route to in situ synthesis of epitaxial Pr2Ir2O7 thin films guided by thermodynamic calculations

In situ growth of pyrochlore iridate thin films has been a long-standing challenge due to the low reactivity of Ir at low temperatures and the vaporization of volatile gas species such as IrO3(g) and IrO2(g) at high temperatures and high oxygen partial pressures. To address this challenge, we combine thermodynamic analysis of the Pr-Ir-O2 system with experimental results from the conventional physical vapor deposition (PVD) technique of co-sputtering. Our results indicate that only high growth temperatures yield films with crystallinity sufficient for utilizing and tailoring the desired topological electronic properties. Thermodynamic calculations indicate that high deposition temperatures and high partial pressures of gas species O2(g) and IrO3(g), are required to stabilize Pr2Ir2O7. We further find that the gas species partial pressure requirements are beyond that achievable by any conventional PVD technique. We experimentally show that conventional PVD growth parameters produce exclusively Pr3IrO7, which conclusion we reproduce with theoretical calculations. Our findings provide solid evidence that in situ synthesis of Pr2Ir2O7 thin films is fettered by the inability to grow with oxygen partial pressure on the order of 10 Torr, a limitation inherent to the PVD process. Thus, we suggest high-pressure techniques, in particular chemical vapor deposition (CVD), as a route to synthesis of Pr2Ir2O7, as this can support thin film deposition under the high pressure needed for in situ stabilization of Pr2Ir2O7.

preprint2020arXiv

Scheil-Gulliver simulations for the design of functionally graded alloys by additive manufacturing using pycalphad

Additive manufacturing (AM), through directed energy deposition, supports planned composition changes between locations within a single component, allowing for functionally graded materials (FGMs) to be developed and fabricated. The formation of deleterious phases along a particular composition path can cause significant cracking during the AM build process that makes the composition path unviable to produce these FGMs, but it is challenging to predict which phases will be present in as-built additively manufactured parts by analyzing only equilibrium phase relations. Solute segregation during solidification can lead to the formation of non-equilibrium phases that are stable at compositions far from the nominal composition of the melt, leading to crack formation. In this work, we developed a Scheil-Gulliver simulation tool based on pycalphad. We used this tool to compare the non-equilibrium phases predicted to form during solidification using the Scheil-Gulliver model with experimentally measured phases at several locations with different composition in a Ti-6Al-4V to Invar-36 FGM and a commercially pure Ti to Invar-36 FGM. We showed that the phases predicted to form by the Scheil-Gulliver model outperform the predictions made by assuming equilibrium solidification. Further, we demonstrated the use of our Scheil-Gulliver simulation tool as a method of screening potential FGM pathways by calculating the solidification phase fractions along the experimental gradient path in composition space.

preprint2020arXiv

Selection of strain and fitting schemes for calculating higher-order elastic constants

Criteria of selecting strain and fitting schemes are proposed for the calculation of higher-order elastic constants more efficiently, robustly and accurately. As demonstrated by the third-order elastic constants (TOECs) of diamond, the proposed method is 3-5 times faster than existing methods, and the range of strain for getting correct TOECs is expanded. In addition, our result provides an evidence for the inaccuracy of some previous experiments caused by higher-order effect, and the difference among experiments and several different theoretical methods is resolved. Finally, we give the recommend TOECs values for diamond.

preprint2020arXiv

Sensitivity estimation for calculated phase equilibria

The development of a consistent framework for Calphad model sensitivity is necessary for the rational reduction of uncertainty via new models and experiments. In the present work, a sensitivity theory for Calphad was developed, and a closed-form expression for the log-likelihood gradient and Hessian of a multi-phase equilibrium measurement was presented. The inherent locality of the defined sensitivity metric was mitigated through the use of Monte Carlo averaging. A case study of the Cr-Ni system was used to demonstrate visualizations and analyses enabled by the developed theory. Criteria based on the classical Cramér-Rao bound were shown to be a useful diagnostic in assessing the accuracy of parameter covariance estimates from Markov Chain Monte Carlo. The developed sensitivity framework was applied to estimate the statistical value of phase equilibria measurements in comparison with thermochemical measurements, with implications for Calphad model uncertainty reduction.

preprint2019arXiv

Quantified Uncertainty in Thermodynamic Modeling for Materials Design

Phase fractions, compositions and energies of the stable phases as a function of macroscopic composition, temperature, and pressure (X-T-P) are the principle correlations needed for the design of new materials and improvement of existing materials. They are the outcomes of thermodynamic modeling based on the CALculation of PHAse Diagrams (CALPHAD) approach. The accuracy of CALPHAD predictions vary widely in X-T-P space due to experimental error, model inadequacy and unequal data coverage. In response, researchers have developed frameworks to quantify the uncertainty of thermodynamic property model parameters and propagate it to phase diagram predictions. In previous studies, uncertainty was represented as intervals on phase boundaries (with respect to composition) or invariant reactions (with respect to temperature) and was unable to represent the uncertainty in eutectoid reactions or in the stability of phase regions. In this work, we propose a suite of tools that leverages samples from the multivariate model parameter distribution to represent uncertainty in forms that surpass previous limitations and are well suited to materials design. These representations include the distribution of phase diagrams and their features, as well as the dependence of phase stability and the distributions of phase fraction, composition activity and Gibbs energy on X-T-P location - irrespective of the total number of components. Most critically, the new methodology allows the material designer to interrogate a certain composition and temperature domain and get in return the probability of different phases to be stable, which can positively impact materials design.

preprint2019arXiv

Realization of Epitaxial Thin Films of the Topological Crystalline Insulator Sr$_3$SnO

Topological materials are derived from the interplay between symmetry and topology. Advances in topological band theories have led to the prediction that the antiperovskite oxide Sr$_3$SnO is a topological crystalline insulator, a new electronic phase of matter where the conductivity in its (001) crystallographic planes is protected by crystallographic point group symmetries. Realization of this material, however, is challenging. Guided by thermodynamic calculations we design and implement a deposition approach to achieve the adsorption-controlled growth of epitaxial Sr$_3$SnO single-crystal films by molecular-beam epitaxy (MBE). In-situ transport and angle-resolved photoemission spectroscopy measurements reveal the metallic and non-trivial topological nature of the as-grown samples. Compared with conventional MBE, the synthesis route used results in superior sample quality and is readily adapted to other topological systems with antiperovskite structures. The successful realization of thin films of topological crystalline insulators opens opportunities to manipulate topological states by tuning symmetries via epitaxial strain and heterostructuring.