arXiv:2610.08497v1 Announce Type: cross Abstract: Transfer functions of Partial Differential Equations (PDEs) are irrational and difficult to obtain. Unlike for rational transfer functions, only a few methods are available for the control and analysis of irrational transfer functions. Thus, for robust analysis and control of PDEs, we need to use a rational approximation of the transfer function, preferably with provable $H_\infty$ error bounds. In this paper, we extend Krylov subspace model reduction to infinite-dimensional systems modeled by Partial Integral Equations (PIEs) -- an equivalent representation of PDEs. First, we propose a state transformation of the PIE system that allows us to compute a controllability matrix. Then, to obtain a rational approximation, we project the PIE states onto the vector space spanned by the truncated controllability matrix. Next, we show that sufficient conditions for exponential stability and well-posedness of the PIE system guarantee the exponential stability of the reduced-order system and provide constructive finite-frequency $H_\infty$ error bounds. Finally, we provide four illustrative numerical examples to demonstrate the performance of the proposed method.