• This paper is concerned with networks of interacting linear systems at sites of a multidimensional lattice. The systems are governed by linear ODEs with constant coefficients driven by external inputs, and their internal dynamics and coupling with the other component systems are translation invariant. Such systems occur, for example, in finite-difference models of large-scale flexible structures manufactured from homogeneous materials. Using the spatio-temporal transfer function of this translation invariant network, we establish conditions for its positive realness in the sense of energy dissipation. The latter is formulated in terms of block Toeplitz bilinear forms of the input and output variables of the composite system. We also discuss quadratic stability of the network in isolation from the environment and phonon theoretic dispersion relations.
  • This paper considers links between the original risk-sensitive performance criterion for quantum control systems and its recent quadratic-exponential counterpart. We discuss a connection between the minimization of these cost functionals and robustness with respect to uncertainty in system-environment quantum states whose deviation from a nominal state is described in terms of the quantum relative entropy. These relations are similar to those in minimax LQG control for classical systems. The results of the paper can be of use in providing a rational choice of the risk-sensitivity parameter in the context of robust quantum control with entropy theoretic quantification of statistical uncertainty in the system-field state.
  • This paper is concerned with risk-sensitive performance analysis for linear quantum stochastic systems interacting with external bosonic fields. We consider a cost functional in the form of the exponential moment of the integral of a quadratic polynomial of the system variables over a bounded time interval. An integro-differential equation is obtained for the time evolution of this quadratic-exponential functional, which is compared with the original quantum risk-sensitive performance criterion employed previously for measurement-based quantum control and filtering problems. Using multi-point Gaussian quantum states for the past history of the system variables and their first four moments, we discuss a quartic approximation of the cost functional and its infinite-horizon asymptotic behaviour. The computation of the asymptotic growth rate of this approximation is reduced to solving two algebraic Lyapunov equations. We also outline further approximations of the cost functional, based on higher-order cumulants and their growth rates, together with large deviations estimates. For comparison, an auxiliary classical Gaussian Markov diffusion process is considered in a complex Euclidean space which reproduces the quantum system variables at the level of covariances but has different higher-order moments relevant to the risk-sensitive criteria. The results of the paper are also demonstrated by a numerical example and may find applications to coherent quantum risk-sensitive control problems, where the plant and controller form a fully quantum closed-loop system, and other settings with nonquadratic cost functionals.
  • This paper is concerned with the generation of Gaussian invariant states in cascades of open quantum harmonic oscillators governed by linear quantum stochastic differential equations. We carry out infinitesimal perturbation analysis of the covariance matrix for the invariant Gaussian state of such a system and the related purity functional subject to inaccuracies in the energy and coupling matrices of the subsystems. This leads to the problem of balancing the state-space realizations of the component oscillators through symplectic similarity transformations in order to minimize the mean square sensitivity of the purity functional to small random perturbations of the parameters. This results in a quadratic optimization problem with an effective solution in the case of cascaded one-mode oscillators, which is demonstrated by a numerical example. We also discuss a connection of the sensitivity index with classical statistical distances and outline infinitesimal perturbation analysis for translation invariant cascades of identical oscillators. The findings of the paper are applicable to robust state generation in quantum stochastic networks.
  • This paper is concerned with coherent quantum linear quadratic Gaussian (CQLQG) control. The problem is to find a stabilizing measurement-free quantum controller for a quantum plant so as to minimize a mean square cost for the fully quantum closed-loop system. The plant and controller are open quantum systems interconnected through bosonic quantum fields. In comparison with the observation-actuation structure of classical controllers, coherent quantum feedback is less invasive to the quantum dynamics. The plant and controller variables satisfy the canonical commutation relations (CCRs) of a quantum harmonic oscillator and are governed by linear quantum stochastic differential equations (QSDEs). In order to correspond to such oscillators, these QSDEs must satisfy physical realizability (PR) conditions in the form of quadratic constraints on the state-space matrices, reflecting the CCR preservation in time. The symmetry of the problem is taken into account by introducing equivalence classes of coherent quantum controllers generated by symplectic similarity transformations. We discuss a modified gradient flow, which is concerned with norm-balanced realizations of controllers. A line-search gradient descent algorithm with adaptive stepsize selection is proposed for the numerical solution of the CQLQG control problem. The algorithm finds a local minimum of the LQG cost over the parameters of the Hamiltonian and coupling operators of a stabilizing coherent quantum controller, thus taking the PR constraints into account. A convergence analysis of the algorithm is presented. Numerical examples of designing locally optimal CQLQG controllers are provided in order to demonstrate the algorithm performance.
  • The paper is concerned with interacting particle systems on the simple cubic lattice obeying the nearest neighbour exclusion (NNE). This constraint forbids any two neighbouring sites of the lattice to be simultaneously occupied, thus reducing the set of admissible configurations for the cubic cell and its subclusters such as edges and faces. This reduction extends applicability of Kikuchi's Cluster Variation Method (CVM) with higher-order clusters to systems with complex site configurations and short-range ordering, which would be impractical beyond the NNE framework because of the "curse of dimensionality". For edges of the cubic cell, which are the operational clusters of the Bethe-Peierls entropy approximation in the CVM hierarchy, the edge-to-site reduction of the entropy cumulants was studied previously. In extending the earlier results, we develop a face-to-diagonal reduction of the Kramers-Wannier entropy approximation of the CVM in the NNE setting. We also outline an application of the resulting approximation to thermodynamic modeling of disordered condensed media, such as liquid silicates, and discuss combinatorial and numerical aspects of the implementation of this approach.
  • This paper is concerned with a filtering problem for a class of nonlinear quantum stochastic systems with multichannel nondemolition measurements. The system-observation dynamics are governed by a Markovian Hudson-Parthasarathy quantum stochastic differential equation driven by quantum Wiener processes of bosonic fields in vacuum state. The Hamiltonian and system-field coupling operators, as functions of the system variables, are represented in a Weyl quantization form. Using the Wigner-Moyal phase-space framework, we obtain a stochastic integro-differential equation for the posterior quasi-characteristic function (QCF) of the system conditioned on the measurements. This equation is a spatial Fourier domain representation of the Belavkin-Kushner-Stratonovich stochastic master equation driven by the innovation process associated with the measurements. We also discuss a more specific form of the posterior QCF dynamics in the case of linear system-field coupling and outline a Gaussian approximation of the posterior quantum state.
  • This paper is concerned with quantum harmonic oscillators consisting of a quantum plant and a directly coupled coherent quantum observer. We employ discounted quadratic performance criteria in the form of exponentially weighted time averages of second-order moments of the system variables. A coherent quantum filtering (CQF) problem is formulated as the minimization of the discounted mean square of an estimation error, with which the dynamic variables of the observer approximate those of the plant. The cost functional also involves a quadratic penalty on the plant-observer coupling matrix in order to mitigate the back-action of the observer on the covariance dynamics of the plant. For the discounted mean square optimal CQF problem with penalized back-action, we establish first-order necessary conditions of optimality in the form of algebraic matrix equations. By using the Hamiltonian structure of the Heisenberg dynamics and related Lie-algebraic techniques, we represent this set of equations in a more explicit form in the case of equally dimensioned plant and observer.
  • This paper considers open quantum systems whose dynamic variables satisfy canonical commutation relations and are governed by Markovian Hudson-Parthasarathy quantum stochastic differential equations driven by external bosonic fields. The dependence of the Hamiltonian and the system-field coupling operators on the system variables is represented using the Weyl functional calculus. This leads to an integro-differential equation (IDE) for the evolution of the quasi-characteristic function (QCF) which encodes the dynamics of mixed moments of the system variables. Unlike quantum master equations for reduced density operators, this IDE involves only complex-valued functions on finite-dimensional Euclidean spaces and extends the Wigner-Moyal phase-space approach for quantum stochastic systems. The dynamics of the QCF and the related Wigner quasi-probability density function (QPDF) are discussed in more detail for the case when the coupling operators depend linearly on the system variables and the Hamiltonian has a nonquadratic part represented in the Weyl quantization form. For this class of quantum stochastic systems, we also consider an approximate computation of invariant states and discuss the deviation from Gaussian quantum states in terms of the $\chi^2$-divergence (or the second-order Renyi relative entropy) applied to the QPDF. The results of the paper may find applications to investigating different aspects of the moment stability, relaxation dynamics and invariant states in open quantum systems.
  • This paper is concerned with coherent quantum control design for translation invariant networks of identical quantum stochastic systems subjected to external quantum noise. The network is modelled as an open quantum harmonic oscillator and is governed by a set of linear quantum stochastic differential equations. The dynamic variables of this quantum plant satisfy the canonical commutation relations. Similar large-scale systems can be found, for example, in quantum metamaterials and optical lattices. The problem under consideration is to design a stabilizing decentralized coherent quantum controller in the form of another translation invariant quantum system, directly coupled to the plant, so as to minimize a weighted mean square functional of the dynamic variables of the interconnected networks. We consider this problem in the thermodynamic limit of infinite network size and present first-order necessary conditions for optimality of the controller.
  • This paper outlines an approach to the approximation of probability density functions by quadratic forms of weighted orthonormal basis functions with positive semi-definite Hermitian matrices of unit trace. Such matrices are called stochastic density matrices in order to reflect an analogy with the quantum mechanical density matrices. The SDM approximation of a PDF satisfies the normalization condition and is nonnegative everywhere in contrast to the truncated Gram-Charlier and Edgeworth expansions. For bases with an algebraic structure, such as the Hermite polynomial and Fourier bases, the SDM approximation can be chosen so as to satisfy given moment specifications and can be optimized using a quadratic proximity criterion. We apply the SDM approach to the Fokker-Planck-Kolmogorov PDF dynamics of Markov diffusion processes governed by nonlinear stochastic differential equations. This leads to an ordinary differential equation for the SDM dynamics of the approximating PDF. As an example, we consider the Smoluchowski SDE on a multidimensional torus.
  • This paper is concerned with variational methods for nonlinear open quantum systems with Markovian dynamics governed by Hudson-Parthasarathy quantum stochastic differential equations. The latter are driven by quantum Wiener processes of the external boson fields and are specified by the system Hamiltonian and system-field coupling operators. We consider the system response to perturbations of these energy operators and introduce a transverse Hamiltonian which encodes the propagation of the perturbations through the unitary system-field evolution. This provides a tool for the infinitesimal perturbation analysis and development of optimality conditions for coherent quantum control problems. We apply the transverse Hamiltonian variational technique to a mean square optimal coherent quantum filtering problem for a measurement-free cascade connection of quantum systems.
  • This paper is concerned with the coherent quantum filtering (CQF) problem, where a quantum observer is cascaded in a measurement-free fashion with a linear quantum plant so as to minimize a mean square error of estimating the plant variables of interest. Both systems are governed by Markovian Hudson-Parthasarathy quantum stochastic differential equations driven by bosonic fields in vacuum state. These quantum dynamics are specified by the Hamiltonians and system-field coupling operators. We apply a recently proposed transverse Hamiltonian variational method to the development of first-order necessary conditions of optimality for the CQF problem in a larger class of observers. The latter is obtained by perturbing the Hamiltonian and system-field coupling operators of a linear coherent quantum observer along linear combinations of unitary Weyl operators, whose role here resembles that of the needle variations in the Pontryagin minimum principle. We show that if the observer is a stationary point of the performance functional in the class of linear observers, then it is also a stationary point with respect to the Weyl variations in the larger class of nonlinear observers.
  • This paper is concerned with a translation invariant network of identical quantum stochastic systems subjected to external quantum noise. Each node of the network is directly coupled to a finite number of its neighbours. This network is modelled as an open quantum harmonic oscillator and is governed by a set of linear quantum stochastic differential equations. The dynamic variables of the network satisfy the canonical commutation relations. Similar large-scale networks can be found, for example, in quantum metamaterials and optical lattices. Using spatial Fourier transform techniques, we obtain a sufficient condition for stability of the network in the case of finite interaction range, and consider a mean square performance index for the stable network in the thermodynamic limit. The Peres-Horodecki-Simon separability criterion is employed in order to obtain sufficient and necessary conditions for quantum entanglement of bipartite systems of nodes of the network in the Gaussian invariant state. The results on stability and entanglement are extended to the infinite chain of the linear quantum systems by letting the number of nodes go to infinity. A numerical example is provided to illustrate the results.
  • This paper is concerned with application of the classical Youla-Ku\v{c}era parameterization to finding a set of linear coherent quantum controllers that stabilize a linear quantum plant. The plant and controller are assumed to represent open quantum harmonic oscillators modelled by linear quantum stochastic differential equations. The interconnections between the plant and the controller are assumed to be established through quantum bosonic fields. In this framework, conditions for the stabilization of a given linear quantum plant via linear coherent quantum feedback are addressed using a stable factorization approach. The class of stabilizing quantum controllers is parameterized in the frequency domain. Also, this approach is used in order to formulate coherent quantum weighted $H_2$ and $H_\infty$ control problems for linear quantum systems in the frequency domain. Finally, a projected gradient descent scheme is proposed to solve the coherent quantum weighted $H_2$ control problem.
  • This paper is concerned with open quantum systems whose dynamic variables satisfy canonical commutation relations and are governed by quantum stochastic differential equations. The latter are driven by quantum Wiener processes which represent external boson fields. The system-field coupling operators are linear functions of the system variables. The Hamiltonian consists of a nominal quadratic function of the system variables and an uncertain perturbation which is represented in a Weyl quantization form. Assuming that the nominal linear quantum system is stable, we develop sufficient conditions on the perturbation of the Hamiltonian which guarantee robust mean square stability of the perturbed system. Examples are given to illustrate these results for a class of Hamiltonian perturbations in the form of trigonometric polynomials of the system variables.
  • This paper is concerned with the Coherent Quantum Linear Quadratic Gaussian (CQLQG) control problem of finding a stabilizing measurement-free quantum controller for a quantum plant so as to minimize an infinite-horizon mean square performance index for the fully quantum closed-loop system. In comparison with the observation-actuation structure of classical controllers, the coherent quantum feedback is less invasive to the quantum dynamics and quantum information. Both the plant and the controller are open quantum systems whose dynamic variables satisfy the canonical commutation relations (CCRs) of a quantum harmonic oscillator and are governed by linear quantum stochastic differential equations (QSDEs). In order to correspond to such oscillators, these QSDEs must satisfy physical realizability (PR) conditions, which are organised as quadratic constraints on the controller matrices and reflect the preservation of CCRs in time. The CQLQG problem is a constrained optimization problem for the steady-state quantum covariance matrix of the plant-controller system satisfying an algebraic Lyapunov equation. We propose a gradient descent algorithm equipped with adaptive stepsize selection for the numerical solution of the problem. The algorithm finds a local minimum of the LQG cost over the parameters of the Hamiltonian and coupling operators of a stabilizing PR quantum controller, thus taking the PR constraints into account. A convergence analysis of the proposed algorithm is presented. A numerical example of a locally optimal CQLQG controller design is provided to demonstrate the algorithm performance.
  • The roundoff errors in computer simulations of continuous dynamical systems, caused by finiteness of machine arithmetic, can lead to qualitative discrepancies between phase portraits of the resulting spatially discretized systems and the original systems. These effects can be modelled on a multidimensional integer lattice by using a dynamical system obtained by composing the transition operator of the original system with a quantizer. Such models manifest pseudorandomness which can be studied using a rigorous probability theoretic approach. To this end, the lattice $\mathbb{Z}^n$ is endowed with a class of frequency measurable subsets and a spatial frequency functional as a finitely additive probability measure on them. Using a multivariate version of Weyl's equidistribution criterion, we introduce an algebra of frequency measurable quasiperiodic subsets of the lattice. This approach is applied to quantized linear systems with the transition operator $R \circ L$, where $L$ is a nonsingular matrix of the original linear system in $\mathbb{R}^n$, and the map $R$ commutes with the additive group of translations of the lattice. For almost every $L$, the events associated with the deviation of trajectories of the quantized and original systems are frequency measurable quasiperiodic subsets of the lattice whose frequencies involve geometric probabilities on finite-dimensional tori. Using the skew products of measure preserving toral automorphisms, we prove mutual independence and uniform distribution of the quantization errors and investigate statistical properties of invertibility loss for the quantized linear system, extending V.V.Voevodin's results. When $L$ is similar to an orthogonal matrix, we establish a functional central limit theorem for the deviations of trajectories of the quantized and original systems. These results are demonstrated for rounded-off planar rotations.
  • This paper is concerned with a problem of robust filtering for a finite-dimensional linear discrete time invariant system with two output signals, one of which is directly observed while the other has to be estimated. The system is assumed to be driven by a random disturbance produced from the Gaussian white noise sequence by an unknown shaping filter. The worst-case performance of an estimator is quantified by the maximum ratio of the root-mean-square (RMS) value of the estimation error to that of the disturbance over stationary Gaussian disturbances whose mean anisotropy is bounded from above by a given parameter $a \ge 0$. The mean anisotropy is a combined entropy theoretic measure of temporal colouredness and spatial "nonroundness" of a signal. We construct an $a$-anisotropic estimator which minimizes the worst-case error-to-noise RMS ratio. The estimator retains the general structure of the Kalman filter, though with modified state-space matrices. Computing the latter is reduced to solving a set of two coupled algebraic Riccati equations and an equation involving the determinant of a matrix. In two limiting cases, where $a = 0$ or $a \to +\infty$, the $a$-anisotropic estimator leads to the standard steady-state Kalman filter or the $H_{\infty}$-optimal estimator, respectively.
  • This paper is concerned with integro-differential identities which are known in statistical signal processing as Price's theorem for expectations of nonlinear functions of jointly Gaussian random variables. We revisit these relations for classical variables by using the Frechet differentiation with respect to covariance matrices, and then show that Price's theorem carries over to a quantum mechanical setting. The quantum counterpart of the theorem is established for Gaussian quantum states in the framework of the Weyl functional calculus for quantum variables satisfying the Heisenberg canonical commutation relations. The quantum mechanical version of Price's theorem relates the Frechet derivative of the generalized moment of such variables with respect to the real part of their quantum covariance matrix with other moments. As an illustrative example, we consider these relations for quadratic-exponential moments which are relevant to risk-sensitive quantum control.
  • This paper is concerned with translation invariant networks of linear quantum stochastic systems with nearest neighbour interaction mediated by boson fields. The systems are associated with sites of a one-dimensional chain or a multidimensional lattice and are governed by coupled linear quantum stochastic differential equations (QSDEs). Such interconnections of open quantum systems are relevant, for example, to the phonon theory of crystalline solids, atom trapping in optical lattices and quantum metamaterials. In order to represent a large-scale open quantum harmonic oscillator, the coefficients of the coupled QSDEs must satisfy certain physical realizability conditions. These are established in the form of matrix algebraic equations for the parameters of an individual building block of the network and its interaction with the neighbours and external fields. We also discuss the computation of mean square performance functionals with block Toeplitz weighting matrices for such systems in the thermodynamic limit per site for unboundedly increasing fragments of the lattice.
  • A Modified Associate Formalism is proposed for thermodynamic modelling of solution phases. The approach is free from the entropy paradox described by L\"{u}ck et al. (Z. Metallkd. 80 (1989) pp. 270--275). The model is considered in its general form for an arbitrary number of solution components and an arbitrary size of associates. Asymptotic behaviour of chemical activities of solution components in binary dilute solutions is also investigated.
  • The paper is concerned with a problem of coherent (measurement-free) filtering for physically realizable (PR) linear quantum plants. The state variables of such systems satisfy canonical commutation relations and are governed by linear quantum stochastic differential equations, dynamically equivalent to those of an open quantum harmonic oscillator. The problem is to design another PR quantum system, connected unilaterally to the output of the plant and playing the role of a quantum filter, so as to minimize a mean square discrepancy between the dynamic variables of the plant and the output of the filter. This coherent quantum filtering (CQF) formulation is a simplified feedback-free version of the coherent quantum LQG control problem which remains open despite recent studies. The CQF problem is transformed into a constrained covariance control problem which is treated by using the Frechet differentiation of an appropriate Lagrange function with respect to the matrices of the filter.
  • This paper is concerned with a dissipativity theory for dynamical systems governed by linear Ito stochastic differential equations driven by random noise with an uncertain drift. The deviation of the noise from a standard Wiener process in the nominal model is quantified by relative entropy. We discuss a dissipation inequality for the noise relative entropy supply. The problem of minimizing the supply required to drive the system between given Gaussian state distributions over a specified time horizon is considered. This problem, known in the literature as the Schroedinger bridge, was treated previously in the context of reciprocal processes. A closed-form smooth solution is obtained for a Hamilton-Jacobi equation for the minimum required relative entropy supply by using nonlinear algebraic techniques.
  • We consider a finite horizon linear discrete time varying system whose input is a random noise with an imprecisely known probability law. The statistical uncertainty is described by a nonnegative parameter a which constrains the anisotropy of the noise as an entropy theoretic measure of deviation of the actual noise distribution from Gaussian white noise laws with scalar covariance matrices. The worst-case disturbance attenuation capabilities of the system with respect to the statistically uncertain random inputs are quantified by the a-anisotropic norm which is an appropriately constrained operator norm of the system. We establish an anisotropic norm bounded real lemma which provides a state-space criterion for the a-anisotropic norm of the system not to exceed a given threshold. The criterion is organized as an inequality on the determinants of matrices associated with a difference Riccati equation and extends the Bounded Real Lemma of the H-infinity-control theory. We also provide a necessary background on the anisotropy-based robust performance analysis.