Search

The kagome lattice, known for its strong frustration in two dimensions, hosts a variety of exotic magnetic and electronic states. A variation of this geometry, where the triangular motifs are twisted to further reduce symmetry, has recently revealed even more complex physics. HoAgGe exemplifies such a structure, with magnetic and electronic properties believed to be driven by strong in-plane anisotropy of the Ho spins, effectively acting as a two-dimensional spin ice. In this study, using a combination of magnetization, Hall conductivity measurements, and density functional theory calculations, we demonstrate how various spin-ice states, stabilized by external magnetic fields, influence the Fermi surface topology. More interestingly, we observe sharp transitions in Hall conductivity without concurrent changes in magnetization when an external magnetic field is applied along a particular crystallographic direction, underscoring the role of strong magnetic frustration and providing a new platform for exploring the interplay between magnetic frustration, electronic topology, and crystalline symmetry. These results also highlight the limitations of a simple spin-ice model, suggesting that a more sophisticated framework is necessary to capture the subtle experimental nuances observed., Comment: 17 pages, 11 figures

Mn-based transition metal Zintl compounds in the 14-1-11 phase are known to host complex atomic and magnetic structures owing to their intricate crystal structure. Among them, Ba14MnBi11 stands out as one of the least understood compounds, with experimental measurements and theoretical findings largely inconsistent. Following up on the earlier attempt [D. Sanchez-Portal et al., PRB 65, 144414 (2002)] at establishing a connection between metallicity and magnetism through a DFT-based analysis, our work aims to provide additional insights to resolve the existing contradictions. Our key finding is that the magnetic ground state is very susceptible to charge doping. DFT calculations for stoichiometric Ba14MnBi11 give a rather stable ferromagnetic metallic ground state. However, by adding exactly one additional electron per Mn, the system becomes semiconducting and the magnetic ground state becomes weakly antiferromagnetic (AF). On the other hand, upon small hole doping the system transitions to a special type of AF state known as altermagnetic ordering. The observed trends suggest that hole and electron doping-induced phase transitions likely result from different underlying mechanisms, influencing various exchange pathways. Additionally, our projected density-of-states along with bandstructure analyses indicate that, besides the largest hole contribution coming from the tetrahedral unit of Bi, the isolated Bi sites also play a substantial role and the dispersive bands near VBM suggest a rather complex hybridization network involving both Bi band characters. Through a comprehensive comparison of available data and our analysis, we propose that the inconsistency in magnetic states between experimental findings and DFT calculations is due to nonstoichiometric effects, likely impurities or defects., Comment: 9 pages, 9 figures

The classical Ising model on the triangular lattice (we will call it I-3 model below), while simple in the nearest-neighbors (NN) only approximation, becomes increasingly richer and more complex when further interactions are incorporated. However, the studies so far have not been exhaustive, nor have any attempts been made to estimate how realistic are the parameter ranges that generate strong metamagnetism with a large number of magnetization steps. In this study, we identify one such candidate, ErGa$_{2}$, a material known to have one strong magnetization step, albeit some narrow steps below and above cannot be confidently excluded. It has been established, and we can confirm the same computationally, to have an easy axis perpendicular to the triangular Er plane, with a strong anisotropy and with a large magnetic moment of $9.5\:\mu_{B}$, making it a perfect implementation of the classical I-3 model. In the first part of the analysis, we present the I-3 model with up to the third nearest-neighbors in a range of parameters $J_{2}$ and $J_{3}$ ($J_{1}$ in this part is set to 1), and in some cases adding a rather small $J_{4}$ in order to reveal new phases otherwise degenerate with some others. The richest phase diagram is, not surprisingly, observed when all interactions are antiferromagnetic (AF). Subsequently, a more realistic case, inspired by RKKY and by our calculations for ErGa$_{2}$, where $J_{1},J_{2}>0$ (antiferromagnetic) and $J_{3},J_{4}<0$ (ferromagnetic), is presented. Finally, we report our first-principles calculations of $J_{1-4}$ in ErGa$_{2}$ and compared the phase diagram in the regime corresponding to the calculated values with the experiment., Comment: 7 pages, 5 figures

Recently, the pressure has been used as external stimuli to induce structural and magnetic phase transitions in many layered quantum materials whose layers are linked by van der Waals forces. Such materials with weakly held layers allow relatively easy manipulation of the superexchange mechanism and lead to novel magnetic behavior. Using the hydrostatic pressure as a disorderless means to manipulate the interlayer coupling, we applied pressure on two quasi-2D sister compounds, namely, Cr2Si2Te6 (CST) and Mn3Si2Te6 (MST), up to ~1 GPa. Magnetic property measurements with the application of pressure revealed that the ferromagnetic transition temperature decreases in CST while the opposite trend occurs for the ferrimagnet MST. In MST, the magnetization decreases with the increase in the pressure, and such trend is not clearly noticed in CST, within the pressure range studied. Theoretical calculations showed the overall pressure effect on layer separation, bond angle, and exchange coupling, strongly influencing the change in subsequent magnetic characteristics. Exchange coupling in Mn3Si2Te6 is strongly frustrated and the first nearest neighbor interaction is the most dominant of the components with the strongest pressure dependence. Whereas, in Cr2Si2Te6, the exchange coupling parameters exhibit very little dependence on the pressure. This combined experimental and theoretical work has the potential to expand to other relevant quantum materials.

Recently proposed local self-interaction correction (LSIC) method [Zope, R. R. et al., J. Chem. Phys. 151, 214108 (2019)] is a one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange-correlation and Coulomb energy densities. The LSIC method is exact for the one-electron densities, also recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to the well-known Perdew-Zunger SIC (PZSIC) method as a special case. This article presents the self-consistent implementation of the LSIC method using the ratio of Weizs\"acker and Kohn-Sham kinetic energy densities as an iso-orbital indicator. The atomic forces as well as the forces on the Fermi-L\"owdin orbitals are also implemented for the LSIC energy functional. Results show that LSIC with the simplest local spin density functional predicts atomization energies of AE6 dataset better than some of the most widely used GGA functional (e.g. PBE) and barrier heights of BH6 database better than some of the most widely used hybrid functionals (e.g. PBE0 and B3LYP). The LSIC method [mean absolute error (MAE) of 0.008 \r{A}] predicts bond lengths of a small set of molecules better than the PZSIC-LSDA (MAE 0.042 \r{A}) and LSDA (0.011 \r{A}). This work shows that accurate results can be obtained from the simplest density functional by removing the self-interaction-errors using an appropriately designed SIC method., Comment: 34 pages, 4 figures

We review two common numerical schemes for Coulomb potential evaluation that differ only in their radial part of the solutions in the spherical harmonic expansion (SHE). One is based on finite-difference method (FDM) while the other is based on the Green's function (GF) solution to the radial part of the Poisson equation. We analyze the methods and observe that the FDM-based approach appears to be more efficient in terms of the convergence with the number of radial points, particularly for monopole (l=0). However, as a known issue, it suffers from error accumulation as the system size increases. We identify the source of error that comes mainly from l=1 (and sometimes l=2) contribution of SHE induced by the charge partitioning. We then propose a hybrid scheme by combining the two methods, where the radial solution for l=0 is obtained using the FDM method and treating the remaining terms using GF approach. The proposed hybrid method is subsequently applied to a variety of systems to examine its performance. The results show improved accuracy than earlier numerical schemes in all cases. We also show that, even with a generic set of radial grid parameters, accurate energy differences can be obtained using a numerical Coulomb solver in standard density functional studies.

We examine the role of self-interaction errors (SIE) removal on the evaluation of magnetic exchange coupling constants. In particular we analyze the effect of scaling down the self-interaction-correction (SIC) for three {\em non-empirical} density functional approximations (DFAs) namely, the local spin density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation, and recent SCAN family of meta-GGA functionals. To this end, we employ three one-electron SIC methods: Perdew-Zunger [Perdew, J. P.; Zunger, A. \textit{Phys. Rev. B}, {\bf 1981}, \textit{23}, 5048] SIC, the orbitalwise scaled SIC method [Vydrov, O. A. \textit{et al.}, \textit{J. Chem. Phys.} {\bf 2006,} \textit{124}, 094108], and the recent {local} scaling method [Zope, R. R. \textit{et al.}, \textit{J. Chem. Phys.} {\bf 2019}, \textit{151}, 214108]. We compute the magnetic exchange coupling constants using the spin projection and non projection approaches for sets of molecules composed of dinuclear and polynuclear H--He models, organic radical molecules, and chlorocuprate, and compare these results against accurate theories and experiment. Our results show that for the systems that mainly consist of single electron regions, PZSIC performs well but for more complex organic systems and the chlorcuprates, an overcorrecting tendency of PZSIC combined with the DFAs utilized in this work is more pronounced, and in such cases LSIC with kinetic energy density ratio performs better than PZSIC. Analysis of the results in terms of SIC corrections to the density and to the total energy shows that both density and energy correction are required to obtain an improved prediction of magnetic exchange couplings., Comment: 40 pages, 5 figures

The van der Waals class of materials offer an approach to two-dimensional magnetism as their spin fluctuations can be tuned upon exfoliation of layers. Moreover, it has recently been shown that spin-lattice coupling and long-range magnetic ordering can be modified with pressure in van der Waals materials. In this work, the magnetic properties of quasi two-dimensional CrBr3 are reported applying hydrostatic pressure. The application of pressure (0 - 0.8 GPa) shows a 72 % decrease in saturation magnetization with small decrease in the Curie temperature from 33 to 29 K. Density functional theory calculations with pressure up to 1 GPa show a reduction in volume and interplanar distance as pressure increases. To further understand magnetic properties with applied pressure, the magnetocrystalline anisotropy energy (MAE) and exchange coupling parameter (J) are calculated. There is minimal decrease in MAE and the first nearest neighbor interaction (J1) (U = 2.7 eV and J = 0.7 eV), shows an increase in J1 with respect to pressure. Overall, CrBr3 displays ferromagnetic interlayer coupling and the calculated exchange coupling and MAE parameters match well with the observations from the experimental work.

We study realizations of spirals and skyrmions in two-dimensional antiferromagnets with a triangular lattice on an inversion-symmetry-breaking substrate. As a possible material realization, we investigate the adsorption of transition-metal atoms (Cr, Mn, Fe, or Co) on a monolayer of MoS$_2$, WS$_2$, or WSe$_2$ and obtain the exchange, anisotropy, and Dzyaloshinskii-Moriya interaction parameters using first-principles calculations. Using energy minimization and parallel-tempering Monte-Carlo simulations, we determine the magnetic phase diagrams for a wide range of interaction parameters. We find that skyrmion lattices can appear even with weak Dzyaloshinskii-Moriya interactions, but their stability is hindered by magnetic anisotropy. However, a weak easy plane magnetic anisotropy can be beneficial for stabilizing the skyrmion phase. Our results suggest that Cr$/$MoS$_2$, Fe$/$MoS$_2$, and Fe$/$WSe$_2$ interfaces can host spin spirals formed from the 120$^{\circ}$ antiferromagnetic states. Our results further suggests that for other interfaces, such as Fe$/$MoS$_2$, the Dzyaloshinskii-Moriya interaction is strong enough to drive the system into a three-sublattice skyrmion lattice in the presence of experimentally feasible external magnetic field., Comment: 11 pages, 10 figures

A study of the Fermi surface of the putative topological semimetal Pd$_3$Pb has been carried out using Shubnikov-de Haas (SdH) oscillations measured in fields of up to 60 T. Pd$_3$Pb has garnered attention in the community due to a peculiar Fermi surface that has been proposed theoretically by Ahn, Pickett, and Lee, [Phys. Rev. B 98, 035130 (2018)] to host a dispersion-less band along $X-\Gamma$ as well as multiple triply-degenerate band crossings that, under the influence of spin-orbit coupling, lead to ten four-fold degenerate Dirac points. Analysis of the SdH oscillation data verifies the calculated multi-sheet Fermi surface, revealing a $\Gamma$ centered spheroid that had not been resolved experimentally in prior studies. A comprehensive, angle-dependent analysis of the phase of the SdH oscillations convincingly demonstrates a non-trivial Berry phase for two bands along $\Gamma-R$, supporting the theoretical predictions, while simultaneously evidencing interference between extremal orbits that mimics a trivial Berry phase at intermediate angles., Comment: 9 pages, 8 figures (including supplemental document)

We present an {\it ab initio}-based theoretical framework which elucidates the origin of the spin-orbit torque (SOT) in Normal-Metal(NM)/Ferromagnet(FM) heterostructures. The SOT is decomposed into two contributions, namely, {\it spin-Hall} and the {\it spin-orbital} components. We find that {\it (i)} the Field-Like (FL) SOT is dominated by the spin-orbital component and {\it (ii)} both components contribute to the damping-like torque with comparable magnitude in the limit of thick Pt film. The contribution of the spin-orbital component to the DL-SOT is present only for NMs with strong SOC coupling strength. We demonstrate that the FL-SOT can be expressed in terms of the non-equilibrium spin-resolved orbital moment accumulation. The calculations reveal that the experimentally reported oxygen-induced sign-reversal of the FL-SOT in Pt/Co bilayers is due to the significant reduction of the majority-spin orbital moment accumulation on the interfacial NM atoms., Comment: 7 pages, 3 figures

We predict that unpolarized charge current injected into a ballistic thin film of prototypical topological insulator (TI) Bi$_2$Se$_3$ will generate a {\it noncollinear spin texture} $\mathbf{S}(\mathbf{r})$ on its surface. Furthermore, the nonequilibrium spin texture will extend into $\simeq 2$ nm thick layer below the TI surfaces due to penetration of evanescent wavefunctions from the metallic surfaces into the bulk of TI. Averaging $\mathbf{S}(\mathbf{r})$ over few \AA{} along the longitudinal direction defined by the current flow reveals large component pointing in the transverse direction. In addition, we find an order of magnitude smaller out-of-plane component when the direction of injected current with respect to Bi and Se atoms probes the largest hexagonal warping of the Dirac-cone dispersion on TI surface. Our analysis is based on an extension of the nonequilibrium Green functions combined with density functional theory (NEGF+DFT) to situations involving noncollinear spins and spin-orbit coupling. We also demonstrate how DFT calculations with properly optimized local orbital basis set can precisely match putatively more accurate calculations with plane-wave basis set for the supercell of Bi$_2$Se$_3$., Comment: 5 pages, 4 figures, PDFLaTeX

Using first-principles quantum transport simulations, based on the nonequilibrium Green function formalism combined with density functional theory (NEGF+DFT), we examine changes in the total and local electronic currents within the plane of graphene nanoribbon with zigzag edges (ZGNR) hosting a nanopore which are induced by inserting a DNA nucleobase into the pore. We find a sizable change of the zero-bias conductance of two-terminal ZGNR + nanopore device after the nucleobase is placed into the most probable position (according to molecular dynamics trajectories) inside the nanopore of a small diameter \mbox{$D=1.2$ nm}. Although such effect decreases as the nanopore size is increased to \mbox{$D=1.7$ nm}, the contrast between currents in ZGNR + nanopore and ZGNR + nanopore + nucleobase systems can be enhanced by applying a small bias voltage $V_b \lesssim 0.1$ V. This is explained microscopically as being due to DNA nucleobase-induced modification of spatial profile of local current density around the edges of ZGNR. We repeat the same analysis using NEGF combined with self-consistent charge density functional tight-binding (NEGF+SCC-DFTB) or self-consistent extended H\"{u}ckel (NEGF+SC-EH) semi-empirical methodologies. The large discrepancy we find between the results obtained from NEGF+DFT vs. those obtained from NEGF+SCC-DFTB or NEGF+SC-EH approaches could be of great importance when selecting proper computational algorithms for {\em in silico} design of optimal nanoelectronic sensors for rapid DNA sequencing., Comment: 16 pages, 7 figures, PDFLaTeX; prepared for the special issue of the Journal of Computational Electronics on "Solid State Membranes for Biomolecule Sensing and Manipulation"

Designing thermoelectric materials with high figure of merit $ZT=S^2 G T/\kappa$ requires fulfilling three often irreconcilable conditions, i.e., the high electrical conductance $G$, small thermal conductance $\kappa$ and high Seebeck coefficient $S$. Nanostructuring is one of the promising ways to achieve this goal as it can substantially suppress lattice contribution to $\kappa$. However, it may also unfavorably influence the electronic transport in an uncontrollable way. Here we theoretically demonstrate that this issue can be ideally solved by fabricating graphene nanoribbons with heavy adatoms and nanopores. These systems, acting as a two-dimensional topological insulator with robust helical edge states carrying electrical current, yield a highly optimized power factor $S^2G$ per helical conducting channel. Concurrently, their array of nanopores impedes the lattice thermal conduction through the bulk. Using quantum transport simulations coupled with first-principles electronic and phononic band structure calculations, the thermoelectric figure of merit is found to reach its maximum $ZT \simeq 3$ at $T \simeq 40$ K. This paves a way to design high-$ZT$ materials by exploiting the nontrivial topology of electronic states through nanostructuring., Comment: 7 pages, 4 figures; PDFLaTeX

We study the transverse spin-Seebeck effect (SSE) on the surface of a three-dimensional topological insulator (TI) thin film, such as Bi$_2$Se$_3$, which is sandwiched between two normal metal leads. The temperature bias $\Delta T$ applied between the leads generates surface charge current which becomes spin-polarized due to strong spin-orbit coupling on the TI surface, with polarization vector acquiring a component $P_x \simeq 60%$ {\em parallel to the direction of transport}. When the third nonmagnetic voltage probe is attached to the portion of the TI surface across its width $L_y$, pure spin current will be injected into the probe where the inverse spin Hall effect (ISHE) converts it into a voltage signal \mbox{$|V_\mathrm{ISHE}|^\mathrm{max}/\Delta T \simeq 2.5$ $\mu$V/K} (assuming the SH angle of Pt voltage probe and $L_y=1$ mm). The existence of predicted nonequilibrium spin-polarization parallel to the direction of electronic transport and the corresponding electron-driven SSE crucially relies on orienting quintuple layers (QLs) of Bi$_2$Se$_3$ {\em orthogonal} to the TI surface and {\em tilted} by $45^\circ$ with respect to the direction of transport. Our analysis is based on the Landauer-B\"{u}ttiker-type formula for spin currents in the leads of a multi-terminal quantum-coherent junction, which is constructed using nonequilibrium Green function formalism within which we show how to take into account arbitrary orientation of QLs via the self-energy describing coupling between semi-infinite normal metal leads and TI., Comment: 8 pages, 4 color figures; PDFLaTeX

We analyze electronic and phononic quantum transport through zigzag or chiral graphene nanoribbons (GNRs) perforated with an array of nanopores. Since local charge current profiles in these GNRs are peaked around their edges, drilling nanopores in their interior does not affect such edge charge currents while drastically reducing heat current carried by phonons in sufficiently long wires. The combination of these two effects can yield highly efficient thermoelectric devices with maximum $ZT \simeq 11$ at liquid nitrogen temperature and $ZT \simeq 4$ at room temperature achieved in $\sim 1$ $\mu$m long zigzag GNRs with nanopores of variable diameter and spacing between them. Our analysis is based on the $\pi$-orbital tight-binding Hamiltonian with up to third nearest-neighbor hopping for electronic subsystem, the empirical fourth-nearest-neighbor model for phononic subsystem, and nonequilibrium Green function formalism to study quantum transport in both of these models., Comment: 5 pages, 5 figures, PDFLaTeX