Source author record

Alberto Ramos

Alberto Ramos appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2022arXiv

Toward the End-to-End Optimization of Particle Physics Instruments with Differentiable Programming: a White Paper

The full optimization of the design and operation of instruments whose functioning relies on the interaction of radiation with matter is a super-human task, given the large dimensionality of the space of possible choices for geometry, detection technology, materials, data-acquisition, and information-extraction techniques, and the interdependence of the related parameters. On the other hand, massive potential gains in performance over standard, "experience-driven" layouts are in principle within our reach if an objective function fully aligned with the final goals of the instrument is maximized by means of a systematic search of the configuration space. The stochastic nature of the involved quantum processes make the modeling of these systems an intractable problem from a classical statistics point of view, yet the construction of a fully differentiable pipeline and the use of deep learning techniques may allow the simultaneous optimization of all design parameters. In this document we lay down our plans for the design of a modular and versatile modeling tool for the end-to-end optimization of complex instruments for particle physics experiments as well as industrial and medical applications that share the detection of radiation as their basic ingredient. We consider a selected set of use cases to highlight the specific needs of different applications.

preprint2020arXiv

An analysis of systematic effects in finite size scaling studies using the gradient flow

We propose a new strategy for the determination of the step scaling function $σ(u)$ in finite size scaling studies using the Gradient Flow. In this approach the determination of $σ(u)$ is broken in two pieces: a change of the flow time at fixed physical size, and a change of the size of the system at fixed flow time. Using both perturbative arguments and a set of simulations in the pure gauge theory we show that this approach leads to a better control over the continuum extrapolations. Following this new proposal we determine the running coupling at high energies in the pure gauge theory and re-examine the determination of the $Λ$-parameter, with special care on the perturbative truncation uncertainties.

preprint2020arXiv

Light quark masses in N_f = 2+1 lattice QCD with Wilson fermions

We present a lattice QCD determination of light quark masses with three sea-quark flavours ($N_f = 2+1$). Bare quark masses are known from PCAC relations in the framework of CLS lattice computations with a non-perturbatively improved Wilson-Clover action and a tree-level Symanzik improved gauge action. They are fully non-perturbatively improved, including the recently computed Symanzik counter-term $b_{\rm A} - b_{\rm P}$. The mass renormalisation at hadronic scales and the renormalisation group running over a wide range of scales are known non-perturbatively in the Schrödinger functional scheme. In the present paper we perform detailed extrapolations to the physical point, obtaining (for the four-flavour theory) $m_{u/d}(2{\rm GeV}) = 3.54(12)(9)$ MeV and $m_s(2{\rm GeV}) = 95.7(2.5)(2.4)$ MeV in the $\bar{MS}$ scheme. For the mass ratio we have $m_s/m_{u/d} = 27.0(1.0)(0.4)$. The RGI values in the three-flavour theory are $M_{u/d} = 4.70(15)(12)$ MeV and $M_s = 127.0(3.1)(3.2)$ MeV.

preprint2020arXiv

The twisted gradient flow running coupling in SU(3): a non-perturbative determination

We report some preliminary results of our ongoing non-perturbative computation of the twisted 't Hooft running coupling in a particular set-up, using the gradient flow to define the coupling and step scaling techniques to compute it. For the computation we considered a pure gauge SU(3) theory in four dimensions, defined on the lattice on an asymmetrical torus endowed with twisted boundary conditions in a single plane, and related the energy scale of the coupling to an effective size combining the size of the torus and the rank of the gauge group. Additionally, we explore some of the effects of the freezing of the topology on the computation of the coupling.

preprint2019arXiv

Non-perturbative renormalization by decoupling

We propose a new strategy for the determination of the QCD coupling. It relies on a coupling computed in QCD with $N_{\rm f} \geq 3$ degenerate heavy quarks at a low energy scale $μ_{\rm dec}$, together with a non-perturbative determination of the ratio $Λ/μ_{\rm dec}$ in the pure gauge theory. We explore this idea using a finite volume renormalization scheme for the case of $N_{\rm f} = 3$ QCD, demonstrating that a precise value of the strong coupling $α_s$ can be obtained. The idea is quite general and can be applied to solve other renormalization problems, using finite or infinite volume intermediate renormalization schemes.

preprint2019arXiv

openQ*D code: a versatile tool for QCD+QED simulations

We present the open-source package openQ*D-1.0, which has been primarily, but not uniquely, designed to perform lattice simulations of QCD+QED and QCD, with and without C* boundary conditions, and O(a) improved Wilson fermions. The use of C* boundary conditions in the spatial direction allows for a local and gauge-invariant formulation of QCD+QED in finite volume, and provides a theoretically clean setup to calculate isospin-breaking and radiative corrections to hadronic observables from first principles. The openQ*D code is based on openQCD-1.6 and NSPT-1.4. In particular it inherits from openQCD-1.6 several core features, e.g. the highly optimized Dirac operator, the locally deflated solver, the frequency splitting for the RHMC, or the 4th order OMF integrator.

preprint2016arXiv

Determination of the QCD $Λ$-parameter and the accuracy of perturbation theory at high energies

We discuss the determination of the strong coupling $α_\mathrm{\overline{MS}}^{}(m_\mathrm{Z})$ or equivalently the QCD $Λ$-parameter. Its determination requires the use of perturbation theory in $α_s(μ)$ in some scheme, $s$, and at some energy scale $μ$. The higher the scale $μ$ the more accurate perturbation theory becomes, owing to asymptotic freedom. As one step in our computation of the $Λ$-parameter in three-flavor QCD, we perform lattice computations in a scheme which allows us to non-perturbatively reach very high energies, corresponding to $α_s = 0.1$ and below. We find that (continuum) perturbation theory is very accurate there, yielding a three percent error in the $Λ$-parameter, while data around $α_s \approx 0.2$ is clearly insufficient to quote such a precision. It is important to realize that these findings are expected to be generic, as our scheme has advantageous properties regarding the applicability of perturbation theory.

preprint2016arXiv

Non-perturbative running of quark masses in three-flavour QCD

We present our preliminary results for the computation of the non-perturbative running of renormalized quark masses in $N_f = 3$ QCD, between the electroweak and hadronic scales, using standard finite-size scaling techniques. The computation is carried out to very high precision, using massless $\mathcal{O}(a)$-improved Wilson quarks. Following the strategy adopted by the ALPHA Collaboration for the running coupling, different schemes are used above and below a scale $μ_0 \sim m_b$, which differ by using either the Schrödinger Functional or Gradient Flow renormalized coupling. We discuss our results for the running in both regions, and the procedure to match the two schemes.

preprint2016arXiv

The determination of $α_s$ by the ALPHA collaboration

We review the ALPHA collaboration strategy for obtaining the QCD coupling at high scale. In the three-flavor effective theory it avoids the use of perturbation theory at $α> 0.2$ and at the same time has the physical scales small compared to the cutoff $1/a$ in all stages of the computation. The result $Λ_\overline{MS}^{(3)}=332(14)$~MeV is translated to $α_\overline{MS}(m_Z)=0.1179(10)(2)$ by use of (high order) perturbative relations between the effective theory couplings at the charm and beauty quark "thresholds". The error of this perturbative step is discussed and estimated as $0.0002$.

preprint2015arXiv

(Dimensional) twisted reduction in large N gauge theories

We show that the spontaneous breaking of center symmetry can be avoided on a $L^2\times 1^2$ lattice with the appropriate choice of twisted boundary conditions. In order for this to work it is crucial that the twisted boundary conditions are chosen in the reduced plane. This suggests that the choice of twist tensor can influence the directions in which color and space degrees of freedom become indistinguishable. We also present some preliminary quantitative data comparing the value of the plaquette for different forms of reduction.

preprint2015arXiv

A comparison of updating algorithms for large N reduced models

We investigate Monte Carlo updating algorithms for simulating $SU(N)$ Yang-Mills fields on a single-site lattice, such as for the Twisted Eguchi-Kawai model (TEK). We show that performing only over-relaxation (OR) updates of the gauge links is a valid simulation algorithm for the Fabricius and Haan formulation of this model, and that this decorrelates observables faster than using heat-bath updates. We consider two different methods of implementing the OR update: either updating the whole $SU(N)$ matrix at once, or iterating through $SU(2)$ subgroups of the $SU(N)$ matrix, we find the same critical exponent in both cases, and only a slight difference between the two.

preprint2015arXiv

A status update on the determination of $Λ_{\overline{\rm MS}}^{N_{\rm f}=3}$ by the ALPHA collaboration

The ALPHA collaboration aims to determine $α_s(m_Z)$ with a total error below the percent level. A further step towards this goal can be taken by combining results from the recent simulations of 2+1-flavour QCD by the CLS initiative with a number of tools developed over the years: renormalized couplings in finite volume schemes, recursive finite size techniques, two-loop renormalized perturbation theory and the (improved) gradient flow on the lattice. We sketch the strategy, which involves both the standard SF coupling in the high energy regime and a gradient flow coupling at low energies. This implies the need for matching both schemes at an intermediate switching scale, $L_{\rm swi}$, which we choose roughly in the range 2-4 GeV. In this contribution we present a preliminary result for this matching procedure, and we then focus on our almost final results for the scale evolution of the SF coupling from $L_{\rm swi}$ towards the perturbative regime, where we extract the $N_{\rm f} = 3$ $Λ$-parameter, $Λ_{\overline{\rm MS}}^{N_{\rm f}=3}$, in units of $L_{\rm swi}$ . Connecting $L_{\rm swi}$ and thus the $Λ$-parameter to a hadronic scale such as $F_K$ requires 2 further ingredients: first, the connection of $L_{\rm swi}$ to $L_{\rm max}$ using a few steps with the step-scaling function of the gradient flow coupling, and, second, the continuum extrapolation of $L_{\rm max} F_K$.

preprint2015arXiv

Charged hadrons in local finite-volume QED+QCD with C* boundary conditions

In order to calculate QED corrections to hadronic physical quantities by means of lattice simulations, a coherent description of electrically-charged states in finite volume is needed. In the usual periodic setup, Gauss's law and large gauge transformations forbid the propagation of electrically-charged states. A possible solution to this problem, which does not violate the axioms of local quantum field theory, has been proposed by Wiese and Polley, and is based on the use of C* boundary conditions. We present a thorough analysis of the properties and symmetries of QED in isolation and QED coupled to QCD, with C* boundary conditions. In particular we learn that a certain class of electrically-charged states can be constructed in this setup in a fully consistent fashion, without relying on gauge fixing. We argue that this class of states covers most of the interesting phenomenological applications in the framework of numerical simulations. We also calculate finite-volume corrections to the mass of stable charged particles and show that these are much smaller than in non-local formulations of QED.

preprint2015arXiv

Lattice study for conformal windows of SU(2) and SU(3) gauge theories with fundamental fermions

We present our investigation of SU(2) gauge theory with 8 flavours, and SU(3) gauge theory with 12 flavours. For the SU(2) case, at strong bare coupling, $β\lesssim 1.45$, the distribution of the lowest eigenvalue of the Dirac operator can be described by chiral random matrix theory for the Gaussian symplectic ensemble. Our preliminary result indicates that the chiral phase transition in this theory is of bulk nature. For the SU(3) theory, we use high-precision lattice data to perform the step-scaling study of the coupling, $g_{\rm GF}$, in the Gradient Flow scheme. We carefully examine the reliability of the continuum extrapolation in the analysis, and conclude that the scaling behaviour of this SU(3) theory is not governed by possible infrared conformality at $g_{\rm GF}^{2} \lesssim 6$.

preprint2015arXiv

On O($a^2$) effects in gradient flow observables

In lattice gauge theories, the gradient flow has been used extensively both, for scale setting and for defining finite volume renormalization schemes for the gauge coupling. Unfortunately, rather large cutoff effects have been observed in some cases. We here investigate these effects to leading order in perturbation theory, considering various definitions of the lattice observable, the lattice flow equation and the Yang Mills lattice action. These considerations suggest an improved set- up for which we perform a scaling test in the pure SU(3) gauge theory, demonstrating strongly reduced cutoff effects. We then attempt to obtain a more complete understanding of the structure of O($a^2$) effects by applying Symanzik's effective theory approach to the 4+1 dimensional local field theory with flow time as the fifth dimension. From these considerations we are led to a fully O($a^2$) improved set-up the study of which is left to future work.

preprint2015arXiv

The Yang-Mills gradient flow and SU(3) gauge theory with 12 massless fundamental fermions in a colour-twisted box

We perform the step-scaling investigation of the running coupling constant, using the gradient-flow scheme, in SU(3) gauge theory with twelve massless fermions in the fundamental representation. The Wilson plaquette gauge action and massless unimproved staggered fermions are used in the simulations. Our lattice data are prepared at high accuracy, such that the statistical error for the renormalised coupling, g_GF, is at the subpercentage level. To investigate the reliability of the continuum extrapolation, we employ two different lattice discretisations to obtain g_GF. For our simulation setting, the corresponding gauge-field averaging radius in the gradient flow has to be almost half of the lattice size, in order to have this extrapolation under control. We can determine the renormalisation group evolution of the coupling up to g^2_GF ~ 6, before the onset of the bulk phase structure. In this infrared regime, the running of the coupling is significantly slower than the two-loop perturbative prediction, although we cannot draw definite conclusion regarding possible infrared conformality of this theory. Furthermore, we comment on the issue regarding the continuum extrapolation near an infrared fixed point. In addition to adopting the fit ansatz a'la Symanzik for performing this task, we discuss a possible alternative procedure inspired by properties derived from low-energy scale invariance at strong coupling. Based on this procedure, we propose a finite-size scaling method for the renormalised coupling as a means to search for infrared fixed point. Using this method, it can be shown that the behaviour of the theory around g^2_GF ~ 6 is still not governed by possible infrared conformality.

preprint2014arXiv

Form factors for $\mathrm B_\mathrm s \to \mathrm K \ell ν$ decays in Lattice QCD

We present the current status of the computation of the form factor $f_+ (q^2)$ for the semi-leptonic decay $\mathrm B_\mathrm s \to \mathrm K \ell ν$ by the ALPHA collaboration. We use gauge configurations which were generated as part of the Coordinated Lattice Simulations (CLS) effort. They have $N_\mathrm f=2$ non-perturbatively $O(a)$ improved Wilson fermions, and pion masses down to $\approx 250 \,\mathrm {MeV}$ with $m_πL \geq 4$. The heavy quark is treated in non-perturbative Heavy Quark Effective Theory (HQET). We discuss how to extract the form factors from the correlation functions and present first results for the form factor at $q^2 = 21.23\,\mathrm{GeV}^2$ extrapolated to the continuum. Next-to-leading order terms in HQET and the chiral extrapolation still need to be included in the analysis.

preprint2014arXiv

SU(3) gauge theory with 12 flavours in a twisted box

We present preliminary result for the step-scaling study of the coupling constant with the Yang-Mills gradient flow, in the twelve-favour SU(3) gauge theory. In this work, the lattice simulation is performed using unimproved staggered fermions and the Wilson plaquette gauge action, from which the gradient flow is also implemented. Imposing twisted boundary condition a'la t'Hooft and Parisi, our calculation is performed at zero fermion mass. The renormalised coupling constant is extracted via the computation of the energy density. In order to examine the reliability of the continuum extrapolation, we investigate this coupling constant using two different lattice discretisations. Our result shows that in order to control the systematic effects in the continuum extrapolation, it is necessary to implement a large enough gradient-flow time. In the current calculation, the gauge-field averaging radius corresponding to the flow time has to be as large as 40% of the lattice size.

preprint2012arXiv

B $\to$ $π$ form factor with 2 flavours of $O(a)$ improved Wilson quarks

The determinations of $|V_{\rm ub}|$ from the exclusive branching ratios of $B\to τν$ and $B \to πl ν$ tend to show a tension at the level of $3σ$ \cite{Beringer:1900zz}. On the theoretical side they depend on the lattice computation of the hadronic matrix elements $f_{\rm B}$ and the $B\to π$ form factor $f_+(q^2)$. To understand the tension, improved precision and a careful analysis of the systematics involved are necessary. Working towards this goal, we present preliminary lattice results of the ALPHA collaboration for the $B\to π$ form factor $f_+(q^2)$ with $N_{\rm f}=2$ flavours of $O(a)$-improved Wilson fermions. Our computation uses HQET in the static limit, pion masses ranging down to $\sim250$ MeV, large volumes with $m_πL >4$, three lattice spacings, and non-perturbative renormalization. We describe the techniques adopted to reduce the statistical noise (stochastic all-to-all with full time dilution) and the contamination from excited states (smearing for the B and the pion). We estimate the size of the chiral and continuum extrapolations. We discuss the impact our result could have to clarify the above mentioned discrepancy in the determination of $|V_{\rm ub}|$.

preprint2012arXiv

Playing with the kinetic term in the HMC

The HMC algorithm, combining the advantages of molecular dynamics and Monte-Carlo methods, is the most efficient algorithm to simulate QCD including the effects of sea quarks. In the standard approach momentum fields are generated with a Gaussian probability density. In this work I will explore another possibility. By using a Lorentz distribution one can dynamically impose a cutoff in the rate of change of the coordinates that potentially could have a better behavior. I will present some results in pure SU(2) gauge theory.

preprint2012arXiv

Systematic errors in partially-quenched QCD plus QED lattice simulations

At the precision reached in current lattice QCD calculations, electromagnetic effects are becoming numerically relevant. Here, electromagnetic effects are included by superimposing $\mathrm{U}(1)$ degrees of freedom on $N_f = 2+1$ QCD configurations from the Budapest-Marseille-Wuppertal Collaboration. We present preliminary results for the electromagnetic corrections to light pseudoscalars mesons masses and discuss some of the associated systematic errors.

preprint2006arXiv

Dynamics of critical vortices on the torus and on the plane

We use the Bradlow parameter expansion to construct the metric tensor in the space of solutions of the Bogomolny equations for the Abelian Higgs model on a two-dimensional torus. Using this metric we study the dynamics and scattering of vortices on the torus within the geodesic approximation. For small torus volumes the metric is determined in terms of a small number of parameters. For large volumes the results provide a very precise approximation to the metric and dynamics on the plane.