We present a novel general Boltzmann machine with continuous visible

and discrete integer valued hidden states, yielding a parametric

density function involving a ratio of Riemann-Theta functions. After a

brief overview of the theory required to define this new ML

architecture, we show how the conditional expectation of a hidden

state for given visible states can be used as activation function...

We report on multi-loop integral computations executed on a

PEZY/Exascaler large-scale (immersion cooling) computing system.

The programming model requires a host program written in C++

with an OpenCL kernel. However the kernel can be generated by

the Goose compiler interface, which allows parallelizing loops

according to compiler directives. As an advantage, the executable

derived from a...

The high-energy community recently witnessed the first attempts at leveraging machine (deep) learning techniques for improving the efficiency of the numerical Monte-Carlo integrations that lie at the core of most high-energy physics simulations.

The first part of my talk will characterise the various type of integrations necessary in these simulations as well as the type of improvements that...

Massively parallel simulations generate increasing volumes of large data, whose exploitation requires large storage resources, efficient network and increasingly large post-processing facilities. In the coming era of exascale computations, there is an emerging need for new data analysis and visualization strategies.

Data manipulation, during the simulation and after, considerably slows down...

A key ingredient in an automated evaluation of two-loop multileg processes is a

fast and numerically stable evaluation of scalar Feynman integrals. In this respect, the calculation of two-loop three- and four-point functions in the general complex mass case so far relies on multidimensional numerical integration through sector decomposition whereby a reliable result has a high computing cost,...

Complete one-loop electroweak radiative corrections to polarized Bhabha

scattering are presented. Higher order QED effects are evaluated in the leading

logarithmic approximation. Numerical results are shown for the conditions of future

circular and linear electron-positron colliders with polarized beams. Theoretical

uncertainties are estimated.

A new Monte Carlo event generator MCSANCee for simulation of processes at future e^+e^- colliders is presented. Complete one-loop electroweak radiative corrections and polarization of the initial beams are taken into account. The present generator includes the following processes: e^+e^- \to e^+e^- (mu^+mu^-, tau^+tau^-, ZH, Z\gamma, \gamma\gamma). Numerical results for all of these processes...

The Grace system is an automatic system to calculate cross sections based on the standard model and MSSM including one-loop corrections.

I would like to report recent progress of the GRACE system including optimization of generated codes.

The Monte Carlo generator to simulate events of single-photon

annihilation to hadrons at center-of-mass energies below 2.5 GeV

is described. The generator is based on existing data on cross sections

of various exclusive channels of e+e- annihilation obtained in various

e+e- experiments by the scan and ISR methods. It is extensively used

in the software packages for analysis of experiments at...

A high-precision calculation of the electron anomalous magnetic moment requires an evaluation of QED Feynman diagrams up to five independent loops. To make this calculation practically feasible it is necessary to remove all infrared and ultraviolet divergences before integration. A procedure of removing both infrared and ultraviolet divergences in each individual Feynman diagram will be...

While the Higgs boson couplings to other particles are increasingly well-measured by LHC experiments, it has proven difficult to set constraints on the Higgs trilinear self-coupling $\lambda$, principally due to the very low cross-section of Higgs boson pair production. We present the results of NLO QCD corrections to Higgs pair production with full top-quark mass dependence, where the...

In this talk, we consider some of the computational aspects encountered in recent computations of double Higgs boson production in gluon fusion. We consider the NLO virtual amplitude in the high-energy limit, and the NNLO virtual amplitude in the low-energy (or large top quark mass) limit. We discuss various optimizations which were necessary to produce our results.

We present an algorithm which allows to solve analytically linear systems of differential equations which factorize to first order. The solution is given in terms of iterated integrals over an alphabet where its structure is implied by the coefficient matrix of the differential equations. These systems appear in a large variety of higher order calculations in perturbative Quantum Field...

In this contribution I will discuss the practicalities of storing events from a NNLO calculation on disk with the view of "replaying" the simulation for a different analysis and under different conditions, such as a different PDF fit or a different scale setting.

We present the HepMC3 library designed to perform manipulations with

event records of High Energy Physics Monte Carlo Event Generators

(MCEGs). The library is a natural successor of HepMC and HepMC2

libraries used in the present and in the past. HepMC3 supports all

functionality of previous versions and significantly extends them.

In comparison to the previous versions, the default event...

I briefly review the recently finished 5-loop renormalization program of QCD, and explain the status and prospects of the computer-algebraic techniques involved.

We propose an algorithm to find a solution to an integro-differential equation of the DGLAP type for all the orders in the running coupling α with splitting functions given at a fixed order in α. Complex analysis is significantly used in the construction of the algorithm, we found a way to calculate the involved integrals over contours in the complex planes in more simple way than by any of...

I give an update on recent developments in FeynArts, FormCalc, and LoopTools, and show how the new features were used in making the latest version of FeynHiggs.

The software framework SModelS, which has already been presented at the ACAT 2016 conference, allows for a very fast confrontation of arbitrary BSM models exhibiting a Z2 symmetry with an ever growing database of simplified models results from CMS and ATLAS. In this talk we shall present its newest features, like the extension to include searches for heavy stable charged particles (HSCPs), or...

FORM is a symbolic manipulation system, which is especially advantageous for handling gigantic expressions with many small terms. Because FORM has been developed in tackling real problems in perturbative quantum field theory, it has some features useful in such problems, although FORM applications are not restricted to any specific research field. In this talk, we discuss recent developments...

The talk is devoted to the overview of *Rings* — an efficient lightweight library for commutative algebra written in Java and Scala languages. Polynomial arithmetic, GCDs, polynomial factorization and Gröbner bases are implemented with the use of modern asymptotically fast algorithms. *Rings* can be easily interacted or embedded in applications in high-energy physics and other research areas...

Over the last few years manipulating expressions with millions of terms has become common in particle physics. Form is the de facto tool for manipulations of extremely large expressions, but it comes with some downsides. In this talk I will discuss an effort to modernize aspects of Form, such as the language and workflow, and the introduction of bindings to C and Python. This new tool is...

The Mathematica package STR (Star-Triangle Relations) is a recently developed tool designed to solve Feynman diagrams by means of the method of uniqueness in any (Euclidean) spacetime dimension D. The method of uniqueness is a powerful technique to solve multi-loop Feynman integrals in theories with conformal symmetry imposing some relations between D and the powers of propagators. In our...

Tensor calculations are an important case in many natural sciences like mathematics and physics. To simplify such expressions, computer algebra is widely used. There are a number of approaches for solving this problem, namely, the component calculations, the calculations when tensor is considered as an abstract symbol with indices possessing some symmetry properties, and finally a pure...

Decoding the nature of Dark Matter (DM) is one of the most important problems of particle physics. DM can potentially provide unique signatures at collider and non-collider experiments. Details of these signatures which we expect to observe in the near future would allow us to delineate the properties of DM and the respective underlying theory Beyond the Standard Model (BSM). While there...

The recent years have shown an exciting development in the scientific commmunity due to the interplay between new methods from data science and artificial intelligence, increasing computational resources and physics. The fundamental object of our theories of nature is the Lagrangian whose form is determined by the symmetries found already. A famous and well-motivated extension of the SM...

The set of the four-loop numerical results for the relation between pole and running heavy quarks masses in QCD at fixed number of lighter flavors $3\leq n_l\leq 15$, which was obtained in Ref.[1] with help of the Lomonosov Supercomputer of Moscow State University, is analysed by the ordinary method of the least squares. We use its variant which allows to solve the overdetermined system of 13...

Amplitude analysis is an important tool for the research of the hadron spectrum, in which the maximum likelihood method is used to estimate the parameters of a probability density function. In each optimization step, the likelihood values of a huge number of events from both data and Monte-Carlo simulations are calculated and summed, which is the most time-consuming part of the whole...

micrOMEGAs is a package for the calculation of the relic density of Dark Matter and of different observables related with Dark Matter searches. The talk will present the general structure of the package and several recent developments including freeze-in relic abundance calculation, interface with different packages that compute collider observables, and recent improvements in direct detection signals.

An efficient phase space integration is important for most calculations for collider experiments. We are developing a phase space integration that distribute phase space points according to the singular limit of QCD, Using the Altarelli-Parisi splitting functions as the underlying probability for a splitting, by developing and applying theoretical and computational tools

We describe a little in detail some techniques used in the calculation of

4-loop QED contributions to some quantities like g-2,

slope of the Dirac form factor, renormalization constants;

in particular, some different approaches to the parallelization of some

parts of the calculations. Some recent results will be also presented.