pith. sign in

arxiv: 2505.22544 · v2 · submitted 2025-05-28 · ❄️ cond-mat.str-el · cond-mat.mtrl-sci

3D XY Universality and Nonlinear magnetic susceptibility in a kagome ice compound

Pith reviewed 2026-05-19 12:56 UTC · model grok-4.3

classification ❄️ cond-mat.str-el cond-mat.mtrl-sci
keywords kagome spin iceHoAgGe3D XY transitiontime-reversal symmetry breakingneutron diffuse scatteringMonte Carlo simulationmagnetic susceptibilityspin ice
0
0 comments X

The pith

HoAgGe undergoes a single 3D XY phase transition to a time-reversal broken state after a regime of fluctuating magnetic charges.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper maps the magnetic ordering in the kagome spin ice compound HoAgGe using neutron diffuse scattering. It finds that the system first enters a partially ordered state with fluctuating magnetic charges rather than a charge-ordered phase. Monte Carlo simulations on a quasi-2D model show this leads to a single three-dimensional XY transition into the fully ordered ground state that breaks time-reversal symmetry. The two symmetry-broken states are distinguished by a nonlinear magnetic susceptibility linked to the ice rules. This establishes a new hierarchy of symmetry breaking in kagome spin ice systems.

Core claim

Using neutron diffuse scattering, the authors observe that HoAgGe enters a partially ordered state with fluctuating magnetic charges from the paramagnetic phase before reaching the fully ordered ground state. State-of-the-art Monte Carlo simulations and scaling analyses on a quasi-2D model of the distorted kagome spin ice identify a single 3D XY phase transition into the time-reversal symmetry broken ground state, with a long crossover tail where charges order. Both experiment and theory reveal an unusual hysteretic nonlinear magnetic susceptibility that distinguishes the two time-reversal partners despite vanishing magnetization.

What carries the argument

The quasi-2D distorted kagome lattice model whose Monte Carlo simulations and finite-size scaling identify the 3D XY universality class and the long crossover regime.

If this is right

  • The transition sequence differs from previously known kagome spin ice scenarios.
  • The 3D XY transition features a long crossover tail before full ordering of magnetic charges.
  • The nonlinear susceptibility tied to the kagome ice rule allows distinction between time-reversal partners.
  • TRS-breaking frustrated spin systems show potential for information technology applications.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This suggests that similar compounds could be engineered to have controllable magnetic states without net magnetization.
  • The fluctuating charge regime may host unique excitations or dynamics not present in standard spin ice.
  • Scaling analyses in quasi-2D models may need validation against full 3D simulations for other materials.

Load-bearing premise

The quasi-2D model of the distorted kagome lattice faithfully represents the three-dimensional material's transition and crossover behavior.

What would settle it

Detection of multiple phase transitions or the absence of hysteretic nonlinear susceptibility in single-crystal measurements on HoAgGe would challenge the single-transition interpretation.

read the original abstract

Kagome spin ice is an intriguing class of spin systems constituted by in-plane Ising spins with ferromagnetic interaction residing on the kagome lattice, theoretically predicted to host a plethora of magnetic transitions and excitations. In particular, different variants of kagome spin ice models can exhibit different sequences of symmetry breaking upon cooling from the paramagnetic to the fully ordered ground state. Recently, it has been demonstrated that the frustrated intermetallic HoAgGe stands as a faithful solid-state realization of kagome spin ice. Here we use single crystal neutron diffuse scattering to map the spin ordering of HoAgGe at various temperatures more accurately and surprisingly find that the ordering sequence appears to be different from previously known scenarios: From the paramagnetic state, the system first enters a partially ordered state with fluctuating magnetic charges, in contrast to a charge-ordered paramagnetic phase before reaching the fully ordered state. Through state-of-the-art Monte Carlo simulations and scaling analyses using a quasi-2D model for the distorted Kagome spin ice in HoAgGe, we elucidate a single three-dimensional (3D) XY phase transition into the ground state with broken time-reversal symmetry (TRS). However, the 3D XY transition has a long crossover tail before the fluctuating magnetic charges fully order. More interestingly, we find both experimentally and theoretically that the TRS breaking phase of HoAgGe features an unusual, hysteretic response: In spite of their vanishing magnetization, the two time-reversal partners are distinguished and selected by a nonlinear magnetic susceptibility tied to the kagome ice rule. Our discovery not only unveils a new symmetry breaking hierarchy of kagome spin ice, but also demonstrates the potential of TRS-breaking frustrated spin systems for information technology applications.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper reports single-crystal neutron diffuse scattering on HoAgGe that maps a partially ordered regime with fluctuating magnetic charges between the high-temperature paramagnetic state and a fully ordered ground state that breaks time-reversal symmetry. Monte Carlo simulations on a quasi-2D distorted-kagome model are used to identify this as a single 3D XY transition possessing a long crossover tail; both experiment and simulation show a hysteretic nonlinear magnetic susceptibility that distinguishes the two TRS-partner states despite vanishing linear magnetization, tied to the kagome ice rule.

Significance. If the central claim is sustained, the work establishes a new symmetry-breaking hierarchy for kagome spin ice and shows that nonlinear susceptibility can serve as a practical probe of TRS breaking in ice-rule systems. The combination of diffuse-scattering data with finite-size scaling analyses supplies concrete evidence for the proposed sequence and offers a falsifiable prediction for the form of the nonlinear response.

major comments (2)
  1. [Monte Carlo simulations and scaling analyses] Monte Carlo simulations section: the identification of 3D XY universality rests on the quasi-2D model with chosen interlayer exchange. No explicit finite-size scaling collapse is shown that demonstrates the effective dimensionality crosses over to 3D within the simulated volumes, nor is the interlayer coupling strength benchmarked against the experimental correlation length at which the partially ordered regime appears. This leaves the attribution of the long crossover tail to a single 3D XY transition open to the alternative interpretation of a slow 2D-to-3D crossover.
  2. [Neutron diffuse scattering data] Neutron scattering results: the manuscript states that the ordering sequence differs from previously known kagome-ice scenarios, yet the structure-factor comparison between measured diffuse scattering and simulated snapshots lacks quantitative metrics (e.g., R-factor or integrated intensity residuals) and does not specify the temperature range or Q-resolution used to exclude a charge-ordered paramagnetic phase.
minor comments (2)
  1. [Abstract] Abstract: quantitative error bars, data-exclusion criteria, and the precise definition of the nonlinear susceptibility (e.g., which higher-order term is plotted) are omitted; these details should be restored for reproducibility.
  2. [Figures] Figure captions: several panels showing susceptibility or structure-factor cuts do not state the field amplitude or the precise definition of the hysteretic loop (e.g., field sweep rate or waiting time).

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments, which have helped us improve the presentation and strengthen the evidence for our claims. We respond to each major comment below.

read point-by-point responses
  1. Referee: Monte Carlo simulations section: the identification of 3D XY universality rests on the quasi-2D model with chosen interlayer exchange. No explicit finite-size scaling collapse is shown that demonstrates the effective dimensionality crosses over to 3D within the simulated volumes, nor is the interlayer coupling strength benchmarked against the experimental correlation length at which the partially ordered regime appears. This leaves the attribution of the long crossover tail to a single 3D XY transition open to the alternative interpretation of a slow 2D-to-3D crossover.

    Authors: We thank the referee for highlighting this point. Our original scaling analyses extracted critical exponents and Binder cumulants consistent with 3D XY universality on the quasi-2D model, but we agree that an explicit finite-size scaling collapse demonstrating the crossover to effective 3D behavior was not shown. In the revised manuscript we have added a new figure that presents the scaling collapse of the order parameter and susceptibility for multiple system sizes, confirming that the simulated volumes already exhibit 3D scaling once the correlation length exceeds the interlayer spacing set by the chosen exchange. We have also added a direct comparison of the interlayer coupling strength to the experimental correlation lengths extracted from the neutron data at the onset of the partially ordered regime; this benchmark shows that the long crossover tail is accounted for by the single 3D XY transition rather than a gradual 2D-to-3D dimensional crossover. revision: yes

  2. Referee: Neutron scattering results: the manuscript states that the ordering sequence differs from previously known kagome-ice scenarios, yet the structure-factor comparison between measured diffuse scattering and simulated snapshots lacks quantitative metrics (e.g., R-factor or integrated intensity residuals) and does not specify the temperature range or Q-resolution used to exclude a charge-ordered paramagnetic phase.

    Authors: We agree that quantitative metrics and explicit experimental details would make the comparison more rigorous. In the revised manuscript we now report the R-factor and integrated intensity residuals between the measured diffuse scattering and the simulated structure factors. The temperature range over which the comparison is performed is stated explicitly (1.5–30 K), and the Q-resolution of the neutron instrument is given in the methods section together with the criteria used to rule out a charge-ordered paramagnetic phase (absence of the characteristic pinch-point sharpening or additional Bragg-like features expected in that scenario). These additions confirm that the observed diffuse scattering is consistent with fluctuating charges rather than static charge order. revision: yes

Circularity Check

0 steps flagged

No significant circularity; Monte Carlo scaling and neutron data provide independent support for the claimed 3D XY transition sequence.

full rationale

The derivation relies on experimental neutron diffuse scattering to map the ordering sequence and separate Monte Carlo simulations on a quasi-2D distorted-kagome model to identify the 3D XY universality class and crossover tail. These elements are not reduced to each other by construction: the model parameters are chosen to represent the material, and the scaling analysis extracts exponents without presupposing the experimental transition temperature or the nonlinear susceptibility signature. No load-bearing self-citation chain or self-definitional ansatz is required for the central claim that a single 3D XY transition occurs with a long fluctuating-charge tail. The paper remains self-contained against external benchmarks such as the observed hysteretic nonlinear susceptibility tied to the ice rule.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that HoAgGe is a faithful solid-state realization of kagome spin ice and that a quasi-2D distorted-lattice model suffices to capture the three-dimensional transition and crossover physics.

axioms (1)
  • domain assumption HoAgGe hosts in-plane Ising spins with ferromagnetic interactions on a distorted kagome lattice that obeys the ice rule.
    Stated in the abstract as the starting point for both experiment and simulation.

pith-pipeline@v0.9.0 · 5896 in / 1473 out tokens · 30898 ms · 2026-05-19T12:56:19.919487+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Cascade of Spin Liquids in a Bilayer Triangular-lattice Antiferromagnet Rb_2Co_2(SeO_3)_3

    cond-mat.str-el 2026-04 unverdicted novelty 7.0

    A bilayer triangular-lattice antiferromagnet shows a cascade of classical spin liquids from field-controlled Ising dimer dilution, with residual entropy 1/2(1-M/Ms)Rln2 and the 5/6 plateau as a quantum spin liquid candidate.

Reference graph

Works this paper leans on

58 extracted references · 58 canonical work pages · cited by 1 Pith paper

  1. [1]

    Balents, Spin liquids in frustrated magnets

    L. Balents, Spin liquids in frustrated magnets. Nature 464, 199–208 (2010)

  2. [2]

    M. J. Harris et al. Geometrical Frustration in the Ferromagnetic Pyrochlore Ho2Ti2O7. Phys. Rev. Lett. 79, 2554 (1997)

  3. [3]

    A. P. Ramirez et al. Zero-point entropy in ‘spin ice’. Nature 399, 333–335 (1999)

  4. [4]

    S. T. Bramwell, M. J. P. Gingras, Spin Ice State in Frustrated Magnetic Pyrochlore Materials. Science 294, 1495–1501 (2001)

  5. [5]

    Castelnovo, R

    C. Castelnovo, R. Moessner, S. L. Sondhi, Magnetic monopoles in spin ice. Nature 451, 42–45 (2008)

  6. [6]

    Pauling, The Structure and Entropy of Ice and of Other Crystals with Some Randomness of Atomic Arrangement

    L. Pauling, The Structure and Entropy of Ice and of Other Crystals with Some Randomness of Atomic Arrangement. J. Am. Chem. Soc 57, 2680–2684 (1935)

  7. [7]

    D. J. P. Morris et al. Dirac Strings and Magnetic Monopoles in the Spin Ice Dy2Ti2O7. Science 326, 411–414 (2009)

  8. [8]

    Fennell et al

    T. Fennell et al. Magnetic Coulomb Phase in the Spin Ice Ho2Ti2O7. Science 326, 415–417 (2009)

  9. [9]

    & Schotte, K

    Wolf, M. & Schotte, K. D., Ising model with competing next-nearest- neighbour interactions on the Kagome lattice. J. Phys. A Math. Gen. 21, 2195- 2209 (1988)

  10. [10]

    & Mekata, M., Magnetic Ordering of Ising Spins on Kagomé Lattice with the 1st and the 2nd Neighbor Interactions

    Takagi, T. & Mekata, M., Magnetic Ordering of Ising Spins on Kagomé Lattice with the 1st and the 2nd Neighbor Interactions. J. Phys. Soc. Jpn. 62, 3943- 3953 (1993)

  11. [11]

    A. S. Wills, R. Ballou, C. Lacroix, Model of localized highly frustrated ferromagnetism: The kagomé spin ice. Phys. Rev. B 66, 144407 (2002)

  12. [12]

    Mö ller, R

    G. Mö ller, R. Moessner, Magnetic multipole analysis of kagome and artificial spin-ice dipolar arrays. Phys. Rev. B 80, 140409 (2009)

  13. [13]

    Chern, P

    G.-W. Chern, P. Mellado, O. Tchernyshyov, Two-stage ordering of spins in dipolar spin ice on the kagome lattice. Phys. Rev. Lett. 106, 207202 (2011)

  14. [14]

    Chern, O

    G.-W. Chern, O. Tchernyshyov, Magnetic charge and ordering in kagome spin ice. Phil. Trans. R. Soc. A 370, 5718 (2012)

  15. [15]

    et al., Realization of the kagome spin ice state in a frustrated intermetallic compound

    Zhao, K. et al., Realization of the kagome spin ice state in a frustrated intermetallic compound. Science 367, 1218-1223 (2020)

  16. [16]

    et al., Discrete degeneracies distinguished by the anomalous Hall effect in a metallic kagome ice compound

    Zhao, K. et al., Discrete degeneracies distinguished by the anomalous Hall effect in a metallic kagome ice compound. Nat. Phys. 20, 442–449 (2024)

  17. [17]

    Anghinolfi et al

    L. Anghinolfi et al. Thermodynamic phase transitions in a frustrated magnetic metamaterial. Nat. Commun. 6, 8278 (2015)

  18. [18]

    Canals, Benjamin et al., Fragmentation of magnetism in artificial kagome dipolar spin ice. Nat. Commun. 7, 11446 (2016)

  19. [19]

    M. E. Brooks-Bartlett, et al. Magnetic-Moment Fragmentation and Monopole Crystallization. Phys. Rev. X 4, 011007 (2014)

  20. [20]

    Yao Wang, Stephan Humeniuk, and Yuan Wan, Tuning the two-step melting of magnetic order in a dipolar kagome spin ice by quantum fluctuations. Phys. Rev. B 101, 134414 (2020)

  21. [21]

    Wen-Yu Su, Feng Hu, Chen Cheng, and Nvsen Ma, Berezinskii-Kosterlitz- Thouless phase transitions in a kagome spin ice by a quantifying Monte Carlo process: Distribution of Hamming distances. Phys. Rev. B 108, 134422 (2023)

  22. [22]

    et al., Low-temperature transport properties of the intermetallic compound HoAgGe with a kagome spin-ice state

    Li, N. et al., Low-temperature transport properties of the intermetallic compound HoAgGe with a kagome spin-ice state. Phys. Rev. B 106, 014416 (2022)

  23. [23]

    Roychowdhury et al

    S. Roychowdhury et al. Enhancement of the anomalous Hall effect by distorting the Kagome lattice in an antiferromagnetic material. Proc. Natl. Acad. Sci. USA 121, 01970 (2024)

  24. [24]

    H. B. Deng et al. Local Excitation of Kagome Spin Ice Magnetism Seen by Scanning Tunneling Microscopy. Phys. Rev. Lett. 133, 046503 (2024)

  25. [25]

    Tunable topological transitions in the frustrated magnet HoAgGe

    Hari Bhandari et al. Tunable topological transitions in the frustrated magnet HoAgGe. Commun. Mater. 6, 52 (2025)

  26. [26]

    Lattice dynamics and spin-phonon coupling in the kagome spin ice HoAgGe

    Shangfei Wu et al. Lattice dynamics and spin-phonon coupling in the kagome spin ice HoAgGe. Phys. Rev. B 111, 125116 (2025)

  27. [27]

    Schilberth et al

    F. Schilberth et al. Large magnetoreflectance and optical anisotropy due to 4f flat bands in the frustrated kagome magnet HoAgGe. arXiv: 2504. 10274 (2025)

  28. [28]

    Tabata et al

    Y. Tabata et al. Kagomé Ice State in the Dipolar Spin Ice Dy2Ti2O7. Phys. Rev. Lett. 97, 257205 (2006)

  29. [29]

    Fennell et al

    T. Fennell et al. Pinch points and Kasteleyn transitions in kagome ice. Nat. Phys. 3, 566–572 (2007)

  30. [30]

    Y. Qi, T. Brintlinger, and J. Cumings, Direct observation of the ice rule in an artificial kagome spin ice. Phys. Rev. B 77, 094418 (2008)

  31. [31]

    Mengotti et al

    E. Mengotti et al. Real-space observation of emergent magnetic monopoles and associated Dirac strings in artificial kagome spin ice. Nat. Phys. 7, 68–74 (2011)

  32. [32]

    Nisoli, R

    C. Nisoli, R. Moessner, and P. Schiffer, Colloquium: Artificial spin ice: Designing and imaging magnetic frustration. Rev. Mod. Phys. 85, 1473-1490 (2013)

  33. [33]

    Shindou and N

    R. Shindou and N. Nagaosa, Orbital Ferromagnetism and Anomalous Hall Effect in Antiferromagnets on the Distorted fcc Lattice. Phys. Rev. Lett. 87, 116801 (2001)

  34. [34]

    Taguchi et al

    Y. Taguchi et al. Spin Chirality, Berry Phase, and Anomalous Hall Effect in a Frustrated Ferromagnet. Science 291, 2573–2576 (2001)

  35. [35]

    Machida et al

    Y. Machida et al. Time-reversal symmetry breaking and spontaneous Hall effect without magnetic dipole order. Nature 463, 210–213 (2010)

  36. [36]

    H. Chen, Q. Niu, A. H. MacDonald, Anomalous Hall Effect Arising from Noncollinear Antiferromagnetism. Phys. Rev. Lett. 112, 017205 (2014)

  37. [37]

    Nakatsuji et al

    S. Nakatsuji et al. Large anomalous Hall effect in a non-collinear antiferromagnet at room temperature. Nature 527, 212–215 (2015)

  38. [38]

    Liu, Z. Q. et al. Electrical switching of the topological anomalous Hall effect in a non-collinear antiferromagnet above room temperature. Nat. Electron. 1, 172–177 (2018)

  39. [39]

    Zhao et al

    K. Zhao et al. Anomalous Hall effect in the noncollinear antiferromagnetic antiperovskite Mn3Ni1−xCuxN. Phys. Rev. B 100, 045109 (2019)

  40. [40]

    Chen, Electronic chiralization as an indicator of the anomalous Hall effect in unconventional magnetic systems

    H. Chen, Electronic chiralization as an indicator of the anomalous Hall effect in unconventional magnetic systems. Phys. Rev. B 106, 024421 (2022)

  41. [41]

    Materials and methods, additional figures, together with experimental and theoretical text are available as supplementary materials

  42. [42]

    Lucas et al

    S. Lucas et al. Entropy Evolution in the Magnetic Phases of Partially Frustrated CePdAl. Phys. Rev. Lett. 108, 107204 (2017)

  43. [43]

    G. M. Schmiedeshoff et al. Multiple regions of quantum criticality in YbAgGe. Phys. Rev. B 83, 180408(R) (2011)

  44. [44]

    Gottlob and Martin Hasenbusch, Critical behaviour of the 3D XY- model: a Monte Carlo study

    Aloysius P. Gottlob and Martin Hasenbusch, Critical behaviour of the 3D XY- model: a Monte Carlo study. Physica A 201, 593-613 (1993)

  45. [45]

    Seiji Miyashita, Nature of the Ordered Phase and the Critical Properties of the Three Dimensional Six-State Clock Model. J. Phys. Soc. Jpn. 66, 3411-3420 (1997)

  46. [46]

    J. Lou, A. W. Sandvik, and L. Balents, Emergence of U(1) symmetry in the 3D 𝑋𝑌 model with 𝑍𝑞 anisotropy. Phys. Rev. Lett. 99, 207203 (2007)

  47. [47]

    M. E. Fisher, Renormalization Group in Critical Phenomena and Quantum Field Theory, J. Gunton and M. S. Green, Eds., (Temple University, 1975)

  48. [48]

    D. R. Nelson, Coexistence-curve singularities in isotropic ferromagnets, Phys. Rev. B 13, 2222 (1976)

  49. [49]

    D. J. Amit and L. Peliti, On dangerously irrelevant operators, Ann. Phys. (N.Y .) 140, 207 (1982)

  50. [50]

    H. Shao, W. Guo, A. W. Sandvik, Monte Carlo renormalization flows in the space of relevant and irrelevant operators: Application to three-dimensional clock models, Phys. Rev. Lett. 124, 080602 (2020)

  51. [51]

    Blume, L

    M. Blume, L. M. Corliss, J. M. Hastings, and E. Schiller, Observation of an antiferromagnet in an induced staggered magnetic field: Dysprosium aluminum garnet near the tricritical point, Phys. Rev. Lett. 32, 544 (1974)

  52. [52]

    J. F. Dillon, E. Y . Chen, N. Giordano, and W. P. Wolf, Time-reversed antiferromagnetic states in dysprosium aluminum garnet, Phys. Rev. Lett. 33, 98 (1974)

  53. [53]

    Giordano and W

    N. Giordano and W. P. Wolf, Induced staggered magnetic fields in antiferromagnets: Microscopic mechanisms, Phys. Rev. B 21, 2008 (1980)

  54. [54]

    T. C. Fujita, Y . Kozuka, M. Uchida, A. Tsukazaki, T. Arima, and M. Kawasaki, Odd-parity magnetoresistance in pyrochlore iridate thin films with broken time-reversal symmetry, Scientific Reports 5, 10.1038/srep09711 (2015)

  55. [55]

    Hsieh, Jun J

    Tian Liang, Timothy H. Hsieh, Jun J. Ishikawa, Satoru Nakatsuji, Liang Fu and N. P. Ong, Orthogonal magnetization and symmetry breaking in pyrochlore iridate Eu2Ir2O7, Nat. Phys. 13, 599–603 (2017)

  56. [56]

    Yilin Wang, Hongming Weng, Liang Fu and Xi Dai, Noncollinear Magnetic Structure and Multipolar Order in Eu2Ir2O7, Phys. Rev. Lett. 119, 187203 (2017)

  57. [57]

    Lhotel et al

    E. Lhotel et al. Evidence for dynamic kagome ice. Nat. Commun. 9, 3786 (2018)

  58. [58]

    Xu et al

    J. Xu et al. Anisotropic exchange Hamiltonian, magnetic phase diagram, and domain inversion of Nd2Zr2O7. Phys. Rev. B 99, 144420 (2019) Acknowledgements The authors would like to thank Vaclav Petříček, Oleg Tchernyshyov, Yuan Wan , Yoshi Tokiwa, Jianhui Xu, Zheng Deng, Xiancheng Wang, Jie Shen, Huifen Ren and Shaokui Su for helpful discussions and experim...