Related papers: Calculation of kinetic parameters $\beta_{\mathit{…
Recent calculations of nuclear structure corrections to the Lamb shift in light muonic atoms are based on an expansion in a parameter eta, where only terms up to second order are retained. The parameter eta can be shown to be proportional…
In this paper we present a new approach to control variates for improving computational efficiency of Ensemble Monte Carlo. We present the approach using simulation of paths of a time-dependent nonlinear stochastic equation. The core idea…
We discuss the current development of MC/DC (Monte Carlo Dynamic Code). MC/DC is primarily designed to serve as an exploratory Python-based MC transport code. However, it seeks to offer improved performance, massive scalability, and backend…
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…
This MATLAB program calculates the dynamics of the reduced density matrix of an open quantum system modeled by the Feynman-Vernon model. The user gives the program a vector describing the coordinate of an open quantum system, a hamiltonian…
Neutrinoless double-beta decay ($0\nu\beta\beta$) is a rare hypothesised process that, if discovered, would establish that the neutrino is Majorana, that is, it is its own antiparticle. Interpretation of experimental results relies on…
We present Path Integral Monte Carlo C code for calculation of quantum mechanical transition amplitudes for 1D models. The SPEEDUP C code is based on the use of higher-order short-time effective actions and implemented to the maximal order…
This paper presents the benchmarking and scaling studies of a GPU accelerated three dimensional compressible magnetohydrodynamic code. The code is developed keeping an eye to explain the large and intermediate scale magnetic field…
In this work, we propose a new flow-matching Markov chain Monte Carlo (FM-MCMC) algorithm for estimating the orbital parameters of exoplanetary systems, especially for those only one exoplanet is involved. Compared to traditional methods…
In this paper, we evaluate the performance of the multilevel Monte Carlo method (MLMC) for deterministic and uncertain hyperbolic systems, where randomness is introduced either in the modeling parameters or in the approximation algorithms.…
We compare our Monte Carlo reaction rates (see Paper II of this series) to previous results that were obtained by using the classical method of computing thermonuclear reaction rates. For each reaction, the comparison is presented using two…
The recently developed auxiliary field diffusion Monte Carlo method is applied to compute the equation of state and the compressibility of neutron matter. By combining diffusion Monte Carlo for the spatial degrees of freedom and auxiliary…
In this paper, we put constraints on neutrino properties such as mass $m_{\nu}$ and degeneracy parameters $\xi_i$ from WMAP5 data and light element abundances by using a Markov chain Monte Carlo (MCMC) approach. In order to take…
Determinant Quantum Monte Carlo (DQMC) provides numerically exact solutions for strongly correlated fermionic systems but faces significant computational challenges with increasing system size. While submatrix updates were originally…
The one-loop contribution of scalar and vector leptoquarks (LQs) to the electromagnetic properties (NEPs) of massive Dirac neutrinos is presented via an effective Lagrangian approach, with emphasis on the effective neutrino charge radius…
We present a deep neural network (DNN) accelerated Hamiltonian Monte Carlo (HMC) algorithm called DeepHMC for the inference of binary neutron star systems. The HMC is a non-random walk sampler that uses background gradient information to…
To better understand the capture process by a nanopore, we introduce an efficient Kinetic Monte Carlo (KMC) algorithm that can simulate long times and large system sizes by mapping the dynamic of a point-like particle in a 3D spherically…
This work covers methodology of solving QCD evolution equation of the parton distribution using Markovian Monte Carlo (MMC) algorithms in a class of models ranging from DGLAP to CCFM. One of the purposes of the above MMCs is to test the…
In many applications, such as plasma edge simulation of a nuclear fusion reactor, a coupled PDE/kinetic description is required, which is usually solved with a coupled finite-volume/Monte-Carlo method. Different procedures have been…
Energy integrated neutrino cross sections are integrals of nuclear responses weighted with kinematic prefactors. We decompose the prefactors into a limited set of functions of energy transfer and show the relevant integrals are the moments…