Related papers: Testing the stochastic LapH method in the twisted …
There is a great need for accurate and efficient computational approaches that can account for both the discrete and stochastic nature of chemical interactions as well as spatial inhomogeneities and diffusion. This is particularly true in…
We modify the coupling method established in [22, 20] and develop a technique to prove the exponential mixing of a 2D stochastic system forced by degenerate Levy noises. In particular, these Levy noises include $\alpha$-stable noises (0 <…
We investigate the utility of partially twisted boundary conditions in lattice calculations of meson observables. For dynamical simulations, we show that the pion dispersion relation is modified by volume effects. In the isospin limit, we…
We introduce a $Z_2$ noise for the stochastic estimation of matrix inversion and discuss its superiority over other noises including the Gaussian noise. This algorithm is applied to the calculation of quark loops in lattice quantum…
In this paper, we investigate the use of a mass lumped fully explicit time stepping scheme for the discretisation of the wave equation with underlying material parameters that vary at arbitrarily fine scales. We combine the leapfrog scheme…
We present a first scaling test of twisted mass QCD with pure Wilson quarks for a twisting angle of pi/2. We have computed the vector meson mass and the pseudoscalar decay constant for different values of beta at fixed value of r_0 m_PS.…
We consider the problem of efficiently simulating stochastic models of chemical kinetics. The Gillespie Stochastic Simulation algorithm (SSA) is often used to simulate these models, however, in many scenarios of interest, the computational…
In this paper we analyze a mixed displacement-pseudostress formulation for the elasticity eigenvalue problem. We propose a finite element method to approximate the pseudostress tensor with Raviart-Thomas elements and the displacement with…
In this paper, we extend the dynamical low-rank approximation method to the space of finite signed measures. Under this framework, we derive stochastic low-rank dynamics for stochastic differential equations (SDEs) coming from classical…
The first part of the cumulative thesis contains the numerical analysis of different $hp$-finite element discretizations related to two different weak formulations of a model problem in elastoplasticity with linearly kinematic hardening.…
In this article, we analyze the SPICE method developed in [1], and establish its connections with other standard sparse estimation methods such as the Lasso and the LAD-Lasso. This result positions SPICE as a computationally efficient…
We calculate finite volume effects on the pion masses and decay constant in twisted mass lattice QCD (tmLQCD) at finite lattice spacing. We show that the lighter neutral pion in tmLQCD gives rise to finite volume effects that are…
This study provides an abstract framework to analyze mixed formulations in viscoelasticity, in the classic saddle point form. Standard hypothesis for mixed methods are adapted to the Volterra type equations in order to obtain stability of…
We propose an alternative method to generate samples of a spatially correlated random field with applications to large-scale problems for forward propagation of uncertainty. A classical approach for generating these samples is the…
The truncated Lanczos method using a variational scheme based on Hilbert space reduction as well as a local basis change is re-examined. The energy is extrapolated as a power law function of the Hamiltonian variance. This systematic…
We investigate the mass differences of charm and bottom heavy baryons, using MILC lattice gauge configurations with 2+1 dynamical sea quark flavors. We extend our previous work to three lattice spacings: fine (a~0.09), coarse (a~0.12), and…
In this paper, a systematic approach of constructing modified equations for weak stochastic symplectic methods of stochastic Hamiltonian systems is given via using the generating functions of the stochastic symplectic methods. This approach…
This paper proposes a novel way to solve transient linear, and non-linear solid dynamics for compressible, nearly incompressible, and incompressible material in the updated Lagrangian framework for tetrahedral unstructured finite elements.…
We determine the light quark masses from lattice QCD simulations incorporating the electromagnetic interaction of valence quarks. The meson masses are calculated on lattice QCD configurations generated by the RBC Collaboration for two…
Simulations that use the clover action in quenched QCD calculations have a lower limit to the quark mass that can be reached, because of the fluctuations caused by exceptional configurations. From this low statistics study, we find that the…