Critical behavior of the pure and random-bond two-dimensional triangular Ising ferromagnet

Critical dynamics of cluster algorithms in the random-bond Ising model

In the present work, we present an extensive Monte Carlo simulation study on the dynamical properties of the two-dimensional random-bond Ising model. The correlation time τ of the Swendsen-Wang and Wolff cluster algorithms is calculated at the critical point. The dynamic critical exponent z of both algorithms is also measured by using the numerical data for several lattice sizes up to L=512. It is found for both algorithms that the autocorrelation time decreases considerably and the critical slowing-down effect reduces upon the introduction of bond disorder. Additionally, simulations with the Metropolis algorithm are performed, and the critical slowing-down effect is observed to be more pronounced in the presence of disorder, confirming the previous findings in the literature. Moreover, the existence of the non-self-averaging property of the model is demonstrated by calculating the scaled form of the standard deviation of autocorrelation times. Finally, the critical exponent ratio of the magnetic susceptibility is estimated by using the average cluster size of the Wolff algorithm.

Dynamic phase transitions on the kagome Ising ferromagnet

We perform extensive Monte Carlo simulations to investigate the dynamic phase transition properties of the two-dimensional kinetic Ising model on the kagome lattice in the presence of square-wave oscillating magnetic field. Through detailed finite-size scaling analysis, we study universality aspects of the nonequilibrium phase transition. Obtained critical exponents indicate that the two-dimensional kagome-lattice kinetic Ising model belongs to the same universality class with the corresponding Ising model in equilibrium. Moreover, dynamic critical exponent of the local moves used in simulations is determined with high precision. Our numerical results are compatible with the previous ones on kinetic Ising models.

First-principles calculation of the configurational energy density of states for a solid-state ion conductor with a variant of the Wang and Landau algorithm

In this work, a variant of the Wang and Landau algorithm for calculation of the configurational energy density of states is proposed. The algorithm was developed for the purpose of using first-principles simulations, such as density functional theory, to calculate the partition function of disordered sublattices in crystal materials. The expensive calculations of first-principles methods make a parallel algorithm necessary for a practical computation of the configurational energy density of states within a supercell approximation of a solid-state material. The algorithm developed in this work is tested with the two-dimensional (2d) Ising model to bench mark the algorithm and to help provide insight for implementation to a materials science application. Tests with the 2d Ising model revealed that the algorithm has good performance compared to the original Wang and Landau algorithm and the 1/t algorithm, in particular the short iteration performance. A proof of convergence is presented within an adiabatic assumption, and the analysis is able to correctly predict the time dependence of the modification factor to the density of states. The algorithm was then applied to the lithium and lanthanum sublattice of the solid-state lithium ion conductor Li_{0.5}La_{0.5}TiO_{3}. This was done to help understand the disordered nature of the lithium and lanthanum. The results find, overall, that the algorithm performs very well for the 2d Ising model and that the results for Li_{0.5}La_{0.5}TiO_{3} are consistent with experiment while providing additional insight into the lithium and lanthanum ordering in the material. The primary result is that the lithium and lanthanum become more mixed between layers along the c axis for increasing temperature. In part, the simulation of the disordered Li_{0.5}La_{0.5}TiO_{3} system serves as a benchmark for what size systems are currently and in the near future practical to calculate with density functional theory methods.

Violation of the Harris-Barghathi-Vojta Criterion

In 1974, Harris proposed his celebrated criterion: Continuous phase transitions in d-dimensional systems are stable against quenched spatial randomness whenever dν>2, where ν is the clean critical exponent of the correlation length. Forty years later, motivated by violations of the Harris criterion for certain lattices such as Voronoi-Delaunay triangulations of random point clouds, Barghathi and Vojta put forth a modified criterion for topologically disordered systems: aν>1, where a is the disorder decay exponent, which measures how fast coordination number fluctuations decay with increasing length scale. Here we present a topologically disordered lattice with coordination number fluctuations that decay as slowly as those of conventional uncorrelated randomness, but for which the clean universal behavior is preserved, hence violating even the modified criterion.

Two-dimensional Ising model on random lattices with constant coordination number

We study the two-dimensional Ising model on networks with quenched topological (connectivity) disorder. In particular, we construct random lattices of constant coordination number and perform large-scale Monte Carlo simulations in order to obtain critical exponents using finite-size scaling relations. We find disorder-dependent effective critical exponents, similar to diluted models, showing thus no clear universal behavior. Considering the very recent results for the two-dimensional Ising model on proximity graphs and the coordination number correlation analysis suggested by Barghathi and Vojta [Phys. Rev. Lett. 113, 120602 (2014)PRLTAO0031-900710.1103/PhysRevLett.113.120602], our results indicate that the planarity and connectedness of the lattice play an important role on deciding whether the phase transition is stable against quenched topological disorder.

Self-averaging in the random two-dimensional Ising ferromagnet

We study sample-to-sample fluctuations in a critical two-dimensional Ising model with quenched random ferromagnetic couplings. Using replica calculations in the renormalization group framework we derive explicit expressions for the probability distribution function of the critical internal energy and for the specific heat fluctuations. It is shown that the disorder distribution of internal energies is Gaussian, and the typical sample-to-sample fluctuations as well as the average value scale with the system size L like ∼Llnln(L). In contrast, the specific heat is shown to be self-averaging with a distribution function that tends to a δ peak in the thermodynamic limit L→∞. While previously a lack of self-averaging was found for the free energy, we here obtain results for quantities that are directly measurable in simulations, and implications for measurements in the actual lattice system are discussed.

Breathing modes of Kolumbo submarine volcano (Santorini, Greece)

Submarine volcanoes, such as Kolumbo (Santorini, Greece) are natural laboratories for fostering multidisciplinary studies. Their investigation requires the most innovative marine technology together with advanced data analysis. Conductivity and temperature of seawater were recorded directly above Kolumbo’s hydrothermal vent system. The respective time series have been analyzed in terms of non–equilibrium techniques. The energy dissipation of the volcanic activity is monitored by the temperature variations of seawater. The venting dynamics of chemical products is monitored by water conductivity. The analysis of the time series in terms of stochastic processes delivers scaling exponents with turning points between consecutive regimes for both conductivity and temperature. Changes of conductivity are shown to behave as a universal multifractal and their variance is subdiffusive as the scaling exponents indicate. Temperature is constant over volcanic rest periods and a universal multifractal behavior describes its changes in line with a subdiffusive character otherwise. The universal multifractal description illustrates the presence of non–conservative conductivity and temperature fields showing that the system never retains a real equilibrium state. The existence of a repeated pattern of the combined effect of both seawater and volcanic activity is predicted. The findings can shed light on the dynamics of chemical products emitted from the vents and point to the presence of underlying mechanisms that govern potentially hazardous, underwater volcanic environments.

Nonconvergence of the Wang-Landau algorithms with multiple random walkers

This paper discusses some convergence properties in the entropic sampling Monte Carlo methods with multiple random walkers, particularly in the Wang-Landau (WL) and 1/t algorithms. The classical algorithms are modified by the use of m-independent random walkers in the energy landscape to calculate the density of states (DOS). The Ising model is used to show the convergence properties in the calculation of the DOS, as well as the critical temperature, while the calculation of the number π by multiple dimensional integration is used in the continuum approximation. In each case, the error is obtained separately for each walker at a fixed time, t; then, the average over m walkers is performed. It is observed that the error goes as 1/sqrt[m]. However, if the number of walkers increases above a certain critical value m>m_{x}, the error reaches a constant value (i.e., it saturates). This occurs for both algorithms; however, it is shown that for a given system, the 1/t algorithm is more efficient and accurate than the similar version of the WL algorithm. It follows that it makes no sense to increase the number of walkers above a critical value m_{x}, since it does not reduce the error in the calculation. Therefore, the number of walkers does not guarantee convergence.

Universality in disordered systems: the case of the three-dimensional random-bond Ising model

We study the critical behavior of the d=3 Ising model with bond randomness through extensive Monte Carlo simulations and finite-size scaling techniques. Our results indicate that the critical behavior of the random-bond model is governed by the same universality class as the site- and bond-diluted models, clearly distinct from that of the pure model, thus providing a complete set of universality in disordered systems.