Reciprocal-space correlation functions of nearest-neighbour spin ice

<p>Reciprocal-space correlation functions&nbsp;<span class="math-tex">\(\langle \sigma_\mu(k)\sigma_\nu(-k)\rangle\)</span> of the nearest-neighbour spin ice Hamiltonian <span class="math-tex">\(H = J \sum_{\langle ij\rangle} \sigma_i\sigma_j\)</span> on the pyrochlore lattice, from Monte Carlo simulations on 128&times;128&times;128 cubic unit cells.</p> <p><strong>File formats</strong></p> <p>The file name indicates the temperature (in units of J) where the measurement was taken.</p> <p>Files with extension .dat contain the correlation functions as double-precision complex floats.<br> Files with extension .err contain the standard errors as single-precision real floats.</p> <p>Both file types contain an array of shape (4, 4, 128, 256, 256). The meaning of the indices (from major to minor):</p> <ul> <li>sublattice index &mu;</li> <li>sublattice index &nu;</li> <li>wave vector components in units of <span class="math-tex">\(2\pi/(128a_0)\)</span></li> </ul> <p>The sublattice indices correspond to the site offsets <span class="math-tex">\(r_0 = [111]/8, r_1=[1\bar1\bar1]/8, r_2=[\bar11\bar1]/8,r_3=[\bar1\bar1 1]/8\)</span> from the nearest FCC lattice point. For the Fourier transforms, all sublattices are shifted to these lattice points, so the correlators stored in the files are periodic with respect to the FCC reciprocal lattice.</p> <p>The wave vector range covered is <span class="math-tex">\(0\le k_x&lt; 2\pi/a_0, 0\le k_y,k_z&lt; 4\pi/a_0\)</span>, an (unconventional) reciprocal-space unit cell of the FCC pyrochlore lattice.</p> <p>The utility loader.py loads the files, shapes them in the correct array format, and extracts single k-points.</p> <p><strong>Details of data generation</strong></p> <p>We performed Monte Carlo simulations of nearest-neighbour spin ice using the efficient loop-string algorithm introduced in <a href="https://doi.org/10.1103/PhysRevB.90.220406">Phys. Rev. B 90, 220406(R)</a>. For nonzero temperatures, correlations between different strings were excluded, which effectively averages all spin configurations compatible with a given loop graph. We ran 32 independent Markov chains at each temperature point and obtained 4096 Monte Carlo samples in each. The reported standard errors are the error on the mean of the 32 Markov chains.</p>