Checking Continuous Stochastic Logic against Quantum Continuous-Time Markov Chains
Ming Xu, Jingyi Mei, Ji Guan, Yuxin Deng, Nengkun Yu
TL;DR
The paper addresses verifying CSL properties of quantum continuous-time Markov chains, models that couple a classical subsystem with a quantum subsystem under Lindblad dynamics. It introduces an algebraic, symbolically computable procedure using projection, matrix exponentiation, and definite integration to evaluate multiphase until path probabilities, proving CSL decidability and achieving polynomial-time performance in the model encoding size while acknowledging exponential dependence on qubit count. A numerical speed-up based on scaling and squaring with Padé approximations and Riemann sums is proposed to make the approach practical for larger systems. An Apollonian-network running example demonstrates the method, and an implementation in Mathematica/MATLAB showcases the feasibility of exact and approximate CSL checking on quantum CTMCs, bridging formal verification with open quantum dynamics. The work lays groundwork for future model checking of quantum CTMDPs and highlights the need to address optimization over density operators in quantum settings.
Abstract
Verifying quantum systems has attracted a lot of interest in the last decades.In this paper, we study the quantitative model-checking of quantum continuous-time Markov chains (quantum CTMCs). The branching-time properties of quantum CTMCs are specified by continuous stochastic logic (CSL), which is well-known for verifying real-time systems, including classical CTMCs. The core of checking the CSL formulas lies in tackling multiphase until formulas. We develop an algebraic method using proper projection, matrix exponentiation, and definite integration to symbolically calculate the probability measures of path formulas. Thus the decidability of CSL is established. To be efficient, numerical methods are incorporated to guarantee that the time complexity is polynomial in the encoding size of the input model and linear in the size of the input formula. A running example of Apollonian networks is further provided to demonstrate our method.
