---
title: Hybrid Quantum-Classical PINNs
url: https://www.emergentmind.com/topics/hybrid-quantum-classical-physics-informed-neural-network-hqpinn
type: topic
---

# Hybrid Quantum-Classical PINNs

Searching arXiv for recent HQPINN and related QPINN papers to ground the article.
Hybrid Quantum-Classical Physics-Informed Neural Network (HQPINN) denotes a class of physics-informed learning architectures in which a physics-informed neural network is coupled to a quantum component, typically a parameterized quantum circuit or a quantum feature map, while residual construction, loss assembly, and optimization remain partly or predominantly classical. Across the literature, the term covers several distinct patterns: fixed quantum feature augmentation for classical PINNs, late-stage quantum latent processors embedded inside otherwise classical PINNs, quantum surrogate models trained through classical physics-informed losses, and trainable-embedding quantum-assisted PINNs for time-dependent PDEs. What unifies these variants is that the learned approximation is constrained by governing equations rather than labeled solution data alone, and that a quantum subroutine contributes either representational structure, latent transformation, or trainable nonlinear features within a hybrid differentiable workflow [2509.16247][2304.11247][2602.14596].

## 1. Conceptual scope and defining characteristics

HQPINNs are most accurately understood as a family rather than a single architecture. In the minimal form, a classical PINN is augmented with quantum-circuit-generated features and trained using an ordinary residual loss. A proof-of-concept example solves the initial value problem
\[
\frac{dy}{dx}+2y=0,\qquad y(0)=1,
\]
by sending each scalar collocation point \(x\in[0,1]\) through a two-qubit circuit, measuring computational-basis probabilities, concatenating the resulting four-dimensional quantum feature vector with the original input, and learning the solution with a classical feed-forward network under a hard initial-condition ansatz [2509.16247]. In that setting, the quantum part is fixed feature extraction rather than a trainable variational layer.

A stronger notion of hybridization inserts a trainable quantum circuit into the interior of the neural approximator. In a 3D computational fluid dynamics application, a classical multilayer perceptron maps spatial coordinates \((x,y,z)\) to a 16-dimensional latent representation, then a parallel hybrid head routes 15 units to a variational quantum circuit and one unit to a classical branch before a final dense layer predicts \((u,v,w,p)\) for the steady incompressible Navier–Stokes equations in a parameterized Y-shaped mixer [2304.11247]. In this design, the PINN remains classical in its loss construction and PDE differentiation, while the latent representation is partly processed quantumly.

Another line replaces the usual classical surrogate altogether with a quantum model while retaining a classical physics-informed training loop. Continuous-variable QPINNs for Poisson and heat equations use qumode circuits as the trainable function approximator and compute PDE residuals classically through autodiff or consistency-based derivative constructions [2209.14754][2503.12244]. A plausible implication is that the boundary between “HQPINN” and “QPINN” is architectural rather than categorical: some works use the quantum model as one module inside a classical network, whereas others use the quantum model as the surrogate and preserve hybridization at the training-loop level.

More recent work frames HQPINNs explicitly as hybrid classical-to-quantum-to-classical solvers. For nonlinear PDEs such as Burgers, Allen–Cahn, and Korteweg–de Vries, a classical backbone first extracts latent features from \((x,t)\), a parameterized quantum circuit transforms a bounded latent vector through angle encoding and entangling layers, and a final classical dense layer maps quantum measurements to the scalar field \(u(x,t)\) [2606.04679]. Closely related trainable-embedding quantum-assisted PINNs for parabolic and reaction-diffusion PDEs study whether the embedding from coordinates to quantum angles should be classical or quantum, identifying the embedding stage itself as a primary design variable rather than mere preprocessing [2602.14596][2602.09291].

## 2. Architectural patterns

The literature supports several recurring HQPINN motifs.

The first is the **quantum feature map plus classical PINN** pattern. In the ODE proof of concept, the trial solution is
\[
\hat y(x)=1+x\,g_\theta([x,q(x)]),
\]
where \(g_\theta\) is a classical MLP and \(q(x)\) is a quantum feature vector obtained from a two-qubit circuit
\[
U(t)=\big(Rx(-2t)\otimes Rx(-2t)\big)(H\otimes I),\qquad t=x.
\]
The measurement output is
\[
q(x)=p_Z(x)=\big[P_Z(00),P_Z(01),P_Z(10),P_Z(11)\big],
\]
and the fused input is \(z(x)=[x,q(x)]\in\mathbb R^5\) [2509.16247]. The initial condition is enforced exactly by construction.

The second is the **late latent-space quantum processor**. In the Y-mixer CFD architecture, the classical front-end has 5 hidden layers of 64 neurons, followed by a 16-dimensional penultimate layer. The hybrid head splits those 16 units into 15 quantum inputs and 1 classical scalar. The quantum branch uses repeated encoding blocks with depth \(d=5\), followed by variational layers repeated \(m=2\), and returns expectation values such as
\[
\langle \sigma_z\rangle=\bra{\psi}\sigma_z\ket{\psi}
\]
that are fused with the classical branch before a final dense map to \((u,v,w,p)\) [2304.11247]. This architecture treats the PQC as a latent feature transformer rather than as the entire solution network.

The third is the **classical–quantum–classical serial architecture**. Work on explaining the advantage of quantum-enhanced PINNs embeds a three-qubit variational circuit between a classical encoding network and a classical decoding network. The encoder maps \((t,x)\) into features \(i_j\), repeated data re-uploading encodes them through alternating \(R_Y\) and \(R_X\) rotations before each variational layer, expectation values of Pauli-\(Z\) are measured on each qubit, and a classical decoder maps the resulting vector to the scalar solution \(u(t,x)\) [2601.15046].

The fourth is the **classical backbone with quantum hidden layer**. For 2D time-dependent Maxwell equations in TE\(_z\) form, a classical PINN backbone performs periodic coordinate preprocessing, learned periodic time mapping, random Fourier feature embedding, and several classical hidden layers. The second-last classical layer is then replaced by a seven-qubit PQC; Pauli-\(Z\) expectation values from each qubit are fed into a final classical layer predicting \((E_z,H_x,H_y)\) [2506.23246]. This model is explicitly hybrid and fully differentiable in simulation.

The fifth is the **trainable-embedding quantum-assisted PINN**. In FNN-TE-QPINN, a classical feed-forward network generates trainable angles
\[
\Gamma(\mathbf x,t)=[\phi_1(\mathbf x,t),\ldots,\phi_n(\mathbf x,t)]^\top
\]
for the encoding unitary
\[
U_{\mathrm{enc}}(\mathbf x,t;\theta_{\mathrm{emb}})=\prod_{k=1}^{n}R_y(\phi_k(\mathbf x,t)),
\]
after which a hardware-efficient variational circuit acts and the final scalar output is read out as
\[
\tilde u(\mathbf x,t)=\langle \psi_{\mathrm{var}}(\mathbf x,t)|\mathcal O|\psi_{\mathrm{var}}(\mathbf x,t)\rangle,\qquad
\mathcal O=\bigotimes_{i=1}^{n}\sigma_z^{(i)}.
\]
In QNN-TE-QPINN, the embedding stage is itself quantum [2602.14596]. An analogous extended framework for reaction-diffusion systems supports both classical and quantum embeddings while keeping the variational ansatz and loss fixed [2602.09291].

The sixth is the **additive–multiplicative coupling architecture**. A purely theoretical proposal, QPINN-MAC, couples a classical MLP output \(\widehat y_j(t,\vec{\mathcal W})\) with a QNode expectation value \(\braket{\hat O}_{\vec\Theta_j}\) through
\[
y_{MAC_j}(\hat O,t,\vec{\mathcal W},\vec\Theta_j)=
\widehat y_j(t,\vec{\mathcal W})\braket{\hat O}_{\vec\Theta_j}+\braket{\hat O}_{\vec\Theta_j}.
\]
This retains explicit classical and quantum branches while making the quantum term both additive and multiplicative [2511.07216].

## 3. Physics-informed objectives and constraint enforcement

Across HQPINNs, the physics-informed component is structurally orthodox: PDE or ODE residuals are computed at collocation points, and the total objective combines residual, boundary, and initial-condition terms when needed.

For the minimal ODE HQPINN, the residual is
\[
r(x_i;\theta)=\frac{d\hat y}{dx}(x_i)+2\hat y(x_i),
\]
and the loss is the mean-squared residual
\[
\mathcal L_{\mathrm{ODE}}(\theta)=\frac{1}{N}\sum_{i=1}^{N}r(x_i;\theta)^2.
\]
There is no supervised data term and no initial-condition penalty because the architecture satisfies \(\hat y(0)=1\) exactly [2509.16247].

For steady incompressible CFD in the Y-shaped mixer, the governing equations are
\[
-(\mathbf v\cdot\nabla)\mathbf v+\nu\Delta\mathbf v-\frac{1}{\rho}\nabla p=0,\qquad
\nabla\cdot\mathbf v=0,
\]
which are enforced through separate residuals for the three momentum equations and continuity, summed into
\[
\mathcal L_{\mathrm{PDE}}=\langle r_u^2\rangle_{\Omega_f}+\langle r_v^2\rangle_{\Omega_f}+\langle r_w^2\rangle_{\Omega_f}+\langle r_c^2\rangle_{\Omega_f},
\]
plus wall, inlet, and outlet boundary losses [2304.11247].

For parabolic PDEs, the generic residual formulation is
\[
\mathcal R_{\mathrm{PDE}}(\mathbf x,t):=\mathcal D(u(\mathbf x,t);\beta)-q(\mathbf x,t)=0,
\]
with corresponding boundary and initial residuals
\[
\mathcal R_{\mathrm{BC}}(\mathbf x,t):=\mathcal B(u(\mathbf x,t))-b(\mathbf x,t)=0,\qquad
\mathcal R_{\mathrm{IC}}(\mathbf x):=u(\mathbf x,0)-u_0(\mathbf x)=0.
\]
The trainable-embedding QPINNs optimize
\[
\mathcal L(\theta_{\mathrm{var}},\theta_{\mathrm{emb}})
=
\mathcal L_{\mathrm{PDE}}
+\lambda_{\mathrm{BC}}\mathcal L_{\mathrm{BC}}
+\lambda_{\mathrm{IC}}\mathcal L_{\mathrm{IC}},
\]
where the surrogate itself is a quantum expectation value [2602.14596].

For nonlinear PDE benchmarks in Burgers, Allen–Cahn, and KdV, the HQPINN uses the unweighted sum
\[
L_T=L_{IC}+L_{BC}+L_{PDE},
\]
with derivatives such as \(u_t\), \(u_x\), \(u_{xx}\), and \(u_{xxx}\) obtained by autodiff through the full hybrid model [2606.04679].

In Maxwell QPINNs, the physics loss combines residual, initial, symmetry, and energy-conservation terms:
\[
\mathcal L_{\mathrm{tot}}=\mathcal L_{\mathrm{res}}+10\,\mathcal L_{\mathrm{IC}}+10\,\mathcal L_{\mathrm{sym}}+10\,\mathcal L_{\mathrm{energy}}.
\]
The added energy term is derived from Poynting’s theorem and is used to stabilize training in the source-free vacuum case [2506.23246].

Not all HQPINNs are PDE solvers in the narrow numerical-analysis sense. Seismic inversion HQ-PINNs use a physics-based decoder rather than a differential operator residual: a hybrid quantum encoder predicts elastic parameters, a geophysical forward model computes reflectivity via Aki–Richards or normal-incidence relations, and the network is trained by seismic reconstruction misfit plus low-frequency regularization [2503.05009]. Hydrological HQC-PINNs similarly combine classification loss with differentiable PDE-constrained hydrological regularizers rather than solving a PDE field directly. Their objective is
\[
\mathcal L_{\mathrm{HQC}}
=
\mathcal L_{\mathrm{data}}
+\lambda_{\mathrm{SV}}\mathcal L_{\mathrm{SV}}
+\lambda_{\mathrm{M}}\mathcal L_{\mathrm{Manning}},
\]
where \(\mathcal L_{\mathrm{SV}}\) is a Saint-Venant continuity residual and \(\mathcal L_{\mathrm{Manning}}\) enforces Manning consistency [2604.09374].

## 4. Quantum representations, embeddings, and differentiability

A central distinction across HQPINNs is whether the quantum component is fixed, trainable, or merely quantum-inspired.

The simplest case is fixed quantum feature extraction. In the ODE proof of concept, the circuit is not variational, has no trainable quantum parameters, and gradients do not pass through the quantum branch. Quantum outputs are generated first, standardized as a feature matrix \(Q\), and then treated as fixed classical inputs [2509.16247]. This is hybrid in the sense of “quantum feature map plus classical PINN,” not in the sense of joint quantum-classical optimization.

Trainable variational circuits appear in most later HQPINNs. The Maxwell architecture uses a seven-qubit PQC with 4 variational layers, angle-embedded latent activations, and several candidate ansätze, including Basic Entangling Layers, Strongly Entangling Layers, Cross-Mesh variants, and a No Entanglement control [2506.23246]. The nonlinear-PDE HQPINN uses a hardware-efficient ansatz in which a bounded classical latent vector
\[
z=(z_1,\dots,z_{n_q})\in[-1,1]^{n_q}
\]
is angle-encoded through
\[
R_x(z_i)=\exp\!\left(-\frac{i z_i\sigma_x}{2}\right),
\]
followed by \(m\) layers of trainable single-qubit \(R_x(\theta_{\ell,i})\) rotations and a ring of CNOT gates [2606.04679].

Trainable embeddings have become a major research direction. FNN-TE-QPINN generates encoding angles with a small classical network, whereas QNN-TE-QPINN uses an auxiliary quantum model to generate the embedding itself [2602.14596]. For multi-species reaction–diffusion systems, the trainable embedding is
\[
\Gamma(x,t)=\big[\alpha_1(\tilde x,\tilde t),\ldots,\alpha_{N_q}(\tilde x,\tilde t)\big]^\mathsf T,
\]
and the encoding operator is
\[
U_{\mathrm{enc}}(x,t;\theta_{\mathrm{emb}})
=
\bigotimes_{i=1}^{N_q}
R_y\!\big(\alpha_i(\tilde x,\tilde t;\theta_{\mathrm{emb}})\big).
\]
The resulting state is then processed by a shared variational circuit, and multiple PDE fields are read out through species-specific observables
\[
\mathcal O_i=\sum_{j\in\mathcal Q_i}\sigma_z^{(j)}.
\]
This suggests that in multi-field HQPINNs the measurement operator itself can encode task structure, not just the circuit architecture [2602.09291].

Continuous-variable variants use qumodes rather than qubits. In these models, inputs are encoded by displacement gates, Gaussian gates perform affine-like transformations, and non-Gaussian Kerr gates provide the analogue of activation functions. The output is typically a quadrature expectation value such as
\[
u_\theta(x)=\langle \psi_\theta(x)|\hat X|\psi_\theta(x)\rangle
\]
or multi-output variants for \(u\) and its derivatives [2209.14754][2503.12244][2404.15015].

Differentiability strategy is a major practical separator. Some architectures rely on classical autodiff through a simulator, including the Poisson and heat QPINNs implemented with Strawberry Fields and TensorFlow [2209.14754][2503.12244], the CFD HQPINN using adjoint differentiation on a classical simulator [2304.11247], and the Maxwell QPINN using a custom GPU-accelerated differentiable simulator [2506.23246]. Other works explicitly use parameter-shift for quantum parameters while using classical backpropagation or hybrid autodiff for embedding and post-processing, as in the trainable-embedding parabolic and reaction–diffusion frameworks [2602.14596][2602.09291] and the hydrological HQC-PINN [2604.09374].

A nearby but distinct line uses quantum-inspired rather than genuinely quantum computation. Battery SOH prognosis with a “QPINN” uses a fixed quantum feature mapping defined through an \(n\)-qubit kernel
\[
K(x_i,x_j)=|\langle\psi(x_i)|\psi(x_j)\rangle|^2,
\]
approximated by a Nyström embedding
\[
\psi_q(x)=K(x,S)\,K(S,S)^{-1/2}.
\]
This feature map is detached and consumed by a classical physics-informed model; there is no explicit variational circuit, no quantum hardware backend, and no trainable quantum parameters [2604.10362]. By strict usage, this is adjacent to HQPINNs rather than a canonical instance.

## 5. Empirical regimes, reported performance, and interpretation

Published HQPINN results do not support a single blanket claim of superiority. Instead, they suggest regime dependence.

For very simple ODEs, feasibility is established but advantage is not. The two-qubit ODE model “successfully reproduced the analytical solution \(y(x)=e^{-2x}\)” with only minor deviations at larger \(x\), but no classical baseline, no numerical error table, and no scaling study are reported [2509.16247]. This is best interpreted as proof of concept.

For CFD in complex shapes, the 3D Y-mixer HQPINN reports a **21%** improvement over a matched classical mini-batch model in terms of loss value or “21% higher accuracy,” depending on phrasing, while preserving the standard Navier–Stokes PINN formulation [2304.11247]. However, the hybrid model still exhibits asymmetry between physically symmetric pipes and is trained only on a simulator, so the result indicates potential rather than CFD-grade validation.

For trainable-embedding parabolic PDEs, the evidence is more granular. On the 1D heat equation, after 150 epochs the reported losses are
\[
4.8670\times 10^{-6}\quad\text{(FNN-TE-QPINN)},
\]
\[
8.9515\times 10^{-6}\quad\text{(PINN)},
\]
\[
1.3406\times 10^{-4}\quad\text{(QNN-TE-QPINN)}.
\]
The best maximum absolute error is
\[
4.5431\times 10^{-7}\quad\text{(FNN-TE-QPINN)},
\]
and the 1D inference metrics favor the classical-embedding hybrid over both the classical PINN and the fully quantum embedding [2602.14596]. In 2D heat, FNN-TE-QPINN again dominates training loss and error, while QNN-TE-QPINN is far more parameter-efficient but notably less accurate:
\[
271 \text{ parameters (PINN)},\qquad
8114 \text{ (FNN-TE-QPINN)},\qquad
96 \text{ (QNN-TE-QPINN)}.
\]
A plausible implication is that embedding design governs trainability more strongly than the mere presence of a quantum variational block.

The multi-species reaction–diffusion study reaches a similar conclusion. FNN-TE-QPINN generally achieves the lowest final loss, whereas QNN-TE-QPINN often converges faster and more smoothly but saturates earlier at a higher error floor [2602.09291]. This reinforces the view that fully quantum embeddings can be compact and stable, but classical trainable embeddings remain stronger in final approximation quality under current simulation constraints.

For nonlinear PDEs, the case for HQPINNs is more substantial. In Burgers’ equation, relative \(L^2\) error decreases from
\[
8.46\times 10^{-2}\quad\text{(PINN)}
\]
to
\[
2.12\times 10^{-2}\quad\text{(HQPINN)},
\]
about a fourfold reduction. In Allen–Cahn, the reduction is from
\[
7.04\times 10^{-2}
\]
to
\[
1.35\times 10^{-2},
\]
about fivefold. In KdV, the gain is more modest, from
\[
2.85\times 10^{-2}
\]
to
\[
2.07\times 10^{-2}.
\]
These results are accompanied by smoother training dynamics and reduced loss oscillations, with the most consistent benefit appearing when the PDE exhibits stiffness, sharp gradients, localized interfaces, or multiscale structure [2606.04679]. The same study reports that output-stage placement of the PQC is consistently best, and that optimal qubit count and depth are problem dependent rather than monotone with model size.

Maxwell QPINNs show that hybridization can improve both accuracy and parameter efficiency on a coupled unsteady PDE. The best settings yield up to **19%** lower \(L^2\) error than the classical PINN while using about **19% fewer trainable parameters**, with the quantum model inserted as a hidden layer rather than replacing the full network [2506.23246]. However, this advantage depends strongly on input-angle scaling, circuit ansatz, and the inclusion of a global energy term in the vacuum case.

Evidence from scientific inverse problems is mixed but informative. The seismic HQ-PINN estimates impedances accurately on synthetic and field data, with ansatz choice \(R_x\), \(R_y\), or \(R_z\) making little difference and even a “No Ansatz” configuration remaining effective [2503.05009]. By contrast, the high-speed-flow HQPINN shows clear benefits only for harmonic or smooth regimes. On a smooth 1D Euler traveling wave, the quantum and hybrid architectures are strongly parameter-efficient and accurate; on a moving contact discontinuity and on 2D transonic aerofoil flow, classical PINNs remain more reliable, and under-parameterized quantum branches degrade trainability [2503.02202]. This suggests that the Fourier-like inductive bias of small PQCs is advantageous for harmonic structure but poorly matched to discontinuities.

A separate hydrological HQC-PINN reports convergence in roughly **3x** fewer epochs and about **44%** fewer trainable parameters than a classical PINN while maintaining competitive classification accuracy. The best \(8\)-qubit, \(3\)-layer model reaches \(71.8\%\) accuracy with \(18{,}944\) parameters versus \(69.4\%\) and \(33{,}793\) parameters for the classical PINN, and a transfer-learning variant reaches \(73.6\%\) [2604.09374]. Because this model is a physics-constrained classifier rather than a direct PDE field solver, it broadens the operational meaning of HQPINN beyond forward PDE approximation.

## 6. Theory, trainability, and open issues

Theoretical work has begun to formalize hybridization, but the literature remains uneven between rigorous theory and empirical validation.

QPINN-MAC provides one of the few explicit universal-approximation and trainability theorems for a hybrid ansatz. Its defining componentwise form,
\[
y_{MAC_j}(t)=\widehat y_j(t,\vec{\mathcal W})\,\braket{\hat O}_{\vec\Theta_j}+\braket{\hat O}_{\vec\Theta_j},
\]
is proved dense in \(L_{\mathbb R^M}^p(\mathcal K)\) under the stated assumptions, and the paper derives an algebraic gradient bound
\[
\|\nabla_{\vec\Xi^\ast}\mathfrak L\|_{p,M}\le
\mathcal O\!\left(\frac{1}{\sqrt{\mathcal N N}}\right)
\]
as a mitigation against barren-plateau-style gradient collapse [2511.07216]. This establishes a theoretical pathway for coupling classical and quantum branches so that the classical factor modulates the quantum gradient path. Yet the paper is theory-only: it contains no experiments, no optimizer comparisons, and no PDE benchmarks.

A complementary explanatory study argues that the practical advantage of quantum-enhanced PINNs lies primarily in faster convergence rather than markedly better asymptotic accuracy. On a family of 1D diffusion–convection PDEs, serial classical–quantum–classical qPINNs often require one to two orders of magnitude fewer epochs than parameter-matched classical PINNs to reach a given MSE, particularly when the attainable error floor is low and collocation data are sparse [2601.15046]. This suggests that optimization efficiency rather than raw expressivity may be the principal near-term value proposition of HQPINNs.

Trainability pathologies remain a recurrent theme. Continuous-variable QPINNs for Poisson report that plain SGD outperforms Adam, Nadam, RMSprop, Adadelta, SPSA, and SGD+L-BFGS-B, reversing common classical PINN heuristics and highlighting that optimizer behavior does not transfer cleanly from classical to quantum settings [2209.14754]. The Maxwell work identifies a distinct “black hole” loss-landscape phenomenon, in which the quantum model initially learns useful dynamics but later collapses to an almost-zero solution for \(t>0\); adding the energy conservation term removes this failure mode in the vacuum case [2506.23246]. A plausible implication is that global physical invariants may serve both as regularizers and as trainability stabilizers in HQPINNs.

Several misconceptions are addressed by the current literature. First, hybridization does not imply a trainable quantum model; fixed quantum feature maps already satisfy a weak notion of HQPINN [2509.16247]. Second, more quantum content is not uniformly better. In trainable-embedding studies, the most quantum-centric architectures are more compact but usually less accurate than hybrids with classical embeddings [2602.14596][2602.09291]. Third, quantum circuits are not automatically beneficial on hard PDEs. High-speed-flow experiments show that Fourier-like quantum inductive bias can be actively harmful near discontinuities or shocks [2503.02202]. Fourth, most reported gains are simulator-only and do not establish hardware-era advantage. Nearly all cited studies use noiseless or classical-simulator settings; shot noise, hardware noise, calibration drift, and finite-shot gradient estimation remain largely unexplored [2304.11247][2606.04679][2506.23246][2604.09374].

The field also has persistent reproducibility limitations. Common omissions include explicit qubit counts, gate-level circuit definitions, optimizer hyperparameters, batch sizes, collocation-point counts, loss-weight schedules, software versions, hardware environments, and full baseline protocols. This is particularly pronounced in early proof-of-concept papers [2509.16247][2304.11247][2404.15015]. As a result, many current HQPINN results should be interpreted as architectural evidence rather than mature numerical benchmarks.

Taken together, the literature indicates that HQPINNs are most credible when the quantum component is tightly co-designed with the classical backbone, the encoding interface, and the physics-informed objective. The strongest empirical support currently favors hybrid latent processing and trainable classical embeddings for smooth, stiff, or multiscale PDEs, while discontinuous and shock-dominated regimes remain challenging and often still favor classical PINNs [2606.04679][2503.02202]. Future work will likely determine whether these gains persist under realistic hardware constraints, whether quantum modules can outperform simply enlarging the classical backbone, and whether physics-based global constraints can systematically regularize hybrid optimization landscapes [2506.23246][2511.07216].

Source: https://www.emergentmind.com/topics/hybrid-quantum-classical-physics-informed-neural-network-hqpinn