LHLeon Hostetler

Leon Hostetler

*This is an algorithmically generated author record based on publicly available information.

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


Leon Hostetler has a total of 260 co-authors. Here are the top co-authors by number of shared publications:


Latest Publications

  • Top Cited Papers
    Open Access

    The anomalous magnetic moment of the muon in the Standard Model: an update

    by R. Aliberti, T. Aoyama, E. Balzani, A. Bashir, G. Benton, J. Bijnens, V. Biloshytskyi, T. Blum, D. Boito, M. Bruno, et al.


  • Open Access

    Hadronic Vacuum Polarization for the Muon g−2 from Lattice QCD: Long-Distance and Full Light-Quark Connected Contribution

    by Alexei Bazavov, Claude W. Bernard, David A. Clarke, Christine Davies, Carleton DeTar, Aida X. El-Khadra, Elvira Gámiz, Steven Gottlieb, Anthony V. Grebe, Leon Hostetler, et al.

    We 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

    Hadronic vacuum polarization for the muon g−2 from lattice QCD: Complete short and intermediate windows

    by Alexei Bazavov, David A. Clarke, Christine T. H. Davies, Carleton DeTar, Aida X. El-Khadra, Elvira Gámiz, Steven Gottlieb, Anthony V. Grebe, Leon Hostetler, William I. Jay, et al.

    We 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

    Efficient state preparation for the Schwinger model with a theta term

    by Alexei Bazavov, Brandon Henke, Leon Hostetler, Dean Lee, Huey-Wen Lin, Giovanni Pederiva, Andrea Shindler

    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

    Improving HISQ propagator solves using deflation

    by Leon Hostetler, M. A. Clark, Carleton DeTar, Steven Gottlieb, Evan Weinberg

    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

    Symmetry breaking in an extended O(2) model

    by Leon Hostetler, Ryo Sakai, Jin Zhang, Alexei Bazavov, Yannick Meurice

    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 −hqcos(qφ). Its coupling hq allows to smoothly interpolate between the O(2) model (hq=0) and a q-state clock model (hq→∞). 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 (hq→∞). 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 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

    by Leon Hostetler, Ryo Sakai, Jin Zhang, Alexei Bazavov, Yannick Meurice

    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

    Clock model interpolation and symmetry breaking in O(2) models

    by Leon Hostetler, Jin Zhang, Ryo Sakai, Judah Unmuth-Yockey, Alexei Bazavov, Yannick Meurice

    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 γcos(qφ) term to the ordinary O(2) model with angular values restricted to a 2π interval. In the γ→∞ limit, the model becomes an extended q-state clock model that reduces to the ordinary q-state clock model when q is an integer and otherwise is a continuation of the clock model for noninteger q. By shifting the 2π 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 q. 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 q is close to an integer by an amount Δq and the system is close to the small-β Berezinskii-Kosterlitz-Thouless transition, the system has a magnetic susceptibility that scales as ∼1/(Δq)1−1/δ′ with δ′ estimates consistent with the magnetic critical exponent δ=15. 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 (β,q) 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.