Search

We present a lattice QCD calculation of the hadronic light-by-light scattering contribution to the anomalous magnetic moment of the muon using $N_f=2+1+1$ flavors of staggered quarks with masses tuned to their physical values. Our final result, in the continuum limit, reads $a_{\mu}^{\mathrm{hlbl}} = 125.5(11.6)_{\mathrm{stat}}(0.4)_{\mathrm{syst}} \times 10^{-11}$ where the first error is statistical and the second is systematic. Light, strange and charm-quark contributions are considered. In addition to the connected and leading disconnected contributions, we also include an estimate of the sub-leading disconnected diagrams. Our result is compatible with previous lattice QCD and data-driven dispersive determinations., Comment: 28 pages, 18 figures

Hadronic contributions dominate the uncertainty of the Standard Model prediction for the anomalous magnetic moment of the muon. In this work, we present results on the hadronic light-by-light contribution obtained from the evaluation of the hadronic four-point function of electromagnetic currents using the position-space formalism developed by the Mainz group. The simulations are performed with staggered fermions directly at the physical point. Several physical volumes are used to estimate finite volume effects. This direct lattice study is supplemented by considering the contribution of the light pseudoscalar pole in both finite and infinite volumes, where we reuse the pseudoscalar transition form factors that have been evaluated in previous simulations on the same ensembles., Comment: 7 pages, 6 figures, LATTICE2024 conference proceedings, submitted to PoS

We present a lattice calculation of the $\eta$ and $\eta^{\prime}$ meson masses at the physical point and in the continuum limit, based on $N_f = 2+1+1$ flavors of rooted staggered quarks. Our analysis includes gauge ensembles at the physical pion and kaon masses spread over six lattice spacings in the range [0.064-0.1315]~fm. Our main results read $m_{\eta} = 543.5(5.6)~$MeV and $m_{\eta^{\prime}} = 986(38)~$MeV, consistent with the experimental values. This is an important numerical test that supports the validity of the fourth root procedure used in the staggered quark formalism. This calculation was the first step towards extracting the pseudoscalar transition form factors of the $\eta$ and $\eta^{\prime}$ mesons that play a crucial role in the hadronic light-by-light contribution to the muon $g-2$., Comment: 19 pages, 12 figures

We present a new lattice QCD calculation of the leading order hadronic vacuum polarization contribution to the muon anomalous magnetic moment $a_\mu$. We reduce uncertainties compared to our earlier computation by $40\%$, arXiv:2002.12347. We perform simulations on finer lattices allowing for an even more accurate continuum extrapolation. We also include a small, long-distance contribution obtained using input from experiments in a low-energy regime where they all agree. Combined with other standard model contributions our result leads to a prediction that differs from the measurement of $a_\mu$ by only 0.9 standard deviations. This provides a remarkable validation of the standard model to 0.37ppm., Comment: 55 pages, 31 figures

The goal of the PolarKID project is testing a new method for the measurement of polarized sources, in order to identify all the possible instrumental systematic effects that could impact the detection of CMB B-modes of polarization. It employs the KISS (KIDs Interferometer Spectrum Survey) instrument coupled to a sky simulator and to sources such as point-like black bodies (simulating planets), a dipole (extended source) and a polarizer. We use filled-arrays Lumped Element Kinetic Inductance Detectors (LEKIDs) since they have multiple advantages when observing both in a photometry and in a polarimetry configuration, Comment: 8 pages, 1 figure, Proceeding of the SPIE conference Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy XII, SPIE Astronomical Telescopes + Instrumentation 2024

Many laboratories routinely cool samples to 10 mK, but relatively few can cool condensed matter below 1 mK. Easy access to the microkelvin range would propel fields such as quantum sensors and quantum materials. Such temperatures are achieved with adiabatic nuclear demagnetization. Existing nuclear demagnetization refrigerators (NDR) are ``single-shot'', and the recycling time is incompatible with some sub-mK experiments. Furthermore, a high cooling power is required to overcome the excess heat load of order nW on NDR pre-cooled by cryogen-free dilution refrigerators. We report the performance of an aluminum NDR designed for powerful cooling when part of a dual stage continuous NDR (CNDR). Its thermal resistance is minimized to maximize the cycling rate of the CNDR and consequently its cooling power. At the same time, its susceptibility to eddy current heating is minimized. A CNDR based on two of the aluminum NDR presented here would achieve a cooling power of approximately 40 nW at 560 $\mu$K less than six days after cooling from room temperature, with a small offset in electronic temperature that decreases as the time-dependent heat load decays.

In a somewhat forgotten paper [1] it was shown how to perform interpolations between relativistic and static computations in order to obtain results for heavy-light observables for masses from, say, $m_{\rm charm}$ to $m_{\rm bottom}$. All quantities are first continuum extrapolated and then interpolated in $1/m_h=1/m_{\rm heavy}$. Large volume computations are combined with finite volume ones where a relativistic bottom quark is accessible with small $am_{\rm bottom}$. We discuss how this strategy is extended to semi-leptonic form factors and other quantities of phenomenological interest. The essential point is to form quantities where the limit $m_h\to\infty$ is approached with power corrections O$(1/m_h)$ only. Perturbative corrections $\sim\alpha_s(m_h)^{\gamma+n}$ are cancelled in the construction of the observables. We also point out how such an approach can help to control systematics in semi-leptonic decays with just large volume data. First numerical results with $N_f = 2 + 1$ and lattice spacings down to 0.039 fm are presented in [2]., Comment: Lattice 2023 talk

We present preliminary results for B-physics from a combination of non-perturbative results in the static limit with relativistic computations satisfying $am_{\mathrm{heavy}}\ll 1$. Relativistic measurements are carried out at the physical b-quark mass using the Schr\"{o}dinger Functional in a $0.5 \ \mathrm{fm}$ box. They are connected to large volume observables through step scaling functions that trace the mass dependence between the physical charm region and the static limit, such that B-physics results can be obtained by interpolation; the procedure is designed to exactly cancel the troublesome $\alpha_s(m_{\mathrm{heavy}})^{n+\gamma}$ corrections to large mass scaling. Large volume computations for both static and relativistic quantities use CLS $N_f=2+1$ ensembles at $m_u=m_d=m_s$, and with five values of the lattice spacing down to $0.039$ fm. Our preliminary results for the b-quark mass and leptonic decay constants have competitive uncertainties, which are furthermore dominated by statistics, allowing for substantial future improvement. Here we focus on numerical results, while the underlying strategy is discussed in a companion contribution., Comment: Lattice 2023 talk

Hadronic contributions dominate the uncertainty of the standard model prediction for the anomalous magnetic moment of the muon. In this work, we describe an ongoing lattice calculation of the hadronic light-by-light contribution, performed with staggered fermions. The presence of quarks with different tastes complicates the analysis of the position-space correlation function. We present a suitable adaption of the "Mainz method". As a first numerical test, we reproduce the well-known lepton-loop contribution. Results at a single lattice spacing for the light quark contribution, using two volumes, are then discussed. Our study of the long distance behavior and finite-volume effects is supplemented by considering the contribution of the light pseudoscalar-pole. The corresponding transition form factors have been evaluated in previous simulations on the same ensembles., Comment: 7 pages, 6 figures, LATTICE2023 conference proceedings, submitted to PoS

In this report we present the status of the Mainz group's lattice QCD calculation of the pion transition form factor, which describes the interaction of an on-shell pion with two off-shell photons. This form factor is the main ingredient in the calculation of the pion-pole contribution to hadronic light-by-light scattering in the muon $g-2$. We use the $N_f = 2 + 1$ CLS gauge ensembles, and we update our previous work by including a physical pion mass ensemble (E250). We compute the transition form factor in a moving frame as well as in the pion rest frame in order to have access to a wider range of photon virtualities. In addition to the quark-line connected correlator we also compute the quark-line disconnected diagrams that contribute to the form factor. At the final stage of the analysis, the result on E250 will be combined with the previous work published in 2019 to extrapolate the form factor to the continuum and to physical quark masses., Comment: Proceedings of the 40th International Symposium on Lattice Field Theory (Lattice 2023), July 31st - August 4th, 2023, Fermilab, Batavia, Illinois, USA

The NSF-funded Robust Epidemic Surveillance and Modeling (RESUME) project successfully convened a workshop entitled "High-performance computing and large-scale data management in service of epidemiological modeling" at the University of Chicago on May 1-2, 2023. This was part of a series of workshops designed to foster sustainable and interdisciplinary co-design for predictive intelligence and pandemic prevention. The event brought together 31 experts in epidemiological modeling, high-performance computing (HPC), HPC workflows, and large-scale data management to develop a shared vision for capabilities needed for computational epidemiology to better support pandemic prevention. Through the workshop, participants identified key areas in which HPC capabilities could be used to improve epidemiological modeling, particularly in supporting public health decision-making, with an emphasis on HPC workflows, data integration, and HPC access. The workshop explored nascent HPC workflow and large-scale data management approaches currently in use for epidemiological modeling and sought to draw from approaches used in other domains to determine which practices could be best adapted for use in epidemiological modeling. This report documents the key findings and takeaways from the workshop.

The precision with which hadronic vacuum polarization (HVP) is obtained determines how accurately important observables, such as the muon anomalous magnetic moment, a_\mu, or the low-energy running of the electromagnetic coupling, \alpha, are predicted. The two most precise approaches for determining HVP are: dispersive relations combined with e+e- to hadrons cross-section data, and lattice QCD. However, the results obtained in these two approaches display significant tensions, whose origins are not understood. Here we present a framework that sheds light on this issue and, if the two approaches can be reconciled, allows them to be combined. Via this framework, we test the hypothesis that the tensions can be explained by modifying the R-ratio in different intervals of center-of-mass energy sqrt(s). As ingredients, we consider observables that have been precisely determined in both approaches. These are the leading hadronic contributions to a_\mu, to the so-called intermediate window observable and to the running of \alpha between spacelike virtualities 1GeV^2 and 10GeV^2 (for which only a preliminary lattice result exists). Our tests take into account all uncertainties and correlations, as well as uncertainties on uncertainties in the lattice results. Among our findings, the most striking is that results obtained in the two approaches can be made to agree for all three observables by modifying the \rho peak in the experimental spectrum. In particular, we find that this requires a common ~5\% increase in the contributions of the peak to each of the three observables. This finding is robust against the presence or absence of one of the constraining observables. However, such an increase is much larger than the uncertainties on the measured R-ratio. We also discuss a variety of generalizations of the methods used here, as well as the limits in the information that can be extracted..., Comment: 38 pages, 8 figures

Nested pairwise frames is a method for relative benchmarking of cell or tissue digital pathology models against manual pathologist annotations on a set of sampled patches. At a high level, the method compares agreement between a candidate model and pathologist annotations with agreement among pathologists' annotations. This evaluation framework addresses fundamental issues of data size and annotator variability in using manual pathologist annotations as a source of ground truth for model validation. We implemented nested pairwise frames evaluation for tissue classification, cell classification, and cell count prediction tasks and show results for cell and tissue models deployed on an H&E-stained melanoma dataset., Comment: 10 pages, 7 figures

The objective of our study is to determine whether using English tools to extract and normalize French medical concepts on translations provides comparable performance to French models trained on a set of annotated French clinical notes. We compare two methods: a method involving French language models and a method involving English language models. For the native French method, the Named Entity Recognition (NER) and normalization steps are performed separately. For the translated English method, after the first translation step, we compare a two-step method and a terminology-oriented method that performs extraction and normalization at the same time. We used French, English and bilingual annotated datasets to evaluate all steps (NER, normalization and translation) of our algorithms. Concerning the results, the native French method performs better than the translated English one with a global f1 score of 0.51 [0.47;0.55] against 0.39 [0.34;0.44] and 0.38 [0.36;0.40] for the two English methods tested. In conclusion, despite the recent improvement of the translation models, there is a significant performance difference between the two approaches in favor of the native French method which is more efficient on French medical texts, even with few annotated documents., Comment: 26 pages, 2 figures, 5 tables

Objective:Develop and validate an algorithm for analyzing the layout of PDF clinical documents to improve the performance of downstream natural language processing tasks. Materials and Methods: We designed an algorithm to process clinical PDF documents and extract only clinically relevant text. The algorithm consists of several steps: initial text extraction using a PDF parser, followed by classification into categories such as body text, left notes, and footers using a Transformer deep neural network architecture, and finally an aggregation step to compile the lines of a given label in the text. We evaluated the technical performance of the body text extraction algorithm by applying it to a random sample of documents that were annotated. Medical performance was evaluated by examining the extraction of medical concepts of interest from the text in their respective sections. Finally, we tested an end-to-end system on a medical use case of automatic detection of acute infection described in the hospital report. Results:Our algorithm achieved per-line precision, recall, and F1 score of 98.4, 97.0, and 97.7, respectively, for body line extraction. The precision, recall, and F1 score per document for the acute infection detection algorithm were 82.54 (95CI 72.86-91.60), 85.24 (95CI 76.61-93.70), 83.87 (95CI 76, 92-90.08) with exploitation of the results of the advanced body extraction algorithm, respectively. Conclusion:We have developed and validated a system for extracting body text from clinical documents in PDF format by identifying their layout. We were able to demonstrate that this preprocessing allowed us to obtain better performances for a common downstream task, i.e., the extraction of medical concepts in their respective sections, thus proving the interest of this method on a clinical use case., Comment: 22 pages, 5 figures

In this paper we present a first ab-initio calculation of the $\pi^0$, $\eta$ and $\eta^{\prime}$ transition form factors performed with physical light-quark masses. We provide a complete parametrization of the form factors that includes both single and double-virtual kinematics. Our results are compared with experimental measurements of the form factors in the space-like region and with the measured two-photon decay widths. In a second step, our parametrizations of the transition form factors are used to compute the dominant pseudoscalar-pole contributions to the hadronic light-by-light scattering in the muon $g-2$. Our final result reads $a_{\mu}^{\rm hlbl, ps-pole} = (85.1 \pm 5.2) \times 10^{-11}$. Although the pion-pole is dominant, we confirm that, together, the $\eta$ and $\eta^{\prime}$ provide roughly half of its contribution., Comment: 29 pages, 16 figures. Typos corrected. Add figure 6

We present results for the hadronic running of the electromagnetic coupling and the weak mixing angle from simulations of lattice QCD with $N_f=2+1$ flavours of $O(a)$-improved Wilson fermions. Using two different discretisations of the vector current, we compute the quark-connected and -disconnected contributions to the hadronic vacuum polarisation (HVP) functions $\bar{\Pi}^{\gamma\gamma}$ and $\bar{\Pi}^{Z\gamma}$ for spacelike squared momenta $Q^2\leq 7$ $\mathrm{GeV}^2$. Our results are extrapolated to the physical point using ensembles at four lattice spacings, with pion masses ranging from 130 to 420 MeV. We observe a tension of up to 3.5 standard deviations between our lattice results for $\Delta\alpha_{\rm had}^{(5)}(-Q^2)$ and estimates based on the $\textit{R}$-ratio for space-like momenta in the range $Q^2=3-7\,\rm GeV^2$. To obtain an estimate for $\Delta\alpha_\mathrm{had}^{(5)}(M_Z^2)$, we employ the Euclidean split technique. The implications for comparison with global electroweak fits are assessed., Comment: 16 pages, 7 figures, proceedings of the 39th International Symposium on Lattice Field Theory, 8th-13th August 2022, Bonn, Germany

Following the publication of the new measurement of the anomalous magnetic moment of the muon, the discrepancy between experiment and the theory prediction from the $g-2$ theory initiative has increased to $4.2\,\sigma$. Recent lattice QCD calculations predict values for the hadronic vacuum polarization contribution that are larger than the data-driven estimates, bringing the Standard Model prediction closer to the experimental measurement. Euclidean time windows in the time-momentum representation of the hadronic vacuum polarization contribution to the muon $g-2$ can help clarify the discrepancy between the phenomenological and lattice predictions. We present our calculation of the intermediate distance window contribution using $N_\mathrm{f}=2+1$ flavors of O$(a)$ improved Wilson quarks. We employ ensembles at six lattice spacings below $0.1\,$fm and pion masses down to the physical value. We present a detailed study of the continuum limit, using two discretizations of the vector current and two independent sets of improvement coefficients. Our result at the physical point displays a tension of $3.9\,\sigma$ with a recent evaluation of the intermediate window based on the data-driven method., Comment: 10 pages, 5 figures, Proceedings of the 39th International Symposium on Lattice Field Theory, 8th-13th August 2022, Bonn, Germany

The energy dependency (running) of the strength of electromagnetic interactions $\alpha$ plays an important role in precision tests of the Standard Model. The running of the former to the $Z$ pole is an input quantity for global electroweak fits, while the running of the mixing angle is susceptible to the effects of Beyond Standard Model physics, particularly at low energies. We present a computation of the hadronic vacuum polarization (HVP) contribution to the running of these electroweak couplings at the non-perturbative level in lattice QCD, in the space-like regime up to $Q^2$ momentum transfers of $7\,\mathrm{GeV}^2$. This quantity is also closely related to the HVP contribution to the muon $g-2$. We observe a tension of up to $3.5$ standard deviation between our lattice results for $\Delta\alpha^{(5)}_{\mathrm{had}}(-Q^2)$ and estimates based on the $R$-ratio for $Q^2$ in the $3$ to $7\,\mathrm{GeV}^2$ range. The tension is, however, strongly diminished when translating our result to the $Z$ pole, by employing the Euclidean split technique and perturbative QCD, which yields $\Delta\alpha^{(5)}_{\mathrm{had}}(M_Z^2)=0.027\,73(15)$. This value agrees with results based on the $R$-ratio within the quoted uncertainties, and can be used as an alternative to the latter in global electroweak fits., Comment: 6 pages, 2 figures, v2: accepted version with minor changes, talk presented at the 41st International Conference on High Energy physics - ICHEP2022

We report on our progress toward the computation of the $\pi^0$, $\eta$ and $\eta^{\prime}$ transition form factors using staggered quarks on $N_f=2+1+1$ gauge ensembles generated by the Budapest-Marseille-Wuppertal collaboration. These form factors are essential ingredients to evaluate the pseudoscalar-pole contributions to the hadronic light-by-light scattering in the muon $g-2$. Preliminary results for the pseudoscalar-pole contributions are presented, at finite lattice spacing, for all three light mesons., Comment: 10 pages, 7 figures, contribution to the 39th International Symposium on Lattice Field Theory

Hadronic light-by-light scattering is one of the virtual processes that causes the gyromagnetic factor $g$ of the muon to deviate from the value of two predicted by Dirac's theory. This process makes one of the largest contributions to the uncertainty of the Standard Model prediction for the muon $(g-2)$. Lattice QCD allows for a first-principles approach to computing this non-perturbative effect. In order to avoid power-law finite-size artifacts generated by virtual photons in lattice simulations, we follow a coordinate-space approach involving a weighted integral over the vertices of the QCD four-point function of the electromagnetic current carried by the quarks. Here we present in detail the semi-analytical calculation of the QED part of the amplitude, employing position-space perturbation theory in continuous, infinite four-dimensional Euclidean space. We also provide some useful information about a computer code for the numerical implementation of our approach that has been made public at https://github.com/RJHudspith/KQED., Comment: 84 pages, 13 figures; text matches the published version

Lattice QCD calculations of the hadronic vacuum polarization (HVP) have reached a precision where the electromagnetic (e.m.) correction can no longer be neglected. This correction is both computationally challenging and hard to validate, as it leads to ultraviolet (UV) divergences and to sizeable infrared (IR) effects associated with the massless photon. While we precisely determine the UV divergence using the operator-product expansion, we propose to introduce a separation scale $\Lambda\sim400\;$MeV into the internal photon propagator, whereby the calculation splits into a short-distance part, regulated in the UV by the lattice and in the IR by the scale $\Lambda$, and a UV-finite long-distance part to be treated with coordinate-space methods, thereby avoiding power-law finite-size effects altogether. In order to predict the long-distance part, we express the UV-regulated e.m. correction to the HVP via the forward hadronic light-by-light (HLbL) scattering amplitude and relate the latter via a dispersive sum rule to $\gamma^*\gamma^*$ fusion cross-sections. Having tested the relation by reproducing the two-loop QED vacuum polarization (VP) from the tree-level $\gamma^*\gamma^*\to e^+e^-$ cross-section, we predict the expected lattice-QCD integrand resulting from the $\gamma^*\gamma^*\to\pi^0$ process., Comment: 29 pages, 6 figures; text matches the published version

Euclidean time windows in the integral representation of the hadronic vacuum polarization contribution to the muon $g-2$ serve to test the consistency of lattice calculations and may help in tracing the origins of a potential tension between lattice and data-driven evaluations. In this paper, we present results for the intermediate time window observable computed using O($a$) improved Wilson fermions at six values of the lattice spacings below 0.1\,fm and pion masses down to the physical value. Using two different sets of improvement coefficients in the definitions of the local and conserved vector currents, we perform a detailed scaling study which results in a fully controlled extrapolation to the continuum limit without any additional treatment of the data, except for the inclusion of finite-volume corrections. To determine the latter, we use a combination of the method of Hansen and Patella and the Meyer-Lellouch-L\"uscher procedure employing the Gounaris-Sakurai parameterization for the pion form factor. We correct our results for isospin-breaking effects via the perturbative expansion of QCD+QED around the isosymmetric theory. Our result at the physical point is $a_\mu^{\mathrm{win}}=(237.30\pm0.79_{\rm stat}\pm1.22_{\rm syst})\times10^{-10}$, where the systematic error includes an estimate of the uncertainty due to the quenched charm quark in our calculation. Our result displays a tension of 3.9$\sigma$ with a recent evaluation of $a_\mu^{\mathrm{win}}$ based on the data-driven method., Comment: 43 pages, 9 figures, 10 tables; version accepted for publication: extended discussion of finite-volume corrections. Results and conclusions unchanged

We compute the hadronic light-by-light scattering contribution to the muon $g-2$ from the charm quark using lattice QCD. The calculation is performed on ensembles generated with dynamical $(u,d,s)$ quarks at the SU(3)$_{\rm f}$ symmetric point with degenerate pion and kaon masses of around 415 MeV. It includes the connected charm contribution, as well as the leading disconnected Wick contraction, involving the correlation between a charm and a light-quark loop. Cutoff effects turn out to be sizeable, which leads us to use lighter-than-physical charm masses, to employ a broad range of lattice spacings reaching down to 0.039 fm and to perform a combined charm-mass and continuum extrapolation. We use the $\eta_c$ meson to define the physical charm-mass point and obtain a final value of $a_\mu^{\rm HLbL,c} = (2.8\pm 0.5) \times 10^{-11}$, whose uncertainty is dominated by the systematics of the extrapolation. Our result is consistent with the estimate based on a simple charm-quark loop, whilst being free of any perturbative scheme dependence on the charm mass. The mixed charm-light disconnected contraction contributes a small negative amount to the final value., Comment: 21 pages, 8 figures, 9 tables

We discuss the prospects for improving the precision on the hadronic corrections to the anomalous magnetic moment of the muon, and the plans of the Muon $g-2$ Theory Initiative to update the Standard Model prediction., Comment: Contribution to the US Community Study on the Future of Particle Physics (Snowmass 2021)

In a heterogeneous medium, the wavefield can be decomposed as an infinite series known as the Born expansion. Each term of the Born expansion corresponds to a scattering order, it is thus theoretically possible to discriminate single and multiple scattering field components. Experimentally, what is actually measured is the total field in which all scattering orders interfere. Conventional imaging methods usually rely on the assumption that the multiple scattering contribution can be disregarded. In a back-scattering configuration, this assumption is valid for small depths, and begins to fail for depths larger than the scattering mean-free path $\ell_s$. It is therefore a key issue to estimate the relative amount of single and multiple scattering in experimental data. To this end, a single scattering estimator $\hat{\rho}$ computed from the reflection matrix has been introduced in order to assess the weight of single scattering in the backscattered wavefield. In this article, the meaning of this estimator is investigated and a particular attention is given to recurrent scattering. In a diffraction-limited experiment, a multiple scattering sequence is said to be recurrent if the first and last scattering events occur in the same resolution cell. Recurrent scattering is shown to be responsible for correlations between single scattering and higher scattering orders of the Born expansion, inducing a bias to the estimator $\hat{\rho}$ that should rather be termed confocal scattering ratio. Interestingly, a more robust estimator is built by projecting the reflection matrix in a focused basis. The argument is sustained by numerical simulations as well as ultrasonic data obtained around 1.5~MHz in a model medium made of nylon rods immersed in water. From a more general perspective, this work raises fundamental questions about the impact of recurrent scattering on wave imaging., Comment: 17 pages, 11 figures

We compute the hadronic running of the electromagnetic and weak couplings in lattice QCD with $N_{\mathrm{f}}=2+1$ flavors of $\mathcal{O}(a)$ improved Wilson fermions. Using two different discretizations of the vector current, we compute the quark-connected and -disconnected contributions to the hadronic vacuum polarization (HVP) functions $\bar{\Pi}^{\gamma\gamma}$ and $\bar{\Pi}^{\gamma Z}$ for Euclidean squared momenta $Q^2\leq 7\,\mathrm{GeV}^2$. Gauge field ensembles at four values of the lattice spacing and several values of the pion mass, including its physical value, are used to extrapolate the results to the physical point. The ability to perform an exact flavor decomposition allows us to present the most precise determination to date of the $\mathrm{SU}(3)$-flavor-suppressed HVP function $\bar{\Pi}^{08}$ that enters the running of $\sin^2\theta_{\mathrm{W}}$. Our results for $\bar{\Pi}^{\gamma\gamma}$, $\bar{\Pi}^{\gamma Z}$ and $\bar{\Pi}^{08}$ are presented in terms of rational functions for continuous values of $Q^2$ below $7 \,\mathrm{GeV}^2$. We observe a tension of up to $3.5$ standard deviation between our lattice results for $\Delta\alpha^{(5)}_{\mathrm{had}}(-Q^2)$ and estimates based on the $R$-ratio for space-like momenta in the range $3$--$7\,\mathrm{GeV}^2$. The tension is, however, strongly diminished when translating our result to the $Z$ pole, by employing the Euclidean split technique and perturbative QCD, which yields $\Delta\alpha^{(5)}_{\mathrm{had}}(M_Z^2)=0.027\,73(15)$ and agrees with results based on the $R$-ratio within the quoted uncertainties., Comment: 63 pages, 17 figures, 11 tables. v2: version accepted for publication in JHEP