Search

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..

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).

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 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 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

Power and energy consumption is becoming key challenges to deploy the first exascale supercomputer successfully. Large-scale HPC applications waste a significant amount of power in communication and synchronization-related idle times. However, due to the time scale at which communication happens, transitioning in low power states during communication's idle times may introduce unacceptable overhead in applications' execution time. In this paper, we present COUNTDOWN, a runtime library, supported by a methodology and analysis tool for identifying and automatically reducing the power consumption of the computing elements during communication and synchronization. COUNTDOWN saves energy without imposing significant time-to-completion increase by lowering CPUs power consumption only during idle times for which power state transition overhead are negligible. This is done transparently to the user, without requiring labor-intensive and error-prone application code modifications, nor requiring recompilation of the application. We test our methodology in a production Tier-0 system. For the NAS benchmarks, COUNTDOWN saves between 6% and 50% energy, with a time-to-solution penalty lower than 5%. In a complete production --- Quantum ESPRESSO --- for a 3.5K cores run, COUNTDOWN saves 22.36% energy, with a performance penalty below 3%. Energy saving increases to 37% with a performance penalty of 6.38%, if the application is executed without communication tuning., Comment: 14 pages, 11 figures

In positive muon spin rotation and relaxation spectroscopy it is becoming nowadays customary to take advantage of Density Functional Theory (DFT) based computational methods to aid the experimental data analysis. DFT aided muon site determination is especially useful for measurements performed in magnetic materials, where large contact hyperfine interactions may arise. Here we present a systematic analysis of the accuracy of the ab initio estimation of muon's hyperfine contact field on elemental transition metals, performing state of the art spin-polarized plane wave DFT and using the projector augmented pseudopotential approach, which allows to include the core state effects due to the spin ordering. We further validate this method in not-so-simple, non-centrosymmetric metallic compounds, presently of topical interest for their spiral magnetic structure giving rise to skyrmion phases, such as MnSi and MnGe. The calculated hyperfine fields agree with experimental values in all cases, provided the spontaneous spin magnetization of the metal is well reproduced within the approach. To overcome the known limits of the conventional mean field approximation of DFT on itinerant magnets, we adopt the so-called reduced Stoner theory [L. Ortenzi et al.,Phys. Rev. B 86, 064437 (2012)]. We establish the accuracy of the estimated muon contact field in metallic compounds with DFT and our results show improved agreement with experiments compared to those of earlier publications., Comment: 8 pages, 4 figures

Muon spin rotation and relaxation studies have been performed on a "111" family of iron-based superconductors NaFe_1-xNi_xAs. Static magnetic order was characterized by obtaining the temperature and doping dependences of the local ordered magnetic moment size and the volume fraction of the magnetically ordered regions. For x = 0 and 0.4 %, a transition to a nearly-homogeneous long range magnetically ordered state is observed, while for higher x than 0.4 % magnetic order becomes more disordered and is completely suppressed for x = 1.5 %. The magnetic volume fraction continuously decreases with increasing x. The combination of magnetic and superconducting volumes implies that a spatially-overlapping coexistence of magnetism and superconductivity spans a large region of the T-x phase diagram for NaFe_1-xNi_xAs . A strong reduction of both the ordered moment size and the volume fraction is observed below the superconducting T_C for x = 0.6, 1.0, and 1.3 %, in contrast to other iron pnictides in which one of these two parameters exhibits a reduction below TC, but not both. The suppression of magnetic order is further enhanced with increased Ni doping, leading to a reentrant non-magnetic state below T_C for x = 1.3 %. The reentrant behavior indicates an interplay between antiferromagnetism and superconductivity involving competition for the same electrons. These observations are consistent with the sign-changing s-wave superconducting state, which is expected to appear on the verge of microscopic coexistence and phase separation with magnetism. We also present a universal linear relationship between the local ordered moment size and the antiferromagnetic ordering temperature TN across a variety of iron-based superconductors. We argue that this linear relationship is consistent with an itinerant-electron approach, in which Fermi surface nesting drives antiferromagnetic ordering., Comment: 20 pages, 14 figures, Correspondence should be addressed to Prof. Yasutomo Uemura: tomo@phys.columbia.edu

The temperature-pressure phase diagram of the ferromagnet LaCrGe$_3$ is determined for the first time from a combination of magnetization, muon-spin-rotation and electrical resistivity measurements. The ferromagnetic phase is suppressed near $2.1$~GPa, but quantum criticality is avoided by the appearance of a magnetic phase, likely modulated, AFM$_Q$. Our density functional theory total energy calculations suggest a near degeneracy of antiferromagnetic states with small magnetic wave vectors $Q$ allowing for the potential of an ordering wave vector evolving from $Q=0$ to finite $Q$, as expected from the most recent theories on ferromagnetic quantum criticality. Our findings show that LaCrGe$_3$ is a very simple example to study this scenario of avoided ferromagnetic quantum criticality and will inspire further study on this material and other itinerant ferromagnets.

The rapid developments of computational quantum chemistry methods and supercomputing facilities motivate the renewed interest in the analysis of the muon/electron interactions in $\mu$SR experiments with \emph{ab initio} approaches. Modern simulation methods seem to be able to provide the answers to the frequently asked questions of many $\mu$SR experiments: where is the muon? Is it a passive probe? What are the interaction parameters governing the muon-sample interaction? In this review we describe some of the approaches used to provide quantitative estimations of the aforementioned quantities and we provide the reader with a short discussion on the current developments in this field.

The discovery of a superconducting dome in the proximity of the magnetic to paramagnetic transition in the electronic phase diagram of MnP as a function of hydrostatic pressure has renewed the interest in the magnetic properties of this binary pnictide. Here we present an ab initio study obtained with Density Functional Theory (DFT) simulations as a function of pressure. We show that the itinerant-electron magnetism of MnP is well characterized by the mean-field Kohn-Sham method which correctly describes the ambient pressure magnetic states, the anomalous trend of the lattice parameters as a function of pressure and the critical pressure for the disappearance of the magnetic order. We are finally able to confirm the nature of the new helical structure observed at high pressure.

Neutron diffraction and muon spin relaxation measurements are used to obtain a detailed phase diagram of Pr(Fe,Ru)AsO. The isoelectronic substitution of Ru for Fe acts effectively as spin dilution, suppressing both the structural and magnetic phase transitions. The temperature of the tetragonal-orthorhombic structural phase transition decreases gradually as a function of x. Slightly below the transition temperature coherent precessions of the muon spin are observed corresponding to static magnetism, possibly reflecting a significant magneto-elastic coupling in the FeAs layers. Short range order in both the Fe and Pr moments persists for higher levels of x. The static magnetic moments disappear at a concentration coincident with that expected for percolation of the J1-J2 square lattice model.

We investigate the effect of external pressure on magnetic order in undoped LnFeAsO (Ln = La, Ce, Pr, La) by using muon-spin relaxation measurements and ab-initio calculations. Both magnetic transition temperature $T_m$ and Fe magnetic moment decrease with external pressure. The effect is observed to be lanthanide dependent with the strongest response for Ln = La and the weakest for Ln = Sm. The trend is qualitatively in agreement with our DFT calculations. The same calculations allow us to assign a value of 0.68(2) $\mu_B$ to the Fe moment, obtained from an accurate determination of the muon sites. Our data further show that the magnetic lanthanide order transitions do not follow the simple trend of Fe, possibly as a consequence of the different $f$-electron overlap., Comment: 16 pages, 11 figures

The effect of isoelectronic substitutions on the microscopic properties of LaFe1-xRuxAsO, for 0< x< 0.8, has been investigated by means of muSR and 139La NMR. It was found that Ru substitution causes a progressive reduction of the N\`eel temperature (T_N) and of the magnetic order parameter without leading to the onset of superconductivity. The temperature dependence of 139La nuclear spin-lattice relaxation rate 1/T_1 can be suitably described within a two-band model. One band giving rise to the spin density wave ground-state, while the other one is characterized by weakly correlated electrons. Fe for Ru substitution yields to a progressive decrease of the density of states at the Fermi level close to the one derived from band structure calculations. The reduction of T_N with doping follows the predictions of the J_1-J_2 model on a square lattice, which appears to be an effective framework to describe the magnetic properties of the spin density wave ground-state., Comment: 6 pages, 8 figures