Source author record

Milan Žukovič

Milan Žukovič appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2021arXiv

XY model with competing higher-order interactions

We study effects of competing pairwise higher-order interactions (HOI) with alternating signs and exponentially decreasing intensity on critical behavior of the XY model. It is found that critical properties of such a generalized model can be very different from the standard XY model and can strongly depend on whether the number of HOI terms is odd or even. Inclusion of any odd number of HOI terms results in two consecutive phase transitions to distinct ferromagnetic quasi-long-range order phases. Even number of HOI terms leads to two phase transitions only if the decay of the HOI intensities is relatively slow. Then the high-temperature transition to the ferromagnetic phase is followed by another transition to a peculiar competition-induced canted ferromagnetic phase. In the limit of an infinite number of HOI terms only one phase transition is confirmed, and under the conditions of fierce competition between the even and odd terms the transition temperature can be suppressed practically to zero.

preprint2019arXiv

Efficient and Scalable Approach to Equilibrium Conditional Simulation of Gibbs Markov Random Fields

We study the performance of an automated hybrid Monte Carlo (HMC) approach for conditional simulation of a recently proposed, single-parameter Gibbs Markov random field (Gibbs MRF). The MRF is based on a modified version of the planar rotator (MPR) model and is used for efficient gap filling in gridded data. HMC combines the deterministic over-relaxation method and the stochastic Metropolis update with dynamically adjusted restriction and performs automatic detection of the crossover to the targeted equilibrium state. We focus on the ability of the algorithm to efficiently drive the system to equilibrium at very low temperatures even with sparse conditioning data. These conditions are the most challenging computationally, requiring extremely long relaxation times if simulated by means of the standard Metropolis algorithm. We demonstrate that HMC has considerable benefits in terms of both computational efficiency and prediction performance of the MPR method.

preprint2019arXiv

Global Thermodynamic Properties of Complex Spin Systems Calculated from Density of States and Indirectly by Thermodynamic Integration Method

Evaluation of global thermodynamic properties, such as the entropy or the free energy, of complex systems featuring a high degree of frustration or disorder is often desirable. Nevertheless, they cannot be measured directly in standard Monte Carlo simulation. Therefore, they are either evaluated indirectly from the directly measured quantities, for example by the thermodynamic integration method (TIM), or by applying more sophisticated simulation methods, such as the Wang-Landau (WL) algorithm, which can directly sample density of states. In the present investigation we compare the performance of the WL and TIM methods in terms of calculation of the entropy of an Ising antiferromagnetic system on a Kagome lattice - a typical example of a complex spin system with high geometrical frustration resulting in a non-zero residual entropy the value of which is exactly known. It is found that in terms of accuracy the implementationally simpler TIM can deliver results comparable with the more involved WL method.

preprint2019arXiv

Numerical Precision Effects on GPU Simulation of Massive Spatial Data, Based on the Modified Planar Rotator Model

The present research builds on a recently proposed spatial prediction method for discretized two-dimensional data, based on a suitably modified planar rotator (MPR) spin model from statistical physics. This approach maps the measured data onto interacting spins and, exploiting spatial correlations between them, which are similar to those present in geostatistical data, predicts the data at unmeasured locations. Due to the short-range nature of the spin pair interactions in the MPR model, parallel implementation of the prediction algorithm on graphical processing units (GPUs) is a natural way of increasing its efficiency. In this work we study the effects of reduced computing precision as well as GPU-based hardware intrinsic functions on the speedup and accuracy of the MPR-based prediction and explore which aspects of the simulation can potentially benefit the most from the reduced precision. It is found that, particularly for massive data sets, a thoughtful precision setting of the GPU implementation can significantly increase the computational efficiency, while incurring little to no degradation in the prediction accuracy.

preprint2019arXiv

Phase Diagram of a Generalized XY Model with Geometrical Frustration

In the present study we investigate the effects of geometrical frustration on the XY model with antiferromagnetic (AFM) coupling on a triangular lattice, generalized by the inclusion of a third-order antinematic term (AN3). We demonstrate that at non-zero temperatures such a generalization leads to a phase diagram consisting of three different quasi-long-range ordered (QLRO) phases. Compared to the model with the second-order AN coupling (AN2), it includes besides the AFM and AN3 phases which appear in the limits of relatively strong AFM and AN3 interactions, respectively, an additional complex noncollinear QLRO phase at lower temperatures wedged between the AFM and AN3 phases. This new phase originates from the competition between the AFM and AN3 couplings, which is absent in the model with the AN2 coupling.

preprint2016arXiv

Frustrated ground states of a generalized XY model and their mapping to non-magnetic structural analogues

Ground state phases of a generalized XY model with magnetic and generalized nematic couplings on a non-bipartite triangular lattice are investigated in the exchange interactions parameter space. We demonstrate that the model displays a number of ordered and quasi-ordered phases as a result of geometrical frustration and/or competition between the magnetic and generalized nematic interactions. The nature and the extent of the respective phases depend on the parameter $q$, characterizing the higher-order harmonics term in the Hamiltonian. Motivated by the recent discovery of the experimental realization of the model with $q=2$ in the seemingly unrelated field of the systems chemistry [A.B. Cairns et al., Nature Chemistry 8, 442 (2016)], the model for $q>2$ is discussed in the context of the prediction of structural phases of a class of bimetalic cyanides based on a mapping between the two systems.

preprint2016arXiv

Thermodynamic properties of a Hubbard model on a cubic cluster - exact diagonalization study at quarter filling

We study the thermodynamics of a zero-dimensional, cubic cluster described with a Hubbard Hamiltonian, focusing our interest on the magnetic properties. The range in which the studied cluster is paramagnetic is considered. The results are obtained by means of exact numerical diagonalization. Such thermodynamic quantities as entropy, specific heat, magnetic susceptibility, spin-spin correlations and double occupancy are discussed. Particular emphasis is put on the behaviour of local maxima of specific heat and susceptibility, which are analysed in terms of Schottky anomalies.

preprint2015arXiv

Mixed spin-1/2 and spin-1 Ising ferromagnets on a triangular lattice

By Monte Carlo simulations we study critical properties of the mixed spin-1/2 and spin-1 Ising model on a triangular lattice, considering two different ways of the spin-value distributions on the three sublattices: $(1/2,1/2,1)$ and $(1/2,1,1)$. In the former case, we find standard Ising universality class phase transitions between paramagnetic and magnetic phases but no phase transitions between two distinct magnetic phases $(\pm 1/2,\pm 1/2,\pm 1)$ and $(\pm 1/2,\pm 1/2,0)$ at any finite temperature, except for some interesting non-critical anomalies displayed by response functions. On the other hand, the latter case turns out to be a rare (or perhaps the only) example of a two-dimensional mixed spin-1/2 and spin-1 Ising model on a standard lattice that displays tricritical behavior.

preprint2014arXiv

Incommensurate Short-Range Order in $S=1$ Triangular Lattice Ising Antiferromagnet

We study the $S = 1$ triangular lattice Ising antiferromagnet by Monte Carlo simulations. Frustrations between a major antiferromagnetic third-neighbor interaction $J_3$ and a minor ferromagnetic nearest-neighbor interaction $J_1$ cause incommensurate short-range orders at intermediate temperatures. At low temperatures (below $T/J_3 \lesssim 0.2$ for $J_1/J_3 = -1/3$), the system exhibits fourfold periodic ordered state. In the short-range order phase, the system shows a glassy two-step relaxation. We demonstrate that the features of the short-range order are attributed to the cooperation between the frustrations and the nonmagnetic Ising spin states which is a particular feature of the integer spin systems.

preprint2014arXiv

Thermodynamic and magnetocaloric properties of geometrically frustrated Ising nanoclusters

Thermodynamic and magnetocaloric properties of geometrically frustrated Ising spin clusters of selected shapes and sizes are studied by exact enumeration. In the ground state the magnetization and the entropy show step-wise variations with an applied magnetic field. The number of steps, their widths and heights depend on the cluster shape and size. While the character of the magnetization plateau heights is always increasing, the entropy is not necessarily decreasing function of the field, as one would expect. For selected clusters showing some interesting ground-state properties, the calculations are extended to finite temperatures by exact enumeration of densities of states in the energy-magnetization space. In zero field the focus is laid on a peculiar behavior of some thermodynamic quantities, such as the entropy, the specific heat and the magnetic susceptibility. In finite fields various thermodynamic functions are studied in the temperature-field parameter plane and particular attention is paid to the cases showing an enhanced magnetocaloric effect. The exact results on the finite clusters are compared with the thermodynamic limit behavior obtained from Monte Carlo simulations.

preprint2013arXiv

A Directional Gradient-Curvature Method for Gap Filling of Gridded Environmental Spatial Data with Potentially Anisotropic Correlations

We introduce the Directional Gradient-Curvature (DGC) method, a novel approach for filling gaps in gridded environmental data. DGC is based on an objective function that measures the distance between the directionally segregated normalized squared gradient and curvature energies of the sample and entire domain data. DGC employs data-conditioned simulations, which sample the local minima configuration space of the objective function instead of the full conditional probability density function. Anisotropy and non-stationarity can be captured by the local constraints and the direction-dependent global constraints. DGC is computationally efficient and requires minimal user input, making it suitable for automated processing of large (e.g., remotely sensed) spatial data sets. Various effects are investigated on synthetic data. The gap-filling performance of DGC is assessed in comparison with established classification and interpolation methods using synthetic and real satellite data, including a skewed distribution of daily column ozone values. It is shown that DGC is competitive in terms of cross validation performance.

preprint2013arXiv

Classification of missing values in spatial data using spin models

A problem of current interest is the estimation of spatially distributed processes at locations where measurements are missing. Linear interpolation methods rely on the Gaussian assumption, which is often unrealistic in practice, or normalizing transformations, which are successful only for mild deviations from the Gaussian behavior. We propose to address the problem of missing values estimation on two-dimensional grids by means of spatial classification methods based on spin (Ising, Potts, clock) models. The "spin" variables provide an interval discretization of the process values, and the spatial correlations are captured in terms of interactions between the spins. The spins at the unmeasured locations are classified by means of the "energy matching" principle: the correlation energy of the entire grid (including prediction sites) is estimated from the sample-based correlations. We investigate the performance of the spin classifiers in terms of computational speed, misclassification rate, class histogram and spatial correlations reproduction using simulated realizations of spatial random fields, real rainfall data, and a digital test image. We also compare the spin-based methods with standard classifiers such as the k-nearest neighbor, the fuzzy k-nearest neighbor, and the Support Vector Machine. We find that the spin-based classifiers provide competitive choices.

preprint2013arXiv

Dynamics of episodic transient correlations in currency exchange rate returns and their predictability

We study the dynamics of the linear and non-linear serial dependencies in financial time series in a rolling window framework. In particular, we focus on the detection of episodes of statistically significant two- and three-point correlations in the returns of several leading currency exchange rates that could offer some potential for their predictability. We employ a rolling window approach in order to capture the correlation dynamics for different window lengths and analyze the distributions of periods with statistically significant correlations. We find that for sufficiently large window lengths these distributions fit well to power-law behavior. We also measure the predictability itself by a hit rate, i.e. the rate of consistency between the signs of the actual returns and their predictions, obtained from a simple correlation-based predictor. It is found that during these relatively brief periods the returns are predictable to a certain degree and the predictability depends on the selection of the window length.

preprint2012arXiv

Critical properties of a spin-1 triangular lattice Ising antiferromagnet

We employ Monte Carlo simulations in order to investigate critical behavior of a geometrically frustrated spin-1 Ising antiferromagnet on a triangular lattice in the presence of a single-ion anisotropy. It has been previously found that long-range order can exist in the isotropic system with a spin larger than some critical value estimated as 11/2. We show that the presence of the single-ion anisotropy can lead to a partial long-range order in the low-temperature region even below this critical value, namely for the spin 1, within a certain range of the anisotropy strength. At higher temperatures we identify another phase of the Berezinsky-Kosterlitz-Thouless type and using a finite-size scaling analysis evaluate the correlation decay exponent. We also study densities of various local spin patterns in the respective phases.

preprint2008arXiv

Multilevel Discretized Random Field Models with "Spin" Correlations for the Simulation of Environmental Spatial Data

A problem of practical significance is the analysis of large, spatially distributed data sets. The problem is more challenging for variables that follow non-Gaussian distributions. We show that the spatial correlations between variables can be captured by interactions between "spins". The spins represent multilevel discretizations of the initial field with respect to a number of pre-defined thresholds. The spatial dependence between the "spins" is imposed by means of short-range interactions. We present two approaches, inspired by the Ising and Potts models, that generate conditional simulations from samples with missing data. The simulations of the "spin system" are forced to respect locally the sample values and the system statistics globally. We compare the two approaches in terms of their ability to reproduce the sample statistical properties, to predict data at unsampled locations, as well as in terms of their computational complexity. We discuss the impact of relevant simulation parameters, such as the domain size, the number of discretization levels, and the initial conditions.