
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 finitedifference models of largescale flexible
structures manufactured from homogeneous materials. Using the spatiotemporal
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 risksensitive performance
criterion for quantum control systems and its recent quadraticexponential
counterpart. We discuss a connection between the minimization of these cost
functionals and robustness with respect to uncertainty in systemenvironment
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 risksensitivity parameter in the context of
robust quantum control with entropy theoretic quantification of statistical
uncertainty in the systemfield state.

This paper is concerned with risksensitive 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 integrodifferential equation is obtained for the time evolution
of this quadraticexponential functional, which is compared with the original
quantum risksensitive performance criterion employed previously for
measurementbased quantum control and filtering problems. Using multipoint
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 infinitehorizon 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 higherorder 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 higherorder moments relevant to the
risksensitive criteria. The results of the paper are also demonstrated by a
numerical example and may find applications to coherent quantum risksensitive
control problems, where the plant and controller form a fully quantum
closedloop 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 statespace 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 onemode 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 measurementfree quantum
controller for a quantum plant so as to minimize a mean square cost for the
fully quantum closedloop system. The plant and controller are open quantum
systems interconnected through bosonic quantum fields. In comparison with the
observationactuation 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 statespace 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 normbalanced realizations of controllers. A linesearch
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 higherorder clusters to systems
with complex site configurations and shortrange 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
BethePeierls entropy approximation in the CVM hierarchy, the edgetosite
reduction of the entropy cumulants was studied previously. In extending the
earlier results, we develop a facetodiagonal reduction of the KramersWannier
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
systemobservation dynamics are governed by a Markovian HudsonParthasarathy
quantum stochastic differential equation driven by quantum Wiener processes of
bosonic fields in vacuum state. The Hamiltonian and systemfield coupling
operators, as functions of the system variables, are represented in a Weyl
quantization form. Using the WignerMoyal phasespace framework, we obtain a
stochastic integrodifferential equation for the posterior quasicharacteristic
function (QCF) of the system conditioned on the measurements. This equation is
a spatial Fourier domain representation of the BelavkinKushnerStratonovich
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 systemfield 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 secondorder 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 plantobserver coupling matrix in order to
mitigate the backaction of the observer on the covariance dynamics of the
plant. For the discounted mean square optimal CQF problem with penalized
backaction, we establish firstorder necessary conditions of optimality in the
form of algebraic matrix equations. By using the Hamiltonian structure of the
Heisenberg dynamics and related Liealgebraic 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
HudsonParthasarathy quantum stochastic differential equations driven by
external bosonic fields. The dependence of the Hamiltonian and the systemfield
coupling operators on the system variables is represented using the Weyl
functional calculus. This leads to an integrodifferential equation (IDE) for
the evolution of the quasicharacteristic 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 complexvalued
functions on finitedimensional Euclidean spaces and extends the WignerMoyal
phasespace approach for quantum stochastic systems. The dynamics of the QCF
and the related Wigner quasiprobability 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 secondorder 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 largescale 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 firstorder 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 semidefinite 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 GramCharlier 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 FokkerPlanckKolmogorov 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 HudsonParthasarathy 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 systemfield 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 systemfield 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
measurementfree cascade connection of quantum systems.

This paper is concerned with the coherent quantum filtering (CQF) problem,
where a quantum observer is cascaded in a measurementfree 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
HudsonParthasarathy quantum stochastic differential equations driven by
bosonic fields in vacuum state. These quantum dynamics are specified by the
Hamiltonians and systemfield coupling operators. We apply a recently proposed
transverse Hamiltonian variational method to the development of firstorder
necessary conditions of optimality for the CQF problem in a larger class of
observers. The latter is obtained by perturbing the Hamiltonian and
systemfield 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 largescale
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 PeresHorodeckiSimon 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 YoulaKu\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 systemfield 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 measurementfree quantum
controller for a quantum plant so as to minimize an infinitehorizon mean
square performance index for the fully quantum closedloop system. In
comparison with the observationactuation 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 steadystate quantum covariance matrix of the plantcontroller
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 finitedimensional 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 roundedoff planar rotations.

This paper is concerned with a problem of robust filtering for a
finitedimensional 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 worstcase
performance of an estimator is quantified by the maximum ratio of the
rootmeansquare (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
worstcase errortonoise RMS ratio. The estimator retains the general
structure of the Kalman filter, though with modified statespace 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 steadystate Kalman filter or the $H_{\infty}$optimal
estimator, respectively.

This paper is concerned with integrodifferential 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 quadraticexponential
moments which are relevant to risksensitive 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 onedimensional 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
largescale 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. 270275). 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 (measurementfree)
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
feedbackfree 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 closedform smooth solution is obtained for
a HamiltonJacobi 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 worstcase disturbance attenuation capabilities of the system
with respect to the statistically uncertain random inputs are quantified by the
aanisotropic norm which is an appropriately constrained operator norm of the
system. We establish an anisotropic norm bounded real lemma which provides a
statespace criterion for the aanisotropic 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 Hinfinitycontrol theory. We also
provide a necessary background on the anisotropybased robust performance
analysis.