Abstract
The accelerated discovery of high-performance lead-free piezoelectric ceramics is hindered by the vast compositional space and the limited interpretability of conventional machine learning (ML) models. Here, we propose a physics-informed and interpretable ML framework with integrated uncertainty quantification to predict and understand the piezoelectric coefficient d(33) of (K(0.5)Na(0.5)) NbO(3) (KNN)-based ceramics. A curated dataset of 1113 experimental samples is used to construct 65 descriptors by decoupling A-site and B-site ionic contributions. Pearson correlation analysis reduces these to an optimized 11-dimensional feature set for training deep neural networks, Wide & Deep networks, and residual networks. A Bayesian neural network further provides predictive uncertainty, which quantitatively reflects the confidence of machine-learning-based d(33) predictions rather than experimental measurement uncertainty. To achieve physical interpretability, SHapley Additive exPlanations (SHAP) are combined with the Sure Independence Screening and Sparsifying Operator (SISSO) to derive a compact analytical descriptor revealing that sintering temperature, B-site electronic anisotropy, and A-site ionic displacement jointly govern d(33). The proposed framework achieves high accuracy (R(2) ≈ 0.81) while offering transparent design rules for next-generation lead-free piezoelectrics.