## Altermagnetism on the Shastry-Sutherland lattice

Francesco Ferrari and Roser Valent´ ı Institute for Theoretical Physics, Goethe University Frankfurt, Max-von-Laue-Straße 1, D-60438 Frankfurt a.M., Germany

∗ †

(Dated: November 21, 2024)

Motivated by recent developments on altermagnetism, we investigate the Hubbard model on the Shastry-Sutherland lattice, where the onsite repulsion between electrons induces the onset of a staggered magnetic order in which opposite magnetic sublattices are related to each other only by rotations and glide reflections, and not by inversion nor translations. As a consequence, the magnetic phase displays an altermagnetic character, i.e. a finite (non-relativistic) Zeeman splitting of the electronic excitations, despite the absence of a net magnetization. By means of a variational Monte Carlo approach based on Jastrow-Slater wave functions, which accounts for electronic correlations beyond the single-particle approximation, we study how d -wave altermagnetism shows up in the ground state properties and in the spectral function of the system. We characterize the metalinsulator transition between altermagnetic phases at half-filling and the stability of altermagnetism as a function of doping. The calculation of the single-particle spectral function displays the spinsplit nature of the electronic excitations, also within the Mott insulating regime. We discuss possible realizations of the model in the context of organic α -or κ -(BEDT-TTF) X charge transfer salts. 2

## I. INTRODUCTION

Altermagnetism has been recently framed as a new form of magnetism [1, 2] which shows no net magnetization, like in antiferromagnets, but fosters a (nonrelativistic) Zeeman splitting of the electronic bands, analogously to ferromagnets [3-6]. This phenomenological distinction is a consequence of the symmetries of the magnetic order and is formalized in the framework of spin groups [2]. Concretely, the key aspect to discriminate between antiferromagnetic and altermagnetic (collinear) orders is identifying the set of symmetries which transform sublattices of opposite spin polarizations into each other. In usual antiferromagnets, this is achieved by inversion or translations. Therefore, the combination of such lattice symmetries and time reversal leaves the magnetic order invariant and implies a spin-degeneracy of the electronic bands. This is not the case of altermagnets, where, instead, the effect of time reversal is compensated only by the application of lattice rotations and/or reflections. As a consequence, the altermagnetic character manifests itself in a momentum-dependent spin splitting of the band structure, with an even-parity order parameter with d , g or i -wave symmetry, which is compatible with the presence of a vanishing net magnetization [1, 2, 6-9].

Altermagnetism is presently a topic of intensive scrutiny, not only for the possible application of nonrelativistic-induced spin splittings for spintronics purposes [10] but, on a more fundamental level, for its potential towards unveiling unexplored phases [11-13]. In this respect, a number of studies have been devoted to the investigation of interacting models where altermagnetism is either imprinted in the local crystallographic symmetries of the model and triggered by the electronic interac-

∗ ferrari@itp.uni-frankfurt.de

† valenti@itp.uni-frankfurt.de

0

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00841v2-figure-6e049153-2c65-4fb2-84c2-eb97e91ec5ac.jpg)

1

Figure 1. Illustration of the Shastry-Sutherland lattice. The yellow box denotes the four sites forming the unit cell ( η = 0 1 2 3 sublattices) and the orange arrows represent the , , , Bravais unit vectors a 1 = (1 0) and , a 2 = (0 1). , Full/empty dots indicate sites of even/odd sublattices. The hoppings of Eq. (1), t and t ′ , are located on the solid and dashed bonds, respectively. The red/blue arrows depict the Q = 0 staggered magnetic order. The symmetry operations that connect sublattices of opposite spin-polarization, i.e. C 4 rotations and G ,G x y glide reflections, are marked by the mauve and green arrows.

tions [14-16], or, as recently proposed [17], spontaneously generated by orbital ordering. While the above studies suggest the realization of interesting altermagnetic phases, most calculations are performed at the singleparticle level without including correlation effects (beyond mean-field) that could modify or even prevent the formation of the conjectured order, and some of the proposed models might be somewhat speculative to have them achieved in real materials.

In this context, a showcase model in quantum magnetism that has received a large amount of attention in the last decades is the Shastry-Sutherland model [1821] which, in the strong coupling limit, consists of in-

teracting S = 1 2 spins on a checkerboard square lat-/ tice, where every other square contains an interaction along alternating diagonals (see Fig. 1); this results in a competition between nearest-neighbor and diagonal exchange interactions, which leads to the formation of distinct ground states depending on the degree of frustration, such as long-range magnetic ordered phases, dimer and plaquette orders, and possibly spin-liquid states [22-26]. One of the most studied material realizations of this model is SrCu 2 (BO ) 3 2 [27], whose behavior under applied magnetic field and pressure illustrates the complex phase diagram of the Shastry-Sutherland model [25, 28-32]. In view of the richness of phases that the Shastry-Sutherland model entails at the level of the Heisenberg Hamiltonian, a question that arises is which phases emerge away from the Mott limit of localized electrons [33], especially considering that the symmetries of the underlying lattice could potentially bear altermagnetic order. In this regard, a distorted version of the Shastry-Sutherland lattice is realized by the molecular arrangements of BEDT-TTF molecules in organic charge-transfer salts [34-36], such as κ -(BEDTTTF) Cu[N(CN) ]Cl, which displays a momentum de2 2 pendent spin-splitting of the band structure in its magnetically ordered phase [37-39].

Motivated by these observations, in this work we investigate the ground state and excitation spectra of the fermionic Hubbard model on the Shastry-Sutherland lattice (Sec. II A). We make use of a variational Monte Carlo approach that allows us to accurately treat the effects of electronic correlations (Secs. II B and II C), which have not been considered in previous mean-field studies of altermagnetic models. Starting from a symmetry analysis of the non-interacting band structure (Sec. III A), we investigate the onset of altermagnetic order at half-filling as a function of the Hubbard repulsion, uncovering both metallic and insulating altermagnetic phases (Sec. III B). We characterize the altermagnetic order in these two phases and discuss the impact of electronic correlations on the metal-insulator transition. We then assess the stability of altermagnetism as a function of electron and hole doping (Sec. III C). Finally, we compute the singleparticle spectral function, which shows signs of altermagnetism, also in the Mott insulating regime (Sec. III D).

## II. MODEL AND METHOD

## A. Hubbard Model

The focus of this work is on the fermionic Hubbard model [40, 41]

$$\mathcal { H } & = - \, t \sum _ { \langle i, j \rangle, \sigma } c ^ { \dagger } _ { i \sigma } c _ { j \sigma } - t ^ { \prime } \sum _ { ( i, j ) \in d, \sigma } c ^ { \dagger } _ { i \sigma } c _ { j \sigma } + h. c. \quad \text{ in w} \} \\ & \quad + U \sum _ { i } n _ { i \uparrow } n _ { i \downarrow }$$

on the Shastry-Sutherland lattice, whose periodicity can be represented by a square lattice with a four-site unit cell, as depicted in Fig. 1. The operator c † iσ ( c iσ ) creates (annihilates) an electron at site i with spin σ ∈ {↑ ↓} , and n iσ = c † iσ c iσ . The Hamiltonian contains two (spinisotropic) hopping terms, t at first neighbors, and t ′ residing on a subset d of the square-lattice diagonals. The Hubbard repulsion is denoted by U . The model of Eq. (1) is invariant under the p g 4 (two-dimensional) space group, which features C 4 rotations around the center of the unit cell and reflections with respect to the axes parallel to the diagonal bonds d . Additional glide reflections are present, e.g. with respect to the axes of first-neighbor bonds [25], as illustrated in Fig. 1. As an alternative to the above notation, it is sometimes convenient to denote the lattice sites as ( R,η ), i.e. specifying the Bravais vector R of the unit cell and the sublattice kind η = 0 1 2 3. , , , This is particularly useful for the definition of Fourier transforms, e.g. in the case of the annihilation operators c k,η,σ = 1 √ N c ∑ R e -ikR c R,η,σ , where N c denotes the number of unit cells.

In absence of diagonal hoppings ( t ′ = 0), the Hamiltonian of Eq. (1) reduces to the standard square-lattice Hubbard model, whose ground state at half-filling is insulating and antiferromagnetically ordered for any U &gt; 0. The full model with both hopping terms, instead, has been widely investigated in its strong coupling regime ( U ≫ t, t ′ ), where it reduces to the Shastry-Sutherland Heisenberg model [18]. Here, tuning the ratio of the resulting exchange couplings, the staggered magnetic order melts due to frustration and leaves room for nonmagnetic ground states [23-26], among which an exact product state of singlets on the diagonal bonds is realized [18, 20]. In this work, on the other hand, we explore the ground state phase diagram of the Hubbard Hamiltonian (1) for small and moderate values of U/t , investigating the metal-insulator transition at half-filling and the effects of doping, with a special focus on the altermagnetic properties of the magnetically ordered phase at finite t ′ . We limit our study to t /t ′ ≤ 1 (and t, t ′ &gt; 0), such that frustration effects are not strong enough to melt magnetic order in the strong coupling regime [23, 33].

## B. Variational Monte Carlo: ground state

We tackle the Hubbard Hamiltonian of Eq. (1) by means of a variational Monte Carlo approach built on Jastrow-Slater wave functions of the form

$$| \Psi _ { 0 } \rangle = \mathcal { P J } _ { s } \mathcal { J } _ { n } | \Phi \rangle \,,$$

in which long-range density

$$\mathcal { J } _ { \mathbf n } = \exp \left ( - \, \frac { 1 } { 2 } \sum _ { i, j } v _ { i, j } n _ { i } n _ { j } \right ) \, \quad \quad ( 3 )$$

and spin

$$\mathcal { J } _ { s } = \exp \left ( \, - \, \frac { 1 } { 2 } \sum _ { i, j } u _ { i, j } S _ { i } ^ { z } S _ { j } ^ { z } \right ) \, \quad \ \ ( 4 ) \, \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Jastrow factors are applied to a fermionic uncorrelated state | Φ , ⟩ i.e. a Slater determinant [42]. The variational Ansatz (2) provides a flexible and efficient tool to account for electronic correlations, going beyond singleparticle approximations (i.e., mean-field methods). In this regard, for instance, the density Jastrow factor J n can be viewed as a generalization of the Gutzwiller correlator [43] and, contrary to the latter, is capable of describing genuine Mott insulators with finite charge fluctuations [44, 45]. The uncorrelated part of the variational state, | Φ , can be chosen to encode fundamental ingredi-⟩ ents of the physical state, e.g. the presence of symmetrybreaking orders. We take | Φ to be the ground state of an ⟩ auxiliary quadratic Hamiltonian, which contains hopping terms and a staggered magnetic field (with pitch vector Q = 0)

$$\mathcal { H } _ { \text{aux} } = & - \sum _ { \langle i, j \rangle, \sigma } c _ { i \sigma } ^ { \dagger } c _ { j \sigma } - \chi ^ { \prime } \sum _ { ( i, j ) \in d, \sigma } c _ { i \sigma } ^ { \dagger } c _ { j \sigma } + h. c. \quad \text{$\text{$\text{$\text{$\text{eter $\xi$}$}$}$} \\ & - h \sum _ { i } ( - ) ^ { \eta _ { i } } ( c _ { i \uparrow } ^ { \dagger } c _ { i \uparrow } - c _ { i \downarrow } ^ { \dagger } c _ { i \downarrow } ) \,. \quad \text{$(5)} \text{$\text{appr}$} \\ & \text{$\text{ized} \,v$} }$$

̸

Here, η i denotes the sublattice kind of site i in the convention of Fig. 1. The diagonal hopping χ ′ and the Zeeman field h are variational parameters. At half-filling, if h = 0 the ground state of H aux is a paramagnetic metal (as long as χ ′ ≤ 2). On the other hand, when h is finite, | Φ is an altermagnetic metallic state if ⟩ χ ′ ≤ | h | ( χ ′ = 0) and an altermagnetic insulator if χ &gt; h ′ | | .

The calculations are carried out on a lattice of N c = 8 × 8 unit cells ( N = 256 sites), with periodicantiperiodic boundary conditions in the directions of a 1 and a 2 , respectively. The projector P included in the variational Ansatz [Eq. (2)] constrains the wave function to the subspace with total number of electrons N e and S z tot = ∑ i ( n i, ↑ -n i, ↓ ) / 2 = 0. Calculations are performed at half-filling, n = N /N e = 1, and in presence of finite doping δ = n -1. The variational parameters of the auxiliary Hamiltonian ( χ ′ and h ) and the Jastrow factors ( v i,j and u i,j [46]) are numerically optimized by the stochastic reconfiguration method [42, 47] to obtain the lowest variational energy

$$E _ { 0 } = \frac { \langle \Psi _ { 0 } | \mathcal { H } | \Psi _ { 0 } \rangle } { \langle \Psi _ { 0 } | \Psi _ { 0 } \rangle } \quad \quad \ ( 6 ) \quad \frac { \text{ac} } { \text{(es)} }$$

for the Hubbard model (1).

## C. Variational Monte Carlo: spectral functions

In addition to ground state properties, we employ the method introduced in Ref. [48] to compute the singleparticle spectral function of the Hubbard model (1). This approach relies on the definition of suitable variational excitations in the sectors with N e -1 and N e + 1 electrons, which can be exploited to obtain an approximate expression for the single-particle Green's functions at momentum k , energy ω and spin σ :

$$\text{varia-} \\ \text{ool to} \\ \text{single-} \\ \text{). \ \, \text{In} \\ \text{or.\ }$$

$$\text{$r$ cor-} \\ \text{$of$ de-} \\ \text{$e$ flux-} \\ \text{total} \\ \text{credit} \\ \text{$f$or$} \, \mathcal { C } ^ { h } \, \text{and} \, \mathcal { C } ^ { e } \, \text{donate to the hole and also return one} \, \text{$coron} \, \text{$coron} \, \text{$e$}$$

Here G h and G e denote the holes and electrons Green's functions, respectively, which are 4 × 4 matrices in the sublattice indices η and η ′ . The sum over n runs over a finite set of variational excitations Ψ h k,σ,n (Ψ e k,σ,n ) living in the sector with N e -1 ( N e +1) electrons; the respective variational energies are E h k,σ,n ( E e k,σ,n ). The real parameter ε introduces a small Lorentzian broadening.

The approximate excited states entering Eqs. (7) and (8) are constructed by a Rayleigh-Ritz variational approach in which the Hubbard Hamiltonian is diagonalized within a restricted subspace of variational states [48]. We adopt the basis set proposed in Ref. [49] and generalize it to the case of a non-Bravais lattice. For the sake of clarity, let us consider the example of hole excitations. The variational subspace at momentum k and spin σ = ↑ is formed by states of the form

$$\frac { 1 } { \sqrt { N _ { c } } } \sum _ { R } e ^ { - i k R } \, T _ { R } c _ { 0, \eta, \uparrow } B _ { 0, \eta } | \Psi _ { 0 } \rangle \,. \quad \quad ( 9 )$$

Here, B 0 ,η is a certain local operator that conserves the number of electrons and is centered at site (0 , η ); T R is the translation operator for the Bravais vector R . The various states forming the variational basis are then defined by considering different operators B 0 ,η . The simplest excitations in our basis are the ones in which B 0 ,η = 1, i.e. the bare excitations c k,η, ↑ | Ψ 0 ⟩ , and the ones with B 0 ,η = n 0 ,η, ↓ . In addition to these, we include a set of states with B 0 ,η = n i 0 ,η , ↑ n j 0 ,η , ↓ , where i 0 ,η and j 0 ,η denote two distinct sites in the neighborhood of site (0 , η ). Werestrict them to first-neighbors of (0 , η ) and sites lying across the square lattice diagonals with respect to (0 , η ) (excluding the diagonals perpendicular to the d bonds). The resulting subspace contains 176 variational states for each momentum and spin. The basis for electron excitations is defined analogously by applying a (uniform) particle-hole transformation to the operators in Eq. (9).

The solution of the generalized eigenvalue problem for the Hamiltonian (1) in the variational subspace of hole and electron excitations provides the energies and matrix elements necessary to compute the approximate Green's functions of Eqs. (7) and (8) (for details see Ref. [48]).

Figure 2. Panels (a-d): band structure of the non-interacting model [Eq. (1) with U = 0] for t = 1 , t ′ = 0 4, plotted along . the high symmetry lines of the Brillouin zone, shown in panel (f). The size of the dots represents the weights of the Bloch states onto the orbitals of Eqs. (11-14), e.g. |⟨ 0 | s k d † k,ν | 0 ⟩| 2 for the case of the s -orbital. Projections onto s , p x , p y and d orbitals are displayed in panels (a),(b),(c) and (d), respectively. In all panels the Fermi level at half-filling is represented by the horizontal grey dotted line. Panel (e): sublattice-resolved band structure. The color of the points corresponds to the value of e ν ( k ) [colorbar in panel (f)], with red and blue colors indicating localization of the Bloch states on even and odd sublattices, respectively. Panel (f): Brillouin zone and sublattice-resolved Fermi surface at half-filling. The Bloch states of the two lines forming the Fermi surface are localized on even/odd sublattice sites.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00841v2-figure-419f2616-6a1a-40ef-8e0e-5dd15411f68d.jpg)

Then, the spectral function can be obtained as

$$A _ { \sigma } ( k, \omega ) = - \frac { 1 } { \pi } \sum _ { \eta } \text{Im} [ G ^ { h } _ { \eta, \eta } ( k, \sigma, \omega ) + G ^ { e } _ { \eta, \eta } ( k, \sigma, \omega ) ] \,, \quad \text{itively} \quad \text{$\nu$ indi$} \\ \text{when} \quad \text{$line$} \quad \text{(} \quad \text{$10$)} \quad \text{general} \quad \text{$the$ $k$}$$

where we have taken the trace over the sublattice indices. The spectral function is thus approximated by a discrete set of delta functions (Lorentzians for finite ε ) at each k -point. In this respect, the present numerical approach is expected to be less effective for the description of incoherent continua, due to the finiteness of the variational basis. Still, a comparison with exact diagonalization calculations on small clusters showed that the finite set of variational excited states is capable of capturing the main spectral features, also within Mott insulating phases [48].

the bands ϵ k,ν shown in Fig. 2, with k = ( k , k x y ) and ν indicating the momentum and the band index, respectively. On the edges of the Brillouin zone, e.g. along the line ( π, π ) → ( π, 0), the two upper (lower) bands are degenerate. The third band from the bottom is flat along the k x = ± k y lines, i.e. the diagonals of the Brillouin zone, and displays a quadratic band touching with the second (fourth) band at the zone center (zone corners). For the case t ′ = 0 (square lattice limit) the partiallyflat bands can be unfolded and correspond to the nested Fermi surface at half-filling. For the remainder of the discussion, we assume that t ′ = 0.

## III. RESULTS

## A. Non-interacting band structure

We begin the discussion by looking at the band structure of the non-interacting model [Eq. (1) with U = 0]. For simplicity, let us drop the spin label at this stage, as the single-particle energies do not depend on σ . The diagonalization of the non-interacting Hamiltonian yields

̸

It is worth taking a look at the character of the Bloch states. We express their annihilation operators d k,ν as linear combinations of c k,η , i.e. d k,ν = ∑ η φ ν k,η c k,η . We are interested in projecting the various Bloch states at momentum k onto the s , p x , p y and d orbitals formed by

combining the different sublattices [15]

$$s _ { k } = \frac { 1 } { 2 } ( c _ { k, 0 } + c _ { k, 1 } + c _ { k, 2 } + c _ { k, 3 } ) \,, \quad \ \ ( 1 1 )$$

$$p _ { k } ^ { x } = \frac { \mu _ { 1 } } { \sqrt { 2 } } ( c _ { k, 0 } - c _ { k, 2 } ) \,, \quad \quad ( 1 2 ) \quad \underset { \mathfrak { C } ^ { \prime } } { \theta }$$

$$p _ { k } ^ { y } = \frac { \cdot _ { 1 } } { \sqrt { 2 } } ( c _ { k, 1 } - c _ { k, 3 } ) \,, \quad \quad \quad ( 1 3 ) \quad \ \ _ { \cdot }$$

$$d _ { k } = \frac { 1 } { 2 } ( c _ { k, 0 } - c _ { k, 1 } + c _ { k, 2 } - c _ { k, 3 } ) \,. \quad \ \ ( 1 4 ) \quad \ \ _ { \, }$$

The projected band structure of the non-interacting Hamiltonian is displayed in Fig. 2(a-d), with the size of the dots in each panel representing the weights of the Bloch states onto the different orbitals (11-14). We note that the lower-most and upper-most bands are mostly localized on the s and d orbitals, respectively. The second and third bands, instead, show contributions also from the p x and p y orbitals. Most importantly, the partially flat band along the Brillouin zone diagonals, which is cut by the Fermi level at half-filling, is entirely due to the p x ( p y ) orbital on the k x = k y ( k x = -k y ) line; this fact implies that the corresponding Bloch states are localized on the even (odd) sublattice, as shown also in Fig. 2 (e). Here the bands are plotted with a color code that reflects the value of e ν ( k ) = | φ ν k, 0 | 2 + | φ ν k, 2 | 2 , i.e. the degree of localization of the Bloch eigenvectors on the sites of the even sublattice.

The Fermi surface at half-filling coincides with the two diagonals of the Brillouin zone, k x = ± k y , each of them having an opposite sublattice character, as shown in Fig. 2(f). This aspect provides a simple argument for the onset of the altermagnetic band splitting in the ordered phase, as discussed in the next section. We remark that the partially flat band along the diagonal lines is present for any value of t /t ′ . Furthermore, as long as | t /t ′ | ≤ 2 the system remains metallic at half-filling and the Fermi surface is always the one shown in Fig. 2(f). For | t /t ′ | &gt; 2, a gap opens and the system becomes a band insulator.

## B. Altermagnetism at half-filling: metallic and insulating phases

̸

We turn to the discussion of the phase diagram of the Hubbard model Eq. (1) as obtained by variational Monte Carlo. At half-filling ( n = 1), the presence of a finite repulsive interaction U induces the onset of a staggered magnetic order (see Fig. 1), which is signalled by the stabilization of a finite value of the Zeeman field parameter h upon energy minimization. In the well-known square lattice limit ( t ′ = 0), any finite U value is sufficient to induce N´ eel order and turn the system into an insulator [50-52]. In this case, the optimization of the variational parameters yields χ ′ = 0 and h = 0, i.e. the auxiliary Hamiltonian of Eq. (5) is the one of an antiferromagnetic insulator.

Figure 3. Spin-resolved band structure ϵ σ k,ν of the auxiliary Hamiltonian of Eq. (5). Examples for an altermagnetic metal ( χ ′ = 0 4 . , h = 0 3) and an altermagnetic insulator . ( χ ′ = 0 4 . , h = 0 8) are shown in panels (a) and (b), . respectively. The Fermi surface of the metallic case of panel (a) is shown in (c-e), for different fillings. Red and blue colors indicate single-particle energies of ↑ and ↓ spins, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00841v2-figure-e84303ac-4f7a-412f-a6c2-f6802aa6a938.jpg)

The situation changes when t ′ is finite. Indeed, the inclusion of the Shastry-Sutherland diagonal hoppings reduces the symmetry of the Hamiltonian from p m 4 (square lattice) to p g 4 [54][55]. Therefore, even and odd sublattices, which correspond to opposite spin orientations of the staggered magnetic order, cease to be connected by translations or inversion. The remaining symmetries which transform the two magnetic sublattices into each other are C 4 rotations and the glide reflections with respect to the axes of first-neighbor bonds ( G x and G y shown in Fig. 1). Under these conditions, the staggered magnetic order induced by the Hubbard repulsion becomes altermagnetic. This is visible in the band structure of the auxiliary Hamiltonian (5) where, in addition to a non-zero Zeeman field h , a finite χ ′ is stabilized: the ↑ and ↓ bands split symmetrically in k -space, in a d -wave fashion, as shown in Fig. 3. The magnitude of the splitting can be traced back to the degree of localization of the paramagnetic bands on the different sublattices [cf Fig. 2, in particular panel (e)]. In this respect, a large splitting is observed for the second and third bands along the k x = ± k y lines, where the non-magnetic Bloch states are entirely localized on the p x and p y orbitals. Here, the effect of the Zeeman field is to rigidly shift the bands with opposite sublattice character in opposite directions. On the other hand, no splitting is present for the momenta which are mapped into themselves by rotations and glide

Figure 4. Schematic phase diagram of the Hubbard model (1) at half-filling. The hatched region delimits the metallic phase in which the magnetic order parameter m is very small and hard to distinguish from zero. The colormap shows the value of lim q → 0 q 2 /N q ( ), which is zero (finite) in a metallic (insulating) phase, i.e. can be taken as a proxy for the charge gap of the system [53].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00841v2-figure-23867fb2-9cf1-4f19-a90f-a00db69d78d6.jpg)

symmetries, namely the k x = 0 and k y = 0 lines and the zone boundaries (i.e., k x = π and k y = π lines). In this case, the combination of the aforementioned point group symmetries with time reversal protects the Kramers' degeneracy of the bands.

As a function of increasing U/t , the ground state first becomes an altermagnetic metal and then an altermagnetic insulator, as shown in the phase diagram of Fig. 4. The metallic/insulating nature of the ground state can be directly inferred by looking at the spectrum of the auxiliary Hamiltonian. More quantitatively, one can observe the smallq behavior of the densitydensity correlation function N q ( ) = ∑ η,η ′ ⟨ n -q,η ′ n q,η ⟩ / 4, where n q,η = N -1 / 2 c ∑ R n R,η e -iqR is the Fourier transform of the electron density on the sublattice η . If N q ( ) ∼ q ( q 2 ) for q → 0, the charge excitations are gapless (gapped) [53, 56]. To characterize the phase transition, we compute the staggered magnetization

$$m = \frac { 1 } { 2 N } \sum _ { i } ( - 1 ) ^ { \eta _ { i } } \langle ( n _ { i, \uparrow } - n _ { i, \downarrow } ) \rangle \,, \quad ( 1 5 ) \quad \text{$$

and the average double occupancy ⟨ n i, ↑ n i, ↓ ⟩ for different U/t and t /t ′ values, as shown in Fig. 5. Here ⟨·⟩ indicates the expectation value over the variational ground state | Ψ 0 ⟩ .

For small values of the diagonal hopping, i.e. t /t ′ ≤ 0 6, . the magnetization grows steadily as a function of U/t . Also the average double occupancy decreases smoothly, suggesting a continuous metal-insulator transition. While the opening of the charge gap can be reliably located, the critical U/t value necessary to induce mag-

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00841v2-figure-5c7d2708-4e0d-4e3f-93ef-6f3141a98198.jpg)

Figure 5. Staggered magnetization ( m , upper panel) and average double occupancy ( ⟨ n i, ↑ n i, ↓ ⟩ , lower panel) as a function of the Hubbard interaction U/t at half-filling. The various sets of data correspond to different values of t /t ′ . Empty (full) dots indicate that the ground state is metallic (insulating).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00841v2-figure-700ded79-f394-461a-b556-7731eac27a4c.jpg)

h

Figure 6. Variational energy landscape as a function of the Zeeman field parameter h . Results for t /t ′ = 1, at U/t = 6 1 . [panel (a)] and U/t = 6 3 [panel (b)]. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00841v2-figure-73a4a504-a350-406d-8a87-bcd326ba8e3e.jpg)

netic order is harder to estimate. We observe vanishingly small magnetic order at small U/t , i.e. a paramagnetic metallic region [hatched area in Fig. 4]. However, distinguishing a small m from zero is numerically challenging, and finite-size effects can play a relevant role, as can be seen in the t ′ = 0 case, where the value of m seems to vanish at U/t = 1. Increasing the hopping ratio t /t ′ induces magnetic frustration in the system and implies a slower growth of m as a function of U , which makes the numerical analysis harder. For this reason, the region of very small magnetic order becomes wider for larger val-

ues of t /t ′ and the latter behavior is compatible with the finite-temperature results of Ref. [33].

Interestingly, the transition between the altermagnetic metal and the insulator becomes of the first order for larger values of the diagonal hoppings ( t /t &gt; ′ 0 8), where . we observe a sudden jump in the magnetization and in the average double occupancy. This is further corroborated by the calculation of the energy landscape reported in Fig. 6, namely the variational ground state energy E 0 obtained by optimizing the wave function for different fixed values of the Zeeman field parameter h . The energy landscape shows two minima, at small and large h values, which swap their relative positions across the phase transition, showcasing its first-order character. It is worth noting that within a pure mean-field approach the metal-insulator transition is found to be continuous for all values of t /t ′ considered here. Hence, the Jastrow correlators play a relevant role in determining the firstorder behavior at large t /t ′ . These results are thus suggestive of the following scenario: for small t /t ′ , a meanfield-like metal-insulator transition takes place, with a smooth increase of the Zeeman field inducing altermagnetic order in the variational state and gradually leading to a Slater insulating state [57], which then crosses over to the Mott regime at large U ; for large t /t ′ , instead, the slow increase of the magnetic order at small U/t is not sufficient to open a charge gap before the Mott regime kicks in, and thus a first-order transition is observed, with an abrupt increase of m and a jump in the average double occupancy.

To conclude this section, we provide further evidence of the altermagnetic nature of the ordered phases. At the mean-field level, altermagnetism is evidenced by the compensated k -dependent splitting of ↑ and ↓ bands in the single-particle spectrum. Therefore, one may expect an analogous splitting of the kinetic contributions in the correlated ground state wave function. We thus compute the kinetic terms ⟨ c † iσ c jσ + c † jσ c iσ ⟩ at first neighbors, for the optimal variational wave functions at different t /t ′ and U/t values. The results are summarized in Fig. 7. We observe a different behavior in the two sets of square plaquettes that contain a diagonal bond. One set of plaquettes is characterized by strong (weak) kinetic terms in the σ = ↑ ( ↓ ) channel. The other set of plaquettes, which is connected to the latter by C 4 rotations, or G ,G x y glide reflections, displays the exact opposite behavior. Overall, this gives rise to a d -wave-like spin-polarization of the kinetic contributions in real space (around the center of C 4 rotations), analogous to the d -wave spin-splitting of the mean-field bands observed in momentum space. We can then take the ratio between the strong and weak kinetic terms on each bond as a real-space marker for altermagnetism, i.e.

$$R = \max ( r, 1 / r ) \ \text{ with } r = \frac { \langle c _ { i \uparrow } ^ { \dagger } c _ { j \uparrow } + c _ { j \uparrow } ^ { \dagger } c _ { i \uparrow } \rangle } { \langle c _ { i \downarrow } ^ { \dagger } c _ { j \downarrow } + c _ { j \downarrow } ^ { \dagger } c _ { i \downarrow } \rangle } \,. \quad ( 1 6 ) \quad \text{notat} \\ \text{and}$$

Indeed, for t ′ = 0, R = 1 for all U values, i.e. the kinetic terms are isotropic in spin and space, and the system is not altermagnetic, as expected. For finite t ′ , the altermagnetic ratio R first increases and then decreases as a function of U/t , denoting the onset of altermagnetism. We again observe a smooth behaviour for small values of t /t ′ , indicative of a continuous transition, with the maximum of R being located in the insulating phase. On the other hand, for large t /t ′ , the ratio jumps when the discontinuous metal-insulator transition is hit.

Figure 7. Altermagnetic ratio R [Eq. (16)] as a function of U/t , for different values of t /t ′ . Empty (full) dots indicate that the ground state is metallic (insulating). In the upper panel, a sketch of the ↑ and ↓ kinetic terms is shown, displaying a splitting between strong and weak contributions. A specular behavior is observed in the two sets of square plaquettes with oppositely oriented diagonals. R is defined as the ratio between the stronger and the weaker kinetic contributions on each first-neighbor bond.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00841v2-figure-0faecbb2-1d9d-4b92-8775-5d20b400f1a5.jpg)

## C. The effect of doping

We assess the stability of the altermagnetic phases when the system is doped away from half-filling. We consider the cases with U/t = 4 and U/t = 8, and track the evolution of the altermagnetic order as a function of δ = n -1. The results are summarized in Fig. 8, where the staggered magnetization m and the altermagnetic ratio R are plotted as a function of doping. In our notation, δ &lt; 0 ( &gt; 0) indicates hole (electron) doping. We note that for t ′ = 0 the system reduces to the standard Hubbard model on the square lattice and it is thus particle-hole symmetric [41]. For finite t ′ , instead, the

)

(

S Q glyph[triangleleft] N

2

80

60

40

20

0

3

1

100

80

t glyph[triangleleft]t

′

t glyph[triangleleft]t

′

t glyph[triangleleft]t

′

= 1 0

glyph[triangleright]

glyph[triangleright]

= 0 8

glyph[triangleright]

= 0 6

t glyph[triangleleft]t

′

t glyph[triangleleft]t

′

100

t glyph[triangleleft]t

′

= 0 8

glyph[triangleright]

glyph[triangleright]

= 0 6

glyph[triangleright]

= 0 4

100

t glyph[triangleleft]t

′

= 1 0

glyph[triangleright]

t glyph[triangleleft]t

′

= 1 0

glyph[triangleright]

t glyph[triangleleft]t

′

= 0 8

glyph[triangleright]

t glyph[triangleleft]t

100= 0 6

t glyph[triangleleft]t

′

= 0 4

glyph[triangleright]

′

= 0 2

glyph[triangleright]

t glyph[triangleleft]t

′

= 1 0

glyph[triangleright]

t glyph[triangleleft]t

′

= 0 8

glyph[triangleright]

t glyph[triangleleft]t

′

glyph[triangleright]

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00841v2-figure-509b4e5c-bc05-42a1-bf99-902c64c29b84.jpg)

80

Figure 8. Staggered magnetization m and altermagnetic ratio R at U/t = 4 [panels (a),(b)] and U/t = 8 [panels (c),(d)], as a function of doping δ . Results for different values of t /t ′ are shown. Negative (positive) values of δ refer to hole (electron) doping. The different sets of data of each panel are plotted on the same scale, which is specified in the upper-most plot (i.e., for t /t ′ = 1). The solid horizontal lines mark the ground reference values, m = 0 in panels (a),(c) and R = 1 in panel (b),(d). The dashed lines delimit the region with finite magnetic order. The regimes in which system is prone to phase separation are indicated by the hatched grey areas in panels (a),(c).

↦

↦

staggered particle-hole transformation, c i,σ → -( ) η i c † i,σ , maps the model at doping δ to the model at doping -δ , but changes t ′ → -t ′ . Thus the system ceases to be particle-hole symmetric, as clearly shown in the results of Fig. 8.

Let us first discuss the results at large U , i.e. U/t = 8 [Fig. 8(c,d)], where the system is in the altermagnetic insulating phase at half-filling for all t /t ′ values considered here. Away from half-filling, the system immediately becomes metallic and the staggered magnetization m starts decreasing as a function of doping, see Fig. 8(c). The region of stability of the altermagnetic order to hole doping shrinks when t /t ′ is increased, while in the case of electron doping it remains almost constant. For hole doping, the altermagnetic ratio R , shown in Fig. 8(d), first increases and then decreases, reaching R = 1 when m vanishes. The disappearence of the altermagnetic order with hole doping is smooth for small t /t ′ and abrupt for large t /t ′ . In presence of electron doping, we observe a similar behavior for the magnetization m , which jumps to zero for large t /t ′ . On the other hand, the altermagnetic ratio displays a smooth behavior for all t /t ′ values.

phase separation in the presence of hole doping, marked by the grey hatched areas in Fig. 8(c). On the contrary, for electron doping phase separation is found only for t /t ′ ≤ 0 2. . We point out, however, that variational approaches like the one employed in this work tend to overestimate the tendency to phase separate, as the accuracy of the variational Ansatz can differ for different electron densities n , thus affecting the above analysis based on the calculation of e δ ( ). In this respect, for instance, the authors of Ref. [59] showed that in the square lattice limit ( t ′ = 0), for relatively large U values (e.g., U/t = 8 considered here) the region of phase separation shrinks considerably when improving Jastrow-Slater variational results by performing more accurate fixed-node Monte Carlo calculations.

It is important to note that our variational results indicate that the altermagnetic metal at finite doping is unstable to phase separation. Considering as an example the case of hole doping, phase separation can be detected by inspecting the behavior of the energy per hole, e δ ( ) = [ E 0 ( δ ) -E 0 ( δ = 0)] / δ | | , as a function of δ [58, 59]. For a stable phase, e δ ( ) is monotonic, while in presence of phase separation (in a finite-sized system) it displays a minimum at a certain value of δ . This indicates that the system can gain energy by forming hole-rich and holepoor regions. At U/t = 8, we find that large regions of

We now turn to the results at smaller U values, i.e. U/t = 4, shown in Fig. 8(a,b). In this case, the system is insulating at half-filling for t /t ′ ≤ 0 6 and metallic for . for larger values of the diagonal hopping, as reported in Fig. 5. In presence of finite doping, the insulating phase at t /t ′ ≤ 0 6 immediately turns to a metallic altermag-. net. Then the altermagnetic order decreases both for hole and electron dopings, displaying a larger stability in the latter case. An interesting effect is observed when doping the altermagnetic metal found at half-filling for t /t ′ = 0 8 and . t /t ′ = 1. In case of electron doping, the altermagnetic order is enhanced, especially in the case t /t ′ = 1, where the magnetization m is vanishingly small at half-filling and clearly becomes finite for δ &gt; 0. Both m and R first increase and then smoothly decrease to zero as a function of δ . Finally, we note that phase separation at U/t = 4 is found only at small t /t ′ and, except for the case t ′ = 0, only in the case of hole doping.

Figure 9. Spectral function A σ ( k, ω ) at half-filling ( n = 1) for t /t ′ = 0 4. . Results for the altermagnetic metal [ U/t = 3, panel (a)] and the altermagnetic insulator [ U/t = 7 and 12, panels (b) and (c)] are shown. The panels are symmetrically split to show the σ = ↓ (left side) and σ = ↑ (right side) contributions along the same k -path, with blue and red color schemes, respectively. The dashed lines represent the spin-resolved bands of the optimal auxiliary Hamiltonian [Eq. (5)]. A Lorentzian broadening ε = 0 2 . t is used. Panels (d-f) show A ↓ ( k, ω ) and A ↑ ( k, ω ) at k = ( π/ , π/ 2 2), for the same U/t values of panels (a-c), respectively. The Fermi level is at ω = 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00841v2-figure-04e270fd-0865-45ce-9f6d-71b451111fef.jpg)

## D. Spectral functions

We conclude the results section by presenting the spectral function of the Hubbard model as obtained by variational Monte Carlo, as it provides the most direct evidence of altermagnetism [60]. Indeed, in the altermagnetic phases a splitting of the A ↑ ( k, ω ) and A ↓ ( k, ω ) components along certain k -directions is observed, analogously to the spin-split band structure of the auxiliary Hamiltonian of Eq. (5) (see Fig. 3). The symmetry of the altermagnetic order on the Shastry-Sutherland lattice implies that the spectral features for σ = ↑ and ↓ components are related to each other by, e.g., a C 4 rotation in k -space.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00841v2-figure-b468e401-f504-4958-961f-75f5f307f906.jpg)

In Fig. 9, we show the spectral function at half-filling for t /t ′ = 0 4 and different . values of U/t . The spinresolved A ↑ ( k, ω ) and A ↓ ( k, ω ) are plotted along the k x = k y diagonal of the Brillouin zone. In the metallic regime at U/t = 3, the spectrum closely resembles the one of the auxiliary Hamiltonian, with well-defined quasiparticle bands, especially close to the Fermi level ( ω = 0). The altermagnetic splitting is showcased by, e.g., the second and third bands, one dispersive and one flat, which are pushed in opposite energy directions in the ↑ and ↓ channels by the presence of finite altermagnetic order. Increasing U/t and entering the insulating regime, we observe the gap opening and a progressive broadening of the spectral features into the upper and lower Hubbard bands. The altermagnetic splitting is visible also within the Mott insulating phase, where the excitations are incoherent, e.g. for U/t = 12. In Fig 9(d-f), we report a comparison of ↑ and ↓ spectral functions at fixed momentum k = ( π/ , π/ 2 2) which helps visualizing the splitting for all U/t values.

Figure 10. The same as in Fig. 9, for t /t ′ = 0 4, . U/t = 8 and two different fillings in the electron doped case. The results of panel (a) are obtained in the altermagnetic metallic phase at small doping ( N e = 282, δ ≈ 0 1); . the results in panel (b) are in the paramagnetic metallic phase at large doping ( N e = 326, δ ≈ 0 27). .

Finally, in Fig. 10, spectral functions at finite electron doping are presented. For small δ , the system is in the altermagnetic metallic phase and the spectral function displays the splitting between the spin channels [see Fig. 10(a)]. While close to the Fermi energy the dispersion is well-captured by the band structure of the auxiliary Hamiltonian, away from ω = 0 the signal becomes incoherent and the altermagnetic splitting is seemingly reduced. On the other hand, at large electron dopings altermagnetism is suppressed. As shown in Fig. 10(b), the A ↑ ( k, ω ) and A ↓ ( k, ω ) spectral functions are symmetric and no spin-splitting is observed.

## IV. CONCLUSIONS

In this work, we investigated the properties of the Hubbard model on the Shastry-Sutherland lattice in the regime of moderate frustration, as it provides a minimal setting to study altermagnetism in a correlated system. We made use of a variational Monte Carlo approach, based on Jastrow-Slater wave functions, to investigate ground state (and spectral) properties beyond mean-field-like approximations.

At half-filling, when the Shastry-Sutherland diagonal hopping t ′ is finite (i.e., away from the square-lattice limit), the system is found to be first metallic for a finite range of U/t , and then insulating for larger values of the Hubbard interaction. Most importantly, the onsite repulsion induces the onset of a staggered magnetic order, which bears altermagnetic properties due to the symmetry of the underlying lattice. We studied the metalinsulator transition between altermagnetic phases at halffilling, as a function of U/t . While for small values of the diagonal hoppings, the transition is found to be continuous and Slater-like, for sufficiently large t /t ′ a first-order jump of the magnetization is detected in correspondence of the gap opening. The change in the nature of the transition testifies the importance of electronic correlations, as it is not captured by a standard mean-field approach. In addition to the magnetic order parameter, we provided a static measure of altermagnetism by computing the difference of the first-neighbors kinetic terms in the two spin channels (parallel to the magnetic order direction). We investigated also the effects of doping, assessing the stability of the altermagnetic order away from half-filling.

- [1] L. Smejkal, R. Gonz´lez-Hern´ndez, T. Jungwirth, and ˇ a a J. Sinova, Crystal time-reversal symmetry breaking and spontaneous hall effect in collinear antiferromagnets, Science Advances 6 , eaaz8809 (2020).
- [2] L. Smejkal, J. Sinova, and T. Jungwirth, Emerging reˇ search landscape of altermagnetism, Phys. Rev. X 12 , 040501 (2022).
- [3] Y. Noda, K. Ohno, and S. Nakamura, Momentumdependent band spin splitting in semiconducting mno2: a density functional calculation, Phys. Chem. Chem. Phys. 18 , 13294 (2016).
- [4] S. Hayami, Y. Yanagi, and H. Kusunose, Momentumdependent spin splitting by collinear antiferromagnetic ordering, Journal of the Physical Society of Japan 88 , 123702 (2019).
- [5] L.-D. Yuan, Z. Wang, J.-W. Luo, E. I. Rashba, and A. Zunger, Giant momentum-dependent spin splitting in centrosymmetric lowz antiferromagnets, Phys. Rev. B 102 , 014422 (2020).
- [6] I. I. Mazin, K. Koepernik, M. D. Johannes, R. Gonz´leza Hern´ andez, and L. ˇ mejkal, Prediction of unconventional S magnetism in doped FeSb , Proceedings of the National 2 Academy of Sciences 118 , e2108924118 (2021).

Finally, we computed the spin-resolved spectral function, providing evidence of the altermagnetic splitting of the spectral features, both in the metallic ground state with well-defined quasiparticle excitations, and in the highlycorrelated Mott regime.

These results present a new avenue to explore correlation-induced altermagnetic phases. Possible material candidates where such phases could be realized are, for instance, α - or κ -(BEDT-TTF) X [34-36] charge2 transfer salts, whose underlying lattice structure can give rise to a distorted version of the Shastry-Sutherland lattice. Although κ -(BEDT-TTF) X salts are actually 2 closer to the regime of strong dimer hopping (i.e., t ′ ≫ t ) and away from half-filling [61, 62], the versatility of these materials as a function of pressure, strain and doping may open up the possibility to tune their properties towards altermagnetic phases related to the ones observed in this work. On the theoretical side, instead, the possible presence of superconductivity in the doped regime warrants further investigation, also in view of the possibility of realizing pair-density wave states [12, 13, 38]. A more accurate variational approach including backflow correlations could be suitable for this purpose [63].

## ACKNOWLEDGMENTS

We would like to thank F. Becca, J. Knolle, I. I. Mazin, J. Profe, J. Sinova, L. ˇ mejkal S and A. Valadkhani for useful discussions, and P. P. Stavropoulos for providing mean-field results for comparison. We gratefully acknowledge support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) for funding through TRR 288-422213477 (project B05).

- [7] L. ˇ mejkal, J. Sinova, and T. Jungwirth, Beyond convenS tional ferromagnetism and antiferromagnetism: A phase with nonrelativistic spin and crystal rotation symmetry, Phys. Rev. X 12 , 031042 (2022).
- [8] S. Bhowal and N. A. Spaldin, Ferroically ordered magnetic octupoles in d -wave altermagnets, Phys. Rev. X 14 , 011019 (2024).
- [9] R. M. Fernandes, V. S. de Carvalho, T. Birol, and R. G. Pereira, Topological transition from nodal to nodeless zeeman splitting in altermagnets, Phys. Rev. B 109 , 024404 (2024).
- [10] L. Bai, W. Feng, S. Liu, L. Smejkal, Y. Mokrousov, and ˇ Y. Yao, Altermagnetism: Exploring new frontiers in magnetism and spintronics, arXiv preprint arXiv:2406.02123 (2024).
- [11] D. Zhu, Z.-Y. Zhuang, Z. Wu, and Z. Yan, Topological superconductivity in two-dimensional altermagnetic metals, Phys. Rev. B 108 , 184505 (2023).
- [12] D. Chakraborty and A. M. Black-Schaffer, Zero-field finite-momentum and field-induced superconductivity in altermagnets (2024), arXiv:2309.14427 [cond-mat.suprcon].
- [13] G. Sim and J. Knolle, Pair density waves and supercurrent diode effect in altermagnets (2024),

- arXiv:2407.01513 [cond-mat.supr-con].
- [14] T. A. Maier and S. Okamoto, Weak-coupling theory of neutron scattering as a probe of altermagnetism, Phys. Rev. B 108 , L100402 (2023).
- [15] A. Bose, S. Vadnais, and A. Paramekanti, Altermagnetism and superconductivity in a multiorbital t-j model (2024), arXiv:2403.17050 [cond-mat.str-el].
- [16] P. Das, V. Leeb, J. Knolle, and M. Knap, Realizing altermagnetism in fermi-hubbard models with ultracold atoms, Phys. Rev. Lett. 132 , 263402 (2024).
- [17] V. Leeb, A. Mook, L. Smejkal, and J. Knolle, Spontaˇ neous formation of altermagnetism from orbital ordering, Phys. Rev. Lett. 132 , 236701 (2024).
- [18] B. Sriram Shastry and B. Sutherland, Exact ground state of a quantum mechanical antiferromagnet, Physica B+C 108 , 1069 (1981).
- [19] M. Albrecht and F. Mila, First-order transition between magnetic order and valence bond order in a 2d frustrated heisenberg model, Europhysics Letters 34 , 145 (1996).
- [20] S. Miyahara and K. Ueda, Exact dimer ground state of the two dimensional heisenberg spin system srcu 2 (bo ) , 3 2 Phys. Rev. Lett. 82 , 3701 (1999).
- [21] A. Koga and N. Kawakami, Quantum phase transitions in the shastry-sutherland model for srcu (bo ) , 2 3 2 Phys. Rev. Lett. 84 , 4461 (2000).
- [22] Z. Weihong, C. J. Hamer, and J. Oitmaa, Series expansions for a heisenberg antiferromagnetic model for srcu 2 (bo ) , Phys. Rev. B 3 2 60 , 6608 (1999).
- [23] P. Corboz and F. Mila, Tensor network study of the shastry-sutherland model in zero magnetic field, Phys. Rev. B 87 , 115144 (2013).
- [24] J. Yang, A. W. Sandvik, and L. Wang, Quantum criticality and spin liquid phase in the shastry-sutherland model, Phys. Rev. B 105 , L060409 (2022).
- [25] J. Y. Lee, Y.-Z. You, S. Sachdev, and A. Vishwanath, Signatures of a deconfined phase transition on the shastrysutherland lattice: Applications to quantum critical srcu 2 (bo ) , Phys. Rev. X 3 2 9 , 041037 (2019).
- [26] L. L. Viteritti, R. Rende, A. Parola, S. Goldt, and F. Becca, Transformer wave function for the shastrysutherland model: emergence of a spin-liquid phase (2024), arXiv:2311.16889 [cond-mat.str-el].
- [27] H. Kageyama, K. Yoshimura, R. Stern, N. V. Mushnikov, K. Onizuka, M. Kato, K. Kosuge, C. P. Slichter, T. Goto, and Y. Ueda, Exact dimer ground state and quantized magnetization plateaus in the two-dimensional spin system srcu (bo ) , Phys. Rev. Lett. 2 3 2 82 , 3168 (1999).
- [28] S. Miyahara and K. Ueda, Theory of the orthogonal dimer heisenberg spin model for srcu2 (bo3)2, Journal of Physics: Condensed Matter 15 , R327 (2003).
- [29] M. E. Zayed, C. R¨ uegg, J. Larrea J., A. M. L¨ auchli, C. Panagopoulos, S. S. Saxena, M. Ellerby, D. F. McMorrow, T. Str¨ssle, a S. Klotz, G. Hamel, R. A. Sadykov, V. Pomjakushin, M. Boehm, M. Jim´ enezRuiz, A. Schneidewind, E. Pomjakushina, M. Stingaciu, K. Conder, and H. M. Rønnow, 4-spin plaquette singlet state in the shastry-sutherland compound srcu2(bo3)2, Nature Physics 13 , 962 (2017).
- [30] J. Guo, G. Sun, B. Zhao, L. Wang, W. Hong, V. A. Sidorov, N. Ma, Q. Wu, S. Li, Z. Y. Meng, A. W. Sandvik, and L. Sun, Quantum phases of srcu (bo ) 2 3 2 from high-pressure thermodynamics, Phys. Rev. Lett. 124 , 206602 (2020).
- [31] J. L. Jim´ enez, S. P. G. Crone, E. Fogh, M. E. Zayed, R. Lortz, E. Pomjakushina, K. Conder, A. M. L¨ auchli, L. Weber, S. Wessel, A. Honecker, B. Normand, C. R¨ uegg, P. Corboz, H. M. Rønnow, and F. Mila, A quantum magnetic analogue to the critical point of water, Nature 592 , 370 (2021).
- [32] Z. Shi, S. Dissanayake, P. Corboz, W. Steinhardt, D. Graf, D. M. Silevitch, H. A. Dabkowska, T. F. Rosenbaum, F. Mila, and S. Haravifard, Discovery of quantum phases in the shastry-sutherland compound srcu2(bo3)2 under extreme conditions of field and pressure, Nature Communications 13 , 2301 (2022).
- [33] H.-D. Liu, Y.-H. Chen, H.-F. Lin, H.-S. Tao, and W.M. Liu, Antiferromagnetic metal and mott transition on shastry-sutherland lattice, Scientific Reports 4 , 4829 (2014).
- [34] N. Tajima, M. Tamura, Y. Nishio, K. Kajita, and Y. Iye, Transport property of an organic conductor α -(bedt-ttf) 2 i 3 under high pressure-discovery of a novel type of conductor, Journal of the Physical Society of Japan 69 , 543 (2000).
- [35] K. Kanoda, Electron correlation, metal-insulator transition and superconductivity in quasi-2d organic systems, (et)2x, Physica C: Superconductivity 282-287 , 299 (1997), materials and Mechanisms of Superconductivity High Temperature Superconductors V.
- [36] K. Riedl, E. Gati, and R. Valent´ ı, Ingredients for generalized models of κ -phase organic charge-transfer salts: A review, Crystals 12 , 10.3390/cryst12121689 (2022).
- [37] M. Naka, S. Hayami, H. Kusunose, Y. Yanagi, Y. Motome, and H. Seo, Spin current generation in organic antiferromagnets, Nature Communications 10 , 4305 (2019).
- [38] S. Sumita, M. Naka, and H. Seo, Fulde-ferrell-larkinovchinnikov state induced by antiferromagnetic order in κ -type organic conductors, Phys. Rev. Res. 5 , 043171 (2023).
- [39] T. Misawa and M. Naka, Correlated zak insulator in organic antiferromagnets, Phys. Rev. B 108 , L081120 (2023).
- [40] J. Hubbard, Electron correlations in narrow energy bands, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 276 , 238 (1963).
- [41] D. P. Arovas, E. Berg, S. A. Kivelson, and S. Raghu, The hubbard model, Annual Review of Condensed Matter Physics 13 , 239 (2022).
- [42] F. Becca and S. Sorella, Quantum Monte Carlo Approaches for Correlated Systems (Cambridge University Press, 2017).
- [43] M. C. Gutzwiller, Effect of correlation on the ferromagnetism of transition metals, Phys. Rev. Lett. 10 , 159 (1963).
- [44] M. Capello, F. Becca, M. Fabrizio, S. Sorella, and E. Tosatti, Variational description of mott insulators, Phys. Rev. Lett. 94 , 026406 (2005).
- [45] M. Capello, F. Becca, S. Yunoki, and S. Sorella, Unconventional metal-insulator transition in two dimensions, Phys. Rev. B 73 , 245116 (2006).
- [46] The Jastrow factors pseudopotentials, v i,j and u i,j , are assumed to depend only on the Euclidean distance between sites. For this purpose, we adopt the 'noncommensurate' geometry of the Shastry-Sutherland lattice, i.e. the arrangement of magnetic Cu 2+ ions in SrCu (BO ) 2 3 2 [27].

- [47] S. Sorella, Wave function optimization in the variational monte carlo method, Phys. Rev. B 71 , 241103 (2005).
- [48] M. Charlebois and M. Imada, Single-particle spectral function formulated and calculated by variational monte carlo method with application to d -wave superconducting state, Phys. Rev. X 10 , 041023 (2020).
- [49] P. Rosenberg, D. S´ en´chal, e A.-M. S. Tremblay, and M. Charlebois, Fermi arcs from dynamical variational monte carlo, Phys. Rev. B 106 , 245132 (2022).
- [50] J. E. Hirsch, Two-dimensional hubbard model: Numerical simulation study, Phys. Rev. B 31 , 4403 (1985).
- [51] S. R. White, D. J. Scalapino, R. L. Sugar, E. Y. Loh, J. E. Gubernatis, and R. T. Scalettar, Numerical study of the two-dimensional hubbard model, Phys. Rev. B 40 , 506 (1989).
- [52] J. P. F. LeBlanc, A. E. Antipov, F. Becca, I. W. Bulik, G. K.-L. Chan, C.-M. Chung, Y. Deng, M. Ferrero, T. M. Henderson, C. A. Jim´ enez-Hoyos, E. Kozik, X.-W. Liu, A. J. Millis, N. V. Prokof'ev, M. Qin, G. E. Scuseria, H. Shi, B. V. Svistunov, L. F. Tocchio, I. S. Tupitsyn, S. R. White, S. Zhang, B.-X. Zheng, Z. Zhu, and E. Gull (Simons Collaboration on the Many-Electron Problem), Solutions of the two-dimensional hubbard model: Benchmarks and results from a wide range of numerical algorithms, Phys. Rev. X 5 , 041041 (2015).
- [53] L. F. Tocchio, F. Becca, and C. Gros, Backflow correlations in the hubbard model: An efficient tool for the study of the metal-insulator transition and the largeu limit, Phys. Rev. B 83 , 195138 (2011).
- [54] K. Parshukov, R. Wiedmann, and A. P. Schnyder, Topological responses from gapped weyl points in 2d altermagnets (2024), arXiv:2403.09520 [cond-mat.mes-hall].
- [55] The spatial symmetries of the problem are those of the wallpaper group p g 4 since the Hamiltonian under in-
- vestigation is strictly two-dimensional. A spin-wallpaper group describing the altermagnetically ordered phases can be defined by supplementing C 4 rotations and G ,G x y glide symmetries with a spin-flip operation (cf. group #12 3 of Ref. [54]). .
- [56] R. P. Feynman and M. Cohen, Energy spectrum of the excitations in liquid helium, Phys. Rev. 102 , 1189 (1956).
- [57] J. C. Slater, Magnetic effects and the hartree-fock equation, Phys. Rev. 82 , 538 (1951).
- [58] V. J. Emery, S. A. Kivelson, and H. Q. Lin, Phase separation in the t-j model, Phys. Rev. Lett. 64 , 475 (1990).
- [59] L. F. Tocchio, F. Becca, and S. Sorella, Hidden mott transition and largeu superconductivity in the twodimensional hubbard model, Phys. Rev. B 94 , 195126 (2016).
- [60] For the calculations of the spectral function we use periodic boundary conditions both in the a 1 and a 2 directions, in order to preserve the point group symmetry of the Shastry-Sutherland lattice. In the case t /t ′ = 0 4 and . U/t = 3, the system is metallic at half-filling and the spectrum of H aux is open-shell. For computational reasons, we split the degeneracy at the Fermi level by adding a tiny hopping of 10 -6 across the diagonals of the empty squares (in H aux ). This has no visible effects on the physical results.
- [61] T. Koretsune and C. Hotta, Evaluating model parameters of the κ - and β ′ -type mott insulating organic solids, Phys. Rev. B 89 , 045102 (2014).
- [62] D. Guterding, M. Altmeyer, H. O. Jeschke, and R. Valent´ ı, Near-degeneracy of extended s + d x 2 -y 2 and d xy order parameters in quasi-two-dimensional organic superconductors, Phys. Rev. B 94 , 024515 (2016).
- [63] L. F. Tocchio, F. Becca, A. Parola, and S. Sorella, Role of backflow correlations for the nonmagnetic phase of the t -t ′ hubbard model, Phys. Rev. B 78 , 041101 (2008).
Francesco Ferrari Roser Valenti
Altermagnetism on the Shastry-Sutherland lattice
Motivated by recent developments on altermagnetism, we investigate the
Hubbard model on the Shastry-Sutherland lattice, where the onsite repulsion
between electrons induces the onset of a staggered magnetic order in which
opposite magnetic sublattices are related to each other only by rotations and
glide reflections, and not by inversion nor translations. As a consequence, the
magnetic phase displays an altermagnetic character, i.e. a finite
(non-relativistic) Zeeman splitting of the electronic excitations, despite the
absence of a net magnetization. By means of a variational Monte Carlo approach
based on Jastrow-Slater wave functions, which accounts for electronic
correlations beyond the single-particle approximation, we study how $d$-wave
altermagnetism shows up in the ground state properties and in the spectral
function of the system. We characterize the metal-insulator transition between
altermagnetic phases at half-filling and the stability of altermagnetism as a
function of doping. The calculation of the single-particle spectral function
displays the spin-split nature of the electronic excitations, also within the
Mott insulating regime. We discuss possible realizations of the model in the
context of organic $\alpha$- or $\kappa$-(BEDT-TTF)$_2$X charge transfer salts.