Related papers: Truncation scheme of time-dependent density-matrix…
In this paper, the concept of matrix splitting is introduced to solve a large sparse ill-posed linear system via Tikhonov's regularization. In the regularization process, we convert the ill-posed system to a well-posed system. The…
We consider model order reduction of a nonlinear cable-mass system modeled by a 1D wave equation with interior damping and dynamic boundary conditions. The system is driven by a time dependent forcing input to a linear mass-spring system at…
In this paper we propose and analyze an energy stable numerical scheme for the Cahn-Hilliard equation, with second order accuracy in time and the fourth order finite difference approximation in space. In particular, the truncation error for…
We analyse two practical aspects that arise in the numerical solution of Hamilton-Jacobi-Bellman (HJB) equations by a particular class of monotone approximation schemes known as semi-Lagrangian schemes. These schemes make use of a wide…
Explicit discretizations of stochastic differential equations often encounter instability when the coefficients are not globally Lipschitz. The truncated schemes and tamed schemes have been proposed to handle this difficulty, but truncated…
A reaction-diffusion problem with a Caputo time derivative is considered. An integral discretization scheme on a graded mesh along with a decomposition of the exact solution is proposed. The truncation error estimate of the discretization…
We propose a density matrix renormalization group approach to tackle a two-state system coupled to a bosonic bath with continuous spectrum. In this approach, the optimized phonon scheme is applied to several hundred phonon modes which are…
In this study, we consider the numerical solution of large systems of linear equations obtained from the stochastic Galerkin formulation of stochastic partial differential equations. We propose an iterative algorithm that exploits the…
The higher-order tensor renormalization group is a tensor-network method providing estimates for the partition function and thermodynamical observables of classical and quantum systems in thermal equilibrium. At every step of the iterative…
We investigate linear dynamical systems consisting of ordinary differential equations with high dimensionality. Model order reduction yields alternative systems of much lower dimensions. However, a reduced system may be unstable, although…
Balanced truncation is one of the most common model order reduction schemes. In this paper, we study finite-frequency model order reduction (FF-MOR) problems of linear continuous-time systems within the framework of balanced truncation…
We have extended the density matrix renormalization group (DMRG) approach to two-fluid open many-fermion systems governed by complex-symmetric Hamiltonians. The applications are carried out for three- and four-nucleon (proton-neutron)…
Given a renormalization scheme, we show how to formulate a tractable convex relaxation of the set of feasible local density matrices of a many-body quantum system. The relaxation is obtained by introducing a hierarchy of constraints between…
A one-dimensional quantum many-body system consisting of particles confined in a harmonic potential and subject to finite-range two-body and three-body inverse-square interactions is introduced. The range of the interactions is set by…
We consider special upwinding Petrov-Galerkin discretizations for convection-diffusion problems. For the one dimensional case with a standard continuous linear element as the trial space and a special exponential bubble test space, we prove…
The aim of this paper is to develop fast second-order accurate difference schemes for solving one- and two-dimensional time distributed-order and Riesz space fractional diffusion equations. We adopt the same measures for one- and…
Novel constructions of empirical controllability and observability gramians for nonlinear systems for subsequent use in a balanced truncation style of model reduction are proposed. The new gramians are based on a generalisation of the…
Accurately simulating long-time dynamics of many-body systems is a challenge in both classical and quantum computing due to the accumulation of Trotter errors. While low-order Trotter-Suzuki decompositions are straightforward to implement,…
Based on recently derived exact stochastic Liouville-von Neumann equations, several strategies for the efficient simulation of open quantum systems are developed and tested on the spin-boson model. The accuracy and efficiency of these…
We present a mass lumping approach based on an isogeometric Petrov-Galerkin method that preserves higher-order spatial accuracy in explicit dynamics calculations irrespective of the polynomial degree of the spline approximation. To…