Search

Kagome materials manifest rich physical properties due to the emergence of abundant electronic phases. Here, we carry out a high-throughput first-principles study of the kagome 1:6:6 family MT$_6$Z$_6$ materials in space group 191, focusing on their phonon instability and electronic flat bands. Different MT$_6$Z$_6$ kagome candidates reveal a remarkable variety of kagome flat bands ranging from unfilled, partially filled, to fully filled. Notably, the Mn/Fe-166 compounds exhibit partially filled flat bands with a pronounced sharp peak in the density of states near the Fermi level, leading to magnetic orders that polarize the bands and stabilize the otherwise unstable phonon. When the flat bands are located away from the Fermi level, we find a large number of phonon instabilities, which can be classified into three types, based on the phonon dispersion and vibrational modes. Type-I instabilities involve the in-plane distortion of kagome nets, while type-II and type-III present out-of-plane distortion of trigonal M and Z atoms. We take MgNi$_6$Ge$_6$ and HfNi$_6$In$_6$ as examples to illustrate the possible CDW structures derived from the emergent type-I and type-II instabilities. The type-I instability in MgNi$_6$Ge$_6$ suggests a nematic phase transition, governed by the local twisting of kagome nets. The type-II instability in HfNi$_6$In$_6$ may result in a hexagonal-to-orthorhombic transition, offering insight into the formation of MT$_6$Z$_6$ in other space groups. Additionally, the predicted ScNb$_6$Sn$_6$ is analyzed as an example of the type-III instability. Our predictions suggest a vast kagome family with rich properties induced by the flat bands, possible CDW transitions, and their interplay with magnetism., Comment: 14 pages, 7 figures, with 1000 pages of additional supplemental materials

Electronic ordering is prevalent in correlated systems, which commonly exhibit competing interactions. Here, we use x-ray diffraction to show that the charge density wave transition temperature of FeGe increases with pressure and evolves towards a $\sqrt{3}\times\sqrt{3}$ periodic lattice modulation, $\mathbf{q}$$^*$=$\left(\frac{1}{3}\ \frac{1}{3}\ \frac{1}{2}\right)$. In the pressure interval between 4$<$$p$$<$12 GPa both orders coexist and the spatial extent of the $\sqrt{3}\times\sqrt{3}$ order at high pressure becomes nearly long-range, $\sim$30 unit cells, while the correlation length of the 2$\times$2 phase remains shorter-ranged. The $\sqrt{3}\times\sqrt{3}$ phase is the ground state above 15 GPa, consistent with harmonic DFT calculations that predict a dimerization induced $\sqrt{3}\times\sqrt{3}$ order without phonon softening. The pressure dependence of the integrated intensities of $\mathbf{q}$$_\mathrm{CDW}=\left(\frac{1}{2}\ 0\ \frac{1}{2}\right)$ and $\mathbf{q}$$^*$ indicates a competition between the 2$\times$2 and $\sqrt{3}\times\sqrt{3}$ and demonstrates that the ground state of FeGe is characterized by a rich landscape of metastable/fragile phases. We discuss possible scenarios based on an order-disorder transformation and the formation of Friedel oscillations.

The intrinsic frustrated nature of a kagome lattice is amenable to the realization of exotic phases of matter, such as quantum spin liquids or spin ices, and more recently the multiple-$\mathrm{\textbf{q}}$ charge density waves (CDW) in the kagome metals. Despite intense efforts to understand the mechanism driving the electronic modulations, its origin is still unknown and hindered by competing interactions and intertwined orders. Here, we identify a dimerization-driven 2D hexagonal charge-diffuse precursor in the antiferromagnetic kagome metal FeGe and demonstrate that the fraction of dimerized/undimerized states is the relevant order parameter of the multiple-$\mathrm{\textbf{q}}$ CDW of a continuous phase transition. The pretransitional charge fluctuations with propagation vector $\mathrm{\textbf{q}=\textbf{q}_M}$ at T$_{\mathrm{CDW}}$$<$T$<$T$^*$(125 K) are anisotropic, hence holding a quasi-long-range bond-orientational order. The broken translational symmetry emerges from the anisotropic diffuse precursor, akin to the Ising scenario of antiferromagnetic triangular lattices. The temperature and momentum dependence of the critical scattering show parallels to the stacked hexatic $\mathrm{B}$-phases reported in liquid crystals and transient states of CDWs and highlight the key role of the topological defect-mediated melting of the CDW in FeGe.

Meeting of non-trivial topology with magnetism results in novel phases of matter, such as Quantum Anomalous Hall (QAH) or axion insulator phases. Even more exotic states with high and tunable Chern numbers are expected at the contact of intrinsic magnetic topological insulators (IMTIs) and 2D topological insulators (TIs).Here we synthesize a heterostructures composed of 2D TI and 3D IMTIs, specifically of bismuth bilayer on top of MnBi$_2$Te$_4$-family of compounds and study their electronic properties by means of angle-resolved photoelectron spectroscopy (ARPES) and density functional theory (DFT). The epitaxial interface is characterized by hybridized Bi and IMTI electronic states. The Bi bilayer-derived states on different members of MnBi$_2$Te$_4$-family of materials are similar, except in the region of mixing with the topological surface states of the substrate. In that region, the new, substrate dependent interface Dirac state is observed. Our \emph{ab initio} calculations show rich interface phases with emergence of exchange split 1D edge states, making the Bi/IMTI heterostructures promising playground for observation of novel members in the family of quantum Hall effects.

The emergence of correlated phenomena arising from the combination of 1$\mathrm{T}$ and 1$\mathrm{H}$ van der Waals layers is the focus of intense research. Here, we synthesize a novel self-stacked 6$\mathrm{R}$ phase in NbSeTe, showing a perfect alternating 1T and 1H layers that grow coherently along the c-direction, as revealed by scanning transmission electron microscopy. Angle resolved photoemission spectroscopy shows a mixed contribution of the trigonal and octahedral Nb bands to the Fermi level. Diffuse scattering reveals temperature-independent short-range charge fluctuations with propagation vector $\mathrm{q_{CO}}$=(0.25,0), derived from the condensation of a longitudinal mode in the 1T layer. We observe that ligand disorder quenches the formation of a charge density wave. Magnetization measurements suggest the presence of an inhomogeneous, short-range magnetic order, further supported by the absence of a clear phase transition in the specific heat. These experimental analyses in combination with \textit{ab initio} calculations indicate that the ground state of 6$\mathrm{R}$-NbSeTe is described by a statistical distribution of short-range charge-modulated and spin-correlated regions driven by ligand disorder. Our results devise a route to synthesize 1$\mathrm{T}$-1$\mathrm{H}$ self-stacked bulk heterostructures to study emergent phases of matter., Comment: 12 pages, including Supplementary Information. 4 figures + 6 supplementary figures

The electronic structure and interactions of kagome materials such as 1:1 (FeGe class) and 1:6:6 (MgFe$_6$Ge$_6$ class) are complicated and involve many orbitals and bands at the Fermi level. Current theoretical models treat the systems in an $s$-orbital kagome representation, unsuited and incorrect both quantitatively and qualitatively to the material realities. In this work, we lay the basis of a faithful framework of the electronic model for this large class of materials. We show that the complicated ``spaghetti" of electronic bands near the Fermi level can be decomposed into three groups of $d$-Fe orbitals coupled to specific Ge orbitals. Such decomposition allows for a clear analytical understanding (leading to different results than the simple $s$-orbital kagome models) of the flat bands in the system based on the $S$-matrix formalism of generalized bipartite lattices. Our three minimal Hamiltonians can reproduce the quasi-flat bands, van Hove singularities, topology, and Dirac points close to the Fermi level, which we prove by extensive ab initio studies. We also obtain the interacting Hamiltonian of $d$ orbitals in FeGe using the constraint random phase approximation (cRPA) method. We then use this as a fundamental ``LEGO"-like building block for a large family of 1:6:6 kagome materials, which can be obtained by doubling and perturbing the FeGe Hamiltonian. We applied the model to its kagome siblings FeSn and CoSn, and also MgFe$_6$Ge$_6$. Our work serves as the first complete framework for the study of the interacting phase diagram of kagome compounds., Comment: 5+50 pages, 3+16 figures, previously submitted. This is the second paper of a series on kagome materials. See also the first paper: arXiv:2305.15469

We derive a general method for determining the necessary and sufficient conditions for exact factorization $|\Psi\rangle=\otimes_p |\psi_p\rangle$ of an eigenstate of a many-body Hamiltonian $H$, based on the quantum covariance matrix of the relevant local operators building the Hamiltonian. The "site" $p$ can be either a single component or a group of subsystems. The formalism is then used to derive exact dimerization and clusterization conditions in spin systems, covering from spin-$s$ singlets and clusters coupled to $0$ total spin to general nonmaximally entangled spin-$s$ dimers (generalized singlets). New results for field induced dimerization in anisotropic $XYZ$ arrays under a magnetic field are obtained., Comment: 16 pages, 6 figures

In a series of recent reports, doped lead apatite (LK-99) has been proposed as a candidate ambient temperature and pressure superconductor. However, from both an experimental and theoretical perspective, these claims are largely unsubstantiated. To this end, our synthesis and subsequent analysis of an LK-99 sample reveals a multiphase material that does not exhibit high-temperature superconductivity. We study the structure of this phase with single-crystal X-ray diffraction (SXRD) and find a structure consistent with doped $\text{Pb}_{10}(\text{PO}_4)_6(\text{OH})_2$. However, the material is transparent which rules out a superconducting nature. From ab initio defect formation energy calculations, we find that the material likely hosts $\text{OH}^-$ anions, rather than divalent $\text{O}^{2-}$ anions, within the hexagonal channels and that Cu substitution is highly thermodynamically disfavored. Phonon spectra on the equilibrium structures reveal numerous unstable phonon modes. Together, these calculations suggest it is doubtful that Cu enters the structure in meaningful concentrations, despite initial attempts to model LK-99 in this way. However for the sake of completeness, we perform ab initio calculations of the topology, quantum geometry, and Wannier function localization in the Cu-dominated flat bands of four separate doped structures. In all cases, we find they are atomically localized by irreps, Wilson loops, and the Fubini-Study metric. It is unlikely that such bands can support strong superfluidity, and instead are susceptible to ferromagnetism (or out-of-plane antiferromagnetism) at low temperatures, which we find in ab initio studies. In sum, $\text{Pb}_{9}\text{Cu}(\text{PO}_4)_6(\text{OH})_2$ could more likely be a magnet, rather than an ambient temperature and pressure superconductor., Comment: 39 pages including appendices. Updated defect calculations and energy-dispersive X-ray spectroscopy data

In order to characterize in detail the charge density wave (CDW) transition of 1$T$-VSe$_2$, its electronic structure and lattice dynamics are comprehensively studied by means of x-ray diffraction, angle resolved photoemission (ARPES), diffuse and inelastic x-ray scattering (IXS), and state-of-the-art first principles density functional theory calculations. Resonant elastic x-ray scattering (REXS) does not show any resonant enhancement at either V or Se K-edges, indicating that the CDW peak describes a purely structural modulation of the electronic ordering. ARPES identifies (i) a pseudogap at T$>$T$_{CDW}$, which leads to a depletion of the density of states in the $ML-M'L'$ plane at T$<$T$_{CDW}$, and (ii) anomalies in the electronic dispersion reflecting a sizable impact of phonons on it. A diffuse scattering precursor, characteristic of soft phonons, is observed at room temperature (RT) and leads to the full collapse of the low-energy phonon ($\omega_1$) with propagation vector (0.25 0 -0.3) r.l.u. We show that the frequency and linewidth of this mode are anisotropic in momentum space, reflecting the momentum dependence of the electron-phonon interaction (EPI), hence demonstrating that the origin of the CDW is, to a much larger extent, due to the momentum dependence EPI with a small contribution from nesting. The pressure dependence of the $\omega_1$ soft mode remains nearly constant up to 13 GPa at RT, with only a modest softening before the transition to the high-pressure monoclinic $C2/m$ phase. The wide set of experimental data are well captured by our state-of-the art first-principles anharmonic calculations with the inclusion of van der Waals (vdW) corrections in the exchange-correlation functional. The description of the electronics and dynamics of VSe$_2$ reported here adds important pieces of information to the understanding of the electronic modulations of TMDs.

The emergence of symmetry-breaking orders such as ferromagnetism and the weak interlayer bonding in van der Waals materials, offers a unique platform to engineer novel heterostructures and tune transport properties like thermal conductivity. Here, we report the experimental and theoretical study of the cross-plane thermal conductivity, $\kappa_\perp$, of the van der Waals 2D ferromagnet Fe$_3$GeTe$_2$. We observe a non-monotonic increase of $\kappa_\perp$ with the thickness and a large suppression in artificially-stacked layers, indicating a diffusive transport regime with ballistic contributions. These results are supported by the theoretical analyses of the accumulated thermal conductivity, which show an important contribution of phonons with mean free paths between 10 and 200 nm. Moreover, our experiments show a reduction of the $\kappa_\perp$ in the low-temperature ferromagnetic phase occurring at the magnetic transition. The calculations show that this reduction in $\kappa_\perp$ is associated with a decrease in the group velocities of the acoustic phonons and an increase in the phonon-phonon scattering of the Raman modes that couple to the magnetic phase. These results demonstrate the potential of van der Waals ferromagnets for thermal transport engineering.

Kagome Materials with flat bands exhibit wildly different physical properties depending on symmetry group, and electron number. For the case of ScV$_6$Sn$_6$ in space group 191, we investigate the existence of a charge density wave (CDW) at vector $\bar{K}=(\frac{1}{3},\frac{1}{3},\frac{1}{3})$ and its relationship with the phonon behavior. The experimental findings reveal a $\sim$95K CDW without nesting/peaks in the electron susceptibility at $\bar{K}$. Notably, ScV$_6$Sn$_6$ exhibits a collapsed phonon mode at $H=(\frac{1}{3},\frac{1}{3},\frac{1}{2})$ and an imaginary flat phonon band in the vicinity of $H$. The soft phonon is attributed to triangular Sn ($Sn^T$) mirror-even vibrations along the $z$-direction. We develop a simple force constant model to describe the entire soft phonon dispersion. By employing a new (Gaussian) approximation of the hopping parameter, we demonstrate the renormalization of the phonon frequency and the consequent collapse of the $H$ phonon. Additionally, we propose an effective model with two order parameters (OPs) to explain the appearance of the CDW at $\bar{K}$, which competes with the collapsed phonon at $H$. Through comparisons with experimental data, we show that the $H$ OP undergoes a second-order phase transition while exhibiting substantial fluctuations, ultimately inducing the first-order transition of the $\bar{K}$ OP. Furthermore, we extend our analysis to the similar compound YV$_6$Sn$_6$, which lacks a CDW phase, attributing this difference to the participation of the heavier Y atom in the out-of-plane phonon behavior. Our comprehensive study not only elucidates the CDW in ScV$_6$Sn$_6$ but also presents a significant advancement in modeling complex electronic systems, fostering collaborations between ab-initio simulations and analytical approaches., Comment: 5+108 pages, 3+42 figures, previously submitted. See also the related experimental study arXiv:2304.09173

The long range electronic modulations recently discovered in the geometrically frustrated kagome lattice have opened new avenues to explore the effect of correlations in materials with topological electron flat bands. The observation of the lattice response to the emergent new phases of matter, a soft phonon mode, has remained elusive and the microscopic origin of charge density waves (CDWs) is still unknown. Here, we show, for the first time, a complete melting of the ScV$_ 6$Sn$_ 6$ (166) kagome lattice. The low energy phonon with propagation vector $\frac{1}{3} \frac{1}{3} \frac{1}{2}$ collapses at 98 K, without the emergence of long-range charge order, which sets in with a propagation vector $\frac{1}{3} \frac{1}{3} \frac{1}{3}$. The CDW is driven (but locks at a different vector) by the softening of an overdamped phonon flat plane at k$_z$=$\pi$. We observe broad phonon anomalies in momentum space, pointing to (1) the existence of approximately flat phonon bands which gain some dispersion due to electron renormalization, and (2) the effects of the momentum dependent electron-phonon interaction in the CDW formation. Ab initio and analytical calculations corroborate the experimental findings to indicate that the weak leading order phonon instability is located at the wave vector $\frac{1}{3} \frac{1}{3} \frac{1}{2}$ of a rather flat collapsed mode. We analytically compute the phonon frequency renormalization from high temperatures to the soft mode, and relate it to a peak in the orbital-resolved susceptibility, obtaining an excellent match with both ab initio and experimental results, and explaining the origin of the approximately flat phonon dispersion. Our data report the first example of the collapse of a softening of a flat phonon plane and promote the 166 compounds of the kagome family as primary candidates to explore correlated flat phonon-topological flat electron physics., Comment: 10 pages, 4 figures

In the study of dynamic charge order correlations in the cuprates, most high energy-resolution resonant inelastic x-ray scattering (RIXS) measurements have focused on momenta along the high-symmetry directions of the copper oxide plane. However, electron scattering along other in-plane directions should not be neglected as they may contain information relevant, for example, to the origin of charge order correlations or to our understanding of the isotropic scattering responsible for strange metal behavior in cuprates. We report high-resolution resonant inelastic x-ray scattering (RIXS) experiments that reveal the presence of dynamic electron correlations over the $q_x$-$q_y$ scattering plane in underdoped $\textrm{Bi}_2\textrm{Sr}_2\textrm{Ca}\textrm{Cu}_2\textrm{O}_{8+\delta}$ with $T_c=54$ K. We use the softening of the RIXS-measured bond stretching phonon line as a marker for the presence of charge-order-related dynamic electron correlations. The experiments show that these dynamic correlations exist at energies below approximately $70$ meV and are centered around a quasi-circular manifold in the $q_x$-$q_y$ scattering plane with radius equal to the magnitude of the charge order wave vector, $q_{CO}$. We also demonstrate how this phonon-tracking procedure provides the necessary experimental precision to rule out fluctuations of short-range directional charge order (i.e. centered around $[q_x=\pm q_{CO}, q_y=0]$ and $[q_x=0, q_y=\pm q_{CO}]$) as the origin of the observed correlations., Comment: Submitted and under review

We investigate the electronic structure of Cs$_2$CuCl$_4$, a material discussed in the framework of a frustrated quantum antiferromagnet, by means of resonant inelastic x-ray scattering (RIXS) and Density Functional Theory (DFT). From the non-dispersive highly localized dd excitations, we resolve the crystal field splitting of the Cu$^{2+}$ ions in a strongly distorted tetrahedral coordination. This allows us to model the RIXS spectrum within the Crystal Field Theory (CFT), assign the dd orbital excitations and retrieve experimentally the values of the crystal field splitting parameters D$_q$, D$_s$ and D$_{\tau}$. The electronic structure obtained ab-initio agrees with the RIXS spectrum and modelled by CFT, highlighting the potential of combined spectroscopic, cluster and DFT calculations to determine the electronic ground state of complex materials., Comment: 7 pages, 4 figures

VSe$_2$ is a layered compound that has attracted great attention due to its proximity to a ferromagnetic state that is quenched by the presence of a charge density wave (CDW) phase. In the monolayer limit, unrelated experiments have reported different CDW orders with transition temperatures in the range of 130 to 220 K, making this monolayer very controversial. Here we perform first-principles non-perturbative anharmonic phonon calculations in monolayer VSe$_2$ in order to estimate the CDW order and the corresponding transition temperature. Our analysis solves previous experimental contradictions as we reveal that monolayer VSe$_2$ develops two independent charge density wave orders associated to $\sqrt{3} \times \sqrt{7}$ and $4 \times 4$ modulations that compete as a function of strain. In fact, tiny variations of only 1.5% in the lattice parameter are enough to stabilize one order or the other, which underlines that the CDW order becomes substrate-dependent. The predicted CDW temperature is strain-dependent and has a value of around 220 K, in good agreement with experiments. Moreover, we analyze the impact of external Lennard-Jones interactions on the CDW. We show that these can act together with the anharmonicity to suppress the CDW orders. In the particular case of monolayer VSe$_2$, this may give rise to the emergence of a ferromagnetic order., Comment: 8 pages, 4 figures, and supplemental material

The nature of the charge density wave phases in the kagome metal compound AV$_3$Sb$_5$ has raised many questions and their origin is still under debate. Here, we combine diffuse scattering (DS) and inelastic x-ray scattering (IXS) to identify a 3-dimensional (3D) precursor of the charge order with propagation vector (0.5 0 0.5), which condenses into a 3D-CDW through a first order phase transition. The quasi-elastic critical scattering indicates that the dominant contribution to the diffuse precursor is the elastic central peak without phonon softening. However, the inelastic spectra show a small broadening of the Einstein-type phonon mode on approaching the phase transition. Our results point to the situation where the Fermi surface instability at the \textit{L} point is of order-disorder type with a critical growth of quasi-static domains. The results go beyond the classical weak-coupling Peierls transition dynamics and are discussed within models including strong-electron phonon coupling and non-adiabaticity.

We analyze the application of the history state formalism to quantum walks. The formalism allows one to describe the whole walk through a pure quantum history state, which can be derived from a timeless eigenvalue equation. It naturally leads to the notion of system-time entanglement of the walk, which can be considered as a measure of the number of orthogonal states visited in the walk. We then focus on one-dimensional discrete quantum walks, where it is shown that such entanglement is independent of the initial spin orientation for real Hadamard-type coin operators and real initial states (in the standard basis) with definite site parity. Moreover, in the case of an initially localized particle it can be identified with the entanglement of the unitary global operator that generates the whole history state, which is related to its entangling power and can be analytically evaluated. Besides, it is shown that the evolution of the spin subsystem can also be described through a spin history state with an extended clock. A connection between its average entanglement (over all initial states) and that of the operator generating this state is also derived. A quantum circuit for generating the quantum walk history state is provided as well., Comment: 12 pages, 7 figures, final version