-
Electrical conductivity of the Quark-Gluon Plasma in the presence of strong magnetic fields
Authors:
Giorgio Almirante,
Nikita Astrakhantsev,
V. V. Braguta,
Massimo D'Elia,
Lorenzo Maio,
Manuel Naviglio,
Francesco Sanfilippo,
Anton Trunin
Abstract:
We compute the electrical conductivity of the strongly interacting medium in the presence of strong magnetic background fields, $eB=4,9~GeV^2$, and for different values of the temperature, both in the confined and in the deconfined Quark-Gluon Plasma (QGP) phase. The conductivity is obtained from the Euclidean lattice time correlator of the electrical current, computed on gauge configurations samp…
▽ More
We compute the electrical conductivity of the strongly interacting medium in the presence of strong magnetic background fields, $eB=4,9~GeV^2$, and for different values of the temperature, both in the confined and in the deconfined Quark-Gluon Plasma (QGP) phase. The conductivity is obtained from the Euclidean lattice time correlator of the electrical current, computed on gauge configurations sampled from Monte-Carlo simulations of an improved staggered discretization of $N_f = 2+1$ QCD. We perform the inverse Laplace transform of the correlator adopting a recently-proposed version of the standard Backus--Gilbert procedure for the inversion. The results obtained in the QGP phase show a sizable enhancement of the conductivity in the direction parallel to the magnetic field, as well as a suppression in the direction orthogonal to it. Such enhancement could be attributed to the manifestation of the Chiral Magnetic Effect (CME): following this guess, we extract the behaviour of the relaxation time of this process, extrapolate it to the continuum limit and compare it to previous results, finding it lower than expected in the explored range of temperatures.
△ Less
Submitted 26 June, 2024;
originally announced June 2024.
-
Non-perturbative determination of the $N_f=2+1$ QCD sphaleron rate
Authors:
Claudio Bonanno,
Francesco D'Angelo,
Massimo D'Elia,
Manuel Naviglio,
Lorenzo Maio
Abstract:
The strong sphaleron rate, i.e., the rate of real time QCD topological transitions, is a key phenomenological quantity, playing a fundamental role in several physical contexts. In heavy-ion collisions, a non-vanishing rate can lead to the so-called Chiral Magnetic Effect. In early-Universe cosmology, instead, it can be related to the rate of thermal production of QCD axions. In this talk, we prese…
▽ More
The strong sphaleron rate, i.e., the rate of real time QCD topological transitions, is a key phenomenological quantity, playing a fundamental role in several physical contexts. In heavy-ion collisions, a non-vanishing rate can lead to the so-called Chiral Magnetic Effect. In early-Universe cosmology, instead, it can be related to the rate of thermal production of QCD axions. In this talk, we present the first reliable fully non-perturbative computation of the strong sphaleron rate in $N_f=2+1$ QCD at the physical point by means of lattice simulations, in a range of temperatures going from 200 MeV to 600 MeV. Our strategy is based on the inversion of lattice correlators via a recently-proposed modified version of the Backus-Gilbert method.
△ Less
Submitted 9 May, 2024;
originally announced May 2024.
-
ML-based identification of the interface regions for coupling local and nonlocal models
Authors:
Noujoud Nader,
Patrick Diehl,
Marta D'Elia,
Christian Glusa,
Serge Prudhomme
Abstract:
Local-nonlocal coupling approaches combine the computational efficiency of local models and the accuracy of nonlocal models. However, the coupling process is challenging, requiring expertise to identify the interface between local and nonlocal regions. This study introduces a machine learning-based approach to automatically detect the regions in which the local and nonlocal models should be used i…
▽ More
Local-nonlocal coupling approaches combine the computational efficiency of local models and the accuracy of nonlocal models. However, the coupling process is challenging, requiring expertise to identify the interface between local and nonlocal regions. This study introduces a machine learning-based approach to automatically detect the regions in which the local and nonlocal models should be used in a coupling approach. This identification process uses the loading functions and provides as output the selected model at the grid points. Training is based on datasets of loading functions for which reference coupling configurations are computed using accurate coupled solutions, where accuracy is measured in terms of the relative error between the solution to the coupling approach and the solution to the nonlocal model. We study two approaches that differ from one another in terms of the data structure. The first approach, referred to as the full-domain input data approach, inputs the full load vector and outputs a full label vector. In this case, the classification process is carried out globally. The second approach consists of a window-based approach, where loads are preprocessed and partitioned into windows and the problem is formulated as a node-wise classification approach in which the central point of each window is treated individually. The classification problems are solved via deep learning algorithms based on convolutional neural networks. The performance of these approaches is studied on one-dimensional numerical examples using F1-scores and accuracy metrics. In particular, it is shown that the windowing approach provides promising results, achieving an accuracy of 0.96 and an F1-score of 0.97. These results underscore the potential of the approach to automate coupling processes, leading to more accurate and computationally efficient solutions for material science applications.
△ Less
Submitted 23 April, 2024;
originally announced April 2024.
-
Full QCD with milder topological freezing
Authors:
Claudio Bonanno,
Giuseppe Clemente,
Massimo D'Elia,
Lorenzo Maio,
Luca Parente
Abstract:
We simulate $N_f=2+1$ QCD at the physical point combining open and periodic boundary conditions in a parallel tempering framework, following the original proposal by M. Hasenbusch for $2d$ $\mathrm{CP}^{N-1}$ models, which has been recently implemented and widely employed in $4d$ $\mathrm{SU}(N)$ pure Yang-Mills theories too. We show that using this algorithm it is possible to achieve an impressiv…
▽ More
We simulate $N_f=2+1$ QCD at the physical point combining open and periodic boundary conditions in a parallel tempering framework, following the original proposal by M. Hasenbusch for $2d$ $\mathrm{CP}^{N-1}$ models, which has been recently implemented and widely employed in $4d$ $\mathrm{SU}(N)$ pure Yang-Mills theories too. We show that using this algorithm it is possible to achieve an impressive reduction of the auto-correlation time of the topological charge in dynamical fermions simulations both at zero and finite temperature up to two orders of magnitude, allowing to avoid topology freezing down to lattice spacings as fine as $a \sim 0.02$ fm. Therefore, this implementation of the Parallel Tempering on Boundary Conditions algorithm has the potential to substantially push forward the investigation of the QCD vacuum properties by means of lattice simulations.
△ Less
Submitted 22 April, 2024;
originally announced April 2024.
-
Dirac Spectral Density in $N_f = 2+1$ QCD at $T = 230$ MeV
Authors:
Andrei Alexandru,
Claudio Bonanno,
Massimo D'Elia,
Ivan Horváth
Abstract:
We compute the renormalized Dirac spectral density in $N_f = 2+1$ QCD at physical quark masses, temperature $T = 230$ MeV and system size $L_s = 3.4$ fm. To that end, we perform a point-wise continuum limit of the staggered density in lattice QCD with staggered quarks. We find, for the first time, that a clear infrared structure (IR peak) emerges in the density of Dirac operator describing dynamic…
▽ More
We compute the renormalized Dirac spectral density in $N_f = 2+1$ QCD at physical quark masses, temperature $T = 230$ MeV and system size $L_s = 3.4$ fm. To that end, we perform a point-wise continuum limit of the staggered density in lattice QCD with staggered quarks. We find, for the first time, that a clear infrared structure (IR peak) emerges in the density of Dirac operator describing dynamical quarks. Features of this structure are consistent with those previously attributed to the recently-proposed IR phase of thermal QCD. Our results (i) provide solid evidence that these IR features are stable and physical; (ii) improve the upper bound for IR-phase transition temperature $T_{\mathrm{IR}}$ so that the new window is $200 < T_{\mathrm{IR}} < 230\,$MeV; (iii) support non-restoration of anomalous U$_{\mathrm A}$(1) symmetry (chiral limit) below $T \!=\! 230\,$MeV.
△ Less
Submitted 18 April, 2024;
originally announced April 2024.
-
The $\mathrm{SU}(3)$ twisted gradient flow strong coupling without topology freezing
Authors:
Claudio Bonanno,
Jorge Luis Dasilva Golán,
Massimo D'Elia,
Margarita García Pérez,
Andrea Giorgieri
Abstract:
We investigate the role of topology in the lattice determination of the renormalized strong coupling via the gradient flow. This is done adopting the Parallel Tempering on Boundary Conditions to deal with the infamous topological freezing problem that affects standard local updating algorithms. We show that, even in the presence of severe topological freezing, both algorithms yield the same strong…
▽ More
We investigate the role of topology in the lattice determination of the renormalized strong coupling via the gradient flow. This is done adopting the Parallel Tempering on Boundary Conditions to deal with the infamous topological freezing problem that affects standard local updating algorithms. We show that, even in the presence of severe topological freezing, both algorithms yield the same strong coupling once projected onto a fixed topological sector. Moreover, we show that using a non-projected definition of the coupling leads to the same step-scaling sequence. This means that projecting the coupling onto a fixed topological sector does not affect the determination of the dynamically-generated scale of the theory $Λ$.
△ Less
Submitted 20 March, 2024;
originally announced March 2024.
-
Chiral condensate from the spectrum of the staggered Dirac operator
Authors:
Claudio Bonanno,
Francesco D'Angelo,
Massimo D'Elia
Abstract:
The chiral condensate is computed from the mode number of the staggered Dirac operator. This result is compared with those obtained with other approaches, based on the quark mass dependence of the topological susceptibility and of the pion mass.
The chiral condensate is computed from the mode number of the staggered Dirac operator. This result is compared with those obtained with other approaches, based on the quark mass dependence of the topological susceptibility and of the pion mass.
△ Less
Submitted 15 January, 2024;
originally announced January 2024.
-
Sphaleron rate as an inverse problem: a novel lattice approach
Authors:
Claudio Bonanno,
Francesco D'Angelo,
Massimo D'Elia,
Lorenzo Maio,
Manuel Naviglio
Abstract:
We compute the sphaleron rate on the lattice. We adopt a novel strategy based on the extraction of the spectral density via a modified version of the Backus-Gilbert method from finite-lattice-spacing and finite-smoothing-radius Euclidean topological charge density correlators. The physical sphaleron rate is computed by performing controlled continuum limit and zero-smoothing extrapolations both in…
▽ More
We compute the sphaleron rate on the lattice. We adopt a novel strategy based on the extraction of the spectral density via a modified version of the Backus-Gilbert method from finite-lattice-spacing and finite-smoothing-radius Euclidean topological charge density correlators. The physical sphaleron rate is computed by performing controlled continuum limit and zero-smoothing extrapolations both in pure gauge and, for the first time, in full QCD.
△ Less
Submitted 24 December, 2023;
originally announced December 2023.
-
The $θ$-dependence of the $\mathrm{SU}(N)$ critical temperature at large $N$
Authors:
Claudio Bonanno,
Massimo D'Elia,
Lorenzo Verzichelli
Abstract:
We investigate, by means of numerical lattice simulations, the $θ$-dependence of the critical deconfinement temperature of $\mathrm{SU}(N)$ gauge theories at large $N$: $T_c(θ) = T_c(0)[1-Rθ^2+O(θ^4)]$, with $R\sim O(1/N^2)$. We follow two different strategies to determine $R$, one based on the calculation of the latent heat of the transition and on the jump of the topological susceptibility at th…
▽ More
We investigate, by means of numerical lattice simulations, the $θ$-dependence of the critical deconfinement temperature of $\mathrm{SU}(N)$ gauge theories at large $N$: $T_c(θ) = T_c(0)[1-Rθ^2+O(θ^4)]$, with $R\sim O(1/N^2)$. We follow two different strategies to determine $R$, one based on the calculation of the latent heat of the transition and on the jump of the topological susceptibility at the $θ=0$ critical point, the other relying on a direct probe of $T_c(θ)$ by means of imaginary-$θ$ Monte Carlo simulations. Our results show that $R$ follows the expected large-$N$ scaling.
△ Less
Submitted 8 February, 2024; v1 submitted 19 December, 2023;
originally announced December 2023.
-
The Twisted Gradient Flow strong coupling with Parallel Tempering on Boundary Conditions
Authors:
Claudio Bonanno,
Jorge Luis Dasilva Golán,
Massimo D'Elia,
Margarita García Pérez,
Andrea Giorgieri
Abstract:
We present a proposal for calculating the running of the coupling constant of the $\mathrm{SU}(3)$ pure-gauge theory, which combines the Twisted Gradient Flow (TGF) renormalization scheme with Parallel Tempering on Boundary Conditions (PTBC). The TGF is a gradient flow-based renormalization scheme formulated in an asymmetric lattice with twisted boundary conditions. Combined with step scaling, it…
▽ More
We present a proposal for calculating the running of the coupling constant of the $\mathrm{SU}(3)$ pure-gauge theory, which combines the Twisted Gradient Flow (TGF) renormalization scheme with Parallel Tempering on Boundary Conditions (PTBC). The TGF is a gradient flow-based renormalization scheme formulated in an asymmetric lattice with twisted boundary conditions. Combined with step scaling, it has been successfully used to calculate the $\mathrm{SU}(3)$ $Λ$ parameter. As with all gradient flow-based schemes, the coupling constant is highly correlated with the topological charge and affected by topology freezing, an issue addressed by projecting the determination of the coupling onto the zero topological sector. As an alternative to the zero-charge projection, we combine TGF with PTBC, by replicating multiple copies of the same lattice, interpolating between periodic and open boundary conditions in a parallel-tempered manner. We present a first exploration of these ideas by analyzing specific ensembles of $\mathrm{SU}(3)$ lattices with and without PTBC.
△ Less
Submitted 15 December, 2023; v1 submitted 14 December, 2023;
originally announced December 2023.
-
Sphaleron rate from lattice QCD
Authors:
Claudio Bonanno,
Francesco D'Angelo,
Massimo D'Elia,
Lorenzo Maio,
Manuel Naviglio
Abstract:
We compute the sphaleron rate on the lattice from the inversion of the Euclidean time correlators of the topological charge density, performing also controlled continuum and zero-smoothing extrapolations. The correlator inversion is performed by means of a recently-proposed modification of the Backus-Gilbert method.
We compute the sphaleron rate on the lattice from the inversion of the Euclidean time correlators of the topological charge density, performing also controlled continuum and zero-smoothing extrapolations. The correlator inversion is performed by means of a recently-proposed modification of the Backus-Gilbert method.
△ Less
Submitted 23 September, 2023;
originally announced September 2023.
-
Quantum Computation of Thermal Averages for a Non-Abelian $D_4$ Lattice Gauge Theory via Quantum Metropolis Sampling
Authors:
Edoardo Ballini,
Giuseppe Clemente,
Massimo D'Elia,
Lorenzo Maio,
Kevin Zambello
Abstract:
In this paper, we show the application of the Quantum Metropolis Sampling (QMS) algorithm to a toy gauge theory with discrete non-Abelian gauge group $D_4$ in (2+1)-dimensions, discussing in general how some components of hybrid quantum-classical algorithms should be adapted in the case of gauge theories. In particular, we discuss the construction of random unitary operators which preserve gauge i…
▽ More
In this paper, we show the application of the Quantum Metropolis Sampling (QMS) algorithm to a toy gauge theory with discrete non-Abelian gauge group $D_4$ in (2+1)-dimensions, discussing in general how some components of hybrid quantum-classical algorithms should be adapted in the case of gauge theories. In particular, we discuss the construction of random unitary operators which preserve gauge invariance and act transitively on the physical Hilbert space, constituting an ergodic set of quantum Metropolis moves between gauge invariant eigenspaces, and introduce a protocol for gauge invariant measurements. Furthermore, we show how a finite resolution in the energy measurements distorts the energy and plaquette distribution measured via QMS, and propose a heuristic model that takes into account part of the deviations between numerical results and exact analytical results, whose discrepancy tends to vanish by increasing the number of qubits used for the energy measurements.
△ Less
Submitted 13 September, 2023;
originally announced September 2023.
-
The chiral condensate of $N_f=2+1$ QCD from the spectrum of the staggered Dirac operator
Authors:
Claudio Bonanno,
Francesco D'Angelo,
Massimo D'Elia
Abstract:
We compute the chiral condensate of $2+1$ QCD from the mode number of the staggered Dirac operator, performing controlled extrapolations to both the continuum and the chiral limit. We consider also alternative strategies, based on the quark mass dependence of the topological susceptibility and of the pion mass, and obtain consistent results within errors. Results are also consistent with phenomeno…
▽ More
We compute the chiral condensate of $2+1$ QCD from the mode number of the staggered Dirac operator, performing controlled extrapolations to both the continuum and the chiral limit. We consider also alternative strategies, based on the quark mass dependence of the topological susceptibility and of the pion mass, and obtain consistent results within errors. Results are also consistent with phenomenological expectations and with previous numerical determinations obtained with different lattice discretizations.
△ Less
Submitted 17 October, 2023; v1 submitted 2 August, 2023;
originally announced August 2023.
-
Sphaleron rate of $N_f=2+1$ QCD
Authors:
Claudio Bonanno,
Francesco D'Angelo,
Massimo D'Elia,
Lorenzo Maio,
Manuel Naviglio
Abstract:
We compute the sphaleron rate of $N_f=2+1$ QCD at the physical point for a range of temperatures $200$ MeV $\lesssim T \lesssim 600$ MeV. We adopt a strategy recently applied in the quenched case, based on the extraction of the rate via a modified version of the Backus-Gilbert method from finite-lattice-spacing and finite-smoothing-radius Euclidean topological charge density correlators. The physi…
▽ More
We compute the sphaleron rate of $N_f=2+1$ QCD at the physical point for a range of temperatures $200$ MeV $\lesssim T \lesssim 600$ MeV. We adopt a strategy recently applied in the quenched case, based on the extraction of the rate via a modified version of the Backus-Gilbert method from finite-lattice-spacing and finite-smoothing-radius Euclidean topological charge density correlators. The physical sphaleron rate is finally computed by performing a continuum limit at fixed physical smoothing radius, followed by a zero-smoothing extrapolation. Dynamical fermions were discretized using the staggered formulation, which is known to yield large lattice artifacts for the topological susceptibility. However, we find them to be rather mild for the sphaleron rate.
△ Less
Submitted 20 January, 2024; v1 submitted 2 August, 2023;
originally announced August 2023.
-
Quantum Algorithms for the computation of quantum thermal averages at work
Authors:
Riccardo Aiudi,
Claudio Bonanno,
Claudio Bonati,
Giuseppe Clemente,
Massimo D'Elia,
Lorenzo Maio,
Davide Rossini,
Salvatore Tirone,
Kevin Zambello
Abstract:
Recently, a variety of quantum algorithms have been devised to estimate thermal averages on a genuine quantum processor. In this paper, we consider the practical implementation of the so-called Quantum-Quantum Metropolis algorithm. As a testbed for this purpose, we simulate a basic system of three frustrated quantum spins and discuss its systematics, also in comparison with the Quantum Metropolis…
▽ More
Recently, a variety of quantum algorithms have been devised to estimate thermal averages on a genuine quantum processor. In this paper, we consider the practical implementation of the so-called Quantum-Quantum Metropolis algorithm. As a testbed for this purpose, we simulate a basic system of three frustrated quantum spins and discuss its systematics, also in comparison with the Quantum Metropolis Sampling algorithm.
△ Less
Submitted 2 August, 2023;
originally announced August 2023.
-
Spectral Observables and Gauge Field Couplings in Causal Dynamical Triangulations
Authors:
Giuseppe Clemente,
Massimo D'Elia
Abstract:
In the first part of this Chapter, we discuss the role of spectral observables, describing possible ways to build them from discretizations of the Laplace--Beltrami operator on triangulations, and how to extract useful geometric information. In the second part, we discuss how to simulate the composite system of gauge fields coupled to CDT for generic groups and dimensions, showing results in some…
▽ More
In the first part of this Chapter, we discuss the role of spectral observables, describing possible ways to build them from discretizations of the Laplace--Beltrami operator on triangulations, and how to extract useful geometric information. In the second part, we discuss how to simulate the composite system of gauge fields coupled to CDT for generic groups and dimensions, showing results in some specific case and pointing out current challenges.
△ Less
Submitted 10 July, 2023;
originally announced July 2023.
-
A scalable domain decomposition method for FEM discretizations of nonlocal equations of integrable and fractional type
Authors:
Manuel Klar,
Giacomo Capodaglio,
Marta D'Elia,
Christian Glusa,
Max Gunzburger,
Christian Vollmann
Abstract:
Nonlocal models allow for the description of phenomena which cannot be captured by classical partial differential equations. The availability of efficient solvers is one of the main concerns for the use of nonlocal models in real world engineering applications. We present a domain decomposition solver that is inspired by substructuring methods for classical local equations. In numerical experiment…
▽ More
Nonlocal models allow for the description of phenomena which cannot be captured by classical partial differential equations. The availability of efficient solvers is one of the main concerns for the use of nonlocal models in real world engineering applications. We present a domain decomposition solver that is inspired by substructuring methods for classical local equations. In numerical experiments involving finite element discretizations of scalar and vectorial nonlocal equations of integrable and fractional type, we observe improvements in solution time of up to 14.6x compared to commonly used solver strategies.
△ Less
Submitted 31 May, 2023;
originally announced June 2023.
-
Sphaleron rate from a modified Backus-Gilbert inversion method
Authors:
Claudio Bonanno,
Francesco D'Angelo,
Massimo D'Elia,
Lorenzo Maio,
Manuel Naviglio
Abstract:
We compute the sphaleron rate in quenched QCD for a temperature $T \simeq 1.24~T_c$ from the inversion of the Euclidean lattice time correlator of the topological charge density. We explore and compare two different strategies: one follows a new approach proposed in this study and consists in extracting the rate from finite lattice spacing correlators, and then in taking the continuum limit at fix…
▽ More
We compute the sphaleron rate in quenched QCD for a temperature $T \simeq 1.24~T_c$ from the inversion of the Euclidean lattice time correlator of the topological charge density. We explore and compare two different strategies: one follows a new approach proposed in this study and consists in extracting the rate from finite lattice spacing correlators, and then in taking the continuum limit at fixed smoothing radius followed by a zero-smoothing extrapolation; the other follows the traditional approach of extracting the rate after performing such double extrapolation directly on the correlator. In both cases the rate is obtained from a recently-proposed modification of the standard Backus-Gilbert procedure. The two strategies lead to compatible estimates within errors, which are then compared to previous results in the literature at the same or similar temperatures; the new strategy permits to obtain improved results, in terms of statistical and systematic uncertainties.
△ Less
Submitted 29 August, 2023; v1 submitted 26 May, 2023;
originally announced May 2023.
-
GNN-based physics solver for time-independent PDEs
Authors:
Rini Jasmine Gladstone,
Helia Rahmani,
Vishvas Suryakumar,
Hadi Meidani,
Marta D'Elia,
Ahmad Zareei
Abstract:
Physics-based deep learning frameworks have shown to be effective in accurately modeling the dynamics of complex physical systems with generalization capability across problem inputs. However, time-independent problems pose the challenge of requiring long-range exchange of information across the computational domain for obtaining accurate predictions. In the context of graph neural networks (GNNs)…
▽ More
Physics-based deep learning frameworks have shown to be effective in accurately modeling the dynamics of complex physical systems with generalization capability across problem inputs. However, time-independent problems pose the challenge of requiring long-range exchange of information across the computational domain for obtaining accurate predictions. In the context of graph neural networks (GNNs), this calls for deeper networks, which, in turn, may compromise or slow down the training process. In this work, we present two GNN architectures to overcome this challenge - the Edge Augmented GNN and the Multi-GNN. We show that both these networks perform significantly better (by a factor of 1.5 to 2) than baseline methods when applied to time-independent solid mechanics problems. Furthermore, the proposed architectures generalize well to unseen domains, boundary conditions, and materials. Here, the treatment of variable domains is facilitated by a novel coordinate transformation that enables rotation and translation invariance. By broadening the range of problems that neural operators based on graph neural networks can tackle, this paper provides the groundwork for their application to complex scientific and industrial settings.
△ Less
Submitted 27 March, 2023;
originally announced March 2023.
-
Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures
Authors:
Huaiqian You,
Xiao Xu,
Yue Yu,
Stewart Silling,
Marta D'Elia,
John Foster
Abstract:
Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. Herein, we propose a learning framework to extract a peridynamic model as a mesoscale continuum surrogate from MD simulated material fracture datas…
▽ More
Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. Herein, we propose a learning framework to extract a peridynamic model as a mesoscale continuum surrogate from MD simulated material fracture datasets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in MD displacement dataset. Inspired by the Weighted Essentially Non-Oscillatory scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from datasets without material damage, to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. Our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for single layer graphene. Our tests show that the proposed data-driven model is robust and generalizable: it is capable in modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.
△ Less
Submitted 11 January, 2023;
originally announced January 2023.
-
Phase Transitions in Particle Physics -- Results and Perspectives from Lattice Quantum Chromo-Dynamics
Authors:
Gert Aarts,
Joerg Aichelin,
Chris Allton,
Andreas Athenodorou,
Dimitrios Bachtis,
Claudio Bonanno,
Nora Brambilla,
Elena Bratkovskaya,
Mattia Bruno,
Michele Caselle,
Costanza Conti,
Roberto Contino,
Leonardo Cosmai,
Francesca Cuteri,
Luigi Del Debbio,
Massimo D'Elia,
Petros Dimopoulos,
Francesco Di Renzo,
Tetyana Galatyuk,
Jana N. Guenther,
Rachel Houtz,
Frithjof Karsch,
Andrey Yu. Kotov,
Maria Paola Lombardo,
Biagio Lucini
, et al. (16 additional authors not shown)
Abstract:
Phase transitions in a non-perturbative regime can be studied by ab initio Lattice Field Theory methods. The status and future research directions for LFT investigations of Quantum Chromo-Dynamics under extreme conditions are reviewed, including properties of hadrons and of the hypothesized QCD axion as inferred from QCD topology in different phases. We discuss phase transitions in strong interact…
▽ More
Phase transitions in a non-perturbative regime can be studied by ab initio Lattice Field Theory methods. The status and future research directions for LFT investigations of Quantum Chromo-Dynamics under extreme conditions are reviewed, including properties of hadrons and of the hypothesized QCD axion as inferred from QCD topology in different phases. We discuss phase transitions in strong interactions in an extended parameter space, and the possibility of model building for Dark Matter and Electro-Weak Symmetry Breaking. Methodological challenges are addressed as well, including new developments in Artificial Intelligence geared towards the identification of different phases and transitions.
△ Less
Submitted 11 July, 2023; v1 submitted 11 January, 2023;
originally announced January 2023.
-
Phase diagram of QCD in strong background magnetic field
Authors:
Massimo D'Elia,
Lorenzo Maio,
Francesco Sanfilippo,
Alfredo Stanzione
Abstract:
We discuss the phase diagram of QCD in the presence of a strong background magnetic field, providing numerical evidence, based on lattice simulations of QCD with $2+1$ flavours and physical quark masses, that the QCD crossover turns into a first order phase transition for large enough magnetic field, with a critical endpoint located between $eB=4$ GeV$^2$ (where we found an analytic crossover at a…
▽ More
We discuss the phase diagram of QCD in the presence of a strong background magnetic field, providing numerical evidence, based on lattice simulations of QCD with $2+1$ flavours and physical quark masses, that the QCD crossover turns into a first order phase transition for large enough magnetic field, with a critical endpoint located between $eB=4$ GeV$^2$ (where we found an analytic crossover at a pseudo-critical temperature $T_c=(98\pm3)$ MeV) and $eB=9$ GeV$^2$ (where the measured critical temperature is $T_c=(63\pm5)$ MeV).
△ Less
Submitted 22 November, 2022;
originally announced November 2022.
-
Localisation of Dirac eigenmodes and confinement in gauge theories: the Roberge-Weiss transition
Authors:
Marco Cardinali,
Massimo D'Elia,
Francesco Garosi,
Matteo Giordano
Abstract:
Ample numerical evidence from lattice calculations shows a strong connection between the confining properties of gauge theories at finite temperature and the localisation properties of the low-lying Dirac eigenmodes. In this contribution we discuss recent progress on this topic, focussing on results for QCD at imaginary chemical potential $μ_I/T = π$ at temperatures above the Roberge-Weiss transit…
▽ More
Ample numerical evidence from lattice calculations shows a strong connection between the confining properties of gauge theories at finite temperature and the localisation properties of the low-lying Dirac eigenmodes. In this contribution we discuss recent progress on this topic, focussing on results for QCD at imaginary chemical potential $μ_I/T = π$ at temperatures above the Roberge-Weiss transition temperature. These confirm the general picture of low modes turning from delocalised to localised at the deconfinement transition, in a previously unexplored setup with a genuine, physical transition in the presence of dynamical fermions. This further supports the use of Dirac eigenmodes as a tool to investigate the mechanisms behind confinement and the deconfinement transition.
△ Less
Submitted 21 November, 2022;
originally announced November 2022.
-
The QCD topological susceptibility at high temperatures via staggered fermions spectral projectors
Authors:
Andreas Athenodorou,
Claudio Bonanno,
Claudio Bonati,
Giuseppe Clemente,
Francesco D'Angelo,
Massimo D'Elia,
Lorenzo Maio,
Guido Martinelli,
Francesco Sanfilippo,
Antonino Todaro
Abstract:
The QCD topological observables are essential inputs to obtain theoretical predictions about axion phenomenology, which are of utmost importance for current and future experimental searches for this particle. Among them, we focus on the topological susceptibility, related to the axion mass. We present lattice results for the topological susceptibility in QCD at high temperatures obtained by discre…
▽ More
The QCD topological observables are essential inputs to obtain theoretical predictions about axion phenomenology, which are of utmost importance for current and future experimental searches for this particle. Among them, we focus on the topological susceptibility, related to the axion mass. We present lattice results for the topological susceptibility in QCD at high temperatures obtained by discretizing this observable via spectral projectors on eigenmodes of the staggered Dirac operator, and we compare them with those obtained with the standard gluonic definition. The adoption of the spectral discretization is motivated by the large lattice artifacts affecting the standard gluonic susceptibility, related to the choice of non-chiral fermions in the lattice action.
△ Less
Submitted 3 November, 2022;
originally announced November 2022.
-
Bayesian Nonlocal Operator Regression (BNOR): A Data-Driven Learning Framework of Nonlocal Models with Uncertainty Quantification
Authors:
Yiming Fan,
Marta D'Elia,
Yue Yu,
Habib N. Najm,
Stewart Silling
Abstract:
We consider the problem of modeling heterogeneous materials where micro-scale dynamics and interactions affect global behavior. In the presence of heterogeneities in material microstructure it is often impractical, if not impossible, to provide quantitative characterization of material response. The goal of this work is to develop a Bayesian framework for uncertainty quantification (UQ) in materia…
▽ More
We consider the problem of modeling heterogeneous materials where micro-scale dynamics and interactions affect global behavior. In the presence of heterogeneities in material microstructure it is often impractical, if not impossible, to provide quantitative characterization of material response. The goal of this work is to develop a Bayesian framework for uncertainty quantification (UQ) in material response prediction when using nonlocal models. Our approach combines the nonlocal operator regression (NOR) technique and Bayesian inference. Specifically, we use a Markov chain Monte Carlo (MCMC) method to sample the posterior probability distribution on parameters involved in the nonlocal constitutive law, and associated modeling discrepancies relative to higher fidelity computations. As an application, we consider the propagation of stress waves through a one-dimensional heterogeneous bar with randomly generated microstructure. Several numerical tests illustrate the construction, enabling UQ in nonlocal model predictions. Although nonlocal models have become popular means for homogenization, their statistical calibration with respect to high-fidelity models has not been presented before. This work is a first step towards statistical characterization of nonlocal model discrepancy in the context of homogenization.
△ Less
Submitted 6 October, 2022;
originally announced November 2022.
-
Control of fractional diffusion problems via dynamic programming equations
Authors:
Alessandro Alla,
Marta D'Elia,
Christian Glusa,
Hugo Oliveira
Abstract:
We explore the approximation of feedback control of integro-differential equations containing a fractional Laplacian term.
To obtain feedback control for the state variable of this nonlocal equation we use the Hamilton--Jacobi--Bellman equation.
It is well-known that this approach suffers from the curse of dimensionality, and to mitigate this problem we couple semi-Lagrangian schemes for the d…
▽ More
We explore the approximation of feedback control of integro-differential equations containing a fractional Laplacian term.
To obtain feedback control for the state variable of this nonlocal equation we use the Hamilton--Jacobi--Bellman equation.
It is well-known that this approach suffers from the curse of dimensionality, and to mitigate this problem we couple semi-Lagrangian schemes for the discretization of the dynamic programming principle with the use of Shepard approximation. This coupling enables approximation of high dimensional problems. Numerical convergence toward the solution of the continuous problem is provided together with linear and nonlinear examples. The robustness of the method with respect to disturbances of the system is illustrated by comparisons with an open-loop control approach.
△ Less
Submitted 18 October, 2022;
originally announced October 2022.
-
Towards glueball masses of large-$N~\mathrm{SU}(N)$ Yang-Mills theories without topological freezing via parallel tempering on boundary conditions
Authors:
Claudio Bonanno,
Massimo D'Elia,
Biagio Lucini,
Davide Vadacchino
Abstract:
Standard local updating algorithms experience a critical slowing down close to the continuum limit, which is particularly severe for topological observables. In practice, the Markov chain tends to remain trapped in a fixed topological sector. This problem further worsens at large $N$, and is known as $\mathit{topological}~\mathit{freezing}$. To mitigate it, we adopt the parallel tempering on bound…
▽ More
Standard local updating algorithms experience a critical slowing down close to the continuum limit, which is particularly severe for topological observables. In practice, the Markov chain tends to remain trapped in a fixed topological sector. This problem further worsens at large $N$, and is known as $\mathit{topological}~\mathit{freezing}$. To mitigate it, we adopt the parallel tempering on boundary conditions proposed by M. Hasenbusch. This algorithm allows to obtain a reduction of the auto-correlation time of the topological charge up to several orders of magnitude. With this strategy we are able to provide the first computation of low-lying glueball masses at large $N$ free of any systematics related to topological freezing.
△ Less
Submitted 14 October, 2022;
originally announced October 2022.
-
Probabilistic partition of unity networks for high-dimensional regression problems
Authors:
Tiffany Fan,
Nathaniel Trask,
Marta D'Elia,
Eric Darve
Abstract:
We explore the probabilistic partition of unity network (PPOU-Net) model in the context of high-dimensional regression problems and propose a general framework focusing on adaptive dimensionality reduction. With the proposed framework, the target function is approximated by a mixture of experts model on a low-dimensional manifold, where each cluster is associated with a local fixed-degree polynomi…
▽ More
We explore the probabilistic partition of unity network (PPOU-Net) model in the context of high-dimensional regression problems and propose a general framework focusing on adaptive dimensionality reduction. With the proposed framework, the target function is approximated by a mixture of experts model on a low-dimensional manifold, where each cluster is associated with a local fixed-degree polynomial. We present a training strategy that leverages the expectation maximization (EM) algorithm. During the training, we alternate between (i) applying gradient descent to update the DNN coefficients; and (ii) using closed-form formulae derived from the EM algorithm to update the mixture of experts model parameters. Under the probabilistic formulation, step (ii) admits the form of embarrassingly parallelizable weighted least-squares solves. The PPOU-Nets consistently outperform the baseline fully-connected neural networks of comparable sizes in numerical experiments of various data dimensions. We also explore the proposed model in applications of quantum computing, where the PPOU-Nets act as surrogate models for cost landscapes associated with variational quantum circuits.
△ Less
Submitted 11 June, 2023; v1 submitted 6 October, 2022;
originally announced October 2022.
-
Topological susceptibility of $N_f=2+1$ QCD from staggered fermions spectral projectors at high temperatures
Authors:
Andreas Athenodorou,
Claudio Bonanno,
Claudio Bonati,
Giuseppe Clemente,
Francesco D'Angelo,
Massimo D'Elia,
Lorenzo Maio,
Guido Martinelli,
Francesco Sanfilippo,
Antonino Todaro
Abstract:
We compute the topological susceptibility of $N_f=2+1$ QCD with physical quark masses in the high-temperature phase, using numerical simulations of the theory discretized on a space-time lattice. More precisely we estimate the topological susceptibility for five temperatures in the range from $\sim200$ MeV up to $\sim600$ MeV, adopting the spectral projectors definition of the topological charge b…
▽ More
We compute the topological susceptibility of $N_f=2+1$ QCD with physical quark masses in the high-temperature phase, using numerical simulations of the theory discretized on a space-time lattice. More precisely we estimate the topological susceptibility for five temperatures in the range from $\sim200$ MeV up to $\sim600$ MeV, adopting the spectral projectors definition of the topological charge based on the staggered Dirac operator. This strategy turns out to be effective in reducing the large lattice artifacts which affect the standard gluonic definition, making it possible to perform a reliable continuum extrapolation. Our results for the susceptibility in the explored temperature range are found to be partially in tension with previous determinations in the literature.
△ Less
Submitted 31 October, 2022; v1 submitted 18 August, 2022;
originally announced August 2022.
-
Topological susceptibility of $2d~\mathrm{CP}^1$ or $\mathrm{O}(3)$ non-linear $σ$-model: is it divergent or not?
Authors:
Claudio Bonanno,
Massimo D'Elia,
Francesca Margari
Abstract:
The topological susceptibility of $2d$ $\mathrm{CP}^{N-1}$ models is expected, based on perturbative computations, to develop a divergence in the limit $N \to 2$, where these models reduce to the well-known non-linear $\mathrm{O}(3)$ $σ$-model. The divergence is due to the dominance of instantons of arbitrarily small size and its detection by numerical lattice simulations is notoriously difficult,…
▽ More
The topological susceptibility of $2d$ $\mathrm{CP}^{N-1}$ models is expected, based on perturbative computations, to develop a divergence in the limit $N \to 2$, where these models reduce to the well-known non-linear $\mathrm{O}(3)$ $σ$-model. The divergence is due to the dominance of instantons of arbitrarily small size and its detection by numerical lattice simulations is notoriously difficult, because it is logarithmic in the lattice spacing. We approach the problem from a different perspective, studying the behavior of the model when the volume is fixed in dimensionless lattice units, where perturbative predictions are turned into more easily checkable behaviors. After testing this strategy for $N = 3$ and $4$, we apply it to $N = 2$, adopting at the same time a multicanonic algorithm to overcome the problem of rare topological fluctuations on asymptotically small lattices. Our final results fully confirm, by means of purely non-perturbative methods, the divergence of the topological susceptibility of the $2d$ $\mathrm{CP}^1$ model.
△ Less
Submitted 31 January, 2023; v1 submitted 30 July, 2022;
originally announced August 2022.
-
Towards glueball masses of large-$N$ $\mathrm{SU}(N)$ pure-gauge theories without topological freezing
Authors:
Claudio Bonanno,
Massimo D'Elia,
Biagio Lucini,
Davide Vadacchino
Abstract:
In commonly used Monte Carlo algorithms for lattice gauge theories the integrated autocorrelation time of the topological charge is known to be exponentially-growing as the continuum limit is approached. This $\mathit{topological}\,\,\textit{freezing}$, whose severity increases with the size of the gauge group, can result in potentially large systematics. To provide a direct quantification of the…
▽ More
In commonly used Monte Carlo algorithms for lattice gauge theories the integrated autocorrelation time of the topological charge is known to be exponentially-growing as the continuum limit is approached. This $\mathit{topological}\,\,\textit{freezing}$, whose severity increases with the size of the gauge group, can result in potentially large systematics. To provide a direct quantification of the latter, we focus on $\mathrm{SU}(6)$ Yang--Mills theory at a lattice spacing for which conventional methods associated to the decorrelation of the topological charge have an unbearable computational cost. We adopt the recently proposed $\mathit{parallel}\,\,\mathit{tempering}\,\,\mathit{on}\,\,\mathit{boundary}\,\,\mathit{conditions}$ algorithm, which has been shown to remove systematic effects related to topological freezing, and compute glueball masses with a typical accuracy of $2-5\%$. We observe no sizeable systematic effect in the mass of the first lowest-lying glueball states, with respect to calculations performed at nearly-frozen topological sector.
△ Less
Submitted 8 July, 2022; v1 submitted 12 May, 2022;
originally announced May 2022.
-
Nonlocal Machine Learning of Micro-Structural Defect Evolutions in Crystalline Materials
Authors:
Eduardo Augusto Barros de Moraes,
Marta D'Elia,
Mohsen Zayernouri
Abstract:
The presence and evolution of defects that appear in the manufacturing process play a vital role in the failure mechanisms of engineering materials. In particular, the collective behavior of dislocation dynamics at the mesoscale leads to avalanche, strain bursts, intermittent energy spikes, and nonlocal interactions producing anomalous features across different time- and length-scales, directly af…
▽ More
The presence and evolution of defects that appear in the manufacturing process play a vital role in the failure mechanisms of engineering materials. In particular, the collective behavior of dislocation dynamics at the mesoscale leads to avalanche, strain bursts, intermittent energy spikes, and nonlocal interactions producing anomalous features across different time- and length-scales, directly affecting plasticity, void and crack nucleation. Discrete Dislocation Dynamics (DDD) simulations are often used at the meso-level, but the cost and complexity increase dramatically with simulation time. To further understand how the anomalous features propagate to the continuum, we develop a probabilistic model for dislocation motion constructed from the position statistics obtained from DDD simulations. We obtain the continuous limit of discrete dislocation dynamics through a Probability Density Function for the dislocation motion, and propose a nonlocal transport model for the PDF. We develop a machine-learning framework to learn the parameters of the nonlocal operator with a power-law kernel, connecting the anomalous nature of DDD to the origin of its corresponding nonlocal operator at the continuum, facilitating the integration of dislocation dynamics into multi-scale, long-time material failure simulations.
△ Less
Submitted 11 May, 2022;
originally announced May 2022.
-
An asymptotically compatible coupling formulation for nonlocal interface problems with jumps
Authors:
Christian Glusa,
Marta D'Elia,
Giacomo Capodaglio,
Max Gunzburger,
Pavel B. Bochev
Abstract:
We introduce a mathematically rigorous formulation for a nonlocal interface problem with jumps and propose an asymptotically compatible finite element discretization for the weak form of the interface problem. After proving the well-posedness of the weak form, we demonstrate that solutions to the nonlocal interface problem converge to the corresponding local counterpart when the nonlocal data are…
▽ More
We introduce a mathematically rigorous formulation for a nonlocal interface problem with jumps and propose an asymptotically compatible finite element discretization for the weak form of the interface problem. After proving the well-posedness of the weak form, we demonstrate that solutions to the nonlocal interface problem converge to the corresponding local counterpart when the nonlocal data are appropriately prescribed. Several numerical tests in one and two dimensions show the applicability of our technique, its numerical convergence to exact nonlocal solutions, its convergence to the local limit when the horizons vanish, and its robustness with respect to the patch test.
△ Less
Submitted 14 March, 2022;
originally announced March 2022.
-
Machine Learning in Heterogeneous Porous Materials
Authors:
Marta D'Elia,
Hang Deng,
Cedric Fraces,
Krishna Garikipati,
Lori Graham-Brady,
Amanda Howard,
George Karniadakis,
Vahid Keshavarzzadeh,
Robert M. Kirby,
Nathan Kutz,
Chunhui Li,
Xing Liu,
Hannah Lu,
Pania Newell,
Daniel O'Malley,
Masa Prodanovic,
Gowri Srinivasan,
Alexandre Tartakovsky,
Daniel M. Tartakovsky,
Hamdi Tchelepi,
Bozo Vazic,
Hari Viswanathan,
Hongkyu Yoon,
Piotr Zarzycki
Abstract:
The "Workshop on Machine learning in heterogeneous porous materials" brought together international scientific communities of applied mathematics, porous media, and material sciences with experts in the areas of heterogeneous materials, machine learning (ML) and applied mathematics to identify how ML can advance materials research. Within the scope of ML and materials research, the goal of the wor…
▽ More
The "Workshop on Machine learning in heterogeneous porous materials" brought together international scientific communities of applied mathematics, porous media, and material sciences with experts in the areas of heterogeneous materials, machine learning (ML) and applied mathematics to identify how ML can advance materials research. Within the scope of ML and materials research, the goal of the workshop was to discuss the state-of-the-art in each community, promote crosstalk and accelerate multi-disciplinary collaborative research, and identify challenges and opportunities. As the end result, four topic areas were identified: ML in predicting materials properties, and discovery and design of novel materials, ML in porous and fractured media and time-dependent phenomena, Multi-scale modeling in heterogeneous porous materials via ML, and Discovery of materials constitutive laws and new governing equations. This workshop was part of the AmeriMech Symposium series sponsored by the National Academies of Sciences, Engineering and Medicine and the U.S. National Committee on Theoretical and Applied Mechanics.
△ Less
Submitted 4 February, 2022;
originally announced February 2022.
-
Efficient optimization-based quadrature for variational discretization of nonlocal problems
Authors:
Marco Pasetto,
Zhaoxiang Shen,
Marta D'Elia,
Xiaochuan Tian,
Nathaniel Trask,
David Kamensky
Abstract:
Casting nonlocal problems in variational form and discretizing them with the finite element (FE) method facilitates the use of nonlocal vector calculus to prove well-posedeness, convergence, and stability of such schemes. Employing an FE method also facilitates meshing of complicated domain geometries and coupling with FE methods for local problems. However, nonlocal weak problems involve the comp…
▽ More
Casting nonlocal problems in variational form and discretizing them with the finite element (FE) method facilitates the use of nonlocal vector calculus to prove well-posedeness, convergence, and stability of such schemes. Employing an FE method also facilitates meshing of complicated domain geometries and coupling with FE methods for local problems. However, nonlocal weak problems involve the computation of a double-integral, which is computationally expensive and presents several challenges. In particular, the inner integral of the variational form associated with the stiffness matrix is defined over the intersections of FE mesh elements with a ball of radius $δ$, where $δ$ is the range of nonlocal interaction. Identifying and parameterizing these intersections is a nontrivial computational geometry problem. In this work, we propose a quadrature technique where the inner integration is performed using quadrature points distributed over the full ball, without regard for how it intersects elements, and weights are computed based on the generalized moving least squares method. Thus, as opposed to all previously employed methods, our technique does not require element-by-element integration and fully circumvents the computation of element-ball intersections. This paper considers one- and two-dimensional implementations of piecewise linear continuous FE approximations, focusing on the case where the element size h and the nonlocal radius $δ$ are proportional, as is typical of practical computations. When boundary conditions are treated carefully and the outer integral of the variational form is computed accurately, the proposed method is asymptotically compatible in the limit of $h \sim δ\to 0$, featuring at least first-order convergence in L^2 for all dimensions, using both uniform and nonuniform grids.
△ Less
Submitted 28 January, 2022;
originally announced January 2022.
-
Machine-learning of nonlocal kernels for anomalous subsurface transport from breakthrough curves
Authors:
Xiao Xu,
Marta D'Elia,
Christian Glusa,
John T. Foster
Abstract:
Anomalous behavior is ubiquitous in subsurface solute transport due to the presence of high degrees of heterogeneity at different scales in the media. Although fractional models have been extensively used to describe the anomalous transport in various subsurface applications, their application is hindered by computational challenges. Simpler nonlocal models characterized by integrable kernels and…
▽ More
Anomalous behavior is ubiquitous in subsurface solute transport due to the presence of high degrees of heterogeneity at different scales in the media. Although fractional models have been extensively used to describe the anomalous transport in various subsurface applications, their application is hindered by computational challenges. Simpler nonlocal models characterized by integrable kernels and finite interaction length represent a computationally feasible alternative to fractional models; yet, the informed choice of their kernel functions still remains an open problem. We propose a general data-driven framework for the discovery of optimal kernels on the basis of very small and sparse data sets in the context of anomalous subsurface transport. Using spatially sparse breakthrough curves recovered from fine-scale particle-density simulations, we learn the best coarse-scale nonlocal model using a nonlocal operator regression technique. Predictions of the breakthrough curves obtained using the optimal nonlocal model show good agreement with fine-scale simulation results even at locations and time intervals different from the ones used to train the kernel, confirming the excellent generalization properties of the proposed algorithm. A comparison with trained classical models and with black-box deep neural networks confirms the superiority of the predictive capability of the proposed model.
△ Less
Submitted 28 January, 2022; v1 submitted 26 January, 2022;
originally announced January 2022.
-
Materials and devices for fundamental quantum science and quantum technologies
Authors:
Marco Polini,
Francesco Giazotto,
Kin Chung Fong,
Ioan M. Pop,
Carsten Schuck,
Tommaso Boccali,
Giovanni Signorelli,
Massimo D'Elia,
Robert H. Hadfield,
Vittorio Giovannetti,
Davide Rossini,
Alessandro Tredicucci,
Dmitri K. Efetov,
Frank H. L. Koppens,
Pablo Jarillo-Herrero,
Anna Grassellino,
Dario Pisignano
Abstract:
Technologies operating on the basis of quantum mechanical laws and resources such as phase coherence and entanglement are expected to revolutionize our future. Quantum technologies are often divided into four main pillars: computing, simulation, communication, and sensing & metrology. Moreover, a great deal of interest is currently also nucleating around energy-related quantum technologies. In thi…
▽ More
Technologies operating on the basis of quantum mechanical laws and resources such as phase coherence and entanglement are expected to revolutionize our future. Quantum technologies are often divided into four main pillars: computing, simulation, communication, and sensing & metrology. Moreover, a great deal of interest is currently also nucleating around energy-related quantum technologies. In this Perspective, we focus on advanced superconducting materials, van der Waals materials, and moiré quantum matter, summarizing recent exciting developments and highlighting a wealth of potential applications, ranging from high-energy experimental and theoretical physics to quantum materials science and energy storage.
△ Less
Submitted 23 January, 2022;
originally announced January 2022.
-
Nonlocal Kernel Network (NKN): a Stable and Resolution-Independent Deep Neural Network
Authors:
Huaiqian You,
Yue Yu,
Marta D'Elia,
Tian Gao,
Stewart Silling
Abstract:
Neural operators have recently become popular tools for designing solution maps between function spaces in the form of neural networks. Differently from classical scientific machine learning approaches that learn parameters of a known partial differential equation (PDE) for a single instance of the input parameters at a fixed resolution, neural operators approximate the solution map of a family of…
▽ More
Neural operators have recently become popular tools for designing solution maps between function spaces in the form of neural networks. Differently from classical scientific machine learning approaches that learn parameters of a known partial differential equation (PDE) for a single instance of the input parameters at a fixed resolution, neural operators approximate the solution map of a family of PDEs. Despite their success, the uses of neural operators are so far restricted to relatively shallow neural networks and confined to learning hidden governing laws. In this work, we propose a novel nonlocal neural operator, which we refer to as nonlocal kernel network (NKN), that is resolution independent, characterized by deep neural networks, and capable of handling a variety of tasks such as learning governing equations and classifying images. Our NKN stems from the interpretation of the neural network as a discrete nonlocal diffusion reaction equation that, in the limit of infinite layers, is equivalent to a parabolic nonlocal equation, whose stability is analyzed via nonlocal vector calculus. The resemblance with integral forms of neural operators allows NKNs to capture long-range dependencies in the feature space, while the continuous treatment of node-to-node interactions makes NKNs resolution independent. The resemblance with neural ODEs, reinterpreted in a nonlocal sense, and the stable network dynamics between layers allow for generalization of NKN's optimal parameters from shallow to deep networks. This fact enables the use of shallow-to-deep initialization techniques. Our tests show that NKNs outperform baseline methods in both learning governing equations and image classification tasks and generalize well to different resolutions and depths.
△ Less
Submitted 6 January, 2022;
originally announced January 2022.
-
Connections between nonlocal operators: from vector calculus identities to a fractional Helmholtz decomposition
Authors:
Marta D'Elia,
Mamikon Gulian,
Tadele Mengesha,
James M. Scott
Abstract:
Nonlocal vector calculus, which is based on the nonlocal forms of gradient, divergence, and Laplace operators in multiple dimensions, has shown promising applications in fields such as hydrology, mechanics, and image processing. In this work, we study the analytical underpinnings of these operators. We rigorously treat compositions of nonlocal operators, prove nonlocal vector calculus identities,…
▽ More
Nonlocal vector calculus, which is based on the nonlocal forms of gradient, divergence, and Laplace operators in multiple dimensions, has shown promising applications in fields such as hydrology, mechanics, and image processing. In this work, we study the analytical underpinnings of these operators. We rigorously treat compositions of nonlocal operators, prove nonlocal vector calculus identities, and connect weighted and unweighted variational frameworks. We combine these results to obtain a weighted fractional Helmholtz decomposition which is valid for sufficiently smooth vector fields. Our approach identifies the function spaces in which the stated identities and decompositions hold, providing a rigorous foundation to the nonlocal vector calculus identities that can serve as tools for nonlocal modeling in higher dimensions.
△ Less
Submitted 9 December, 2021;
originally announced December 2021.
-
Coupling Yang-Mills with Causal Dynamical Triangulations
Authors:
Alessandro Candido,
Giuseppe Clemente,
Massimo D'Elia,
Federico Rottoli
Abstract:
We discuss the algorithmic problem of minimal coupling gauge fields of the Yang-Mills type to Quantum Gravity in the approach known as Causal Dynamical Triangulations (CDT) as a step towards studying, ultimately, systems of gravity coupled with bosonic and fermionic matter. We first describe the algorithm for general dimensions and gauge groups and then focus on the results obtained from simulatio…
▽ More
We discuss the algorithmic problem of minimal coupling gauge fields of the Yang-Mills type to Quantum Gravity in the approach known as Causal Dynamical Triangulations (CDT) as a step towards studying, ultimately, systems of gravity coupled with bosonic and fermionic matter. We first describe the algorithm for general dimensions and gauge groups and then focus on the results obtained from simulations of 2d CDT coupled to gauge fields with U{(1)} and SU{(2)} gauge groups, where we studied both observables related to gravity and gauge fields, and compared them with analogous simulations in the static flat case.
△ Less
Submitted 7 December, 2021; v1 submitted 6 December, 2021;
originally announced December 2021.
-
Topology in high-$T$ QCD via staggered spectral projectors
Authors:
Andreas Athenodorou,
Claudio Bonanno,
Claudio Bonati,
Giuseppe Clemente,
Francesco D'Angelo,
Massimo D'Elia,
Lorenzo Maio,
Guido Martinelli,
Francesco Sanfilippo,
Antonino Todaro
Abstract:
We present preliminary lattice results for the topological susceptibility in high-$T$ $N_f=2+1$ QCD obtained discretizing this observable via spectral projectors on eigenmodes of the staggered operator, and we compare them with those obtained with the standard gluonic definition. The adoption of the spectral discretization is motivated by the large lattice artifacts affecting the continuum scaling…
▽ More
We present preliminary lattice results for the topological susceptibility in high-$T$ $N_f=2+1$ QCD obtained discretizing this observable via spectral projectors on eigenmodes of the staggered operator, and we compare them with those obtained with the standard gluonic definition. The adoption of the spectral discretization is motivated by the large lattice artifacts affecting the continuum scaling of the gluonic susceptibility at high $T$, related to the choice of non-chiral fermions in the action.
△ Less
Submitted 6 December, 2021;
originally announced December 2021.
-
Phase diagram of QCD in a magnetic background
Authors:
Massimo D'Elia,
Lorenzo Maio,
Francesco Sanfilippo,
Alfredo Stanzione
Abstract:
We provide numerical evidence that the thermal QCD crossover turns into a first order transition in the presence of large enough magnetic background fields. The critical endpoint is found to be located between $eB = 4$ GeV$^2$ (where the pseudocritical temperature is $T_c = (98 \pm 3)$ MeV) and $eB = 9$ GeV$^2$ (where the critical temperature is $T_c = (63 \pm 5)$ MeV). Results are based on the an…
▽ More
We provide numerical evidence that the thermal QCD crossover turns into a first order transition in the presence of large enough magnetic background fields. The critical endpoint is found to be located between $eB = 4$ GeV$^2$ (where the pseudocritical temperature is $T_c = (98 \pm 3)$ MeV) and $eB = 9$ GeV$^2$ (where the critical temperature is $T_c = (63 \pm 5)$ MeV). Results are based on the analysis of quark condensates and number susceptibilities, determined by lattice simulations of $N_f = 2+1$ QCD at the physical point, discretized with three different lattice spacings, $a = 0.114, 0.086$ and $0.057$ fm, via rooted stout staggered fermions and a Symanzik tree level improved pure gauge action. We also present preliminary results regarding the confining properties of the thermal theory, suggesting that they could change drastically going across the phase transition
△ Less
Submitted 22 November, 2021;
originally announced November 2021.
-
Fractional Modeling in Action: A Survey of Nonlocal Models for Subsurface Transport, Turbulent Flows, and Anomalous Materials
Authors:
Jorge Suzuki,
Mamikon Gulian,
Mohsen Zayernouri,
Marta D'Elia
Abstract:
Modeling of phenomena such as anomalous transport via fractional-order differential equations has been established as an effective alternative to partial differential equations, due to the inherent ability to describe large-scale behavior with greater efficiency than fully-resolved classical models. In this review article, we first provide a broad overview of fractional-order derivatives with a cl…
▽ More
Modeling of phenomena such as anomalous transport via fractional-order differential equations has been established as an effective alternative to partial differential equations, due to the inherent ability to describe large-scale behavior with greater efficiency than fully-resolved classical models. In this review article, we first provide a broad overview of fractional-order derivatives with a clear emphasis on the stochastic processes that underlie their use. We then survey three exemplary application areas - subsurface transport, turbulence, and anomalous materials - in which fractional-order differential equations provide accurate and predictive models. For each area, we report on the evidence of anomalous behavior that justifies the use of fractional-order models, and survey both foundational models as well as more expressive state-of-the-art models. We also propose avenues for future research, including more advanced and physically sound models, as well as tools for calibration and discovery of fractional-order models.
△ Less
Submitted 21 October, 2021;
originally announced October 2021.
-
Electromagnetic conductivity of quark-gluon plasma at non-zero baryon density
Authors:
N. Astrakhantsev,
V. V. Braguta,
M. Cardinali,
M. D'Elia,
L. Maio,
F. Sanfilippo,
A. Trunin,
A. Vasiliev
Abstract:
In this preprint we present our results on the study of the electromagnetic conductivity in dense quark-gluon plasma obtained within lattice simulations with $N_f = 2 + 1$ dynamical quarks. We employ stout improved rooted staggered quarks at the physical point and the tree-level Symanzik improved gauge action. The simulations are performed at imaginary baryon chemical potential, and the Tikhonov r…
▽ More
In this preprint we present our results on the study of the electromagnetic conductivity in dense quark-gluon plasma obtained within lattice simulations with $N_f = 2 + 1$ dynamical quarks. We employ stout improved rooted staggered quarks at the physical point and the tree-level Symanzik improved gauge action. The simulations are performed at imaginary baryon chemical potential, and the Tikhonov regularisation method is used to extract the conductivity from current-current correlators. Our results indicate an increase of QGP electromagnetic conductivity with real baryon density, and this dependence is quite strong.
△ Less
Submitted 20 October, 2021;
originally announced October 2021.
-
Localization properties of Dirac modes at the Roberge-Weiss phase transition
Authors:
Marco Cardinali,
Massimo D'Elia,
Francesco Garosi,
Matteo Giordano
Abstract:
We study the localization properties of the low-lying Dirac eigenmodes in QCD at imaginary chemical potential $\hatμ_I=π$ at temperatures above the Roberge-Weiss transition temperature $T_{\rm RW}$. We find that modes are localized up to a temperature-dependent "mobility edge" and delocalized above it, and that the mobility edge extrapolates to zero at a temperature compatible with $T_{\rm RW}$. T…
▽ More
We study the localization properties of the low-lying Dirac eigenmodes in QCD at imaginary chemical potential $\hatμ_I=π$ at temperatures above the Roberge-Weiss transition temperature $T_{\rm RW}$. We find that modes are localized up to a temperature-dependent "mobility edge" and delocalized above it, and that the mobility edge extrapolates to zero at a temperature compatible with $T_{\rm RW}$. This supports the existence of a strong connection between localization of the low Dirac modes and deconfinement, studied here for the first time in a model with a genuine deconfinement transition in the continuum limit in the presence of dynamical fermions.
△ Less
Submitted 19 October, 2021;
originally announced October 2021.
-
An optimization-based strategy for peridynamic-FEM coupling and for the prescription of nonlocal boundary conditions
Authors:
Marta D'Elia,
David Littlewood,
Jeremy Trageser,
Mauro Perego,
Pavel Bochev
Abstract:
We develop and analyze an optimization-based method for the coupling of a static peri-dynamic (PD) model and a static classical elasticity model. The approach formulates the coupling as a control problem in which the states are the solutions of the PD and classical equations, the objective is to minimize their mismatch on an overlap of the PD and classical domains, and the controls are virtual vol…
▽ More
We develop and analyze an optimization-based method for the coupling of a static peri-dynamic (PD) model and a static classical elasticity model. The approach formulates the coupling as a control problem in which the states are the solutions of the PD and classical equations, the objective is to minimize their mismatch on an overlap of the PD and classical domains, and the controls are virtual volume constraints and boundary conditions applied at the local-nonlocal interface. Our numerical tests performed on three-dimensional geometries illustrate the consistency and accuracy of our method, its numerical convergence, and its applicability to realistic engineering geometries. We demonstrate the coupling strategy as a means to reduce computational expense by confining the nonlocal model to a subdomain of interest, and as a means to transmit local (e.g., traction) boundary conditions applied at a surface to a nonlocal model in the bulk of the domain.
△ Less
Submitted 8 October, 2021;
originally announced October 2021.
-
Confining and chiral properties of QCD in extremely strong magnetic fields
Authors:
Massimo D'Elia,
Lorenzo Maio,
Francesco Sanfilippo,
Alfredo Stanzione
Abstract:
We investigate, by numerical lattice simulations, the static quark-antiquark potential, the flux tube properties and the chiral condensate for $N_f = 2+1$ QCD with physical quark masses in the presence of strong magnetic fields, going up to $eB = 9$ GeV$^2$, with continuum extrapolated results. The string tension for quark-antiquark separations longitudinal to the magnetic field is suppressed by o…
▽ More
We investigate, by numerical lattice simulations, the static quark-antiquark potential, the flux tube properties and the chiral condensate for $N_f = 2+1$ QCD with physical quark masses in the presence of strong magnetic fields, going up to $eB = 9$ GeV$^2$, with continuum extrapolated results. The string tension for quark-antiquark separations longitudinal to the magnetic field is suppressed by one order of magnitude at the largest explored magnetic field with respect to its value at zero magnetic background, but is still non-vanishing; in the transverse direction, instead, the string tension is enhanced but seems to reach a saturation at around 50 % of its value at $B = 0$. The flux tube shows a consistent suppression/enhancement of the overall amplitude, with mild modifications of its profile. Finally, we observe magnetic catalysis in the whole range of explored fields with a behavior compatible with a lowest Landau level approximation, in particular with a linear dependence of the chiral condensate on $B$ which is in agreement, within errors, with that already observed for $eB \sim 1$ GeV$^2$.
△ Less
Submitted 21 September, 2021; v1 submitted 15 September, 2021;
originally announced September 2021.
-
Peridynamic Model for Single-Layer Graphene Obtained from Coarse Grained Bond Forces
Authors:
Stewart A. Silling,
Marta D'Elia,
Yue Yu,
Huaiqian You,
Muge Fermen-Coker
Abstract:
An ordinary state-based peridynamic material model is proposed for single sheet graphene. The model is calibrated using coarse grained molecular dynamics simulations. The coarse graining method allows the dependence of bond force on bond length to be determined, including the horizon. The peridynamic model allows the horizon to be rescaled, providing a multiscale capability and allowing for substa…
▽ More
An ordinary state-based peridynamic material model is proposed for single sheet graphene. The model is calibrated using coarse grained molecular dynamics simulations. The coarse graining method allows the dependence of bond force on bond length to be determined, including the horizon. The peridynamic model allows the horizon to be rescaled, providing a multiscale capability and allowing for substantial reductions in computational cost compared with molecular dynamics. The calibrated peridynamic model is compared to experimental data on the deflection and perforation of a graphene monolayer by an atomic force microscope probe.
△ Less
Submitted 8 September, 2021;
originally announced September 2021.
-
On the fractional Laplacian of variable order
Authors:
Eric Darve,
Marta D'Elia,
Roberto Garrappa,
Andrea Giusti,
Natalia L. Rubio
Abstract:
We present a novel definition of variable-order fractional Laplacian on Rn based on a natural generalization of the standard Riesz potential. Our definition holds for values of the fractional parameter spanning the entire open set (0, n/2). We then discuss some properties of the fractional Poisson's equation involving this operator and we compute the corresponding Green function, for which we prov…
▽ More
We present a novel definition of variable-order fractional Laplacian on Rn based on a natural generalization of the standard Riesz potential. Our definition holds for values of the fractional parameter spanning the entire open set (0, n/2). We then discuss some properties of the fractional Poisson's equation involving this operator and we compute the corresponding Green function, for which we provide some instructive examples for specific problems.
△ Less
Submitted 2 September, 2021;
originally announced September 2021.
-
A data-driven peridynamic continuum model for upscaling molecular dynamics
Authors:
Huaiqian You,
Yue Yu,
Stewart Silling,
Marta D'Elia
Abstract:
Nonlocal models, including peridynamics, often use integral operators that embed lengthscales in their definition. However, the integrands in these operators are difficult to define from the data that are typically available for a given physical system, such as laboratory mechanical property tests. In contrast, molecular dynamics (MD) does not require these integrands, but it suffers from computat…
▽ More
Nonlocal models, including peridynamics, often use integral operators that embed lengthscales in their definition. However, the integrands in these operators are difficult to define from the data that are typically available for a given physical system, such as laboratory mechanical property tests. In contrast, molecular dynamics (MD) does not require these integrands, but it suffers from computational limitations in the length and time scales it can address. To combine the strengths of both methods and to obtain a coarse-grained, homogenized continuum model that efficiently and accurately captures materials' behavior, we propose a learning framework to extract, from MD data, an optimal Linear Peridynamic Solid (LPS) model as a surrogate for MD displacements. To maximize the accuracy of the learnt model we allow the peridynamic influence function to be partially negative, while preserving the well-posedness of the resulting model. To achieve this, we provide sufficient well-posedness conditions for discretized LPS models with sign-changing influence functions and develop a constrained optimization algorithm that minimizes the equation residual while enforcing such solvability conditions. This framework guarantees that the resulting model is mathematically well-posed, physically consistent, and that it generalizes well to settings that are different from the ones used during training. We illustrate the efficacy of the proposed approach with several numerical tests for single layer graphene. Our two-dimensional tests show the robustness of the proposed algorithm on validation data sets that include thermal noise, different domain shapes and external loadings, and discretizations substantially different from the ones used for training.
△ Less
Submitted 4 August, 2021;
originally announced August 2021.
