2026-07-30 | | Total: 33
Materials behavior is often treated as a deterministic mapping from structure to properties, yet many important phenomena emerge from the conditional activation of multiple mechanisms across scales. This is especially evident in fatigue of metals, where crack growth is typically modeled as monotonic and irreversible process, despite evidence that local microstructure, loading history, and competing unit processes can shift the balance among propagation, arrest, and self-healing. Here we present a probabilistic framework that describes materials behavior as an ensemble of constituent mechanisms whose activation, interaction, and evolution determine emergent outcomes. The framework connects mechanism activation, state evolution, and macroscopic observables in a probabilistic way. In the case of fatigue crack propagation, it reframes damage tolerance as an inference problem over mechanism competition and provides a basis for integrating multiscale simulation, multimodal characterization, and machine learning. The same logic extends to other physical and chemical systems suggesting a portable framework for any system in which emergent behavior reflects mechanism competition under changing conditions. The broader ambition of this perspective review is a shift from correlating structure and performance after the fact to identifying, in advance, the conditions that make desired emergent behavior probable.
Precipitation-strengthened Cu-Ni-Al alloys are of interest for technological applications because coherent, $\mathrm{L}1_2$-ordered $γ'$-$\mathrm{Ni}_3\mathrm{Al}$ precipitates can confer high mechanical strength while allowing the material to retain many of the good transport properties characteristic of elemental Cu. In this work, we study the thermodynamics and phase stability of the pseudobinary $\textrm{Cu}_x (\textrm{Ni}_{3/4} \textrm{Al}_{1/4})_{1-x}$ system, $0 \leq x \leq 1$. We use a computational modelling framework combining first-principles electronic structure calculations with a concentration wave analysis from which atom-atom effective pair interactions are extracted for use in atomistic Monte Carlo simulations. Our modelling reveals three distinct, composition-dependent regimes of phase behaviour, in qualitative agreement with the experimentally determined phase diagram. At low Cu content, Cu is soluble in the $\mathrm{L}1_2$-ordered $\mathrm{Ni}_3\mathrm{Al}$ phase, with a single identifiable phase transition corresponding to chemical ordering between Ni and Al. At intermediate compositions, this high-temperature ordering is followed at lower temperatures by phase separation of Cu and $\mathrm{L}1_2$-ordered $\mathrm{Ni}_3\mathrm{Al}$. Finally, at high Cu content, $\mathrm{L}1_2$-ordered $\mathrm{Ni}_3\mathrm{Al}$ precipitates directly from the solid solution, with no clearly identifiable secondary transition. We relate these phase transformations to features of the underlying electronic structures of the considered alloys. Overall, this work demonstrates a computationally efficient workflow capturing both chemical ordering and coherent precipitation in multicomponent substitutional alloys, with relevance to the study of phenomena such as precipitation strengthening.
Density functional theory (DFT) underestimates the experimental band gap---the infamous band gap problem. As the band gap defines the energy scale of defect levels, this complicates computation of charge transition energies for atomic defects. In the extreme case of narrow-gap semiconductors, the DFT band gap collapses to zero, seemingly precluding quantitative predictions of defect levels. We present a band-avoiding occupation-constrained DFT (ba-occ-DFT) approach that prevents spurious occupation of band-edge states and enables reliable total energy calculations of atomic defects. Application to indium arsenide (InAs) shows that ba-occ-DFT circumvents the band gap problem, separates band-edge errors from defect level calculations, and enables rigorous defect level predictions in a narrow-gap semiconductor despite a zero DFT band gap.
Although dense sphere packings serve as the initial state for simulations in powder metallurgy, additive manufacturing and granular physics, the choice of a packing generator is rarely guided by a systematic benchmark. A representative packing must be (i)-(ii) physically admissible (non-overlapping particles in gravitational equilibrium); (iii) faithful to the target particle size distribution (PSD); (iv) representative in relative density Phi; and (v) computationally affordable. Four open-source tools have been benchmarked, meeting (i)-(ii) by construction: the sequential DR (dropping-and-rolling) and its densified variant DR-ME, and the discrete element method (DEM) codes LAMMPS (gravity) and dp3D (isostatic compression). Across four configurations (2D/3D lognormal, 3D bimodal, and a 3D domain-size study), they are compared against an industrial MIM-grade powder, with PSD fidelity measured by the bin-width-independent Hellinger distance and Phi against the feedstock solid loading (phi_exp = 0.62, by Archimedes' method). In 3D, the DEM codes reach the densest packings but run more than three orders of magnitude slower: for approximatively 20 000 particles, DR shows a 9% phi shortfall relative to dp3D while running 1800x faster. These idealised model packings yield application-driven tool-selection guidelines.
In magnetic materials angular momentum can be mediated by different carriers, including electrons, magnons, and phonons. The magnons can interact with circularly polarized phonons which are close in energy, provided specific symmetry conditions are met. If the interaction is strong enough the phonons and magnons shift in frequency and start to mix to form hybrid magneto-elastic quasi-particles. In this paper, we develop a constrained Hamiltonian framework which incorporates hybrid stiffness matrices and Berry curvatures. We quantify the degree of hybridization between phonons and magnons (up to 8% in CrI$_3$ and 25% in CrBr$_3$) using a decomposition of the total energy, which is a generalization of the norm decomposition for atomic contributions to phonons used in the literature. We also explore how the total angular momentum is conserved but shared between the phononic and magnonic subsystems upon hybridization.
Unconventional collinear magnets with almost zero magnetization but prominent nonrelativistic spin-splitting, such as altermagnets, can inherit the advantages of both ferromagnets and antiferromagnets. By incorporating more degrees of freedom such as ferroelectricity and charge ordering, these unconventional magnets can be even more interesting and functionalized. With this design principle, the Fe$_3$O$_5$ monolayer is predicted to exhibit a hybrid spin-splitting mechanism, with the superposition of the altermagnetic-like $k$-path alternating splitting and ferrimagnet-like Zeeman splitting. Benefiting from the hidden magnetoelectricity based on the spin-charge coupling, such spin-splitting can be fully switched by an electric field. Its conductivity is highly spin-polarized, with a polarization ratio above $99\%$, comparable to half-metals but with zero magnetization.
High-entropy alloys (HEAs), with their unique compositional-complexity and tunable surface chemistry, have emerged as promising electrocatalysts for energy conversion. The catalytic activity of HEA often arises from the interplay between the intrinsic activity of the individual elements and the synergistic effects generated at the interfaces. Even a single-element substitution in a multicomponent HEA can substantially alter the surface-chemistry and electrochemical kinetics of the active surface. Here we investigated the structure-property relationship of CrMnFeCoNi (HEA-Cr) and MnFeCoNiCu (HEA-Cu) HEAs by comparing the alkaline hydrogen evolution reaction (HER) and oxygen evolution reaction (OER) activity. Under the same experimental condition, HEA-Cu outperforms HEA-Cr towards both HER and OER. Substituting Cr with Cu significantly enhances the bifunctional activity, where HEA-Cu achieved a lower overpotential (538 mV) and Tafel slope (165 mVdec-1) when compared with HEA-Cr. Computational analysis corroborates these findings, showing that Cu substitution modulates the electronic-structure to provide favorable binding energies for reaction intermediates (H*, O*, OH*, and OOH*). Interestingly, unlike HER, recovered HEA-Cu after OER showed migration of Cu forming a Cu-rich outer layer shell with a multimetallic core. These findings demonstrate the potential of single-element substitution in HEAs as a strategy for designing high-performance, cost-effective catalysts for efficient water electrolysis.
Amorphous refractory metal coatings combine high hardness with chemical inertness, yet their metastability makes them prone to mechanically induced crystallisation (devitrification) under contact loading, and how readily such a glass re-orders to its parent body centred cubic (bcc) crystal is unknown across the refractory series. We prepared amorphous V, Nb, Mo, Ta and W by melt quenching to 300 K, validated each interatomic potential against ab initio liquid radial distribution functions, and probed them by large scale molecular dynamics nanoindentation. Indentation drives a localised amorphous to bcc transformation by bulk nucleation, growth and coalescence. A crystal forming ability (CFA), the maximum slope of the sigmoidal bcc fraction versus depth curve, spans about a factor of four and decreases as V > Mo > Nb > Ta > W; it falls with indenter velocity as CFA $\propto v^{-m}$ with a mean exponent of 1.08, so CFA$\cdot v$ is nearly constant and the transformation rate is almost velocity independent. Transforming atoms carry excess non affine displacement and local shear strain, not hydrostatic pressure, marking a shear associated displacive pathway; the bulk driving force is largest for the most resistant element, W, so resistance tracks cohesive bond strength rather than the thermodynamic driving force. Early bcc like order sets the sharpness and depth of the transition, whereas the mechanical work to half transformation sets the persistent nucleus density. A grain population balance links nucleation, growth and coalescence to a terminal microstructure whose completeness does not follow the CFA order.
Hopfions are three-dimensional (3D) topological solitons predicted to exist in diverse magnetic systems, yet their practical utility has been largely restricted to cryogenic environments. Here, we overcome this temperature constraint by demonstrating stable magnetic hopfions in the chiral magnet Co8Zn8Mn4 at and above room temperature. Using a transmission electron microscope equipped for in situ optical excitation, we generate magnetic hopfions with femtosecond laser pulses. Long-term observations further reveal Brownian-like motion at room temperature and thermally activated collapse upon approaching the high-temperature regime. Together with micromagnetic simulations and homotopy group analysis, our experimental observations uncover the hopfion formation mechanism through the fusion of bimeron pairs. These findings establish room-temperature magnetic hopfions and provide a framework for their further studies under technologically relevant conditions.
Magnetic material parameters such as the exchange stiffness and magnetic anisotropy govern the behavior and functionality of magnetic systems, yet their local inference from magnetization data remains challenging, particularly in strongly fluctuating regimes with polycrystalline or multiphase microstructure, where conventional texture-based methods become unreliable. We introduce a magnetization-only framework for inferring material parameters from thermally driven magnetization dynamics. Using micromagnetic simulations, we extract statistical quantities such as temporal mean and latent entropy from the magnetization dynamics, fit models to these descriptors, and invert the models to infer material parameters. We show that this framework enables material-parameter inference as well as grain-boundary detection in a heterogeneous sample. Among the descriptors considered, latent entropy yields more accurate parameter estimates than the temporal mean. Our results establish latent entropy as an efficient descriptor for inferring magnetic material parameters from dynamical magnetization data and point toward its use for experimental parameter extraction at high temperatures and, more broadly, under strongly fluctuating conditions.
Memristors are promising devices for applications such as non-volatile memory, neuromorphic computing, logic circuits, and analog signal processing. The development of such systems requires accurate simulations based on models that reproduce the electrical behavior of real devices under both continuous and pulsed excitation. This work presents the development of a simulation environment for a volatile TiO2-based memristor. Experimental measurement data are analyzed to verify the memristive behavior of the device and to identify a suitable model. The model parameters are then optimized to match the measured characteristics. The resulting model is implemented in SPICE and validated by comparing simulation results with measurement data. The comparison shows a good agreement between simulation and experiment, demonstrating that the developed model is suitable for reproducing the electrical behavior of the investigated memristor and can be applied in circuit-level simulations, as demonstrated by a leaky integrate-and-fire neuron.
In the quest for sustainable energy solutions, water splitting emerges as a crucial process for generating clean hydrogen a versatile and renewable fuel essential for energy storage, emissions reduction, and achieving sustainability goals. This study employs a comprehensive computational approach, utilizing atomistic simulations to systematically investigate the effects of water absorption on the electronic and optical properties of g-C3N4 nanosheets. Our methodology integrates ab initio computations grounded in density functional theory (DFT), which allows for a detailed characterization of the nanosheet and serves as a benchmark for self-consistent charge density functional tight binding (SCC-DFTB) simulations. This approach provides valuable insights into the behavior of the nanosheet under the influence of absorbed OH and H2O molecules by considering calculated parameters for photocatalytic efficiency. Additionally, we extend our investigation to classical molecular dynamics simulations within the ReaxFF framework, modeling the emission of multiple H2O molecules and assessing the subsequent rate of H2 evolution. A key finding of our study reveals that the dissociation of H2O into HO and O molecules significantly enhances both the optical absorbance and conductance of the nanosheet compared to its pristine state. These results underscore the potential of g-C3N4 nanosheets as effective materials for water splitting applications.
Hafnium zirconium oxide (HZO) has attracted significant attention as a core material for next-generation non-volatile memories due to its excellent ferroelectricity in the ultra-thin film regime and its CMOS process compatibility. However, the exploration of its polarization switching mechanism has predominantly relied on static energy barrier analyses, leaving the transient bond formation and cooperative dynamic mechanisms under actual electric field driving unresolved. In this study, we performed Electric-Field-Induced MD simulations on a defect-free ideal HZO lattice using a fine-tuned machine learning force field (MACEField). As a result, we successfully reproduced the P-E hysteresis loop dynamically and demonstrated that the polarization switching in HZO is driven not by conventional simple displacement models (S:N/S:T models), but by the dynamic mutual exchange of 3-coordinated oxygen (O3c) and 4-coordinated oxygen (O4c). Analysis of the oxygen atom displacement trajectories revealed that this pathway is accompanied by a unique "detour" behavior originating from transient cation-oxygen bond formation. Furthermore, we identified an "internal self-compensation mechanism" in which the local volumetric expansion and contraction accompanying the coordination number changes are effectively offset within the cell. These findings provide, from a dynamic perspective, a microscopic physical origin of for HZO's exceptional ability to sustain stable polarization switching without macroscopic strain, a property that has long distinguished HZO from conventional perovskite ferroelectrics yet lacked atomistic explanation. These findings suggest that preserving the integrity of cooperative O3c/O4c exchange pathways, rather than minimizing individual atomic displacements, is the key design principle for endurance and scalability in next-generation ferroelectric memories.
Giant piezoelectricity in ferroelectrics is commonly associated with phase-boundary instabilities, among which the polymorphic phase boundary (PPB) is a prominent example conventionally attributed to the coexistence of ferroelectric phases. Here, using large-scale molecular dynamics simulations of the lead-free (K,Na)NbO3-(Bi,Na)ZrO3 solid solutions, we show that the PPB hosts a percolating multifractal polar domain which governs the dielectric and piezoelectric responses. By quantifying the global fractal dimension and multifractal spectrum width of this polar network, we identify fractal connectivity and multiscale heterogeneity as microstructural order parameters for the PPB. The maximum reversible piezoelectric response occurs when the fractal-domain volume fraction approaches the three-dimensional percolation threshold, suggesting that near-critical polar connectivity enables giant reversible electromechanical coupling. In this mechanism, the fractal backbone preserves polar memory and provides the restoring force required for reversibility, while the surrounding nonfractal regions supply the polar compliance needed for large polarization rotation and strain. These results establish percolating multifractal polar domains as a microscopic mechanism for PPB-enhanced piezoelectricity and suggest fractal connectivity as a design parameter for high-performance piezoelectrics.
Altermagnets combine compensated magnetic order with momentum-dependent spin splitting, offering a route to spintronic functionality without the stray fields of conventional ferromagnets. Mechanical strain provides a promising means of controlling their Néel order, yet the microscopic pathway by which strain reorganizes an altermagnetic texture remains unresolved. Here, we integrate a piezo-driven uniaxial strain cell with scanning nitrogen-vacancy magnetometry to image the magnetic domains of bulk α-MnTe during in situ compression at room temperature. We find that compression reorganizes the magnetic texture through domain coalescence, increasing the size of the largest connected domain while reducing the domain-wall density. Upon unloading, however, the strain-formed domain network does not retrace the loading pathway. Instead, the large connected regions fragment into a new metastable configuration, producing pronounced hysteresis in the maximum domain size and stray-field distribution. These results identify domain connectivity and topology as key carriers of strain-induced magnetic memory. Our work reveals domain coalescence and hysteretic fragmentation as the microscopic pathway of strain control in α-MnTe and establishes a route toward strain-programmable altermagnetic textures and reconfigurable spintronic devices.
Physical Unified Device Architecture (PUDA) is an AI-native hardware harness for self-driving laboratories (SDLs). Rather than building a human-centered graphical user interface (GUI) orchestration layer, PUDA creates a command-line runtime environment that lets agents observe, orient, decide, and act over experiments while hardware execution remains deterministic, atomic, and auditable. Headless by design, devices appear through discoverable command-line interfaces, JSON protocols are routed through a distributed messaging system, and command responses, data products, and reports are preserved as structured records. PUDA organizes protocols, runs, samples, measurements, and command logs into an AI-native data structure linked by run identifiers and timestamps, preserving provenance from submitted protocol through hardware response to resulting data products. PUDA separates scientific orchestration from physical operation and data telemetry: agents choose experiments, while PUDA executes validated commands and captures provenance-linked state, responses, and data. The contribution is not another optimizer, orchestrator, or recipe language. It is a practical execution and data environment for agentic SDLs; the broader physical AI implication is that PUDA provides an AI-native hardware harness for AI systems to interact with physical tools.
Accurate determination of Hubbard interaction parameters is essential for beyond-DFT approaches such as DFT+$U$, DFT+DMFT, and DFT+$U$+$V$ in correlated materials. In practice, however, these parameters are often chosen empirically, limiting their transferability across materials. Advanced computational approaches such as the constrained random-phase approximation (cRPA) provide a rigorous route for evaluating Hubbard interactions, but their computational cost remains a bottleneck for large-scale materials screening. Here, we present machine-learning (ML) models for predicting cRPA-derived Hubbard interaction parameters: effective on-site $U_{\rm eff}$, inter-site $V$, and Hund's coupling $J$ for transition-metal oxides (TMOs). We combine ensemble-learning models with a regression-based brute-force search (BFS) approach to achieve both predictive accuracy and explicit analytical expressions. We construct features that capture electronic, structural, and atomic properties, including the TM-$d$ bandwidth and TM-$d$/O-$p$ band-center separation, as physically motivated descriptors of localization and screening. Our ensemble models achieve RMSEs of 0.148 eV, 0.062 eV, and 0.007 eV for $U_{\rm eff}$, $V$, and $J$, respectively. The derived analytical forms directly relate $U_{\rm eff}$ to electron localization and TM-$d$/O-$p$ hybridization, suggest the importance of hybridization and structural compactness in determining $V$, and indicate that $J$ is governed primarily by elemental descriptors of the TM ion. Together, the present study provides an efficient approach for predicting cRPA-derived $U_{\rm eff}$, $V$, and $J$, while offering physical insight into the factors underlying these Hubbard interactions.
Machine learning (ML) is transforming materials discovery by enabling rapid prediction of properties that previously required computationally expensive first-principles calculations. Yet most current ML models remain fundamentally limited to zero-temperature descriptions, learning static lattice energies while neglecting the thermodynamic effects that govern materials behaviour at finite temperature. Because phase stability, functional response, and performance are governed by free-energy landscapes rather than static energies alone, this limitation represents a major barrier to predictive materials design under realistic operating conditions. In this Perspective, we argue that developing thermodynamics-informed ML constitutes one of the most important and least explored frontiers in materials discovery. We examine the fundamental shortcomings of energy-based models, highlighting the essential roles of entropy and anharmonicity in determining free energies and materials functionality. We review emerging strategies, including machine-learned interatomic potentials and hybrid ML-statistical mechanics frameworks, while identifying key challenges related to data availability, transferability, and thermodynamic consistency. Building on these advances, we outline a roadmap for thermodynamics-informed ML centred on direct free-energy learning, entropy-aware representations, and adaptive sampling across temperature. We highlight the transformative opportunities this paradigm offers for energy materials and argue that the next generation of ML models must move beyond static energy predictions towards a thermodynamic description of materials behaviour under realistic operating conditions.
Fracture and calving of glaciers are key contributors to ice-mass loss and sea-level rise, yet predictive modeling remains challenging. Fracture in grounded glaciers is driven by gravitational forces and is typically studied within the framework of quasi-static linear elastic fracture mechanics. In this work, we show that purely quasistatic brittle fracture simulations within the phase field fracture framework under fixed self-weight can become strongly overdriven after crevasse initiation, producing unphysical thickening of the diffusive crack band and diffuse damage patterns. This pathology arises because gravity drives a growing region ahead of the crack tip beyond the strength surface. To resolve this, we show that the post-nucleation propagation is fundamentally a dynamic instability rather than a quasistatic process and propose the use of dynamic formulations of fracture. We demonstrate that accounting for inertia results in sharp, localized cracks that propagate through the ice thickness. As a second objective, this paper introduces a new dynamic formulation of the phase-field fracture model of Kumar et al. (J. Mech. Phys. Solids 2018) in which elastic, inertial, and gravitational contributions are degraded consistently in fractured regions.
This work is devoted to the mechanisms of dielectric response and electric conductivity of suspensions consisting of the nematic liquid crystal 5CB with different concentrations (from 0 to 10 wt.%) of ferroelectric BaTiO_3 nanoparticles with an average size of 24 nm. We revealed that the incorporation of nanoparticles influences significantly the dielectric permittivity magnitude and anisotropy, as well as dielectric losses of the suspension. A pronounced temperature dependence of the anisotropic dielectric permittivity of the suspensions was found at lower temperatures corresponding to the mesophase state; but it is also present at higher temperatures corresponding to the isophase. The dependence of the mesophase-isophase transition temperature on the concentration of BaTiO_3 nanoparticles appeared nonmonotonic. With increasing temperature, both the capacitance and the electrical resistance of the pure liquid crystal increase, as well as it increases in the suspensions with small concentration of BaTiO_3 nanoparticles. Due to space charge accumulation in the shells of nanoparticles, larger concentrations of BaTiO_3 nanoparticles influence strongly the ionic transport by promoting the formation of ionic-electronic screening. This effect modifies the dielectric properties and conduction mechanisms of the suspension, leading to the nonmonotonic dependence of the mesophase - isophase transition temperature versus the nanoparticle concentration.
Crystalline symmetries stabilize topological states with distinct electronic properties, while altermagnets exhibit momentum-dependent spin splitting without net magnetization. Here, we combine first-principles calculations with a minimal tight-binding model to realize $C$-paired spin-valley-locked nodal-line fermions in the $d$-wave altermagnet $\mathrm{AV_2X_2O}$ (A = Rb, Cs, or K; X = Te, Se, or S). The low-energy electronic structure hosts coexisting spin-degenerate and spin-polarized nodal lines around $C_{4z}$-paired valleys near the Fermi level. The spin-polarized nodal lines are protected by the out-of-plane mirror symmetry $\mathcal{M}_z$ and remain robust against spin-orbit coupling. The minimal model reveals their microscopic origin and establishes a general design principle for their isolation. Layer engineering and electronic correlations serve as material-specific knobs for realizing these isolated spin-valley-locked nodal lines near the Fermi level. Our results establish the $\mathrm{AV_2X_2O}$ family as a versatile platform for exploring topological spin-valley locking in $d$-wave altermagnets.
We investigate magnetic breakdown near Van Hove singularities (VHSs) in the electron-doped rhombohedral tetralayer graphene, where chiral superconductivity has recently been reported. Using the noninteracting band structure and Kubo formula, we identify anomalous Shubnikov-de Haas effects: Ring-like structures in the Landau fan and anomalous high-frequency peaks in the frequency spectra. These anomalous quantum oscillations can be understood by the reconstruction from magnetic breakdown among three nearby Fermi pockets separated by VHSs. Remarkably, these qualitative anomalous features persist into a stronger-VHS regime, where the semiclassical picture breaks down, and chiral superconductivity emerges. The temperature and (weak) disorder dependence of the oscillations are also investigated. Our results establish that the magnetic-breakdown-induced anomalous quantum oscillation provides a general distinctive probe for the underlying Fermi-surface geometry associated with VHSs and may explain the recent quantum oscillation experiment in rhombohedral tetralayer graphene [arXiv:2606.05356].
In this work, we investigate the problem of electromagnetic angular momentum within the framework of topological materials, described by axion electrodynamics. Specifically, we calculate the electromagnetic angular momentum for a spherical sample of magnetic Weyl semimetal in the presence of a point charge. We then analyze how, by moving the point charge quasi-statically, the angular momentum stored in the electromagnetic field can be converted into mechanical angular momentum of the sphere. Finally, we derive analytical expressions for the resulting angular velocity and angular displacement, demonstrating that the effect is enhanced in the Weyl semimetal compared with a topological insulator.
Controlling magnetic losses requires identifying the microscopic origin of extrinsic linewidth broadening. We introduce a Lippmann--Schwinger framework for two-magnon scattering from weak defect potentials, showing that scattering is set by the overlap between the perturbation of the effective field and degenerate spin-wave states. Applied to iron garnet films, YIG/GGG(111) and YIG/GGG(110), we identify crystallographic directions along which these perturbation are elongated. This approach provides a link between the crystallography, effective field on the sample scale and two-magnon scattering.
Deployable shape-morphing structures that transform from flat sheets into stable three-dimensional configurations are highly desirable for applications ranging from soft robotics and biomedical devices to adaptive architecture and aerospace systems. Existing kirigami-based morphing systems primarily rely on isotropic deployment, compliant soft materials, or external constraints to maintain deployed shapes, which limits geometric programmability, structural integrity, and applicability in rigid-material systems. Here, we present an inverse design framework for anisotropic bistable kirigami structures that enables programmable shape morphing through controlled geometric frustration and instability-induced deployment. The framework combines a semi-analytical mechanical model with geometry to establish a direct connection between geometric transformation and the underlying energy landscape. We show that instability-induced shape morphing leads to tunable bistability and directional deployment in anisotropic kirigami structures. The results are validated through finite element simulations and experiments, demonstrating stable deployed configurations and programmable anisotropic morphing. The proposed framework further provides a general design strategy that can be integrated with various active actuation systems, enabling broader engineering applications.