Markov Chain Mote Carlo solution of BK equation through Newton-Kantorovich method

Abstract

We propose a new method for Monte Carlo solution of non-linear integral equations by combining the Newton-Kantorovich method for solving non-linear equations with the Markov Chain Monte Carlo (MCMC) method for solving linear equations. The Newton-Kantorovich method allows to express the non-linear equation as a system of the linear equations which then can be treated by the MCMC (random walk) algorithm. We apply this method to the Balitsky-Kovchegov (BK) equation describing evolution of gluon density at low x. Results of numerical computations show that the MCMC method is both precise and efficient. The presented algorithm may be particularly suited for solving more complicated and higher-dimensional non-linear integral equation, for which traditional methods become unfeasible.

10 pages, 2 figures, typos corrected, accepted by JHEP

References in corpus (9)

• PYTHIA 6.4 Physics and Manual
• Herwig++ Physics and Manual
• The ridge in proton-proton collisions at the LHC
• The CCFM Monte Carlo generator CASCADE 2.2.0
• Azimuthal correlations of forward di-hadrons in d+Au collisions at RHIC in the Color Glass Condensate
• Evidence for BFKL and saturation dynamics from di-hadron spectra at the LHC
• Gluon saturation in dijet production in p-Pb collisions at Large Hadron Collider
• Nonlinear equation for coherent gluon emission
• Resummation in nonlinear equation for high energy factorizable gluon density and its extension to include coherence