Short-time critical dynamics of damage spreading in the two-dimensional Ising model

ISING MODEL OF CHORIOCAPILLARIS FLOW

PURPOSE

To develop a mathematical model of local blood flow in the choriocapillaris using an Ising model.

METHODS

A JavaScript Ising model was used to create images that emulated the development of signal voids as would be seen in optical coherence tomography angiography of the choriocapillaris. The model was produced by holding the temperature near criticality and varying the field strength. Individual frames were evaluated, and a movie video was created to show the hypothetical development of flow-related signal voids over a lifetime.

RESULTS

Much the same as actual choriocapillaris images in humans, the model of flow-related signal voids followed a power-law distribution. The slope and intercept both decreased with age, as is seen in human subjects.

CONCLUSION

This model is a working hypothesis, and as such can help predict system characteristics, evaluate conclusions drawn from studies, suggest new research questions, and provide a way of obtaining an estimate of behavior in which experimental data are not yet available. It may be possible to understand choriocapillaris blood flow in health and disease states by determining by observing deviations from an expected model.

Generalized Metropolis dynamics with a generalized master equation: an approach for time-independent and time-dependent Monte Carlo simulations of generalized spin systems

The extension of Boltzmann-Gibbs thermostatistics, proposed by Tsallis, introduces an additional parameter q to the inverse temperature β. Here, we show that a previously introduced generalized Metropolis dynamics to evolve spin models is not local and does not obey the detailed energy balance. In this dynamics, locality is only retrieved for q=1, which corresponds to the standard Metropolis algorithm. Nonlocality implies very time-consuming computer calculations, since the energy of the whole system must be reevaluated when a single spin is flipped. To circumvent this costly calculation, we propose a generalized master equation, which gives rise to a local generalized Metropolis dynamics that obeys the detailed energy balance. To compare the different critical values obtained with other generalized dynamics, we perform Monte Carlo simulations in equilibrium for the Ising model. By using short-time nonequilibrium numerical simulations, we also calculate for this model the critical temperature and the static and dynamical critical exponents as functions of q. Even for q≠1, we show that suitable time-evolving power laws can be found for each initial condition. Our numerical experiments corroborate the literature results when we use nonlocal dynamics, showing that short-time parameter determination works also in this case. However, the dynamics governed by the new master equation leads to different results for critical temperatures and also the critical exponents affecting universality classes. We further propose a simple algorithm to optimize modeling the time evolution with a power law, considering in a log-log plot two successive refinements.