Researcher profile

M. Žukovič

M. Žukovič contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
32works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

32 published item(s)

preprint2022arXiv

Absence of long-range order in a three-dimensional stacked Ising antiferromagnet on kagome lattice

We study the possibility of long-range ordering (LRO) in a 3D system of vertically stacked layers of Ising antiferromagnet on a kagome lattice (SIAKL). Monte Carlo simulations are carried out for a varying interlayer coupling strength and a finite-size scaling analysis is performed for selected cases. Unlike in the related Ising system on a triangular lattice, in which even a finite number of the layers can stabilize LRO, or the Heisenberg system on the same kagome lattice, in which LRO emerged in 3D for a sufficient strength of the interlayer coupling, no LRO could be observed in the present model. This makes SIAKL a rare example of a 3D Ising paramagnet composed of frustrated layers coupled by unfrustrated interaction.

preprint2022arXiv

Road to zero-field antiferromagnetic skyrmions in a frustrated AFM/FM heterostructure

We demonstrate a mechanism of significant reduction, including complete elimination, of the external magnetic field required for the stabilization of a skyrmion lattice (SkX) phase in a frustrated triangular Heisenberg antiferromagnet (AFM) with the Dzyaloshinskii-Moriya interaction. It is achieved by coupling of such a AFM plane to a reference ferromagnetic (FM) layer, which generates an effective field cooperating with the external magnetic field. If the FM layer shows some axial single-ion anisotropy then the effective field can also be generated in zero external field due to a spontaneous FM long-range ordering. Then a sufficiently large interlayer coupling can fully substitute the external magnetic field and the SkX phase in the AFM layer can be stabilized even in zero external field.

preprint2020arXiv

Absence of long-range order in a general spin-$S$ kagome lattice Ising antiferromagnet

The possibility of the emergence of some kind of long-range ordering (LRO) due to the increase of multiplicity of the local degrees of freedom (spin value $S$) is studied in an Ising antiferromagnet on a kagome lattice (IAKL) by Monte Carlo simulation. In particular, the critical exponent of the spin correlation function, obtained from a finite-size scaling analysis, is evaluated for various values of $S$, including $S=\infty$, with the goal to determine whether there exists some threshold value of the spin $S_C$ above which the system would show true or quasi-LRO, similar to a related model on a triangular lattice (IATL). It is found that, unlike in the IATL case, the IAKL model remains disordered for any spin value and any finite temperature.

preprint2020arXiv

New ordered phase in geometrically frustrated generalized $XY$ model

Critical properties of a geometrically frustrated generalized $XY$ model with antiferromagnetic (AFM) and third-order antinematic (AN3) couplings on a triangular lattice are studied by Monte Carlo simulation. It is found that 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 antinematic (AN2) coupling, besides the AFM and AN3 phases which appear in the limits of relatively strong AFM and AN3 interactions, respectively, it includes an additional complex canted antiferromagnetic (CAFM) phase. It emerges at lower temperatures, wedged between the AFM and AN3 phases, as a result of the competition between the AFM and AN3 couplings, which is absent in the model with the AN2 coupling. The AFM-CAFM and AN3-CAFM phase transitions are concluded to belong to the weak Ising and weak three-state Potts universality classes, respectively. Additionally, all three QLRO phases also feature true LRO of the standard and generalized chiralities, which both vanish simultaneously at second-order phase transitions with non-Ising critical exponents and the critical temperatures slightly higher than the magnetic and nematic order-disorder transition temperatures.

preprint2019arXiv

Emergence of a Skyrmion Phase in a Frustrated Heisenberg Antiferromagnet with Dzyaloshinskii-Moriya Interaction

We study formation of a periodical structure of quasiparticle-like magnetic configurations with non-trivial topological charge, known as a magnetic skyrmion phase, on a frustrated triangular lattice antiferromagnetic Heisenberg model with Dzyaloshinskii-Moriya interaction (DMI) by means of Monte Carlo simulations. The existence of such a phase, formed on the three interpenetrating sublattices, has been demonstrated for a sufficiently large strength of DMI. The goal of the present investigation is to establish the minimum values of both the DM interaction strength as well as the magnetic field intensity at which such a skyrmion phase can exist at very low temperatures (close to zero). We find that the skyrmion phase emerges above the DMI parameter value of $D = 0.2$ and persists within the field intensity $2 \lesssim h < 4.5$.

preprint2019arXiv

Magnetism of the diluted Ising antiferromagnet in a magnetic field on the Kagome lattice: single-spin cluster approximation

An effective-field theory based on the single-spin cluster has been used to study of a diluted spin-$1/2$ Ising antiferromagnet on the Kagome lattice with nearest-neighbor interactions. We observe five plateaus in the magnetization curve of the diluted antiferromagnet when a magnetic field is applied which is in agreement with Monte Carlo calculation. The effect of the site-dilution on the magnetic susceptibility is also investigated and discussed. In particular, we have found that the frustrated Kagome lattice inverse susceptibility fall to zero at $0$K.

preprint2019arXiv

Magnetocaloric properties of an Ising antiferromagnet on a kagome lattice

Owing to a high degree of geometrical frustration an Ising antiferromagnet on a kagome lattice (IAKL) is known to exhibit no long-range ordering at any temperature, including the ground state. Nevertheless, at low temperatures it shows a strongly correlated, highly fluctuating regime known as a cooperative paramagnet or classical spin liquid. In the ground state it is characterized by a macroscopic degeneracy which translates to a relatively large value of the residual entropy. It has been shown that the presence of a macroscopic degeneracy associated with geometrical frustration below the saturation field can facilitate an enhanced magnetocaloric effect (MCE), which can exceed that of an ideal paramagnet with equivalent spin by more than an order of magnitude. In the present study we investigate magnetic and magnetocaloric properties of IAKL by Monte Carlo simulation. In particular, we calculate the entropy of the system using the thermodynamic integration method and evaluate quantities which characterize MCE, such as the isothermal entropy and adiabatic temperature changes in a varying magnetic field. It is found that IAKL shows the most interesting magnetocaloric properties at low temperatures and moderate magnetic fields, suggesting its potential to be used in technological applications for low-temperature magnetic refrigeration.

preprint2019arXiv

Spin-glass-like ordering in a frustrated $J_1-J_2$ Ising antiferromagnet on a honeycomb lattice

We study the nature of a low-temperature phase in the frustrated honeycomb-lattice Ising model with first- and second-neighbor antiferromagnetic (AF) interactions, $J_1$ and $J_2$, respectively, for $R = J_2/J_1 > 1/4$. It is known that for $R < 1/4$ there is a phase transition at low temperatures to the AF phase. Nevertheless, little is known about the critical behavior of the model for $R > 1/4$, except for recent effective field results which detected no phase transition down to zero temperature. Our Monte Carlo results suggest that for $R > 1/4$ there is at least one peculiar phase transition accompanied by a spin-glass-like freezing to a highly degenerate state consisting of frozen domains with stripe-type AF ordering separated by zero-energy domain walls. In spite of the local ordering within the respective domains there is no ordering among them and thus, unlike in the corresponding square-lattice model with $R > 1/2$, there is no conventional magnetic long-range ordering spanning the entire system.

preprint2015arXiv

Frustrated mixed spin-1/2 and spin-1 Ising ferrimagnets on a triangular lattice

Mixed spin-1/2 and spin-1 Ising ferrimagnets on a triangular lattice with sublattices A, B and C are studied for two spin value distributions $(S_{\rm A},S_{\rm B},S_{\rm C})=(1/2,1/2,1)$ and $(1/2,1,1)$ by Monte Carlo simulations. The non-bipartite character of the lattice induces geometrical frustration in both systems, which leads to the critical behavior rather different from their ferromagnetic counterparts. We confirm second-order phase transitions belonging to the standard Ising universality class occurring at higher temperatures, however, in both models these change at tricritical points (TCP) to first-order transitions at lower temperatures. In the model $(1/2,1/2,1)$, TCP occurs on the boundary between paramagnetic and ferrimagnetic $(\pm 1/2,\pm 1/2,\mp 1)$ phases. The boundary between two ferrimagnetic phases $(\pm 1/2,\pm 1/2,\mp 1)$ and $(\pm 1/2,\mp 1/2,0)$ at lower temperatures is always first order and it is joined by a line of second-order phase transitions between the paramagnetic and the ferrimagnetic $(\pm 1/2,\mp 1/2,0)$ phases at a critical endpoint. The tricritical behavior is also confirmed in the model $(1/2,1,1)$ on the boundary between the paramagnetic and ferrimagnetic $(0,\pm 1,\mp 1)$ phases.

preprint2015arXiv

Phase transitions in a frustrated Ising antiferromagnet on a square lattice

The phase transitions occurring in the frustrated Ising square antiferromagnet with first- ($J_1 < 0$) and second- ($J_2 < 0$) nearest-neighbor interactions are studied within the framework of the effective-field theory with correlations based on the different cluster sizes and for a wide range of $R = J_2/|J_1|$. Despite the simplicity of the model, it has proved difficult to precisely determine the order of the phase transitions. In contrast to the previous effective-field study, we have found a first-order transition line in the region close to $R = -0.5$ not only between the superantiferromagnetic and paramagnetic ($R<-0.5$) but also between antiferromagnetic and paramagnetic ($R>-0.5$) phases.

preprint2015arXiv

Short-range correlations in modified planar rotator model

We introduce a model inspired from statistical physics that is shown to display flexible short-range spatial correlations which are potentially useful in geostatistical modeling. In particular, we consider a suitably modified planar rotator or XY model, traditionally used for modeling continuous spin systems in magnetism, and we demonstrate that it can capture spatial correlations typically present in geostatistical data. The empirical study of the spin configurations produced by Monte Carlo simulations at various temperatures and stages in the nonequilibrium regime shows that their spatial variability can be modeled by the flexible class of Matérn covariance functions. The correlation range and the smoothness of these functions vary significantly in the parameter space that consists of the temperature and the simulation time. We briefly discuss the potential of the model for efficient and automatic prediction of spatial data with short-range correlations, such as commonly encountered in geophysical and environmental applications.

preprint2015arXiv

Tricritical Properties of Antiferromagnetic Ising Model on the Square Lattice

The Ising square lattice model with nearest-neighbor (nn) interactions ($J_1$) is one of the few exactly solvable models [1]. Adding next-neareast- neighbor (nnn) interactions ($J_2$) or a magnetic field (or both) leads to the non solvability of the model and only some approximate solutions are possible. In this brief report we will review some results obtained within effective field theory. We will show that besides second-order transitions there are also lines of first-order transitions and the coordinates of tricritical points are calculated.

preprint2014arXiv

Thermodynamic description of the Ising antiferromagnet on triangular lattice with selective dilution by modified Pair Approximation method

The Pair Approximation method is modified in order to describe the systems with geometrical frustration. The Ising antiferromagnet on triangular lattice with selective dilution (Kaya-Berker model) is considered and a self-consistent thermodynamic description of this model is obtained. For this purpose, the Gibbs free-energy as a function of temperature, concentration of magnetic atoms on the selected sublattice, and external magnetic field is derived. In particular, the phase diagram is constructed and a comparison of different methods is presented. The thermodynamic quantities are discussed in the context of their physical validity and the improvement in the description introduced by modified method is emphasized.

preprint2013arXiv

Chiral universality class behavior of a non-chiral antiferroquadrupole system

The planar Heisenberg system with antiferroquadrupolar exchange on 3D stacked triangular lattice is shown to belong to the new &#34;chiral&#34; universality class, predicted for chiral antiferromagnets. The present system, however, displays no such chirality, which, according to the currently widely accepted concept, is supposed to be a key ingredient for the new critical behavior. Our claim that the new universality class should not be limited only to the chiral antiferromagnets is based on a simple mapping between the chiral antiferromagnetic and antiferroquadrupolar systems and supported by actual calculations of the transition temperatures and critical indices from finite-size scaling analysis of data issued from Histogram Monte Carlo simulations. In order to demonstrate the existence of two different universality classes in the behavior of quadrupolar systems, the finite-size scaling is also performed for the system with ferroquadrupolar exchange, which, in contrast to the antiferroquadrupolar exchange system case, produced standard critical behavior, as could also be anticipated from the mapping.

preprint2013arXiv

Effect of biquadratic exchange on phase transitions of a planar classical Heisenberg ferromagnet

Effect of biquadratic exchange on phase transitions of a planar classical Heisenberg (or XY) ferromagnet on a stacked triangular lattice is investigated by Standard Monte Carlo and Histogram Monte Carlo simulations in the region of a bilinear to biquadratic exchange interaction ratio $J_{1}/J_{2} \leq 1$. The biquadratic exchange is found to cause separate second-order phase transitions in a strong biquadratic exchange limit, followed by simultaneous dipole and quadrupole ordering, which is of first order for an intermediate range of the exchange ratio and changes to a second-order one again as $J_{1}/J_{2}$ is further increased. Thus, a phase diagram featuring both triple and tricritical points is obtained. Furthermore, a finite-size scaling analysis is used to calculate the critical indices for both dipole and quadrupole kinds of ordering.

preprint2013arXiv

Entropy of spin clusters with frustrated geometry

Geometrically frustrated clusters of Ising spins of different shapes on a triangular lattice are studied by exact enumeration and Monte Carlo simulation. The focus is laid on the ground-state energy and residual entropy behaviors as functions of the cluster shape and size, as well as the spin value. Depending on the cluster shape, the residual entropy density on approach to the thermodynamic limit can either vanish or remain finite and the dependence can be decreasing, increasing or non-monotonic. Nevertheless, the relative entropies normalized by the respective thermodynamic limit values turn out to be little sensitive to the spin value. Some interesting finite-temperature properties of the clusters are also discussed and attention is drawn to their magnetocaloric properties.

preprint2013arXiv

Environmental Time Series Interpolation Based on Spartan Random Processes

In many environmental applications, time series are either incomplete or irregularly spaced. We investigate the application of the Spartan random process to missing data prediction. We employ a novel modified method of moments (MMoM) for parameter inference. The CPU time of MMoM is shown to be much faster than that of maximum likelihood estimation and almost independent of the data size. We formulate an explicit Spartan interpolator for estimating missing data. The model validation is performed on both synthetic data and real time series of atmospheric aerosol concentrations. The prediction performance is shown to be comparable with that attained by the best linear unbiased (Kolmogorov-Wiener) predictor at reduced computational cost.

preprint2013arXiv

Histogram Monte Carlo Simulation of the Geometrically Frustrated XY Antiferromagnet with Biquadratic Exchange

Histogram Monte Carlo simulation is used to investigate effects of biquadratic exchange $J_{2}$ on phase transitions of a 3D classical XY antiferromagnet with frustration induced by the antiferromagnetic exchange $J_{1}$ and the stacked triangular lattice geometry. The biquadratic exchange is considered negative (antiferroquadrupolar) within the triangular planes and positive (ferroquadrupolar) between the planes. The phase diagram obtained features a variety of interesting phenomena arising from the presence of both the biquadratic exchange and frustration. In a strong biquadratic exchange limit ($|J_{1}|/|J_{2}| \leq 0.25$), the antiferroquadrupolar phase transition which is of second order is followed by the antiferromagnetic one which can be either first or second order. The separate antiferroquadrupolar and antiferromagnetic second-order transitions are found to belong to the chiral XY and Ising universality classes, respectively. If the biquadratic exchange is reduced both transitions are found to be first order and occur simultaneously in a wide region of $|J_{1}|/|J_{2}|$. However, if $|J_{2}| \rightarrow 0$ the transition changes to the second-order one with the chiral universality class critical behavior.

preprint2013arXiv

Magnetocaloric properties of a frustrated Blume-Capel antiferromagnet

Low-temperature magnetization processes and magnetocaloric properties of a geometrically frustrated spin-1 Blume-Capel model on a triangular lattice are studied by Monte Carlo simulations. The model is found to display qualitatively different behavior depending on the sign of the single-ion anisotropy D. For positive values of D we observe two magnetization plateaus, similar to the spin-1/2 Ising antiferromagnet, and negative isothermal entropy changes for any field intensity. For a range of small negative values of D there are four magnetization plateaus and the entropy changes can be either negative or positive, depending on the field. If D is negative but large in absolute value then the entropy changes are solely positive.

preprint2013arXiv

Monte Carlo investigation of the tricritical point stability in a three-dimensional Ising metamagnet

We use Monte Carlo simulations to study multicritical properties of an Ising metamagnet in an external field. According to the mean field theory predictions, a three-dimensional layered metamagnet is expected to display a tricritical point decomposition to a critical endpoint and a bicritical endpoint, when a ratio between intralayer ferromagnetic and interlayer antiferromagnetic couplings becomes sufficiently small. Our simulations show no evidence of such a decomposition and produce a tricritical behaviour even for a coupling ratio as small as R=0.01.

preprint2013arXiv

Monte Carlo simulation of the three-dimensional XY model with bilinear-biquadratic exchange interaction

The three-dimensional XY model with bilinear-biquadratic exchange interactions $J$ and $J&#39;$, respectively, has been studied by Monte Carlo simulations. From the detailed analysis of the thermal variation of various physical quantities, as well as the order parameter and energy histogram analysis, the phase diagram including two different ordered phases has been determined. There is a single phase boundary from a paramagnetic to a dipole-quadrupole ordered phase, which is of second order in a high $J/J&#39;$ ratio region, changing to a first-order one for $0.35 \leq J/J&#39; \leq 0.5$. Below $J/J&#39;=0.35$ there are two separate transitions: the first one to the quadrupole long-range order (QLRO) phase at higher temperatures, followed by another one to the dipole-quadrupole long-range order (DLRO) phase at lower temperatures. The finite-size scaling analysis yields values of the critical exponents for both the DLRO and QLRO transitions close to the values for the conventional XY model which includes no biquadratic exchange.

preprint2013arXiv

Phase boundaries of a spin-3/2 Blume-Emery-Griffiths model on a honeycomb lattice

The spin-3/2 Blume-Emery-Griffiths model on a honeycomb lattice is studied by Monte Carlo simulations with the goal to determine phase diagrams for a range of the model parameters and to investigate the nature of the phase transitions between the respective phases. For positive values of the biquadratic to bilinear interaction ratio $α$, we find two ferromagnetically ordered phases, F1 and F2, with the sublattice magnetizations (1/2,1/2) and (3/2,3/2), respectively, and our results confirm the discontinuous character of the order-disorder critical line as a function of the single-ion anisotropy strength, predicted by the effective-field theory (EFT). For negative values of $α$, there is another ferrimagnetic (FRM) phase of the type (1/2,3/2), located between F1 and F2. However, the step-like variation of the order-disorder critical frontier obtained from EFT for large negative $α$ is not reproduced and thus deemed artifact of the EFT approximation. Finite-size scaling analysis performed at various points between the respective identified phases gave the ratio of critical exponents $γ/ν$ consistent with the 2D Ising universality class, except in the vicinity of the boundary intersection, where the results deviated from the standard values beyond the measurement errors.

preprint2013arXiv

Residual entropy of spin-s triangular Ising antiferromagnet

We employ a thermodynamic integration method (TIM) to establish the values of the residual entropy for the geometrically frustrated spin-s triangular Ising antiferromagnet, with the spin values s = 1/2, 1, 3/2, 2 and 5/2. The case of s = 1/2, for which the exact value is known, is used to assess the TIM performance. We also obtain an analytical formula for the lower bound in a general spin-s model and conjecture that it should reasonably approximate the true residual entropy for sufficiently large $s$. Implications of the present results in relation to reliability of the TIM as an indirect method for calculating global thermodynamic quantities, such as the free energy and the entropy, in similar systems involving frustration and/or higher spin values by standard Monte Carlo sampling are briefly discussed.

preprint2012arXiv

Critical and compensation phenomena in a mixed-spin ternary alloy: a Monte Carlo study

By means of standard and histogram Monte Carlo simulations, we investigate the critical and compensation behaviour of a ternary mixed spin alloy of the type $AB_pC_{1-p}$ on a cubic lattice. We focus on the case with the parameters corresponding to the Prussian blue analog $({\rm Ni}^{\rm II}_p {\rm Mn}^{\rm II}_{1-p})_{1.5}[{\rm Cr^{III}(CN)}_6] \cdot n {\rm H_2O}$ and confront our findings with those obtained by some approximative approaches and the experiments.

preprint2012arXiv

Critical and tricritical behavior of a selectively diluted triangular Ising antiferromagnet in a field

We study a geometrically frustrated triangular Ising antiferromagnet in an external magnetic field which is selectively diluted with nonmagnetic impurities employing an effective-field theory with correlations and Monte Carlo simulations. We focus on the frustration-relieving effects of such a selective dilution on the phase diagram and find that it can lead to rather intricate phase diagrams in the dilution-field parameters space. In particular, in a highly (weakly) diluted system the frustration is greatly (little) relieved and such a system is found to display only the second(first)-order phase transitions at any field. On the other hand, for a wide interval of intermediate dilution values the transition remains second order at low fields but it changes to first order at higher fields and the system displays a tricritical behavior. The existence of the first-order transition in the region of intermediate dilution and high fields is verified by Monte Carlo simulations.

preprint2012arXiv

Effects of selective dilution on phase diagram and ground-state magnetizations of an Ising antiferromagnet on triangular and honeycomb lattices

We employ an effective-field theory with correlations in order to study the phase diagram and ground-state magnetizations of a selectively diluted Ising antiferromagnet on triangular and honeycomb lattices. Dilution of different sublattices with generally unequal probabilities results in a rather intricate phase diagram in the sublattice dilution parameters space. In the case of the frustrated triangular lattice antiferromagnet the selective dilution affects the degree of frustration which can lead to some peculiar phenomena, such as reentrant behavior of long-range order or unsaturated sublattice magnetizations at zero temperature. The selectively diluted Ising antiferromagnet on the honeycomb lattice is obtained as a special case when one sublattice of the triangular lattice is completely removed by dilution.

preprint2012arXiv

Ground states of the frustrated Blume-Emery-Griffiths model in a field

Ground-state properties of the Blume-Emery-Griffiths model with antiferromagnetic nearest-neighbor interactions on a triangular lattice are investigated in the presence of an external magnetic field. In particular, we explore the model&#39;s parameter space and identify regions with different degenerate ground states that may give rise to different magnetic phases also at finite temperatures. We demonstrate the presence of such phases by Monte Carlo simulations of magnetization processes for selected values of parameters.

preprint2012arXiv

Low-temperature metastable states in a stacked triangular Ising antiferromagnet

We study low-temperature magnetization processes in a stacked triangular Ising antiferromagnet by Monte Carlo simulations. In increasing and decreasing magnetic fields we observe multiple steps and hysteresis corresponding to formation of different metastable states. Besides the equidistant threefold splitting of the 1/3 ferrimagnetic plateau, we additionally confirm a fourth plateau in the field-increasing branch and a sizable remanence when the field is decreased to zero. The newly observed plateau only appears at sufficiently low temperature and sufficiently large exchange interaction in the stacking direction. These observations reasonably reproduce low-temperatures measurements on the spin-chain compound $\rm{Ca}_3\rm{Co}_2\rm{O}_6$.

preprint2012arXiv

Phase Diagram of a Diluted Triangular Lattice Ising Antiferromagnet in a Field

Magnetization processes and phase transitions in a geometrically frustrated triangular lattice Ising antiferromagnet in the presence of an external magnetic field and a random site dilution are studied by the use of an effective-field theory with correlations. We find that the interplay between the applied field and the frustration-relieving dilution results in peculiar phase diagrams in the temperature-field-dilution parameter space.

preprint2012arXiv

Phase diagram of a mixed spin-1 and spin-3/2 Ising ferrimagnet

Critical and compensation properties of a mixed spin-1 and spin-3/2 Ising ferrimagnet on a square lattice are investigated by standard and histogram Monte Carlo simulations. The critical temperature is studied as a function of a single-ion anisotropy strength. The second order of the phase transition is established by finite-size scaling for the entire boundary. Some previously obtained results, such as a tricritical point, predicted by the mean field theory (MFT) and effective field theory (EFT), or a first-order transition line separating two different ordered phases, obtained by the cluster variational theory (CVT), are deemed artifacts of the respective approximations. So is a reentrant phenomenon produced by CVT. Nevertheless, the multicompensation behavior predicted by MFT and EFT was confirmed.

preprint2012arXiv

Phase transitions in a triangular Blume-Capel antiferromagnet

We study the critical behavior of a frustrated Blume-Capel (BC) antiferromagnet on a triangular lattice by Monte Carlo simulations. For a reduced single-ion anisotropy strength $-1.47 \lesssim d < 0$ we find two phase transitions. The low-temperature phase is characterized by the antiferromagnetic long-range ordering (LRO) on two sublattices with the third one remaining in a non-magnetic state. At higher temperatures there is a critical region of the Berezinskii-Kosterlitz-Thouless (BKT) type with a power-law decaying spin-correlation function. For $-1.5 \le d \lesssim -1.47$, there is only one phase transition from the LRO to the paramagnetic region and the transition is of first order. The presence of the BKT phase in the current frustrated BC model is a new feature not observed in its non-frustrated counterparts. The values of the decay exponent $η$ of the BKT phase corresponding to upper and lower temperatures appear to be consistent with the theoretical predictions for the six-state clock model.