Junyuan He

Junyuan He, Zhonghao Sun, Jizu Huang

Tree tensor networks (TTNs) provide a compact and structured representation of high-dimensional data, making them valuable in various areas of computational mathematics and physics. In this paper, we present a rigorous mathematical framework for expressing high-order derivatives of functional TTNs, both with or without constraints. Our framework decomposes the total derivative of a given TTN into a summation of TTNs, each corresponding to the partial derivatives of the original TTN. Using this decomposition, we derive the Taylor expansion of vector-valued functions subject to ordinary differential equation constraints or algebraic constraints imposed by Runge--Kutta (RK) methods. As a concrete application, we employ this framework to construct order conditions for RK methods. Due to the intrinsic tensor properties of partial derivatives and the separable tensor structure in RK methods, the Taylor expansion of numerical solutions can be obtained in a manner analogous to that of exact solutions using tensor operators. This enables the order conditions of RK methods to be established by directly comparing the Taylor expansions of the exact and numerical solutions, eliminating the need for mathematical induction. For a given function $\vector{f}$, we derive sharper order conditions that go beyond the classical ones, enabling the identification of situations where a standard RK scheme of order $p$ achieves unexpectedly higher convergence order for the particular function. These results establish new connections between tensor network theory and classical numerical methods, potentially opening new avenues for both analytical exploration and practical computation.

Junyuan He, Jizu Huang

Constructing explicit Runge--Kutta (ERK) methods with as few stages as possible for a given order is a classical problem in numerical analysis. In this work, we introduce a $Q$/$D$-space framework of sufficient order conditions for ERK methods. This framework generalizes Butcher's classical simplifying assumptions by reformulating them in terms of simplified $Q$- and $D$-spaces defined through their residual vectors. It yields sufficient conditions which, together with $B(p)$, ensure order $p$. It also leads to a recursive construction procedure for ERK methods of arbitrary even order, in which the Butcher coefficients are obtained from two structured linear systems. For every even order $p\ge 4$, the construction produces ERK methods with stage number $s(p)=(p^2-2p+8)/4$. This stage count has the same leading term as that of the classical Gragg families, while improving the linear term. The free parameters retained by the construction further provide a systematic framework for designing methods with enhanced stability and short-time accuracy.