šŸ“£ Help Shape the Future of UKRI's Gateway to Research (GtR)

We're improving UKRI's Gateway to Research and are seeking your input! If you would be interested in being interviewed about the improvements we're making and to have your say about how we can make GtR more user-friendly, impactful, and effective for the Research and Innovation community, please email gateway@ukri.org.

New Physics meets the Strong Force

Lead Research Organisation: University of Edinburgh
Department Name: Sch of Physics and Astronomy

Abstract

The Standard Model (SM) is the (admittedly dull) name given to one of most profound theoretical achievements of science. It is an overwhelmingly successful description, not just of fundamental particles, but of the forces that govern their interactions. Despite its success, the SM is incomplete and leaves us with deep open questions. One such puzzle is the nature of dark matter. This mysterious material fills our universe, in greater abundance than the well-understood material forming stars and planets. To date it is only understood through its gravitational signature and its relation to standard material is unknown. A second question follows from the well-established fact that every particle has a corresponding anti-particle and that, upon contact, matter and anti-matter annihilate in an explosion of light. Given this remarkable correspondence, it is unknown how matter came to dominate anti-matter, leading to the formation of galaxies, stars and planets.

This proposal will combine advanced theoretical ideas with high-performance computing to investigate these two fundamental issues. The project distinguishes itself from other on-going research in this area by exploring the connections between dark matter, matter/anti-matter imbalance and a specific part of the Standard Model known as the strong force. The strong force binds fundamental particles called quarks and gluons into protons and neutrons, and the energy of these interactions generates nearly all of the mass of everyday matter. It also binds protons and neutrons together into nuclei, giving rise to all known elements.

But the same properties that make the strong-force so rich also make it an incredible challenge for theoretical physicists. The best known methods involve combining powerful theoretical tools with numerical calculations driven by world-class supercomputers. The first aspect of the proposed project will use these methods to understand how strongly-interacting dark matter particles might interact with each other, and investigate the experimental implications. The second part will shed light on exactly how heavy particles produced in collider experiments decay into the lighter particles that are detected and explore the remarkable fact that such decays violate symmetries relating matter and anti-matter. It is currently unknown if new physics is playing a role but, if it is, this could well be related to the mechanism that led to matter dominance in the early universe.

Planned Impact

Despite its incredible success in describing a wide variety of phenomena, both in earth-based experiments and astro-physical observations, the Standard Model of particle physics is unambiguously incomplete. This proposal will consider two key gaps in the current understanding: (i) the nature of the dark matter that is known to fill the universe but has escaped direct detection and (ii) the dominance of matter over anti-matter and its apparent contradiction with the approximate symmetries of the microscopic world. The proposed research will address the interplay of these issues with the strong force, an essential sector of the Standard Model. As a result of the theoretically challenging nature of the strong force, theoretical calculations are often either impossible or rely on uncontrolled approximations.

To overcome this, I will lead a team in combining advanced theoretical methods with large scale numerical calculations using a framework called lattice QCD. The proposal requires high-performance computing as well as advanced reconstruction techniques and the development required for the scientific investigation is expected to have consequences reaching well outside of academia. In addition, the fundamental nature of the questions elevate them to a level of interest for anyone with an underlying scientific curiosity.

The software and data analysis techniques, in particular to be developed in the component concerning advanced reconstruction techniques, are expected to have a high impact outside the physics community. The reconstruction of spectral information from data correlations is a ubiquitous problem, appearing, for example, in seismology and medical imagining. The tools we develop will be provided in an open source package and a special effort will be made to disseminate the progress outside the scientific community. Further, the ambitious nature of the calculations considered will require pre-exascale computing and the interplay between the scientific and computational development will serve as a powerful driving force on both sides. In this respect the University of Edinburgh already has a strong track record that the successful proposal will push to a new level.

Finally, this proposal will undoubtedly find resonance with the wider public. The nature of dark matter, a mysterious weakly-interacting background density with a clear gravitational signature, is a mystery so profound and simple that a wider audience can appreciate the intrigue. The same is true for the excess of matter over anti-matter, the contradiction between the nearly symmetric microscopic world and the massive imbalance that defines our daily experience. With this in mind an important output of this proposal will be dedicated public outreach, through public lectures, blog posts and online videos, describing the progress of this work to the general public. One cannot overestimate the long-term positive effect, on society and industry, of funding such fundamental research and presenting the progress in an open and accessible way.
 
Description Four years after the start of the funding period, the first work package of the research grant continues to make strong progress. In WP 1.1 "D decays at the SU(3)F point", code was developed, computational time was secured, and numerical data were generated in collaboration with OpenLatt. The resulting data included finite-volume energies that were related to the scattering amplitudes for the final-state particles of the decay. This is a necessary first step in predicting the decay amplitudes. The progress was summarized at the Lattice Conference 2022 by Fabian Joswig and at the Lattice Conference 2023 by the PI, Maxwell T. Hansen.

Most recently, the newly hired PDRA, Rajnandini Mukherjee, has joined the team to bring the calculation to its final form, with publications planned for the middle of 2025. The progress has further been accelerated through collaboration with researchers at the University of Edinburgh and CERN. This includes the development of OpenSource code in the Hadrons library to calculate the final matrix elements, which will be combined with previously determined information to reach the final experimental prediction.

In addition, the project has been extended to consider not only D decays but also B decays into multi-hadron final states. The key deliverable on this front is the completion of two publications predicting the interactions of the decay products (one to appear in PRD and the other in PRL). Again, this will be a necessary input for the full decay but represents an ambitious calculation in and of itself. In collaboration with the teams in Edinburgh and CERN, computing resources have been secured to perform the full calculation, with the allocation beginning in April 2025.

Concerning WP 1.2 'Formal mapping to physical final states', the PI and collaborators completed a theoretical derivation making it possible to extract decay amplitudes into final states with three or more particles (of direct relevance for D decays), resulting in a publication in the Journal for High Energy Physics. Concerning WP 1.2 "Formal mapping to physical final states", the PI and collaborators completed a theoretical derivation making it possible to extract decay amplitudes into final states with three or more particles (of direct relevance for D decays), resulting in a publication in the Journal for High Energy Physics. Progress has also been made in the implementation of these formulas, including the development of an open-source Python library. This work was presented by the PI's student at the Lattice Conference 2024 and formalized in a lattice proceedings posted on arXiv.org and published in the Proceedings of Science.

Work Package 2, concerning the interaction of glueballs as a dark matter candidate, is progressing well with publications expected at year-end, while significant progress has been made in Work Package 3, resulting in a publication in the Journal of High Energy Physics. The publication shows that Backus-Gilbert-like algorithms for extracting dynamical information can be applied systematically in realistic numerical calculations in a two-dimensional toy theory, as planned in WP 3.1. This publication was well received in the community and led to one of the PI's collaborators presenting a plenary at the 39th International Symposium on Lattice Field Theory.
Exploitation Route In particular, the results of Work Package 3 already show clear potential to find applications throughout the field of lattice QCD and possibly in a broader context beyond. The key findings are that spectral functions, which encode dynamical information in systems ranging from medical imaging to atmospheric science, can be extracted in a controlled manner from a numerical estimate of their Laplace transform. The key features of the method are that one has precise knowledge and partial control over the resolution function that is convolved with the desired spectral function. This means that all sources of uncertainty can be reliably estimated. In addition, the work considers strategies for removing the effects of the resolution function by exploiting analytic properties of the system under consideration.
Sectors Digital/Communication/Information Technologies (including Software)

Other

 
Title Lattice dataset for the paper "Physical-mass calculation of $\rho(770)$ and $K^*(892)$ resonance parameters via $\pi \pi$ and $K \pi$ scattering amplitudes from lattice QCD" 
Description Release for publications, currently available as arxiv preprints: arXiv:2406.19193: Physical-mass calculation of $\rho(770)$ and $K^*(892)$ resonance parameters via $\pi \pi$ and $K \pi$ scattering amplitudes from lattice QCD arXiv:2406.19194: Light and strange vector resonances from lattice QCD at physical quark masses ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ This dataset contains the Wick contraction diagrams for the correlators used in the articles mentioned above. They were computed via multiplication of the distillation meson fields with the appropriate momentum, quark content and time source, followed by a trace over dilution indices (see Appendix B of arxiv:2406.19193). Data is provided as gzip tarballs which extract into a directory tree. The name of the tarballs is hinting at the directory structure. E.g. the file K-Pi.Kpi_Kpi.P0.tar.gz will extract into K-Pi/Kpi_Kpi/P0 and all containing subdirectories. The final extracted structure is designed as follows: The data used in the GEVP (before any lattice-irrep projection) is organized in the folder tree     //P/data./..h5 where     : K-Pi, Pi-Pi     : V_V, _V, _ (in correspondence to eqs. 4,5,6 of arxiv:2406.19193, where =Kpi,Pipi, and V generically denotes the vector gamma matrices contained in the files within)     : 0,1,2,3,4     : trajectory number in MD units (for a total of 90 trajectories, ranging from 420 to 2200 in steps of 20 MD units)     : filename prefix indicating the input meson fields, with their respective gamma matrices and momentum projections, as follows:        V_V (as a Kpi subfolder)  : rpl__p__rps__p<-P>..h5, with =GammaX,GammaY,GammaZ        V_V (as a Pipi subfolder) : rpl__p__rpl__p<-P>..h5, with =GammaX,GammaY,GammaZ                Kpi_V  : rpl_Gamma5_p__rpl_Gamma5_p__rps__p<-P>..h5        Pipi_V : rpl_Gamma5_p__rpl_Gamma5_p__rpl__p<-P>..h5         Kpi_Kpi   : rpl_Gamma5_p__rpl_Gamma5_p__rps_Gamma5_p__rpl_Gamma5_p..h5        Pipi_Pipi : rpl_Gamma5_p__rpl_Gamma5_p__rpl_Gamma5_p__rpl_Gamma5_p..h5 The rpl and rps stand for the 'rhophi'-like meson fields with 'light' and 'strange' content, used in arxiv:2406.19193. The ,,,,,<-P> are in the format 'px_py_pz', and denote the 3-vectors used in the momentum-projection of the accompanying meson field, exactly matching the LHS of the equations (in the order they appear) in Appendix B of arxiv:2406.19193. The momentum conservation implies that these vectors sum to zero in any given diagram. We only compute the upper triangular part of the GEVP matrix, and enforce it to be symmetric. Each hdf5 file has the structure     /DistillationContraction/Correlators///double_array(96,2) where identifies the possible diagram topologies:     'connected' for V_V correlators (eqs. A1, A4)    'triangle' for _V correlators (eqs. A2, A5)    'direct' for _ correlators (1st diagram of Eq. A3, 1st and 2nd diagrams of A6)    'rectangle' for _ correlators (2nd diagram of Eq. A3, 3rd, 4th, 5th and 6th diagrams of A6)    'cross' for _ correlators (3rd diagram of Eq. A3) To recover the diagrams where the arrows are swapped or reversed, the appropriate momentum combinations need to be chosen, as indicated in the manuscript. In this way, for Kpi each file corresponds to exactly one correlator, but for Pipi the 2nd diagram of a given correlator will come from a file whose momenta are swapped in relation to the 1st diagram. The =0,1,...,95 denotes the time source (corresponding to 't=0' in Appendix B of arxiv:2406.19193), but where no time translation was done at this stage. The h5 dataset has shape (96,2) corresponding to 96 time slices, and their real("re") and imaginary("im") parts, respectively, in double precision. To obtain a full correlator, one sums the datasets coming from the diagrams in Appendix B, including the prefactors indicated there. The pion- and kaon-like correlators (using the pi+ and K+ interpolators from eq. arxiv:2406.19193) are under     /P/data./..h5 where =Pion,Kaon, with h5 structure      /DistillationContraction/Correlators/connected//double_array(96,2) where all naming convention is the same as for the files described above. 
Type Of Material Database/Collection of data 
Year Produced 2024 
Provided To Others? Yes  
Impact This open data allows interested researchers to reproduce the results of the paper "Physical-mass calculation of $\rho(770)$ and $K^*(892)$ resonance parameters via $\pi \pi$ and $K \pi$ scattering amplitudes from lattice QCD". 
URL https://repository.cern/doi/10.17181/zbrqg-z0m06
 
Title Lattice dataset for the paper "Physical-mass calculation of $\rho(770)$ and $K^*(892)$ resonance parameters via $\pi \pi$ and $K \pi$ scattering amplitudes from lattice QCD" 
Description Release for publications, currently available as arxiv preprints: arXiv:2406.19193: Physical-mass calculation of $\rho(770)$ and $K^*(892)$ resonance parameters via $\pi \pi$ and $K \pi$ scattering amplitudes from lattice QCD arXiv:2406.19194: Light and strange vector resonances from lattice QCD at physical quark masses ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ This dataset contains the Wick contraction diagrams for the correlators used in the articles mentioned above. They were computed via multiplication of the distillation meson fields with the appropriate momentum, quark content and time source, followed by a trace over dilution indices (see Appendix B of arxiv:2406.19193). Data is provided as gzip tarballs which extract into a directory tree. The name of the tarballs is hinting at the directory structure. E.g. the file K-Pi.Kpi_Kpi.P0.tar.gz will extract into K-Pi/Kpi_Kpi/P0 and all containing subdirectories. The final extracted structure is designed as follows: The data used in the GEVP (before any lattice-irrep projection) is organized in the folder tree     //P/data./..h5 where     : K-Pi, Pi-Pi     : V_V, _V, _ (in correspondence to eqs. 4,5,6 of arxiv:2406.19193, where =Kpi,Pipi, and V generically denotes the vector gamma matrices contained in the files within)     : 0,1,2,3,4     : trajectory number in MD units (for a total of 90 trajectories, ranging from 420 to 2200 in steps of 20 MD units)     : filename prefix indicating the input meson fields, with their respective gamma matrices and momentum projections, as follows:        V_V (as a Kpi subfolder)  : rpl__p__rps__p<-P>..h5, with =GammaX,GammaY,GammaZ        V_V (as a Pipi subfolder) : rpl__p__rpl__p<-P>..h5, with =GammaX,GammaY,GammaZ                Kpi_V  : rpl_Gamma5_p__rpl_Gamma5_p__rps__p<-P>..h5        Pipi_V : rpl_Gamma5_p__rpl_Gamma5_p__rpl__p<-P>..h5         Kpi_Kpi   : rpl_Gamma5_p__rpl_Gamma5_p__rps_Gamma5_p__rpl_Gamma5_p..h5        Pipi_Pipi : rpl_Gamma5_p__rpl_Gamma5_p__rpl_Gamma5_p__rpl_Gamma5_p..h5 The rpl and rps stand for the 'rhophi'-like meson fields with 'light' and 'strange' content, used in arxiv:2406.19193. The ,,,,,<-P> are in the format 'px_py_pz', and denote the 3-vectors used in the momentum-projection of the accompanying meson field, exactly matching the LHS of the equations (in the order they appear) in Appendix B of arxiv:2406.19193. The momentum conservation implies that these vectors sum to zero in any given diagram. We only compute the upper triangular part of the GEVP matrix, and enforce it to be symmetric. Each hdf5 file has the structure     /DistillationContraction/Correlators///double_array(96,2) where identifies the possible diagram topologies:     'connected' for V_V correlators (eqs. A1, A4)    'triangle' for _V correlators (eqs. A2, A5)    'direct' for _ correlators (1st diagram of Eq. A3, 1st and 2nd diagrams of A6)    'rectangle' for _ correlators (2nd diagram of Eq. A3, 3rd, 4th, 5th and 6th diagrams of A6)    'cross' for _ correlators (3rd diagram of Eq. A3) To recover the diagrams where the arrows are swapped or reversed, the appropriate momentum combinations need to be chosen, as indicated in the manuscript. In this way, for Kpi each file corresponds to exactly one correlator, but for Pipi the 2nd diagram of a given correlator will come from a file whose momenta are swapped in relation to the 1st diagram. The =0,1,...,95 denotes the time source (corresponding to 't=0' in Appendix B of arxiv:2406.19193), but where no time translation was done at this stage. The h5 dataset has shape (96,2) corresponding to 96 time slices, and their real("re") and imaginary("im") parts, respectively, in double precision. To obtain a full correlator, one sums the datasets coming from the diagrams in Appendix B, including the prefactors indicated there. The pion- and kaon-like correlators (using the pi+ and K+ interpolators from eq. arxiv:2406.19193) are under     /P/data./..h5 where =Pion,Kaon, with h5 structure      /DistillationContraction/Correlators/connected//double_array(96,2) where all naming convention is the same as for the files described above. 
Type Of Material Database/Collection of data 
Year Produced 2024 
Provided To Others? Yes  
Impact This open data allows interested researchers to reproduce the results of the paper "Physical-mass calculation of $\rho(770)$ and $K^*(892)$ resonance parameters via $\pi \pi$ and $K \pi$ scattering amplitudes from lattice QCD". 
URL https://repository.cern/doi/10.17181/vy9x7-bzn92
 
Title pyerrors: A python framework for error analysis of Monte Carlo data 
Description We present the pyerrors python package for statistical error analysis of Monte Carlo data. Linear error propagation using automatic differentiation in an object oriented framework is combined with the G-method for a reliable estimation of autocorrelation times. Data from different sources can easily be combined, keeping the information on the origin of error components intact throughout the analysis. pyerrors can be smoothly integrated into the existing scientific python ecosystem which allows for efficient and compact analyses. 
Type Of Material Database/Collection of data 
Year Produced 2023 
Provided To Others? Yes  
Impact This open source code has been used by many lattice QCD collaborations to efficiently analyse data. 
URL https://data.mendeley.com/datasets/7ncw242ymh
 
Description RBC/UKQCD 
Organisation Brookhaven National Laboratory
Country United States 
Sector Public 
PI Contribution My research team and I participate in weekly collaborative conference calls, participate in requests for computing resources, manage the efficient use of allocated computing resources, and collaborate on software development and maintenance. I am also a contributing author to a major collaborative publication, which is summarized below.
Collaborator Contribution The collaboration works closely together to develop software, apply for high-performance computing resources, perform intensive numerical calculations, analyse the resulting data and disseminate the results through publications and seminars.
Impact Publication of "Isospin-breaking corrections to light-meson leptonic decays from lattice simulations at physical quark masses", P. Boyle, M. Di Carlo, F. Erben, V. Gülpers, M. T. Hansen et al. (Nov 23, 2022) in JHEP 02 (2023) 242
Start Year 2022
 
Description RBC/UKQCD 
Organisation University of Southampton
Country United Kingdom 
Sector Academic/University 
PI Contribution My research team and I participate in weekly collaborative conference calls, participate in requests for computing resources, manage the efficient use of allocated computing resources, and collaborate on software development and maintenance. I am also a contributing author to a major collaborative publication, which is summarized below.
Collaborator Contribution The collaboration works closely together to develop software, apply for high-performance computing resources, perform intensive numerical calculations, analyse the resulting data and disseminate the results through publications and seminars.
Impact Publication of "Isospin-breaking corrections to light-meson leptonic decays from lattice simulations at physical quark masses", P. Boyle, M. Di Carlo, F. Erben, V. Gülpers, M. T. Hansen et al. (Nov 23, 2022) in JHEP 02 (2023) 242
Start Year 2022