Leon Hostetler
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Leon Hostetler has a total of 260 co-authors. Here are the top co-authors by number of shared publications:
- Carleton DeTar
University of Utah
4 shared publications - Aida X. El-Khadra
University of Illinois, Department of Physics, Urbana, IL, USA
3 shared publications - Alejandro Vaquero
University of Utah
3 shared publications Alexei Bazavov
Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, Michigan 48824, USA
3 shared publications- Alexei Bazavov
Michigan State University
3 shared publications - Andreas S. Kronfeld
Fermi National Accelerator Laboratory
3 shared publications
Latest Publications
- Top Cited PapersOpen Access
- 10 September 2025
- Published by Elsevier in Physics Reports
- Vol. 1143, 1-158
The anomalous magnetic moment of the muon in the Standard Model: an update
- 10 September 2025
- Published by Elsevier in Physics Reports
- Vol. 1143, 1-158
Top Cited PapersOpen AccessThe anomalous magnetic moment of the muon in the Standard Model: an update
- Open Access
- 1 July 2025
- Published by American Physical Society (APS) in Physical Review Letters
- Vol. 135 (1) , 011901
Hadronic Vacuum Polarization for the Muon from Lattice QCD: Long-Distance and Full Light-Quark Connected ContributionWe present results for the dominant light-quark connected contribution to the long-distance window of the hadronic vacuum polarization (HVP) contribution to the muon g − 2 from lattice quantum chromodynamics. Specifically, with a new determination of the lattice scale on MILC’s physical-mass HISQ ensembles, using the Ω − baryon mass, we obtain a result of a μ l l , LD ( conn ) = 400.2 ( 2.3 ) stat ( 3.7 ) syst [ 4.3 ] total × 10 − 10 . Summing this result with our recent determinations of the light-quark connected contributions to the short- and intermediate-distance windows, we obtain a subpercent precision determination of the light-quark-connected contribution to HVP of a μ l l ( conn ) = 655.2 ( 2.3 ) stat ( 3.9 ) syst [ 4.5 ] total × 10 − 10 . Finally, as a consistency check, we verify that an independent analysis of the full contribution is in agreement with the sum of individual windows. We discuss our future plans for improvements of our HVP calculations to meet the target precision of the Fermilab g − 2 experiment.
- 1 July 2025
- Published by American Physical Society (APS) in Physical Review Letters
- Vol. 135 (1) , 011901
Open AccessHadronic Vacuum Polarization for the Muon from Lattice QCD: Long-Distance and Full Light-Quark Connected ContributionWe present results for the dominant light-quark connected contribution to the long-distance window of the hadronic vacuum polarization (HVP) contribution to the muon g − 2 from lattice quantum chromodynamics. Specifically, with a new determination of the lattice scale on MILC’s physical-mass HISQ ensembles, using the Ω − baryon mass, we obtain a result of a μ l l , LD ( conn ) = 400.2 ( 2.3 ) stat ( 3.7 ) syst [ 4.3 ] total × 10 − 10 . Summing this result with our recent determinations of the light-quark connected contributions to the short- and intermediate-distance windows, we obtain a subpercent precision determination of the light-quark-connected contribution to HVP of a μ l l ( conn ) = 655.2 ( 2.3 ) stat ( 3.9 ) syst [ 4.5 ] total × 10 − 10 . Finally, as a consistency check, we verify that an independent analysis of the full contribution is in agreement with the sum of individual windows. We discuss our future plans for improvements of our HVP calculations to meet the target precision of the Fermilab g − 2 experiment.
- Open Access
- 12 May 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (9)
Hadronic vacuum polarization for the muon from lattice QCD: Complete short and intermediate windowsWe present complete results for the hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment a μ in the short- and intermediate-distance window regions, which account for roughly 10% and 35% of the total HVP contribution to a μ , respectively. In particular, we perform lattice-QCD calculations for the isospin-symmetric connected and disconnected contributions, as well as corrections due to strong-isospin breaking. For the short-distance window observables, we investigate the so-called log-enhancement effects as well as the significant oscillations associated with staggered quarks in this region. For the dominant, isospin-symmetric light-quark-connected contribution, we obtain a μ l l , SD ( conn ) = 48.139 ( 11 ) stat ( 91 ) syst [ 92 ] total × 10 − 10 and a μ l l , W ( conn ) = 206.90 ( 14 ) stat ( 61 ) syst [ 63 ] total × 10 − 10 . We use Bayesian model averaging to fully estimate the covariance matrix between the individual contributions. Our determinations of the complete window contributions are a μ SD = 69.05 ( 1 ) stat ( 21 ) syst [ 21 ] total × 10 − 10 and a μ W = 236.45 ( 17 ) stat ( 83 ) syst [ 85 ] total × 10 − 10 . This work is part of our ongoing effort to compute all contributions to HVP with an overall uncertainty at the few-permille level. Published by the American Physical Society 2025
- 12 May 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (9)
Open AccessHadronic vacuum polarization for the muon from lattice QCD: Complete short and intermediate windowsWe present complete results for the hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment a μ in the short- and intermediate-distance window regions, which account for roughly 10% and 35% of the total HVP contribution to a μ , respectively. In particular, we perform lattice-QCD calculations for the isospin-symmetric connected and disconnected contributions, as well as corrections due to strong-isospin breaking. For the short-distance window observables, we investigate the so-called log-enhancement effects as well as the significant oscillations associated with staggered quarks in this region. For the dominant, isospin-symmetric light-quark-connected contribution, we obtain a μ l l , SD ( conn ) = 48.139 ( 11 ) stat ( 91 ) syst [ 92 ] total × 10 − 10 and a μ l l , W ( conn ) = 206.90 ( 14 ) stat ( 61 ) syst [ 63 ] total × 10 − 10 . We use Bayesian model averaging to fully estimate the covariance matrix between the individual contributions. Our determinations of the complete window contributions are a μ SD = 69.05 ( 1 ) stat ( 21 ) syst [ 21 ] total × 10 − 10 and a μ W = 236.45 ( 17 ) stat ( 83 ) syst [ 85 ] total × 10 − 10 . This work is part of our ongoing effort to compute all contributions to HVP with an overall uncertainty at the few-permille level. Published by the American Physical Society 2025
- Open Access
- 24 April 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (7)
Efficient state preparation for the Schwinger model with a theta term
We present a comparison of different quantum state preparation algorithms and their overall efficiency for the Schwinger model with a theta term. While adiabatic state preparation is proved to be effective, in practice it leads to large gate counts to prepare the ground state. The quantum approximate optimization algorithm (QAOA) provides excellent results while keeping the counts small by design, at the cost of an expensive classical minimization process. We introduce a “blocked” modification of the Schwinger Hamiltonian to be used in the QAOA that further decreases the length of the algorithms as the size of the problem is increased. The rodeo algorithm (RA) provides a powerful tool to efficiently prepare any eigenstate of the Hamiltonian, as long as its overlap with the initial guess is large enough. We obtain the best results when combining the blocked QAOA ansatz and the RA, as this provides an excellent initial state with a relatively short algorithm without the need to perform any classical steps for large problem sizes.
- 24 April 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (7)
Open AccessEfficient state preparation for the Schwinger model with a theta term
We present a comparison of different quantum state preparation algorithms and their overall efficiency for the Schwinger model with a theta term. While adiabatic state preparation is proved to be effective, in practice it leads to large gate counts to prepare the ground state. The quantum approximate optimization algorithm (QAOA) provides excellent results while keeping the counts small by design, at the cost of an expensive classical minimization process. We introduce a “blocked” modification of the Schwinger Hamiltonian to be used in the QAOA that further decreases the length of the algorithms as the size of the problem is increased. The rodeo algorithm (RA) provides a powerful tool to efficiently prepare any eigenstate of the Hamiltonian, as long as its overlap with the initial guess is large enough. We obtain the best results when combining the blocked QAOA ansatz and the RA, as this provides an excellent initial state with a relatively short algorithm without the need to perform any classical steps for large problem sizes.
- Open Access
- 19 February 2025
- Published by Sissa Medialab Srl
- Vol. 466, 046
Improving HISQ propagator solves using deflation
Typically, the conjugate gradient (CG) algorithm employs mixed precision and even-odd preconditioning to compute propagators for highly improved staggered quarks (HISQ). This approach suffers from critical slowing down as the light quark mass is decreased to its physical value. Multigrid is one alternative to combat critical slowing down; however, it involves setup costs that are not always easy to amortize. We consider deflation, which can also remove critical slowing down, but incurs its own setup cost to compute eigenvectors. Results using the MILC and QUDA software libraries to generate eigenvectors and to perform deflated solves on lattices up to $144^3 \times 288$ (with lattice spacing 0.04 fm) and with a range of quark masses from the physical strange down to the physical light quark values will be presented. We compare with CG and comment on deflation versus multigrid.- 19 February 2025
- Published by Sissa Medialab Srl
- Vol. 466, 046
Open AccessImproving HISQ propagator solves using deflation
Typically, the conjugate gradient (CG) algorithm employs mixed precision and even-odd preconditioning to compute propagators for highly improved staggered quarks (HISQ). This approach suffers from critical slowing down as the light quark mass is decreased to its physical value. Multigrid is one alternative to combat critical slowing down; however, it involves setup costs that are not always easy to amortize. We consider deflation, which can also remove critical slowing down, but incurs its own setup cost to compute eigenvectors. Results using the MILC and QUDA software libraries to generate eigenvectors and to perform deflated solves on lattices up to $144^3 \times 288$ (with lattice spacing 0.04 fm) and with a range of quark masses from the physical strange down to the physical light quark values will be presented. We compare with CG and comment on deflation versus multigrid. - Open Access
- 29 March 2024
- Published by American Physical Society (APS) in Physical Review D
- Vol. 109 (5) , 054514
Symmetry breaking in an extended O(2) model
Motivated by attempts to quantum simulate lattice models with continuous Abelian symmetries using discrete approximations, we study an extended-O(2) model in two dimensions that differs from the ordinary O(2) model by the addition of an explicit symmetry breaking term . Its coupling allows to smoothly interpolate between the O(2) model () and a -state clock model (). In the latter case, a -state clock model can also be defined for noninteger values of . Thus, such a limit can also be considered as an analytic continuation of an ordinary -state clock model to noninteger . In previous work, we established the phase diagram for noninteger in the infinite coupling limit (). We showed that there is a second-order phase transition at low temperature and a crossover at high temperature. In this work, we seek to establish the phase diagram at finite values of the coupling using Monte Carlo and tensor methods. We show that for noninteger , the second-order phase transition at low temperature and crossover at high temperature persist to finite coupling. For integer , 3, 4, we know there is a second-order phase transition at infinite coupling (i.e. the well-known clock models). At finite coupling, we find that the critical exponents for , 4 vary with the coupling, and for the transition may turn into a Berezinskii-Kosterlitz-Thouless transition at small coupling. We comment on the similarities and differences of the phase diagrams with those of quantum simulators of the Abelian-Higgs model based on ladder-shaped arrays of Rydberg atoms.- 29 March 2024
- Published by American Physical Society (APS) in Physical Review D
- Vol. 109 (5) , 054514
Open AccessSymmetry breaking in an extended O(2) model
Motivated by attempts to quantum simulate lattice models with continuous Abelian symmetries using discrete approximations, we study an extended-O(2) model in two dimensions that differs from the ordinary O(2) model by the addition of an explicit symmetry breaking term . Its coupling allows to smoothly interpolate between the O(2) model () and a -state clock model (). In the latter case, a -state clock model can also be defined for noninteger values of . Thus, such a limit can also be considered as an analytic continuation of an ordinary -state clock model to noninteger . In previous work, we established the phase diagram for noninteger in the infinite coupling limit (). We showed that there is a second-order phase transition at low temperature and a crossover at high temperature. In this work, we seek to establish the phase diagram at finite values of the coupling using Monte Carlo and tensor methods. We show that for noninteger , the second-order phase transition at low temperature and crossover at high temperature persist to finite coupling. For integer , 3, 4, we know there is a second-order phase transition at infinite coupling (i.e. the well-known clock models). At finite coupling, we find that the critical exponents for , 4 vary with the coupling, and for the transition may turn into a Berezinskii-Kosterlitz-Thouless transition at small coupling. We comment on the similarities and differences of the phase diagrams with those of quantum simulators of the Abelian-Higgs model based on ladder-shaped arrays of Rydberg atoms. - Open Access
Symmetry Breaking in an Extended O(2) Model
Motivated by attempts to quantum simulate lattice models with continuous Abelian symmetries using discrete approximations, we study an extended-O(2) model in two dimensions that differs from the ordinary O(2) model by the addition of an explicit symmetry breaking term $-h_q\cos(q\varphi)$. Its coupling $h_q$ allows to smoothly interpolate between the O(2) model ($h_q=0$) and a $q$-state clock model ($h_q\rightarrow\infty$). In the latter case, a $q$-state clock model can also be defined for noninteger values of $q$. Thus, such a limit can also be considered as an analytic continuation of an ordinary $q$-state clock model to noninteger $q$. In previous work, we established the phase diagram for noninteger $q$ in the infinite coupling limit ($h_q\rightarrow\infty$). We showed that there is a second-order phase transition at low temperature and a crossover at high temperature. In this work, we seek to establish the phase diagram at finite values of the coupling using Monte Carlo and tensor methods. We show that for noninteger $q$, the second-order phase transition at low temperature and crossover at high temperature persist to finite coupling. For integer $q=2,3,4$, we know there is a second-order phase transition at infinite coupling (i.e. the well-known clock models). At finite coupling, we find that the critical exponents for $q=3,4$ vary with the coupling, and for $q=4$ the transition may turn into a BKT transition at small coupling. We comment on the similarities and differences of the phase diagrams with those of quantum simulators of the Abelian-Higgs model based on ladder-shaped arrays of Rydberg atoms.Open AccessSymmetry Breaking in an Extended O(2) Model
Motivated by attempts to quantum simulate lattice models with continuous Abelian symmetries using discrete approximations, we study an extended-O(2) model in two dimensions that differs from the ordinary O(2) model by the addition of an explicit symmetry breaking term $-h_q\cos(q\varphi)$. Its coupling $h_q$ allows to smoothly interpolate between the O(2) model ($h_q=0$) and a $q$-state clock model ($h_q\rightarrow\infty$). In the latter case, a $q$-state clock model can also be defined for noninteger values of $q$. Thus, such a limit can also be considered as an analytic continuation of an ordinary $q$-state clock model to noninteger $q$. In previous work, we established the phase diagram for noninteger $q$ in the infinite coupling limit ($h_q\rightarrow\infty$). We showed that there is a second-order phase transition at low temperature and a crossover at high temperature. In this work, we seek to establish the phase diagram at finite values of the coupling using Monte Carlo and tensor methods. We show that for noninteger $q$, the second-order phase transition at low temperature and crossover at high temperature persist to finite coupling. For integer $q=2,3,4$, we know there is a second-order phase transition at infinite coupling (i.e. the well-known clock models). At finite coupling, we find that the critical exponents for $q=3,4$ vary with the coupling, and for $q=4$ the transition may turn into a BKT transition at small coupling. We comment on the similarities and differences of the phase diagrams with those of quantum simulators of the Abelian-Higgs model based on ladder-shaped arrays of Rydberg atoms. - Open Access
- 14 September 2021
- Published by American Physical Society (APS) in Physical Review D
- Vol. 104 (5) , 054505
Clock model interpolation and symmetry breaking in O(2) models
Motivated by recent attempts to quantum simulate lattice models with continuous Abelian symmetries using discrete approximations, we define an extended-O(2) model by adding a term to the ordinary O(2) model with angular values restricted to a interval. In the limit, the model becomes an extended -state clock model that reduces to the ordinary -state clock model when is an integer and otherwise is a continuation of the clock model for noninteger . By shifting the integration interval, the number of angles selected can change discontinuously and two cases need to be considered. What we call case 1 has one more angle than what we call case 2. We investigate this class of clock models in two space-time dimensions using Monte Carlo and tensor renormalization group methods. Both the specific heat and the magnetic susceptibility show a double-peak structure for fractional . In case 1, the small- peak is associated with a crossover, and the large- peak is associated with an Ising critical point, while both peaks are crossovers in case 2. When is close to an integer by an amount and the system is close to the small- Berezinskii-Kosterlitz-Thouless transition, the system has a magnetic susceptibility that scales as with estimates consistent with the magnetic critical exponent . The crossover peak and the Ising critical point move to Berezinskii-Kosterlitz-Thouless transition points with the same power-law scaling. A phase diagram for this model in the plane is sketched. These results are possibly relevant for configurable Rydberg-atom arrays where the interpolations among phases with discrete symmetries can be achieved by varying continuously the distances among atoms and the detuning frequency.- 14 September 2021
- Published by American Physical Society (APS) in Physical Review D
- Vol. 104 (5) , 054505
Open AccessClock model interpolation and symmetry breaking in O(2) models
Motivated by recent attempts to quantum simulate lattice models with continuous Abelian symmetries using discrete approximations, we define an extended-O(2) model by adding a term to the ordinary O(2) model with angular values restricted to a interval. In the limit, the model becomes an extended -state clock model that reduces to the ordinary -state clock model when is an integer and otherwise is a continuation of the clock model for noninteger . By shifting the integration interval, the number of angles selected can change discontinuously and two cases need to be considered. What we call case 1 has one more angle than what we call case 2. We investigate this class of clock models in two space-time dimensions using Monte Carlo and tensor renormalization group methods. Both the specific heat and the magnetic susceptibility show a double-peak structure for fractional . In case 1, the small- peak is associated with a crossover, and the large- peak is associated with an Ising critical point, while both peaks are crossovers in case 2. When is close to an integer by an amount and the system is close to the small- Berezinskii-Kosterlitz-Thouless transition, the system has a magnetic susceptibility that scales as with estimates consistent with the magnetic critical exponent . The crossover peak and the Ising critical point move to Berezinskii-Kosterlitz-Thouless transition points with the same power-law scaling. A phase diagram for this model in the plane is sketched. These results are possibly relevant for configurable Rydberg-atom arrays where the interpolations among phases with discrete symmetries can be achieved by varying continuously the distances among atoms and the detuning frequency.