Ordered Ground State Research Articles

Scaling Behavior in the Spectral and Power Density Dependent Photovoltaic Response of Hot Polaronic Heterojunctions

AbstractStrongly correlated manganites can be considered as model systems for the study of photovoltaic harvesting of hot polarons that can be excited from the electronically ordered ground state. In order to gain basic understanding of hot polaron harvesting, the deviations of the photovoltaic response of a heterojunction with polaronic absorber from a conventional semiconductor are analyzed. Specifically, the spectral and photon power density dependence of the open circuit voltage Uoc and the short circuit current density Jsc in heterojunctions consisting of orbital ordered Pr1‐xCaxMnO3 (x = 0.1, PCMO) thin films epitaxially grown on single crystalline (100) SrTiO3 (STO) and Nb‐doped (100) SrTiO3 (STNO) substrates are investigated. The observed behavior is fundamentally different from conventional solar cells, in which Uoc is limited by fast carrier relaxation to the band edges. Whereas the spectral and photon power dependence of Uoc of conventional semiconductor junctions is well described by the Shockley‐Queisser (SQ) theory, the hot polaron junctions surprisingly show a scaling law behavior of Uoc. Such scaling laws otherwise apply to equilibrium order parameters in second‐order phase transitions. It is concluded that its physical origin is the unique dependence of the quasi‐Fermi level splitting on temperature, photon energy, and power density in a hot polaron system.

Quantum Spin Dynamics Due to Strong Kitaev Interactions in the Triangular-Lattice Antiferromagnet CsCeSe_{2}.

The extraordinary properties of the Kitaev model have motivated an intense search for new physics in materials that combine geometrical and bond frustration. In this Letter, we employ inelastic neutron scattering, spin wave theory, and exact diagonalization to study the spin dynamics in the perfect triangular-lattice antiferromagnet (TLAF) CsCeSe_{2}. This material orders into a stripe phase, which is demonstrated to arise as a consequence of the off-diagonal bond-dependent terms in the spin Hamiltonian. By studying the spin dynamics at intermediate fields, we identify an interaction between the single-magnon state and the two-magnon continuum that causes decay of coherent magnon excitations, level repulsion, and transfer of spectral weight to the continuum that are controlled by the strength of the magnetic field. Our results provide a microscopic mechanism for the stabilization of the stripe phase in TLAF and show how complex many-body physics can be present in the spin dynamics in a magnet with strong Kitaev coupling even in an ordered ground state.

Band Structure Calculations, Magnetic Properties and Magnetocaloric Effect of GdCo1.8M0.2 Compounds with M = Fe, Mn, Cu, Al

The magnetic properties, band structure results, and magnetocaloric effect of GdCo1.8M0.2 with M = Fe, Mn, Cu, and Al are reported. The band structure calculations demonstrate that all the samples have a ferrimagnetically ordered ground state, in perfect agreement with the magnetic measurements. Calculated magnetic moments and variation with the alloy composition are strongly influenced by hybridisation mechanisms as sustained by an analysis of the orbital projected local density of states. The XPS measurements reveal no significant shift in the binding energy of the investigated Co core levels with a change in the dopant element. The Co 3s core-level spectra gave us direct evidence of the local magnetic moments on Co sites and an average magnetic moment of 1.3 µB /atom was found, being in good agreement with the theoretical estimation and magnetic measurements. From the Mn 3s core-level spectra, a value of 2.1 µB/Mn was obtained. The symmetric shapes of magnetic entropy changes, the Arrott plots, and the temperature dependence of Landau coefficients clearly indicate a second-order phase transition. The relative cooling power, RCP(S), normalized relative cooling power, RCP(∆S)/∆B, and temperature-averaged entropy change values indicate that these compounds could be promising candidates for applications in magnetic refrigeration devices.

Effect of the dipolar interaction on the dynamic hysteresis of properties of 2D-nanodisks: In-plane driving field case

The effect of the dipolar interaction strength on the dynamical hysteresis behavior of the nanodisk is investigated within the Heisenberg model and Monte Carlo simulation. The system is under the effect of the dynamical periodic sinusoidal magnetic field. It has been shown that dipolar interaction strength can induce dynamical phase transition and the transition from the fully dynamically ordered ground state to the weak dynamically ordered ground state. Besides, it has been shown that rising dipolar interaction raises the hysteresis loop area for higher temperatures and a linear relationship between the peak positions of the hysteresis loop area and the temperature locations of the peaks obtained.

Flat-band hybridization between f and d states near the Fermi energy of SmCoIn5

We present high-quality angle-resolved photoemission (ARPES) and density functional theory calculations (DFT+U) of SmCoIn5. We find broad agreement with previously published studies of LaCoIn5 and CeCoIn51,2, confirming that the Sm 4f electrons are mostly localized. Nevertheless, our model is consistent with an additional delocalized Sm component, stemming from hybridization between the 4f electrons and the metallic bands at “hot spot” positions in the Brillouin zone. The dominant hot spot, called γZ, is similar to a source of delocalized f states found in previous experimental and theoretical studies of CeCoIn51,3. In this work, we identify and focus on the role of the Co d states in exploring the relationship between heavy quasiparticles and the magnetic interactions in SmCoIn5, which lead to a magnetically ordered ground state from within an intermediate valence scenario4–6. Specifically, we find a globally flat band consisting of Co d states near E = − 0.7 eV, indicating the possibility of enhanced electronic and magnetic interactions in the “115” family of materials through localization in the Co layer, and we discuss a possible origin in geometric frustration. We also show that the delocalized Sm 4f states can hybridize directly with the Co 3dxz/3dyz orbitals, which occurs in our model at the Brillouin zone boundary point R in a band that is locally flat and touches the Fermi level from above. Our work identifies microscopic ingredients for additional magnetic interactions in the “115” materials beyond the RKKY mechanism, and strongly suggests that the Co d bands are an important ingredient in the formation of both magnetic and superconducting ground states.

Ground-state phase diagram and universality of sequential topological valence-bond-solid quantum transitions in a mixed tetramer chain

The ground-state ordering of a quantum mixed-spin Heisenberg tetramer chain composed of an alternate sequence of s = 1 and S=3/2 dimers is studied in detail as a function of two considered exchange interactions ascribed to similar and dissimilar spin pairs. At zero magnetic field, the ferrimagnetic mixed spin-(1, 1, 3/2, 3/2) Heisenberg tetramer chain displays, depending on a mutual interplay between two considered exchange interactions, three distinct gapped valence-bond-solid phases separated by gap-closing quantum critical points. Using density-matrix renormalization group calculations we construct the full ground-state phase diagram as a function of the interaction ratio and magnetic field, which exhibits besides three gapped valence-bond-solid phases special Kosterlitz–Thouless and topological quantum critical points. A tangential finite-size scaling analysis is employed to obtain precise estimates of the zero-field valence-bond-solid transitions and unveil their common logarithmic correction to a power-law scaling of the correlation length.

Resilience of the Aurivillius structure upon La and Cr doping in a Bi5Ti3FeO15 multiferroic.

Combining the experimental techniques of high-resolution X-ray diffraction, magnetometry, specific heat measurement, and X-ray photoelectron, Raman and dielectric spectroscopy techniques, we have studied the influence of La and Cr doping on the crystal structure and magnetism of the room temperature Aurivillius multiferroic Bi5Ti3FeO15 by investigating the physical properties of (Bi4La)Ti3FeO15 and Bi5Ti3 (Fe0.5Cr0.5)O15. The parent (Bi5Ti3FeO15) and the doped ((Bi4La)Ti3FeO15 and Bi5Ti3(Fe0.5Cr0.5)O15) compounds crystallize in the A21am space group, which is confirmed through our analysis of high-resolution synchrotron X-ray diffraction data obtained on phase-pure polycrystalline powders. We determined the oxidation states of the metal atoms in the studied compounds as Fe3+, Ti4+, Cr3+, and La3+ through the analysis of X-ray photoelectron spectroscopy data. The magnetic susceptibilities of the three compounds are marked by the absence of a long-range ordered ground state, but dominated by superparamagnetic clusters with dominant antiferromagnetic interactions. This signature of short-range magnetism is also seen in specific heat as a low temperature enhancement which is suppressed upon the application of external magnetic fields up to 8 T. Our dielectric spectroscopy experiments showed that the three studied compounds display similar features in the dielectric constant measured as a function of frequency. However, upon doping La at the Bi site, the width of the ferroelectric hysteresis loop increases for (Bi4La)Ti3FeO15 compared to that of the parent compound Bi5Ti3FeO15, and with Cr doping, Bi5Ti3(Fe0.5Cr0.5)O15 becomes a leaky dielectric. The resilience of the Aurivillius crystal structure towards doping of La at the Bi site and Cr at the Fe site is clearly seen in the bulk properties of magnetic susceptibility, specific heat and the average crystal structure. The relevance of changes in the local structure is evident from our Raman spectroscopy and X-ray pair distribution function studies.

Ordered ground state configurations of the asymmetric Wigner bilayer system-Revisited with unsupervised learning.

We have reanalyzed the rich plethora of ground state configurations of the asymmetric Wigner bilayer system that we had recently published in a related diagram of states [Antlanger et al., Phys. Rev. Lett. 117, 118002 (2016)], comprising roughly 60 000 state points in the phase space spanned by the distance between the plates and the charge asymmetry parameter of the system. In contrast to this preceding contribution where the classification of the emerging structures was carried out "by hand," we have used for the present contribution machine learning concepts, notably based on a principal component analysis and a k-means clustering approach: using a 30-dimensional feature vector for each emerging structure (containing relevant information, such as the composition of the configuration as well as the most relevant order parameters), we were able to reanalyze these ground state configurations in a considerably more systematic and comprehensive manner than we could possibly do in the previously published classification scheme. Indeed, we were now able to identify new structures in previously unclassified regions of the parameter space and could considerably refine the previous classification scheme, thereby identifying a rich wealth of new emerging ground state configurations. Thorough consistency checks confirm the validity of the newly defined diagram of states.

Robust Spin Liquidity in 2D Metal-Organic Framework Cu3 (HHTP)2 with S=1 /2 Kagome Lattice.

On one hand electron or hole doping of quantum spin liquid (QSL) may unlock high-temperature superconductivity and on the other hand it can disrupt the spin liquidity, giving rise to a magnetically ordered ground state. Recently, a 2D MOF, Cu3 (HHTP)2 (HHTP - 2,3,6,7,10,11-hexahydroxytriphenylene), containing Cu(II) S= frustrated spins in the Kagome lattice is emerging as a promising QSL candidate. Herein, we present an elegant in situ redox-chemistry strategy of anchoring Cu3 (HHTP)2 crystallites onto diamagnetic reduced graphene oxide (rGO) sheets, resulting in the formation of electron-doped Cu3 (HHTP)2 -rGO composite which exhibited a characteristic semiconducting behavior (5 K to 300 K) with high electrical conductivity of 70 S ⋅ m-1 and a carrier density of ~1.1×1018 cm-3 at 300 K. Remarkably, no magnetic transition in the Cu3 (HHTP)2 -rGO composite was observed down to 1.5 K endorsing the robust spin liquidity of the 2D MOF Cu3 (HHTP)2 . Specific heat capacity measurements led to the estimation of the residual entropy values of 28 % and 34 % of the theoretically expected value for the pristine Cu3 (HHTP)2 and Cu3 (HHTP)2 -rGO composite, establishing the presence of strong quantum fluctuations down to 1.5 K (two times smaller than the value of the exchange interaction J).

Weyl Fermions and broken symmetry phases of laterally confined 3He films

Broken symmetries in topological condensed matter systems have implications for the spectrum of Fermionic excitations confined on surfaces or topological defects. The Fermionic spectrum of confined (quasi-2D) 3He-A consists of branches of chiral edge states. The negative energy states are related to the ground-state angular momentum, , for Cooper pairs. The power law suppression of the angular momentum, for , in the fully gapped 2D chiral A-phase reflects the thermal excitation of the chiral edge Fermions. We discuss the effects of wave function overlap, and hybridization between edge states confined near opposing edge boundaries on the edge currents, ground-state angular momentum and ground-state order parameter of superfluid 3He thin films. Under strong lateral confinement, the chiral A phase undergoes a sequence of phase transitions, first to a pair density wave (PDW) phase with broken translational symmetry at . The PDW phase is described by a periodic array of chiral domains with alternating chirality, separated by domain walls. The period of PDW phase diverges as the confinement length . The PDW phase breaks time-reversal symmetry, translation invariance, but is invariant under the combination of time-reversal and translation by a one-half period of the PDW. The mass current distribution of the PDW phase reflects this combined symmetry, and originates from the spectra of edge Fermions and the chiral branches bound to the domain walls. Under sufficiently strong confinement a second-order transition occurs to the non-chiral ‘polar phase’ at , in which a single p-wave orbital state of Cooper pairs is aligned along the channel.

Effect of neighbouring molecules on ground-state properties of many-body polar linear rotor systems

A path integral ground state approach has been used to estimate the ground-state energy and structural properties of hydrogen fluoride molecules pinned to a one-dimensional lattice. In the simulations, the molecules are assumed to be rigid, and only the continuous rotational degrees of freedom are considered. The constituents of a many-body system interact through the dipole-dipole interaction because the molecules have a permanent dipole moment. The workability of our approach has been demonstrated by estimating the ground-state energy, order parameter and nearest neighbour correlation for the systems of 2 and 3 HF molecules using quantum Monte Carlo simulations based on the path integral ground state methodology. The results agree satisfactorily with those obtained from exact Hamiltonian matrix diagonalisation. In addition, the effect of neighbours on the ground state properties has been investigated for larger systems. The converged ground-state energy per neighbour as a function of inter-nuclear separation is considered to be the equation of state for the system. The chemical potential of the rotor chains also supports the convergence test.

Kondo Lattice Model of Magic-Angle Twisted-Bilayer Graphene: Hund's Rule, Local-Moment Fluctuations, and Low-Energy Effective Theory.

We apply a generalized Schrieffer-Wolff transformation to the extended Anderson-like topological heavy fermion (THF) model for the magic-angle (θ=1.05°) twisted bilayer graphene (MATBLG) [Phys. Rev. Lett. 129, 047601 (2022)PRLTAO0031-900710.1103/PhysRevLett.129.047601], to obtain its Kondo lattice limit. In this limit localized f electrons on a triangular lattice interact with topological conduction c electrons. By solving the exact limit of the THF model, we show that the integer fillings ν=0,±1,±2 are controlled by the heavy f electrons, while ν=±3 is at the border of a phase transition between two f-electron fillings. For ν=0,±1,±2, we then calculate the Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions between the f moments in the full model and analytically prove the SU(4) Hund's rule for the ground state which maintains that two f electrons fill the same valley-spin flavor. Our (ferromagnetic interactions in the) spin model dramatically differ from the usual Heisenberg antiferromagnetic interactions expected at strong coupling. We show the ground state in some limits can be found exactly by employing a positive semidefinite "bond-operators" method. We then compute the excitation spectrum of the f moments in the ordered ground state, prove the stability of the ground state favored by RKKY interactions, and discuss the properties of the Goldstone modes, the (reason for the accidental) degeneracy of (some of) the excitation modes, and the physics of their phase stiffness. We develop a low-energy effective theory for the f moments and obtain analytic expressions for the dispersion of the collective modes. We discuss the relevance of our results to the spin-entropy experiments in TBG.

High Proton Conductivity of Magnetically Ordered 2D Oxalate-Bridged [MnIICrIII] Coordination Polymers with Irregular Topology.

Two heterometallic coordination polymers {[NH(CH3)2(C2H5)]8[Mn4Cl4Cr4(C2O4)12]}n (1) and {[NH(CH3)-(C2H5)2]8[Mn4Cl4Cr4(C2O4)12]}n (2) were obtained by slow evaporation of an aqueous solution containing the building block [A]3[Cr(C2O4)3] [A = (CH3)2(C2H5)NH+ or (CH3)(C2H5)2NH+] and MnCl2·2H2O. The isostructural compounds comprise irregular two-dimensional (2D) oxalate-bridged anionic layers [Mn4Cl4Cr4(C2O4)12]n8n- with a Shubnikov plane net fes topology designated as (4·82), interleaved by the hydrogen-bonded templating cations (CH3)2(C2H5)NH+ (1) or (CH3)(C2H5)2NH+ (2). They exhibit remarkable humidity-sensing properties and very high proton conductivity at room temperature [1.60 × 10-3 (Ω·cm)-1 at 90% relative humidity (RH) of 1 and 9.6 × 10-4 (Ω·cm)-1 at 94% RH of 2]. The layered structure facilitates the uptake of water molecules, which contributes to the enhancement of proton conductivity at high RH. The better proton transport observed in 1 compared to that in 2 can be tentatively attributed to the higher hydrophilicity of the cations (CH3)2(C2H5)NH+, which is closely related to their affinity for water molecules. The original topology of the anionic networks in both compounds leads to the development of interesting magnetic phases upon cooling. The magnetically ordered ground state can be described as the coupling of ferromagnetic spin chains in which Mn2+ and Cr3+ ions are bridged by bis(bidentate) oxalate groups into antiferromagnetic planes through monodentate-bidentate oxalate bridges in the layers, which are triggered to long-range order below temperature 4.45 K via weaker interlayer interactions.

High-pressure order-disorder transition in Mg2SiO4 : Implications for super-Earth mineralogy

(Mg, Fe)SiO$_3$ post-perovskite is the highest pressure silicate mineral phase in the Earth's interior. The extreme pressure and temperature conditions inside large extrasolar planets will likely lead to phase transitions beyond pPv. In this work we have explored the high-pressure phase relations in Mg$_2$SiO$_4$ using computations based on density functional theory. We find that a partially disordered I-42d type structure would be stable in the interiors of these super-Earth planets. The discovery of a structure where two very dissimilar cations, Mg$^{2+}$ and Si$^{4+}$ occupy the same crystallographic site opens up a domain of interesting crystal chemistry and provides a foundation for other silicates and oxides with mixed occupancy. We have explored the mechanism of the phase transition from the ordered ground state and the effect of the disordering on electronic properties of the silicate phase.

Tunable magnetism and large inverse magnetocaloric effect in Shastry-Sutherland compounds Tb2Be2Si1-xGexO7 (0 ≤ x ≤ 1)

The structure, magnetism and magnetocaloric effect (MCE) have been systematically investigated on the series of Tb2Be2Si1-xGexO7 (0 ≤ x ≤ 1) compounds, where the magnetic Tb3+ ions are located on the quasi-two-dimensional Shastry−Sutherland-Lattice (SSL). The X-ray diffraction measurements reveal that all compounds crystallize into the tetragonal structure with space group P4¯21m, the lattice parameters and nearest neighbor Tb-Tb distances increase monotonically with the variation of doping content x. All samples show the antiferromagnetic (AFM) ordered ground state and field (μ0H) induced spin-flop (SP) transitions below Néel temperature (TN). Moreover, the constructed field-temperature (μ0H−T) phase diagrams of Tb2Be2Si1-xGexO7 (0 ≤ x ≤ 1) reveal that the AFM phase shifts to the low temperature and low field sides with increased x. Below TN, we can observe a crossover behavior of MCE from the low-field (Δμ0H<2.58 T) inverse MCE (IMCE) to the CMCE at high fields (Δμ0H>2.58 T). Among these serial compounds, Tb2Be2SiO7 exhibits the maximal IMCE with a value of 8.30 J kg−1 K−1 under Δμ0H= 0.7 T and 2 K, comparable to the MCE (8.6 J kg−1 K−1) of the benchmark material Tb3Ga5O12 at Δμ0H = 2 T. Additionally, it has a slightly better MCE performance than the materials RE3Ga5O12 (RE = Dy, Tb, Gd) under Δμ0H = 2 T. The realization of large IMCE performance for Tb2Be2SiO7 provides an alternative route for magnetic refrigeration at low field and liquid Helium temperatures.

Order-disorder transition in the zero-temperature Ising model on random graphs.

The zero-temperature Ising model is known to reach a fully ordered ground state in sufficiently dense random graphs. In sparse random graphs, the dynamics gets absorbed in disordered local minima at magnetization close to zero. Here, we find that the nonequilibrium transition between the ordered and the disordered regime occurs at an average degree that slowly grows with the graph size. The system shows bistability: The distribution of the absolute magnetization in the reached absorbing state is bimodal, with peaks only at zero and unity. For a fixed system size, the average time to absorption behaves nonmonotonically as a function of average degree. The peak value of the average absorption time grows as a power law of the system size. These findings have relevance for community detection, opinion dynamics, and games on networks.

Charge order in the kagome lattice Holstein model: a hybrid Monte Carlo study

The Holstein model is a paradigmatic description of the electron-phonon interaction, in which electrons couple to local dispersionless phonon modes, independent of momentum. The model has been shown to host a variety of ordered ground states such as charge density wave (CDW) order and superconductivity on several geometries, including the square, honeycomb, and Lieb lattices. In this work, we study CDW formation in the Holstein model on the kagome lattice, using a recently developed hybrid Monte Carlo simulation method. We present evidence for sqrt{3}times sqrt{3} CDW order at an average electron filling of 〈n〉 = 2/3 per site, with an ordering wavevector at the K-points of the Brillouin zone. We estimate a phase transition occurring at Tc ≈ t/18, where t is the nearest-neighbor hopping parameter. Our simulations find no signature of CDW order at other electron fillings or ordering momenta for temperatures T ≥ t/20.

Charged magnons on the surface of a topological insulator

We study a system of two-dimensional Dirac electrons (as is realized on the surface of a 3D topological insulator) coupled to an array of localized spins. The spins are coupled ferromagnetically to each other, forming an ordered ground state with low-energy spin-wave excitations (magnons). The Dirac electrons couple to the spins through a spin-dependent effective Zeeman field. The out-of-plane effective Zeeman field therefore serves as a Dirac mass that gaps the electronic spectrum. Once a spin is flipped, it creates a surrounding domain in which the sign of the Dirac mass is opposite to that of the rest of the sample. Therefore, an electronic bound state appears on the domain wall, as predicted by Jackiw and Rebbi. However, in a quantum magnet, a localized spin flip does not produce an eigenstate. Instead, the eigenstates correspond to delocalized spin waves (magnons). As in the case of the single flipped spin, the delocalized magnon also binds an in-gap electronic state. We name this excitation a `Jackiw-Rebbi-Magnon' (JRM) and study its signature in the dynamic spin susceptibility. When the sample is tunnel-coupled to an electronic reservoir, a magnon produced in a system without any electrons hybridizes with a JRM (which binds a single electron), producing magnon-JRM polaritons. For such a system, we identify a quantum phase transition when the magnon-JRM polariton energy falls below that of the fully polarized ferromagnetic ground state.

Theory of the collective behavior of two-dimensional periodic ensembles of dipole-coupled magnetic nanoparticles

The magnetic behavior of two-dimensional (2D) periodic ensembles of dipole-coupled nanoparticles (NPs) is investigated theoretically by considering all possible structural organizations in 2D Bravais lattices, namely, square, triangular, rectangular, rhombic, and oblique lattices. The different interaction-energy landscapes (ELs) are characterized by determining the local minima that define the stable and metastable magnetic configurations as well as the transition states connecting them. The topology of the resulting ergodic ensemble of stationary states is analyzed from both local and energy perspectives by calculating the corresponding kinetic networks and disconnectivity graphs. For all lattices, the magnetic orders of the ground-state and low-lying metastable configurations are identified, including the elementary relaxation processes between them. A remarkable profound dependence of the collective magnetic behavior on the structural arrangement of the NPs is revealed. Square and triangular nanostructures are extremely good structure seekers, showing starlike kinetic networks and disconnectivity graphs with palm-tree form. Their continuously degenerate ground states are throughout-reaching hubs connected to nearly all excited magnetic configurations over a single first-order saddle point. Rhombic ensembles have double-funnel ELs in which the time-inversion-related ground states play the role of hubs with exceedingly high connectivity densities. These systems need to undergo very few elementary transitions with small downward energy barriers to relax from any configuration towards one of the ground states. However, the energy barriers between the two funnels are quite large and, furthermore, they increase as the number of particles in the unit cell is increased. Therefore, ergodicity breaking is expected in the thermodynamic limit. Finally, the rectangular and oblique lattices show contrasting ground-state orders, the former consists of an antiferromagnetic alternation of head-to-tail chains of spins while the latter is ferromagnetic. Nevertheless, the ELs of these nanostructures are qualitatively very similar. In both cases, all excited magnetic configurations consist of independent flipping of the chains of spins that are formed along the direction of the shorter Bravais vector. The whole energy spectrum of metastable magnetic configurations can thus be mapped to a one-dimensional Ising model. As a result, the kinetic network of stationary points is latticelike and the disconnectivity graphs have a willow-tree form. In conclusion, the collective magnetic behaviors of nanostructures having different point-group symmetries are contrasted and related by varying the parameters of the lattices that interpolate between them.

Double dome structure of the Bose–Einstein condensation in diluted S = 3/2 quantum magnets

Bose–Einstein condensation (BEC) in quantum magnets, where bosonic spin excitations condense into ordered ground states, is a realization of BEC in a thermodynamic limit. Although previous magnetic BEC studies have focused on magnets with small spins of S ≤ 1, larger spin systems potentially possess richer physics because of the multiple excitations on a single site level. Here, we show the evolution of the magnetic phase diagram of S = 3/2 quantum magnet Ba2CoGe2O7 when the averaged interaction J is controlled by a dilution of magnetic sites. By partial substitution of Co with nonmagnetic Zn, the magnetic order dome transforms into a double dome structure, which can be explained by three kinds of magnetic BECs with distinct excitations. Furthermore, we show the importance of the randomness effects induced by the quenched disorder: we discuss the relevance of geometrical percolation and Bose/Mott glass physics near the BEC quantum critical point.