-
Exchange interactions and intermolecular hybridization in a spin-1/2 nanographene dimer
Authors:
N. Krane,
E. Turco,
A. Bernhardt,
D. Jacob,
G. Gandus,
D. Passerone,
M. Luisier,
M. Juríček,
R. Fasel,
J. Fernández-Rossier,
P. Ruffieux
Abstract:
Phenalenyl is a radical nanographene with triangular shape that hosts an unpaired electron with spin S = 1/2. The open-shell nature of phenalenyl is expected to be retained in covalently bonded networks. Here, we study a first step in that direction and report the synthesis of the phenalenyl dimer by combining in-solution synthesis and on-surface activation and its characterization both on Au(111)…
▽ More
Phenalenyl is a radical nanographene with triangular shape that hosts an unpaired electron with spin S = 1/2. The open-shell nature of phenalenyl is expected to be retained in covalently bonded networks. Here, we study a first step in that direction and report the synthesis of the phenalenyl dimer by combining in-solution synthesis and on-surface activation and its characterization both on Au(111) and on a monolayer of NaCl on top of Au(111) by means of inelastic electron tunneling spectroscopy (IETS). IETS shows inelastic steps that, together with a thorough theoretical analysis, are identified as the singlet-triplet excitation arising from interphenalenyl exchange. Two prominent features of our data permit to shed light on the nature of spin interactions in this system. First, the excitation energies with and without the NaCl decoupling layer are 48 and 41 meV, respectively, indicating a significant renormalization of the spin excitation energies due to exchange with the Au(111) electrons. Second, a position-dependent bias-asymmetry of the height of the inelastic steps is accounted for by an interphenalenyl hybridization of the singly occupied phenalenyl orbitals that is only possible via third neighbor hopping. This hybridization is also essential to activate kinetic interphenalenyl exchange. Our results set the stage for future work on the bottom-up synthesis of spin S = 1/2 spin lattices with large exchange interaction.
△ Less
Submitted 19 July, 2023;
originally announced July 2023.
-
Platinum contacts for 9-atom-wide armchair graphene nanoribbons
Authors:
Chunwei Hsu,
Michael Rohde,
Gabriela Borin Barin,
Guido Gandus,
Daniele Passerone,
Mathieu Luisier,
Pascal Ruffieux,
Roman Fasel,
Herre S. J. van der Zant,
Maria El Abbassi
Abstract:
Creating a good contact between electrodes and graphene nanoribbons (GNRs) has been a longstanding challenge in searching for the next GNR-based nanoelectronics. This quest requires the controlled fabrication of sub-20 nm metallic gaps, a clean GNR transfer minimizing damage and organic contamination during the device fabrication, as well as work function matching to minimize the contact resistanc…
▽ More
Creating a good contact between electrodes and graphene nanoribbons (GNRs) has been a longstanding challenge in searching for the next GNR-based nanoelectronics. This quest requires the controlled fabrication of sub-20 nm metallic gaps, a clean GNR transfer minimizing damage and organic contamination during the device fabrication, as well as work function matching to minimize the contact resistance. Here, we transfer 9-atom-wide armchair-edged GNRs (9-AGNRs) grown on Au(111)/mica substrates to pre-patterned platinum electrodes, yielding polymer-free 9-AGNR field-effect transistor devices. Our devices have a resistance in the range of $10^6$ to $10^8$ $Ω$ in the low-bias regime, which is 2 to 4 orders of magnitude lower than previous reports. Density functional theory (DFT) calculations combined with the non-equilibrium Green's function method (NEGF) explain the observed p-type electrical characteristics and further demonstrate that platinum gives strong coupling and higher transmission in comparison to other materials such as graphene.
△ Less
Submitted 31 January, 2023;
originally announced January 2023.
-
Strongly correlated physics in organic open-shell quantum systems
Authors:
G. Gandus,
D. Passerone,
R. Stadler,
M. Luisier,
A. Valli
Abstract:
Strongly correlated physics arises due to electron-electron scattering within partially-filled orbitals, and in this perspective, organic molecules in open-shell configuration are good candidates to exhibit many-body effects. With a focus on neutral organic radicals with a molecular orbital hosting a single unpaired electron (SOMO) we investigate many-body effects on electron transport in a single…
▽ More
Strongly correlated physics arises due to electron-electron scattering within partially-filled orbitals, and in this perspective, organic molecules in open-shell configuration are good candidates to exhibit many-body effects. With a focus on neutral organic radicals with a molecular orbital hosting a single unpaired electron (SOMO) we investigate many-body effects on electron transport in a single-molecule junction setup. Within a combination of density functional theory and many-body techniques, we perform numerical simulations for an effective model for which all the parameters, including the Coulomb tensor, are derived ab-initio. We demonstrate that the SOMO resonance is prone towards splitting, and identify a giant electronic scattering rate as the driving many-body mechanism, akin to a Mott metal-to-insulator transition. The nature of the splitting, and thus of the resulting gap, as well as the spatial distribution of the SOMO and its coupling to the electrodes, have dramatic effects on the transport properties of the junction. We argue that the phenomenon and the underlying microscopic mechanism are general, and apply to a wide family of open-shell molecular systems.
△ Less
Submitted 31 December, 2022;
originally announced January 2023.
-
Efficient and accurate defect level modelling in monolayer MoS$_2$ via GW+DFT with open boundary conditions
Authors:
Guido Gandus,
Youseung Lee,
Leonard Deuschle,
Mathieu Luisier,
Daniele Passerone
Abstract:
Within the framework of many-body perturbation theory (MBPT) integrated with density functional theory (DFT), a novel defect-subspace projection GW method, the so-called p-GW, is proposed. By avoiding the periodic defect interference through open boundary self-energies, we show that the p-GW can efficiently and accurately describe quasi-particle correlated defect levels in two-dimensional (2D) mon…
▽ More
Within the framework of many-body perturbation theory (MBPT) integrated with density functional theory (DFT), a novel defect-subspace projection GW method, the so-called p-GW, is proposed. By avoiding the periodic defect interference through open boundary self-energies, we show that the p-GW can efficiently and accurately describe quasi-particle correlated defect levels in two-dimensional (2D) monolayer MoS$_2$. By comparing two different defect states originating from sulfur vacancy and adatom to existing theoretical and experimental works, we show that our GW correction to the DFT defect levels is precisely modelled. Based on these findings, we expect that our method can provide genuine trap states for various 2D transition-metal dichalcogenide (TMD) monolayers, thus enabling the study of defect-induced effects on the device characteristics of these materials via realistic simulations.
△ Less
Submitted 29 June, 2022;
originally announced June 2022.
-
Efficient partitioning of surface Green's function: toward ab initio contact resistance study
Authors:
Guido Gandus,
Youseung Lee,
Daniele Passerone,
Mathieu Luisier
Abstract:
In this work, we propose an efficient computational scheme for first-principle quantum transport simulations to evaluate the open-boundary conditions. Its partitioning differentiates from conventional methods in that the contact self-energy matrices are constructed on smaller building blocks, principal layers (PL), while conventionally it was restricted to have the same lateral dimensions of the a…
▽ More
In this work, we propose an efficient computational scheme for first-principle quantum transport simulations to evaluate the open-boundary conditions. Its partitioning differentiates from conventional methods in that the contact self-energy matrices are constructed on smaller building blocks, principal layers (PL), while conventionally it was restricted to have the same lateral dimensions of the adjoining atoms in a channel region. Here, we obtain the properties of bulk electrodes through non-equilibrium Green's function (NEGF) approach with significant improvements in the computational efficiency without sacrificing the accuracy of results. To exemplify the merits of the proposed method we investigate the carrier density dependency of contact resistances in silicon nanowire devices connected to bulk metallic contacts.
△ Less
Submitted 29 December, 2020;
originally announced December 2020.
-
Smart local orbitals for efficient calculations within density functional theory and beyond
Authors:
G. Gandus,
A. Valli,
D. Passerone,
R. Stadler
Abstract:
Localized basis sets in the projector augmented wave formalism allow for computationally efficient calculations within density functional theory (DFT). However, achieving high numerical accuracy requires an extensive basis set, which also poses a fundamental problem for the interpretation of the results. We present a way to obtain a reduced basis set of atomic orbitals through the subdiagonalizati…
▽ More
Localized basis sets in the projector augmented wave formalism allow for computationally efficient calculations within density functional theory (DFT). However, achieving high numerical accuracy requires an extensive basis set, which also poses a fundamental problem for the interpretation of the results. We present a way to obtain a reduced basis set of atomic orbitals through the subdiagonalization of each atomic block of the Hamiltonian. The resulting local orbitals (LOs) inherit the information of the local crystal field. In the LO basis, it becomes apparent that the Hamiltonian is nearly block-diagonal, and we demonstrate that it is possible to keep only a subset of relevant LOs which provide an accurate description of the physics around the Fermi level. This reduces to some extent the redundancy of the original basis set, and at the same time it allows one to perform post-processing of DFT calculations, ranging from the interpretation of electron transport to extracting effective tight-binding Hamiltonians, very efficiently and without sacrificing the accuracy of the results.
△ Less
Submitted 17 November, 2020;
originally announced November 2020.
-
AiiDAlab -- an ecosystem for developing, executing, and sharing scientific workflows
Authors:
Aliaksandr V. Yakutovich,
Kristjan Eimre,
Ole Schütt,
Leopold Talirz,
Carl S. Adorf,
Casper W. Andersen,
Edward Ditler,
Dou Du,
Daniele Passerone,
Berend Smit,
Nicola Marzari,
Giovanni Pizzi,
Carlo A. Pignedoli
Abstract:
Cloud platforms allow users to execute tasks directly from their web browser and are a key enabling technology not only for commerce but also for computational science. Research software is often developed by scientists with limited experience in (and time for) user interface design, which can make research software difficult to install and use for novices. When combined with the increasing comple…
▽ More
Cloud platforms allow users to execute tasks directly from their web browser and are a key enabling technology not only for commerce but also for computational science. Research software is often developed by scientists with limited experience in (and time for) user interface design, which can make research software difficult to install and use for novices. When combined with the increasing complexity of scientific workflows (involving many steps and software packages), setting up a computational research environment becomes a major entry barrier. AiiDAlab is a web platform that enables computational scientists to package scientific workflows and computational environments and share them with their collaborators and peers. By leveraging the AiiDA workflow manager and its plugin ecosystem, developers get access to a growing range of simulation codes through a python API, coupled with automatic provenance tracking of simulations for full reproducibility. Computational workflows can be bundled together with user-friendly graphical interfaces and made available through the AiiDAlab app store. Being fully compatible with open-science principles, AiiDAlab provides a complete infrastructure for automated workflows and provenance tracking, where incorporating new capabilities becomes intuitive, requiring only Python knowledge.
△ Less
Submitted 29 September, 2020;
originally announced October 2020.
-
Mapping the Structure of Oxygen-Doped Wurtzite Aluminum Nitride Coatings From Ab Initio Random Structure Search and Experiments
Authors:
Piero Gasparotto,
Maria Fischer,
Daniele Scopece,
Maciej Oskar Liedke,
Maik Butterling,
Andreas Wagner,
Oguz Yildirim,
Mathis Trant,
Daniele Passerone,
Hans J. Hug,
Carlo Antonio Pignedoli
Abstract:
Machine learning is changing how we design and interpret experiments in materials science. In this work, we show how unsupervised learning, combined with ab initio modeling, improves our understanding of structural metastability in multicomponent alloys. We use the example case of Al-O-N alloys where the formation of aluminum vacancies in wurtzite AlN upon the incorporation of substitutional oxyge…
▽ More
Machine learning is changing how we design and interpret experiments in materials science. In this work, we show how unsupervised learning, combined with ab initio modeling, improves our understanding of structural metastability in multicomponent alloys. We use the example case of Al-O-N alloys where the formation of aluminum vacancies in wurtzite AlN upon the incorporation of substitutional oxygen can be seen as a general mechanism of solids where crystal symmetry is reduced to stabilize defects. The ideal AlN wurtzite crystal structure occupation cannot be matched due to the presence of an aliovalent hetero-element into the structure. The traditional interpretation of the c-lattice shrinkage in sputter-deposited Al-O-N films from X-ray diffraction (XRD) experiments suggests the existence of a solubility limit at 8at.% oxygen content. Here we show that such naive interpretation is misleading. We support XRD data with a machine learning analysis of ab initio simulations and positron annihilation lifetime spectroscopy data, revealing no signs of a possible solubility limit. Instead, the presence of a wide range of non-equilibrium oxygen-rich defective structures emerging at increasing oxygen contents suggests that the formation of grain boundaries is the most plausible mechanism responsible for the lattice shrinkage measured in Al-O-N sputtered films.
△ Less
Submitted 12 January, 2021; v1 submitted 28 September, 2020;
originally announced September 2020.
-
Metallic Carbon Nanotube Quantum Dots with Broken Symmetries as a Platform for Tunable Terahertz Detection
Authors:
G. Buchs,
M. Marganska,
J. W. González,
K. Eimre,
C. A. Pignedoli,
D. Passerone,
A. Ayuela,
O. Gröning,
D. Bercioux
Abstract:
Generating and detecting radiation in the technologically relevant range of the so-called terahertz gap ($0.1 - 10$ THz) is challenging because of a lack of efficient sources and detectors. Quantum dots in carbon nanotubes have shown great potential to build sensitive terahertz detectors usually based on photon-assisted tunnelling. A recently reported mechanism combining resonant quantum dot trans…
▽ More
Generating and detecting radiation in the technologically relevant range of the so-called terahertz gap ($0.1 - 10$ THz) is challenging because of a lack of efficient sources and detectors. Quantum dots in carbon nanotubes have shown great potential to build sensitive terahertz detectors usually based on photon-assisted tunnelling. A recently reported mechanism combining resonant quantum dot transitions and tunnelling barriers asymmetries results in a narrow linewidth photocurrent response with a large signal-to-noise ratio under weak THz radiation. That device was sensitive to one frequency, corresponding to transitions between equidistant quantized states. In this work we show, using numerical together with scanning tunnelling spectroscopy studies of a defect-induced metallic zigzag single-walled carbon nanotube quantum dot that simultaneously breaking various symmetries in metallic nanotube quantum dots of arbitrary chirality strongly relaxes the selection rules in the electric dipole approximation, and removes energy degeneracies. This leads to a richer set of allowed optical transitions spanning frequencies from 1 THz to several tens of THz, for a $\sim$10 nm quantum dot. Based on these findings, we propose a terahertz detector device based on a metallic single-walled carbon nanotube quantum dot defined by artificial defects. Depending on its length and contacts transparency, the operating regimes range from a high-resolution gate-tunable terahertz sensor to a broadband terahertz detector. Our calculations indicate that the device is largely unaffected by temperatures up to 100 K, making carbon nanotube quantum dots with broken symmetries a promising platform to design tunable terahertz detectors that could operate at liquid nitrogen temperatures.
△ Less
Submitted 14 May, 2021; v1 submitted 18 June, 2020;
originally announced June 2020.
-
Electronic structure of pristine and Ni-substituted LaFeO$_3$ from near edge x-ray absorption fine structure experiments and first-principles simulations
Authors:
Iurii Timrov,
Piyush Agrawal,
Xinyu Zhang,
Selma Erat,
Riping Liu,
Artur Braun,
Matteo Cococcioni,
Matteo Calandra,
Nicola Marzari,
Daniele Passerone
Abstract:
We present a joint theoretical and experimental study of the oxygen $K$-edge spectra for LaFeO$_3$ and homovalent Ni-substituted LaFeO$_3$ (LaFe$_{0.75}$Ni$_{0.25}$O$_3$), using first-principles simulations based on density-functional theory with extended Hubbard functionals and x-ray absorption near edge structure (XANES) measurements. Ground-state and excited-state XANES calculations employ Hubb…
▽ More
We present a joint theoretical and experimental study of the oxygen $K$-edge spectra for LaFeO$_3$ and homovalent Ni-substituted LaFeO$_3$ (LaFe$_{0.75}$Ni$_{0.25}$O$_3$), using first-principles simulations based on density-functional theory with extended Hubbard functionals and x-ray absorption near edge structure (XANES) measurements. Ground-state and excited-state XANES calculations employ Hubbard on-site $U$ and inter-site $V$ parameters determined from first principles and the Lanczos recursive method to obtain absorption cross sections, which allows for a reliable description of XANES spectra in transition-metal compounds in a very broad energy range, with an accuracy comparable to that of hybrid functionals but at a substantially lower cost. We show that standard gradient-corrected exchange-correlation functionals fail in capturing accurately the electronic properties of both materials. In particular, for LaFe$_{0.75}$Ni$_{0.25}$O$_3$ they do not reproduce its semiconducting behaviour and provide a poor description of the pre-edge features at the O $K$ edge. The inclusion of Hubbard interactions leads to a drastic improvement, accounting for the semiconducting ground state of LaFe$_{0.75}$Ni$_{0.25}$O$_3$ and for a good agreement between calculated and measured XANES spectra. We show that the partial substitution of Fe for Ni affects the conduction-band bottom by generating a strongly hybridized O($2p$)-Ni($3d$) minority-spin empty electronic state. The present work, based on a consistent correction of self-interaction errors, outlines the crucial role of extended Hubbard functionals to describe the electronic structure of complex transition-metal oxides such as LaFeO$_3$ and LaFe$_{0.75}$Ni$_{0.25}$O$_3$ and paves the way to future studies on similar systems.
△ Less
Submitted 24 August, 2020; v1 submitted 8 April, 2020;
originally announced April 2020.
-
Materials Cloud, a platform for open computational science
Authors:
Leopold Talirz,
Snehal Kumbhar,
Elsa Passaro,
Aliaksandr V. Yakutovich,
Valeria Granata,
Fernando Gargiulo,
Marco Borelli,
Martin Uhrin,
Sebastiaan P. Huber,
Spyros Zoupanos,
Carl S. Adorf,
Casper W. Andersen,
Ole Schütt,
Carlo A. Pignedoli,
Daniele Passerone,
Joost VandeVondele,
Thomas C. Schulthess,
Berend Smit,
Giovanni Pizzi,
Nicola Marzari
Abstract:
Materials Cloud is a platform designed to enable open and seamless sharing of resources for computational science, driven by applications in materials modelling. It hosts 1) archival and dissemination services for raw and curated data, together with their provenance graph, 2) modelling services and virtual machines, 3) tools for data analytics, and pre-/post-processing, and 4) educational material…
▽ More
Materials Cloud is a platform designed to enable open and seamless sharing of resources for computational science, driven by applications in materials modelling. It hosts 1) archival and dissemination services for raw and curated data, together with their provenance graph, 2) modelling services and virtual machines, 3) tools for data analytics, and pre-/post-processing, and 4) educational materials. Data is citable and archived persistently, providing a comprehensive embodiment of the FAIR principles that extends to computational workflows. Materials Cloud leverages the AiiDA framework to record the provenance of entire simulation pipelines (calculations performed, codes used, data generated) in the form of graphs that allow to retrace and reproduce any computed result. When an AiiDA database is shared on Materials Cloud, peers can browse the interconnected record of simulations, download individual files or the full database, and start their research from the results of the original authors. The infrastructure is agnostic to the specific simulation codes used and can support diverse applications in computational science that transcend its initial materials domain.
△ Less
Submitted 27 March, 2020;
originally announced March 2020.
-
Band Gap of Atomically Precise Graphene Nanoribbons as a Function of Ribbon Length and Termination
Authors:
Leopold Talirz,
Hajo Söde,
Shigeki Kawai,
Pascal Ruffieux,
Ernst Meyer,
Xinliang Feng,
Klaus Müllen,
Roman Fasel,
Carlo A. Pignedoli,
Daniele Passerone
Abstract:
We study the band gap of finite $N_A=7$ armchair graphene nanoribbons (7-AGNRs) on Au(111) through scanning tunneling microscopy/spectroscopy combined with density functional theory calculations. The band gap of 7-AGNRs with lengths of 6 nm and more is converged to within 0.1 eV of its bulk value of 2.3 eV, while the band gap opens by several hundred meV in very short 7-AGNRs. The termination has…
▽ More
We study the band gap of finite $N_A=7$ armchair graphene nanoribbons (7-AGNRs) on Au(111) through scanning tunneling microscopy/spectroscopy combined with density functional theory calculations. The band gap of 7-AGNRs with lengths of 6 nm and more is converged to within 0.1 eV of its bulk value of 2.3 eV, while the band gap opens by several hundred meV in very short 7-AGNRs. The termination has a significant effect on the band gap, doubly hydrogenated termini yielding a lower band gap than singly hydrogenated ones.
△ Less
Submitted 3 October, 2019;
originally announced October 2019.
-
Superlubricity of Graphene Nanoribbons on Gold Surfaces
Authors:
Shigeki Kawai,
Andrea Benassi,
Enrico Gnecco,
Hajo Soede,
Remy Pawlak,
Xinliang Feng,
Klaus Muellen,
Daniele Passerone,
Carlo A. Pignedoli,
Pascal Ruffieux,
Roman Fasel,
Ernst Meyer
Abstract:
The state of vanishing friction known as superlubricity has important applications for energy saving and increasing the lifetime of devices. Superlubricity detected with atomic force microscopy appears in examples like sliding large graphite flakes or gold nanoclusters across surfaces. However, the origin of the behavior is poorly understood due to the lack of a controllable nano-contact. We demon…
▽ More
The state of vanishing friction known as superlubricity has important applications for energy saving and increasing the lifetime of devices. Superlubricity detected with atomic force microscopy appears in examples like sliding large graphite flakes or gold nanoclusters across surfaces. However, the origin of the behavior is poorly understood due to the lack of a controllable nano-contact. We demonstrate graphene nanoribbons superlubricity when sliding on gold with a joint experimental and computational approach. The atomically well-defined contact allows us to trace the origin of superlubricity, unravelling the role played by edges, surface reconstruction and ribbon elasticity. Our results pave the way to the scale-up of superlubricity toward the realization of frictionless coatings.
△ Less
Submitted 13 March, 2016;
originally announced March 2016.
-
On-surface synthesis of graphene nanoribbons with zigzag edge topology
Authors:
Pascal Ruffieux,
Shiyong Wang,
Bo Yang,
Carlos Sanchez,
Jia Liu,
Thomas Dienel,
Leopold Talirz,
Prashant Shinde,
Carlo A. Pignedoli,
Daniele Passerone,
Tim Dumslaff,
Xinliang Feng,
Klaus Muellen,
Roman Fasel
Abstract:
Graphene-based nanostructures exhibit a vast range of exciting electronic properties that are absent in extended graphene. For example, quantum confinement in carbon nanotubes and armchair graphene nanoribbons (AGNRs) leads to the opening of substantial electronic band gaps that are directly linked to their structural boundary conditions. Even more intriguing are nanostructures with zigzag edges,…
▽ More
Graphene-based nanostructures exhibit a vast range of exciting electronic properties that are absent in extended graphene. For example, quantum confinement in carbon nanotubes and armchair graphene nanoribbons (AGNRs) leads to the opening of substantial electronic band gaps that are directly linked to their structural boundary conditions. Even more intriguing are nanostructures with zigzag edges, which are expected to host spin-polarized electronic edge states and can thus serve as key elements for graphene-based spintronics. The most prominent example is zigzag graphene nanoribbons (ZGNRs) for which the edge states are predicted to couple ferromagnetically along the edge and antiferromagnetically between them. So far, a direct observation of the spin-polarized edge states for specifically designed and controlled zigzag edge topologies has not been achieved. This is mainly due to the limited precision of current top-down approaches, which results in poorly defined edge structures. Bottom-up fabrication approaches, on the other hand, were so far only successfully applied to the growth of AGNRs and related structures. Here, we describe the successful bottom-up synthesis of ZGNRs, which are fabricated by the surface-assisted colligation and cyclodehydrogenation of specifically designed precursor monomers including carbon groups that yield atomically precise zigzag edges. Using scanning tunnelling spectroscopy we prove the existence of edge-localized states with large energy splittings. We expect that the availability of ZGNRs will finally allow the characterization of their predicted spin-related properties such as spin confinement and filtering, and ultimately add the spin degree of freedom to graphene-based circuitry.
△ Less
Submitted 16 November, 2015;
originally announced November 2015.
-
Does Rotational Melting Make Molecular Crystal Surfaces More Slippery?
Authors:
Andrea Benassi,
Andrea Vanossi,
Carlo A. Pignedoli,
Daniele Passerone,
Erio Tosatti
Abstract:
The surface of a crystal made of roughly spherical molecules exposes, above its bulk rotational phase transition at T= T$_r$, a carpet of freely rotating molecules, possibly functioning as "nanobearings" in sliding friction. We explored by extensive molecular dynamics simulations the frictional and adhesion changes experienced by a sliding C$_{60}$ flake on the surface of the prototype system C…
▽ More
The surface of a crystal made of roughly spherical molecules exposes, above its bulk rotational phase transition at T= T$_r$, a carpet of freely rotating molecules, possibly functioning as "nanobearings" in sliding friction. We explored by extensive molecular dynamics simulations the frictional and adhesion changes experienced by a sliding C$_{60}$ flake on the surface of the prototype system C$_{60}$ fullerite. At fixed flake orientation both quantities exhibit only a modest frictional drop of order 20% across the transition. However, adhesion and friction drop by a factor of $\sim$ 2 as the flake breaks its perfect angular alignment with the C$_{60}$ surface lattice suggesting an entropy-driven aligned-misaligned switch during pull-off at T$_r$. The results can be of relevance for sliding Kr islands, where very little frictional differences were observed at T$_r$, but also to the sliding of C$_{60}$ -coated tip, where a remarkable factor $\sim$ 2 drop has been reported.
△ Less
Submitted 10 October, 2014;
originally announced October 2014.
-
Microscale Motion Control Through Ferromagnetic Films
Authors:
Andrea Benassi,
Johannes Schwenk,
Miguel A. Marioni,
Hans J. Hug,
Daniele Passerone
Abstract:
Actuation and control of motion in micro-mechanical systems are technological challenges, since they are accompanied by mechanical friction and wear, principal and well known sources of device lifetime reduction. In this theoretical work we propose a non-contact motion control technique based on the introduction of a tunable magnetic interaction. The latter is realized by coating two non-touching…
▽ More
Actuation and control of motion in micro-mechanical systems are technological challenges, since they are accompanied by mechanical friction and wear, principal and well known sources of device lifetime reduction. In this theoretical work we propose a non-contact motion control technique based on the introduction of a tunable magnetic interaction. The latter is realized by coating two non-touching sliding bodies with ferromagnetic films. The resulting dynamics is determined by shape, size and ordering of magnetic domains arising in the films below the Curie temperature. We demonstrate that the domain behavior can be tailored by acting on handles like ferromagnetic coating preparation, external magnetic fields and the finite distance between the plates. In this way, motion control can be achieved without mechanical contact. Moreover, we discuss how such handles can disclose a variety of sliding regimes. Finally, we propose how to practically implement the proposed model sliding system.
△ Less
Submitted 22 July, 2014;
originally announced July 2014.
-
Role of interface coupling inhomogeneity in domain evolution in exchange bias
Authors:
Andrea Benassi,
Miguel A. Marioni,
Daniele Passerone,
Hans J. Hug
Abstract:
Models of exchange-bias in thin films have been able to describe various aspects of this technologically relevant effect. Through appropriate choices of free parameters the modelled hysteresis loops adequately match experiment, and typical domain structures can be simulated. However, the use of these parameters, notably the coupling strength between the systems' ferromagnetic (F) and antiferromagn…
▽ More
Models of exchange-bias in thin films have been able to describe various aspects of this technologically relevant effect. Through appropriate choices of free parameters the modelled hysteresis loops adequately match experiment, and typical domain structures can be simulated. However, the use of these parameters, notably the coupling strength between the systems' ferromagnetic (F) and antiferromagnetic (AF) layers, obscures conclusions about their influence on the magnetization reversal processes. Here we develop a 2D phase-field model of the magnetization process in exchange-biased CoO/(Co/Pt)xn that incorporates the 10 nm-resolved measured local biasing characteristics of the antiferromagnet. Just three interrelated parameters set to measured physical quantities of the ferromagnet and the measured density of uncompensated spins thus suffice to match the experiment in microscopic and macroscopic detail. We use the model to study changes in bias and coercivity caused by different distributions of pinned uncompensated spins of the antiferromagnet, in application-relevant situations where domain wall motion dominates the ferromagnetic reversal. We show the excess coercivity can arise solely from inhomogeneity in the density of biasing- and anti-biasing pinned uncompensated spins in the antiferromagnet. Counter to conventional wisdom, irreversible processes in the latter are not essential.
△ Less
Submitted 29 March, 2014;
originally announced April 2014.
-
Stable ferromagnetism and doping induced half-metallicity in asymmetric graphene nanoribbons
Authors:
Donat J. Adams,
Oliver Gröning,
Carlo A. Pignedoli,
Pascal Ruffieux,
Roman Fasel,
Daniele Passerone
Abstract:
We propose a class of graphene nanoribbons showing strong intrinsic ferromagnetic behavior due to their asymmetry. Such ribbons are based on a zig-zag edged backbone surmounted by a periodic, triangular notched region of variable size. The electronic properties as a function of the topology are investigated. Interestingly, substitutional doping by boron or nitrogen induces half-metallicity. The mo…
▽ More
We propose a class of graphene nanoribbons showing strong intrinsic ferromagnetic behavior due to their asymmetry. Such ribbons are based on a zig-zag edged backbone surmounted by a periodic, triangular notched region of variable size. The electronic properties as a function of the topology are investigated. Interestingly, substitutional doping by boron or nitrogen induces half-metallicity. The most effective doping sites can be inferred from the band structure. Given the present rapid development of bottom-up strategies for the synthesis of atomically precise carbon nanostructures the proposed class of nanoribbons emerges as a real candidate for spintronic applications at ambient temperature.
△ Less
Submitted 23 January, 2012;
originally announced January 2012.
-
A Concerted Variational Strategy for Investigating Rare Events
Authors:
Daniele Passerone,
Matteo Ceccarelli,
Michele Parrinello
Abstract:
A strategy for finding transition paths connecting two stable basins is presented. The starting point is the Hamilton principle of stationary action; we show how it can be transformed into a minimum principle through the addition of suitable constraints like energy conservation. Methods for improving the quality of the paths are presented: for example, the Maupertuis principle can be used for de…
▽ More
A strategy for finding transition paths connecting two stable basins is presented. The starting point is the Hamilton principle of stationary action; we show how it can be transformed into a minimum principle through the addition of suitable constraints like energy conservation. Methods for improving the quality of the paths are presented: for example, the Maupertuis principle can be used for determining the transition time of the trajectory and for coming closer to the desired dynamic path. A saddle point algorithm (conjugate residual method) is shown to be efficient for reaching a ``true'' solution of the original variational problem.
△ Less
Submitted 30 December, 2002;
originally announced December 2002.
-
Bending strain-driven modification of surface reconstructions: Au(111)
Authors:
U. Tartaglino,
E. Tosatti,
D. Passerone,
F. Ercolessi
Abstract:
Strain can affect the morphology of a crystal surface, and cause modifications of its reconstruction even when weak, as in the case of mechanical bending. We carried out calculations of strain-dependent surface free energy and direct bending simulations demonstrating the change of incommensurate reconstruction in Au(111) under strain, in good agreement with recent data. Time-dependent strain sho…
▽ More
Strain can affect the morphology of a crystal surface, and cause modifications of its reconstruction even when weak, as in the case of mechanical bending. We carried out calculations of strain-dependent surface free energy and direct bending simulations demonstrating the change of incommensurate reconstruction in Au(111) under strain, in good agreement with recent data. Time-dependent strain should cause a sliding of the topmost layer over the second, suggesting an interesting case of nanofriction. Bending strain could also be used to fine tune the spacing of selectively absorbed nanoclusters.
△ Less
Submitted 17 July, 2002; v1 submitted 27 October, 2001;
originally announced October 2001.
-
Action-derived molecular dynamics in the study of rare events
Authors:
D. Passerone,
M. Parrinello
Abstract:
We present a practical method to generate classical trajectories with fixed initial and final boundary conditions. Our method is based on the minimization of a suitably defined discretized action. The method finds its most natural application in the study of rare events. Its capabilities are illustrated by non-trivial examples. The algorithm lends itself to straightforward parallelization, and w…
▽ More
We present a practical method to generate classical trajectories with fixed initial and final boundary conditions. Our method is based on the minimization of a suitably defined discretized action. The method finds its most natural application in the study of rare events. Its capabilities are illustrated by non-trivial examples. The algorithm lends itself to straightforward parallelization, and when combined with molecular dynamics (MD) it promises to offer a powerful tool for the study of chemical reactions.
△ Less
Submitted 11 April, 2001;
originally announced April 2001.
-
Two-Stage Rotational Disordering of a Molecular Crystal Surface: C60
Authors:
C. Laforge,
D. Passerone,
A. B. Harris,
Ph. Lambin,
E. Tosatti
Abstract:
We propose a two-stage mechanism for the rotational surface disordering phase transition of a molecular crystal, as realized in C$_{60}$ fullerite. Our study, based on Monte Carlo simulations, uncovers the existence of a new intermediate regime, between a low temperature ordered $(2 \times 2)$ state, and a high temperature $(1 \times 1)$ disordered phase. In the intermediate regime there is part…
▽ More
We propose a two-stage mechanism for the rotational surface disordering phase transition of a molecular crystal, as realized in C$_{60}$ fullerite. Our study, based on Monte Carlo simulations, uncovers the existence of a new intermediate regime, between a low temperature ordered $(2 \times 2)$ state, and a high temperature $(1 \times 1)$ disordered phase. In the intermediate regime there is partial disorder, strongest for a subset of particularly frustrated surface molecules. These concepts and calculations provide a coherent understanding of experimental observations, with possible extension to other molecular crystal surfaces.
△ Less
Submitted 23 January, 2001;
originally announced January 2001.
-
Bent surface free energy differences from simulation
Authors:
U. Tartaglino,
D. Passerone,
E. Tosatti,
F. Di Tolla
Abstract:
We present a calculation of the change of free energy of a solid surface upon bending of the solid. It is based on extracting the surface stress through a molecular dynamics simulation of a bent slab by using a generalized stress theorem formula, and subsequent integration of the stress with respect to strain as a function of bending curvature. The method is exemplified by obtaining and comparin…
▽ More
We present a calculation of the change of free energy of a solid surface upon bending of the solid. It is based on extracting the surface stress through a molecular dynamics simulation of a bent slab by using a generalized stress theorem formula, and subsequent integration of the stress with respect to strain as a function of bending curvature. The method is exemplified by obtaining and comparing free energy changes with curvature of various reconstructed Au(001) surfaces.
△ Less
Submitted 28 November, 2000;
originally announced November 2000.
-
Surface molecular dynamics simulation with two orthogonal surface steps: how to beat the particle conservation problem
Authors:
D. Passerone,
U. Tartaglino,
F. Ercolessi,
E. Tosatti
Abstract:
Due to particle conservation, Canonical Molecular Dynamics (MD) simulations fail in the description of surface phase transitions involving coverage or lateral density changes. However, a step on the surface can act effectively as a source or a sink of atoms, in the simulation as well as in real life. A single surface step can be introduced by suitably modifying planar Periodic Boundary Condition…
▽ More
Due to particle conservation, Canonical Molecular Dynamics (MD) simulations fail in the description of surface phase transitions involving coverage or lateral density changes. However, a step on the surface can act effectively as a source or a sink of atoms, in the simulation as well as in real life. A single surface step can be introduced by suitably modifying planar Periodic Boundary Conditions (PBC), to accommodate the generally inequivalent stacking of two adjacent layers. We discuss here how, through the introduction of two orthogonal surface steps, particle number conservation may no longer represent a fatal constraint for the study of these surface transitions. As an example, we apply the method for estimating temperature-induced lateral density increase of the reconstructed
Au (001) surface; the resulting anisotropic cell change is consistent with experimental observations. Moreover, we implement this kind of scheme in conjunction with the variable curvature MD method, recently introduced by our group.
△ Less
Submitted 25 October, 2000;
originally announced October 2000.
-
(Meta-)stable reconstructions of the diamond(111) surface: interplay between diamond- and graphite-like bonding
Authors:
A. V. Petukhov,
D. Passerone,
F. Ercolessi,
E. Tosatti,
A. Fasolino
Abstract:
Off-lattice Grand Canonical Monte Carlo simulations of the clean diamond (111) surface, based on the effective many-body Brenner potential, yield the $(2\times1)$ Pandey reconstruction in agreement with \emph{ab-initio} calculations and predict the existence of new meta-stable states, very near in energy, with all surface atoms in three-fold graphite-like bonding. We believe that the long-standi…
▽ More
Off-lattice Grand Canonical Monte Carlo simulations of the clean diamond (111) surface, based on the effective many-body Brenner potential, yield the $(2\times1)$ Pandey reconstruction in agreement with \emph{ab-initio} calculations and predict the existence of new meta-stable states, very near in energy, with all surface atoms in three-fold graphite-like bonding. We believe that the long-standing debate on the structural and electronic properties of this surface could be solved by considering this type of carbon-specific configurations.
△ Less
Submitted 18 March, 2000; v1 submitted 27 January, 2000;
originally announced January 2000.
-
Realistic simulations of Au(100): Grand Canonical Monte Carlo and Molecular Dynamics
Authors:
D. Passerone,
F. Ercolessi,
F. Celestini,
E. Tosatti
Abstract:
The large surface density changes associated with the (100) noble metals surface hex-reconstruction suggest the use of non-particle conserving simulation methods. We present an example of a surface Grand Canonical Monte Carlo applied to the transformation of a square non reconstructed surface to the hexagonally covered low temperature stable Au(100). On the other hand, classical Molecular Dynami…
▽ More
The large surface density changes associated with the (100) noble metals surface hex-reconstruction suggest the use of non-particle conserving simulation methods. We present an example of a surface Grand Canonical Monte Carlo applied to the transformation of a square non reconstructed surface to the hexagonally covered low temperature stable Au(100). On the other hand, classical Molecular Dynamics allows to investigate microscopic details of the reconstruction dynamics, and we show, as an example, retraction of a step and its interplay with the surface reconstruction/deconstruction mechanism.
△ Less
Submitted 19 January, 2000;
originally announced January 2000.
-
Islands, craters, and a moving surface step on a hexagonally reconstructed (100) noble metal surface
Authors:
D. Passerone,
F. Ercolessi,
E. Tosatti
Abstract:
Deposition/removal of metal atoms on the hex reconstructed (100) surface of Au, Pt and Ir should present intriguing aspects, since a new island implies hex -> square deconstruction of the substrate, and a new crater the square -> hex reconstruction of the uncovered layer. To obtain a microscopic understanding of how islands/craters form in these conditions, we have conducted simulations of islan…
▽ More
Deposition/removal of metal atoms on the hex reconstructed (100) surface of Au, Pt and Ir should present intriguing aspects, since a new island implies hex -> square deconstruction of the substrate, and a new crater the square -> hex reconstruction of the uncovered layer. To obtain a microscopic understanding of how islands/craters form in these conditions, we have conducted simulations of island and crater growth on Au(100), whose atomistic behavior, including the hex reconstruction on top of the square substrate, is well described by mean s of classical many-body forces. By increasing/decreasing the Au coverage on Au(100), we find that island/craters will not grow unless they exceed a critical size of about 8-10 atoms. This value is close to that which explains the nonlinear coverage dependence observed in molecular adsorption on the closely related surface Pt (100). This threshold size is rationalized in terms of a transverse step correlation length, measuring the spatial extent where reconstruction of a given plane is disturbed by the nearby step.
△ Less
Submitted 14 January, 2000;
originally announced January 2000.
-
Reentrant layering in rare gas adsorption: preroughening or premelting?
Authors:
F. Celestini,
D. Passerone,
F. Ercolessi,
E. Tosatti
Abstract:
The reentrant layering transitions found in rare gas adsorption on solid substrates have conflictually been explained either in terms of preroughening (PR), or of top layer melting-solidification phenomena. We obtain adsorption isotherms of Lennard-Jones particles on an attractive substrate by off lattice Grand Canonical Monte Carlo (GCMC) simulation, and reproduce reentrant layering. Microscopi…
▽ More
The reentrant layering transitions found in rare gas adsorption on solid substrates have conflictually been explained either in terms of preroughening (PR), or of top layer melting-solidification phenomena. We obtain adsorption isotherms of Lennard-Jones particles on an attractive substrate by off lattice Grand Canonical Monte Carlo (GCMC) simulation, and reproduce reentrant layering. Microscopic analysis, including layer-by-layer occupancies, surface diffusion and pair correlations, confirms the switch of the top surface layer from solid to quasi-liquid across the transition temperature. At the same time, layer occupancy is found at each jump to switch from close to full to close to half, indicating a disordered flat (DOF) surface and establishing preroughening as the underlying mechanism. Our results suggest that top layer melting is essential in triggering preroughening, which thus represents the threshold transition to surface melting in rare gas solids.
△ Less
Submitted 10 January, 2000;
originally announced January 2000.
-
Variable Curvature Slab Molecular Dynamics as a Method to Determine Surface Stress
Authors:
Daniele Passerone,
Erio Tosatti,
Guido L. Chiarotti,
Furio Ercolessi
Abstract:
A thin plate or slab, prepared so that opposite faces have different surface stresses, will bend as a result of the stress difference. We have developed a classical molecular dynamics (MD) formulation where (similar in spirit to constant-pressure MD) the curvature of the slab enters as an additional dynamical degree of freedom. The equations of motion of the atoms have been modified according to…
▽ More
A thin plate or slab, prepared so that opposite faces have different surface stresses, will bend as a result of the stress difference. We have developed a classical molecular dynamics (MD) formulation where (similar in spirit to constant-pressure MD) the curvature of the slab enters as an additional dynamical degree of freedom. The equations of motion of the atoms have been modified according to a variable metric, and an additional equation of motion for the curvature is introduced. We demonstrate the method to Au surfaces, both clean and covered with Pb adsorbates, using many-body glue potentials. Applications to stepped surfaces, deconstruction and other surface phenomena are under study.
△ Less
Submitted 11 June, 1998;
originally announced June 1998.
-
Realistic grand canonical Monte Carlo surface simulation: application to Ar(111)
Authors:
Franck Celestini,
Daniele Passerone,
Furio Ercolessi,
Erio Tosatti
Abstract:
Most realistic, off-lattice surface simulations are done canonically--- conserving particles. For some applications, however, such as studying the thermal behavior of rare gas solid surfaces, these constitute bad working conditions. Surface layer occupancies are believed to change with temperature, particularly at preroughening, and naturally call for a grand canonical approach, where particle n…
▽ More
Most realistic, off-lattice surface simulations are done canonically--- conserving particles. For some applications, however, such as studying the thermal behavior of rare gas solid surfaces, these constitute bad working conditions. Surface layer occupancies are believed to change with temperature, particularly at preroughening, and naturally call for a grand canonical approach, where particle number is controlled by a chemical potential. We report preliminary results of novel realistic grand canonical Monte Carlo simulations of the Lennard-Jones (LJ) fcc(111) surface, believed to represent a quantitative model of e.g. Ar(111). The results are successful and highly informative for temperatures up to roughly 0.8 T_m, where clear precursor signals of preroughening are found. At higher temperatures, convergence to equilibrium is hampered by large fluctuations.
△ Less
Submitted 8 October, 1997;
originally announced October 1997.