Dynamical low-rank approximation in the Tucker tensor format of given large time-dependent tensors and of tensor differential equations is the subject of this paper. In particular, a discrete time integration method for rank-constrained Tucker tensors is presented and analyzed. It extends the known projector-splitting integrator for dynamical low-rank approximation of matrices to Tucker tensors and is shown to inherit the same favorable properties. The integrator is based on iteratively applying the matrix projector-splitting integrator to tensor unfoldings but with inexact solution in a substep. It has the property that it reconstructs time-dependent Tucker tensors of the given rank exactly. The integrator is also shown to be robust to the presence of small singular values in the tensor unfoldings. Numerical examples with time-dependent problems from quantum physics and tensor optimization methods illustrate our theoretical results.
Differential-algebraic equations (DAEs) arise naturally in constrained dynamical systems, but their algebraic constraints and hidden compatibility conditions make them more subtle than standard ordinary differential equations. This paper initiates a quantum-algorithmic study of constrained linear DAEs. We introduce a dilation framework that embeds the generally non-Hermitian constrained evolution into a projected Schrödinger-type dynamics on an enlarged Hilbert space, idtdΨ(t)=PHPΨ(t), where H is Hermitian and P is the orthogonal projector onto the lifted constraint subspace. This identifies the DAE evolution with a quantum Zeno-type dynamics and enables the use of block encodings, QSVT-based projector construction, and Hamiltonian simulation.
We apply the framework to structure-preserving discretizations of the unsteady Stokes equations, where the pressure enforces the discrete incompressibility constraint. For Stokes, the Zeno-reduced generator has the projected square factorization Sh=−ΠhΔhΠh=(GhΠh)†(GhΠh), which can be represented through a Gaussian moment dilation and implemented as a Gaussian superposition of unitary Zeno evolutions generated by a first-order square-root Hamiltonian. In the generic sparse-access model, this gives a simulation-stage cost O(h−2t), up to the usual postselection factor for preparing the normalized dissipative state. The results provide a first step toward understanding the intersection of quantum algorithms, DAEs, constrained PDE dynamics, and square-root Gaussian dilations.
In quantum theory, equilibrium statistical mechanics is usually formulated through the canonical ensemble, whose privileged status is tied to the Euclidean continuation of time evolution. The microcanonical ensemble, by contrast, is commonly introduced as a separate spectral construction. In this work we show that this asymmetry is representational rather than structural. We formulate the system in an extended Hilbert space in which time is promoted to an auxiliary clock degree of freedom and physical states are selected by a reparametrization-invariant constraint operator C^=P^T+H^. The corresponding projector δ(C^) provides a single unified object from which both canonical and microcanonical ensembles emerge as complementary projections in the clock sector. In the clock-time representation, a purely imaginary clock separation yields the Euclidean kernel and the canonical partition function. In the conjugate clock-energy representation, the same projector reduces to the spectral operator δ(H^−E) and hence to the microcanonical density of states. The main consequence is structural: canonical and microcanonical statistics need not be introduced as independent constructions, since both are already encoded in the same constrained quantum dynamics.