Source author record

Erio Tosatti

Erio Tosatti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

66works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

66 published item(s)

preprint2022arXiv

Critical Peeling of Tethered Nanoribbons

The peeling of an immobile adsorbed membrane is a well known problem in engineering and macroscopic tribology. In the classic setup, picking up at one extreme and pulling off results in a peeling force that is a decreasing function of the pickup angle. As one end is lifted, the detachment front retracts to meet the immobile tail. At the nanoscale, interesting situations arise with the peeling of graphene nanoribbons (GNRs) on gold, as realized, e.g., by atomic force microscopy. The nanosized system shows a constant-force steady peeling regime, where the tip lifting h produces no retraction of the ribbon detachment point, and just an advancement h of the free tail end. This is opposite to the classic case, where the detachment point retracts and the tail end stands still. Here we characterise, by analytical modeling and numerical simulations, a third, experimentally relevant, setup where the nanoribbon, albeit structurally lubric, does not have a freely moving tail end, which is instead elastically tethered. Surprisingly, novel nontrivial scaling exponents appear that regulate the peeling evolution. As the detachment front retracts and the tethered tail is stretched, power laws of h characterize the shrinking of the adhered length the growth of peeling force and the peeling angle. These exponents precede the final total detachment as a critical point, where the entire ribbon eventually hangs suspended between the tip and tethering spring. These analytical predictions are confirmed by realistic MD simulations, retaining the full atomistic description, also confirming their survival at finite experimental temperatures.

preprint2022arXiv

Multiwalled nanotube faceting unravelled

Nanotubes show great promise for miniaturizing advanced technologies. Their exceptional physical properties are intimately related to their morphological and crystal structure. Circumferential faceting of multiwalled nanotubes reinforces their mechanical strength and alters their tribological and electronic properties. Here, the nature of this important phenomenon is fully rationalized in terms of interlayer registry patterns. Regardless of the nanotube identity (that is, diameter, chirality, chemical composition), faceting requires the matching of the chiral angles of adjacent layers. Above a critical diameter that corresponds well with experimental results, achiral multiwalled nanotubes display evenly spaced extended axial facets whose number equals the interlayer difference in circumferential unit cells. Elongated helical facets, commonly observed in experiment, appear in nanotubes that exhibit small interlayer chiral angle mismatch. When the wall chiralities are uncorrelated, faceting is suppressed and outer layer corrugation, which is induced by the moiré superlattice, is obtained in agreement with experiments. Finally, we offer an explanation for the higher incidence of faceting in multiwalled boron nitride nanotubes with respect to their carbon-based counterparts.

preprint2022arXiv

The Smallest Archimedean Screw: Facet Dynamics and Friction in Multi-Walled Nanotubes

We identify a new material phenomenon, where minute mechanical manipulations induce pronounced global structural reconfigurations in faceted multi-walled nanotubes. This behavior has strong implications on the tribological properties of these systems and may be the key to understand the enhanced inter-wall friction recently measured for boron-nitride nanotubes with respect to their carbon counterparts. Notably, the fast rotation of helical facets in these systems upon coaxial sliding may serve as a nanoscale Archimedean screw for directional transport of physisorbed molecules.

preprint2021arXiv

Amplitude Nanofriction Spectroscopy

Atomic scale friction, an indispensable element of nanotechnology, requires a direct access to, under actual growing shear stress, its successive live phases: from static pinning, to depinning and transient evolution, eventually ushering in steady state kinetic friction. Standard tip-based atomic force microscopy generally addresses the steady state, but the prior intermediate steps are much less explored. Here we present an experimental and simulation approach, where an oscillatory shear force of increasing amplitude leads to a one-shot investigation of all these successive aspects. Demonstration with controlled gold nanocontacts sliding on graphite uncovers phenomena that bridge the gap between initial depinning and large speed sliding, of potential relevance for atomic scale time and magnitude dependent rheology.

preprint2020arXiv

Proton strings and rings in atypical nucleation of ferroelectricity in ice

Ordinary ice has a proton-disordered phase which is kinetically metastable, unable to reach, spontaneously, the ferroelectric (FE) ground state at low temperature where a residual Pauling entropy persists. Upon light doping with KOH at low temperature, the transition to FE ice takes place, but its microscopic mechanism still needs clarification. We introduce a lattice model based on dipolar interactions plus a competing, frustrating term that enforces the ice rule (IR). In the absence of IR-breaking defects, standard Monte Carlo (MC) simulation leaves this ice model stuck in a state of disordered proton ring configurations with the correct Pauling entropy. A replica exchange accelerated MC sampling strategy succeeds, without open path moves, interfaces, or off-lattice configurations, in equilibrating this defect-free ice, reaching its low-temperature FE order through a well-defined first-order phase transition. When proton vacancies mimicking the KOH impurities are planted into the IR-conserving lattice, they enable standard MC simulation to work, revealing the kinetics of evolution of ice from proton disorder to partial FE order below the transition temperature. Replacing ordinary nucleation, each impurity opens up a proton ring generating a linear string, an actual FE hydrogen bond wire that expands with time. Reminiscent of those described for spin ice, these impurity-induced strings are proposed to exist in doped water ice too, where IRs are even stronger. The emerging mechanism yields a dependence of the long-time FE order fraction upon dopant concentration, and upon quenching temperature, that compares favorably with that known in real-life KOH doped ice.

preprint2016arXiv

CDW slips and giant frictional dissipation peaks at the NbSe$_2$ surface

Accessing, controlling and understanding nanoscale friction and dissipation is a crucial issue in nanotechnology, where moving elements are central. Recently, ultra-sensitive noncontact pendulum Atomic Force Microscope (AFM) succeeded in detecting the electronic friction drop caused by the onset of superconductivity in Nb, raising hopes that a wider variety of mechanisms of mechanical dissipation arising from electron organization into different collective phenomena will become accessible through this unconventional surface probe. Among them, the driven phase dynamics of charge-density-waves (CDWs) represents an outstanding challenge as a source of dissipation. Here we report a striking multiplet of AFM dissipation peaks arising at nanometer distances above the surface of NbSe$_2$ - a layered compound exhibiting an incommensurate CDW. Each peak appears at a well defined tip-surface interaction force of the order of a nN, and persists until T=70K where CDW short-range order is known to disappear. A theoretical model is presented showing that the peaks are connected to tip-induced local 2$π$ CDW phase slips. Under the attractive potential of the approaching tip, the local CDW surface phase landscape deforms continuously until a series of 2$π$ jumps occur between different values of the local phase. As the tip oscillates to and fro, each slip gives rise to a hysteresis cycle, appearing at a selected distance, the dissipation corresponding to "pumping" in and out a local slip in the surface CDW phase of NbSe$_2$.

preprint2016arXiv

Creating new layered structures at high pressures: SiS$_2$

Old and novel layered structures are attracting increasing attention for their physical, electronic, and frictional properties. SiS$_2$, isoelectronic to SiO$_2$, CO$_2$ and CS$_2$, is a material whose phases known experimentally up to 6 GPa exhibit 1D chain-like, 2D layered and 3D tetrahedral structures. We present highly predictive $ab$ $initio$ calculations combined with evolutionary structure search and molecular dynamics simulations of the structural and electronic evolution of SiS$_2$ up to 100 GPa. A highly stable CdI$_2$-type layered structure, which is octahedrally coordinated with space group $P\bar{3}m1$ surprisingly appears between 4 and up to at least 100 GPa. The tetrahedral-octahedral switch is naturally expected upon compression, unlike the layered character realized here by edge-sharing SiS$_6$ octahedral units connecting within but not among sheets. The predicted phase is semiconducting with an indirect band gap of about 2 eV at 10 GPa, decreasing under pressure until metallization around 40 GPa. The robustness of the layered phase suggests possible recovery at ambient pressure, where calculated phonon spectra indicate dynamical stability. Even a single monolayer is found to be dynamically stable in isolation, suggesting that it could possibly be sheared or exfoliated from bulk $P\bar{3}m1$-SiS$_2$.

preprint2016arXiv

Metallic, Magnetic and Molecular Nanocontacts

Scanning tunnelling microscopy and break-junction experiments realize metallic and molecular nanocontacts that act as ideal one-dimensional channels between macroscopic electrodes. Emergent nanoscale phenomena typical of these systems encompass structural, mechanical, electronic, transport, and magnetic properties. This Review focuses on the theoretical explanation of some of these properties obtained with the help of first-principles methods. By tracing parallel theoretical and experimental developments from the discovery of nanowire formation and conductance quantization in gold nanowires to recent observations of emergent magnetism and Kondo correlations, we exemplify the main concepts and ingredients needed to bring together ab initio calculations and physical observations. It can be anticipated that diode, sensor, spin-valve and spin-filter functionalities relevant for spintronics and molecular electronics applications will benefit from the physical understanding thus obtained.

preprint2016arXiv

Nanoscale orbital excitations and the infrared spectrum of a molecular Mott insulator: A15-Cs$_{3}$C$_{60}$

The quantum physics of ions and electrons behind low-energy spectra of strongly correlated {\it molecular} conductors, superconductors and Mott insulators is poorly known, yet fascinating especially in orbitally degenerate cases. The fulleride insulator Cs$_{3}$C$_{60}$ (A15), one such system, exhibits infrared (IR) spectra with low temperature peak features and splittings suggestive of static Jahn-Teller distortions with breakdown of orbital symmetry in the molecular site. That is puzzling, for there is no detectable static distortion, and because the features and splittings disappear upon modest heating, which they should not. Taking advantage of the Mott-induced collapse of electronic wavefunctions from lattice-extended to nanoscale localized inside a caged molecular site, we show that unbroken spin and orbital symmetry of the ion multiplets explains the IR spectrum without adjustable parameters. This demonstrates the importance of a fully quantum treatment of nuclear positions and orbital momenta in the Mott insulator sites, dynamically but not statically distorted. The observed demise of these features with temperature is explained by the thermal population of a multiplet term whose nuclear positions are essentially undistorted, but whose energy is very low-lying. That term is in facts a scaled-down orbital excitation analogous to that of other Mott insulators, with the same spin $1/2$ as the ground state, but with a larger orbital moment of two instead of one.

preprint2016arXiv

Non-equilibrium and non-homogeneous phenomena around a first-order quantum phase transition

We consider non-equilibrium phenomena in a very simple model that displays a zero-temperature first-order phase transition. The quantum Ising model with a four-spin exchange is adopted as a general representative of first-order quantum phase transitions that belong to the Ising universality class, such as for instance the order-disorder ferroelectric transitions, and possibly first-order T = 0 Mott transitions. In particular, we address quantum quenches in the exactly solvable limit of infinite connectivity and show that, within the coexistence region around the transition, the system can remain trapped in a metastable phase, as long as it is spatially homogeneous so that nucleation can be ignored. Motivated by the physics of nucleation, we then study in the same model static but inhomogeneous phenomena that take place at surfaces and interfaces. The first order nature implies that both phases remain locally stable across the transition, and with that the possibility of a metastable wetting layer showing up at the surface of the stable phase, even at T=0. We use mean-field theory plus quantum fluctuations in the harmonic approximation to study quantum surface wetting.

preprint2016arXiv

Pressure-induced gap closing and metallization of MoSe$_{2}$ and MoTe$_{2}$

Layered molybdenum dichalchogenides are semiconductors whose gap is controlled by delicate interlayer interactions. The gap tends to drop together with the interlayer distance, suggesting collapse and metallization under pressure. We predict, based on first principles calculations, that layered semiconductors 2H$_c$-MoSe$_2$ and 2H$_c$-MoTe$_2$ should undergo metallization at pressures between 28 and 40 GPa (MoSe$_2$) and 13 and 19 GPa (MoTe$_2$). Unlike MoS$_2$ where a 2H$_c$ $\to$ 2H$_a$ layer sliding transition is known to take place, these two materials appear to preserve the original 2H$_c$ layered structure at least up to 100 GPa and to increasingly resist lubric layer sliding under pressure. Similar to metallized MoS$_2$ they are predicted to exhibit a low density of states at the Fermi level, and presumably very modest superconducting temperatures if any. We also study the $β$-MoTe$_2$ structure, metastable with a higher enthalpy than 2H$_c$-MoTe$_2$. Despite its ready semimetallic and (weakly) superconducting character already at zero pressure, metallicity is not expected to increase dramatically with pressure.

preprint2016arXiv

Slider Thickness Promotes Lubricity: from 2D Islands to 3D Clusters

The sliding of three-dimensional clusters and two-dimensional islands adsorbed on crystal surfaces represent an important test case to understand friction. Even for the same material, monoatomic islands and thick clusters will not as a rule exhibit the same friction, but specific differences have not been explored. Through realistic molecular dynamics simulations of the static friction gold on graphite, an experimentally relevant system, we uncover as a function of gold thickness a progressive drop of static friction from monolayer islands, that are easily pinned, towards clusters, that slide more readily. The main ingredient contributing to this thickness-induced lubricity appears to be the increased effective rigidity of the atomic contact, acting to reduce the cluster interdigitation with the substrate. A second element which plays a role is lateral contact size, which can accommodate the solitons typical of the incommensurate interface only above a critical contact diameter, which is larger for monolayer islands than for thick clusters. The two effects concur to make clusters more lubric than islands, and large sizes more lubric than smaller ones. These conclusions are expected to be of broader applicability in diverse nanotribological systems, where the role played by static, and dynamic, friction is generally quite important.

preprint2015arXiv

Electrical charging effects on the sliding friction of a model nano-confined ionic liquid

Recent measurements suggest the possibility to exploit ionic liquids (ILs) as smart lubricants for nano-contacts, tuning their tribological and rheological properties by charging the sliding interfaces. Following our earlier theoretical study of charging effects on nanoscale confinement and squeezout of a model IL, we present here molecular dynamics simulations of the frictional and lubrication properties of that model under charging conditions.First we describe the case when two equally charged plates slide while being held together to a confinement distance of a few molecular layers.The shear sliding stress is found to rise as the number of IL layers decreases stepwise. However the shear stress shows, within each given number of layers, only a weak dependence upon the precise value of the normal load, a result in agreement with data extracted from recent experiments.We subsequently describe the case of opposite charging of the sliding plates, and follow the shear stress when the charging is slowly and adiabatically reversed in the course of time, under fixed load. Despite the fixed load, the number and structure of the confined IL layers changes with changing charge, and that in turn drives strong friction variations. The latter involve first of all charging-induced freezing of the IL film, followed by a discharging-induced melting, both made possible by the nanoscale confinement. Another mechanism for charging-induced frictional changes is a shift of the plane of maximum shear from mid-film to the plate-film interface, and viceversa. While these occurrences and results invariably depend upon the parameters of the model IL and upon its specific interaction with the plates, the present study helps identifying a variety of possible behavior, obtained under very simple assumptions, while connecting it to an underlying equilibrium thermodynamics picture.

preprint2015arXiv

Friction Boosted by Equilibrium Misalignment of Incommensurate Two-Dimensional Colloid Monolayers

Colloidal 2D monolayers sliding in an optical lattice are of recent importance as a frictional system. In the general case when the monolayer and optical lattices are incommensurate, we predict two important novelties, one in the static equilibrium structure, the other in the frictional behavior under sliding. Structurally, realistic simulations show that the colloid layer should possess in full equilibrium a small misalignment rotation angle relative to the optical lattice, an effect so far unnoticed but visible in some published experimental moiré patterns. Under sliding, this misalignment has the effect of boosting the colloid monolayer friction by a considerable factor over the hypothetical aligned case discussed so far. A frictional increase of similar origin must generally affect other incommensurate adsorbed monolayers and contacts, to be sought out case by case.

preprint2015arXiv

Frictional Transition from Superlubric Islands to Pinned Monolayers

The inertial sliding of physisorbed submonolayer islands on crystal surfaces contains unexpected information on the exceptionally incommensurate smooth sliding state associated with superlubricity and on the mechanisms of its disappearance. Here we show in a joint quartz crystal microbalance (QCM) and molecular dynamics (MD) simulation case study of Xe on Cu(111) how superlubricity emerges in the large size limit of naturally incommensurate Xe islands. Theory also predicts, as coverage approaches a full monolayer, an abrupt adhesion-driven 2D density compression of order several %, implying here a hysteretic jump from superlubric free islands to a pressurized commensurate immobile monolayer. This scenario is fully supported by QCM data showing remarkably large slip times at increasing submonolayer coverages, signaling superlubricity, followed by a dramatic drop to zero for the dense commensurate monolayer. Careful analysis of this variety of island sliding phenomena will be mandatory in future applications of friction of crystal-adsorbate interfaces.

preprint2015arXiv

High pressure layered structure of carbon disulfide

Solid CS$_{2}$ is superficially similar to CO$_{2}$, with the same $Cmca$ molecular crystal structure at low pressures, which has suggested similar phases also at high pressures. We carried out an extensive first principles evolutionary search in order to identify the zero temperature lowest enthalpy structures of CS$_{2}$ for increasing pressure up to 200\,GPa. Surprisingly, the molecular $Cmca$ phase does not evolve into $β$-cristobalite as in CO$_{2}$, but transforms instead into phases HP2 and HP1, both recently described in high pressure SiS$_{2}$. HP1 in particular, with a wide stability range, is a layered $P2_{1}/c$ structure characterized by pairs of edge-sharing tetrahedra, and theoretically more robust than all other CS$_{2}$ phases discussed so far. Its predicted Raman spectrum and pair correlation function agree with experiment better than those of $β$-cristobalite, and further differences are predicted between their respective IR spectra. The band gap of HP1-CS$_{2}$ is calculated to close under pressure yielding an insulator-metal transition near 50 GPa in agreement with experimental observations. However, the metallic density of states remains modest above this pressure, suggesting a different origin for the reported superconductivity.

preprint2015arXiv

Subharmonic Shapiro steps of sliding colloidal monolayers in optical lattices

We investigate theoretically the possibility to observe dynamical mode locking, in the form of Shapiro steps, when a time-periodic potential or force modulation is applied to a two-dimensional (2D) lattice of colloidal particles that are dragged by an external force over an optically generated periodic potential. Here we present realistic molecular dynamics simulations of a 2D experimental setup, where the colloid sliding is realized through the motion of soliton lines between locally commensurate patches or domains, and where the Shapiro steps are predicted and analyzed. Interestingly, the jump between one step and the next is seen to correspond to a fixed number of colloids jumping from one patch to the next, across the soliton line boundary, during each AC cycle. In addition to ordinary "integer" steps, coinciding here with the synchronous rigid advancement of the whole colloid monolayer, our main prediction is the existence of additional smaller "subharmonic" steps due to localized solitonic regions of incommensurate layers executing synchronized slips, while the majority of the colloids remains pinned to a potential minimum. The current availability and wide parameter tunability of colloid monolayers makes these predictions potentially easy to access in an experimentally rich 2D geometrical configuration.

preprint2015arXiv

Superlubric-Pinned Transition in Sliding Incommensurate Colloidal Monolayers

Two-dimensional (2D) crystalline colloidal monolayers sliding over a laser-induced optical lattice recently emerged as a new tool for the study of friction between ideal crystal surfaces. Here we focus in particular on static friction, the minimal sliding force necessary to depin one lattice from the other. If the colloid and the optical lattices are mutually commensurate, the colloid sliding is always pinned by static friction; but when they are incommensurate the presence or absence of pinning can be expected to depend upon the system parameters. If a 2D analogy to the mathematically established Aubry transition of one-dimensional systems were to hold, an increasing periodic corrugation strength $U_0$ should turn an initially free-sliding monolayer into a pinned state through a well-defined dynamical phase transition. We address this problem by the simulated sliding of a realistic model 2D colloidal lattice, confirming the existence of a clear and sharp superlubric-pinned transition for increasing corrugation strength. Unlike the 1D Aubry transition which is continuous, the 2D transition exhibits a definite first-order character. With no change of symmetry, the transition entails a structural character, with a sudden increase of the colloid-colloid interaction energy, accompanied by a compensating downward jump of the colloid-corrugation energy. The transition value for the corrugation amplitude $U_0$ depends upon the misalignment angle $θ$ between the optical and the colloidal lattices, superlubricity surviving until larger corrugations for angles away from the energetically favored orientation, which is itself generally slightly misaligned, as shown in recent work. The observability of the superlubric-pinned colloid transition is proposed and discussed.

preprint2014arXiv

Charge density wave surface phase slips and non-contact nanofriction

Bulk electrical dissipation caused by charge-density-wave (CDW) depinning and sliding is a classic subject. We present a novel local, nanoscale mechanism describing the occurrence of mechanical dissipation peaks in the dynamics of an atomic force microscope tip oscillating above the surface of a CDW material. Local surface 2$π$ slips of the CDW phase are predicted to take place giving rise to mechanical hysteresis and large dissipation at discrete tip surface distances. The results of our static and dynamic numerical simulations are believed to be relevant to recent experiments on NbSe$_2$; other candidate systems in which similar effects should be observable are also discussed.

preprint2014arXiv

Crystal structure search and electronic properties of alkali doped phenanthrene and picene

Alkali doped aromatic compounds have shown evidence of metallic and superconducting phases whose precise nature is still mysterious. In potassium and rubidium doped phenanthrene, superconducting temperatures around 5 K have been detected, but such basic elements as the stoichiometry, crystal structure, and electronic bands are still speculative. We seek to predict the crystal structure of M3-phenanthrene (M = K, Rb) using ab-initio evolutionary simulation in conjunction with density functional theory (DFT), and find metal but also insulator phases with distinct structures. The original P21 herringbone structure of the pristine molecular crystal is generally abandoned in favor of different packing and chemical motifs. The metallic phases are frankly ionic with three electrons acquired by each molecule. In the nonmagnetic insulating phases the alkalis coalesce reducing the donated charge from three to two per phenanthrene molecule. A similar search for K3-picene yields an old and a new structure, with unlike potassium positions and different electronic bands, but both metallic retaining the face-to-edge herringbone structure and the P21 symmetry of pristine picene. Both the new K3-picene and the best metallic M3-phenanthrene are further found to undergo a spontaneous transition from metal to antiferromagnetic insulator when spin polarization is allowed, a transition which is not necessarily real, but which underlines the necessity to include correlations beyond DFT. Features of the metallic phases that may be relevant to phonon-driven superconductivity are underlined.

preprint2014arXiv

Does Rotational Melting Make Molecular Crystal Surfaces More Slippery?

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.

preprint2014arXiv

High Pressure Structure and Decomposition of MoS$_{2}$

The high pressure structural and electronic evolution of bulk MoS$_2$, an important transition metal layered dichalchogenide, is currently under active investigation. Recent theoretical and experimental work predicted and verified a 2H$_c \to$ 2H$_a$ layer sliding structural transition at 20 GPa and a band overlap semiconductor-semimetal transition in the same pressure range. The 2H$_a$ structure is known to persist up to pressure of 81 GPa but properties at higher pressures remain experimentally unknown. Here we predict, with a reliable first-principles evolutionary search, that major structural transformations should take place in equilibrium at higher pressures near 130-140 GPa. The main motif is a decomposition into MoS + S, also heralded in a small bimolecular cell by the appearance of a metastable non-layered metallic MoS$_2$ structure with space group \textit{P4/mmm}. Unlike semimetallic 2H$_a$-MoS$_2$, both this phase and sulphur in the fully phase separated system are fully metallic and superconducting with higher critical temperatures than alkali-intercalated MoS$_2$.

preprint2014arXiv

Non-contact dissipation on the surface of charge density wave materials

Bulk electrical dissipation caused by charge-density-wave (CDW) depinning and sliding is a classic subject. We present a novel local, nanoscale mechanism describing the occurrence of mechanical dissipation peaks in the dynamics of an atomic force microscope tip oscillating above the surface of a CDW material. Local surface 2$π$ slips of the CDW phase are predicted to take place giving rise to mechanical hysteresis and large dissipation at discrete tip surface distances. The results of our static and dynamic numerical simulations are believed to be relevant to recent experiments on NbSe$_2$; other candidate systems in which similar effects should be observable are also discussed.

preprint2014arXiv

Robust $s\pm$ Superconductivity in a Two-Band Hubbard-Fr{ö}hlich Model of Alkali Doped Organics

The damaging effect of strong electron-electron repulsion on regular, electron-phonon %$s$-wave superconductivity is a standard tenet. In spite of that, an increasing number of compounds such as fullerides and more recently alkali-doped aromatics exhibit %$s$-wave or presumably $s$ wave superconductivity despite very narrow bands and very strong electron repulsion. Here, we explore superconducting solutions of a model Hamiltonian inspired by the electronic structure of alkali doped aromatics. The model is a two-site, two-narrow-band metal with a single intersite phonon, leading to attraction-mediated, two-order parameter superconductivity. On top of that, the model includes a repulsive on-site Hubbard $U$, whose effect on the superconductivity we study. Starting within mean field, we find that $s \pm$ superconductivity is the best solution surviving the presence of $U$, whose effect is canceled out by the opposite signs of the two order parameters. The correlated Gutzwiller study that follows is necessary because without electron correlations the superconducting state would in this model be superseded by an antiferromagnetic insulating state with lower energy. The Gutzwiller correlations lower the energy of the metallic state, with the consequence that the $s \pm$ superconducting state is stabilized and even strengthened for small Hubbard $U$.

preprint2014arXiv

Shape and area fluctuation effects on nucleation theory

In standard nucleation theory, the nucleation process is characterized by computing $ΔΩ(V)$, the reversible work required to form a cluster of volume $V$ of the stable phase inside the metastable mother phase. However, other quantities besides the volume could play a role in the free energy of cluster formation, and this will in turn affect the nucleation barrier and the shape of the nucleus. Here we exploit our recently introduced mesoscopic theory of nucleation to compute the free energy cost of a nearly-spherical cluster of volume $V$ and a fluctuating surface area $A$, whereby the maximum of $ΔΩ(V)$ is replaced by a saddle point in $ΔΩ(V,A)$. Compared to the simpler theory based on volume only, the barrier height of $ΔΩ(V,A)$ at the transition state is systematically larger by a few $k_BT$. More importantly, we show that, depending on the physical situation, the most probable shape of the nucleus may be highly non spherical, even when the surface tension and stiffness of the model are isotropic. Interestingly, these shape fluctuations do not influence or modify the standard Classical Nucleation Theory manner of extracting the interface tension from the logarithm of the nucleation rate near coexistence.

preprint2014arXiv

Shot noise and magnetism of Pt atomic chains: accumulation of points at the boundary

Pt is known to show spontaneous formation of monatomic chains upon breaking a metallic contact. From model calculations, these chains are expected to be spin polarized. However, direct experimental evidence for or against magnetism is lacking. Here, we investigate shot noise as a potential source of information on the magnetic state of Pt atomic chains. We observe a remarkable structure in the distribution of measured shot-noise levels, where the data appear to be confined to the region of nonmagnetic states. While this suggests a nonmagnetic ground state for the Pt atomic chains, from calculations we find that the magnetism in Pt chains is due to 'actor' electron channels, which contribute very little to ballistic conductance and noise. On the other hand, there are weakly polarized 'spectator' channels, which carry most of the current and are only slightly modified by the magnetic state.

preprint2014arXiv

Spin-orbit interaction and Dirac cones in d orbital noble metal surface states

Band splittings, chiral spin polarization and topological surface states generated by spin-orbit interactions at crystal surfaces are receiving a lot of attention for their potential device applications as well as fascinating physical properties. Most studies have focused on sp states near the Fermi energy, which are relevant for transport and have long lifetimes. Far less explored, though in principle stronger, are spin-orbit interaction effecs within d states, including those deep below the Fermi energy. Here, we report a joint photoemission/ab initio study of spin-orbit effects in the deep d orbital surface states of a 24-layer Au film grown on Ag(111) and a 24-layer Ag film grown on Au(111), singling out a conical intersection (Dirac cone) between two surface states in a large surface-projected gap at the time-reversal symmetric M points. Unlike the often isotropic dispersion at Gamma point Dirac cones, the M point cones are strongly anisotropic. An effective k.p Hamiltonian is derived to describe the anisotropic band splitting and spin polarization near the Dirac cone.

preprint2013arXiv

A fingerprint of surface-tension anisotropy in the free-energy cost of nucleation

We focus on the Gibbs free energy $ΔG$ for nucleating a droplet of the stable phase (e.g. solid) inside the metastable parent phase (e.g. liquid), close to the first-order transition temperature. This quantity is central to the theory of homogeneous nucleation, since it superintends the nucleation rate. We recently introduced a field theory describing the dependence of $ΔG$ on the droplet volume $V$, taking into account besides the microscopic fuzziness of the droplet-parent interface, also small fluctuations around the spherical shape whose effect, assuming isotropy, was found to be a characteristic logarithmic term. Here we extend this theory, introducing the effect of anisotropy in the surface tension, and show that in the limit of strong anisotropy $ΔG(V)$ once more develops a term logarithmic on $V$, now with a prefactor of opposite sign with respect to the isotropic case. Based on this result, we argue that the geometrical shape that large solid nuclei mostly prefer could be inferred from the prefactor of the logarithmic term in the droplet free energy, as determined from the optimization of its near-coexistence profile.

preprint2013arXiv

Electron-doped organics: Charge-disproportionate insulators and Hubbard-Fröhlich metals

Several examples of metallic electron doped polycyclic aromatic hydrocarbons (PAHs) molecular crystals have recently been experimentally proposed. Some of them have superconducting components, but most other details are still unknown beginning with structure and the nature of metallicity. We carried out ab-initio density functional calculations for La-Phenanthrene (La-PA), here meant to represent a generic case of three-electron doping, to investigate structure and properties of a conceptually simple case. To our surprise we found first of all that the lowest energy state is not metallic but band insulating, with a disproportionation of two inequivalent PA molecular ions and a low P1 symmetry, questioning the common assumption that three electrons will automatically metallize a PAH crystal. Our best metallic structure is metastable and slightly higher in energy, and retains equivalent PA ions and a higher P21 symmetry -- the same generally claimed for metallic PAHs. We show that a "dimerizing" periodic distortion opens very effectively a gap in place of a symmetry related degeneracy of all P21 structures near the Fermi level, foreshadowing a possible role of that special intermolecular phonon in superconductivity of metallic PAHs. A Hubbard-Frohlich model describing that situation is formulated for future studies.

preprint2013arXiv

Ferromagnetic Kondo effect in a triple quantum dot system

We propose that a simple device of three laterally-coupled quantum dots, the central one contacted by metal leads, can realize the ferromagnetic Kondo model, which is characterized by interesting properties like a non-analytic inverted zero-bias anomaly and an extreme sensitivity to a magnetic field. Furthermore, by tuning the gate voltages of the lateral dots, this device may allow to study the transition from ferromagnetic to antiferromagnetic Kondo effect, a simple case of a Berezinskii-Kosterlitz-Thouless transition. We model the device by three coupled Anderson impurities that we study by numerical renormalization group. We calculate the single-particle spectral function of the central dot, which at zero frequency is proportional to the zero-bias conductance, across the transition, both in the absence and in the presence of a magnetic field.

preprint2013arXiv

Kondo conductance across the smallest spin 1/2 radical molecule

Molecular contacts are generally poorly conducting because their energy levels tend to lie far from the Fermi energy of the metal contact, necessitating undesirably large gate and bias voltages in molecular electronics applications. Molecular radicals are an exception because their partly filled orbitals undergo Kondo screening, opening the way to electron passage even at zero bias. While that phenomenon has been experimentally demonstrated for several complex organic radicals, quantitative theoretical predictions have not been attempted so far. It is therefore an open question whether and to what extent an ab initio-based theory is able to make accurate predictions for Kondo temperatures and conductance lineshapes. Choosing nitric oxide NO as a simple and exemplary spin 1/2 molecular radical, we present calculations based on a combination of density functional theory and numerical renormalization group (DFT+NRG) predicting a zero bias spectral anomaly with a Kondo temperature of 15 K for NO/Au(111). A scanning tunneling spectroscopy study is subsequently carried out to verify the prediction, and a striking zero bias Kondo anomaly is confirmed, still quite visible at liquid nitrogen temperatures. Comparison shows that the experimental Kondo temperature of about 43 K is larger than the theoretical one, while the inverted Fano lineshape implies a strong source of interference not included in the model. These discrepancies are not a surprise, providing in fact an instructive measure of the approximations used in the modeling, which supports and qualifies the viability of the DFT+NRG approach to the prediction of conductance anomalies in larger molecular radicals.

preprint2013arXiv

Self-consistent Gutzwiller study of bcc Fe: interplay of ferromagnetic order and kinetic energy

The Gutzwiller technique has long been known as a method to include correlations in electronic structure calculations. Here we implement an ab-initio Gutzwiller+LDA calculation, and apply it to a classic problem, the ferromagnetism of bulk bcc Fe, whose nature has attracted recent interest. In the conventional Stoner-Wohlfarth model, the ferromagnetic ordering of iron sets in so that the electrons can reduce their mutual Coulomb repulsion at the extra cost of some increase of electron kinetic energy. Density functional theory within the spin polarized local density approximation (LDA) has long supported that picture, showing that ferromagnetic alignment causes band narrowing and a corresponding wavefunction localization, whence a kinetic energy increase. However, because of its inadequate treatment of strong intra-site correlations for localized d orbitals, LDA cannot be relied upon, particularly when it comes to separately describing fine potential and kinetic energy imbalances. With ab-initio Gutzwiller+LDA, we indeed find that the effect of correlations is to dramatically reverse the balance, the ferromagnetic ordering of Fe in fact causing a decrease of kinetic energy, at the cost of some increase of potential energy. The underlying physical mechanism, foreshadowed long ago by Goodenough and others, and more recently supported by LDA+DMFT calculations, is that correlations cause eg and t2g 3d orbitals to behave very differently. Weakly dispersive eg states are spin-polarized and almost localized, while, more than half filled, the t2g are broad band, fully delocalized states. Owing to intra-atomic Hund's rule exchange which aligns eg and t2g spins, the propagation of itinerant t2g holes is only allowed when different atomic spins are ferromagnetically aligned. We thus conclude that double exchange is at work already in the most popular ferromagnetic metal.

preprint2013arXiv

Size scaling of static friction

Sliding friction across a thin soft lubricant film typically occurs by stick-slip, the lubricant fully solidifying at stick, yielding and flowing at slip. The static friction force per unit area preceding slip is known from molecular dynamics (MD) simulations to decrease with increasing contact area. That makes the large-size fate of stick-slip unclear and unknown; its possible vanishing is important as it would herald smooth sliding with a dramatic drop of kinetic friction at large size. Here we formulate a scaling law of the static friction force, which for a soft lubricant is predicted to decrease as f_m + Δf /A^gamma for increasing contact area A, with gamma>0. Our main finding is that the value of f_m, controlling the survival of stick-slip at large size, can be evaluated by simulations of comparably small size. MD simulations of soft lubricant sliding are presented, which verify this theory.

preprint2013arXiv

Sliding on a Nanotube: Interplay of Friction, Deformations and Structure

The frictional properties of individual carbon nanotubes (CNTs) are studied by sliding an atomic force microscopy tip across and along its principle axis. This direction-dependent frictional behavior is found to correlate strongly with the presence of structural defects, surface chemistry, and CNT chirality. This study shows that it is experimentally possible to tune the frictional/adhesion properties of a CNT by controlling the CNT structure and surface chemistry, as well as use friction force to predict its structural and chemical properties.

preprint2013arXiv

Stick-Slip Nanofriction in Trapped Cold Ion Chains

Stick-slip -- the sequence of mechanical instabilities through which a slider advances on a solid substrate -- is pervasive throughout sliding friction, from nano to geological scales. Here we suggest that trapped cold ions in an optical lattice can also be of help in understanding stick-slip friction, and also the way friction changes when one of the sliders undergoes structural transitions. For that scope, we simulated the dynamical properties of a 101-ions chain, driven to slide back and forth by a slowly oscillating electric field in an incommensurate periodic "corrugation" potential of increasing magnitude U0. We found the chain sliding to switch, as U0 increases and before the Aubry transition, from a smooth-sliding regime with low dissipation to a stick-slip regime with high dissipation. In the stick-slip regime the onset of overall sliding is preceded by precursor events consisting of partial slips of few ions only, leading to partial depinning of the chain, a nutshell remnant of precursor events at the onset of motion also observed in macroscopic sliders. Seeking to identify the possible effects on friction of a structural transition, we reduced the trapping potential aspect ratio until the ion chain shape turned from linear to zigzag. Dynamic friction was found to rise at the transition, reflecting the opening of newer dissipation channels.

preprint2013arXiv

Structure change, layer sliding, and metallization in high-pressure MoS$_{2}$

Based on ab initio calculations and metadynamics simulations, we predict that 2H-MoS$_2$, a layered insulator, will metallize under pressures in excess of 20-30 GPa. In the same pressure range, simulations and enthalpy optimization predict a structural transition. Free mutual sliding of layers takes place at this transition, the original 2H$_c$ stacking changing to a 2H$_a$ stacking typical of 2H-NbSe$_2$, an event explaining for the first time previously mysterious X-ray diffraction and Raman spectroscopy data. Phonon and electron phonon calculations suggest that pristine MoS$_2$, once metallized, will require ultrahigh pressures in order to develop superconductivity.

preprint2012arXiv

Collective effects at frictional interfaces

We discuss the role of the long-range elastic interaction between the contacts inside an inhomogeneous frictional interface. The interaction produces a characteristic elastic correlation length $λ_c = a^2 E / k_c$ (where $a$ is the distance between the contacts, $k_c$ is the elastic constant of a contact, and $E$ is the Young modulus of the sliding body), below which the slider may be considered as a rigid body. The strong inter-contact interaction leads to a narrowing of the effective threshold distribution for contact breaking and enhances the chances for an elastic instability to appear. Above the correlation length, $r > λ_c$, the interaction leads to screening of local perturbations in the interface, or to appearance of collective modes --- frictional cracks propagating as solitary waves.

preprint2012arXiv

Hindered rolling and friction anisotropy in supported carbon nanotubes

Carbon nanotubes (CNTs) are well known for their exceptional thermal, mechanical and electrical properties. For many CNT applications it is of the foremost importance to know their frictional properties. However, very little is known about the frictional forces between an individual nanotube and a substrate or tip. Here, we present a combined theoretical and experimental study of the frictional forces encountered by a nanosize tip sliding on top of a supported multiwall CNT along a direction parallel or transverse to the CNT axis. Surprisingly, we find a higher friction coefficient in the transverse direction compared with the parallel direction. This behaviour is explained by a simulation showing that transverse friction elicits a soft 'hindered rolling' of the tube and a frictional dissipation that is absent, or partially absent for chiral CNTs, when the tip slides parallel to the CNT axis. Our findings can help in developing better strategies for large-scale CNT assembling and sorting on a surface.

preprint2012arXiv

Orthorhombic fulleride (CH3NH2)K3C60 close to Mott-Hubbard instability: Ab initio study

We study the electronic structure and magnetic interactions in methylamine-intercalated orthorhombic alkali-doped fullerene (CH3NH2)K3C60 within the density functional theory. As in the simpler ammonia intercalated compound (NH3)K3C60, the orthorhombic crystal-field anisotropy Δlifts the t1u triple degeneracy at the Γpoint and drives the system deep into the Mott-insulating phase. However, the computed Δand conduction electron bandwidth W cannot alone account for the abnormally low experimental Néel temperature, T_N = 11 K of the methylamine compound, compared to the much higher value T_N = 40 K of the ammonia one. Significant interactions between CH3NH2 and C60^{3-} are responsible for the stabilization of particular pseudo-Jahn-Teller fullerene-cage distortions and the ensuing low-spin S = 1/2 state. These interactions also seem to affect the magnetic properties, as interfullerene exchange interactions depend on the relative orientation of pseudo-Jahn-Teller distortions of neighboring C60^{3-} molecules. For the ferro-orientational order of CH3NH2-K^+ groups we find an apparent reduced dimensionality in magnetic exchange interactions, which may explain the suppressed Néel temperature. The disorder in exchange interactions caused by orientational disorder of CH3NH2-K^+ groups could further contribute to this suppression.

preprint2012arXiv

Static and dynamic friction in sliding colloidal monolayers

In a pioneer experiment, Bohlein et al. realized the controlled sliding of two-dimensional colloidal crystals over laser-generated periodic or quasi-periodic potentials. Here we present realistic simulations and arguments which besides reproducing the main experimentally observed features, give a first theoretical demonstration of the potential impact of colloid sliding in nanotribology. The free motion of solitons and antisolitons in the sliding of hard incommensurate crystals is contrasted with the soliton-antisoliton pair nucleation at the large static friction threshold Fs when the two lattices are commensurate and pinned. The frictional work directly extracted from particles' velocities can be analysed as a function of classic tribological parameters, including speed, spacing and amplitude of the periodic potential (representing respectively the mismatch of the sliding interface, and the corrugation, or "load"). These and other features suggestive of further experiments and insights promote colloid sliding to a novel friction study instrument.

preprint2012arXiv

Systematic Improvement of Classical Nucleation Theory

We reconsider the applicability of classical nucleation theory (CNT) to the calculation of the free energy of solid cluster formation in a liquid and its use to the evaluation of interface free energies from nucleation barriers. Using two different freezing transitions (hard spheres and NaCl) as test cases, we first observe that the interface-free-energy estimates based on CNT are generally in error. As successive refinements of nucleation-barrier theory, we consider corrections due to a non-sharp solid-liquid interface and to a non-spherical cluster shape. Extensive calculations for the Ising model show that corrections due to a non-sharp and thermally fluctuating interface account for the barrier shape with excellent accuracy. The experimental solid nucleation rates that are measured in colloids are better accounted for by these non-CNT terms, whose effect appears to be crucial in the interpretation of data and in the extraction of the interface tension from them.

preprint2011arXiv

Crossover from adiabatic to antiadiabatic quantum pumping with dissipation

Quantum pumping, in its different forms, is attracting attention from different fields, from fundamental quantum mechanics, to nanotechnology, to superconductivity. We investigate the crossover of quantum pumping from the adiabatic to the anti-adiabatic regime in the presence of dissipation, and find general and explicit analytical expressions for the pumped current in a minimal model describing a system with the topology of a ring forced by a periodic modulation of frequency omega. The solution allows following in a transparent way the evolution of pumped DC current from much smaller to much larger omega values than the other relevant energy scale, the energy splitting introduced by the modulation. We find and characterize a temperature-dependent optimal value of the frequency for which the pumped current is maximal.

preprint2011arXiv

Finite temperature properties of clusters by replica exchange metadynamics: the water nonamer

We introduce an approach for the accurate calculation of thermal properties of classical nanoclusters. Based on a recently developed enhanced sampling technique, replica exchange metadynamics, the method yields the true free energy of each relevant cluster structure, directly sampling its basin and measuring its occupancy in full equilibrium. All entropy sources, whether vibrational, rotational anharmonic and especially configurational -- the latter often forgotten in many cluster studies -- are automatically included. For the present demonstration we choose the water nonamer (H2O)9, an extremely simple cluster which nonetheless displays a sufficient complexity and interesting physics in its relevant structure spectrum. Within a standard TIP4P potential description of water, we find that the nonamer second relevant structure possesses a higher configurational entropy than the first, so that the two free energies surprisingly cross for increasing temperature.

preprint2011arXiv

Nanofriction in Cold Ion Traps

Sliding friction between crystal lattices and the physics of cold ion traps are so far non-overlapping fields. Two sliding lattices may either stick and show static friction or slip with dynamic friction; cold ions are known to form static chains, helices, or clusters, depending on trapping conditions. Here we show, based on simulations, that much could be learnt about friction by sliding, via e.g. an electric field, the trapped ion chains over a periodic corrugated potential. Unlike infinite chains where, according to theory, the classic Aubry transition to free sliding may take place, trapped chains are always pinned. Nonetheless we find that a properly defined static friction still vanishes Aubry-like at a symmetric-asymmetric structural transition, ubiquitous for decreasing corrugation in both straight and zig-zag trapped chains. Dynamic friction can also be addressed by ringdown oscillations of the ion trap. Long theorized static and dynamic one dimensional friction phenomena could thus become exquisitely accessible in future cold ion tribology.

preprint2011arXiv

Sliding Over a Phase Transition

The effects of a displacive structural phase transition on sliding friction are in principle accessible to nanoscale tools such as the Atomic Force Microscopy, yet they are still surprisingly unexplored. We present model simulations demonstrating and clarifying the mechanism and potential impact of these effects. A structural order parameter inside the material will yield a contribution to stick-slip friction that is nonmonotonic as temperature crosses the phase transition, peaking at the critical Tc where critical fluctuations are strongest, and the sliding-induced order parameter local flips from one value to another more numerous. Accordingly, the friction below Tc is larger when the order parameter orientation is such that flips are more effectively triggered by the slider. The observability of these effects and their use for friction control are discussed, for future application to sliding on the surface of and ferro- or antiferro-distortive materials.

preprint2010arXiv

AFM Dissipation Topography of Soliton Superstructures in Adsorbed Overlayers

In the atomic force microscope, the nanoscale force topography of even complex surface superstructures is extracted by the changing vibration frequency of a scanning tip. An alternative dissipation topography with similar or even better contrast has been demonstrated recently by mapping the (x,y)-dependent tip damping but the detailed damping mechanism is still unknown. Here we identify two different tip dissipation mechanisms: local mechanical softness and hysteresis. Motivated by recent data, we describe both of them in a onedimensional model of Moire' superstructures of incommensurate overlayers. Local softness at "soliton" defects yields a dissipation contrast that can be much larger than the corresponding density or corrugation contrast. At realistically low vibration frequencies, however, a much stronger and more effective dissipation is caused by the tip-induced nonlinear jumping of the soliton, naturally developing bistability and hysteresis. Signatures of this mechanism are proposed for experimental identification.

preprint2010arXiv

Atomic spin sensitive dissipation on magnetic surfaces

We identify the mechanism of energy dissipation relevant to spin-sensitive nanomechanics including the recently introduced magnetic exchange force microscopy, where oscillating magnetic tips approach surface atomic spins. The tip-surface exchange couples spin and atom coordinates, leading to a spin-phonon problem with Caldeira-Leggett type dissipation. In the overdamped regime, that can lead to a hysteretic flip of the local spin with a large spin-dependent dissipation, even down to the very low experimental tip oscillation frequencies, describing recent observations for Fe tips on NiO. A phase transition to an underdamped regime with dramatic drop of magnetic tip dissipation should in principle be possible by tuning tip-surface distance.

preprint2010arXiv

Ballistic nanofriction

Sliding parts in nanosystems such as Nano ElectroMechanical Systems (NEMS) and nanomotors, increasingly involve large speeds, and rotations as well as translations of the moving surfaces; yet, the physics of high speed nanoscale friction is so far unexplored. Here, by simulating the motion of drifting and of kicked Au clusters on graphite - a workhorse system of experimental relevance -- we demonstrate and characterize a novel "ballistic" friction regime at high speed, separate from drift at low speed. The temperature dependence of the cluster slip distance and time, measuring friction, is opposite in these two regimes, consistent with theory. Crucial to both regimes is the interplay of rotations and translations, shown to be correlated in slow drift but anticorrelated in fast sliding. Despite these differences, we find the velocity dependence of ballistic friction to be, like drift, viscous.

preprint2010arXiv

Parameter-free dissipation in simulated sliding friction

Non-equilibrium molecular dynamics simulations, of crucial importance in sliding friction, are hampered by arbitrariness and uncertainties in the way Joule heat is removed. We implement in a realistic frictional simulation a parameter-free, non-markovian, stochastic dynamics, which, as expected from theory, absorbs Joule heat precisely as a semi-infinite harmonic substrate would. Simulating stick-slip friction of a slider over a 2D Lennard-Jones solid, we compare our virtually exact frictional results with approximate ones from commonly adopted empirical dissipation schemes. While the latter are generally in serious error, we show that the exact results can be closely reproduced by a viscous Langevin dissipation at the boundary layer, once the back-reflected frictional energy is variationally optimized.

preprint2009arXiv

Kondo Conductance in an Atomic Nanocontact from First Principles

The electrical conductance of atomic metal contacts represents a powerful tool to detect nanomagnetism. Conductance reflects magnetism through anomalies at zero bias -- generally with Fano lineshapes -- due to the Kondo screening of the magnetic impurity bridging the contact. A full atomic-level understanding of this nutshell many-body system is of the greatest importance, especially in view of our increasing need to control nanocurrents by means of magnetism. Disappointingly, zero bias conductance anomalies are not presently calculable from atomistic scratch. In this Letter we demonstrate a working route connecting approximately but quantitatively density functional theory (DFT) and numerical renormalization group (NRG) approaches and leading to a first-principles conductance calculation for a nanocontact, exemplified by a Ni impurity in a Au nanowire. A Fano-like conductance lineshape is obtained microscopically, and shown to be controlled by the impurity s-level position. We also find a relationship between conductance anomaly and geometry, and uncover the possibility of opposite antiferromagnetic and ferromagnetic Kondo screening -- the latter exhibiting a totally different and unexplored zero bias anomaly. The present matching method between DFT and NRG should permit the quantitative understanding and exploration of this larger variety of Kondo phenomena at more general magnetic nanocontacts.

preprint2009arXiv

Strongly correlated metal interfaces in the Gutzwiller approximation

We study the effect of spatial inhomogeneity on the physics of a strongly correlated electron system exhibiting a metallic phase and a Mott insulating phase, represented by the simple Hubbard model. In three dimensions, we consider various geometries, including vacuum-metal-vacuum, a junction between a weakly and a strongly correlated metal, and finally the double junctions metal-Mott insulator-metal and metal-strongly correlated metal- metal. We applied to these problems the self-consistent Gutzwiller technique recently developed in our group, whose approximate nature is compensated by an extreme flexibility,ability to treat very large systems, and physical transparency. The main general result is a clear characterization of the position dependent metallic quasiparticle spectral weight. Its behavior at interfaces reveals the ubiquitous presence of exponential decays and crossovers, with decay lengths of clear physical significance. The decay length of metallic strength in a weakly-strongly correlated metal interface is due to poor screening in the strongly correlated side. The decay length of metallic strength from a metal into a Mott insulator (or into vacuum) is due to tunneling. In both cases, the decay length is a bulk property, and diverges with a critical exponent ($\sim 1/2$ in the present approximation, mean field in character) as the (continuous, paramagnetic) Mott transition is approached.

preprint2009arXiv

Tribology of the lubricant quantized-sliding state

In the framework of Langevin dynamics, we demonstrate clear evidence of the peculiar quantized sliding state, previously found in a simple 1D boundary lubricated model [Phys. Rev. Lett. 97, 056101 (2006)], for a substantially less idealized 2D description of a confined multi-layer solid lubricant under shear. This dynamical state, marked by a nontrivial ``quantized'' ratio of the averaged lubricant center-of-mass velocity to the externally imposed sliding speed, is recovered, and shown to be robust against the effects of thermal fluctuations, quenched disorder in the confining substrates, and over a wide range of loading forces. The lubricant softness, setting the width of the propagating solitonic structures, is found to play a major role in promoting in-registry commensurate regions beneficial to this quantized sliding. By evaluating the force instantaneously exerted on the top plate, we find that this quantized sliding represents a dynamical ``pinned'' state, characterized by significantly low values of the kinetic friction. While the quantized sliding occurs due to solitons being driven gently, the transition to ordinary unpinned sliding regimes can involve lubricant melting due to large shear-induced Joule heating, for example at large speed.

preprint2008arXiv

Magnetic tunnel junctions with ferroelectric barriers: Prediction of four resistance states from first-principles

Magnetic tunnel junctions (MTJs), composed of two ferromagnetic electrodes separated by a thin insulating barrier layer, are currently used in spintronic devices, such as magnetic sensors and magnetic random access memories. Recently, driven by demonstrations of ferroelectricity at the nanoscale, thin-film ferroelectric barriers were proposed to extend the functionality of MTJs. Due to the sensitivity of conductance to the magnetization alignment of the electrodes (tunnelling magnetoresistance) and the polarization orientation in the ferroelectric barrier (tunnelling electroresistance), these multiferroic tunnel junctions (MFTJs) may serve as four-state resistance devices. Based on first-principles calculations we demonstrate four resistance states in SrRuO3/BaTiO3/SrRuO3 MFTJs with asymmetric interfaces. We find that the resistance of such a MFTJ is significantly changed when the electric polarization of the barrier is reversed and/or when the magnetizations of the electrodes are switched from parallel to antiparallel. These results reveal the exciting prospects of MFTJs for application as multifunctional spintronic devices.

preprint2008arXiv

Role of transverse displacements for a quantized-velocity state of the lubricant

Within the idealized scheme of a 1-dimensional Frenkel-Kontorova-like model, a special "quantized" sliding state was found for a solid lubricant confined between two periodic layers [PRL 97, 056101 (2006)]. This state, characterized by a nontrivial geometrically fixed ratio of the mean lubricant drift velocity <v_cm> and the externally imposed translational velocity v_ext, was understood as due to the kinks (or solitons), formed by the lubricant due to incommensuracy with one of the substrates, pinning to the other sliding substrate. A quantized sliding state of the same nature is demonstrated here for a substantially less idealized 2-dimensional model, where atoms are allowed to move perpendicularly to the sliding direction and interact via Lennard-Jones potentials. Clear evidence for quantized sliding at finite temperature is provided, even with a confined solid lubricant composed of multiple (up to 6) lubricant layers. Characteristic backward lubricant motion produced by the presence of "anti-kinks" is also shown in this more realistic context.

preprint2008arXiv

The role of lubricant molecular shape in microscopic friction

With the help of a simple two-dimensional model we simulate the tribological properties of a thin lubricant film consisting of linear (chain) molecules in the ordinary soft-lubricant regime. We find that friction generally increases with chain length, in agreement with their larger bulk viscosity. When comparing the tribological properties of molecules which stick bodily to the substrates with others carrying a single sticking termination, we find that the latter generally yield a larger friction than the former.

preprint2007arXiv

Lubricated friction between incommensurate substrates

This paper is part of a study of the frictional dynamics of a confined solid lubricant film - modelled as a one-dimensional chain of interacting particles confined between two ideally incommensurate substrates, one of which is driven relative to the other through an attached spring moving at constant velocity. This model system is characterized by three inherent length scales; depending on the precise choice of incommensurability among them it displays a strikingly different tribological behavior. Contrary to two length-scale systems such as the standard Frenkel-Kontorova (FK) model, for large chain stiffness one finds that here the most favorable (lowest friction) sliding regime is achieved by chain-substrate incommensurabilities belonging to the class of non-quadratic irrational numbers (e.g., the spiral mean). The well-known golden mean (quadratic) incommensurability which slides best in the standard FK model shows instead higher kinetic-friction values. The underlying reason lies in the pinning properties of the lattice of solitons formed by the chain with the substrate having the closest periodicity, with the other slider.

preprint2007arXiv

Velocity plateaus and jumps in carbon nanotube sliding

The friction between concentric carbon nanotubes sliding one inside the other has been widely studied and simulated, but not so far using external force as the driving variable. Our molecular dynamics (MD) simulations show that as the pulling force grows, the sliding velocity increases by jumps and plateaus rather than continuously as expected. Dramatic friction peaks (similar to that recently noted by Tangney {\it et al.} in Phys. Rev. Lett. 97 (2006) 195901) which develop around some preferential sliding velocities, are at the origin of this phenomenon. The (stable) rising edge of the peak produces a velocity plateau; the (unstable) dropping edge produces a jump to the nearest stable branch. The outcome is reminiscent of conduction in ionized gases, the plateau correspon ding to a current stabilization against voltage variations, the jump corresponding to a discharge or breakdown.

preprint2006arXiv

Hysteresis from dynamically pinned sliding states

We report a surprising hysteretic behavior in the dynamics of a simple one-dimensional nonlinear model inspired by the tribological problem of two sliding surfaces with a thin solid lubricant layer in between. In particular, we consider the frictional dynamics of a harmonic chain confined between two rigid incommensurate substrates which slide with a fixed relative velocity. This system was previously found, by explicit solution of the equations of motion, to possess plateaus in parameter space exhibiting a remarkable quantization of the chain center-of-mass velocity (dynamic pinning) solely determined by the interface incommensurability. Starting now from this quantized sliding state, in the underdamped regime of motion and in analogy to what ordinarily happens for static friction, the dynamics exhibits a large hysteresis under the action of an additional external driving force F_ext. A critical threshold value F_c of the adiabatically applied force F_ext is required in order to alter the robust dynamics of the plateau attractor. When the applied force is decreased and removed, the system can jump to intermediate sliding regimes (a sort of ``dynamic'' stick-slip motion) and eventually returns to the quantized sliding state at a much lower value of F_ext. On the contrary no hysteretic behavior is observed as a function of the external driving velocity.

preprint2006arXiv

Kink plateau dynamics in finite-size lubricant chains

We extend the study of velocity quantization phenomena recently found in the classical motion of an idealized 1D model solid lubricant -- consisting of a harmonic chain interposed between two periodic sliding potentials [Phys. Rev. Lett. 97, 056101 (2006)]. This quantization is due to one slider rigidly dragging the commensurate lattice of kinks that the chain forms with the other slider. In this follow-up work we consider finite-size chains rather than infinite chains. The finite size (i) permits the development of robust velocity plateaus as a function of the lubricant stiffness, and (ii) allows an overall chain-length re-adjustment which spontaneously promotes single-particle periodic oscillations. These periodic oscillations replace the quasi-periodic motion produced by general incommensurate periods of the sliders and the lubricant in the infinite-size model. Possible consequences of these results for some real systems are discussed.

preprint2002arXiv

Coulomb couplings in positively charged fullerene

We compute, based on density-functional electronic-structure calculations, the Coulomb couplings in the h_u highest occupied orbital of molecular C60. We obtain a multiplet-averaged Hubbard U ~ 3 eV, and four Hund-rule-like intra-molecular multiplet-splitting terms, each of the order of few hundreds of meVs. According to these couplings, all C60^n+ ions should possess a high-spin ground state if kept in their rigid, undistorted form. Even after molecular distortions are allowed, however, the Coulomb terms still appear to be somewhat stronger than the previously calculated Jahn-Teller couplings, the latter favoring low-spin states. Thus for example in C60^2+, unlike C60^2-, the balance between Hund rule and Jahn Teller yields, even if marginally, a high-spin ground state. That seems surprising in view of reports of superconductivity in field-doped C60^n+ systems.

preprint1998arXiv

Non-Crystalline Structures of Ultra-Thin Unsupported Nanowires

Computer simulations suggest that ultrathin metal wires should develop exotic, non-crystalline stable atomic structures, once their diameter decreases below a critical size of the order of a few atomic spacings. The new structures, whose details depend upon the material and the wire thickness, may be dominated by icosahedral packings. Helical, spiral-structured wires with multi-atom pitches are also predicted. The phenomenon, analogous to the appearance of icosahedral and other non-crystalline shapes in small clusters, can be rationalized in terms of surface energy anisotropy and optimal packing.

preprint1995arXiv

Low lying excitations of the Dynamical Jahn-Teller ions C_60^- and C_60^2-

We compute the vibronic spectrum of gas phase fullerene negative ions C_{60}^{-} and C_{60}^{2-}. We treat accurately the linear dynamical Jahn-Teller of these ions in their ground state. In particular, coupling of the t_{1u} orbital to the eight H_g vibrational modes is handled by exact diagonalizations including up to six vibrons. The resulting spectrum is characterized by large splittings, which should be readily observable spectroscopically. The lowest excitation symmetry is argued to be related to the Berry phase present in this problem.

preprint1995arXiv

Phase Diagram of a Model of Correlated Hopping of Electrons in a Lattice of Berry Molecules

The $1D$ phase diagram of a model for correlated hopping of electrons in a lattice of Berry phase molecules is presented. Electrons hop in presence of an extra orbital degree of freedom at each site. This is mimicked as a spin-1 variable whose allowed states depend on the electron occupancy so as to take into account the orbital degeneracies of different molecular occupancies. In the $1D$ case we find that at low electron densities $n<<1$ there is a region with dominant superconducting correlations surviving an additional repulsive on-site interaction $U$ as strong as the bandwidth, $W=4t$. The critical value $U_c$ of $U$ below which superconductivity is found to be dominant decreases with increasing density $n$. For $n=1/2$ we find $U_c/t\approx 1$, whereas at $n=1$ (half-filling) our (less accurate) results are compatible with $U_c/t\approx 0$. For $U>U_c(n)$ and away from half-filling ($n\neq 1$) the system is metallic with dominant $2k_F$ charge density wave (CDW) correlations. At half-filling a charge gap opens for $U>U_c$ and the system becomes an insulator. A spin-gap characterizes the phase-diagram for all densities and for all values of $U$, even in the metallic regime $U>U_c$.

preprint1994arXiv

Anomalous low-energy electron attachment in C_{60}

We propose that thermal electron attachment to C$_{60}$ should occur preferentially in the p-wave channel, following an analysis of the vibron excitation spectrum of C$_{60}^-$. A very simple model based on this idea is shown to account very well for recent attachment data. The unexplained activation energy of $\approx$ 0.26 eV found experimentally is attributed to the p-wave centrifugal barrier.