Source author record

Hikaru Kawamura

Hikaru Kawamura 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

46works
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

46 published item(s)

preprint2022arXiv

Emergent skyrmion-based chiral order in zero-field Heisenberg antiferromagnets on the breathing kagome lattice

We show that classical Heisenberg antiferromagnets on the breathing kagome lattice can be a platform to realize a zero-field topological order of the scalar spin chirality which can be viewed as a miniature skyrmion crystal (SkX) of discrete form with a small number of spins in its magnetic unit cell. In the model, a third nearest-neighbor (NN) antiferromagnetic interaction along the bond direction $J_3$ and the breathing bond-alternation characterized by the ratio of the NN interaction for large triangles to that for small ones, $J_1'/J_1$, are essential. It is found by means of Monte Carlo simulations that a commensurate triple-${\bf Q}$ state appearing for relatively strong $J_3$ at zero field is the noncoplanar state with the SkX structure in the breathing case of $J_1'/J_1 \neq 1$, while in the uniform case of $J_1'/J_1 =1$, it is a collinear state favored by thermal fluctuations. The origin of this chiral order and experimental implications of our result are also discussed.

preprint2022arXiv

Hedgehog lattice and field-induced chirality in breathing-pyrochlore Heisenberg antiferromagnets

We theoretically investigate a $J_1$-$J_3$ classical Heisenberg model on the breathing pyrochlore lattice, where the nearest-neighbor (NN) exchange interactions for small and large tetrahedra, $J_1$ and $J_1'$, take different values due to the breathing bond-alternation and $J_3$ is the third NN antiferromagnetic interaction along the bond direction. It is found by means of Monte Carlo simulations that for large $J_3$, a hedgehog lattice, a three-dimensional periodic array of magnetic monopoles and antimonopoles, emerges in the form of a quadruple-${\bf Q}$ state characterized by the ordering vector of ${\bf Q}=(\pm\frac{1}{2},\pm\frac{1}{2},\pm\frac{1}{2})$, being irrespective of the signs of $J_1$ and/or $J_1'$ as long as $J_1\neq J_1'$. It is also found that in an applied magnetic field, there appear six quadruple-${\bf Q}$ states depending on the values of $J_1$ and $J_1'$, among which three phases including the in-field hedgehog-lattice state exhibit nonzero total chirality $\mbox{\boldmath $χ$}^{\rm T}$ associated with the anomalous Hall effect of chirality origin. In the remaining two chiral phases, which are realized in the presence of ferromagnetic $J_1$ and/or $J_1'$, the spin structure is not topologically nontrivial, in spite of the fact that $\mbox{\boldmath $χ$}^{\rm T} \neq 0$. The role of the topological objects of the monopoles in $\mbox{\boldmath $χ$}^{\rm T}$ is also discussed.

preprint2022arXiv

Skyrmion crystal in the RKKY system on the two-dimensional triangular lattice

We study the ordering properties of the isotropic RKKY Heisenberg model on the two-dimensional (2D) triangular lattice by extensive Monte Carlo simulations to get insights into the chiral-degenerate skyrmion crystal (SkX) of metallic magnets. Our Hamiltonian contains only the spin-quadratic RKKY interaction derived from the spherical Fermi surface, containing neither the nesting nor the many-body interaction. The SkX phase is stabilized under applied fields where the frustration associated with the oscillating nature of the RKKY interaction and the emergent many-body interactions generated by thermal fluctuations play important roles. Replica symmetry breaking, reported in our recent study on the 3D RKKY model [Phys. Rev. B 104, 184432 (2021)], turns out to be absent in the present 2D model. Implications to the SkX formation mechanism are discussed.

preprint2020arXiv

Hedgehog-lattice spin texture in classical Heisenberg antiferromagnets on the breathing pyrochlore lattice

The hedgehog lattice, a three-dimensional periodic array of magnetic monopoles and antimonopoles, is known to be realized in the presence of the Dzyaloshinskii-Moriya (DM) interaction. Here, we demonstrate by means of Monte Carlo simulations that the hedgehog lattice is induced by not the DM interaction but frustration in classical Heisenberg antiferromagnets on the breathing pyrochlore lattice. In the model, the breathing bond-alternation is characterized by the ratio of the nearest-neighbor (NN) antiferromagnetic exchange interaction for large tetrahedra to that for small ones, J1'/J1. A quadruple-q state with the ordering vector of q=(\pm 1/2,\pm 1/2,\pm 1/2), which is realized for a large third NN antiferromagnetic interaction along the bond direction J3, turns out to become the hedgehog-lattice state in the breathing case of J1'/J1 <1, while in the uniform case of J1'/J1 =1, it is a collinear state favored by thermal fluctuations. It is also found that in a magnetic field, the structure of the (1/2,1/2,1/2) hedgehog lattice is changed from cubic to tetragonal, resulting in a nonzero net spin chirality which in a metalic system, should yield a characteristic topological Hall effect.

preprint2020arXiv

Monte Carlo studies of the spin-chirality decoupling in the three-dimensional Heisenberg spin glass

An extensive equilibrium Monte Carlo simulation is performed on the 3D isotropic Heisenberg SG model with the random nearest-neighbor Gaussian coupling, with particular interest in its chiral-glass (CG) and spin-glass (SG) orderings. For this model, the possibility of the spin-chirality decoupling, {\it i.e.\}, the CG order setting in at a higher temperature than that of the SG order was suggested earlier, but still remains controversial. We simulate the model up to the maximum size (linear dimension) $L=48$ under both periodic and open boundary conditions (BC). In locating the CG and SG transition temperatures $T_{\rm CG}$ and $T_{\rm SG}$ by the $L\rightarrow \infty$ extrapolation, a variety of independent physical quantities under the both BC are computed and utilized to get larger number of degrees of freedom (NDF). Thanks to the large NDF up to NDF=43, we succeed in obtaining stable and accurate estimates of the CG and SG transition temperatures, $T_{\rm CG}=0.142\pm 0.001$ and $T_{\rm SG}=0.131^{+0.001}_{-0.006}$. No sign of the size crossover is observed. For larger $L$, the CG correlation length progressively outgrows the SG correlation length at low temperatures. These results provide strong numerical support for the spin-chirality decoupling. The critical exponents associated with the CG and SG transitions are evaluated by use of the finite-size scaling with the scaling correction. For the CG transition, we get the CG exponents, $ν_{\rm CG}=1.36\pm 0.10$ and $η_{\rm CG}=0.49\pm 0.10$, consistently with the corresponding experimental exponents of canonical SG. Implications to the chirality scenario of experimental SG ordering is discussed.

preprint2019arXiv

Monte Carlo study of the critical properties of noncollinear Heisenberg magnets: $O(3)\times O(2)$ universality class

The critical properties of the antiferromagnetic Heisenberg model on the three-dimensional stacked-triangular lattice are studied by means of a large-scale Monte Carlo simulation in order to get insight into the controversial issue of the criticality of the noncollinear magnets with the $O(3)\times O(2)$ symmetry. The maximum size studied is $384^3$, considerably larger than the sizes studied by the previous numerical works on the model. Availability of such large-size data enables us to examine the detailed critical properties including the effect of corrections to the leading scaling. Strong numerical evidence of the continuous nature of the transition is obtained. Our data indicates the existence of significant corrections to the leading scaling. Careful analysis by taking account of the possible corrections yield critical exponents estimates, $α=0.44(3)$, $β=0.26(2)$, $γ=1.03(5)$, $ν=0.52(1)$, $η=0.02(5)$, and the chirality exponents $β_κ=0.40(3)$ and $γ_κ=0.77(6)$, supporting the existence of the $O(3)$ chiral (or $O(3)\times O(2)$) universality class governed by a new `chiral' fixed point. We also obtain an indication that the underlying fixed point is of the focus-type, characterized by the complex-valued correction-to-scaling exponent, $ω=0.1^{+0.4}_{-0.05} + i\ 0.7^{+0.1}_{-0.4}$. The focus-like nature of the chiral fixed point accompanied by the spiral-like renormalization-group (RG) flow is likely to be the origin of the apparently complicated critical behavior. The results are compared and discussed in conjunction with the results of other numerical simulations, several distinct types of RG calculations including the higher-order perturbative massive and massless RG calculations and the nonperturbative functional RG calculation, and the conformal-bootstrap program.

preprint2019arXiv

Nature of the randomness-induced quantum spin liquids in two dimensions

The nature of the randomness-induced quantum spin liquid state, the random-singlet state, is investigated in two dimensions (2D) by means of the exact-diagonalization and the Hams-de Raedt methods for several frustrated lattices, e.g., the triangular, the kagome and the J_1-J_2 square lattices. Properties of the ground state, the low-energy excitations and the finite-temperature thermodynamic quantities are investigated. The ground state and the low-lying excited states consist of nearly isolated singlet-dimers, clusters of resonating singlet-dimers, and orphan spins. Low-energy excitations are either singlet-to-triplet excitations, diffusion of orphan spins accompanied by the recombination of nearby singlet-dimers, creation or destruction of resonating singlet-dimers clusters. The latter two excitations give enhanced dynamical `liquid-like' features to the 2D random-singlet state. Comparison is made with the random-singlet state in a 1D chain without frustration, the similarity and the difference between in 1D and in 2D being highlighted. Frustration in a wide sense, not only the geometrical one but also including the one arising from the competition between distinct types of interactions, play an essential role in stabilizing this frustrated random singlet state. Recent experimental situations on both organic and inorganic materials are reviewed and discussed.

preprint2018arXiv

Nature of the high-speed rupture of the two-dimensional Burridge-Knopoff model of earthquakes

The nature of the high-speed rupture or the main shock of the Burridge-Knopoff spring-block model in two dimensions obeying the rate-and-state dependent friction law is studied by means of extensive computer simulations. It is found that the rupture propagation in larger events is highly anisotropic and irregular in shape on longer length scales, although the model is completely uniform and the rupture-propagation velocity is kept constant everywhere at the rupture front. The manner of the rupture propagation sometimes mimics the successive ruptures of neighboring "asperities" observed in real large earthquakes. Large events tend to be unilateral, with its epicenter lying at the rim of its rupture zone. The epicenter site is also located next to the rim of the rupture zone of some past event. Event-size distributions are computed and discussed in comparison with those of the corresponding one-dimensional model. The magnitude distribution exhibits a power-law behavior resembling the Gutenberg-Richter law for smaller magnitudes, which changes over to a more characteristic behavior for larger magnitudes. The behavior of the rupture length for larger events is discussed in terms of the strongly anisotropic rupture propagation of large events reflecting the underlying geometry.

preprint2016arXiv

Finite-temperature crossover phenomenon in the $S=1/2$ antiferromagnetic Heisenberg model on the kagome lattice

Thermal properties of the $S=1/2$ kagome Heisenberg antiferromagnet at low temperatures are investigated by means of the Hams-de Raedt method for clusters of up to 36 sites possessing a full symmetry of the lattice. The specific heat exhibits, in addition to the double peaks, the third and the forth peaks at lower temperatures. With decreasing the temperature, the type of the magnetic short-range order (SRO) changes around the third-peak temperature from the $\sqrt{3} \times \sqrt{3}$ to the $q$=0 states, suggesting that the third peak of the specific heat is associated with a crossover phenomenon between the spin-liquid states with distinct magnetic SRO. Experimental implications are discussed.

preprint2016arXiv

Spin-lattice-coupled order in Heisenberg antiferromagnets on the pyrochlore lattice

Effects of local lattice distortions on the spin ordering are investigated for the antiferromagnetic classical Heisenberg model on the pyrochlore lattice. It is found by Monte Carlo simulations that the spin-lattice coupling (SLC) originating from site phonons induces a first-order transition into two different types of collinear magnetic ordered states. The state realized at stronger SLC is cubic symmetric characterized by the magnetic (1/2,1/2,1/2) Bragg peaks, while that at weaker SLC is tetragonal symmetric characterized by the (1,1,0) ones, each accompanied by the commensurate local lattice distortions. Experimental implications to chromium spinels are discussed.

preprint2015arXiv

Dynamics of earthquake nucleation process represented by the Burridge-Knopoff model

Dynamics of earthquake nucleation process is studied on the basis of the one-dimensional Burridge-Knopoff (BK) model obeying the rate- and state-dependent friction (RSF) law. We investigate the properties of the model at each stage of the nucleation process, including the quasi-static initial phase, the unstable acceleration phase and the high-speed rupture phase or a mainshock. Two kinds of nucleation lengths L_sc and L_c are identified and investigated. The nucleation length L_sc and the initial phase exist only for a weak frictional instability regime, while the nucleation length L_c and the acceleration phase exist for both weak and strong instability regimes. Both L_sc and L_c are found to be determined by the model parameters, the frictional weakening parameter and the elastic stiffness parameter, hardly dependent on the size of an ensuing mainshock. The sliding velocity is extremely slow in the initial phase up to L_sc, of order the pulling speed of the plate, while it reaches a detectable level at a certain stage of the acceleration phase. The continuum limits of the results are discussed. The continuum limit of the BK model lies in the weak frictional instability regime so that a mature homogeneous fault under the RSF law always accompanies the quasi-static nucleation process. Duration times of each stage of the nucleation process are examined. The relation to the elastic continuum model and implications to real seismicity are discussed.

preprint2015arXiv

Low-temperature magnetic properties of the Kondo lattice model in one dimension

Low-temperature magnetic properties of the one-dimensional Kondo-lattice model with classical localized spins are investigated by means of a Monte Carlo simulation combined with an exact diagonalization technique. Comparative simulations are also made for the RKKY classical Heisenberg model known to be perturbatively derived in the weak-coupling limit. In addition to the previously identified antiferromagnetic, ferromagnetic and coplanar helical phases, the chiral noncoplanar phase and the collinear phases of the period-4 and 6 are newly identified, neither of which is realized in the corresponding RKKY model. The period-4 collinear phase persists even in the weak-coupling limit at the 1/4- and 3/4-fillings, invalidating the perturbative approach at these fillings.

preprint2015arXiv

Static and dynamical spin correlations of the S=1/2 random-bond antiferromagnetic Heisenberg model on the triangular and the kagome lattices

Inspired by the recent theoretical suggestion that the random-bond $S=1/2$ antiferromagnetic Heisenberg model on the triangular and the kagome lattices might exhibit a randomness-induced quantum spin liquid (QSL) behavior when the strength of the randomness exceeds a critical value, and that this "random-singlet state" might be relevant to the QSL behaviors experimentally observed in triangular organic salts $κ{\rm -(ET)_2 Cu_2 (CN)_3}$and ${\rm EtMe_3 Sb[Pd(dmit)_2]_2}$ and in kagome herbertsmithite ${\rm CuZn_3(OH)_6Cl_2}$, we further investigate the nature of the static and the dynamical spin correlations of these models. We compute the static and the dynamical spin structure factors, $S({\bf q})$ and $S({\bf q},ω)$, by means of an exact diagonalization method. In both triangular and kagome models, the computed $S({\bf q},ω)$ in the random-singlet state depends on the wavevector ${\bf q}$ only weakly, robustly exhibiting gapless behaviors accompnied by the broad distribution extending to higher energy $ω$. Especially in the strongly random kagome model, $S({\bf q},ω)$ hardly depends on ${\bf q}$, and exhibits an almost flat distribution for a wide range of $ω$, together with a $ω=0$ peak. These features agree semi-quantitatively with the recent neutron-scattering data on a single-crystal herbertsmithite, suggesting that the QSL state observed in herbersmithite might indeed be the randomness-induced QSL state, {\it i.e.\/}, the random-singlet state.

preprint2014arXiv

Quantum spin-liquid behavior in the spin-1/2 random Heisenberg antiferromagnet on the triangular lattice

Experimental quest for the hypothetical "quantum spin liquid" state has recently met a few promising candidate materials including organic salts κ-(ET)2Cu2(CN)3 and EtMe3Sb[Pd(dmit)2]2, S=1/2 triangular-lattice Heisenberg antiferromagnets consisting of molecular dimers. These compounds exhibit no magnetic ordering nor the spin freezing down to very low temperature, while various physical quantities exhibit gapless behaviors. Recent dielectric measurements revealed the glassy dielectric response suggesting the random freezing of the electric polarization degrees of freedom. Inspired by this observation, we propose as a minimal model of the observed quantum spin-liquid behavior the S=1/2 antiferromagnetic Heisenberg on the triangular lattice with a quenched randomness in the exchange interaction. We study both zero- and finite-temperature properties of the model by an exact diagonalization method, to find that when the randomness exceeds a critical value the model exhibits a quantum spin-liquid ground state. This randomness-induced quantum spin-liquid state exhibits gapless behaviors including the temperature-linear specific heat. The results provide a consistent explanation of the recent experimental observations on organic salts.

preprint2014arXiv

Quantum spin-liquid behavior in the spin-1/2 random-bond Heisenberg antiferromagnet on the kagome lattice

The effect of the quenched bond-randomness on the ordering of the $S=1/2$ antiferromagnetic Heisenberg model on the kagome lattice is investigated by means of an exact-diagonalization method. When the randomness exceeds a critical value, the ground state of the model exhibits a transition within the non-magnetic state into the randomness-relevant gapless spin-liquid state. Implications to the S=1/2 kagome-lattice antiferromagnet herbertsmithite is discussed.

preprint2013arXiv

Monte Carlo simulations of the three-dimensional XY spin glass focusing on the chiral and the spin order

The ordering of the three-dimensional isotropic {\it XY} spin glass with the nearest-neighbor random Gaussian coupling is studied by extensive Monte Carlo simulations. To investigate the ordering of the spin and the chirality, we compute several independent physical quantities including the glass order parameter, the Binder parameter, the correlation-length ratio, the overlap distribution and the non-self-averageness parameter, {\it etc}, for both the spin-glass (SG) and the chiral-glass (CG) degrees of freedom. Evidence of the spin-chirality decoupling, {\it i.e.}, the CG and the SG order occurring at two separated temperatures, $0<T_{SG}<T_{CG}$, is obtained from the glass order parameter, which is fully corroborated by the Binder parameter. By contrast, the CG correlation-length ratio yields a rather pathological and inconsistent result in the range of sizes we studied, which may originate from the finite-size effect associated with a significant short-length drop-off of the spatial CG correlations. Finite-size-scaling analysis yields the CG exponents $ν_{CG}=1.36^{+0.15}_{-0.37}$ and $η_{CG}=0.26^{+0.29}_{-0.26}$, and the SG exponents $ν_{SG}=1.22^{+0.26}_{-0.06}$ and $η_{SG}=-0.54^{+0.24}_{-0.52}$. The obtained exponents are close to those of the Heisenberg SG, but are largely different from those of the Ising SG. The chiral overlap distribution and the chiral Binder parameter exhibit the feature of a continuous one-step replica-symmetry breaking (1RSB), consistently with the previous reports. Such a 1RSB feature is again in common with that of the Heisenberg SG, but is different from the Ising one, which may be the cause of the difference in the CG critical properties from the Ising SG ones despite of a common $Z_2$ symmetry.

preprint2012arXiv

Finite-temperature transition of the distorted kagome-lattice Heisenberg antiferromagnet

Motivated by the recent experiment on kagome-lattice antiferromagnets, we study the zero-field ordering behavior of the antiferromagnetic classical Heisenberg model on a uniaxially distorted kagome lattice by Monte Carlo simulations. A first-order transition, which has no counterpart in the corresponding undistorted model, takes place at a very low temperature. Origin of the transition is ascribed to a cooperative proliferation of topological excitations inherent to the model.

preprint2012arXiv

Ordering of the Heisenberg spin glass in four dimensions

Ordering of the Heisenberg spin glass in four dimensions (4D) with the nearest-neighbor Gaussian coupling is investigated by equilibrium Monte Carlo simulations, with particular attention to its spin and chiral orderings. It is found that the spin and the chirality order simultaneously with a common correlation-length exponent $ν_{CG}=ν_{SG}\simeq 1.0$, i.e., the absence of the spin-chirality decoupling in 4D. Yet, the spin-glass ordered state exhibits a nontrivial phase-space structure associated with a continuous one-step-like replica-symmetry breaking, different in nature from that of the Ising spin glass or of the mean-field spin glass. Comparison is made with the ordering of the Heisenberg spin glass in 3D, and with that of the 1D Heisenberg spin glass with a long-range power-law interaction. It is argued that the 4D might be close to the marginal dimension separating the spin-chirality decoupling/coupling regimes.

preprint2012arXiv

Spin and chiral orderings of the antiferromagnetic XY model on the triangular lattice and their critical properties

We study the antiferromagnetic {\it XY} model on a triangular lattice by extensive Monte Carlo simulations, focusing on its ordering and critical properties. Our result clearly shows that two separate transitions occur at two distinct temperatures, the one at a higher temperature is associated with a $Z_2$-symmetry breaking driven by the chirality, and the one at a lower temperature is associated with the onset of the quasi-long-range order of the {\it XY} spin. We carefully examine the critical properties of each transition to find that the criticality of the chiral transition is consistent with the standard two-dimensional Ising universality class, whereas that of the spin transition might differ from the conventional Kosterlitz-Thouless (KT) one. The observed non-KT nature of the spin criticality is consistent with the most recent simulation result on the fully-frustrated {\it XY} model on a square lattice.

preprint2012arXiv

Statistical Physics of Fracture, Friction and Earthquake

We review the present status of our research and understanding regarding the dynamics and the statistical properties of earthquakes, mainly from a statistical physical viewpoint. Emphasis is put both on the physics of friction and fracture, which provides a "microscopic" basis of our understanding of an earthquake instability, and on the statistical physical modelling of earthquakes, which provides "macroscopic" aspects of such phenomena. Recent numerical results on several representative models are reviewed, with attention to both their "critical" and "characteristic" properties. We highlight some of relevant notions and related issues, including the origin of power-laws often observed in statistical properties of earthquakes, apparently contrasting features of characteristic earthquakes or asperities, the nature of precursory phenomena and nucleation processes, the origin of slow earthquakes, etc.

preprint2011arXiv

Cubic and non-cubic multiple-q states in the Heisenberg antiferromagnet on the pyrochlore lattice

The ordering of the classical Heisenberg model on the pyrochlore lattice with the antiferromagnetic nearest-neighbor interaction J1 and the ferromagnetic next-nearest-neighbour interaction J2 is investigated by means of a mean-field analysis and a Monte Carlo simulation. For a moderate J2/J1-value, the model exhibits a first-order transition into an incommensurate multiple-q ordered state where multiple Bragg peaks coexist in the spin structure factor. We show that there are two types of metastable multiple-q states, a cubic symmetric sextuple-q state and a non-cubic symmetric quadruple-q state. Based on a Monte Carlo simulation, we find that the cubic sextuple-q state appears just below the first-order transition temperature, while another transition from the cubic sextuple-q state to the non-cubic quadruple-q state occurs at a lower temperature.

preprint2011arXiv

Multiple-q states and skyrmion lattice of the triangular-lattice Heisenberg antiferromagnet under magnetic fields

Ordering of the frustrated classical Heisenberg model on the triangular-lattice with an incommensurate spiral spin structure is studied under magnetic fields by means of a mean-field analysis and a Monte Carlo simulation. Several types of multiple-q states including the "skyrmion-lattice" state is observed in addition to the standard single-q state. In contrast to the Dzyaloshinskii-Moriya interaction driven system, the present model allows both skyrmions and anti-skyrmions, together with a new thermodynamic phase where skyrmion and anti-skyrmion lattices form a domain state.

preprint2011arXiv

The ordering of XY spin glasses

Ordering properties of XY-like spin-glass magnets with an easy-plane magnetic anisotropy are studied based on a symmetry consideration and the results of recent numerical simulations on the pure Heisenberg and XY spin-glass models. The effects of an easy-plane-type uniaxial anisotropy, a random magnetic anisotropy and an applied magnetic field are investigated. In the XY regime in zero field, the "spin-chirality decoupling" persists even under the random magnetic anisotropy, escaping the "spin-chirality recoupling" phenomenon which inevitably occurs in the Heisenberg regime. Contrast between the scalar chiral order and the vector chiral order is emphasized. Implications to experiments are discussed.

preprint2011arXiv

Z2-vortex order of frustrated Heisenberg antiferromagnets in two dimensions

We discuss the recent experimental data on various frustrated quasi-two-dimensional Heisenberg antiferromagnets from the viewpoint of the Z2-vortex order, which include S=3/2 triangular-lattice antiferromagnet NaCrO2, S=1 triangular-lattice antiferromagnet NiGa2S4, S=1/2 organic triangular-lattice antiferromagnets κ-(BEDT-TTF)2Cu2(CN)3 and EtMe3Sb[Pd(dmit)2]2, and S=1/2 kagome-lattice antiferromagnet volborthite Cu3V2O7(OH)22H2O, etc.

preprint2010arXiv

Asperity characteristics of the Olami-Feder-Christensen model of earthquakes

Properties of the Olami-Feder-Christensen (OFC) model of earthquakes are studied by numerical simulations. The previous study indicated that the model exhibits ``asperity''-like phenomena, {\it i.e.}, the same region ruptures many times near periodically [T.Kotani {\it et al}, Phys. Rev. E {\bf 77}, 010102 (2008)]. Such periodic or characteristic features apparently coexist with power-law-like critical features, {\it e.g.}, the Gutenberg-Richter law observed in the size distribution. In order to clarify the origin and the nature of the asperity-like phenomena, we investigate here the properties of the OFC model with emphasis on its stress distribution. It is found that the asperity formation is accompanied by self-organization of the highly concentrated stress state. Such stress organization naturally provides the mechanism underlying our observation that a series of asperity events repeat with a common epicenter site and with a common period solely determined by the transmission parameter of the model. Asperity events tend to cluster both in time and in space.

preprint2010arXiv

Monte Carlo studies of the ordering of the one-dimensional Heisenberg spin glass with long-range power-law interactions

The nature of the ordering of the one-dimensional Heisenberg spin-glass model with a long-range power-law interaction is studied by extensive Monte Carlo simulations, with particular attention to the issue of the spin-chirality decoupling/coupling. Large system sizes up to $L=4096$ are studied. With varying the exponent $σ$ describing the power-law interaction, we observe three distinct types of ordering regimes. For smaller $σ$, the spin and the chirality order at a common finite temperature with a common correlation-length exponent, exhibiting the standard spin-chirality coupling behavior. For intermediate $σ$, the chirality orders at a temperature higher than the spin, exhibiting the spin-chirality decoupling behavior. For larger $σ$, both the spin and the chirality order at zero temperature. We construct a phase diagram in the $σ$ versus the temperature plane, and discuss implications of the results. Critical properties associated with both the chiral-glass and the spin-glass transitions are also determined.

preprint2010arXiv

Novel spin-liquid states in the frustrated Heisenberg antiferromagnet on the honeycomb lattice

Recent experiment on a honeycomb-lattice Heisenberg antiferromagnet (AF) Bi$_3$Mn$_4$O$_{12}$(NO$_3$) revealed a novel spin-liquid-like behavior down to low temperature, which was ascribed to the frustration effect due to the competition between the AF nearest- and next-nearest-neighbor interactions $J_1$ and $J_2$. Motivated by the experiment, we study the ordering of the $J_1$ -$J_2$ frustrated classical Heisenberg AF on a honeycomb lattice both by a low-temperature expansion and a Monte Carlo simulation. The model has been known to possess a massive degeneracy of the ground state, which, however, might be lifted due to thermal fluctuations leading to a unique ordered state, the effect known as 'order-by-disorder'. We find that the model exhibits an intriguing ordering behavior, particularly near the AF phase boundary. The energy scale of the order-by-disorder is suppressed there down to extremely low temperatures, giving rise to exotic spin-liquid states like a "ring-liquid" or a "pancake-liquid" state accompanied by the characteristic spin structure factor and the field-induced antiferromagnetism. We argue that the recent experimental data are explicable if the system is in such exotic spin-liquid states.

preprint2010arXiv

Phase transition of the three-dimensional chiral Ginzburg-Landau model -- search for the chiral phase

Nature of the phase transition of regularly frustrated vector spin systems in three dimensions is investigated based on a Ginzburg-Landau-type effective Hamiltonian. On the basis of the variational analysis of this model, Onoda et al recently suggested the possible occurrence of a chiral phase, where the vector chirality exhibits a long-range order without the long-range order of the spin [Phys. Rev. Lett. 99, 027206 (2007)]. In the present paper, we elaborate their analysis by considering the possibility of a first-order transition which was not taken into account in their analysis. We find that the first-order transition indeed occurs within the variational approximation, which significantly reduces the stability range of the chiral phase, while the chiral phase still persists in a restricted parameter range. Then, we perform an extensive Monte Carlo simulation focusing on such a parameter range. Contrary to the variational result, however, we do not find any evidence of the chiral phase. The range of the chiral phase, if any, is estimated to be less than 0.1% in the temperature width.

preprint2010arXiv

Signature of a Z_2 vortex in the dynamical correlations of the triangular-lattice Heisenberg antiferromagnet

Dynamical properties of the classical Heisenberg antiferromagnet on the triangular lattice are investigated by means of a spin-dynamics simulation and an analytical calculation. While the model was suggested to exhibit a topological transition driven by the Z_2-vortex binding-unbinding, weakness of the associated thermodynamic singularity has made it difficult to observe the evidence of the Z_2 vortex experimentally so far, only some indirect support. Here, we show that the signature of the Z_2-vortex excitation and the vortex-induced topological transition can be captured from the dynamical spin correlations. In particular, the dynamical spin structure factor exhibits a characteristic central peak around the Z_2-vortex transition, and this central peak will be a fingerprint of the Z_2-vortex excitation.

preprint2010arXiv

Simulation study of the inhomogeneous Olami-Feder-Christensen model of earthquakes

Statistical properties of the inhomogeneous version of the Olami-Feder-Christensen (OFC) model of earthquakes is investigated by numerical simulations. The spatial inhomogeneity is assumed to be dynamical. Critical features found in the original homogeneous OFC model, e.g., the Gutenberg-Richter law and the Omori law are often weakened or suppressed in the presence of inhomogeneity, whereas the characteristic features found in the original homogeneous OFC model, e.g., the near-periodic recurrence of large events and the asperity-like phenomena persist.

preprint2010arXiv

Spin-chirality decoupling and critical properties of a two-dimensional fully frustrated XY model

We study the ordering of the spin and the chirality in the fully frustrated XY model on a square lattice by extensive Monte Carlo simulations. Our results indicate unambiguously that the spin and the chirality exhibit separate phase transitions at two distinct temperatures, i. e. , the occurrence of the spin-chirality decoupling. The chirality exhibits a long-range order at T_c=0.45324(1) via a second-order phase transition, where the spin remains disordered with a finite correlation length ξ_s(T_c) \sim 120. The critical properties of the chiral transition determined from a finite-size scaling analysis for large enough systems of linear size L > ξ_s(T_c) are well compatible with the Ising universality. On the other hand, the spin exhibits a phase transition at a lower temperature T_s=0.4418(5) into the quasi-long-range-ordered phase. We found η(T_s)=0.201(1), suggesting that the universality of the spin transition is different from that of the conventional Kosterlitz-Thouless (KT) transition.

preprint2010arXiv

Spin-chirality decoupling in the one-dimensional Heisenberg spin glass with long-range power-law interactions

We study the issue of the spin-chirality decoupling/coupling in the ordering of the Heisenberg spin glass by performing large-scale Monte Carlo simulations on a one-dimensional Heisenberg spin-glass model with a long-range power-law interaction up to large system sizes. We find that the spin-chirality decoupling occurs for an intermediate range of the power-law exponent. Implications to the corresponding $d$-dimensional short-range model is discussed.

preprint2010arXiv

Two models of spin glasses -- Ising versus Heisenberg

Brief review is given on recent numerical research of the ordering of two typical models of spin glasses (SGs), the three-dimensional (3D) Ising SG and the 3D Heisenberg SG models. Particular attention is paid to the questions of whether there is a thermodynamic transition in zero field, what are the associated critical properties, what is the nature of the ordered state, particularly of a possible replica-symmetry breaking, and whether there is a thermodynamic transition in applied fields. The properties of the two models are contrasted, and possible relation to experiments is discussed.

preprint2010arXiv

Z_2-vortex ordering of the triangular-lattice Heisenberg antiferromagnet

Ordering of the classical Heisenberg antiferromagnet on the triangular lattice is studied by means of a mean-field calculation, a scaling argument and a Monte Carlo simulation, with special attention to its vortex degree of freedom. The model exhibits a thermodynamic transition driven by the Z_2-vortex binding-unbinding, at which various thermodynamic quantities exhibit an essential singularity. The low-temperature state is a "spin-gel" state with a long but finite spin correlation length where the ergodicity is broken topologically. Implications to recent experiments on triangular-lattice Heisenberg antiferromagnets are discussed.

preprint2009arXiv

Monte Carlo studies of the chiral and spin orderings of the three-dimensional Heisenberg spin glass

The nature of the ordering of the three-dimensional isotropic Heisenberg spin glass with nearest-neighbor random Gaussian coupling is studied by extensive Monte Carlo simulations. Several independent physical quantities are measured both for the spin and for the chirality, including the correlation-length ratio, the Binder ratio, the glass order parameter, the overlap distribution function and the non-self-averageness parameter. By controlling the effect of the correction-to-scaling, we have obtained a numerical evidence for the occurrence of successive chiral-glass and spin-glass transitions at nonzero temperatures, T_{CG} > T_{SG} > 0. Hence, the spin and the chirality are decoupled in the ordering of the model. The chiral-glass exponents are estimated to be ν_{CG}=1.4+-0.2 and η_{CG}=0.6+-0.2, indicating that the chiral-glass transition lies in a universality class different from that of the Ising spin glass. The possibility that the spin and chiral sectors undergo a simultaneous Kosterlitz-Thouless-type transition is ruled out. The chiral-glass state turns out to be non-self-averaging, possibly accompanying a one-step-like peculiar replica-symmetry breaking. Implications to the chirality scenario of experimental spin-glass transitions are discussed.

preprint2008arXiv

Simulation study of the two-dimensional Burridge-Knopoff model of earthquakes

Spatiotemporal correlations of the two-dimensional spring-block (Burridge-Knopoff) model of earthquakes are extensively studied by means of numerical computer simulations. The model is found to exhibit either ``subcritical'' or ``supercritical'' behavior, depending on the values of the model parameters. Transition between these regimes is either continuous or discontinuous. Seismic events in the ``subcritical'' regime and those in the ``supercritical'' regime at larger magnitudes exhibit universal scaling properties. In the ``supercritical'' regime, eminent spatiotemporal correlations, {\it e.g.}, remarkable growth of seismic activity preceding the mainshock, arise in earthquake occurrence, whereas such spatiotemporal correlations are significantly suppressed in the ``subcritical'' regime. Seismic activity is generically suppressed just before the mainshock in a close vicinity of the epicenter of the upcoming event while it remains to be active in the surroundings (the Mogi doughnut). It is also observed that, before and after the mainshock, the apparent $B$-value of the magnitude distribution decreases or increases in the ``supercritical'' or ``subcritical'' regimes, respectively. Such distinct precursory phenomena may open a way to the prediction of the upcoming large event.

preprint2006arXiv

Novel ordering of the pyrochlore Heisenberg antiferromagnet with the ferromagnetic next-nearest-neighbor interaction

The ordering property of the classical pyrochlore Heisenberg antiferromagnet with the ferromagnetic next-nearest-neighbor interaction is investigated by means of a Monte Carlo simulation. The model is found to exhibit a first-order transition at a finite temperature into a peculiar ordered state. While the spin structure factor, i.e., the thermal average of the squared Fourier amplitude of the spin, exhibits a finite long-range order characterized by the commensurate spin order of the period four, the thermal average of the spin itself almost vanishes. It means that, although the amplitude of the spin Fourier component is long-range ordered, the associated phase degree of freedom remains to be fluctuating.

preprint2006arXiv

Simulation Study of the One-Dimensional Burridge-Knopoff Model of Earthquakes

Spatio-temporal correlations of the one-dimensional spring-block (Burridge-Knopoff) model of earthquakes are extensively studied by means of numerical computer simulations. Particular attention is paid to clarifying how the statistical properties of earthquakes depend on the frictional and elastic properties of earthquake faults. It is found that, as the velocity-weakening tendency of the friction force gets weaker, the system tends to be more critical, while, as the velocity-weakening tendency gets stronger, the system tends to be more off-critical with enhanced features of a characteristic earthquake. The model exhibits several eminent precursory phenomena prior to the large event in its spatio-temporal correlations. Preceding the mainshock, the frequency of smaller events is gradually enhanced, whereas, just before the mainshock, it is suppressed in a close vicinity of the epicenter of the upcoming event (the Mogi doughnut). The time scale of the onset of the doughnut-like quiescence depends on the extent of the frictional instability. Under certain conditions, preceding the mainshock, the apparent B-value of the magnitude distribution increases significantly. The existence of such distinct precursory phenomena may open a way to the prediction of the time and the position of the upcoming large event.

preprint2006arXiv

Spatiotemporal Correlations of Earthquakes

Statistical properties of earthquakes are studied both by the analysis of real earthquake catalog of Japan and by numerical computer simulations of the spring-block model in both one and two dimensions. Particular attention is paid to the spatiotemporal correlations of earthquakes, e.g., the recurrence-time distribution or the time evolution before and after the mainshock of seismic distribution functions, including the magnitude distribution and the spatial seismic distribution. Certain eminent features of the spatiotemporal correlations, e.g., foreshocks, aftershocks, swarms and doughnut-like seismic pattern, are discussed.

preprint2006arXiv

Spin and chirality orderings of the one-dimensional Heisenberg spin glass with the long-range power-law interaction

The ordering of the one-dimensional Heisenberg spin glass interacting via the long-range power-law interaction is studied by Monte Carlo simulations. Particular attention is paid to the possible occurrence of the ``spin-chirality decoupling'' for appropriate values of the power-law exponent σ. Our result suggests that, for intermediate values of $σ$, the chiral-glass order occurs at finite temperatures while the standard spin-glass order occurs only at zero temperature.

preprint2002arXiv

Ordering of the Antiferromagnetic Heisenberg Model on a Pyrochlore Slab

Ordering of the geometrically frustrated two-dimensional Heisenberg antiferromagnet on a pyrochlore slab is studied by Monte Carlo simulations. The model is expected to serve as a reference system of SrCrGaO compound studied extensively. In sharp contrast to the kagomé Heisenberg antiferromagnet, the model exhibits locally non-coplanar spin structures at low temperatures, bearing nontrivial chiral degrees of freedom. We find that under certain conditions the model exhibits a novel Kosterlitz-Thouless-type transition at a finite temperature associated with these chiral degrees of freedom. Implications to experiments are discussed.

preprint2001arXiv

Spin and Chirality Orderings of Frustrated Magnets -- Stacked-Triangular Antiferromagnets and Spin Glasses

``Chirality'' is a multispin quantity representing the sense or the handedness of the noncollinear spin structures induced by spin frustration. Recent studies have revealed that the chirality often plays an important role in the ordering of certain frustrated magnets. Here I take up two such examples, stacked-triangular antiferromagnets and spin glasses, where the inherent chiral degree of freedom affects underlying physics and might lead to novel ordering phenomena. The first topic is the criticality of the magnetic phase transition of vector (i.e., XY or Heisenberg) antiferromagnets on the three-dimensional stacked-triangular lattice. The second topic is the nature of the spin-glass ordering. I will review the recent theoretical and experimental works on these topics, with particular emphasis on the important role played by the chirality.