Field-induced spin level crossings within a quasi-XY antiferromagnetic state in Ba2FeSi2O7
Minseong Lee,1, ∗Rico Schoenemann,1Hao Zhang,2, 3 David Dahlbom,3Tae-Hwan Jang,4
Seung-Hwan Do,5Andrew D. Christianson,5Sang-Wook Cheong,4, 6 Jae-Hoon Park,4, 7 Eric
Brosha,8Marcelo Jaime,1Kipton Barros,9Cristian D. Batista,3, †and Vivien S. Zapf1, ‡
1National High Magnetic Field Laboratory, Los Alamos National Laboratory, Los Alamos, NM 87545, USA.
2Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, TN, USA.
3Department of Physics and Astronomy, University of Tennessee, Knoxville, TN, USA.
4MPPHC-CPM, Max Planck POSTECH/Korea Research Initiative, Pohang, Republic of Korea.
5Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA.
6Rutgers Center for Emergent Materials and Department of Physics and Astronomy, Rutgers University, Piscataway, NJ, USA.
7Department of Physics, Pohang University of Science and Technology, Pohang, Republic of Korea.
8Materials Synthesis and Integration, Los Alamos National Laboratory, Los Alamos, NM 87545, USA.
9Theory Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA.
(Dated: March 15, 2023)
We present a high-field study of the strongly anisotropic easy-plane square lattice S= 2 quantum mag-
net Ba2FeSi2O7. This compound is a rare high-spin antiferromagnetic system with very strong easy-plane
anisotropy, such that the interplay between spin level crossings and antiferromagnetic order can be studied. We
observe a magnetic field-induced spin level crossing occurring within an ordered state. This spin level crossing
appears to preserve the magnetic symmetry while producing a non-monotonic dependence the order parame-
ter magnitude. The resulting temperature-magnetic field phase diagram exhibits two dome-shaped regions of
magnetic order overlapping around 30 T. The ground state of the lower-field dome is predominantly a linear
combination of |Sz= 0iand |Sz= 1istates, while the ground state of the higher-field dome can be ap-
proximated by a linear combination of |Sz= 1iand |Sz= 2istates. At 30 T, where the spin levels cross,
the magnetization exhibits a slanted plateau, the magnetocaloric effect shows a broad hump, and the electric
polarization shows a weak slope change. We determined the detailed magnetic phase boundaries and the spin
level crossings using measurements of magnetization, electric polarization, and the magnetocaloric effect in
pulsed magnetic fields to 60 T. We calculate these properties using a mean field theory based on direct products
of SU(5) coherent states and find good agreement. Finally, we measure and calculate the magnetically-induced
electric polarization that reflects magnetic ordering and spin level crossings. This multiferroic behavior provides
another avenue for detecting phase boundaries and symmetry changes.
I. INTRODUCTION
Antiferromagnets with strong easy-plane magnetic
anisotropy are a rich source of quantum magnetism and
quantum phase transitions1–3. In the limit of strong easy-
plane single-ion anisotropy, the Hamiltonian contains a
term HD=PiD(Sz
i)2with D > 0that is dominant
over the other interactions. This term splits the 2S+ 1
spin levels into an Sz= 0 ground state and excited states
with Sz=±1,±2,· · · ,±S, separated by energy gaps
of 12D,22D,· · · ,S2Dfor integer spin. The zdirection
coincides with the uniaixal anisotropy direction created
by HD. Since the Zeeman term generated by an external
magnetic field parallel to the z-direction commutes with HD,
the eigenstates of the single-ion Hamiltonian are preserved
and the energy gaps evolve linearly with the field. The energy
gaps between positive Szstates and the ground state are
suppressed by the field, resulting in successive spin level
crossings. Small exchange interactions between neighboring
ions J < D disperse and broaden these levels, producing
a mean field ground state that is direct product of linear
combinations of different |Szistates. The expectation value
of the in-plane component of the local magnetization is finite
in this wave function, leading to quasi-XY antiferromagnetic
long-range ordering occurring in dome-shaped regions of the
temperature-magnetic (T−H) phase diagram4–6. The order
parameter (planar staggered magnetization) can be described
as a complex number, whose amplitude is the magnitude of
the planar magnetization, while the phase determines the
direction of the spins. The pioneering work of Matsubara
and Matsuda showed that interacting spin-1/2systems can
be mapped into a interacting bosonic gas with hard-core
repulsion7. This mapping can be generalized for higher
spin systems8,9, implying that field-induced quantum phase
transitions of magnets with uniaxial symmetry belong to
the Bose-Einstein condensation (BEC) universality class.
In the presence of disorder or competing interactions, Bose
glass10–13 or even super-solidity14–17 can also occur. It is
important to note, however, that the Hamiltonians of real
spin systems always include terms that break the uniaxial
symmetry, but these terms are small for many compounds and
the BEC description becomes an excellent approximation,
particularly of the quantum phase transitions and their
excitations.2
Most work in this field has been done on S= 1/2
dimers6,18 or monoatomic S= 1 systems2,19. Higher spin sys-
tems (S > 1) with strong easy-plane anisotropy DJare
rare. However, when they occur, the energy of more than one
high-spin states can be lowered by magnetic field to become
the ground state in high fields, creating multiple domes of an-
tiferromagnetism in T−Hspace. For example, Ba3Mn2O8
forms exchange-coupled dimers of Mn5+(S= 1) spins5. In
applied magnetic fields, two dome-shaped regions of antifer-
arXiv:2210.14323v2 [cond-mat.str-el] 14 Mar 2023