Search

Half-Heusler compounds are known for their various compositions and multifunctional properties including topological phases. In this study, we investigate the topological classification of this class of materials based on the ordering of the $\Gamma_6$, $\Gamma_8$ states, and a previously overlooked $\Gamma_6^*$ state, during an adiabatic expansion process. Using first-principles calculations based on density functional theory, we observed that the non-trivial topology is governed by a three-band mechanism. We provide a simple model derived from $k\cdot p$ Hamiltonian that interprets the topological phase in half-Heusler systems. Additionally, we explore the robustness of the topological phase under tension and a new perspective on the topological nature of half-Heusler compounds.

Machine-learning models can be trained to predict the converged electron charge density of a density functional theory (DFT) calculation. In general, the value of the density at a given point in space is invariant under global translations and rotations having that point as a centre. Hence, one can construct locally invariant machine-learning density predictors. However, the widely used projector augmented wave (PAW) implementation of DFT requires the evaluation of the one-center augmentation contributions, that are not rotationally invariant. Building on our recently proposed Jacobi-Legendre charge-density scheme, we construct a covariant Jacobi-Legendre model capable of predicting the local occupancies needed to compose the augmentation charge density. Our formalism is then applied to the prediction of the energy barrier for the 1H-to-1T phase transition of two-dimensional MoS$_2$. With extremely modest training, the model is capable of performing a non-self-consistent nudged elastic band calculation at virtually the same accuracy as a fully DFT-converged one, thus saving thousands of self-consistent DFT steps. Furthermore, at variance with machine-learning force fields, the charge density is here available for any nudged elastic band image, so that we can trace the evolution of the electronic structure across the phase transition., Comment: 12 pages, 4 figures

We report experimental investigations of spin-to-charge current conversion and charge transfer dynamics (CT) at the interface of graphene/WS$_2$ van der Waals heterostructure. Pure spin current was produced by the spin precession in the microwave-driven ferromagnetic resonance of a permalloy film (Py-Ni$_{81}$Fe$_{19}$) and injected into the graphene/WS$_2$ heterostructure through the spin pumping process. The observed spin-to-charge current conversion in the heterostructure is attributed to inverse Rashba-Edelstein effect (IREE) at the graphene/WS$_2$ interface. Interfacial CT dynamics in this heterostructure was investigated based on the framework of core-hole-clock (CHC) approach. The results obtained from spin pumping and CHC studies show that the spin-to-charge current conversion and charge transfer process are more efficient in the graphene/WS$_2$ heterostructure compared to isolated WS2 and graphene films. The results show that the presence of WS$_2$ flakes improves the current conversion efficiency. These experimental results are corroborated by density functional theory (DFT) calculations, which reveal (i) Rashba spin-orbit splitting of graphene orbitals and (ii) electronic coupling between graphene and WS$_2$ orbitals. This study provides valuable insights for optimizing the design and performance of spintronic devices., Comment: 25 pages, 7 figures

The discovery of the quantum spin Hall effect led to the exploration of the electronic transport for spintronic devices. Here, we theoretically investigated the electronic conductance in large-gap realistic quantum spin Hall system, Pt$_2$HgSe$_3$ nanoribbons. By an ab initio approach, we found that the edge states present a penetration depth of about $0.9$\,{nm}, which is much smaller than those predicted in other 2D topological systems. Thus, suggesting that Pt$_2$HgSe$_3$ allows the exploitation of topological transport properties in narrow ribbons. Using non-equilibrium Green's functions calculations, we have examined the electron conductivity upon the presence of Se\,$\leftrightarrow$\,Hg antistructure defects randomly distributed in the Pt$_2$HgSe$_3$ scattering region. By considering scattering lengths up to $109$\,nm, we found localization lengths that can surpass $\mu$m sizes for narrow nanoribbons ($<9$\,nm). These findings can contribute to further understanding the behavior of topological insulators under realistic conditions and their integration within electronic, spintronic devices.

In recent years, synthesis and experimental research of fractalized materials has evolved in a paradigmatic crossover with topological phases of matter. We present here a theoretical investigation of the helical edge transport in Sierpinski carpets (SCs), combining the Bernevig-Hughes-Zhang (BHZ) model and the Landauer approach. Starting from a pristine two-dimensional topological insulator (2DTI), according to the BHZ model, our results reveal resonant transport modes when the SC fractal generation reaches the same scale as the space discretization; these modes are analyzed within a contour plot mapping of the local spin-polarized currents, shown spanned and assisted by inner-edge channels. From such a deeply fractalized SC building block, we introduce a rich tapestry formed by superior SC hierarchies, enlightening intricate patterns and unique fingerprints that offer valuable insights into how helical edge transport occurs in these fractal dimensions.

We investigate the topological phase transitions driven by band warping and a transverse magnetic field, for three-dimensional Weyl semimetals. First, we use the Chern number as a mathematical tool to derive the topological phase diagram. Next, we associate each of the topological sectors to a given angular momentum state of a rotating wave packet. Then we show how the position of the Weyl nodes can be manipulated by a transverse external magnetic field that ultimately quenches the wave packet rotation, first partially and then completely, thus resulting in a sequence of field-induced topological phase transitions. Finally, we calculate the current-induced magnetization and the anomalous Hall conductivity of a prototypical warped Weyl material. Both observables reflect the topological transitions associated with the wave packet rotation and can help to identify the elusive 3D quantum anomalous Hall effect in three-dimensional, warped Weyl materials., Comment: 6 pages, 5 figures

Spin-orbital textures in topological insulators due to the spin locking with the electron momentum, play an important role in spintronic phenomena that arise from the interplay between charge and spin degrees of freedom. We have explored interfaces between a ferromagnetic system (CrI$_3$) and a topological insulator (Bi$_2$Se$_3$) that allow the manipulation of spin-orbital textures. Within an {\it ab initio} approach we have extracted the spin-orbital-textures dependence of experimentally achievable interface designs. The presence of the ferromagnetic system introduces anisotropic transport of the electronic spin and charge. From a parameterized Hamiltonian model we capture the anisotropic backscattering behavior, showing its extension to other ferromagnetic/topological insulator interfaces. We verified that the van der Waals TI/MI interface is an excellent platform for controlling the spin degree of freedom arising from topological states, providing a rich family of unconventional spin texture configurations.

The catalytic sites of MoS$_{2}$ monolayers towards hydrogen evolution are well known to be vacancies and edge-like defects. However, it is still very challenging to control the position, size, and defective areas on the basal plane of Mo$S_{2}$ monolayers by most of defect-engineering routes. In this work, the fabrication of etched arrays on ultra-large supported and free-standing MoS$_{2}$ monolayers using focused ion beam (FIB) is reported for the first time. By tuning the Ga+ ion dose, it is possible to confine defects near the etched edges or spread them over ultra-large areas on the basal plane. The electrocatalytic activity of the arrays toward hydrogen evolution reaction (HER) was measured by fabricating microelectrodes using a new method that preserves the catalytic sites. We demonstrate that the overpotential can be decreased up to 290 mV by assessing electrochemical activity only at the basal plane. High-resolution transmission electron microscopy images obtained on FIB patterned freestanding MoS$_{2}$ monolayers reveal the presence of amorphous regions and X-ray photoelectron spectroscopy indicates sulfur excess in these regions. Density-functional theory calculations provide identification of catalytic defect sites. Our results demonstrate a new rational control of amorphous-crystalline surface boundaries and future insight for defect optimization in MoS$_{2}$ monolayers., Comment: 10 pages, 5 figures

Topological insulators have remained as candidates for future electronic devices since their first experimental realization in the past decade. The existence of topologically protected edge states could be exploited to generate a robust platform and develop quantum computers. In this work we explore the role of magnetic impurities in the transport properties of topological insulators, in particular, we study the effect on the edge states conductivity. By means of realistic $\it{ab}$ $\it{initio}$ calculations we simulate the interaction between magnetic adatoms and topological insulators, furthermore, our main goal is to obtain the transport properties for large samples as it would be possible to localize edge states at large scales.

2D materials present an interesting platform for device designs. However, oxidation can drastically change the system's properties, which need to be accounted for. Through {\it ab initio} calculations, we investigated freestanding and SiC-supported As, Sb, and Bi mono-elemental layers. The oxidation process occurs through an O$_2$ spin-state transition, accounted for within the Landau-Zener transition. Additionally, we have investigated the oxidation barriers and the role of spin-orbit coupling. Our calculations pointed out that the presence of SiC substrate reduces the oxidation time scale compared to a freestanding monolayer. We have extracted the energy barrier transition, compatible with our spin-transition analysis. Besides, spin-orbit coupling is relevant to the oxidation mechanisms and alters time scales. The energy barriers decrease as the pnictogen changes from As to Sb to Bi for the freestanding systems, while for SiC-supported, they increase across the pnictogen family. Our computed energy barriers confirm the enhanced robustness against oxidation for the SiC-supported systems.

Two-dimensional arsenic allotropes have been grown on metallic surfaces, while topological properties have been theoretically described on strained structures. Here we experimentally grow arsenene by molecular beam epitaxy over the insulating SiC substrate. The arsenene presents a flat structure with a strain field that follows the SiC surface periodicity. Our ab initio simulations, based on the density functional theory, corroborate the experimental observation. The strained structure presents a new arsenene allotrope with a triangular structure, rather than the honeycomb previously predicted for other pnictogens. This strained structure presents a Peierls-like transition leading to an indirect gap semiconducting behavior.

Kohn-Sham density functional theory (KS-DFT) is a powerful method to obtain key materials' properties, but the iterative solution of the KS equations is a numerically intensive task, which limits its application to complex systems. To address this issue, machine learning (ML) models can be used as surrogates to find the ground-state charge density and reduce the computational overheads. We develop a grid-centred structural representation, based on Jacobi and Legendre polynomials combined with a linear regression, to accurately learn the converged DFT charge density. This integrates into a ML pipeline that can return any density-dependent observable, including energy and forces, at the quality of a converged DFT calculation, but at a fraction of the computational cost. Fast scanning of energy landscapes and producing starting densities for the DFT self-consistent cycle are among the applications of our scheme.

Transition metal dichalcogenides (TMDCs) are promising materials for applications in nanoelectronics and correlated fields, where their metallic edge states play a fundamental role in the electronic transport. In this work, we investigate the transport properties of MoS$_2$ zigzag nanoribbons under a butadiene (C$_4$H$_6$) atmosphere, as this compound has been used to obtain MoS$_2$ flakes by exfoliation. We use density functional theory combined with non-equilibrium Green's function techniques, in a methodology contemplating disorder and different coverages. Our results indicate a strong modulation of the TMDC electronic transport properties driven by butadiene molecules anchored at their edges, producing the suppression of currents due to a backscattering process. Our results indicate a high sensitivity of the TMDC edge states. Thus, the mechanisms used to reduce the dimensionality of MoS$_2$ considerably modify its transport properties.

Transition metal dichalcogenides have been the subject of numerous studies addressing technological applications and fundamental issues. Single-layer PtSe2 is a semiconductor with a trivial bandgap, in contrast, its counterpart with 25% of Se atoms substituted by Hg, Pt2HgSe3 (jacutingaite, a naturally occurring mineral), is a 2D topological insulator with a large bandgap. Based on ab-initio calculations, we investigate the energetic stability, and the topological transition in Pt(HgxSe1-x)2 as a function of alloy concentration, and the distribution of Hg atoms embedded in the PtSe2 host. Our findings reveal the dependence of the topological phase with respect to the alloy concentration and robustness with respect distribution of Hg. Through a combination of our ab-initio results and a defect wave function percolation model, we estimate the random alloy concentration threshold for the topological transition to be only 9%. Our results expand the possible search for non-trivial topological phases in random alloy systems.

The development of electronic devices based on the functionalization of (nano)cellulose platforms relies upon an atomistic understanding of the structural, and electronic properties of the combined system, cellulose/functional element. In this work, we present a theoretical study of the nanocellulose/graphene interface (nCL/G) based on first-principles calculations. We find that the binding energies of both hydrophobic/G (nCL$^{\rm phob}$/G) and hydrophilic/G (nCL$^{\rm phil}$/G) interfaces are primarily dictated by the van der Waals interactions, and are comparable with that of their 2D interface counterparts. We verify that the energetic preference of nCL$^{\rm phob}$/G has been reinforced by the inclusion of an aqueous media via the implicit solvation model. Further structural characterization was carried out using a set of simulations of Carbon K-edge X-ray absorption spectra to identify and distinguish the key absorption features of the nCL$^{\rm phob}$/G and nCL$^{\rm phil}$/G interfaces. The electronic structure calculations reveal that the linear energy bands of graphene lie in the band gap of the nCL, sheet, while depletion/accumulation charge density regions are observed. We show that external agents, i.e. electric field and mechanical strain, allow for tunability of the Dirac cone and the charge density at the interface. The control/maintenance of the Dirac cone states in nCL/G is an important feature for the development of electronic devices based on cellulosic platforms.

Two-dimensional transition metal dichalcogenides (MX$_2$) vacancy formation energetics is extensively investigated. Within an ab-initio approach we study the MX$_2$ systems, with M=Mo, W, Ni, Pd and Pt, and X=S, Se, and Te. Here we classify that chalcogen vacancies are always energetic favorable over the transition metal ones. However, for late transition metals Pd $4d$, and Pt $5d$ the metal vacancy are experimentally achievable, bringing up localized magnetic moments within the semiconducting matrix. By quantifying the localization of the chalcogen vacancy states we evidentiate that it rules the intra- and inter-vacancy interactions that establish both the number of vacancy states neatly lying within the semiconducting gap, as well as its electronic dispersion and SOC splitting. Combining different vacancies and phase variability 1T and 1H of the explored systems allow us to construct a guiding picture for the vacancy states localization.

Monolayers of transition metal dichalcogenides (TMDs) in the 2H structural phase have been recently classified as higher-order topological insulators (HOTI), protected by $C_3$ rotation symmetry. In addition, theoretical calculations show an orbital Hall plateau in the insulating gap of TMDs, characterized by an orbital Chern number. We explore the correlation between these two phenomena in TMD monolayers in two structural phases: the noncentrosymmetric 2H and the centrosymmetric 1T. Using density functional theory, we confirm the characteristics of 2H-TMDs and reveal that 1T-TMDs are identified by a $\mathbb{Z}_4$ topological invariant. As a result, when cut along appropriate directions, they host conducting edge-states, which cross their bulk energy-band gaps and can transport orbital angular momentum. Our linear response calculations thus indicate that the HOTI phase is accompanied by an orbital Hall effect. Using general symmetry arguments, we establish a connection between the two phenomena with potential implications for spin-orbitronics., Comment: 5 pages, 4 figures

Molecular linkers have emerged as an effective strategy to improve electronic transport properties on solution-processed layered materials via defect functionalization. However, a detailed discussion on the microscopic mechanisms behind the beneficial effects of functionalization is still missing. Here, by first-principles calculations based on density functional theory, we investigate the effects on the electronic properties of interconnected MoS$_2$ model flakes systems upon functionalization with different thiol molecule linkers, namely thiophenol, 1,4-benzenedithiol, 1,2-ethanedithiol, and 1,3-propanedithiol. The bonding of benzene- and ethanedithiol bridging adjacent armchair MoS$_2$ nanoflakes leads to electronic states just above or at the Fermi level, thus forming a molecular channel for electronic transport between flakes. In addition, the molecular linker reduces the potential barrier for thermally activated hopping between neighboring flakes, improving the conductivity as verified in experiments. The comprehension of such mechanisms helps in future developments of solution-processed layered materials for use on 2D electronic devices.

The development of spintronic devices demands the existence of materials with some kind of spin splitting (SS). In this Data Descriptor, we build a database of ab initio calculated SS in 2D materials. More than that, we propose a workflow for materials design integrating an inverse design approach and a Bayesian inference optimization. We use the prediction of SS prototypes for spintronic applications as an illustrative example of the proposed workflow. The prediction process starts with the establishment of the design principles (the physical mechanism behind the target properties), that are used as filters for materials screening, and followed by density functional theory (DFT) calculations. Applying this process to the C2DB database, we identify and classify 358 2D materials according to SS type at the valence and/or conduction bands. The Bayesian optimization captures trends that are used for the rationalized design of 2D materials with the ideal conditions of band gap and SS for potential spintronics applications. Our workflow can be applied to any other material property., Comment: Manuscript: 506 pages

Vacancies in materials structure -- lowering its atomic density -- take the system closer to the atomic limit, to which all systems are topologically trivial. Here we show a mechanism of mediated interaction between vacancies inducing a topologically non-trivial phase. Within an {\it ab initio} approach we explore topological transition dependence with the vacancy density in transition metal dichalcogenides. As a case of study, we focus on the PtSe$_2$, to which pristine form is a trivial semiconductor with an energy gap of $1.2$\,eV. The vacancies states lead to a large topological gap of $180$\,meV within the pristine system gap. We derive an effective model describing this topological phase in other transition metal dichalcogenide systems. The mechanism driving the topological phase allows the construction of backscattering protected metallic channels embedded in a semiconducting host.

Crystalline $\rm Bi_2Se_3$ is one of the most explored three-dimensional topological insulator, with a $0.3\;\rm eV$ energy gap making it promising for applications. Its amorphous counterpart could bring to light new possibilities for large scale synthesis and applications. Using ab initio molecular dynamics simulations, we have studied realistic amorphous $\rm Bi_2Se_3$ phases generated by different processes of melting, quenching, and annealing. Extensive structural and electronic characterizations show that the melting process induces an energy gap decrease ruled by growth of the defective local environments. This behavior dictates a weak stability of the topological phase to disorder, characterized by the spin Bott index. Interestingly, we identify the occurrence of topologically trivial surface states in amorphous $\rm Bi_2Se_3$ that show a strong resemblance with standard helical topological states. Our results and methods advance the search of topological phases in three-dimensional amorphous solids.

One of the main goals and challenges of materials discovery is to find the best candidates for each interest property or application. Machine learning rises in this context to efficiently optimize this search, exploring the immense materials space, consisting of simultaneously the atomic, compositional, and structural spaces. Topological insulators, presenting symmetry-protected metallic edge states, are a promising class of materials for different applications. However, further, development is limited by the scarcity of viable candidates. Here we present and discuss machine learning-accelerated strategies for searching the materials space for two-dimensional topological materials. We show the importance of detailed investigations of each machine learning component, leading to different results. Using recently created databases containing thousands of ab initio calculations of 2D materials, we train machine learning models capable of determining the electronic topology of materials, with an accuracy of over 90%. We can then generate and screen thousands of novel materials, efficiently predicting their topological character without the need for a priori structural knowledge. We discover 56 non-trivial materials, of which 17 novel insulating candidates for further investigation, for which we corroborate their topological properties with density functional theory calculations. This strategy is 10$\times$ more efficient than the trial-and-error approach while few orders of magnitude faster and is a proof of concept for guiding improved materials discovery search strategies., Comment: Accepted for publication in Applied Physical Review. For Supplemental Materials please contact one of the authors

Understanding the role of microscopic attributes in nanocomposites allows for a controlled and, therefore, acceleration in experimental system designs. In this work, we extracted the relevant parameters controlling the graphene oxide binding strength to cellulose by combining first-principles calculations and machine learning algorithms. We were able to classify the systems among two classes with higher and lower binding energies, which are well defined based on the isolated graphene oxide features. By a theoretical X-ray photoelectron spectroscopy analysis, we show the extraction of these relevant features. Additionally, we demonstrate the possibilities of a refined control within a machine learning regression between the binding energy values and the system's characteristics. Our work presents a guiding map to the control graphene oxide/cellulose interaction.

Polyphenols are natural molecules of crucial importance in many applications, of which tannic acid (TA) is one of the most abundant and established. Most high-value applications require precise control of TA interactions with the system of interest. However, the molecular structure of TA is still not comprehended at the atomic level, of which all electronic and reactivity properties depend. Here, we combine an enhanced sampling global optimization method with density functional theory (DFT)-based calculations to explore the conformational space of TA assisted by unsupervised machine learning visualization, and then investigate its lowest energy conformers. We study the external environment's effect on the TA structure and properties. We find that vacuum favors compact structures by stabilizing peripheral atoms' weak interactions, while in water, the molecule adopts more open conformations. The frontier molecular orbitals of the conformers with lowest harmonic vibrational free energy have a HOMO-LUMO energy gap of 2.21 (3.27) eV, increasing to 2.82 (3.88) eV in water, at the DFT generalized gradient approximation (and hybrid) level of theory. Structural differences also change the distribution of potential reactive sites. We establish the fundamental importance of accurate structural consideration in determining TA and related polyphenols interactions in relevant technological applications.

Understanding the effects of atomic structure modification in hematite photoanodes is essential for the rational design of high-efficiency functionalizations. Recently it was found that interface modification with Sn/Sb segregates considerably increases hematite photocatalytic efficiency. However, the understanding of the different electronic effects of these modifications at the atomic level is still lacking. This letter describes the segregation effects of two different dopants-Sn and Sb-on both the solid-solid (grain-boundaries) and solid-liquid interfaces (surfaces) of hematite. Within an ab-initio approach, we quantitatively extract the potential barrier reduction on polycrystalline interfaces due to the dopant, which causes an increase in the inter-grain electron transport. Concomitantly, the dopants' segregation on hematite surfaces results in a decrease of the oxygen vacancy formation energy. Such vacancies lead to the experimentally observed rise of the flat-band potential. The comprehension of the electronic effects of dopants on both types of interfaces explains the experimental peak efficiency of interface-modified hematite with dopant segregates, also enabling the control and design of interfaces for different higher-efficiency applications.

In recent years, the notion of Quantum Materials has emerged as a powerful unifying concept across diverse fields of science and engineering, from condensed-matter and cold atom physics to materials science and quantum computing. Beyond traditional quantum materials such as unconventional superconductors, heavy fermions, and multiferroics, the field has significantly expanded to encompass topological quantum matter, two-dimensional materials and their van der Waals heterostructures, Moire materials, Floquet time crystals, as well as materials and devices for quantum computation with Majorana fermions. In this Roadmap collection we aim to capture a snapshot of the most recent developments in the field, and to identify outstanding challenges and emerging opportunities. The format of the Roadmap, whereby experts in each discipline share their viewpoint and articulate their vision for quantum materials, reflects the dynamic and multifaceted nature of this research area, and is meant to encourage exchanges and discussions across traditional disciplinary boundaries. It is our hope that this collective vision will contribute to sparking new fascinating questions and activities at the intersection of materials science, condensed matter physics, device engineering, and quantum information, and to shaping a clearer landscape of quantum materials science as a new frontier of interdisciplinary scientific inquiry., Comment: 93 pages, 27 figures

Recently cycloarene has been experimentally obtained in a self-assembled structure, forming graphene-like monoatomic layered systems. Here, we establish the bandgap engineering/prediction in cycloarene assemblies within a combination of density functional theory and tight-binding Hamiltonians. Our results show a weak dependence of the gap with the assembly geometry, contrasting a strong dependence with the inter-molecule bond density. We derived a effective model that allows the interpretation of the arising energy gap for general particle-hole symmetric molecular arranges based on inter-molecular bond strength.

Material's geometrical structure is a fundamental part of their properties. The honeycomb geometry of graphene is responsible for the arising of its Dirac cone, while the kagome and Lieb lattice hosts flat bands and pseudospin-1 Dirac dispersion. These features seem to be particular for few 2D systems rather than a common occurrence. Given this correlation between structure and properties, exploring new geometries can lead to unexplored states and phenomena. Kepler is the pioneer of the mathematical tiling theory, describing ways of filing the euclidean plane with geometrical forms in its book {\it Harmonices Mundi}. In this letter, we characterize $1255$ lattices composed of the euclidean plane's k-uniform tiling, with its intrinsic properties unveiled - this class of arranged tiles present high-degeneracy points, exotic quasiparticles, and flat bands as a common feature. Here, we present aid for experimental interpretation and prediction of new 2D systems.

Proximity effects are one of the pillars of exotic phenomena and technological applications of two dimensional materials. However, the interactions nature depends strongly on the materials involved, their crystalline symmetries, and interfacial properties. Here we used large-scale first-principle calculations to demonstrate that strain and twist-angle are efficient knobs to tailor the spin-orbit coupling in graphene transition metal dichalcogenide heterobilayers. We found that by choosing a twist-angle of 30 degrees, the spin relaxation times increase by two orders of magnitude, opening a path to improve these heterostructures spin transport capability. Moreover, we demonstrate that strain and twist angle will modify the relative values of valley-Zeeman and Rashba spin-orbit coupling, allowing to tune the system into an ideal Dirac-Rashba regime. These results enable us to envision an answer for the variability of spin-orbit coupling found in different experiments and have significant consequences for applications that depend on polycrystallinity, where grains form at different orientations.

Cellulose disassembly is an important issue in designing nanostructures using cellulose-based materials. In this work, we present a joint of experimental and theoretical study addressing the disassembly of cellulose nanofibrils. Through 2,2,6,6-tetramethylpiperidine-1-oxyl (TEMPO) mediated oxidation processes, combined with atomic force microscopy results, we find the formation of nanofibers with diameters corresponding to a single cellulose polymer chain. The formation of these polymer chains is ruled by repulsive electrostatic interactions between the oxidized chains. Further first-principles calculations have been done in order to provide an atomistic understanding the cellulose disassembling processes, focusing on the balance of the interchain and intersheet interactions upon oxidation. Firstly we analyse these interaction in pristine systems, where we found the intersheet interaction stronger than the interchain one. In the oxidized systems, we have considered the formation of (charged) carboxylate groups along the inner sites of elementary fibrils. We show a net charge concentration on the carboxylate groups, supporting the emergence of repulsive electrostatic interactions between the cellulose nanofibers. Indeed, our total energy results show that the weakening of the binding strength between the fibrils is proportional to the concentration and the net charge density of the carboxylate group. Moreover, by comparing interchain and intersheet binding energies, we found that most of the disassembly processes should take place by breaking the interchain O--H$\cdots$O hydrogen bond interactions, and thus supporting the experimental observation of single and double cellulose polymeric chains., Comment: 9 pages, 5 figures

Jacutingate, a recently discovered Brazilian naturally occurring mineral, has shown to be the first experimental realization of the Kane-Mele topological model. In this letter we have unveiled a class of materials $M_2NX_3$ ($M$=Ni, Pt, Pd; $N$=Zn, Cd, Hg; and $X$=S, Se, Te), sharing jacutingaite's key features, i.e., high stability, and topological phase. By employing first-principles calculations we extensively characterize the energetic stability of this class while showing a common occurrence of the Kane-Mele topological phase. Here we found Pt-based materials surpassing jacutingaite's impressive topological gap and lower exfoliation barrier while retaining its stability.

We investigate the structure and electronic spectra properties of two-dimensional amorphous bismuthene structures and show that these systems are topological insulators. We employ realistic modeling of amorphous geometries together with density functional theory for electronic structure calculations. We investigate the system topological properties throughout the amorphization process and find that the robustness of the topological phase is associated with the spin-orbit coupling strength and size of the pristine topological gap. Using recursive non-equilibrium Green's function, we study the electronic transport properties of nanoribbons devices with lengths comparable to experimentally synthesized materials. We find a $2e^2/h$ conductance plateau within the topological gap and an onset of Anderson localization at the trivial insulator phase.

The robustness of topological materials against disorder and defects is presumed but has not been demonstrated explicitly in realistic systems. In this work, we use state-of-the-art density functional theory and recursive nonequilibrium Green's functions methods to study the effect of disorder in the electronic transport of long nanoribbons, up to 157 nm, as a function of vacancy concentration. In narrow nanoribbons, even for small vacancy concentrations, defect-like localized states give rise to hybridization between the edge states erasing topological protection and enabling backscattering events. We show that the topological protection is more robust for wide nanoribbons, but surprisingly it breaks down at moderate structural disorder. Our study helps to establish some bounds on defective bismuthene nanoribbons as promising candidates for spintronic applications.

It has recently been demonstrated that N-heterocyclic carbenes (NHCs) form self-assembled monolayers (SAMs) on metal surfaces. Consequently, it is important to both characterize and understand their binding modes to fully exploit NHCs in functional surface systems. To assist with this effort, we have performed {\it first-principles} total energy calculations for NHCs on Au(111) and simulations of X-ray absorption near edge structure (XANES). The NHCs we have considered are N,N-dimethyl-, N,N-diethyl-, N,N-diisopropylbenzimidazolylidene ($^B$NHC$^X$, with X=Me, Et, and iPr, respectively) and the bis-$^B$NHC$^X$ complexes with Au derived from these molecules. We present a comprehensive analysis of the energetic stability of both the $^B$NHC$^X$ and the complexes on Au(111) and, for the former, examine the role of the wing group in determining the attachment geometry. Further structural characterization is performed by calculating the nitrogen K-edge X-ray absorption spectra. Our simulated XANES results give insight into (i) the relationship between the $^B$NHC$^X$/Au geometry and the N($1s$) $\rightarrow$ $\pi^\ast/\sigma^\ast$, pre-edge/near-edge, absorption intensities, and (ii) the contributions of the molecular deformation and molecule-surface electronic interaction to the XANES spectrum. Our simulations are compared with recent experimental results.

Two-dimensional Na$_3$Bi is a dual topological insulator protected by time-reversal and mirror symmetry, resulting in a promising platform for devices design. However, in reality, the design of topological devices is hindered by a sensitivity against disorder and temperature. We study the topological properties of Na$_3$Bi in the presence of intrinsic defects, investigating the robustness of the edge states and the resulting transport properties. We apply a recursive Green's function technique enabling the study of disordered systems with lengths comparable to experimentally synthesized materials, in the order of micrometers. We combine our findings to propose a topological insulator device, where intrinsic defects are used to filter the response of trivial bulk states. This results in a stable conductance throughout a large range of electronic temperatures, and controllable by a perpendicular electric field. Our proposal is general, enabling the design of various dual topological insulators devices.

The discovery and realization of topological insulators, a phase of matter which hosts metallic boundary states when the $d$-dimension insulating bulk is confined to ($d-1$)-dimensions, led to several potential applications. Recently, it was shown that protected topological states can manifest in ($d-2$)-dimensions, such as hinge and corner states for three- and two-dimensional systems, respectively. These nontrivial materials are named higher-order topological insulators (HOTIs). Here we show a connection between spin Hall effect and HOTIs using a combination of {\it ab initio} calculations and tight-binding modeling. The model demonstrates how a non-zero bulk midgap spin Hall conductivity (SHC) emerges within the HOTI phase. Following this, we performed high-throughput density functional theory calculations to find unknown HOTIs, using the SHC as a criterion. We calculated the SHC of 693 insulators resulting in seven stable two-dimensional HOTIs. Our work guides novel experimental and theoretical advances towards higher-order topological insulators realization and applications.

The spin-orbit -induced spin splitting of energy bands in low symmetry compounds (the Rashba Effect) has a long-standing relevance to spintronic applications and to the fundamental understanding of symmetry breaking in solids, yet the knowledge of what controls its magnitude in different materials is difficult to anticipate. Indeed, rare discoveries of compounds with large Rashba coefficients are invariably greeted as pleasant surprises. We advance the understanding of the "Rashba Scale" using the "inverse design" approach by formulating theoretically the relevant design principle and then identifying compounds that satisfy it. We show that the presence of energy band anti-crossing provides a causal design principle of compounds with large Rashba coefficients, leading to the identification via first-principles calculations of 34 rationally designed strong-Rashba compounds. Since topological insulators must have band anti crossing, this leads us to establish an interesting cross functionality of "Topological Rashba Insulators" (TRI) that may provide a platform for the simultaneous control of spin splitting and spin-polarization., Comment: 24 pages and 9 figures

The topological properties of materials are, until now, associated with the features of their crystalline structure, although translational symmetry is not an explicit requirement of the topological phases. Recent studies of hopping models on random lattices have demonstrated that amorphous model systems show a nontrivial topology. Using {\it ab initio} calculations we show that two-dimensional amorphous materials can also display topological insulator properties. More specifically, we present a realistic state-of-the-art study of the electronic and transport properties of amorphous bismuthene systems, showing that these materials are topological insulators. These systems are characterized by the topological index $\mathbb{Z}_{2}=1$ and bulk-edge duality, and their linear conductance is quantized, ${\cal G}=2e^{2}/h$, for Fermi energies within the topological gap. Our study opens the path to the experimental and theoretical investigation of amorphous topological insulator materials.

Despite its potential for device application, the non-magnetic Zeeman effect has only been predicted and observed in two-dimensional compounds. We demonstrate that non-centrosymmetric three-dimensional compounds can also exhibit a Zeeman-type spin splitting, allowing the splitting control by changing the growth direction of slabs formed by these compounds. We determine the required conditions for this effect: $i$) non-centrosymmetric including polar and non-polar point groups, $ii$) valence band maximum or conduction band minimum in a generic $k$-point, i.e., non-time-reversal-invariant momentum, and $iii$) zero magnetic moment. Using these conditions as filters, we perform a material screening to systematically search for these systems in the aflow-ICSD database. We find 20 material candidates featuring the Zeeman-type effect. We also found that the spin-splitting in confined systems can be controlled by an external electric field, which in turns can induce a metal-insulator transition. We believe that the Zeeman-type effect in three-dimensional compounds can potentially be used for spin-filtering devices., Comment: 9 pages, 3 figures

Three-dimensional topological insulators protected by both the time reversal (TR) and mirror symmetries were recently predicted and observed. Two-dimensional materials featuring this property and their potential for device applications have been less explored. We find that in these systems, the spin-polarization of edge states can be controlled with an external electric field breaking the mirror symmetry. This symmetry requires that the spin-polarization is perpendicular to the mirror plane, therefore, the electric field induces spin-polarization components parallel to the mirror plane. Since this field preserves the TR topological protection, we propose a transistor model using the spin-direction of protected edge states as a switch. In order to illustrate the generality of the proposed phenomena, we consider compounds protected by mirror planes parallel and perpendicular to the structure, e.g., Na$_3$Bi and half-functionalized (HF) hexagonal compounds, respectively. For this purpose, we first construct a tight-binding effective model for the Na$_3$Bi compound and predict that HF-honeycomb lattice materials are also dual topological insulators., Comment: 6 pages and 3 figures

We study the spin excitation spectra and the dynamical exchange coupling between iron adatoms on a Bi bilayer nanoribbon. We show that the topological character of the edge states is preserved in the presence of the magnetic adatoms. Nevertheless, they couple significantly to the edge spin currents, as witnessed by the large and long-ranged dynamical coupling we obtain in our calcula- tions. The large effective magnetocrystalline anisotropy of the magnetic adatoms combined with the transport properties of the topologically protected edge states make this system a strong candidate for implementation of spintronics devices and quantum information and/or computation protocols.

Quantum spin-Hall insulators (QSHIs), i.e., two-dimensional topological insulators (TIs) with a symmetry-protected band inversion, have attracted considerable scientific interest in recent years. In this work, we have computed the topological Z2 invariant for 220 functionalized honeycomb lattices that are isoelectronic to functionalized graphene. Besides confirming the TI character of well-known materials such as functionalized stanene, our study identifies 45 yet unreported QSHIs. We applied a compressed-sensing approach to identify a physically meaningful descriptor for the Z2 invariant that only depends on the properties of the material's constituent atoms. This enables us to draw a map of materials, in which metals, trivial insulators, and QSHI form distinct regions. This analysis yields fundamental insights in the mechanisms driving topological transitions. The transferability of the identified model is explicitly demonstrated for an additional set of honeycomb lattices with different functionalizations that are not part of the original set of 220 graphene-type materials used to identify the descriptor. In this class, we predict 74 more novel QSHIs that have not been reported in literature yet.

We have performed a systematic investigation of the formation of topologically protected boundary states (TPBS) in topological/normal insulators (TI/NI) heterostructures. Using a recently developed scheme to construct {\it ab-initio} tight-binding Hamiltonian matrices from density functional theory (DFT) calculations, we studied systems of realistic size with high accuracy and control over the relevant parameters such as TI and NI band alignment, NI gap and spin-orbit coupling strength. Our findings point to the existence of an NI critical thickness for the emergence of TPBS and to the importance of the band alignment between the TI and NI for the appearance of the TPBS. We chose Bi$_{2}$Se$_{3}$ as a prototypical case where the topological/normal insulator behavior is modeled by regions with/without spin-orbit coupling. Finally, we validate our approach comparing our model with fully relativistic DFT calculations for TI/NI heterostructures of Bi$_{2}$Se$_{3}$/Sb$_{2}$Se$_{3}$.

We put forward a tight-binding model for rhombohedral topological insulators materials with the space group $D^{5}_{3d}(R\bar{3}m)$. The model describes the bulk band structure of these materials over the whole Brillouin zone. Within this framework, we also describe the topological nature of surface states, characterized by a Dirac cone-like dispersion and the emergence of surface projected bulk states near to the Dirac-point in energy. We find that the breaking of the $R_{3}$ symmetry as one moves away from the $\Gamma$ point has an important role in the hybridization of the $p_x$, $p_y$, and $p_z$ atomic orbitals. In our tight-binding model, the latter leads to a band mixing matrix element ruled by a single parameter. We show that our model gives a good description of the strategies/mechanisms proposed in the literature to eliminate and/or energy shift the bulk states away from the Dirac point, such as stacking faults and the introduction of an external applied electric field.

The electronic and the structural properties of two dimensional van der Waals metal/semiconductor heterostructures have been investigated through first-principles calculations. We have considered the recently synthesized borophene [Science {\bf 350}, 1513 (2015)], and the planar boron sheets (S1 and S2) [Nature Chemistry {\bf 8}, 563 (2016)] as the 2D metal layer, and the transition metal dichalcogenides (TMDCs) \mose, and \wse\ as the semiconductor monolayer. We find that the energetic stability of those 2D metal/semiconductor heterojunctions is mostly ruled by the vdW interactions; however, chemical interactions also take place in borophene/TMDC. The electronic charge transfers at the metal/semiconductor interface has been mapped, where we find a a net charge transfer from the TMDCs to the boron sheets. Further electronic structure calculations reveal that the metal/semiconductor interfaces, composed by planar boron sheets S1 and S2, present a p-type Schottky barrier which can be tuned to a p-type ohmic contact upon an external electric field.

We proposed that the simultaneous presence of both Rashba and band inversion can lead to a Rashba-like spin-splitting formed by two bands with the same in-plane helical spin texture. Because of this unconventional spin texture, the backscattering is forbidden in edge and bulk conductivity channels. We propose a new non-centrosymmetric honeycomb-lattice quantum spin Hall (QSH) insulator family formed by the IV, V, and VII elements with this property. The system formed by Bi, Pb and I atoms is mechanically stable and has both a large Rashba spin-splitting of 60 meV and a large nontrivial band gap of 0.14 eV. Since the edge and the bulk states are protected by the TR symmetry, contrary to what happens in most doped QSH insulators, the bulk states do not contribute to the backscattering in the electronic transport, allowing the construction of a spintronic device with less energy loss.

Phosphorene is a novel two-dimensional material that can be isolated through mechanical exfoliation from layered black phosphorus, but unlike graphene and silicene, monolayer phosphorene has a large band gap. It was thus unsuspected to exhibit band inversion and the ensuing topological insulator behavior. It has recently attracted interest because of its proposed application as field effect transistors. Using first-principles calculations with applied perpendicular electric field F we predict a continuous transition from the normal insulator to a topological insulator and eventually to a metal as a function of F. The continuous tuning of topological behavior with electric field would lead to spin-separated, gapless edge states, i.e., quantum spins Hall effect. This finding opens the possibility of converting normal insulating materials into topological ones via electric field, and making a multi-functional field effect topological transistor that could manipulate simultaneously both spins and charge carrier., Comment: 19 pages, 4 figures

Topological insulators (TIs) are a new class of matter characterized by the unique electronic properties of an insulating bulk and metallic boundaries arising from non-trivial bulk band topology. While the surfaces of TIs have been well studied, the interface between TIs and semiconductors may not only be more technologically relevant but the interaction with non-topological states may fundamentally alter the physics. Here, we present a general model to show that such an interaction can lead to spin-momentum locked non-topological states, the Dirac cone can split in two, and the particle-hole symmetry can be fundamentally broken, along with their possible ramifications. Unlike magnetic doping or alloying, these phenomena occur without topological transitions or the breaking of time reversal symmetry. The model results are corroborated by first-principles calculations of the technologically relevant Bi$_2$Se$_3$ film van der Waals bound to a Se-treated GaAs substrate.

We show that germanene nanoroads embedded in a completely hydrogenated germanene (germanane) exhibits a quantum spin Hall effect (QSHE). These nanoroads can be obtained experimentally by local hydrogen dissociation from germanane. Using first principle calculations we predict that germanene nanoroads with zigzag interfaces show dissipationless conducting channels with in-plane and out-of-plane spin textures., Comment: 5 pages, 5 figures

We have performed an ab initio investigation of the electronic properties of the graphene sheet adsorbed by Ru adatoms (Ru/graphene). For a particular set of triangular arrays of Ru adatoms, we find the formation of four (spin-polarized) Dirac cones attributed to a suitable overlap between two hexagonal lattices: one composed by the C sites of the graphene sheet, and the other formed by the surface potential induced by the Ru adatoms. Upon the presence of spin-orbit coupling (SOC) nontrivial band gaps take place at the Dirac cones promoting several topological phases. Depending on the Ru concentration, the system can be topologically characterized among the phases i) Quantum Spin Hall (QSH), ii) Quantum Anomalous Hall (QAH), iii) metal iv) or trivial insulator. For each concentration, the topological phase is characterized by the ab-initio calculation of the Chern number., Comment: 8 pages, 6 figures