Search

Kagome superconductors AV$_{3}$Sb$_{5}$ provide a unique platform for studying the interplay between a variety of electronic orders, including superconductivity, charge density waves, nematic phases and more. Understanding the evolution of the electronic state from the charge density wave to the superconducting transition is essential for unraveling the interplay of charge, spin, and lattice degrees of freedom giving rise to the unusual magnetic properties of these nonmagnetic metals. Previous zero-field and high-field $\mu$SR studies revealed two anomalies in the muon spin relaxation rate, a first change at $T_{CDW} \sim 100$ K and a second steep increase at $T^{*}\sim 40$ K, further enhanced by an applied magnetic field, thus suggesting a contribution of magnetic origin. In this study, we use the avoided level crossing $\mu$SR technique to investigate charge order in near-zero applied field. By tracking the temperature dependence of quadrupolar level-crossing resonances, we examined the evolution of the electric field gradient at V nuclei in the kagome plane. Our results show a significant rearrangement of the charge density starting at $T^{*}$ indicating a transition in the charge distribution, likely electronic in origin, well below $T_{CDW}$. These findings, combined with previous $\mu$SR, STM, and NMR studies, emphasize the intertwined nature of proximate phases in these systems, with the charge rearrangement dominating the additional increase in $\mu$SR relaxation rate below $T^{*}$.

The discovery of chiral helical magnetism (CHM) in Cr$_{1/3}$NbS$_2$ and the stabilization of a chiral soliton lattice (CSL) has attracted considerable interest in view of their potential technological applications. However, there is an ongoing debate regarding whether the sister compound, Mn$_{1/3}$NbS$_2$, which shares the same crystal structure, exhibits similar nontrivial properties which rely on the stabilization of the lack of inversion symmetry at the magnetic ion. In this study, we conduct a comprehensive investigation of the magnetically ordered states of both compounds, using $^{53}$Cr, $^{55}$Mn and $^{93}$Nb nuclear magnetic resonance. Our results, supported by density functional calculations, detect in a high-quality single crystal of Cr$_{1/3}$NbS$_2$ all the signatures of the monoaxial CHM in a magnetic field, identifying it as a textbook NMR case. The detailed understanding of this prototypic behavior provides a reference for Mn$_{1/3}$NbS$_2$. Despite the much larger density of specific defects in this second single crystal, we confirm the presence of a CHM phase in the Mn compound, characterized by a very large critical field for the forced ferromagnetic phase ($\approx 5$ T for H$\parallel\hat c$)., Comment: 10 pages, 7 figures, Supplemental MAterial

Positive muon spin rotation and relaxation spectroscopy is a well established experimental technique for studying materials. It provides a local probe that generally complements scattering techniques in the study of magnetic systems and represents a valuable alternative for materials that display strong incoherent scattering or neutron absorption. Computational methods can effectively quantify the microscopic interactions underlying the experimentally observed signal, thus substantially boosting the predictive power of this technique. Here, we present an efficient set of algorithms and workflows devoted to the automation of this task. In particular, we adopt the so-called DFT+{\mu} procedure, where the system is characterised in the density functional theory (DFT) framework with the muon modeled as a hydrogen impurity. We devise an automated strategy to obtain candidate muon stopping sites, their dipolar interaction with the nuclei, and hyperfine interactions with the electronic ground state. We validate the implementation on well-studied compounds, showing the effectiveness of our protocol in terms of accuracy and simplicity of use

A comprehensive set of muon spin spectroscopy and neutron scattering measurements supported by ab-initio and model Hamiltonian simulations have been used to investigate the magnetic ground state of Na$_2$PrO$_3$. $\mu$SR reveals N\'eel antiferromagnetic order below $T_{\rm N}\! \sim\! 4.9$ K, with a small static magnetic moment $m_{\rm static}\!\leq \! 0.22$~$\mu_{\rm B}/{\rm Pr}$ collinearly aligned along the $c-$axis. Inelastic neutron measurements reveal the full spectrum of crystal field excitations and confirm that the Pr$^{4+}$ ground state wave function deviates significantly from the $\Gamma_7$ limit relevant to the Kitaev model. Single and two magnon excitations are observed in the ordered state below $T_N=4.6$ K and are well described by non-linear spin wave theory from the N\'eel state using a magnetic Hamiltonian with Heisenberg exchange $J=1$ meV and symmetric anisotropic exchange $\Gamma/J=0.1$, corresponding to an XY model. Intense two magnon excitations are accounted for by $g$-factor anisotropy $g_\mathrm{z}/g_\pm = 1.29$. A fluctuating moment $\delta m^2 = 0.57(22)$ $\mu_{\rm B}^2/{\rm Pr}$ extracted from the energy and momentum integrated inelastic neutron signal is reduced from expectations for a local $J=1/2$ moment with average $g$-factor $g_{\rm avg}\approx 1.1$. Together, the results demonstrate that the small moment in Na$_2$PrO$_3$ arises from crystal field and covalency effects and that the material does not exhibit significant quantum fluctuations..

The breaking of time-reversal symmetry (TRS) in the normal state of kagome superconductors $A$V$_{3}$Sb$_{5}$ stands out as a significant feature. Yet the extent to which this effect can be tuned remains uncertain, a crucial aspect to grasp in light of the varying details of TRS breaking observed through different techniques. Here, we employ the unique low-energy muon spin rotation technique combined with local field numerical analysis to study the TRS breaking response as a function of depth from the surface in single crystals of RbV$_{3}$Sb$_{5}$ with charge order and Cs(V$_{0.86}$Ta$_{0.14}$)$_{3}$Sb$_{5}$ without charge order. In the bulk (i.e., > 33 nm from the surface) of RbV$_{3}$Sb$_{5}$, we have detected a notable increase in the internal magnetic field width experienced by the muon ensemble. This increase occurs only within the charge ordered state. Intriguingly, the muon spin relaxation rate is significantly enhanced near the surface (i.e., < 33 nm from the surface) of RbV$_{3}$Sb$_{5}$, and this effect commences at temperatures significantly higher than the onset of charge order. Conversely, in Cs(V$_{0.86}$Ta$_{0.14}$)$_{3}$Sb$_{5}$, we do not observe a similar enhancement in the internal field width, neither in the bulk nor near the surface. These observations indicate a strong connection between charge order and TRS breaking on one hand, and on the other hand, suggest that TRS breaking can occur prior to long-range charge order. This research offers compelling evidence for depth-dependent magnetism in $A$V$_{3}$Sb$_{5}$ superconductors in the presence of charge order. Such findings are likely to elucidate the intricate microscopic mechanisms that underpin the TRS breaking phenomena in these materials., Comment: 16 pages, 13 figures

Magnetostriction drives a rhombohedral distortion in the cubic rock salt antiferromagnet MnO at the N\'eel temperature $T_{N}=118$ K. As an unexpected consequence we show that this distortion acts to localize the site of an implanted muon due to the accompanying redistribution of electron density. This lifts the degeneracy between equivalent sites, resulting in a single observed muon precession frequency. Above $T_{N}$, the muon instead becomes delocalized around a network of equivalent sites. Our first-principles simulations based on Hubbard-corrected density-functional theory and molecular dynamics are consistent with our experimental data and help to resolve a long-standing puzzle regarding muon data on MnO, as well as having wider applicability to other magnetic oxides., Comment: 7 pages, 3 figures

Mejia-Rodriguez and Trickey recently proposed a procedure for removing the explicit dependence of meta-GGA exchange-correlation energy functionals $E_{\rm xc}$ on the kinetic energy density $\tau$. We present a simple modification to this approach in which the exact Kohn-Sham $\tau$ is used as input for $E_{\rm xc}$ but the functional derivative of $\tau$ with respect to the density $\rho$, required to calculate the potential term $\int d^3r'\,\delta E_{\rm xc}/\delta\tau({\bf r}')|_{\rho}\cdot \delta\tau({\bf r}')/\delta\rho({\bf r})$, is evaluated using an approximate kinetic energy density functional. This ensures that the Kohn-Sham potential is a local multiplicative function as opposed to the non-local potential of a generalized Kohn-Sham approach. Electronic structure codes can be easily modified to use the new method. We validate it by quantifying the accuracy of the predicted lattice parameters, bulk moduli, magnetic moments and cohesive energies of a large set of periodic solids. An unanticipated benefit of this method is to gauge the quality of approximate kinetic energy functionals by checking if the self-consistent solution is indeed at the variational minimum.

Heusler compounds have been intensively studied owing to the important technological advancements that they provide in the field of shape memory, thermomagnetic energy conversion and spintronics. Many of their intriguing properties are ultimately governed by their magnetic states and understanding and possibly tuning them is evidently of utmost importance. In this work we examine the \alloys alloys with Density Functional Theory simulations and $^{55}$Mn Nuclear Magnetic Resonance and combine these two methods to carefully describe their ground state magnetic order. In addition, we compare the results obtained with the conventional generalized gradient approximation with the ones of strongly constrained and appropriately normed (SCAN) semilocal functionals for exchange and correlation. Experimental results eventually allow to discriminate between two different scenarios identified by ab initio simulations.

The recently discovered vanadium-based kagome metals AV$_3$Sb$_5$ (A = K, Rb, Cs) offer the possibility to study the interplay between competing electronic orderings, such as charge density order and superconductivity. We focus on the former and provide a comprehensive set of $^{51}$V, $^{87}$Rb, and $^{121}$Sb magnetic resonance measurements on an RbV$_3$Sb$_5$ single crystal. Elucidating the symmetries and properties of the CDW phase is essential to understanding the unconventional electronic orderings occurring in this material. We establish the structure of the $2\times 2 \times 2$ superlattice that describes the system below the charge density wave transition by combining both experimental and computational methods, with a methodology that can be readily applied to the remaining compounds of the same family. Our results give compelling evidence that the CDW structure occurring below 103 K for RbV$_3$Sb$_5$ is the so-called Inverse Start of David pattern $\pi$-shifted along the c axis (also known as staggered tri-hexagonal).

Magnetic frustration in ferromagnetic metallic systems is unusual due to the long-range and symmetric nature of the exchange interactions. In this work we prove that it is possible to obtain a highly frustrated ferromagnetic phase in a multi-sublattices cubic structure through a fine tuning of the magnetic interactions. This peculiar state is achieved in Ni-Mn-(In, Sn) Heusler alloys and results in the effective decoupling of their two intertwined ferromagnetic sublattices. One sublattice is ferromagnetic long range ordered below the macroscopic Curie temperature (TC ) whereas the second one remains disordered until a crossover to a polarized state occurs at T << TC . This result points out that a fine engineering of the magnetic interactions in metallic systems can lead to interesting novel and emergent phenomena.

We report on the first example of quantum coherence between the spins of muons and quadrupolar nuclei. We observe this effect in vanadium intermetallic compounds which adopt the A15 crystal structure, and whose members include all technologically dominant superconductors. The entangled states are extremely sensitive to the local structural and electronic environments through the electric field gradient at the quadrupolar nuclei. This case-study demonstrates that positive muons can be used as a quantum sensing tool to probe also structural and charge related phenomena in materials, even in the absence of magnetic order.

We report muon spin rotation ($\mu$SR) and neutron diffraction on the rare-earth based magnets (Mo$_{2/3}$RE$_{1/3}$)$_2$AlC, also predicted as parent materials for 2D derivatives, where RE = Nd, Gd (only ($\mu$SR), Tb, Dy, Ho and Er. By crossing information between the two techniques, we determine the magnetic moment ($m$), structure, and dynamic properties of all compounds. We find that only for RE = Nd and Gd the moments are frozen on a microsecond time scale. Out of these two, the most promising compound for a potential 2D high ($m$) magnet is the Gd variant, since the parent crystals are pristine with $m = 6.5 \pm 0.5 \mu_B$, N\'eel temperature of $29 \pm 1$ K, and the magnetic anisotropy between in and out of plane coupling is smaller than $10^{-8}$. This result suggests that magnetic ordering in the Gd variant is dominated by in-plane magnetic interactions and should therefore remain stable if exfoliated into 2D sheets., Comment: Includes Supplementary Information

We present a systematic investigation of muon-stopping states in superconductors that reportedly exhibit spontaneous magnetic fields below their transition temperatures due to time-reversal symmetry breaking. These materials include elemental rhenium, several intermetallic systems and Sr$_2$RuO$_4$. We demonstrate that the presence of the muon leads to only a limited and relatively localized perturbation to the local crystal structure, while any small changes to the electronic structure occur several electron volts below the Fermi energy leading to only minimal changes in the charge density on ions close to the muon. Our results imply that the muon-induced perturbation alone is unlikely to lead to the observed spontaneous fields in these materials, whose origin is more likely intrinsic to the time-reversal symmetry broken superconducting state., Comment: Main text: 6 pages, 3 figures. Supplemental Material: 13 pages, 16 figures

We present a detailed study of the electronic and magnetic interactions of Ba2CoO4, structurally very uncommon because of the isolated CoO4 distorted tetrahedral coordination. We show the presence of Co(d)-O(p) hybridized states characterized by spin polarized oxygen atoms, with their magnetic moments parallel to that on Co. The calculated isotropic exchange interaction parameters, which include the contributions from ligand spins, demonstrate the presence of a 3D network of magnetic couplings, that are partially frustrated in the identified magnetic ground state. Our results indicate that the dominant indirect exchange mechanism responsible for this ground state is mediated by O atoms along the Co-O...O-Co path., Comment: 5 pages, 3 figures, 3 tables and supplemental information

Quantum ESPRESSO is an open-source distribution of computer codes for quantum-mechanical materials modeling, based on density-functional theory, pseudopotentials, and plane waves, and renowned for its performance on a wide range of hardware architectures, from laptops to massively parallel computers, as well as for the breadth of its applications. In this paper we present a motivation and brief review of the ongoing effort to port Quantum ESPRESSO onto heterogeneous architectures based on hardware accelerators, which will overcome the energy constraints that are currently hindering the way towards exascale computing.

We report on the electronic properties of an artificial system obtained by the intercalation of equiatomic FeCo layers under graphene grown on Ir(111). Upon intercalation, the FeCo film grows epitaxially on Ir(111), resulting in a lattice-mismatched system. By performing Density Functional Theory calculations, we show that the intercalated FeCo layer leads to a pronounced corrugation of the graphene film. At the same time, the FeCo intercalated layers induce a clear transition from a nearly undisturbed to a strongly hybridized graphene {\pi}-band, as measured by angle-resolved photoemission spectroscopy. A comparison of experimental results with the computed band structure and the projected density of states unveils a spin-selective hybridization between the {\pi} band of graphene and FeCo-3d states. Our results demonstrate that the reduced dimensionality, as well as the hybridization within the FeCo layers, induce a narrowing and a clear splitting of Fe 3d-up and Fe 3d-down spin bands of the confined FeCo layers with respect to bulk Fe and Co., Comment: 5 pages, 4 figures

Fe$_2$P alloys have been identified as promising candidates for magnetic refrigeration at room-temperature and for custom magnetostatic applications. The intent of this study is to accurately characterize the magnetic ground state of the parent compound, Fe$_2$P, with two spectroscopic techniques, $\mu$SR and NMR, in order to provide solid bases for further experimental analysis of Fe$_2$P-type transition metal based alloys. We perform zero applied field measurements using both techniques below the ferromagnetic transition $T_C=220~\mathrm K$. The experimental results are reproduced and interpreted using first principles simulations validating this approach for quantitative estimates in alloys of interest for technological applications., Comment: 10 pages, 2 figures

In the heavy-fermion system $Yb_2Pd_2In_{1-x}Sn_x$, the interplay of crystal-field splitting, Kondo effect, and Ruderman-Kittel-Kasuya-Yosida interactions leads to complex chemical-, pressure-, and magnetic-field phase diagrams, still to be explored in full detail. By using a series of techniques, we show that even modest changes of parameters other than temperature are sufficient to induce multiple quantum-critical transitions in this highly susceptible heavy-fermion family. In particular, we show that, above $\sim 10$ kbar, hydrostatic pressure not only induces an antiferromagnetic phase at low temperature, but it likely leads to a reorientation of the Yb magnetic moments and/or the competition among different antiferromagnetic configurations.

In muon spin rotation experiments the positive implanted muon vibrates with large zero point amplitude by virtue of its light mass. Quantum mechanical calculations of the host material usually treat the muon as a point impurity, ignoring this large vibrational amplitude. As a first order correction, the muon zero point motion is usually described within the harmonic approximation, despite the large anharmonicity of the crystal potential. Here we apply the stochastic self-consistent harmonic approximation, a quantum variational method devised to include strong anharmonic effects in total energy and vibrational frequency calculations, in order to overcome these limitations and provide an accurate ab initio description of the quantum nature of the muon. We applied this full quantum treatment to the calculation of the muon contact hyperfine field in textbook-case metallic systems, such as Fe, Ni, Co including MnSi and MnGe, significantly improving agreement with experiments. Our results show that muon vibrational frequencies are strongly renormalized by anharmonicity. Finally, in contrast to the harmonic approximation, we show that including quantum anharmonic fluctuations, the muon stabilizes at the octahedral site in bcc Fe., Comment: 10 pages

yambo is an open source project aimed at studying excited state properties of condensed matter systems from first principles using many-body methods. As input, yambo requires ground state electronic structure data as computed by density functional theory codes such as quantum-espresso and abinit. yambo's capabilities include the calculation of linear response quantities (both independent-particle and including electron-hole interactions), quasi-particle corrections based on the GW formalism, optical absorption, and other spectroscopic quantities. Here we describe recent developments ranging from the inclusion of important but oft-neglected physical effects such as electron-phonon interactions to the implementation of a real-time propagation scheme for simulating linear and non-linear optical properties. Improvements to numerical algorithms and the user interface are outlined. Particular emphasis is given to the new and efficient parallel structure that makes it possible to exploit modern high performance computing architectures. Finally, we demonstrate the possibility to automate workflows by interfacing with the yambopy and AiiDA software tools.