This paper introduces a formal end-to-end hybrid classical-AI quantum-computing architecture designed to model, diagnose, and predict the behavior of pressure-tuned moiré quantum materials. By partitioning the problem into a classical AI surrogate stage (handling structural relaxations and phonon dynamics) and a quantum variational stage (solving the low-energy interacting Hamiltonian), we overcome the exponential scaling limits of classical exact diagonalization (ED) while ensuring structural precision over thousands of atoms.
2. Mathematical Framework
-------------------------------------------------------------
| 1. Structural Relaxation (Classical AI) |
| Minimize H_enthalpy via MLFF (DPmoire) under P and s |
-------------------------------------------------------------
|
v
-------------------------------------------------------------
| 2. Multi-Physics Property Extraction |
| - Pass coordinates to PARPHOM for Phonon Dynamics [D(k)] |
| - Fit Interlayer Distances D_z(P) and Hoppings w(P) to DFT |
-------------------------------------------------------------
|
v
-------------------------------------------------------------
| 3. Continuum Electronic Modeling |
| Construct K-Valley H(θ,r,P); compute Chern Numbers & Bands |
-------------------------------------------------------------
|
v
-------------------------------------------------------------
| 4. Quantum Correlated State Engine |
| Project Interactions onto Flat Bands -> H_eff(P) |
| Execute Hardware-Efficient HVA Circuit (PennyLane/Qiskit) |
-------------------------------------------------------------
2.1 Structural Relaxation under Twist, Sliding, and Pressure
To determine the equilibrium atomic configurations of a twisted multilayer assembly under macroscopic twist angles $\theta$, rigid interlayer sliding vectors $\mathbf{s}$, and external pressure
$P$, we define a physics-informed MLFF enthalpy objective function (Liu et al., 2025):
$$H(\{\mathbf{r}_i\}, \theta, \mathbf{s}, P) = E_{\text{MLFF}}(\{\mathbf{r}_i\}, \theta, \mathbf{s}) P \cdot V(\{\mathbf{r}_i\})$$
where $\{\mathbf{r}_i\}$ denotes the full set of relaxed $3\text{D}$ atomic coordinates within the moiré supercell, and
$V(\{\mathbf{r}_i\})$ represents the instantaneous volume of the simulation cell. The machine-learning force field energy term
$E_{\text{MLFF}}$ is evaluated using a deep neural network potential trained on ab initio molecular dynamics (AIMD) datasets generated under varying out-of-plane strains, satisfying:
$$-\nabla_i H(\{\mathbf{r}_i\}, \theta, \mathbf{s}, P) = \mathbf{0}$$
2.2 Pressure-Dependent Effective K-Valley Hamiltonian, Following structural relaxation, the low-energy electronic structure at the $\pm K$ valleys is projected onto a continuum Hamiltonian. Extending the standard Bistritzer-MacDonald framework to encompass pressure dependence (Anfa et al., 2025), the single-particle effective Hamiltonian for a twisted bilayer system is expressed as:
$$\mathcal{H}(\theta, \mathbf{r}, P) = \begin{pmatrix} -\frac{\hbar^2}{2m^*(P)}(\mathbf{k}-\boldsymbol{\kappa}_-)^2 \Delta_t(\mathbf{r},P) & \Delta_T(\mathbf{r},P) \\ \Delta_T^\dagger(\mathbf{r},P) & -\frac{\hbar^2}{2m^*(P)}(\mathbf{k}-\boldsymbol{\kappa}_ )^2 \Delta_b(\mathbf{r},P) \end{pmatrix}$$; where $\boldsymbol{\kappa}_\pm$ denote the shifted moiré Dirac points, and
$m^*(P)$ represents the pressure-dependent effective mass. The interlayer tunneling matrix $\Delta_T(\mathbf{r},P)$ and intra-layer moiré potentials $\Delta_{t/b}(\mathbf{r},P)$-