Reconfigurable computing for Monte Carlo simulations: results and prospects of the Janus project

We describe Janus, a massively parallel FPGA-based computer optimized for the simulation of spin glasses, theoretical models for the behavior of glassy materials. FPGAs (as compared to GPUs or many-core processors) provide a complementary approach to massively parallel computing. In particular, our model problem is formulated in terms of binary variables, and floating-point operations can be (almost) completely avoided. The FPGA architecture allows us to run many independent threads with almost no latencies in memory access, thus updating up to 1024 spins per cycle. We describe Janus in detail and we summarize the physics results obtained in four years of operation of this machine; we discuss two types of physics applications: long simulations on very large systems (which try to mimic and provide understanding about the experimental non equilibrium dynamics), and low-temperature equilibrium simulations using an artificial parallel tempering dynamics. The time scale of our non-equilibrium simulations spans eleven orders of magnitude (from picoseconds to a tenth of a second). On the other hand, our equilibrium simulations are unprecedented both because of the low temperatures reached and for the large systems that we have brought to equilibrium. A finite-time scaling ansatz emerges from the detailed comparison of the two sets of simulations. Janus has made it possible to perform spin glass simulations that would take several decades on more conventional architectures. The paper ends with an assessment of the potential of possible future versions of the Janus architecture, based on state-of-the-art technology.

Static versus dynamic heterogeneities in the D=3 Edwards-Anderson-ising spin glass

We numerically study the aging properties of the dynamical heterogeneities in the Ising spin glass. We find that a phase transition takes place during the aging process. Statics-dynamics correspondence implies that systems of finite size in equilibrium have static heterogeneities that obey finite-size scaling, thus signaling an analogous phase transition in the thermodynamical limit. We compute the critical exponents and the transition point in the equilibrium setting, and use them to show that aging in dynamic heterogeneities can be described by a finite-time scaling ansatz, with potential implications for experimental work.

Spin glasses on the hypercube

We present a mean field model for spin glasses with a natural notion of distance built in, namely, the Edwards-Anderson model on the diluted D-dimensional unit hypercube in the limit of large D. We show that finite D effects are strongly dependent on the connectivity, being much smaller for a fixed coordination number. We solve the nontrivial problem of generating these lattices. Afterward, we numerically study the nonequilibrium dynamics of the mean field spin glass. Our three main findings are the following: i the dynamics is ruled by an infinite number of time sectors, ii the aging dynamics consists of the growth of coherent domains with a nonvanishing surface-volume ratio, and iii the propagator in Fourier space follows the p4 law. We study as well the finite D effects in the nonequilibrium dynamics, finding that a naive finite size scaling ansatz works surprisingly well.

Microcanonical finite-size scaling in second-order phase transitions with diverging specific heat

A microcanonical finite-size ansatz in terms of quantities measurable in a finite lattice allows extending phenomenological renormalization the so-called quotients method to the microcanonical ensemble. The ansatz is tested numerically in two models where the canonical specific heat diverges at criticality, thus implying Fisher renormalization of the critical exponents: the three-dimensional ferromagnetic Ising model and the two-dimensional four-state Potts model (where large logarithmic corrections are known to occur in the canonical ensemble). A recently proposed microcanonical cluster method allows simulating systems as large as L = 1024 Potts or L= 128 (Ising). The quotients method provides accurate determinations of the anomalous dimension, η, and of the (Fisher-renormalized) thermal ν exponent. While in the Ising model the numerical agreement with our theoretical expectations is very good, in the Potts case, we need to carefully incorporate logarithmic corrections to the microcanonical ansatz in order to rationalize our data.

Cosmological perturbations in coherent oscillating scalar field models

The fact that fast oscillating homogeneous scalar fields behave as perfect fluids in average and their intrinsic isotropy have made these models very fruitful in cosmology. In this work we will analyse the perturbations dynamics in these theories assuming general power law potentials V(ϕ) = λ|ϕ|^n /n. At leading order in the wavenumber expansion, a simple expression for the effective sound speed of perturbations is obtained c_eff^ 2  = ω = (n − 2)/(n + 2) with ω the effective equation of state. We also obtain the first order correction in k^ 2/ω_eff^ 2 , when the wavenumber k of the perturbations is much smaller than the background oscillation frequency, ω_eff. For the standard massive case we have also analysed general anharmonic contributions to the effective sound speed. These results are reached through a perturbed version of the generalized virial theorem and also studying the exact system both in the super-Hubble limit, deriving the natural ansatz for δϕ; and for sub-Hubble modes, exploiting Floquet’s theorem.

Integral elastic, electronic-state, ionization, and total cross sections for electron scattering with furfural

We report absolute experimental integral cross sections (ICSs) for electron impact excitation of bands of electronic-states in furfural, for incident electron energies in the range 20-250 eV. Wherever possible, those results are compared to corresponding excitation cross sections in the structurally similar species furan, as previously reported by da Costa et al. [Phys. Rev. A 85, 062706 (2012)] and Regeta and Allan [Phys. Rev. A 91, 012707 (2015)]. Generally, very good agreement is found. In addition, ICSs calculated with our independent atom model (IAM) with screening corrected additivity rule (SCAR) formalism, extended to account for interference (I) terms that arise due to the multi-centre nature of the scattering problem, are also reported. The sum of those ICSs gives the IAM-SCAR+I total cross section for electron-furfural scattering. Where possible, those calculated IAM-SCAR+I ICS results are compared against corresponding results from the present measurements with an acceptable level of accord being obtained. Similarly, but only for the band I and band II excited electronic states, we also present results from our Schwinger multichannel method with pseudopotentials calculations. Those results are found to be in good qualitative accord with the present experimental ICSs. Finally, with a view to assembling a complete cross section data base for furfural, some binary-encounter-Bethe-level total ionization cross sections for this collision system are presented. (C) 2016 AIP Publishing LLC.

Variational mean-field approach to the double-exchange model

It has been recently shown that the double exchange Hamiltonian, with weak antiferromagnetic interactions, has a richer variety of first- and second-order transitions than previously anticipated, and that such transitions are consistent with the magnetic properties of manganites. Here we present a thorough discussion of the variational mean-field approach that leads to these results. We also show that the effect of the Berry phase turns out to be crucial to produce first-order paramagnetic-ferromagnetic transitions near half filling with transition temperatures compatible with the experimental situation. The computation relies on two crucial facts: the use of a mean-field ansatz that retains the complexity of a system of electrons with off-diagonal disorder, not fully taken into account by the mean-field techniques, and the small but significant antiferromagnetic superexchange interaction between the localized spins.