Showing 1–50 of 57 results for author: Hotta, C

Diffusive dynamics of fractionalized particles and the enhanced conductivity at the border of the neutral-ionic transition

Authors: Yuta Sakai, Chisa Hotta

Abstract: We study the diffusive dynamics of the classical one-dimensional lattice model of mobile particles featuring the incoherent metallic state of the organic TTF-CA found in the vicinity of the neutral-ionic(NI)-transition. The particles are strongly correlated and feel the alternating site potentials, exhibiting uniform ionic (I) Mott insulating and neutral band insulating (N) phases when the Coulomb… ▽ More We study the diffusive dynamics of the classical one-dimensional lattice model of mobile particles featuring the incoherent metallic state of the organic TTF-CA found in the vicinity of the neutral-ionic(NI)-transition. The particles are strongly correlated and feel the alternating site potentials, exhibiting uniform ionic (I) Mott insulating and neutral band insulating (N) phases when the Coulomb interaction and site potentials are large, respectively. We focus on the neutral-ionic domain walls (NIDW) activated by their competition, which is typically regarded as fractionalized particles. The finite temperature phase diagram reveals a thermodynamically stable NIDW phase not found in previous literature that emerges at the triple point, where the transition line splits into two first-order lines toward the critical endpoints. We analyze the long-time behavior of dynamics of the NIDWs and find that it captures the incoherent transport of the system. The conductivity derived from the diffusion constant and the number population of NIDWs shows a strong enhancement inside the NIDW phase, which can explain the large conductivity at the NI crossover recently found in TTF-CA. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: 10 pages, 6 figures

Cluster-projected matrix product state: framework for engineering exact quantum many-body ground states in one and two dimensions

Authors: Hidehiro Saito, Chisa Hotta

Abstract: We propose a framework to design concurrently a frustration-free quantum many-body Hamiltonian and its exact ground states in one and two dimensions using an elementary matrix product state (MPS) representation. Our approach strategically chooses a local cluster Hamiltonian, which is arranged to overlap with neighboring clusters on a designed lattice. Ensuring that there exists a state spanned on… ▽ More We propose a framework to design concurrently a frustration-free quantum many-body Hamiltonian and its exact ground states in one and two dimensions using an elementary matrix product state (MPS) representation. Our approach strategically chooses a local cluster Hamiltonian, which is arranged to overlap with neighboring clusters on a designed lattice. Ensuring that there exists a state spanned on the lattice that has its local submanifolds as the lowest-energy eigenstate of every cluster, we can construct the bulk Hamiltonian as the sum of the cluster Hamiltonians. The key to find such a solution is a systematic protocol, which projects out excited states on every cluster using MPS and effectively entangles the cluster states. The protocol offers several advantages, including the ability to achieve exact many-body ground-state solutions at nearly equal cost in one and two dimensions including those with gapless or long-range entangled ground states, flexibility in designing Hamiltonians unbiasedly across various forms of models, and numerically feasible validation through energy calculations. Our protocol offers exact ground state for any given-frustration free Hamiltonian, and enables the exploration of exact phase boundaries and the analysis of even a spatially nonuniform random system, providing platforms for quantum simulations and benchmarks. △ Less

Submitted 19 June, 2024; v1 submitted 18 June, 2024; originally announced June 2024.

Pinch-point spectral singularity from the interference of topological loop states

Authors: Masafumi Udagawa, Hiroki Nakai, Chisa Hotta

Abstract: Pinch point is a spectral discontinuity found in the neutron diffraction image of spin ice. Similar spectral singularity is commonly observed in a broad range of systems that have a close connection with flat bands. We focus on the electron flat band and its two topologically distinct classes of wavefunction: the compact localized state (CLS), and the non-contractible loop state (NLS). We establis… ▽ More Pinch point is a spectral discontinuity found in the neutron diffraction image of spin ice. Similar spectral singularity is commonly observed in a broad range of systems that have a close connection with flat bands. We focus on the electron flat band and its two topologically distinct classes of wavefunction: the compact localized state (CLS), and the non-contractible loop state (NLS). We establish their simple mathematical relationship, showing that different Bloch NLSs can be derived as momentum derivatives of a Bloch CLS, depending on the approaching direction toward the singular point. This CLS-NLS correspondence helps visualize the pinch point as an interference pattern among NLSs through a ``polarizer", which encodes the information about the location of singular momentum and the experimental techniques like spin-polarized photoemission spectroscopy. It helps extract topological information knit to microscopic electronic and magnetic structures. △ Less

Submitted 21 April, 2024; originally announced April 2024.

Comments: 11 pages, 4 figures

Phase diagram of the quantum spin-1/2 Heisenberg-$Γ$ model on a frustrated zigzag chain

Authors: Hidehiro Saito, Chisa Hotta

Abstract: We investigate the quantum spin-1/2 zigzag chain with frustrated $J_1$-$J_2$ Heisenberg interactions, incorporating additional off-diagonal exchange interactions known as the $Γ$ term, both with and without an applied magnetic field. Based on the density-matrix renormalization group calculation, we map out the ground state phase diagram that shows a variety of magnetic and nonmagnetic phases inclu… ▽ More We investigate the quantum spin-1/2 zigzag chain with frustrated $J_1$-$J_2$ Heisenberg interactions, incorporating additional off-diagonal exchange interactions known as the $Γ$ term, both with and without an applied magnetic field. Based on the density-matrix renormalization group calculation, we map out the ground state phase diagram that shows a variety of magnetic and nonmagnetic phases including multicritical points and several exactly solvable points. Upon introducing a finite $Γ$ term, we observe the persistent dimer singlet state of the $J_1$-$J_2$ Heisenberg model, sustaining a nonzero spin gap, while also giving rise to a gapless nonmagnetic excitation, manifesting in the substantial zero-energy peak in the nematic dynamical structure factor. This gapless peak-mode remaining almost as a fluctuation to the ground state, induces dilute but robust concentration of nematicity on top of singlets on dimers, which we call the nematic singlet-dimer phase. When the whole nematic excited mode condenses and replaces the singlet, the nematic-dimer phase transforms to the Ising-type ferromagnetic or antiferromagnetic long-range orders that arise from the $Γ$ term spontaneously selecting magnetic easy axes. Its orientations dictate the type of magnetic order under geometric frustration effects as predicted by Landau's mean-field theory. These theoretical findings provide insights into the exotic low-temperature phase observed in YbCuS$_2$, characterized by gapless excitations and seemingly nonmagnetic behavior accompanied by incommensurate correlations. △ Less

Submitted 22 June, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Sample complexity of matrix product states at finite temperature

Authors: Atsushi Iwaki, Chisa Hotta

Abstract: For quantum many-body systems in one dimension, computational complexity theory reveals that the evaluation of ground-state energy remains elusive on quantum computers, contrasting the existence of a classical algorithm for temperatures higher than the inverse logarithm of the system size. This highlights a qualitative difference between low- and high-temperature states in terms of computational c… ▽ More For quantum many-body systems in one dimension, computational complexity theory reveals that the evaluation of ground-state energy remains elusive on quantum computers, contrasting the existence of a classical algorithm for temperatures higher than the inverse logarithm of the system size. This highlights a qualitative difference between low- and high-temperature states in terms of computational complexity. Here, we describe finite-temperature states using the matrix product state formalism. Within the framework of random samplings, we derive an analytical formula for the required number of samples, which provides both quantitative and qualitative measures of computational complexity. At high and low temperatures, its scaling behavior with system size is linear and quadratic, respectively, demonstrating a distinct crossover between these numerically difficult regimes of quantitative difference. △ Less

Submitted 8 June, 2024; v1 submitted 15 March, 2024; originally announced March 2024.

Comments: 12 pages, 3 figures

Journal ref: Phys. Rev. B 109, 224410 (2024)

Exact matrix product states at the quantum Lifshitz tricritical point in a spin-1/2 zigzag-chain antiferromagnet with anisotropic $Γ$-term

Authors: Hidehiro Saito, Chisa Hotta

Abstract: Quantum anisotropic exchange interactions in magnets can induce competitions between phases in a different manner from those typically driven by geometrically frustrated interactions. We study a one-dimensional spin-1/2 zigzag chain with such an interaction, $Γ$-term, in conjunction with the Heisenberg interactions. We find a ground state phase diagram featuring a multicritical point where five ph… ▽ More Quantum anisotropic exchange interactions in magnets can induce competitions between phases in a different manner from those typically driven by geometrically frustrated interactions. We study a one-dimensional spin-1/2 zigzag chain with such an interaction, $Γ$-term, in conjunction with the Heisenberg interactions. We find a ground state phase diagram featuring a multicritical point where five phases converge: a uniform ferromagnet, two antiferromagnets, Tomonaga-Luttinger liquid and a dimer-singlet coexisting with nematic order. This multicritical point is simultaneously quantum tricritical and Lifshitz, and most remarkably, it hosts multi-degenerate ground state wave functions, whose exact form is obtained in the matrix product form and its degeneracy increases in squares of system size. △ Less

Submitted 12 January, 2024; originally announced January 2024.

Deriving quantum spin model for a zigzag-chain ytterbium magnet with anisotropic exchange interactions

Authors: Hidehiro Saito, Hiroki Nakai, Chisa Hotta

Abstract: We derive a quantum spin Hamiltonian of the spin-1/2 zigzag chain realized in a rare earth ytterbium-based magnetic insulator, YbCuS2. This material undergoes a transition at 0.95K to an incommensurate magnetic phase with small moments, which does not conform to the nonmagnetic singlet ground state of the spin-1/2 Heisenberg model. We take account of octahedral crystal field effect, atomic spin-or… ▽ More We derive a quantum spin Hamiltonian of the spin-1/2 zigzag chain realized in a rare earth ytterbium-based magnetic insulator, YbCuS2. This material undergoes a transition at 0.95K to an incommensurate magnetic phase with small moments, which does not conform to the nonmagnetic singlet ground state of the spin-1/2 Heisenberg model. We take account of octahedral crystal field effect, atomic spin-orbit coupling, and strong Coulomb interactions on Yb ions, and perform a four-order perturbation theory to evaluate the superexchange coupling constants. A small but finite anisotropic exchange coupling called Γ-term appears similarly to the case reported previously in other triangular-based magnets. By varying several material parameters, we figure out two important factors to enhance Γ-term. One is the splitting of excited f-states with two holes, which efficiently selects the perturbation terms associated with the lowest excited state having large total angular momentum, and accordingly with high spatial anisotropy. The other is the tilting of octahedra or distortion of S-Yb-S bond angle, which break the symmetry of the Slater-Koster overlap. These effects are systematically analyzed by the exchange anisotropy in units of pairs of octahedra within the local spin frame, which significantly reduces the complexity of directly referring to the Hamiltonian discussed in the global axis. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Thermal pure matrix product state in two dimensions: tracking thermal equilibrium from paramagnet down to the Kitaev honeycomb spin liquid state

Authors: Matthias Gohlke, Atsushi Iwaki, Chisa Hotta

Abstract: We present the first successful application of the matrix product state (MPS) representing a thermal quantum pure state (TPQ) in equilibrium in two spatial dimensions over almost the entire temperature range. We use the Kitaev honeycomb model as a prominent example hosting a quantum spin liquid (QSL) ground state to target the two specific-heat peaks previously solved nearly exactly using the free… ▽ More We present the first successful application of the matrix product state (MPS) representing a thermal quantum pure state (TPQ) in equilibrium in two spatial dimensions over almost the entire temperature range. We use the Kitaev honeycomb model as a prominent example hosting a quantum spin liquid (QSL) ground state to target the two specific-heat peaks previously solved nearly exactly using the free Majorana fermionic description. Starting from the high-temperature random state, our TPQ-MPS framework on a cylinder precisely reproduces these peaks, showing that the quantum many-body description based on spins can still capture the emergent itinerant Majorana fermions in a ${\mathbb Z}_2$ gauge field. The truncation process efficiently discards the high-energy states, eventually reaching the long-range entangled topological state approaching the exact ground state for a given finite size cluster. An advantage of TPQ-MPS over exact diagonalization or purification-based methods is its lowered numerical cost coming from a reduced effective Hilbert space even at finite temperature. △ Less

Submitted 17 October, 2023; v1 submitted 3 August, 2023; originally announced August 2023.

Comments: 15 pages, 4 figures, submission to SciPost

Journal ref: SciPost Phys. 15, 206 (2023)

Sine-square deformed mean-field theory

Authors: Masataka Kawano, Chisa Hotta

Abstract: We develop a theory that accurately evaluates quantum phases with any large-scale emergent structures including incommensurate density waves or topological textures without {\it a priori} knowing their periodicity. We spatially deform a real-space mean-field Hamiltonian on a finite-size cluster using a sine-squared envelope function with zero energy at system edges. The wave functions become insen… ▽ More We develop a theory that accurately evaluates quantum phases with any large-scale emergent structures including incommensurate density waves or topological textures without {\it a priori} knowing their periodicity. We spatially deform a real-space mean-field Hamiltonian on a finite-size cluster using a sine-squared envelope function with zero energy at system edges. The wave functions become insensitive to the misfit of the lattice and ordering periods. We successfully extract the ordering wave vectors by our deformed Fourier transformation, updating the previous results for hole-doped and spin-orbit coupled Mott insulators. The method further enables the evaluation of a charge gap beyond the mean-field level. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 7 pages, 4 figures

Journal ref: Phys. Rev. Research 4, L012033 (2022)

Spin-orbit coupled Hubbard skyrmions

Authors: Ryo Makuta, Chisa Hotta

Abstract: When the material phases exhibit topological quantum numbers, they host defects protected by the nontrivial topology. Magnetic skyrmions are such ``quantized" objects and although many of them are metals they had been most likely treated in theories as purely classical spin states. Here, we show that the electrons described by the Hubbard model with strong spin-orbit couplings can exhibit various… ▽ More When the material phases exhibit topological quantum numbers, they host defects protected by the nontrivial topology. Magnetic skyrmions are such ``quantized" objects and although many of them are metals they had been most likely treated in theories as purely classical spin states. Here, we show that the electrons described by the Hubbard model with strong spin-orbit couplings can exhibit various nano- or flake-size skyrmions in their ground state. Importantly, the conducting electrons forming Fermi pockets themselves carry textured magnetic moments in both real and momentum space and on top of that, possess Chern numbers in their energy bands. This quantum and conducting skyrmion is related to small skyrmions observed in atomic-layered compounds. We clarify how the effective magnetic interactions and magnetic anisotropies are tied to the spin-orbit coupling, and how they influence the stability of skyrmions beyond the phenomenology. △ Less

Submitted 9 May, 2024; v1 submitted 17 April, 2023; originally announced April 2023.

Comments: 15 pages, 11 figures, 1 table

Journal ref: Phys. Rev. Research 6, 023133 (2024)

Emergent SU(3) magnons and thermal Hall effect in the antiferromagnetic skyrmion lattice

Authors: Hikaru Takeda, Masataka Kawano, Kyo Tamura, Masatoshi Akazawa, Jian Yan, Takeshi Waki, Hiroyuki Nakamura, Kazuki Sato, Yasuo Narumi, Masayuki Hagiwara, Minoru Yamashita, Chisa Hotta

Abstract: Complexity of quantum phases of matter is often understood by the underlying gauge structures, as was recognized by the $\mathbb{Z}_2$ and U(1) gauge theory description of spin liquid in frustrated magnets. Anomalous Hall effect of conducting electrons can intrisically arise from U(1) gauges expressing the spatial modulation of ferromagnetic moments or from SU(2) gauges representing the spin-orbit… ▽ More Complexity of quantum phases of matter is often understood by the underlying gauge structures, as was recognized by the $\mathbb{Z}_2$ and U(1) gauge theory description of spin liquid in frustrated magnets. Anomalous Hall effect of conducting electrons can intrisically arise from U(1) gauges expressing the spatial modulation of ferromagnetic moments or from SU(2) gauges representing the spin-orbit coupling effect. Similarly, in insulating ferro and antiferromagnets, the magnon excitations can contribute to anomalous transports by feeling the U(1) and SU(2) gauges arising from the features of ordered moments or interactions. In this work, we report the emergent higher rank SU(3) gauge structure in the magnon transport based on the thermal conductivity measurements of MnSc$_2$S$_4$ in an applied field up to 14\,T. The thermal Hall coefficient takes a substantial value when the material enters a three-sublattice antiferromagnetic skyrmion phase, which is confirmed by the large-scale spin wave theory. The excited magnons are dressed with SU(3) gauge field, which is a mixture of three species of U(1) gauge fields originating from the slowly varying magnetic moments on these sublattices. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 16 pages, 12 figures

Journal ref: Nat. Commun. 15, 566 (2024)

Quantum critical dynamics in two-dimensional transverse Ising model

Authors: Chisa Hotta, Tempei Yoshida, Kenji Harada

Abstract: In the vicinity of the quantum critical point(QCP), thermodynamic properties diverge toward zero temperature governed by universal exponents. Although this fact is well known, how it is reflected in quantum dynamics has not been addressed. As an ideal experimental platform to test the issue, we consider an organic Mott insulator whose dielectric degrees of freedom, a quantum electric dipole, is de… ▽ More In the vicinity of the quantum critical point(QCP), thermodynamic properties diverge toward zero temperature governed by universal exponents. Although this fact is well known, how it is reflected in quantum dynamics has not been addressed. As an ideal experimental platform to test the issue, we consider an organic Mott insulator whose dielectric degrees of freedom, a quantum electric dipole, is described by the transverse Ising model on a triangular lattice that has a QCP. We track the Glauber-type dynamics of the model by constructing a kinetic protocol based on the quantum Monte Carlo method. The dynamical susceptibility takes the form of the Debye function and shows a significant peak-narrowing in approaching a QCP due to the divergence of the relaxation timescale. It resembles the anomaly of dielectric constants observed in the organic materials \k{appa}-ET 2 X dimer Mott insulating phase, indicating that the material is very near the ferroelectric QCP. △ Less

Submitted 30 November, 2022; v1 submitted 23 September, 2022; originally announced September 2022.

Comments: 12 pages, 8 figures

Phase diagram of the square lattice Hubbard model with Rashba-type antisymmetric spin-orbit coupling

Authors: Masataka Kawano, Chisa Hotta

Abstract: We clarify the ground state phase diagram of the half-filled square-lattice Hubbard model with Rashba spin-orbit coupling (SOC) characterized by the spin-split energy bands due to broken inversion symmetry. Although the Rashba metals and insulating magnets have been studied well, the intermediate interaction strength of the system remained elusive due to the lack of appropriate theoretical tools t… ▽ More We clarify the ground state phase diagram of the half-filled square-lattice Hubbard model with Rashba spin-orbit coupling (SOC) characterized by the spin-split energy bands due to broken inversion symmetry. Although the Rashba metals and insulating magnets have been studied well, the intermediate interaction strength of the system remained elusive due to the lack of appropriate theoretical tools to unbiasedly describe the large-scale magnetic structures. We complementarily apply four different methods; sine square deformed mean-field theory, random phase approximation, Luttinger-Tisza method, and density matrix embedding theory, and succeed in capturing the incommensurate spin density-wave (SDW) phases with very long spatial periods which were previously overlooked. The transition to the SDW phases from the metallic phase is driven by an unprecedented instability that nests the two parts of the Fermi surface carrying opposite spins. For large SOC, the spiral, stipe, and vortex phases are obtained, when the four Dirac points exist near the Fermi level and their whole linear dispersions nest by a wavelength $π$, opening a band gap. These two types of transition provide Fermiology that distinguishes the antisymmetric SOC systems, generating a variety of magnetic phases that start from the relatively weak correlation regime and continue to the strongly interacting limit. △ Less

Submitted 18 January, 2023; v1 submitted 21 August, 2022; originally announced August 2022.

Comments: 20 pages, 10 figures

Journal ref: Phys. Rev. B 107, 045123 (2023)

Emergent chiral symmetry in non-bipartite kagome and pyrochlore lattices with spin-orbit coupling

Authors: Hiroki Nakai, Masataka Kawano, Chisa Hotta

Abstract: Chiral symmetry in energy bands appears as perfectly symmetric anti-bonding and bonding pairs of energy levels. It has only been observed in a few classes of models with a bipartite lattice structure or Bogoliubov-de-Gennes systems having the pairwise basis. We show that the non-bipartite kagome and pyrochlore lattices can host chiral symmetric bands when the strong spin-orbit coupling is introduc… ▽ More Chiral symmetry in energy bands appears as perfectly symmetric anti-bonding and bonding pairs of energy levels. It has only been observed in a few classes of models with a bipartite lattice structure or Bogoliubov-de-Gennes systems having the pairwise basis. We show that the non-bipartite kagome and pyrochlore lattices can host chiral symmetric bands when the strong spin-orbit coupling is introduced. There, the electrons hop to their neighbors by always converting the spin orientation up-side-down, which allows the up and down spin bases to form fictitious bipartite connections. The gauge invariant Wilson loop operator defined on a triangular unit serves as a marker to detect the presence of chiral symmetry, and using this property, the chiral operator is constructed. This allows us to access their topological symmetry classes that can easily change with small perturbations. △ Less

Submitted 15 August, 2022; originally announced August 2022.

Comments: 6 pages, 3 figures, 1 table, supplementary material

Journal ref: Phys. Rev. B 108, L081106 (2023)

Quantum mechanically driven structural-spin glass in two dimensions at finite temperature

Authors: Chisa Hotta, Kazumasa Ueda, Masatoshi Imada

Abstract: In magnetic materials, spins sometimes freeze into spatially disordered glassy states. Glass forming liquids or structural glasses are found very often in three dimensions. However, in two dimensions(2D) it is believed that both spin glass and structural glass can never exist at a finite temperature because they are destroyed by thermal fluctuations. Using a large-scale quantum Monte Carlo simulat… ▽ More In magnetic materials, spins sometimes freeze into spatially disordered glassy states. Glass forming liquids or structural glasses are found very often in three dimensions. However, in two dimensions(2D) it is believed that both spin glass and structural glass can never exist at a finite temperature because they are destroyed by thermal fluctuations. Using a large-scale quantum Monte Carlo simulation, we discover a quantum-mechanically driven 2D glass phase at finite temperatures. Our platform is an Ising spin model with a quantum transverse field on a frustrated triangular lattice. How the present glass phase is formed is understood by the following three steps. First, by the interplay of geometrical frustration and quantum fluctuation, part of the spins spontaneously form an antiferromagnetic honeycomb spin-superstructure. Then, small randomness in the bond interaction works as a relevant perturbation to this superstructure and breaks it up into {\it domains}, making it a structural glass. The glassiness of the superstructure, in turn, generates an emergent random magnetic field acting on the remaining fluctuating spins and freezes them. The shape of domains thus formed depends sensitively on the quenching process, which is one of the characteristic features of glass, originating from a multi-valley free-energy landscape. The present system consists only of {\it a single} bistable Ising degree of freedom, which naturally does not become a structural glass alone nor a spin glass alone. Nevertheless, a glass having both types of nature emerges in the form of coexisting two-component glasses, algebraic structural-glass and long-range ordered spin-glass. This new concept of glass-forming mechanism opens a way to realize functional glasses even in low dimensional systems. △ Less

Submitted 15 July, 2022; originally announced July 2022.

Comments: 21pages, 14figures

Purity of thermal mixed quantum states

Authors: Atsushi Iwaki, Chisa Hotta

Abstract: We develop a formula to evaluate the purity of a series of thermal equilibrium states that can be calculated in numerical experiments without knowing the exact form of the quantum state \textit{a priori}. Canonical typicality guarantees that there are numerous microscopically different expressions of such states, which we call thermal mixed quantum (TMQ) states. Suppose that we construct a TMQ sta… ▽ More We develop a formula to evaluate the purity of a series of thermal equilibrium states that can be calculated in numerical experiments without knowing the exact form of the quantum state \textit{a priori}. Canonical typicality guarantees that there are numerous microscopically different expressions of such states, which we call thermal mixed quantum (TMQ) states. Suppose that we construct a TMQ state by a mixture of $N_\mathrm{samp}$ independent pure states. The weight of each pure state is given by its norm, and the partition function is given by the average of the norms. To qualify how efficiently the mixture is done, we introduce a quantum statistical quantity called "normalized fluctuation of partition function (NFPF)". For smaller NFPF, the TMQ state is closer to the equally weighted mixture of pure states, which means higher efficiency, requiring a smaller $N_\mathrm{samp}$. The largest NFPF is realized in the Gibbs state with purity-0 and exponentially large $N_\mathrm{samp}$, while the smallest NFPF is given for thermal pure quantum state with purity-1 and $N_\mathrm{samp}=1$. The purity is formulated using solely the NFPF and roughly gives $N_\mathrm{samp}^{-1}$. Our analytical results are numerically tested and confirmed by the two random sampling methods built on matrix-product-state-based wave functions. △ Less

Submitted 11 September, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: 20 pages, 9 figures

Journal ref: Phys. Rev. B 106, 094409 (2022)

Supercooled Jahn-Teller Ice

Authors: Kota Mitsumoto, Chisa Hotta, Hajime Yoshino

Abstract: When the spins on the frustrated pyrochlore lattice obey the celebrated 2-$in$-2-$out$ ice rule, they stay in a correlated disordered phase and break the third law of thermodynamics. Similarly, if the atomic ions on the pyrochlore lattice move in and outward of the tetrahedra, they may obey a constraint resembling the ice rule. We discover that a model for pyrochlore molybdates $A_2$Mo$_2$O$_7$ (… ▽ More When the spins on the frustrated pyrochlore lattice obey the celebrated 2-$in$-2-$out$ ice rule, they stay in a correlated disordered phase and break the third law of thermodynamics. Similarly, if the atomic ions on the pyrochlore lattice move in and outward of the tetrahedra, they may obey a constraint resembling the ice rule. We discover that a model for pyrochlore molybdates $A_2$Mo$_2$O$_7$ ($A=$Y, Dy, Tb) exhibits a "supercooled ice" state of the displacement degrees of freedom of Mo$^{4+}$ ions, when we take account of the Jahn-Teller (JT) effect. The JT effect occurs when the lattice distortions reduce the symmetry of the local crystal field, resulting in the orbital-energy-splitting that causes the local energy gain. Unlike the standard JT effect that leads to periodic long range ordering, the displacements of Mo$^{4+}$ ions are disordered following the ice-like rule. We microscopically derive a model that describes this situation by having the 2nd and 3rd neighbor interactions between in-out lattice displacements comparably as strong as the nearest neighbor interactions of standard ice. There, the well-known nearly flat energy landscape of the ice state is altered to a metastable highly quasi-degenerate ice-like liquid state coexisting with a crystalline-like ground state. Our Monte Carlo simulations show that this liquid remains remarkably stable down to low temperatures by avoiding the putative first order transition. The relaxation in the supercooled JT ice state exhibits glassy dynamics with a plateau structure. They fit the feature of a "good glassformer" very often found in molecular liquids, but that has never been observed in material solids. The high glass-forming ability of the interacting lattice degrees of freedom will play a key role in the spin-glass transition of the material. △ Less

Submitted 30 August, 2022; v1 submitted 11 February, 2022; originally announced February 2022.

Comments: 18 pages, 14 figures

Journal ref: Phys. Rev. Research 4, 033157 (2022)

Dimensional reduction in quantum spin-1/2 system on a 1/7-depleted triangular lattice

Authors: Ryo Makuta, Chisa Hotta

Abstract: We study the magnetism of a quantum spin-1/2 antiferromagnet on a maple-leaf lattice which is obtained by regularly depleting 1/7 of the sites of a triangular lattice. Although the interactions are set to be spatially uniform, the ground state shows a stripe Neel order and the temperature dependence of magnetic susceptibility follows that of the one-dimensional XXZ model with a finite spin gap. We… ▽ More We study the magnetism of a quantum spin-1/2 antiferromagnet on a maple-leaf lattice which is obtained by regularly depleting 1/7 of the sites of a triangular lattice. Although the interactions are set to be spatially uniform, the ground state shows a stripe Neel order and the temperature dependence of magnetic susceptibility follows that of the one-dimensional XXZ model with a finite spin gap. We examine the nature of frustration by mapping the low energy degenerate manifold of states to the fully packed loop-string model on a dual cluster-depleted honeycomb lattice, finding that the order-by-disorder due to quantum fluctuation characteristic of highly frustrated magnets is responsible for the emergent stripes. The excited magnons split into two spinons and propagate in the one-dimensional direction along the stripe which is reminiscent of the XXZ or Ising model in one dimension. Unlike most of the previously studied dimensional reduction effects, our case is purely spontaneous as the interactions of the Hamiltonian retains a spatially isotropic two-dimensional structure. △ Less

Submitted 13 December, 2021; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: 14 pages, 14 figures

Journal ref: Phys. Rev. B 104, 224415 (2021)

Sine-square deformation applied to classical Ising models

Authors: Chisa Hotta, Takashi Nakamaniwa, Tota Nakamura

Abstract: Sine-square deformation (SSD) is a treatment proposed in quantum systems, which spatially modifies a Hamiltonian, gradually decreasing the local energy scale from the center of the system toward the edges by a sine-squared envelope function. It is known to serve as a good boundary condition as well as to provide physical quantities reproducing those of the infinite-size systems. We apply the SSD t… ▽ More Sine-square deformation (SSD) is a treatment proposed in quantum systems, which spatially modifies a Hamiltonian, gradually decreasing the local energy scale from the center of the system toward the edges by a sine-squared envelope function. It is known to serve as a good boundary condition as well as to provide physical quantities reproducing those of the infinite-size systems. We apply the SSD to one- and two-dimensional classical Ising models. Based on the analytical calculations and Monte Carlo simulations, we find that the classical SSD system is regarded as an extended canonical ensemble of a local subsystem each characterized by its own effective temperature. This effective temperature is defined by normalizing the system temperature by the deformed local energy scale. A single calculation for a fixed system temperature provides a set of physical quantities of various temperatures that quantitatively reproduce well those of the uniform system. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 15pages 8 figures

Perfect flat band with chirality and charge ordering out of strong spin-orbit interaction

Authors: Hiroki Nakai, Chisa Hotta

Abstract: Spin-orbit interaction established itself as a major role player for emergent phenomena in modern condensed matter including a topological insulator, spin liquid and spin-dependent transports. However, its function is rather limited to adding topological nature to each phases of matter. We prove by our spinor line graph theory that a very strong spin-orbit interaction realized in 5d pyrochlore ele… ▽ More Spin-orbit interaction established itself as a major role player for emergent phenomena in modern condensed matter including a topological insulator, spin liquid and spin-dependent transports. However, its function is rather limited to adding topological nature to each phases of matter. We prove by our spinor line graph theory that a very strong spin-orbit interaction realized in 5d pyrochlore electronic systems generates multiply degenerate perfect flat bands. Unlike any of the previous flat bands, the electrons living there localize in real space by destructively interfering with each other in a spin selective manner ruled by the SU(2) gauge field. These electrons avoid the Coulomb interaction by self-organizing their localized wave functions. This gives rise to the trimerized charge ordering hand in hand with a stiff spin chirality, which may explain the recently found exotic low-temperature insulating phase of CsW2O6. △ Less

Submitted 17 September, 2021; v1 submitted 25 March, 2021; originally announced March 2021.

Comments: 11pages, 7 figures

Comparative study of the density matrix embedding theory for the Hubbard models

Authors: Masataka Kawano, Chisa Hotta

Abstract: We examine the performance of the density matrix embedding theory (DMET) recently proposed in [G. Knizia and G. K.-L. Chan, Phys. Rev. Lett. 109, 186404 (2012)]. The core of this method is to find a proper one-body potential that generates a good trial wave function for projecting a large scale original Hamiltonian to a local subsystem with a small number of basis. The resultant ground state of th… ▽ More We examine the performance of the density matrix embedding theory (DMET) recently proposed in [G. Knizia and G. K.-L. Chan, Phys. Rev. Lett. 109, 186404 (2012)]. The core of this method is to find a proper one-body potential that generates a good trial wave function for projecting a large scale original Hamiltonian to a local subsystem with a small number of basis. The resultant ground state of the projected Hamiltonian can locally approximate the true ground state. However, the lack of the variational principle makes it difficult to judge the quality of the choice of the potential. Here we focus on the entanglement spectrum (ES) as a judging criterion; accurate evaluation of the ES guarantees that the corresponding reduced density matrix well reproduces all physical quantities on the local subsystem. We apply the DMET to the Hubbard model on the one-dimensional chain, zigzag chain, and triangular lattice and test several variants of potentials and cost functions. It turns out that ES serves as a more sensitive quantity than the energy and double occupancy to probe the quality of the DMET outcomes. A symmetric potential reproduces the ES of the phase that continues from a noninteracting limit. The Mott transition as well as symmetry-breaking transitions can be detected by the singularities in the ES. However, the details of the ES in the strongly interacting parameter region depends much on these variants, meaning that the present DMET algorithm allowing for numerous variant is insufficient to fully characterize the particular phases that require characterization by the ES. △ Less

Submitted 3 December, 2020; v1 submitted 23 September, 2020; originally announced September 2020.

Comments: 12 pages, 6 figures

Journal ref: Phys. Rev. B 102, 235111 (2020)

Entanglement spectrum as a marker for phase transitions in the density embedding theory for interacting spinless fermionic models

Authors: Xavier Plat, Chisa Hotta

Abstract: Entanglement related properties work as nice fingerprint of the quantum many-body wave function. However, those of fermionic models are hard to evaluate in standard numerical methods because they suffer from finite size effects. We show that a so-called density embedding theory (DET) can evaluate them without size scaling analysis in comparably high quality with those obtained by the large-size de… ▽ More Entanglement related properties work as nice fingerprint of the quantum many-body wave function. However, those of fermionic models are hard to evaluate in standard numerical methods because they suffer from finite size effects. We show that a so-called density embedding theory (DET) can evaluate them without size scaling analysis in comparably high quality with those obtained by the large-size density matrix renormalization group analysis. This method projects the large scale original many-body Hamiltonian to the small number of basis sets defined on a local cluster, and optimizes the choice of these bases by tuning the local density matrix. The DET entanglement spectrum of one-dimensional interacting fermions perfectly reproduces the exact ones and works as a marker of the phase transition point. It is further shown that the phase transitions in two-dimension could be determined by the entanglement entropy and the fidelity that reflects the change of the structure of the wave function. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 3figures 5pages

Journal ref: Phys. Rev. B 102, 140410 (2020)

Ferromagnetism out of charge fluctuation of strongly correlated electrons in $κ$-(BEDT-TTF)$_2$Hg(SCN)$_2$Br

Authors: Minoru Yamashita, Shiori Sugiura, Akira Ueda, Shun Dekura, Taichi Terashima, Shinya Uji, Yoshiya Sunairi, Hatsumi Mori, Elena I. Zhilyaeva, Svetlana A. Torunova, Rimma N. Lyubovskaya, Natalia Drichko, Chisa Hotta

Abstract: We perform magnetic susceptibility and magnetic torque measurements on the organic $κ$-(BEDT-TTF)$_2$Hg(SCN)$_2$Br, which is recently suggested to host an exotic quantum dipole-liquid in its low-temperature insulating phase. Below the metal-insulator transition temperature, the magnetic susceptibility follows a Curie-Weiss law with a positive Curie-Weiss temperature, and a particular… ▽ More We perform magnetic susceptibility and magnetic torque measurements on the organic $κ$-(BEDT-TTF)$_2$Hg(SCN)$_2$Br, which is recently suggested to host an exotic quantum dipole-liquid in its low-temperature insulating phase. Below the metal-insulator transition temperature, the magnetic susceptibility follows a Curie-Weiss law with a positive Curie-Weiss temperature, and a particular $M\propto \sqrt{H}$ curve is observed. The emergent ferromagnetically interacting spins amount to about 1/6 of the full spin moment of localized charges. Taking account of the possible inhomogeneous quasi-charge-order that forms a dipole-liquid, we construct a model of antiferromagnetically interacting spin chains in two adjacent charge-ordered domains, which are coupled via fluctuating charges on a Mott-dimer at the boundary. We find that the charge fluctuations can draw a weak ferromagnetic moment out of the spin singlet domains. △ Less

Submitted 24 August, 2021; v1 submitted 9 June, 2020; originally announced June 2020.

Comments: 9 pages, 5 figures, Supplementary Material attached

Journal ref: npj Quantum Materials volume 6, 87 (2021)

Thermal Pure Quantum Matrix Product States Recovering a Volume Law Entanglement

Authors: Atsushi Iwaki, Akira Shimizu, Chisa Hotta

Abstract: We propose a way to construct a thermal pure quantum matrix product state (TPQ-MPS) that can simulate finite temperature quantum many-body systems with a minimal numerical cost comparable to the matrix product algorithm for the ground state. The MPS was originally designed for the wave function with area-law entanglement. However, by attaching the auxiliary sites to the edges of the random matrix… ▽ More We propose a way to construct a thermal pure quantum matrix product state (TPQ-MPS) that can simulate finite temperature quantum many-body systems with a minimal numerical cost comparable to the matrix product algorithm for the ground state. The MPS was originally designed for the wave function with area-law entanglement. However, by attaching the auxiliary sites to the edges of the random matrix product state, we find that the degree of entanglement is automatically tuned so as to recover the volume law of the entanglement entropy that characterizes the TPQ state. The finite temperature physical quantities of the transverse Ising and the spin-1/2 Heisenberg chains evaluated by a TPQ-MPS show excellent agreement even for bond dimension $\sim 10$-$20$ with those of the exact results. △ Less

Submitted 28 May, 2021; v1 submitted 14 May, 2020; originally announced May 2020.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. Research 3, 022015 (2021)

Finite temperature thermodynamic properties of the spin-1 nematics in an applied magnetic field

Authors: Katsuhiro Tanaka, Chisa Hotta

Abstract: We study numerically the thermodynamic properties of the spin nematic phases in a magnetic field in the spin-1 bilinear-biquadratic model. When the field is applied, the phase transition temperature once goes up and then decreases rapidly toward zero, which is detected by the peak-shift in the specific heat. The underlying mechanism of the reentrant behavior is the entropic effect. In a weak field… ▽ More We study numerically the thermodynamic properties of the spin nematic phases in a magnetic field in the spin-1 bilinear-biquadratic model. When the field is applied, the phase transition temperature once goes up and then decreases rapidly toward zero, which is detected by the peak-shift in the specific heat. The underlying mechanism of the reentrant behavior is the entropic effect. In a weak field the high temperature paramagnetic phase rapidly loses its entropy while the ferroquadrupolar nematic phase remains robust by modifying the shape of the ferroquadrupolar moment. This feature serves as a fingerprint of generic ferroquadrupolar phases, while it is not observed for the case of antiferroquadrupoles. △ Less

Submitted 16 September, 2020; v1 submitted 22 April, 2020; originally announced April 2020.

Comments: 6 pages, 4 figures + supplemental material

Journal ref: Phys. Rev. B 102, 140401 (2020)

Discovering momentum-dependent magnon spin texture in insulating antiferromagnets: Role of the Kitaev interaction

Authors: Masataka Kawano, Chisa Hotta

Abstract: We theoretically show that the Kitaev interaction generates a novel class of spin texture in the excitation spectrum of the antiferromagnetic insulator found in the Kitaev-Heisenberg-$Γ$ model. In conducting electronic systems, there is a series of vortex type of spin texture along the Fermi surface induced by Rashba and Dresselhaus spin-orbit coupling. Such spin textures are rarely found in magne… ▽ More We theoretically show that the Kitaev interaction generates a novel class of spin texture in the excitation spectrum of the antiferromagnetic insulator found in the Kitaev-Heisenberg-$Γ$ model. In conducting electronic systems, there is a series of vortex type of spin texture along the Fermi surface induced by Rashba and Dresselhaus spin-orbit coupling. Such spin textures are rarely found in magnetic insulators, since there had been no systematic ways to control the kinetics of its quasi-particle called magnon using a magnetic field or spacially asymmetric exchange couplings. Here, we propose a general framework to explore such spin textures in arbitrary insulating antiferromagnets. We introduce an analytical method to transform any complicated Hamiltonian to the simple representation based on pseudo-spin degrees of freedom. The direction of the pseudo-spin on a Bloch sphere describes the degree of contributions from the two magnetic sublattices to the spin moment carried by the magnon. The momentum dependent fictitious "Zeeman field" determines the direction of the pseudo-spin and thus becomes the control parameter of the spin texture, which is explicitly described by the original model parameters. The framework enabled us to clarify the uncovered aspect of the Kitaev interaction, and further provides a tool to easily design or explore materials with intriguing magnetic properties. Since these spin textures can be a source of a pure spin current, the Kitaev materials $A_{2}$PrO$_{3}$ ($A$ =Li, Na) shall become a potential platform of power-saving spintronics devices. △ Less

Submitted 6 November, 2019; v1 submitted 11 September, 2019; originally announced September 2019.

Comments: 19 pages, 6 figures

Journal ref: Phys. Rev. B 100, 174402 (2019)

Abundant quadrupolar or nematic phases driven by the Heisenberg interactions in a spin-1 dimer system forming a bilayer

Authors: Katsuhiro Tanaka, Chisa Hotta

Abstract: We explore several classes of quadrupolar ordering in a system of antiferromagnetically coupled quantum spin-1 dimers, which are stacked in the triangular lattice geometry forming a bilayer. Low-energy properties of this model is described by an $\mathcal{S}=1$ hard-core bosonic degrees of freedom defined on each dimer-bond, where the singlet and triplet states of the dimerized spins are interpret… ▽ More We explore several classes of quadrupolar ordering in a system of antiferromagnetically coupled quantum spin-1 dimers, which are stacked in the triangular lattice geometry forming a bilayer. Low-energy properties of this model is described by an $\mathcal{S}=1$ hard-core bosonic degrees of freedom defined on each dimer-bond, where the singlet and triplet states of the dimerized spins are interpreted as the vacuum and the occupancy of boson, respectively. The number of bosons per dimer and the magnetic and density fluctuations of bosons are controlled by the inter-dimer Heisenberg interactions. In a solid phase where each dimer hosts one boson and the inter-dimer interaction is weak, a conventional spin nematic phase is realized by the pair-fluctuation of bosons. Larger inter-dimer interaction favors Bose Einstein condensates (BEC) carrying quadrupolar moments. Among them, we find one exotic phase where the quadrupoles develop a spatially modulated structure on the top of a uniform BEC, interpreted in the original dimerized spin-1 model as coexistent $p$-type nematic and 120$^\circ$-magnetic correlations. This may explain an intriguing nonmagnetic phase found in Ba$_{3}$ZnRu$_{2}$O$_{9}$. △ Less

Submitted 5 January, 2020; v1 submitted 10 September, 2019; originally announced September 2019.

Comments: 16 pages, 12 figures, manuscript is updated, 3 figures are added

Spin-orbital glass transition in a model of a frustrated pyrochlore magnet without quenched disorder

Authors: Kota Mitsumoto, Chisa Hotta, Hajime Yoshino

Abstract: We show theoretically that spin and orbital degrees of freedom in the pyrochlore oxide Y2Mo2O7, which is free of quenched disorder, can exhibit a simultaneous glass transition, working as dynamical randomness to each other. The interplay of spins and orbitals is mediated by the Jahn-Teller lattice distortion that selects the choice of orbitals, which then generates variant spin exchange interactio… ▽ More We show theoretically that spin and orbital degrees of freedom in the pyrochlore oxide Y2Mo2O7, which is free of quenched disorder, can exhibit a simultaneous glass transition, working as dynamical randomness to each other. The interplay of spins and orbitals is mediated by the Jahn-Teller lattice distortion that selects the choice of orbitals, which then generates variant spin exchange interactions ranging from ferromagnetic to antiferromagnetic ones. Our Monte Carlo simulations detect the power-law divergence of the relaxation times and the negative divergence of both the magnetic and dielectric non-linear susceptibilities, resolving the long-standing puzzle on the origin of the disorder-free spin glass. △ Less

Submitted 27 February, 2020; v1 submitted 14 August, 2019; originally announced August 2019.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. Lett. 124, 087201 (2020)

Magnetic susceptibility of quantum spin systems calculated by sine square deformation: one-dimensional, square lattice, and kagome lattice Heisenberg antiferromagnets

Authors: Chisa Hotta, Kenichi Asano

Abstract: We develop a simple and unbiased numerical method to obtain the uniform susceptibility of quantum many body systems. When a Hamiltonian is spatially deformed by multiplying it with a sine square function that smoothly decreases from the system center toward the edges, the size-scaling law of the excitation energy is drastically transformed to a rapidly converging one. Then, the local magnetization… ▽ More We develop a simple and unbiased numerical method to obtain the uniform susceptibility of quantum many body systems. When a Hamiltonian is spatially deformed by multiplying it with a sine square function that smoothly decreases from the system center toward the edges, the size-scaling law of the excitation energy is drastically transformed to a rapidly converging one. Then, the local magnetization at the system center becomes nearly size independent; the one obtained for the deformed Hamiltonian of a system length as small as L=10 provides the value obtained for the original uniform Hamiltonian of L=100. This allows us to evaluate a bulk magnetic susceptibility by using the magnetization at the center by existing numerical solvers without any approximation, parameter tuning, or the size-scaling analysis. We demonstrate that the susceptibilities of the spin-1/2 antiferromagnetic Heisenberg chain and square lattice obtained by our scheme at L=10 agree within 10 to (-3) with exact analytical and numerical solutions for L=infinite down to temperature of 0.1 times the coupling constant. We apply this method to the spin-1/2 kagome lattice Heisenberg antiferromagnet which is of prime interest in the search of spin liquids. △ Less

Submitted 13 September, 2018; originally announced September 2018.

Comments: 5pages 3figures

Journal ref: Phys. Rev. B 98, 140405 (2018)

From Linear to Nonlinear Responses of Thermal Pure Quantum States

Authors: Hiroyuki Endo, Chisa Hotta, Akira Shimizu

Abstract: We propose a self-validating scheme to calculate the unbiased responses of quantum many-body systems to external fields of arbibraty strength at any temperature. By switching on a specified field to a thermal pure quantum state of an isolated system, and tracking its time evolution, one can observe an intrinsic thermalization process driven solely by many-body effects. The transient behavior befor… ▽ More We propose a self-validating scheme to calculate the unbiased responses of quantum many-body systems to external fields of arbibraty strength at any temperature. By switching on a specified field to a thermal pure quantum state of an isolated system, and tracking its time evolution, one can observe an intrinsic thermalization process driven solely by many-body effects. The transient behavior before thermalization contains rich information on excited states, giving the linear and nonlinear response functions at all frequencies. We uncover the necessary conditions to clarify the applicability of this formalism, supported by a proper definition of the nonlinear response function. The accuracy of the protocol is guaranteed by a rigorous upper bound of error exponentially decreasing with system size, and is well implemented in the simple ferromagnetic Heisenberg chain, whose response at high fields exhibits a nonlinear band deformation. We further extract the characteristic features of excitation of the spin-1/2 kagome antiferromagnet; the wavenumber-insensitive linear responses from the possible spin liquid ground state, and the significantly broad nonlinear peaks which should be generated from numerous collisions of quasi-particles, that are beyond the perturbative description. △ Less

Submitted 25 November, 2018; v1 submitted 6 June, 2018; originally announced June 2018.

Comments: published version; 6 pages with 3 figures, and Supplemental Material

Journal ref: Phys. Rev. Lett. 121, 220601 (2018)

Origin of Biquadratic Exchange Interactions in a Mott Insulator as a Driving Force of Spin Nematic Order

Authors: Katsuhiro Tanaka, Yuto Yokoyama, Chisa Hotta

Abstract: We consider a series of Mott insulators in unit of two orbitals each hosting spin-1/2 electron, and by pairing two spin-1/2 into spin-1 triplet, derive the effective exchange interaction between the adjacent units via fourth order perturbation theory. It turns out that the biquadratic exchange interaction between spin-1, which is one of the origins of the nematic order, arises only in processes wh… ▽ More We consider a series of Mott insulators in unit of two orbitals each hosting spin-1/2 electron, and by pairing two spin-1/2 into spin-1 triplet, derive the effective exchange interaction between the adjacent units via fourth order perturbation theory. It turns out that the biquadratic exchange interaction between spin-1, which is one of the origins of the nematic order, arises only in processes where the four different electrons exchange cyclically along the twisted loop, which we call "twisted ring exchange" processes. We show that the term becomes the same order with the Heisenberg exchange interactions when the on-orbital Coulomb interaction is not too large. Whereas, the inter-orbital Coulomb interactions give rise to additional processes that cancel the twisted ring exchange, and strongly suppresses the biquadratic term. The Mott insulator with two electrons on degenerate two orbitals is thus not an ideal platform to study such nematic orders. △ Less

Submitted 21 May, 2018; originally announced May 2018.

Journal ref: J. Phys. Soc. Jpn. 87, 023702 (2018)

Spin nematics next to spin singlets

Authors: Yuto Yokoyama, Chisa Hotta

Abstract: We provide a route to generate nematic order in a spin-1/2 system. Unlike the well-known magnon-binding mechanism, our spin nematics requires neigher the frustration effect nor a spin polarization in a high field or in the vicinity of a ferromagnet, but instead appears next to the spin singlet phase. We start from a state consisting of a quantum spin-1/2 singlet dimer placed on each site of a tria… ▽ More We provide a route to generate nematic order in a spin-1/2 system. Unlike the well-known magnon-binding mechanism, our spin nematics requires neigher the frustration effect nor a spin polarization in a high field or in the vicinity of a ferromagnet, but instead appears next to the spin singlet phase. We start from a state consisting of a quantum spin-1/2 singlet dimer placed on each site of a triangular lattice, and show that inter-dimer ring exchange interactions efficiently dope the SU(2) triplets that itinerate and interact, easily driving a stable singlet state to either Bose Einstein condensates or a triplet crystal, some hosting a spin nematic order. A variety of roles the ring exchange serves include the generation of a bilinear-biquadratic interaction between nearby triplets, which is responsible for the emergent nematic order separated from the singlet phase by a first order transition. △ Less

Submitted 21 May, 2018; originally announced May 2018.

Comments: 6 pages , 4 figures

Journal ref: Phys. Rev. B 97, 180404R (2018)

Thermal Hall effect and topological edge states in a square lattice antiferromagnet

Authors: Masataka Kawano, Chisa Hotta

Abstract: We show that the two dimensional spatial inversion-symmetry (SIS) broken square lattice antiferromagnet with easy-plane spin anisotropy exhibits a thermal Hall effect and the edge modes characterized by the Z 2 topological invariant. These topological properties require a nonzero Berry curvature, and its origin is ascribed to the Dzyaloshinskii-Moriya (DM) interactions or the noncoplanar magnetic… ▽ More We show that the two dimensional spatial inversion-symmetry (SIS) broken square lattice antiferromagnet with easy-plane spin anisotropy exhibits a thermal Hall effect and the edge modes characterized by the Z 2 topological invariant. These topological properties require a nonzero Berry curvature, and its origin is ascribed to the Dzyaloshinskii-Moriya (DM) interactions or the noncoplanar magnetic ordering generating a U(1) gauge field that couples to the kinetic motion of magnons. Although this picture is established in ferromagnets on the kagome and pyrochlore lattices, it does not apply to our square lattice model since such gauge field cancels out in an edge shared geometry. Instead, our case has an analogy with the anomalous Hall effect of Rashba electronic system where the spin orbit coupling generates an SU(2) gauge field. The two species of magnons defined on antiferromagnetic sublattices can be regarded as the pseudo-spin degrees of freedom of magnons. The DM interactions that emerge due to the SIS breaking serve as a pseudo-spin orbit coupling of magnons and generate a Berry curvature, when the direction of the magnetic moments is properly controlled by the magnetic field. The thermal Hall conductivity of our antiferromagnet shows rapid growth in temperature T beyond the power-law, reflecting the almost gapless low energy branch of the antiferromagnet, which distinctively differs from the T^7/2 -dependence in the pyrochlore ferromagnets. The present system serves as a standard model for the noncentrosymmetric crystals Ba2MnGe2O7 and Ba2CoGe2O7. △ Less

Submitted 21 February, 2019; v1 submitted 15 May, 2018; originally announced May 2018.

Comments: 17 pages, 9 figures

Journal ref: Phys. Rev. B 99, 054422 (2019)

Designing Rashba-Dresselhaus effect in magnetic insulators

Authors: Masataka Kawano, Yoshinori Onose, Chisa Hotta

Abstract: One of the major strategies to control magnetism in spintronics is to utilize the coupling between electron spin and its orbital motion. The Rashba and Dresselhaus spin-orbit couplings induce magnetic textures of band electrons called spin momentum locking, which produces a spin torque by the injection of electric current. However, joule heating had been a bottleneck for device applications. Here,… ▽ More One of the major strategies to control magnetism in spintronics is to utilize the coupling between electron spin and its orbital motion. The Rashba and Dresselhaus spin-orbit couplings induce magnetic textures of band electrons called spin momentum locking, which produces a spin torque by the injection of electric current. However, joule heating had been a bottleneck for device applications. Here, we propose a theory to generate further rich spin textures in insulating antiferromagnets with broken spatial inversion symmetry (SIS), which is easily controlled by a small magnetic field. In antiferromagnets, the ordered moments host two species of magnons that serve as internal degrees of freedom in analogy with electron spins. The Dzyaloshinskii-Moriya interaction introduced by the SIS breaking couples the two-magnon-degrees of freedom with the magnon momentum. We present a systematic way to design such texture and to detect it via magnonic spin current for the realization of antiferromagnetic memory. △ Less

Submitted 8 March, 2019; v1 submitted 10 May, 2018; originally announced May 2018.

Comments: 14 pages, 7 figures

Journal ref: Commun. Phys. 2, 27 (2019)

Kinetic frustration induced supersolid in the $S=1/2$ kagome lattice antiferromagnet in a magnetic field

Authors: Xavier Plat, Tsutomu Momoi, Chisa Hotta

Abstract: We examine instabilities of the plateau phases in the spin-1/2 kagome-lattice antiferromagnet in an applied field by means of degenerate perturbation theory, and find some emergent supersolid phases below the $m=5/9$ plateau. The wave functions of the plateau phases in a magnetic field have the particular construction based on the building blocks of resonating hexagons and their surrounding sites.… ▽ More We examine instabilities of the plateau phases in the spin-1/2 kagome-lattice antiferromagnet in an applied field by means of degenerate perturbation theory, and find some emergent supersolid phases below the $m=5/9$ plateau. The wave functions of the plateau phases in a magnetic field have the particular construction based on the building blocks of resonating hexagons and their surrounding sites. Magnon excitations on each of these blocks suffer from a kinetic frustration effect, namely, they cannot hop easily to the others since the hopping amplitudes through the two paths destructively cancel out with each other. The itineracy is thus weakened, and the system is driven toward the strong coupling regime, which together with the selected paths allowed in real space bears a supersolid phase. This mechanism is contrary to that proposed in lattice-Bose gases, where the strong competing interactions suppress with each other, allowing a small kinetic energy scale to attain the itinerancy. Eventually, we find a supersolid state in which the pattern of resonating hexagons are preserved from the plateau crystal state and only one-third of the originally polarized spins outside the hexagons dominantly join the superfluid component, or equivalently, participate in the magnetization process. △ Less

Submitted 19 July, 2018; v1 submitted 2 April, 2018; originally announced April 2018.

Comments: 16 pages, 11 figures, published version

Journal ref: Phys. Rev. B 98, 014415 (2018)

Pressure-induced two-step spin crossover in double-layered elastic model

Authors: Daisuke Taniguchi, Jun Okabayashi, Chisa Hotta

Abstract: We study the two-step spin crossover in a double-layered elastic model based on transition metal complexes each taking high spin (HS) and low spin (LS) states. Here, only the simplest elastic interactions between adjacent molecules are considered and the system is exposed to the external pressure within the framework of $NPT$-Monte Carlo method. As a certain amount of pressure is applied, the firs… ▽ More We study the two-step spin crossover in a double-layered elastic model based on transition metal complexes each taking high spin (HS) and low spin (LS) states. Here, only the simplest elastic interactions between adjacent molecules are considered and the system is exposed to the external pressure within the framework of $NPT$-Monte Carlo method. As a certain amount of pressure is applied, the first order thermal transition between uniform HS and LS phases transforms to a two-step transition with an emergent intermediate spin (IS) phase, where the HS and LS molecules are paired face to face between layers and form diagonally striped clusters within the layer. The difference in the size of HS and LS molecules is reflected both in the elastic interactions and in the enthalpy, and the IS phase could gain the latter over the loss of the former by significantly reducing its volume. The present pressure effect is interpreted to the chemical one in double-layered transition metal materials, which actually reveals a variety of multistep spin crossover transitions relevant to our numerical result. △ Less

Submitted 1 November, 2017; v1 submitted 9 August, 2017; originally announced August 2017.

Comments: 8 pages, 6 figures, accepted for publication in PRB

Journal ref: Phys. Rev. B 96, 174104 (2017)

Interplay of charge and spin fluctuations of strongly interacting electrons on the kagome lattice

Authors: Frank Pollmann, Krishanu Roychowdhury, Chisa Hotta, Karlo Penc

Abstract: We study electrons hopping on a kagome lattice at third filling described by an extended Hubbard Hamiltonian with on-site and nearest-neighbour repulsions in the strongly correlated limit. As a consequence of the commensurate filling and the large interactions, each triangle has precisely two electrons in the effective low energy description, and these electrons form chains of different lengths. T… ▽ More We study electrons hopping on a kagome lattice at third filling described by an extended Hubbard Hamiltonian with on-site and nearest-neighbour repulsions in the strongly correlated limit. As a consequence of the commensurate filling and the large interactions, each triangle has precisely two electrons in the effective low energy description, and these electrons form chains of different lengths. The effective Hamiltonian includes the ring exchange around the hexagons as well as the nearest-neighbor Heisenberg interaction. Using large scale exact diagonalization, we find that the effective model exhibits two different phases: If the charge fluctuations are small, the magnetic fluctuations confine the charges to short loops around hexagons, yielding a gapped charge ordered phase. When the charge fluctuations dominate, the system undergoes a quantum phase transition to a resonating plaquette phase with ordered spins and gapless spin excitations. We find that a peculiar conservation law is fulfilled: the electron in the chains can be divided into two sublattices, and this division is conserved by the ring exchange term. △ Less

Submitted 31 July, 2014; v1 submitted 20 February, 2014; originally announced February 2014.

Comments: 12 pages, 13 figures, small modifications in the text and in the figures

Journal ref: Phys. Rev. B 90, 035118 (2014)

"Grand Canonical" Finite Size Numerical Approaches : a Route to Measuring Bulk Properties under Applied Field

Authors: Chisa Hotta, Naokazu Shibata

Abstract: We exploit a prescription to observe directly the physical properties of the thermodynamic limit under continuously applied field in one-dimensional quantum finite lattice systems. By systematically scaling down the energy of the Hamiltonian of the open system from center toward both ends, one could adopt the edge sites with negligibly small energy scale as a "grand canonical" small particle bath,… ▽ More We exploit a prescription to observe directly the physical properties of the thermodynamic limit under continuously applied field in one-dimensional quantum finite lattice systems. By systematically scaling down the energy of the Hamiltonian of the open system from center toward both ends, one could adopt the edge sites with negligibly small energy scale as a "grand canonical" small particle bath, and an equilibrium states with non-integer arbitrary conserved numbers, e.g., electron numbers or sz, are realized in the main part of the system. This will enable the evaluation of the response functions under continuously varying external field in a small lattice without any fine tuning or scaling of parameters while keeping the standard numerical accuracy. Demonstrations are given on quantum spin systems as well as on a Hubbard model by the density matrix renormalization group calculation. △ Less

Submitted 14 July, 2013; originally announced July 2013.

Comments: 5pages, 3 figures

Journal ref: Phys. Rev. B 86, 041108 (2012)

Controlling frustrated liquids and solids with an applied field in a kagome Heisenberg antiferromagnet

Authors: Satoshi Nishimoto, Naokazu Shibata, Chisa Hotta

Abstract: Quantum spin-1/2 kagome Heisenberg antiferromagnet is the representative frustrated system possibly hosting a spin liquid. Clarifying the nature of this elusive topological phase is a key challenge in condensed matter, however, even identifying it still remains unsettled. Here, we apply a magnetic field and discover a series of spin gapped phases appearing at five different fractions of magnetizat… ▽ More Quantum spin-1/2 kagome Heisenberg antiferromagnet is the representative frustrated system possibly hosting a spin liquid. Clarifying the nature of this elusive topological phase is a key challenge in condensed matter, however, even identifying it still remains unsettled. Here, we apply a magnetic field and discover a series of spin gapped phases appearing at five different fractions of magnetization by means of grand canonical density matrix renormalization group, an unbiased state-of-art numerical technique. The magnetic field dopes magnons and first gives rise to a possible Z3 spin liquid plateau at 1/9-magnetization. Higher field induces a self-organized super-lattice-unit, a six-membered ring of quantum spins, resembling an atomic orbital structure. Putting magnons into this unit one by one yields three quantum solid plateaus. We thus find that the magnetic field could control the transition between various emergent phases by continuously releasing the frustration. △ Less

Submitted 17 July, 2013; v1 submitted 14 July, 2013; originally announced July 2013.

Comments: 6 pages 4 figures and 4 page supplementary file. Nature Communications, accepted, July 10, 2013

Multiferroicity in an organic charge-transfer salt: Electric-dipole-driven magnetism

Authors: P. Lunkenheimer, J. Müller, S. Krohns, F. Schrettle, A. Loidl, B. Hartmann, R. Rommel, M. de Souza, C. Hotta, J. A. Schlueter, M. Lang

Abstract: Multiferroics, showing simultaneous ordering of electrical and magnetic degrees of freedom, are remarkable materials as seen from both the academic and technological points of view. A prominent mechanism of multiferroicity is the spin-driven ferroelectricity, often found in frustrated antiferromagnets with helical spin order. There, similar to conventional ferroelectrics, the electrical dipoles ar… ▽ More Multiferroics, showing simultaneous ordering of electrical and magnetic degrees of freedom, are remarkable materials as seen from both the academic and technological points of view. A prominent mechanism of multiferroicity is the spin-driven ferroelectricity, often found in frustrated antiferromagnets with helical spin order. There, similar to conventional ferroelectrics, the electrical dipoles arise from an off-centre displacement of ions. However, recently a different mechanism, namely purely electronic ferroelectricity, where charge order breaks inversion symmetry, has attracted considerable interest. Here we provide evidence for this exotic type of ferroelectricity, accompanied by antiferromagnetic spin order, in a two-dimensional organic charge-transfer salt, thus representing a new class of multiferroics. Quite unexpectedly for electronic ferroelectrics, dipolar and spin order arise nearly simultaneously. This can be ascribed to the loss of spin frustration induced by the ferroelectric ordering. Hence, here the spin order is driven by the ferroelectricity, in marked contrast to the spin-driven ferroelectricity in helical magnets. △ Less

Submitted 24 November, 2011; v1 submitted 11 November, 2011; originally announced November 2011.

Comments: 8 pages, 9 figures (including 4 pages and 6 figures in supplementary information). Version 2 with minor errors corrected (legend of Fig. 3c and definition of vectors e and Q)

Journal ref: Nature Mater. 11 (2012) 755

Wilson ratio of a Tomonaga-Luttinger liquid in a spin-1/2 Heisenberg ladder

Authors: K. Ninios, Tao Hong, T. Manabe, C. Hotta, S. N. Herringer, M. M. Turnbull, C. P. Landee, Y. Takano, H. B. Chan

Abstract: Using micromechanical force magnetometry, we have measured the magnetization of the strong-leg spin-1/2 ladder compound (C$_7$H$_{10}$N)$_2$CuBr$_2$ at temperatures down to 45 mK. Low-temperature magnetic susceptibility as a function of field exhibits a maximum near the critical field H_c at which the magnon gap vanishes, as expected for a gapped one-dimensional antiferromagnet. Above H_c a clear… ▽ More Using micromechanical force magnetometry, we have measured the magnetization of the strong-leg spin-1/2 ladder compound (C$_7$H$_{10}$N)$_2$CuBr$_2$ at temperatures down to 45 mK. Low-temperature magnetic susceptibility as a function of field exhibits a maximum near the critical field H_c at which the magnon gap vanishes, as expected for a gapped one-dimensional antiferromagnet. Above H_c a clear minimum appears in the magnetization as a function of temperature as predicted by theory. In this field region, the susceptibility in conjunction with our specific heat data yields the Wilson ratio R_W. The result supports the relation R_W=4K, where K is the Tomonaga-Luttinger-liquid parameter. △ Less

Submitted 18 March, 2012; v1 submitted 25 October, 2011; originally announced October 2011.

Journal ref: Phys. Rev. Lett. 108, 097201 (2012)

Boundary effects in the density-matrix renormalization group calculation

Authors: Naokazu Shibata, Chisa Hotta

Abstract: We investigate the boundary effect of the density matrix renormalization group calculation (DMRG), which is an artifactual induction of symmetry-breaking pseudo-long-range order and takes place when the long-range quantum fluctuation cannot be properly included in the variational wave function due to numerical limitation. The open boundary condition often used in DMRG suffers from the boundary eff… ▽ More We investigate the boundary effect of the density matrix renormalization group calculation (DMRG), which is an artifactual induction of symmetry-breaking pseudo-long-range order and takes place when the long-range quantum fluctuation cannot be properly included in the variational wave function due to numerical limitation. The open boundary condition often used in DMRG suffers from the boundary effect the most severely, which is directly reflected in the distinct spatial modulations of the local physical quantity. By contrast, the other boundary conditions such as the periodic one or the sin^2-deformed interaction [A. Gendiar, R. Krcmar, and T. Nishino, Prog. Theor. Phys. {\bf 122}, 953 (2009)] keep spatial homogeneity, and are relatively free from the boundary effect. By comparing the numerical results of those various boundary conditions, we show that the open boundary condition sometimes gives unreliable results even after the finite-size scaling. We conclude that the examination of the boundary condition dependence is required besides the usual treatment based on the system size or accuracy dependence in cases where the long-range quantum fluctuation is important. △ Less

Submitted 1 September, 2011; v1 submitted 30 June, 2011; originally announced June 2011.

Comments: 9 pages, 8 figures, to appear in PRB

Journal ref: Phys. Rev. B 84, 115116 (2011)

Designing Dirac points in two-dimensional lattices

Authors: Kenichi Asano, Chisa Hotta

Abstract: We present a framework to elucidate the existence of accidental contacts of energy bands, particularly those called Dirac points which are the point contacts with linear energy dispersions in their vicinity. A generalized von-Neumann-Wigner theorem we propose here gives the number of constraints on the lattice necessary to have contacts without fine tuning of lattice parameters. By counting this n… ▽ More We present a framework to elucidate the existence of accidental contacts of energy bands, particularly those called Dirac points which are the point contacts with linear energy dispersions in their vicinity. A generalized von-Neumann-Wigner theorem we propose here gives the number of constraints on the lattice necessary to have contacts without fine tuning of lattice parameters. By counting this number, one could quest for the candidate of Dirac systems without solving the secular equation. The constraints can be provided by any kinds of symmetry present in the system. The theory also enables the analytical determination of k-point having accidental contact by selectively picking up only the degenerate solution of the secular equation. By using these frameworks, we demonstrate that the Dirac points are feasible in various two-dimensional lattices, e.g. the anisotropic Kagome lattice under inversion symmetry is found to have contacts over the whole lattice parameter space. Spin-dependent cases, such as the spin-density-wave state in LaOFeAs with reflection symmetry, are also dealt with in the present scheme. △ Less

Submitted 20 April, 2011; v1 submitted 13 April, 2011; originally announced April 2011.

Comments: 15pages, 9figures (accepted to Phys. Rev. B)

Field-induced Tomonaga-Luttinger liquid phase of a two-leg spin-1/2 ladder with strong leg interactions

Authors: Tao Hong, Y. H. Kim, C. Hotta, Y. Takano, G. Tremelling, M. M. Turnbull, C. P. Landee, H. -J. Kang, N. B. Christensen, K. Lefmann, K. P. Schmidt, G. S. Uhrig, C. Broholm

Abstract: We study the magnetic-field-induced quantum phase transition from a gapped quantum phase that has no magnetic long-range order into a gapless phase in the spin-1/2 ladder compound bis(2,3-dimethylpyridinium) tetrabromocuprate (DIMPY). At temperatures below about 1 K, the specific heat in the gapless phase attains an asymptotic linear temperature dependence, characteristic of a Tomonaga-Luttinger l… ▽ More We study the magnetic-field-induced quantum phase transition from a gapped quantum phase that has no magnetic long-range order into a gapless phase in the spin-1/2 ladder compound bis(2,3-dimethylpyridinium) tetrabromocuprate (DIMPY). At temperatures below about 1 K, the specific heat in the gapless phase attains an asymptotic linear temperature dependence, characteristic of a Tomonaga-Luttinger liquid. Inelastic neutron scattering and the specific heat measurements in both phases are in good agreement with theoretical calculations, demonstrating that DIMPY is the first model material for an S=1/2 two-leg spin ladder in the strong-leg regime. △ Less

Submitted 27 August, 2010; v1 submitted 23 August, 2010; originally announced August 2010.

Comments: 4.1 pages, 4 figures (Fig. 4 updated), to appear in Physical Review Letters

Journal ref: Phys. Rev. Lett. 105, 137207 (2010)

Anomalous dielectric response in the dimer Mott insulator $κ$-(BEDT-TTF)$_2$Cu$_2$(CN)$_3$

Authors: Majed Abdel-Jawad, Ichiro Terasaki, Takahiko Sasaki, Naoki Yoneyama, Norio Kobayashi, Yoshiaki Uesu, Chisa Hotta

Abstract: We have measured and analyzed the dielectric constant of the dimer Mott insulator $κ-$(BEDT-TTF)$_2$Cu$_2$(CN)$_3$, which is known as a playground for a spin-liquid state. Most unexpectedly, this particular organic salt has nontrivial charge degrees of freedom, being characterized by a dipole-glass-like dielectric relaxation below around 60 K. This is ascribed to the charge disproportionation with… ▽ More We have measured and analyzed the dielectric constant of the dimer Mott insulator $κ-$(BEDT-TTF)$_2$Cu$_2$(CN)$_3$, which is known as a playground for a spin-liquid state. Most unexpectedly, this particular organic salt has nontrivial charge degrees of freedom, being characterized by a dipole-glass-like dielectric relaxation below around 60 K. This is ascribed to the charge disproportionation within the dimer due to the intersite Coulomb repulsion. In this respect, the dipole observed here is purely quantum mechanical, and thereby exhibits unexpected responses. △ Less

Submitted 19 July, 2010; v1 submitted 19 March, 2010; originally announced March 2010.

Comments: 5 pages, 5 figures

Interplay of strongly correlated electrons and localized Ising moments in one-dimension

Authors: Chisa Hotta

Abstract: We study the ground state properties of the one-dimensional quarter-filled strongly correlated electronic chain coupled by $J$ to another chain of antiferromagnetic Ising moments. We focus on the case where the large Coulomb interactions localize the charges on every other site. Both the electronic spins and the Ising moments interact antiferromagnetically within each chain by $J_{\rm eff}$ and… ▽ More We study the ground state properties of the one-dimensional quarter-filled strongly correlated electronic chain coupled by $J$ to another chain of antiferromagnetic Ising moments. We focus on the case where the large Coulomb interactions localize the charges on every other site. Both the electronic spins and the Ising moments interact antiferromagnetically within each chain by $J_{\rm eff}$ and $J'$, respectively. Since the number of electrons is half as that of the Ising moments the period of magnetic correlation of these two chains are incommensurate. In the presence of $J$, the frustration of $J_{\rm eff}$ and $J'$ arises, which may lead the system to the intriguing magneto-electric effect. △ Less

Submitted 18 December, 2009; originally announced December 2009.

Comments: 8pages 6figures

"Nonmagnetic" Dimer Mott insulator induced by the coupling of spin with quantum electric dipoles

Authors: Chisa Hotta

Abstract: We present an effective dipolar-spin model based on the strong coupling analysis, which may explain the possible origin of the "spin liquid Mott insulator". The issue is related to the dimer Mott insulator reminiscent of the organic triangular lattice system, $κ$-ET$_2$Cu$_2$(CN)$_3$, whose gapless spin liquid state had been discussed in the context of geometrical frustration of exchange coupling,… ▽ More We present an effective dipolar-spin model based on the strong coupling analysis, which may explain the possible origin of the "spin liquid Mott insulator". The issue is related to the dimer Mott insulator reminiscent of the organic triangular lattice system, $κ$-ET$_2$Cu$_2$(CN)$_3$, whose gapless spin liquid state had been discussed in the context of geometrical frustration of exchange coupling, $J$, between spins on dimer orbitals. It turns out that another degrees of freedom within the insulator, {\it quantum electric-dipoles} on dimers, interacts with each other by the Coulomb interaction and brings about a significant suppression of $J$ through the dipolar-spin coupling. △ Less

Submitted 26 August, 2010; v1 submitted 18 December, 2009; originally announced December 2009.

Comments: 4pages 2figures

Non-local network-excitation in the Ising antiferromagnet on the anisotropic triangular lattice

Authors: Chisa Hotta, Tetsuhiro Kiyota, Nobuo Furukawa

Abstract: We study the Ising antiferromagnet on the triangular lattice which are interacting along three different directions via two $J$-bonds and one anisotropic $J'$-bond with $J' \ge J$. Although its finite temperature state has long been considered as simply disordered, we find a systematic generation of a ferromagnetic $J'$-bonds named "good defects" each carrying the energy $2δ_J\equiv 2(J'-J)$. Th… ▽ More We study the Ising antiferromagnet on the triangular lattice which are interacting along three different directions via two $J$-bonds and one anisotropic $J'$-bond with $J' \ge J$. Although its finite temperature state has long been considered as simply disordered, we find a systematic generation of a ferromagnetic $J'$-bonds named "good defects" each carrying the energy $2δ_J\equiv 2(J'-J)$. They exhibit an eminent correlation which is regarded as a network or a soft lattice. The specific heat shows a universal peak at $T^* \sim δ_J$ and the corresponding energy follows the activation type of temperature-dependence with an excitation gap $2δ_J$. Energetics and the crossover between this particular two-dimensional state and the other disordered or one-dimensional states are discussed. △ Less

Submitted 3 December, 2009; originally announced December 2009.

Comments: 4pages 4figures

Density-matrix renormalization study of the frustrated fermions on the triangular lattice

Authors: S. Nishimoto, C. Hotta

Abstract: We show that the two-dimensional density-matrix renormalization analysis is useful to detect the symmetry breaking in the fermionic model on a triangular lattice. Under the cylindrical boundary conditions with chemical potentials on edge sites, we find that the open edges work as perturbation to select the strongest correlations {\it only in the presence of a long range order}. We also demonstra… ▽ More We show that the two-dimensional density-matrix renormalization analysis is useful to detect the symmetry breaking in the fermionic model on a triangular lattice. Under the cylindrical boundary conditions with chemical potentials on edge sites, we find that the open edges work as perturbation to select the strongest correlations {\it only in the presence of a long range order}. We also demonstrate that the ordinary size scaling analysis on the charge gap as well as that of the local charge density under this boundary condition could determine the metal-insulator phase boundary, which scales almost perfectly with the density of states and the exact solutions in the weak and strong coupling region, respectively. △ Less

Submitted 7 February, 2009; originally announced February 2009.

Comments: 5 pages, 7 figures

Variational Monte Carlo Study of a Spinless Fermion t-V Model on a Triangular Lattice: Formation of a Pinball Liquid

Authors: M. Miyazaki, C. Hotta, S. Miyahara, K. Matsuda, N. Furukawa

Abstract: We analyze a model of spinless fermions on a triangular lattice at half-filling interacting via strong nearest-neighbor repulsive interactions, V, using the variational Monte Carlo simulation technique. The existence of three-sublattice long-range order is confirmed by the finite-size scaling analysis of the charge structural factor at V_c/t > 12. This ordered phase shows characteristics expecte… ▽ More We analyze a model of spinless fermions on a triangular lattice at half-filling interacting via strong nearest-neighbor repulsive interactions, V, using the variational Monte Carlo simulation technique. The existence of three-sublattice long-range order is confirmed by the finite-size scaling analysis of the charge structural factor at V_c/t > 12. This ordered phase shows characteristics expected for a so called "pinball liquid" state, which has the spontaneous separation of fermionic degrees of freedom into coexisting Wigner crystal-like charge order (pin) and a metal (ball). The pins are fixed in order to maximize the kinetic energy gain of balls which move almost freely. The Fermi surface is reconstructed at V=V_c and remains towards the strong coupling limit. These features reminiscent of the strong correlation together with the large value of V_c/t distinguishes the pinball liquid from the conventional charge-density-wave. △ Less

Submitted 23 July, 2008; originally announced July 2008.

Comments: 7 pages, 7 figures