Critical magnetic flux for Weyl points in the three-dimensional Hofstadter model

We investigate the band structure of the three-dimensional Hofstadter model on cubic lattices, with an isotropic magnetic field oriented along the diagonal of the cube with flux Φ =2 π m /n , where m ,n are coprime integers. Using reduced exact diagonalization in momentum space, we show that, at fixed m , there exists an integer n (m ) associated with a specific value of the magnetic flux, that we denote by Φ_c(m ) ≡2 π m /n (m ) , separating two different regimes. The first one, for fluxes Φ <Φ_c(m ) , is characterized by complete band overlaps, while the second one, for Φ >Φ_c(m ) , features isolated band-touching points in the density of states and Weyl points between the m th and the (m +1 ) -th bands. In the Hasegawa gauge, the minimum of the (m +1 ) -th band abruptly moves at the critical flux Φ_c(m ) from k_z=0 to k_z=π . We then argue that the limit for large m of Φ_c(m ) exists and it is finite: limm_→∞Φ_c(m ) ≡Φ_c . Our estimate is Φ_c/2 π =0.1296 (1 ) . Based on the values of n (m ) determined for integers m ≤60 , we propose a mathematical conjecture for the form of Φ_c(m ) to be used in the large-m limit. The asymptotic critical flux obtained using this conjecture is Φ_c^{(conj )}/2 π =7 /54 .