We develop a generalized pseudo-fermion functional renormalization group (PFFRG) approach that can be applied to arbitrary Heisenberg models with spins ranging from the quantum case S = 1/2 to the classical limit S → ∞. Within this framework, spins of magnitude S are realized by implementing M = 2S copies of spin-1/2 degrees of freedom on each lattice site. We confirm that even without explicitly projecting onto the highest spin sector of the Hilbert space, ground states tend to select the largest possible local spin magnitude. This justifies the average treatment of the pseudo fermion constraint in previous spin-1/2 PFFRG studies. We apply this method to the antiferromagnetic J1-J2 honeycomb Heisenberg model with nearest neighbor J1 > 0 and second neighbor J2 > 0 interactions. Mapping out the phase diagram in the J2/J1-S plane we find that upon increasing S quantum fluctuations are rapidly decreasing. In particular, already at S = 1 we find no indication for a magnetically disordered phase. In the limit S → ∞, the known phase diagram of the classical system is exactly reproduced. More generally, we prove that for S → ∞ the PFFRG approach is identical to the Luttinger-Tisza method.

We study the three-dimensional Kasteleyn transition in both nearest neighbours and dipolar spin ice models using an algorithm that conserves the number of excitations. We first limit the interactions range to nearest neighbours to test the method in the presence of a field applied along [Formula: see text], and then focus on the dipolar spin ice model. The effect of dipolar interactions, which is known to be greatly self screened at zero field, is particularly strong near full polarization. It shifts the Kasteleyn transition to lower temperatures, which decreases  ≈0.4 K for the parameters corresponding to the best known spin ice materials, [Formula: see text] and [Formula: see text]. This shift implies effective dipolar fields as big as 0.05 T opposing the applied field, and thus favouring the creation of 'strings' of reversed spins. We compare the reduction in the transition temperature with results in previous experiments, and study the phenomenon quantitatively using a simple molecular field approach. Finally, we relate the presence of the effective residual field to the appearance of string-ordered phases at low fields and temperatures, and we check numerically that for fields applied along [Formula: see text] there are only three different stable phases at zero temperature.

The dynamical structure factor is one of the experimental quantities crucial in scrutinizing the validity of the microscopic description of strongly correlated systems. However, despite its long-standing importance, it is exceedingly difficult in generic cases to numerically calculate it, ensuring that the necessary approximations involved yield a correct result. Acknowledging this practical difficulty, we discuss in what way results on the hardness of classically tracking time evolution under local Hamiltonians are precisely inherited by dynamical structure factors and, hence, offer in the same way the potential computational capabilities that dynamical quantum simulators do: We argue that practically accessible variants of the dynamical structure factors are bounded-error quantum polynomial time (BQP)-hard for general local Hamiltonians. Complementing these conceptual insights, we improve upon a novel, readily available measurement setup allowing for the determination of the dynamical structure factor in different architectures, including arrays of ultra-cold atoms, trapped ions, Rydberg atoms, and superconducting qubits. Our results suggest that quantum simulations employing near-term noisy intermediate-scale quantum devices should allow for the observation of features of dynamical structure factors of correlated quantum matter in the presence of experimental imperfections, for larger system sizes than what is achievable by classical simulation.