Related papers: Efficient Magnus-type integrators for solar energy…
The integrating factor technique is widely used to solve numerically (in particular) the Schr\"odinger equation in the context of spectral methods. Here, we present an improvement of this method exploiting the freedom provided by the gauge…
Among the family of fourth-order time integration schemes, the two-stage Gauss--Legendre method, which is an implicit Runge--Kutta method based on collocation, is the only superconvergent. The computational cost of this implicit scheme for…
We change a previous time-stepping algorithm for solving a multi-scale Vlasov-Poisson system within a Particle-In-Cell method, in order to do accurate long time simulations. As an exponential integrator, the new scheme allows to use large…
A method of numerically solving the Maxwell equations is considered for modeling harmonic electromagnetic fields. The vector finite element method makes it possible to obtain a physically consistent discretization of the differential…
We take advantage of recent improvements in the grand canonical Hybrid Monte Carlo (HMC) algorithm, to perform a precision study of the single-particle gap in the hexagonal Hubbard model, with on-site electron-electron interactions. After…
Many complex systems can be accurately modeled as a set of coupled time-dependent partial differential equations (PDEs). However, solving such equations can be prohibitively expensive, easily taxing the world's largest supercomputers. One…
A fixed time-step variational integrator cannot preserve momentum, energy, and symplectic form simultaneously for nonintegrable systems. This barrier can be overcome by treating time as a discrete dynamic variable and deriving adaptive…
We study the effect of coupling an electronic reservoir to a Hubbard model and to a Dimer Hubbard Model. This is motivated by recent experiments on the effect of illumination on the insulator-metal transition in a vanadium oxides and…
We present a straightforward implementation scheme for solving the time dependent Schr\"odinger equation for systems described by the Hubbard Hamiltonian with time dependent hoppings. The computations can be performed for clusters of up to…
Assume interest is in sampling from a probability distribution $\mu$ defined on $(\mathsf{Z},\mathscr{Z})$. We develop a framework for sampling algorithms which takes full advantage of ODE numerical integrators, say…
A fast and robust Jacobian-free time-integration method - called Minimum-error Adaptation of a Chemical-Kinetic ODE Solver (MACKS) - for solving stiff ODEs pertaining to chemical-kinetics is proposed herein. The MACKS formulation is based…
We study the Hubbard model with time-reversal invariant flux and spin-orbit coupling and position-dependent onsite energies on the kagome lattice, using numerical and analytical methods. In particular, we perform calculations using real…
The dynamical low-rank approximation of time-dependent matrices is a low-rank factorization updating technique. It leads to differential equations for factors of the matrices, which need to be solved numerically. We propose and analyze a…
In this report we have analyzed a simple effective model for a description of magnetically ordered insulators. The Hamiltonian considered consists of the effective on-site interaction (U) and the intersite Ising-like magnetic exchange…
We develop an explicit, second-order, variational time integrator for full body dynamics that preserves the momenta of the continuous dynamics, such as linear and angular momenta, and exhibits near-conservation of total energy over…
Intermittent renewable energy resources like wind and solar pose great uncertainty of multiple time scales, from minutes to years, on the design and operation of power systems. Energy system optimization models have been developed to find…
We provide a unified, comprehensive treatment of all operators that contribute to the anti-ferromagnetic, ferromagnetic, and charge-density-wave structure factors and order parameters of the hexagonal Hubbard Model. We use the Hybrid Monte…
The semiclassical Schr\"odinger equation with time-dependent potentials is an important model to study electron dynamics under external controls in the mean-field picture. In this paper, we propose two multiscale finite element methods to…
It is often the case that, while the numerical solution of the non-linear dispersive equation $\mathrm{i}\partial_t u(t)=\mathcal{H}(u(t),t)u(t)$ represents a formidable challenge, it is fairly easy and cheap to solve closely related linear…
The Hubbard model is a standard theoretical tool for studying materials with strong electron-electron interactions, such as the cuprate superconductors. Unfortunately, interaction-driven phenomena such as the transition into the strongly…