Magnetic systems offer themselves as the ideal benchmark for generic concepts pertaining to collective phenomena in nature. This is due in part to the availability of a large variety of diverse magnetic materials that can be chosen to approximate simple theoretical “toy models” of collective behavior and, in part, to their ease of study by a battery of experimental techniques.^{1}
The MARQOV project focuses on the development of an adaptable framework for the simulation and analysis of equilibrium lattice spin models and its optimization on massively parallel architectures. As this area of research, with origins in statistical physics, particularly aims at improving our understanding of critical phenomena ^{2}, there exists considerable overlap with other research fields, such as solidstate physics, material science, and highenergy physics. In general, phase transitions are characterized by the emergence of longrange collective behavior of the individual microscopic degrees of freedom. This gives rise to the important phenomenon of universality^{3}, according to which the character of a phase transition is solely determined by only few fundamental properties of the system (such as its global symmetries and interaction ranges). Theoretically, the emergence of universality can be motivated by renormalization group approaches^{4}, where at criticality the system forgets about its microscopic realization under repeated scale transformations.
Quenched spatial disorder, however, can present a relevant perturbation to the phase transition, leaving the character of the transition unchanged in some scenarios^{5}^{6}^{7}, but leading to entirely different critical behavior in others^{8}^{9}^{10}. Spatial disorder is not only of theoretical interest though, but in fact intrinsic to most real systems, arising from impurities, defects, and other inhomogeneities.
Figure 1. Voronoi graph generated from a random point cloud.
In principle, the role of uncorrelated randomness (such as in the form of diluted regular lattices) is successfully explained by Harris' seminal criterion^{11}, using a remarkably simple argument of competing fluctuation scales. The role of topological randomness, i.e., lattices with random connectivity (see Fig. 1 for an example), however, has yet to be fully understood. Recent contributions in this field^{12}^{13} still strongly rely on numerical simulations, as topological disorder has so far not been within the reach of analytical methods. Specifically, the MARQOV code was recently used for probing disordered geometries and determine whether the universal character of an Ising transition is broken with respect to the clean model without disorder^{14}^{15}^{16}. Finally, Isinglike systems are also studied on abstract topologies like scalefree^{17}^{18} and smallworld networks^{19}^{20} as well as in curved geometries such as the hyperbolic plane^{21}, see Fig. 2.
Figure 2. Regular heptagonal hyperbolic tessellation.
Apart from being used as a probe to investigate aspects of topology and disorder relevance, the Ising model, due to its conceptual simplicity, is widely applied in different areas of research. Besides uniaxial magnets, it describes, e.g., experimental liquidvapor transitions and binary mixtures, certain transitions in highenergy physics (see [^13] and references therein), voter dynamics^{22} and can even be found in fairly unexpected applications such as the dynamics of traffic signals^{23} or pistachio trees^{24}. Generalizing the Ising model to \( O(N) \) symmetric spin vectors, one obtains the XY model (\( N=2 \)), featuring a Kosterlitz–Thouless transition and the Heisenberg case (\( N=3 \)), which models certain real magnetic materials. Even higher symmetries (\(N=4,5\)) become, for instance, again relevant in particle physics and hightemperature superconductor transitions (see Ref. ^{25} for a comprehensive review).
This large scope of applications justifies the continued relevance of classical spin models through decades of research. However, a comprehensive software package for those models does not exist to this date. Filling this gap, MARQOV offers a software framework combining modern HPCoptimized numerical simulation techniques and a stateoftheart statistical analysis into a convenient workflow. Beyond that, MARQOV’s particular strengths include quenched disordered and frustrated systems, such as spin glasses, where parallelization and automatization aspects become important.

Bramwell, S. T., & Gingras, M. J. (2001). Spin ice state in frustrated magnetic pyrochlore materials. Science, 294(5546), 14951501. ↩︎

Täuber, U. C. (2014). Critical dynamics: a field theory approach to equilibrium and nonequilibrium scaling behavior. Cambridge University Press. ↩︎

Amit, D. J., & MartinMayor, V. (2005). Field Theory, the Renormalization Group, and Critical Phenomena: Graphs to Computers Third Edition. World Scientific Publishing Company. ↩︎

Cardy, J. (1996). Scaling and renormalization in statistical physics (Vol. 5). Cambridge university press. ↩︎

Lima, F. W. S., Moreira, J. E., Andrade Jr, J. S., & Costa, U. M. S. (2000). The ferromagnetic Ising model on a Voronoi–Delaunay lattice. Physica A: Statistical Mechanics and its Applications, 283(12), 100106. ↩︎

Janke, W., & Villanova, R. (2002). Ising model on threedimensional random lattices: A Monte Carlo study. Physical Review B, 66(13), 134208. ↩︎

GordilloGuerrero, A., & RuizLorenzo, J. J. (2007). Selfaveraging in the threedimensional site diluted Heisenberg model at the critical point. Journal of Statistical Mechanics: Theory and Experiment, 2007(06), P06014. ↩︎

Hasenbusch, M., Toldin, F. P., Pelissetto, A., & Vicari, E. (2007). The universality class of 3D sitediluted and bonddiluted Ising systems. Journal of Statistical Mechanics: Theory and Experiment, 2007(02), P02016. ↩︎

Vojta, T., & Lee, M. Y. (2006). Nonequilibrium phase transition on a randomly diluted lattice. Physical Review Letters, 96(3), 035701. ↩︎

Vojta, T., Farquhar, A., & Mast, J. (2009). Infiniterandomness critical point in the twodimensional disordered contact process. Physical Review E, 79(1), 011111. ↩︎

Harris, A. B. (1974). Effect of random defects on the critical behaviour of Ising models. Journal of Physics C: Solid State Physics, 7(9), 1671. ↩︎

Barghathi, H., & Vojta, T. (2014). Phase transitions on random lattices: How random is topological disorder?. Physical Review Letters, 113(12), 120602. ↩︎

Schrauth, M., Portela, J. S., & Goth, F. (2018). Violation of the HarrisBarghathiVojta Criterion. Physical Review Letters, 121(10), 100601. ↩︎

Schrauth, M., Richter, J. A., & Portela, J. S. (2018). Twodimensional Ising model on random lattices with constant coordination number. Physical Review E, 97(2), 022144. ↩︎

Schrauth, M., & Portela, J. S. (2019). Fast algorithm for topologically disordered lattices with constant coordination number. Physical Review Research, 1(3), 033061. ↩︎

Schrauth, M., & Portela, J. S. (2019). Universality of continuous phase transitions on random Voronoi graphs. Physical Review E, 100(6), 062118. ↩︎

Herrero, C. P. (2004). Ising model in scalefree networks: A monte carlo simulation. Physical Review E, 69(6), 067109. ↩︎

Lima, F. W. S., & Sumour, M. A. (2012). Ising model with spins S=1/2 and 1 on directed and undirected Erdös–Rényi random graphs. Physica A: Statistical Mechanics and its Applications, 391(4), 948953. ↩︎

Herrero, C. P. (2002). Ising model in smallworld networks. Physical Review E, 65(6), 066110. ↩︎

Lima, F. W. S., & Plascak, J. A. (2013). Critical behavior of the Ising and BlumeCapel models on directed twodimensional smallworld networks. The European Physical Journal B, 86(7), 300. ↩︎

Baek, S. K., Shima, H., & Kim, B. J. (2009). Curvatureinduced frustration in the X Y model on hyperbolic surfaces. Physical Review E, 79(6), 060106. ↩︎

Dornic, I., Chaté, H., Chave, J., & Hinrichsen, H. (2001). Critical coarsening without surface tension: The universality class of the voter model. Physical Review Letters, 87(4), 045701. ↩︎

Suzuki, H., Imura, J. I., & Aihara, K. (2013). Chaotic Isinglike dynamics in traffic signals. Scientific reports, 3, 1127. ↩︎

Noble, A. E., Rosenstock, T. S., Brown, P. H., Machta, J., & Hastings, A. (2018). Spatial patterns of tree yield explained by endogenous forces through a correspondence between the Ising model and ecology. Proceedings of the National Academy of Sciences, 115(8), 18251830. ↩︎

Pelissetto, A., & Vicari, E. (2002). Critical phenomena and renormalizationgroup theory. Physics Reports, 368(6), 549727. ↩︎