Skip to content
DraconiticResearch
Search Sources API
203,234 studies31 sourcesUpdated 49 min ago
Peer-reviewedCC BY 3.0RSC Advances

Hybrid quantum-classical machine learning potential with variational quantum circuits

4 affiliations
  1. 1Department of Energy Science, Sungkyunkwan University 2066 Seobu-ro Suwon Gyeonggi-do 16419 Republic of Korea
  2. 2Department of Chemistry, Ulsan National Institute of Science and Technology 50 UNIST-gil Ulsan 44919 Republic of Korea
  3. 3Department of Energy, Sungkyunkwan University 2066 Seobu-ro Suwon Gyeonggi-do 16419 Republic of Korea
  4. 4Department of Quantum Information Engineering, Sungkyunkwan University 2066 Seobu-ro Suwon Gyeonggi-do 16419 Republic of Korea

Abstract

Machine learning potentials (MLPs) based on equivariant message-passing neural networks have achieved state-of-the-art accuracy for atomistic simulation, yet integration of variational quantum circuit (VQC) readouts into such frameworks remains unexplored for condensed-matter systems. Here we employ and characterize a hybrid quantum-classical MLP (HQC-MLP) by coupling an E(3)-equivariant message-passing network with VQC readout layers for predicting DFT energies and forces of liquid silicon (64 atoms, periodic cell, 2000–3000 K). We benchmark five VQC Ansatz architectures differing in entangling gate type (controlled-Z (CZ) vs. controlled-X (CX)) and circuit structure against a classical baseline. All HQC-MLP variants reproduce PBE-DFT energies with test-set RMSEs of ∼3.5 meV per atom when evaluated using exact expectation values, while entangled circuits exhibit faster training convergence in epoch count, although not in wall-clock time. CZ-based entangling gates achieve lower validation losses than their CX-based counterparts. This trend is consistent with two factors: the phase-preserving nature of CZ gates for amplitude-encoded features and their reduced effective entangling depth, since the final CZ layer commutes with the Pauli-Z measurement operators. Both molecular dynamics and Monte Carlo simulations accurately reproduce the DFT radial distribution function at 2000 K, despite the models being trained exclusively on 3000 K configurations; molecular dynamics also reproduces the Si–Si–Si bond-angle distribution. Finite-shot analysis further shows that achieving an RMSE of approximately 2 meV per atom relative to shot-noise-free expectation values requires about 5 × 104 shots per circuit execution. This corresponds to approximately 2.56 × 108 shots for a single configuration force evaluation, making direct force calculations prohibitively expensive on current noisy intermediate-scale quantum (NISQ) hardware and thereby motivating the proposed Monte Carlo protocol. Finally, the cutoff-radius sensitivity analysis reveals that HQC-MLP exhibits truncation behavior comparable to that of the classical MLP.

Related

  1. Quantum-Classical Effective Fragment Potential Embedding for Condensed-Phase Quantum Chemistry
    Federico Zahariev, Vassiliki-Alexandra Glezakou, Mark S. Gordon
    arXivPreprintCC BY 4.0Full textPhysical chemistryQuantum physics & quantum informationarXiv
  2. Machine learning for sample-based quantum diagonalization: a review of generative configuration recovery and the classical-simulability frontier
    Nicolás Bonilla Vargas
    arXivPreprintCC BY 4.0Full textPhysical chemistryQuantum physics & quantum informationarXiv
  3. Testing Small Parameterized Quantum Circuits for Molecular Property Prediction
    Bhagya Sri Uddandam, Chirag Girish Mokashi, Jose Rodrigo Arizpe Figueroa +7
    ChemRxivPreprintCC BYFull textGeneral chemistryChemRxiv
  4. Occupation-Space Variational Quantum Optimization
    Edison X. Salazar, Mario Piris
    ChemRxivPreprintCC BY-NC-NDGeneral chemistryChemRxiv
  5. Mitigating Errors in Variational Quantum Eigensolver Calculations through Extended Random Phase Approximation
    Michael Wagner, Pavel Begy, Jiří Brabec +4
    arXivPreprintarXiv licenceFull textPhysical chemistryQuantum physics & quantum informationarXiv
BibTeX
@article{willow2026hybrid,
  title = {{Hybrid quantum-classical machine learning potential with variational quantum circuits}},
  author = {Soohaeng Yoo Willow and D. ChangMo Yang and Kwang S. Kim and Chang Woo Myung},
  journal = {RSC Advances},
  year = {2026},
  month = oct,
  doi = {10.1039/d6ra07180e},
  url = {https://doi.org/10.1039/d6ra07180e},
  note = {CC BY 3.0},
}
RIS
TY  - JOUR
TI  - Hybrid quantum-classical machine learning potential with variational quantum circuits
AU  - Soohaeng Yoo Willow
AU  - D. ChangMo Yang
AU  - Kwang S. Kim
AU  - Chang Woo Myung
PY  - 2026
DA  - 2026/10/26
JO  - RSC Advances
DO  - 10.1039/d6ra07180e
UR  - https://doi.org/10.1039/d6ra07180e
AB  - Machine learning potentials (MLPs) based on equivariant message-passing neural networks have achieved state-of-the-art accuracy for atomistic simulation, yet integration of variational quantum circuit (VQC) readouts into such frameworks remains unexplored for condensed-matter systems. Here we employ and characterize a hybrid quantum-classical MLP (HQC-MLP) by coupling an E(3)-equivariant message-passing network with VQC readout layers for predicting DFT energies and forces of liquid silicon (64 atoms, periodic cell, 2000–3000 K). We benchmark five VQC Ansatz architectures differing in entangling gate type (controlled-Z (CZ) vs. controlled-X (CX)) and circuit structure against a classical baseline. All HQC-MLP variants reproduce PBE-DFT energies with test-set RMSEs of ∼3.5 meV per atom when evaluated using exact expectation values, while entangled circuits exhibit faster training convergence in epoch count, although not in wall-clock time. CZ-based entangling gates achieve lower validation losses than their CX-based counterparts. This trend is consistent with two factors: the phase-preserving nature of CZ gates for amplitude-encoded features and their reduced effective entangling depth, since the final CZ layer commutes with the Pauli-Z measurement operators. Both molecular dynamics and Monte Carlo simulations accurately reproduce the DFT radial distribution function at 2000 K, despite the models being trained exclusively on 3000 K configurations; molecular dynamics also reproduces the Si–Si–Si bond-angle distribution. Finite-shot analysis further shows that achieving an RMSE of approximately 2 meV per atom relative to shot-noise-free expectation values requires about 5 × 104 shots per circuit execution. This corresponds to approximately 2.56 × 108 shots for a single configuration force evaluation, making direct force calculations prohibitively expensive on current noisy intermediate-scale quantum (NISQ) hardware and thereby motivating the proposed Monte Carlo protocol. Finally, the cutoff-radius sensitivity analysis reveals that HQC-MLP exhibits truncation behavior comparable to that of the classical MLP.
SN  - 2046-2069
ER  -