Researcher profile

Hikaru Kawamura

Hikaru Kawamura contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
5topics
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

9 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&#39;/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&#39;/J1 <1, while in the uniform case of J1&#39;/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&#39; 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&#39; 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 &#34;asperities&#34; 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.