Related papers: A new method for probing the late-time dynamics in…
The Gaussian expansion has been developed since early 80s as a powerful analytical method, which enables nonperturbative studies of various systems using `perturbative' calculations. Recently the method has been used to suggest that 4d…
The effective dynamics of the low-frequency modes is derived for the O(N) symmetric scalar field theory in the broken symmetry phase. The effect of the high-frequency fluctuations is taken into account at one-loop level exactly. A new…
We propose a modified power method for computing the subdominant eigenvalue $\lambda_2$ of a matrix or continuous operator. Here we focus on defining simple Monte Carlo methods for its application. The methods presented use random walkers…
We present preliminary numerical results from a lattice study of the two-dimensional O(3) non-linear sigma model. In the continuum this model possesses N=2 supersymmetry. The lattice formulation we use retains an exact (twisted)…
There have been several algorithms designed to optimise matrix multiplication. From schoolbook method with complexity $O(n^3)$ to advanced tensor-based tools with time complexity $O(n^{2.3728639})$ (lowest possible bound achieved), a lot of…
In linear stability analysis of field quantities described by partial differential equations, the well-established classical theory is all but impossible to apply to concrete problems in its entirety even for uniform backgrounds when the…
Micromagnetics depends on high-fidelity numerical methods for magnetization dynamics. This work proposes a third-order temporal accuracy scheme for the Landau-Lifshitz-Gilbert equation, addressing accuracy-efficiency trade-offs in existing…
It has long been conjectured that certain supersymmetric Yang-Mills (SYM) theories provide us with nonperturbative formulations of the string/M-theory. Although the supersymmetry (SUSY) on lattice is notoriously difficult in general, for a…
The second moment method has always been an effective tool to lower bound the satisfiability threshold of many random constraint satisfaction problems. However, the calculation is usually hard to carry out and as a result, only some loose…
The lattice provides a powerful tool to non-perturbatively investigate strongly coupled supersymmetric Yang-Mills (SYM) theories. The pure SU(2) SYM theory with one supercharge is simulated on large lattices with small Majorana gluino…
An important paradigm in the understanding of mixing times of Glauber dynamics for spin systems is the correspondence between spatial mixing properties of the models and bounds on the mixing time of the dynamics. This includes, in…
We apply the bootstrap technique to find the moments of certain multi-trace and multi-matrix random matrix models suggested by noncommutative geometry. Using bootstrapping we are able to find the relationships between the coupling constant…
We modified the recently proposed multicanonical MC algorithm for the case of a magnetic field driven order--order phase transition. We test this {\it multimagnetic} Monte Carlo algorithm for the D=2 Ising model at $\beta=0.5$ and simulate…
This paper introduces a novel method for approximating the dynamics of a large autonomous system projected onto a fixed subspace. The core contribution is a novel recursive algorithm to construct an effective time-dependent generator that…
It is shown that perturbation theory in $2D$ nonlinear $\sigma$-models as well gauge theories in dimension $D\geq 2$ produces answers that depend on boundary conditions even after the infinite volume limit has been taken. This unphysical…
Matrix Models are the most effective way to describe strongly interacting systems with many degrees of freedom. They have proven successful in describing very different settings, from nuclei spectra to conduction in mesoscopic systems, from…
The problem of late time instability in time domain integral equations for electromagnetics is longstanding. While several techniques have been suggested for addressing this problem, they either require impractically high degrees of freedom…
Using Monte Carlo simulations we study two-dimensional prey-predator systems. Measuring the variance of densities of prey and predators on the triangular lattice and on the lattice with eight neighbours, we conclude that temporal…
An efficient Monte Carlo algorithm for the simulation of spin models with long-range interactions is discussed. Its central feature is that the number of operations required to flip a spin is independent of the number of interactions…
This project investigates the applicability of quasi-Monte Carlo methods to Euclidean lattice systems in order to improve the asymptotic error scaling of observables for such theories. The error of an observable calculated by averaging over…