Theoretical and Computational Nanoscience

Group leader: Stephan Roche

Website

Publications

  • Blue emission at atomically sharp 1D heterojunctions between graphene and h-BN

    Kim G., Ma K.Y., Park M., Kim M., Jeon J., Song J., Barrios-Vargas J.E., Sato Y., Lin Y.-C., Suenaga K., Roche S., Yoo S., Sohn B.-H., Jeon S., Shin H.S. Nature Communications; 11 (1, 5359) 2020. 10.1038/s41467-020-19181-2. IF: 12.121

    Atomically sharp heterojunctions in lateral two-dimensional heterostructures can provide the narrowest one-dimensional functionalities driven by unusual interfacial electronic states. For instance, the highly controlled growth of patchworks of graphene and hexagonal boron nitride (h-BN) would be a potential platform to explore unknown electronic, thermal, spin or optoelectronic property. However, to date, the possible emergence of physical properties and functionalities monitored by the interfaces between metallic graphene and insulating h-BN remains largely unexplored. Here, we demonstrate a blue emitting atomic-resolved heterojunction between graphene and h-BN. Such emission is tentatively attributed to localized energy states formed at the disordered boundaries of h-BN and graphene. The weak blue emission at the heterojunctions in simple in-plane heterostructures of h-BN and graphene can be enhanced by increasing the density of the interface in graphene quantum dots array embedded in the h-BN monolayer. This work suggests that the narrowest, atomically resolved heterojunctions of in-plane two-dimensional heterostructures provides a future playground for optoelectronics. © 2020, The Author(s).

  • Canted Persistent Spin Texture and Quantum Spin Hall Effect in WTe2

    Garcia J.H., Vila M., Hsu C.-H., Waintal X., Pereira V.M., Roche S. Physical Review Letters; 125 (25, 256603) 2020. 10.1103/PhysRevLett.125.256603. IF: 8.385

    We report an unconventional quantum spin Hall phase in the monolayer WTe2, which exhibits hitherto unknown features in other topological materials. The low symmetry of the structure induces a canted spin texture in the yz plane, which dictates the spin polarization of topologically protected boundary states. Additionally, the spin Hall conductivity gets quantized (2e2/h) with a spin quantization axis parallel to the canting direction. These findings are based on large-scale quantum simulations of the spin Hall conductivity tensor and nonlocal resistances in multiprobe geometries using a realistic tight-binding model elaborated from first-principle methods. The observation of this canted quantum spin Hall effect, related to the formation of topological edge states with nontrivial spin polarization, demands for specific experimental design and suggests interesting alternatives for manipulating spin information in topological materials. © 2020 American Physical Society.

  • Computation of topological phase diagram of disordered Pb1-xSnxTe using the kernel polynomial method

    Dániel Varjas, Michel Fruchart, Anton R. Akhmerov, Pablo M. Perez-Piskunow Physical Review Research; 2 (13229) 2020. 10.1103/PhysRevResearch.2.013229 .

    We present an algorithm to determine topological invariants of inhomogeneous systems, such as alloys, disordered crystals, or amorphous systems. Based on the kernel polynomial method, our algorithm allows us to study samples with more than 107 degrees of freedom. Our method enables the study of large complex compounds, where disorder is inherent to the system. We use it to analyze Pb1−xSnxTe and tighten the critical concentration for the phase transition. Moreover, we obtain the topological phase diagram for related alloys in the family of three-dimensional mirror Chern insulators.

  • Efficient machine-learning based interatomic potentialsfor exploring thermal conductivity in two-dimensional materials

    Bohayra Mortazavi, V Evgeny Podryabinkin, Ivan S Novikov, Stephan Roche, Timon Rabczuk,Xiaoying Zhuang; V Alexander Shapeev Journal of Physics-Materials; 3 (2, 02LT02) 2020. 10.1088/2515-7639/ab7cbb. IF: 0.000

  • Emergence of intraparticle entanglement and time-varying violation of Bell's inequality in Dirac matter

    De Moraes B.G., Cummings A.W., Roche S. Physical Review B; 102 (4, 041403) 2020. 10.1103/PhysRevB.102.041403. IF: 3.575

    We demonstrate the emergence and dynamics of intraparticle entanglement in massless Dirac fermions. This entanglement, generated by spin-orbit coupling, arises between the spin and sublattice pseudospin of electrons in graphene. The entanglement is a complex dynamic quantity but is generally large, independent of the initial state. Its time dependence implies a dynamical violation of a Bell inequality, while its magnitude indicates that large intraparticle entanglement is a general feature of graphene on a substrate. These features are also expected to impact entanglement between pairs of particles, and may be detectable in experiments that combine Cooper pair splitting with nonlocal measurements of spin-spin correlation in mesoscopic devices based on Dirac materials. © 2020 American Physical Society.

  • Exploring event horizons and Hawking radiation through deformed graphene membranes

    Morresi T., Binosi D., Simonucci S., Piergallini R., Roche S., Pugno N.M., Simone T. 2D Materials; 7 (4, 041006) 2020. 10.1088/2053-1583/aba448. IF: 7.140

    Analogue gravitational systems are becoming an increasing popular way of studying the behaviour of quantum systems in curved spacetime. Setups based on ultracold quantum gases in particular, have been recently harnessed to explore the thermal nature of Hawking's and Unruh's radiation that was theoretically predicted almost 50 years ago. For solid state implementations, a promising system is graphene, in which a link between the Dirac-like low-energy electronic excitations and relativistic quantum field theories has been unveiled soon after its discovery. This link could be extended to the case of curved quantum field theory when the graphene sheet is shaped in a surface of constant negative curvature, known as Beltrami's pseudosphere. Here we provide numerical evidence that energetically stable negative curvature graphene surfaces can be realized. Owing to large-scale simulations, our geometrical realizations are characterized by a ratio between the carbon-carbon bond length and the pseudosphere radius small enough to allow the formation of an analog of a black hole event horizon. Additionally, from the energy dependence of the spatially resolved density of states, we infer some thermal properties of the corresponding gravitational system, which could be investigated using low temperature scanning tunnelling microscopy or optical near field spectroscopy. These findings pave the way to the realization of a solid-state system in which the curved spacetime dynamics of quantum many body systems can be investigated. © 2020 IOP Publishing Ltd.

  • Exploring phononic properties of two-dimensional materials using machine learning interatomic potentials

    Mortazavi B., Novikov I.S., Podryabinkin E.V., Roche S., Rabczuk T., Shapeev A.V., Zhuang X. Applied Materials Today; 20 (100685) 2020. 10.1016/j.apmt.2020.100685. IF: 8.352

    Phononic properties are commonly studied by calculating force constants using the density functional theory (DFT) simulations. Although DFT simulations offer accurate estimations of phonon dispersion relations or thermal properties, but for low-symmetry and nanoporous structures the computational cost quickly becomes very demanding. Moreover, the computational setups may yield nonphysical imaginary frequencies in the phonon dispersion curves, impeding the assessment of phononic properties and the dynamical stability of the considered system. Here, we compute phonon dispersion relations and examine the dynamical stability of a large ensemble of novel materials and compositions. We propose a fast and convenient alternative to DFT simulations which derived from machine-learning interatomic potentials passively trained over computationally efficient ab-initio molecular dynamics trajectories. Our results for diverse two-dimensional (2D) nanomaterials confirm that the proposed computational strategy can reproduce fundamental thermal properties in close agreement with those obtained via the DFT approach. The presented method offers a stable, efficient, and convenient solution for the examination of dynamical stability and exploring the phononic properties of low-symmetry and porous 2D materials. © 2020 Elsevier Ltd

  • Impact of oxidation morphology on reduced graphene oxides upon thermal annealing

    Aleandro Antidormi, Stephan Roche, Luciano Colombo Journal Of Physics-Materials; 3 (1, 15011) 2020. 10.1088/2515-7639/ab5ef2. IF: 0.000

  • Impact of synthetic conditions on the anisotropic thermal conductivity of poly(3,4-ethylenedioxythiophene) (PEDOT): A molecular dynamics investigation

    Cappai A., Antidormi A., Bosin A., Narducci D., Colombo L., Melis C. Physical Review Materials; 4 (3, 035401) 2020. 10.1103/PhysRevMaterials.4.035401. IF: 3.337

    In this work we study the effect of different synthetic conditions on thermal transport properties of poly(3,4-ethylenedioxythiophene) (PEDOT) by focusing in particular on the role of proton scavengers. To this aim, different PEDOT samples were generated in silico using a novel computational algorithm based on a combination of first-principles density functional theory and classical molecular dynamics simulations. The corresponding thermal conductivities were then estimated using the approach to equilibrium molecular dynamics methodology. The results show that the initial synthetic conditions strongly affect the corresponding thermal conductivities, which display variations up to a factor of ∼2 depending on the proton scavenger. By decomposing the thermal conductivity tensor along the direction of maximum chain alignment and the corresponding perpendicular directions, we attribute the thermal conductivity differences to the variations in the average polymer chain length λave. A dependence of the thermal conductivity with the polydispersity index was finally observed, suggesting a possible role of intercrystallite chains in enhancing thermal transport properties. By means of the Green-Kubo modal analysis, we eventually characterize the vibrational modes involved in PEDOT thermal transport and investigate how they are related to the thermal conductivity values of the samples. © 2020 American Physical Society.

  • Machine-learning interatomic potentials enable first-principles multiscale modeling of lattice thermal conductivity in graphene/borophene heterostructures

    Mortazavi B., Podryabinkin E.V., Roche S., Rabczuk T., Zhuang X., Shapeev A.V. Materials Horizons; 7 (9): 2359 - 2367. 2020. 10.1039/d0mh00787k. IF: 12.319

    One of the ultimate goals of computational modeling in condensed matter is to be able to accurately compute materials properties with minimal empirical information. First-principles approaches such as density functional theory (DFT) provide the best possible accuracy on electronic properties but they are limited to systems up to a few hundreds, or at most thousands of atoms. On the other hand, classical molecular dynamics (CMD) simulations and the finite element method (FEM) are extensively employed to study larger and more realistic systems, but conversely depend on empirical information. Here, we show that machine-learning interatomic potentials (MLIPs) trained over short ab initio molecular dynamics trajectories enable first-principles multiscale modeling, in which DFT simulations can be hierarchically bridged to efficiently simulate macroscopic structures. As a case study, we analyze the lattice thermal conductivity of coplanar graphene/borophene heterostructures, recently synthesized experimentally (Sci. Adv., 2019, 5, eaax6444), for which no viable classical modeling alternative is presently available. Our MLIP-based approach can efficiently predict the lattice thermal conductivity of graphene and borophene pristine phases, the thermal conductance of complex graphene/borophene interfaces and subsequently enable the study of effective thermal transport along the heterostructures at continuum level. This work highlights that MLIPs can be effectively and conveniently employed to enable first-principles multiscale modeling via hierarchical employment of DFT/CMD/FEM simulations, thus expanding the capability for computational design of novel nanostructures. © 2020 The Royal Society of Chemistry.

  • Magnetic proximity in a van der Waals heterostructure of magnetic insulator and graphene

    Karpiak B., Cummings A.W., Zollner K., Vila M., Khokhriakov D., Hoque A.M., Dankert A., Svedlindh P., Fabian J., Roche S., Dash S.P. 2D Materials; 7 (1, 015026) 2020. 10.1088/2053-1583/ab5915. IF: 7.140

    Engineering 2D material heterostructures by combining the best of different materials in one ultimate unit can offer a plethora of opportunities in condensed matter physics. Here, in the van der Waals heterostructures of the ferromagnetic insulator Cr2Ge2Te6 and graphene, our observations indicate an out-of-plane proximity-induced ferromagnetic exchange interaction in graphene. The perpendicular magnetic anisotropy of Cr2Ge2Te6 results in significant modification of the spin transport and precession in graphene, which can be ascribed to the proximity-induced exchange interaction. Furthermore, the observation of a larger lifetime for perpendicular spins in comparison to the in-plane counterpart suggests the creation of a proximity-induced anisotropic spin texture in graphene. Our experimental results and density functional theory calculations open up opportunities for the realization of proximity-induced magnetic interactions and spin filters in 2D material heterostructures and can form the basic building blocks for future spintronic and topological quantum devices. © 2019 IOP Publishing Ltd.

  • Magnetism, spin dynamics, and quantum transport in two-dimensional systems

    Savero Torres W., Sierra J.F., Benítez L.A., Bonell F., García J.H., Roche S., Valenzuela S.O. MRS Bulletin; 45 (5): 357 - 365. 2020. 10.1557/mrs.2020.121. IF: 5.061

    Two-dimensional (2D) quantum materials offer a unique platform to explore mesoscopic phenomena driven by interfacial and topological effects. Their tunable electric properties and bidimensional nature enable their integration into sophisticated heterostructures with engineered properties, resulting in the emergence of new exotic phenomena not accessible in other platforms. This has fostered many studies on 2D ferromagnetism, proximity-induced effects, and quantum transport, demonstrating their relevance for fundamental research and future device applications. Here, we review ongoing progress in this lively research field with special emphasis on spin-related phenomena. © Materials Research Society 2020.

  • Nonlocal Spin Dynamics in the Crossover from Diffusive to Ballistic Transport

    Vila M., Garcia J.H., Cummings A.W., Power S.R., Groth C.W., Waintal X., Roche S. Physical Review Letters; 124 (19, 196602) 2020. 10.1103/PhysRevLett.124.196602. IF: 8.385

    Improved fabrication techniques have enabled the possibility of ballistic transport and unprecedented spin manipulation in ultraclean graphene devices. Spin transport in graphene is typically probed in a nonlocal spin valve and is analyzed using spin diffusion theory, but this theory is not necessarily applicable when charge transport becomes ballistic or when the spin diffusion length is exceptionally long. Here, we study these regimes by performing quantum simulations of graphene nonlocal spin valves. We find that conventional spin diffusion theory fails to capture the crossover to the ballistic regime as well as the limit of long spin diffusion length. We show that the latter can be described by an extension of the current theoretical framework. Finally, by covering the whole range of spin dynamics, our study opens a new perspective to predict and scrutinize spin transport in graphene and other two-dimensional material-based ultraclean devices. © 2020 American Physical Society.

  • Optimization of the sensitivity of a double-dot magnetic detector

    Macucci M., Marconcini P., Roche S. Electronics (Switzerland); 9 (7, 1134): 1 - 13. 2020. 10.3390/electronics9071134. IF: 2.412

    We investigate, by means of numerical simulations, the lowest magnetic field level that can be detected with a given relative accuracy with a sensor based on a double-dot device fabricated in a high-mobility two-dimensional electron gas. The double dot consists of a cavity delimited by an input and an output constriction, with a potential barrier exactly in the middle. In conditions of perfect symmetry, a strong conductance enhancement effect appears as a consequence of the constructive interference between symmetric trajectories. When the symmetry is broken, for example by the presence of an applied magnetic field, this enhancement effect is suppressed. We explore the design parameter space and assess the minimum magnetic field value that can be measured with a given accuracy in the presence of flicker noise. © 2020 by the authors. Licensee MDPI, Basel, Switzerland.

  • Tunable room-temperature spin galvanic and spin Hall effects in van der Waals heterostructures

    Benítez L.A., Savero Torres W., Sierra J.F., Timmermans M., Garcia J.H., Roche S., Costache M.V., Valenzuela S.O. Nature Materials; 19 (2): 170 - 175. 2020. 10.1038/s41563-019-0575-1. IF: 38.663

    Spin–orbit coupling stands as a powerful tool to interconvert charge and spin currents and to manipulate the magnetization of magnetic materials through spin-torque phenomena. However, despite the diversity of existing bulk materials and the recent advent of interfacial and low-dimensional effects, control of this interconversion at room temperature remains elusive. Here, we demonstrate strongly enhanced room-temperature spin-to-charge interconversion in graphene driven by the proximity of WS2. By performing spin precession experiments in appropriately designed Hall bars, we separate the contributions of the spin Hall and the spin galvanic effects. Remarkably, their corresponding conversion efficiencies can be tailored by electrostatic gating in magnitude and sign, peaking near the charge neutrality point with an equivalent magnitude that is comparable to the largest efficiencies reported to date. Such electric-field tunability provides a building block for spin generation free from magnetic materials and for ultra-compact magnetic memory technologies. © 2020, The Author(s), under exclusive licence to Springer Nature Limited.

  • Ultralow-dielectric-constant amorphous boron nitride

    Hong S., Lee C.-S., Lee M.-H., Lee Y., Ma K.Y., Kim G., Yoon S.I., Ihm K., Kim K.-J., Shin T.J., Kim S.W., Jeon E.-C., Jeon H., Kim J.-Y., Lee H.-I., Lee Z., Antidormi A., Roche S., Chhowalla M., Shin H.-J., Shin H.S. Nature; 582 (7813): 511 - 514. 2020. 10.1038/s41586-020-2375-9. IF: 42.779

    Decrease in processing speed due to increased resistance and capacitance delay is a major obstacle for the down-scaling of electronics1–3. Minimizing the dimensions of interconnects (metal wires that connect different electronic components on a chip) is crucial for the miniaturization of devices. Interconnects are isolated from each other by non-conducting (dielectric) layers. So far, research has mostly focused on decreasing the resistance of scaled interconnects because integration of dielectrics using low-temperature deposition processes compatible with complementary metal–oxide–semiconductors is technically challenging. Interconnect isolation materials must have low relative dielectric constants (κ values), serve as diffusion barriers against the migration of metal into semiconductors, and be thermally, chemically and mechanically stable. Specifically, the International Roadmap for Devices and Systems recommends4 the development of dielectrics with κ values of less than 2 by 2028. Existing low-κ materials (such as silicon oxide derivatives, organic compounds and aerogels) have κ values greater than 2 and poor thermo-mechanical properties5. Here we report three-nanometre-thick amorphous boron nitride films with ultralow κ values of 1.78 and 1.16 (close to that of air, κ = 1) at operation frequencies of 100 kilohertz and 1 megahertz, respectively. The films are mechanically and electrically robust, with a breakdown strength of 7.3 megavolts per centimetre, which exceeds requirements. Cross-sectional imaging reveals that amorphous boron nitride prevents the diffusion of cobalt atoms into silicon under very harsh conditions, in contrast to reference barriers. Our results demonstrate that amorphous boron nitride has excellent low-κ dielectric characteristics for high-performance electronics. © 2020, The Author(s), under exclusive licence to Springer Nature Limited.