Search

Aggregating large-scale radiotherapy planning and delivery data is crucial for advancing radiation oncology research and improving clinical practice, yet challenges persist due to the diversity of treatment planning systems (TPS), record and verify (R&V) systems, and complex data formats lacking standardized retrieval methods. We developed a robust software framework that automates the collection and integration of multi-institutional radiotherapy data from diverse TPS and R&V systems. By utilizing the unidirectional references of DICOM objects, our framework reconstructs complete patient datasets starting from Radiotherapy Treatment Records (RTRECORDs), managing tasks such as data queries, transfers, verification, and logging. It effectively maps DICOM linkages between RTRECORDs, RTPLANs, RTDOSEs, RTSTRUCTs, planning images, registrations, and associated diagnostic images, incorporating custom modules for data conversion and comprehensive error handling. Implemented across multiple institutions using various systems$-$ including ARIA, Eclipse, MOSAIQ, RayStation, MIM, Pinnacle$-$ the framework successfully collected data from two clinics over an 11-year period, aggregating data from 6,022 patients and 13,871 treatment plans with a success rate of 99.76% and an average processing time of approximately 18 minutes per patient. Ongoing efforts are extending data collection to clinics lacking DICOM Query/Retrieve capabilities, demonstrating the framework's adaptability to various clinical environments. This efficient automation of comprehensive data collection overcomes significant technical barriers, facilitating the creation of large-scale datasets that can accelerate advancements in radiation oncology., Comment: 11 pages, 5 figures

We introduce a class of efficient multiple right-hand side multigrid algorithm for domain wall fermions. The simultaneous solution for a modest number of right hand sides concurrently allows for a significant reduction in the time spent solving the coarse grid operator in a multigrid preconditioner. We introduce a preconditioned block conjuate gradient with a multigrid preconditioner, giving additional algorithmic benefit from the multiple right hand sides. There is also a very significant additional to computation rate benefit to multiple right hand sides. This both increases the arithmetic intensity in the coarse space and increases the amount of work being performed in each subroutine call, leading to excellent performance on modern GPU architectures. Further, the software implementation makes use of vendor linear algebra routines (batched GEMM) that can make use of high throughput tensor hardware on recent Nvidia, AMD and Intel GPUs. The cost of the coarse space is made sub-dominant in this algorithm, and benchmarks from the Frontier supercomputer system show up to a factor of twenty speed up over the standard red-black preconditioned conjugate gradient algorithm on a large system with physical quark masses., Comment: 33 pages

We present our study of the $\rho(770)$ and $K^*(892)$ resonances from lattice quantum chromodynamics (QCD) employing domain-wall fermions at physical quark masses. We determine the finite-volume energy spectrum in various momentum frames and obtain phase-shift parameterizations via the L\"uscher formalism, and as a final step the complex resonance poles of the $\pi \pi$ and $K \pi$ elastic scattering amplitudes via an analytical continuation of the models. By sampling a large number of representative sets of underlying energy-level fits, we also assign a systematic uncertainty to our final results. This is a significant extension to data-driven analysis methods that have been used in lattice QCD to date, due to the two-step nature of the formalism. Our final pole positions, $M+i\Gamma/2$, with all statistical and systematic errors exposed, are $M_{K^{*}} = 893(2)(8)(54)(2)~\mathrm{MeV}$ and $\Gamma_{K^{*}} = 51(2)(11)(3)(0)~\mathrm{MeV}$ for the $K^*(892)$ resonance and $M_{\rho} = 796(5)(15)(48)(2)~\mathrm{MeV}$ and $\Gamma_{\rho} = 192(10)(28)(12)(0)~\mathrm{MeV}$ for the $\rho(770)$ resonance. The four differently grouped sources of uncertainties are, in the order of occurrence: statistical, data-driven systematic, an estimation of systematic effects beyond our computation (dominated by the fact that we employ a single lattice spacing), and the error from the scale-setting uncertainty on our ensemble.

We present the first ab initio calculation at physical quark masses of scattering amplitudes describing the lightest pseudoscalar mesons interacting via the strong force in the vector channel. Using lattice quantum chromodynamics, we postdict the defining parameters for two short-lived resonances, the $\rho(770)$ and $K^*(892)$, which manifest as complex energy poles in $\pi \pi$ and $K \pi$ scattering amplitudes, respectively. The calculation proceeds by first computing the finite-volume energy spectrum of the two-hadron systems, and then determining the amplitudes from the energies using the L\"uscher formalism. The error budget includes a data-driven systematic error, obtained by scanning possible fit ranges and fit models to extract the spectrum from Euclidean correlators, as well as the scattering amplitudes from the latter. The final results, obtained by analytically continuing multiple parameterizations into the complex energy plane, are $M_\rho = 796(5)(50)~\mathrm{MeV}$, $\Gamma_\rho = 192(10)(31)~\mathrm{MeV}$, $M_{K^*} = 893(2)(54)~\mathrm{MeV}$ and $\Gamma_{K^*} = 51(2)(11)~\mathrm{MeV}$, where the subscript indicates the resonance and $M$ and $\Gamma$ stand for the mass and width, respectively, and where the first bracket indicates the statistical and the second bracket the systematic uncertainty.

We present a direct lattice QCD calculation of the $x$-dependence of the pion distribution amplitude (DA), which is performed using the quasi-DA in large momentum effective theory on a domain-wall fermion ensemble at physical quark masses and spacing $a\approx 0.084$ fm. The bare quais-DA matrix elements are renormalized in the hybrid scheme and matched to $\overline{\rm MS}$ with a subtraction of the leading renormalon in the Wilson-line mass. For the first time, we include threshold resummation in the perturbative matching onto the light-cone DA, which resums the large logarithms in the soft gluon limit at next-to-next-to-leading log. The resummed results show controlled scale-variation uncertainty within the range of momentum fraction $x\in[0.25,0.75]$ at the largest pion momentum $P_z\approx 1.85$~GeV. In addition, we apply the same analysis to quasi-DAs from a highly-improved-staggered-quark ensemble at physical pion mass and $a=0.076$ fm. By comparison we find with $2\sigma$ confidence level that the DA obtained from chiral fermions is flatter and lower near $x=0.5$., Comment: update to match the version published in journal

We present non-perturbative results for beyond the standard model kaon mixing matrix elements in the isospin symmetric limit ($m_u=m_d$) of QCD, including a complete estimate of all dominant sources of systematic error. Our results are obtained from numerical simulations of lattice QCD with $N_f = 2+1$ flavours of dynamical domain wall fermions. For the first time, these quantities are simulated directly at the physical pion mass $m_\pi$~$\sim$~$139\,\mathrm{MeV}$ for two different lattice spacings. We include data at three lattice spacings in the range $a = 0.11 $ - $ 0.07\,\mathrm{fm}$ and with pion masses ranging from the physical value up to 450$\,\mathrm{MeV}$. Compared to our earlier work, we have added both direct calculations at physical quark masses and a third lattice spacing making the removal of discretisation effects significantly more precise and eliminating the need for any significant mass extrapolation beyond the range of simulated data. We renormalise the lattice operators non-perturbatively using RI-SMOM off-shell schemes. These schemes eliminate the need to model and subtract non-perturbative pion poles that arises in the RI-MOM scheme and, since the calculations are performed with domain wall fermions, the unphysical mixing between chirality sectors is suppressed. Our results for the bag parameters in the $\overline{\mathrm{MS}}$ scheme at $3\,\mathrm{GeV}$ are $B_K~\equiv~\mathcal{B}_1 = 0.5240(17)(54)$, $\mathcal{B}_2 = 0.4794(25)(35)$, $\mathcal{B}_3 = 0.746(13)(17)$, $\mathcal{B}_4 = 0.897(02)(10)$ and $\mathcal{B}_5 = 0.6882(78)(94)$, where the first error is from lattice uncertainties and the second is the uncertainty due to the perturbative matching to $\overline{\mathrm{MS}}$., Comment: 31 pages, 18 figures, 21 tables, 2 ancillary files

I review recent research and advances in algorithms for solvers and gauge generation, with an emphasis on practical algorithms for four dimensional simulations. Particular consideration is given to advances in multigrid solvers, fourier acceleration and field transformation approaches to accelerating evolution dynamics, and to parallel tempering approaches to solving the topological tunneling problem. Particular consideration is given to the interaction between rapidly evolving computer architecture and optimal algorithms that exploit these. In this conference, nascent machine learning algorithms were separately reviewed, Comment: Lattice 2023

Four-dimensional gauge theories based on symplectic Lie groups provide elegant realisations of the microscopic origin of several new physics models. Numerical studies pursued on the lattice provide quantitative information necessary for phenomenological applications. To this purpose, we implemented Sp(2N) gauge theories using Monte Carlo techniques within Grid, a performant framework designed for the numerical study of quantum field theories on the lattice. We show the first results obtained using this library, focusing on the case-study provided by the Sp(4) theory coupled to Nas = 4 Wilson-Dirac fermions transforming in the 2-index antisymmetric representation. In particular, we discuss preliminary tests of the algorithm and we test some of its main functionalities., Comment: 9 pages, 7 figures, contribution for the 40th International Symposium on Lattice Field Theory (Lattice 2023), July 31st - August 4th, 2023, Fermi National Accelerator Laboratory, presenting the results of the paper: arXiv:2306.11649

Symplectic gauge theories coupled to matter fields lead to symmetry enhancement phenomena that have potential applications in such diverse contexts as composite Higgs, top partial compositeness, strongly interacting dark matter, and dilaton-Higgs models. These theories are also interesting on theoretical grounds, for example in reference to the approach to the large-N limit. A particularly compelling research aim is the determination of the extent of the conformal window in gauge theories with symplectic groups coupled to matter, for different groups and for field content consisting of fermions transforming in different representations. Such determination would have far-reaching implications, but requires overcoming huge technical challenges. Numerical studies based on lattice field theory can provide the quantitative information necessary to this endeavour. We developed new software to implement symplectic groups in the Monte Carlo algorithms within the Grid framework. In this paper, we focus most of our attention on the Sp(4) lattice gauge theory coupled to four (Wilson-Dirac) fermions transforming in the 2-index antisymmetric representation, as a case study. We discuss an extensive catalogue of technical tests of the algorithms and present preliminary measurements to set the stage for future large-scale numerical investigations. We also include the scan of parameter space of all asymptotically free Sp(4) lattice gauge theories coupled to varying number of fermions transforming in the antisymmetric representation., Comment: 42 pages, 16 figures. Version accepted for publication

We present a lattice calculation of the $K\to\pi\pi$ matrix elements and amplitudes with both the $\Delta I = 3/2$ and 1/2 channels and $\varepsilon'$, the measure of direct $CP$ violation. We use periodic boundary conditions (PBC), where the correct kinematics of $K\to\pi\pi$ can be achieved via an excited two-pion final state. To overcome the difficulty associated with the extraction of excited states, our previous work \cite{Bai:2015nea,RBC:2020kdj} successfully employed G-parity boundary conditions, where pions are forced to have non-zero momentum enabling the $I=0$ two-pion ground state to express the on-shell kinematics of the $K\to\pi\pi$ decay. Here instead we overcome the problem using the variational method which allows us to resolve the two-pion spectrum and matrix elements up to the relevant energy where the decay amplitude is on-shell. In this paper we report an exploratory calculation of $K\to\pi\pi$ decay amplitudes and $\varepsilon'$ using PBC on a coarser lattice size of $24^3\times64$ with inverse lattice spacing $a^{-1}=1.023$ GeV and the physical pion and kaon masses. The results are promising enough to motivate us to continue our measurements on finer lattice ensembles in order to improve the precision in the near future.

A study of two-pion scattering for the isospin channels, $I=0$ and $I=2$, using lattice QCD is presented. M\"obius domain wall fermions on top of the Iwasaki-DSDR gauge action for gluons with periodic boundary conditions are used for the lattice computations which are carried out on two ensembles of gauge field configurations generated by the RBC and UKQCD collaborations with physical masses, inverse lattice spacings of 1.023 and 1.378 GeV, and spatial extents of $L=4.63$ and 4.58 fm, respectively. The all-to-all propagator method is employed to compute a matrix of correlation functions of two-pion operators. The generalized eigenvalue problem (GEVP) is solved for a matrix of correlation functions to extract phase shifts with multiple states, two pions with a non-zero relative momentum as well as two pions at rest. Our results for phase shifts for both $I=0$ and $I=2$ channels are consistent with and the Roy Equation and chiral perturbation theory, though at this preliminary stage our errors for $I=0$ are large. An important outcome of this work is that we are successful in extracting two-pion excited states, which are useful for studying $K\to\pi\pi$ decay, on physical-mass ensembles using GEVP.

We construct an approximate trivializing map by using a Schwinger-Dyson equation. The advantage of this method is that: (1) The basis for the flow kernel can be chosen arbitrarily by hand. (2) It can be applied to the general action of interest. (3) The coefficients in the kernel are determined by lattice estimates of the observables, which does not require analytic calculations beforehand. We perform the HMC with the effective action obtained by the Schwinger-Dyson method, and show that we can have better control of the effective action than the known $t$-expansion construction. However, the algorithmic overhead is still large and overwhelming the gain though faster decorrelation is observed for long-range observables in some cases. This contribution reports the preliminary results of this attempt., Comment: 19 pages, 12 figures, talk presented at the 39th International Symposium on Lattice Field Theory (Lattice 2022)

We report on the physical-point RBC/UKQCD calculation of the leading isospin-breaking corrections to light-meson leptonic decays. This is highly relevant for future precision tests in the flavour physics sector, in particular the first-row unitarity of the Cabibbo-Kobayashi-Maskawa matrix containing the elements $V_{us}$ and $V_{ud}$. The simulations were performed using Domain-Wall fermions for $2+1$ flavours, and with isospin-breaking effects included perturbatively in the path integral through order $\alpha$ and $(m_u - m_d)/\Lambda _{\mathrm{QCD}}$. We use QED$_{\mathrm{L}}$ for the inclusion of electromagnetism, and discuss here the non-locality of this prescription which has significant impact on the infinite-volume extrapolation., Comment: Proceedings for The 39th International Symposium on Lattice Field Theory, 8th-13th August, 2022, Rheinische Friedrich-Wilhelms-Universit\"at Bonn, Bonn, Germany

Resonances play an important role in Standard Model phenomenology. In particular, hadronic resonances feature in $B$ and $D$ decays, which can be central for New Physics searches. Lattice QCD simulations combined with the finite-volume method can nowadays be used to reliably study strongly coupled scattering processes such as $K\pi$ and thus the hadronic resonance $K^*$. In this work, we approach $K\pi$ scattering on a domain-wall $N_f = 2+1$ RBC-UKQCD ensemble at a physical pion mass. We use the distillation method within Grid and Hadrons software to compute sets of operator basis. That allows solving an eigenvalue problem to extract the low-energy finite-volume spectra, which are then translated into scattering information. We update the state of the calculation by reviewing the smearing process, outlining the variational analysis and concluding by showing preliminary data., Comment: arXiv admin note: substantial text overlap with arXiv:2112.09804

The decreasing uncertainties in theoretical predictions and experimental measurements of several hadronic observables related to weak processes, which in many cases are now smaller than $\mathrm{O}(1\%)$, require theoretical calculations to include subleading corrections that were neglected so far. Precise determinations of leptonic and semi-leptonic decay rates, including QED and strong isospin-breaking effects, can play a central role in solving the current tensions in the first-row unitarity of the CKM matrix. In this work we present the first RBC/UKQCD lattice calculation of the isospin-breaking corrections to the ratio of leptonic decay rates of kaons and pions into muons and neutrinos. The calculation is performed with $N_\mathrm{f}=2+1$ dynamical quarks close to the physical point and domain wall fermions in the M\"obius formulation are employed. Long-distance QED interactions are included according to the $\mathrm{QED_L}$ prescription and the crucial role of finite-volume electromagnetic corrections in the determination of leptonic decay rates, which produce a large systematic uncertainty, is extensively discussed. Finally, we study the different sources of uncertainty on $|V_\mathrm{us}|/|V_\mathrm{ud}|$ and observe that, if finite-volume systematics can be reduced, the error from isospin-breaking corrections is potentially sub-dominant in the final precision of the ratio of the CKM matrix elements., Comment: 66 pages, 23 figures and 4 tables