ترغب بنشر مسار تعليمي؟ اضغط هنا

Asymptotics of Fixed Point Distributions for Inexact Monte Carlo Algorithms

30   0   0.0 ( 0 )
 نشر من قبل Michael A. Clark
 تاريخ النشر 2007
  مجال البحث
والبحث باللغة English




اسأل ChatGPT حول البحث

We introduce a simple general method for finding the equilibrium distribution for a class of widely used inexact Markov Chain Monte Carlo algorithms. The explicit error due to the non-commutivity of the updating operators when numerically integrating Hamiltons equations can be derived using the Baker-Campbell-Hausdorff formula. This error is manifest in the conservation of a ``shadow Hamiltonian that lies close to the desired Hamiltonian. The fixed point distribution of inexact Hybrid algorithms may then be derived taking into account that the fixed point of the momentum heatbath and that of the molecular dynamics do not coincide exactly. We perform this derivation for various inexact algorithms used for lattice QCD calculations.

قيم البحث

اقرأ أيضاً

99 - Alberto Ramos 2018
Automatic Differentiation (AD) allows to determine exactly the Taylor series of any function truncated at any order. Here we propose to use AD techniques for Monte Carlo data analysis. We discuss how to estimate errors of a general function of measur ed observables in different Monte Carlo simulations. Our proposal combines the $Gamma$-method with Automatic differentiation, allowing exact error propagation in arbitrary observables, even those defined via iterative algorithms. The case of special interest where we estimate the error in fit parameters is discussed in detail. We also present a freely available fortran reference implementation of the ideas discussed in this work.
We present a comprehensive new global QCD analysis of polarized inclusive deep-inelastic scattering, including the latest high-precision data on longitudinal and transverse polarization asymmetries from Jefferson Lab and elsewhere. The analysis is pe rformed using a new iterative Monte Carlo fitting technique which generates stable fits to polarized parton distribution functions (PDFs) with statistically rigorous uncertainties. Inclusion of the Jefferson Lab data leads to a reduction in the PDF errors for the valence and sea quarks, as well as in the gluon polarization uncertainty at $x gtrsim 0.1$. The study also provides the first determination of the flavor-separated twist-3 PDFs and the $d_2$ moment of the nucleon within a global PDF analysis.
We explore connections between two common methods for quantifying the uncertainty in parton distribution functions (PDFs), based on the Hessian error matrix and Monte-Carlo sampling. CT14 parton distributions in the Hessian representation are convert ed into Monte-Carlo replicas by a numerical method that reproduces important properties of CT14 Hessian PDFs: the asymmetry of CT14 uncertainties and positivity of individual parton distributions. The ensembles of CT14 Monte-Carlo replicas constructed this way at NNLO and NLO are suitable for various collider applications, such as cross section reweighting. Master formulas for computation of asymmetric standard deviations in the Monte-Carlo representation are derived. A correction is proposed to address a bias in asymmetric uncertainties introduced by the Taylor series approximation. A numerical program is made available for conversion of Hessian PDFs into Monte-Carlo replicas according to normal, log-normal, and Watt-Thorne sampling procedures.
Supersymmetric Yang-Mills (SYM) theories in four dimensions exhibit many interesting non-perturbative phenomena that can be studied by means of Monte Carlo lattice simulations. However, the lattice regularization breaks supersymmetry explicitly, and in general a fine tuning of a large number of parameters is required to correctly extrapolate the theory to the continuum limit. From this perspective, it is important to preserve on the lattice as many symmetries of the original continuum action as possible. Chiral symmetry for instance prevents an additive renormalization of the fermion mass. A (modified) version of chiral symmetry can be preserved exactly if the Dirac operator fulfills the Ginsparg-Wilson relation. In this contribution, we present an exploratory non-perturbative study of N=1 supersymmetric Yang-Mills theory using the overlap formalism to preserve chiral symmetry at non-zero lattice spacings. N=1 SYM is an ideal benchmark toward the extension of our studies to more complex supersymmetric theories, as the only parameter to be tuned is the gluino mass. Overlap fermions allow therefore to simulate the theory without fine-tuning. We compare our approach to previous investigations of the same theory, and we present clear evidences for gluino condensation.
Hasenbusch has proposed splitting the pseudo-fermionic action into two parts, in order to speed-up Hybrid Monte Carlo simulations of QCD. We have tested a different splitting, also using clover-improved Wilson fermions. An additional speed-up between 5 and 20% over the original proposal was achieved in production runs.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا