Search

We investigated the electronic structure of the Si(111)--7$\times$7 surface below 20 K by scanning tunneling and photoemission spectroscopies and by density functional theory calculations. Previous experimental studies have questioned the ground state of this surface, which is expected to be metallic in a band picture because of the odd number of electrons per unit cell. Our differential conductance spectra instead show the opening of an energy gap at the Fermi level and a significant temperature dependence of the electronic properties, especially for the adatoms at the center of the unfaulted half of the unit cell. Complementary photoemission spectra with improved correction of the surface photovoltage shift corroborate the differential conductance data and demonstrate the absence of surface bands crossing the Fermi level at 17 K. These consistent experimental observations point to an insulating ground state and contradict the prediction of a metallic surface obtained by density functional theory in the generalized gradient approximation. The calculations indicate that this surface has or is near a magnetic instability, but remains metallic in the magnetic phases even including correlation effects at mean-field level. We discuss possible origins of the observed discrepancies between experiments and calculations.

The surprising insulating and superconducting states of narrow-band graphene twisted bilayers have been mostly discussed so far in terms of strong electron correlation, with little or no attention to phonons and electron-phonon effects. We found that, among the 33492 phonons of a fully relaxed $\theta=1.08^\circ$ twisted bilayer, there are few special, hard and nearly dispersionless modes that resemble global vibrations of the moir\'e supercell, as if it were a single, ultralarge molecule. One of them, doubly degenerate at $\Gamma$ with symmetry $A_1+B_1$, couples very strongly with the valley degrees of freedom, also doubly degenerate, realizing a so-called $\text{E}\otimes\text{e}$ Jahn-Teller (JT) coupling. The JT coupling lifts very efficiently all degeneracies which arise from the valley symmetry, and may lead, for an average atomic displacement as small as $0.5~$mA, to an insulating state at charge neutrality. This insulator possesses a non-trivial topology testified by the odd winding of the Wilson loop. In addition, freezing the same phonon at a zone boundary point brings about insulating states at most integer occupancies of the four ultra-flat electronic bands. Following that line, we further study the properties of the superconducting state that might be stabilized by these modes. Since the JT coupling modulates the hopping between AB and BA stacked regions, pairing occurs in the spin-singlet Cooper channel at the inter-(AB-BA) scale, which may condense a superconducting order parameter in the extended $s$-wave and/or $d\pm id$-wave symmetry., Comment: Revised version published in PRX (Open access)

We present a tight-binding calculation of a twisted bilayer graphene at magic angle $\theta\sim 1.08^\circ$, allowing for full, in- and out-of-plane, relaxation of the atomic positions. The resulting band structure displays as usual four narrow mini bands around the neutrality point, well separated from all other bands after the lattice relaxation. A thorough analysis of the mini-bands Bloch functions reveals an emergent $D_6$ symmetry, despite the lack of any manifest point group symmetry in the relaxed lattice. The Bloch functions at the $\Gamma$ point are degenerate in pairs, reflecting the so-called valley degeneracy. Moreover, each of them is invariant under C$_{3z}$, i.e., transforming like one-dimensional, in-plane symmetric irreducible representation of an "emergent" $D_6$ group. Out of plane, the lower doublet is even under C$_{2x}$, while the upper doublet is odd, which implies that at least eight Wannier orbitals, two $s$-like and two $p_z$-like for each of the two supercell sublattices AB and BA are necessary, probably not sufficient, to describe the four mini bands. This unexpected one-electron complexity is likely to play an important role in the still unexplained metal-insulator-superconductor phenomenology of this system., Comment: 9 pages, 8 figures, 1 table

The possibility to achieve entirely frictionless, i.e. superlubric, sliding between solids, holds enormous potential for the operation of mechanical devices. At small length scales, where mechanical contacts are well-defined, Aubry predicted a transition from a superlubric to a pinned state when the mechanical load is increased. Evidence for this intriguing Aubry transition (AT), which should occur in one dimension (1D) and at zero temperature, was recently obtained in few-atom chains. Here, we experimentally and theoretically demonstrate the occurrence of the AT in an extended two-dimensional (2D) system at room temperature using a colloidal monolayer on an optical lattice. Unlike the continuous nature of the AT in 1D, we observe a first-order transition in 2D leading to a coexistence regime of pinned and unpinned areas. Our data demonstrate that the original concept of Aubry does not only survive in 2D but is relevant for the design of nanoscopic machines and devices at ambient temperature., Comment: 12 pages including 4 figures + 9 pages supplemental information

Tribology, which studies surfaces in contact and relative motion, includes friction, wear, and lubrication, straddling across different fields: mechanical engineering, materials science, chemistry, nanoscience, physics. This short review restricts to the last two disciplines, with a qualitative survey of a small number of recent progress areas in the physics of nanofriction., Comment: 17 pages, 7 figures

Non-equilibrium Markov State Modeling (MSM) has recently been proposed [Phys. Rev. E 94, 053001 (2016)] as a possible route to construct a physical theory of sliding friction from a long steady state atomistic simulation: the approach builds a small set of collective variables, which obey a transition-matrix based equation of motion, faithfully describing the slow motions of the system. A crucial question is whether this approach can be extended from the original 1D small size demo to larger and more realistic size systems, without an inordinate increase of the number and complexity of the collective variables. Here we present a direct application of the MSM scheme to the sliding of an island made of over 1000 harmonically bound particles over a 2D periodic potential. Based on a totally unprejudiced phase space metric and without requiring any special doctoring, we find that here too the scheme allows extracting a very small number of slow variables, necessary and sufficient to describe the dynamics of island sliding., Comment: 6 pages, 5 figures

The nonlinear dynamics associated with sliding friction forms a broad interdisciplinary research field that involves complex dynamical processes and patterns covering a broad range of time and length scales. Progress in experimental techniques and computational resources has stimulated the development of more refined and accurate mathematical and numerical models, capable of capturing many of the essentially nonlinear phenomena involved in friction., Comment: 41 pages, 13 figures

Inspired by spring-block models, we elaborate a "minimal" physical model of earthquakes which reproduces two main empirical seismological laws, the Gutenberg-Richter law and the Omori aftershock law. Our new point is to demonstrate that the simultaneous incorporation of ageing of contacts in the sliding interface and of elasticity of the sliding plates constitute the minimal ingredients to account for both laws within the same frictional model.

Markov State Modeling has recently emerged as a key technique for analyzing rare events in thermal equilibrium molecular simulations and finding metastable states. Here we export this technique to the study of friction, where strongly non-equilibrium events are induced by an external force. The approach is benchmarked on the well-studied Frenkel-Kontorova model, where we demonstrate the unprejudiced identification of the minimal basis microscopic states necessary for describing sliding, stick-slip and dissipation. The steps necessary for the application to realistic frictional systems are highlighted., Comment: 9 pages, 11 figures

The quantum well states of a film can be used to sample the electronic structure of the parent bulk material and determine its band parameters. We highlight the benefits of two-dimensional film band mapping, with respect to complex bulk analysis, in an angle-resolved photoemission spectroscopy study of the 5d states of Au(111). Discrete 5d-derived quantum well states of various orbital characters form in Au(111) films and span the width of the corresponding bulk bands. For sufficiently thick films, the dispersion of these states samples the bulk band edges, as confirmed by first-principles calculations, thus providing the positions of the critical points of bulk Au in agreement with previously determined values. In turn, this analysis identifies several d-like surface states and resonances with large spin-splittings that originate from the strong spin-orbit coupling of the Au 5d atomic levels.

The Kondo zero bias anomaly of Co adatoms probed by scanning tunneling microscopy is known to depend on the height of the tip above the surface, and this dependence is different on different low index Cu surfaces. On the (100) surface, the Kondo temperature first decreases then increases as the tip approaches the adatom, while on the (111) surface it is virtually unaffected. These trends are captured by combined density functional theory and numerical renormalization group (DFT+NRG) calculations. The adatoms are found to be described by an S = 1 Anderson model on both surfaces, and ab initio calculations help identify the symmetry of the active d orbitals. We correctly reproduce the Fano lineshape of the zero bias anomaly for Co/Cu(100) in the tunneling regime but not in the contact regime, where it is probably dependent on the details of the tip and contact geometry. The lineshape for Co/Cu(111) is presumably affected by the presence of surface states, which are not included in our method. We also discuss the role of symmetry, which is preserved in our model scattering geometry but most likely broken in experimental conditions.

The critical fluctuations at second order structural transitions in a bulk crystal may affect the dissipation of mechanical probes even if completely external to the crystal surface. Here we show that noncontact force microscope dissipation bears clear evidence of the antiferrodistortive phase transition of SrTiO3, known for a long time to exhibit a unique, extremely narrow neutron scattering "central peak". The noncontact geometry suggests a central peak linear response coupling connected with strain. The detailed temperature dependence reveals for the first time the intrinsic central peak width of order 80 kHz, two orders of magnitude below the established neutron upper bound., Comment: 5 pages, 3 figures

Electrical charging of parallel plates confining a model ionic liquid down to nanoscale distances yields a variety of charge-induced changes in the structural features of the confined film. That includes even-odd switching of the structural layering and charging-induced solidification and melting, with important changes of local ordering between and within layers, and of squeezout behavior. By means of molecular dynamics simulations, we explore this variety of phenomena in the simplest charged Lennard-Jones coarse-grained model including or excluding the effect a neutral tail giving an anisotropic shape to one of the model ions. Using these models and open conditions permitting the flow of ions in and out of the interplate gap, we simulate the liquid squeezout to obtain the distance dependent structure and forces between the plates during their adiabatic appraoch under load. Simulations at fixed applied force illustrate an effective electrical pumping of the ionic liquid, from a thick nearly solid film that withstands the interplate pressure for high plate charge to complete squeezout following melting near zero charge. Effective enthalpy curves obtained by integration of interplate forces versus distance show the local minima that correspond to layering, and predict the switching between one minimum and another under squeezing and charging., Comment: submitted to J. Chem. Phys

Oxygen, one of the most common and important elements in nature, has an exceedingly well explored phase diagram under pressure, up and beyond 100 GPa. At low temperatures, the low pressures antiferromagnetic phases below 8 GPa where O$_2$ molecules have spin S=1 are followed by the broad apparently nonmagnetic $\epsilon$ phase from about 8 to 96 GPa. In this phase which is our focus molecules group structurally together to form quartets while switching, as believed by most, to spin S=0. Here we present theoretical results strongly connecting with existing vibrational and optical evidence, showing that this is true only above 20 GPa, whereas the S=1 molecular state survives up to at about 20 GPa. The $\epsilon$ phase thus breaks up into two: a spinless $\epsilon_0$ (20-96 GPa), and another $\epsilon_1$ (8-20 GPa) where the molecules have S=1 but possess only short range antiferromagnetic correlations. A local spin liquid-like singlet ground state akin to some earlier proposals and whose optical signature we identify in existing data, is proposed for this phase. Our proposed phase diagram thus has a first order phase transition just above 20 GPa, extending at finite temperature and most likely terminating into a crossover with a critical point near 30 GPa and 200 K.

Recent highly idealized model studies of lubricated nanofriction for two crystalline sliding surfaces with an interposed thin solid crystalline lubricant layer showed that the overall relative velocity of the lubricant $v_{\rm lub} / v_{\rm slider}$ depends only on the ratio of the lattice spacings, and retains a strictly constant value even when system parameters are varied within a wide range. This peculiar "quantized" dynamical locking was understood as due to the sliding-induced motion of misfit dislocations, or soliton structures. So far, the practical relevance of this concept to realistic sliding three dimensional crystals has not been demonstrated. In this work, by means of classical molecular dynamics simulations and theoretical considerations, we realize a realistic three-dimensional crystal-lubricant-crystal geometry. Results show that the flux of lubricant particles associated to the advancing soliton lines gives rise here too to a quantized velocity ratio. Moreover, depending on the interface lattice spacing mismatch, both forward and backward quantized motion of the lubricant is predicted. The persistence under realistic conditions of the dynamically pinned state and quantized sliding is further investigated by varying sliding speed, temperature, load, and lubricant film thickness. The possibilities of experimental observation of quantized sliding are also discussed., Comment: 14 pages, 16 figures

Low-temperature electronic conductance in nanocontacts, scanning tunneling microscopy (STM), and metal break junctions involving magnetic atoms or molecules is a growing area with important unsolved theoretical problems. While the detailed relationship between contact geometry and electronic structure requires a quantitative ab initio approach such as density functional theory (DFT), the Kondo many-body effects ensuing from the coupling of the impurity spin with metal electrons are most properly addressed by formulating a generalized Anderson impurity model to be solved with, for example, the numerical renormalization group (NRG) method. Since there is at present no seamless scheme that can accurately carry out that program, we have in recent years designed a systematic method for semiquantitatively joining DFT and NRG. We apply this DFT-NRG scheme to the ideal conductance of single wall (4,4) and (8,8) nanotubes with magnetic adatoms (Co and Fe), both inside and outside the nanotube, and with a single carbon atom vacancy. A rich scenario emerges, with Kondo temperatures generally in the Kelvin range, and conductance anomalies ranging from a single channel maximum to destructive Fano interference with cancellation of two channels out of the total four. The configuration yielding the highest Kondo temperature (tens of Kelvins) and a measurable zero-bias anomaly is that of a Co or Fe impurity inside the narrowest nanotube. The single atom vacancy has a spin, but a very low Kondo temperature is predicted. The geometric, electronic, and symmetry factors influencing this variability are all accessible, which makes this approach methodologically instructive and highlights many delicate and difficult points in the first-principles modeling of the Kondo effect in nanocontacts.

The increase of sliding friction upon increasing load is a classic in the macroscopic world. Here we discuss the possibility that friction rise might sometimes turn into a drop when, at the mesoscale and nanoscale, a confined lubricant film separating crystalline sliders undergoes strong layering and solidification. Under pressure, transitions from N to N-1 layers may imply a change of lateral periodicity of the crystallized lubricant sufficient to alter the matching of crystal structures, influencing the ensuing friction jump. A pressure-induced friction drop may occur as the shear gradient maximum switches from the lubricant middle, marked by strong stick-slip with or without shear melting, to the crystalline slider-lubricant interface, characterized by smooth superlubric sliding. We present high pressure sliding simulations to display examples of frictional drops, suggesting their possible relevance to the local behavior in boundary lubrication.

Non-equilibrium molecular dynamics simulations, of crucial importance in sliding friction, are hampered by arbitrariness and uncertainties in the removal of the frictionally generated Joule heat. Building upon general pre-existing formulation, we implement a fully microscopic dissipation approach which, based on a parameter-free, non-Markovian, stochastic dynamics, absorbs Joule heat equivalently to a semi-infinite solid and harmonic substrate. As a test case, we investigate the stick-slip friction of a slider over a two-dimensional Lennard-Jones solid, comparing our virtually exact frictional results with approximate ones from commonly adopted dissipation schemes. Remarkably, the exact results can be closely reproduced by a standard Langevin dissipation scheme, once its parameters are determined according to a general and self-standing variational procedure.

The effect of magnetic impurities on the ballistic conductance of nanocontacts is, as suggested in recent work, amenable to ab initio study \cite{naturemat}. Our method proceeds via a conventional density functional calculation of spin and symmetry dependent electron scattering phase shifts, followed by the subsequent numerical renormalization group solution of Anderson models -- whose ingredients and parameters are chosen so as to reproduce these phase shifts. We apply this method to investigate the Kondo zero bias anomalies that would be caused in the ballistic conductance of perfect metallic (4,4) and (8,8) single wall carbon nanotubes, ideally connected to leads at the two ends, by externally adsorbed Co and Fe adatoms. The different spin and electronic structure of these impurities are predicted to lead to a variety of Kondo temperatures, generally well below 10 K, and to interference between channels leading to Fano-like conductance minima at zero bias.

Transition metal impurities will yield zero bias anomalies in the conductance of well contacted metallic carbon nanotubes, but Kondo temperatures and geometry dependences have not been anticipated so far. Applying the density functional plus numerical renormalization group approach of Lucignano \textit{et al.} to Co and Fe impurities in (4,4) and (8,8) nanotubes, we discover a huge difference of behaviour between outside versus inside adsorption of the impurity. The predicted Kondo temperatures and zero bias anomalies, tiny outside the nanotube, turn large and strongly radius dependent inside, owing to a change of symmetry of the magnetic orbital. Observation of this Kondo effect should open the way to a host of future experiments.

In this Letter we report a simulation study in which we compare the solid-liquid interfacial free energy of NaCl at coexistence, with the value that follows from the height of the homogeneous nucleation barrier. We find that the two estimates differ by more than 100%. Similar, although smaller discrepancies are found for crystals of hard-sphere colloids and of Lennard-Jones ("argon") particles. We consider a variety of possible causes for this discrepancy and are forced to conclude that it is due to a finite-size effect that cannot be corrected for by any simple thermodynamic procedure. Importantly, we find that the surface free energies that follow from real nucleation experiments should be subject to large finite size effects. Taking this in to account, we obtain quantitative agreement between the simulation data and the surface free energy of NaCl that follows from nucleation experiments. Our finding suggests that most published solid-liquid surface free energies derived from nucleation experiments will have to be revised.

The magnetic ground state of nanosized systems such as Fe double chains, recently shown to form in the early stages of Fe deposition on Ir(100), is generally nontrivial. Using ab initio density functional theory we find that the straight ferromagnetic (FM) state typical of bulk Fe as well as of isolated Fe chains and double chains is disfavored after deposition on Ir(100) for all the experimentally relevant double chain structures considered. So long as spin-orbit coupling (SOC) is neglected, the double chain lowest energy state is generally antiferromagnetic (AFM), a state which appears to prevail over the FM state due to Fe-Ir hybridization. Successive inclusion of SOC adds two further elements, namely a magnetocrystalline anisotropy, and a Dzyaloshinskii-Moriya (DM) spin-spin interaction, the former stabilizing the collinear AFM state, the second favoring a long-period spin modulation. We find that anisotropy is most important when the double chain is adsorbed on the partially deconstructed Ir(100) -- a state which we find to be substantially lower in energy than any reconstructed structure -- so that in this case the Fe double chain should remain collinear AFM. Alternatively, when the same Fe double chain is adsorbed in a metastable state onto the (5x1) fully reconstructed Ir(100) surface, the FM-AFM energy difference is very much reduced and the DM interaction is expected to prevail, probably yielding a helical spin structure., Comment: to appear on PRB

We discovered in simulations of sliding coaxial nanotubes an unanticipated example of dynamical symmetry breaking taking place at the nanoscale. While both nanotubes are perfectly left-right symmetric and nonchiral, a nonzero angular momentum of phonon origin appears spontaneously at a series of critical sliding velocities, in correspondence with large peaks of the sliding friction. The non-linear equations governing this phenomenon resemble the rotational instability of a forced string. However, several new elements, exquisitely "nano" appear here, with the crucial involvement of Umklapp and of sliding nanofriction., Comment: To appear in PRL

The trivalent alkali fullerides A$_3$C$_{60}$, where C$_{60}$ are a well established family of molecular superconductors. The electron pairing has s-wave symmetry and is due to standard electron-phonon coupling, in particular by Jahn-Teller intramolecular C$_{60}$ vibrations. A source of renewed interest in these systems are indications of strong electron-electron repulsion, which emerges especially in compounds where the C$_{60}$-C$_{60}$ distance is expanded. In several compounds after an initial increases T$_c$, further expansions leads to a decline of superconductivity and its eventual disappearance in favor of a Mott insulating state We theoretically study a three-orbital Hubbard model including the phonon-mediated interaction using Dynamical Mean-Field Theory, which is particularly suitable due to the local nature of all the interactions. We studied the system as a function of the ratio of intra-molecular repulsion $U$ over the electron bandwidth $W$, the increase of $U/W$ representing the main effect of lattice expansion. The phase diagram is close to that of actual materials, with a dome-shaped superconducting region preceding the Mott transition. Unconventional properties predicted by this model include: (i) a pseudogap in the normal phase; (ii) a gain of kinetic energy and of d.c. conductivity at the onset of superconductivity; (iii) regular spin susceptibility and specific heat despite strong correlations; (iv) the emergence of more than one energy scale governing the renormalized single particle dispersion. These predictions, if confirmed, would establish fullerides, especially the expanded ones, as members of the wider family of strongly correlated superconductors., Comment: 17 pages, 11 figures, to be published on Rev. Mod. Phys. (Colloquia)

Recent experiments showed that the conductance of Pt nanocontacts and nanowires is measurably reduced by adsorption of CO. We present DFT calculations of the electronic structure and ballistic conductance of a Pt monoatomic chain and a CO molecule adsorbed in an on-top position. We find that the main electronic molecule-chain interaction occurs via the $5\sigma$ and $2\pi^{\star}$ orbitals of the molecule, involved in a donation/back-donation process similar to that of CO on transition-metal surfaces. The ideal ballistic conductance of the monoatomic chain undergoes a moderate reduction by about 1.0 G_0 (from 4 G_0 to 3.1 G_0) upon adsorption of CO. By repeating all calculations with and without spin-orbit coupling, no substantial spin-orbit induced change emerges either in the chain-molecule interaction mechanism or in the conductance., Comment: 4 pages, 2 figures, in proceedings of Frontiers of Fundamental and Computational Physics

We carry out a first-principles density functional study of the interaction between a monatomic Pt wire and a CO molecule, comparing the energy of different adsorption configurations (bridge, on top, substitutional, and tilted bridge) and discussing the effects of spin-orbit (SO) coupling on the electronic structure and on the ballistic conductance of two of these systems (bridge and substitutional). We find that, when the wire is unstrained, the bridge configuration is energetically favored, while the substitutional geometry becomes possible only after the breaking of the Pt-Pt bond next to CO. The interaction can be described by a donation/back-donation process similar to that occurring when CO adsorbs on transition-metal surfaces, a picture which remains valid also in presence of SO coupling. The ballistic conductance of the (tipless) nanowire is not much reduced by the adsorption of the molecule on the bridge and on-top sites, but shows a significant drop in the substitutional case. The differences in the electronic structure due to the SO coupling influence the transmission only at energies far away from the Fermi level so that fully- and scalar-relativistic conductances do not differ significantly., Comment: 12 pages, 12 figures; figure misplacement and minor syntax issues fixed, some references updated and corrected

We consider a model inspired by a metal break-junction hypothetically caught at its breaking point, where the non-adiabatic center-of-mass motion of the bridging atom can be treated as a two-level system. By means of Numerical Renormalization Group (NRG) we calculate the influence of the two level system on the ballistic conductance across the bridge atom. The results are shown to be fully consistent with a conformal field theory treatment. We find that the conductance, calculated by coupling Fermi liquid theory to our NRG is always finite and fractional at zero temperature, but drops quite fast as the temperature increases., Comment: 9 pages 6 figures

A nonlinear model representing the tribological problem of a thin solid lubricant layer between two sliding periodic surfaces is used to analyze the phenomenon of hysteresis at pinning/depinning around a moving state rather than around a statically pinned state. The cycling of an external driving force F_ext is used as a simple means to destroy and then to recover the dynamically pinned state previously discovered for the lubricant center-of-mass velocity. De-pinning to a quasi-freely sliding state occurs either directly, with a single jump, or through a sequence of discontinuous transitions. The intermediate sliding steps are reminiscent of phase-locked states and stick-slip motion in static friction, and can be interpreted in terms of the appearance of travelling density defects in an otherwise regular arrangement of kinks. Re-pinning occurs more smoothly, through the successive disappearance of different travelling defects. The resulting bistability and multistability regions may also be explored by varying mechanical parameters other than F_ext, e.g. the sliding velocity or the corrugation amplitude of the sliders., Comment: 14 pages 10 figures

The dragging velocity of a model solid lubricant confined between sliding periodic substrates exhibits a phase transition between two regimes, respectively with quantized and with continuous lubricant center-of-mass velocity. The transition, occurring for increasing external driving force F_ext acting on the lubricant, displays a large hysteresis, and has the features of depinning transitions in static friction, only taking place on the fly. Although different in nature, this phenomenon appears isomorphic to a static Aubry depinning transition in a Frenkel-Kontorova model, the role of particles now taken by the moving kinks of the lubricant-substrate interface. We suggest a possible realization in 2D optical lattice experiments., Comment: 5 pages, 4 figures, revtex, in print in Phys. Rev. Lett

Alkali halide (100) crystal surfaces are poorly wetted by their own melt at the triple point. We carried out simulations for NaCl(100) within the well tested BMHFT model potential. Calculations of the solid-vapor, solid-liquid and liquid-vapor free energies showed that solid NaCl(100) is a non-melting surface, and explain its bad wetting in detail. The extreme stability of NaCl(100) is ideal for a study of the nanofriction in the high temperature regime, close to and even above the bulk melting temperature (T_M). Our simulations reveal in this regime two distinct and opposite phenomena for plowing and for grazing friction. We found a frictional drop close to T_M for deep ploughing and wear, but on the contrary a frictional rise for grazing, wearless sliding. For both phenomena we obtain a fresh microscopic understanding, relating the former to ``skating'' through a local liquid cloud, the latter to softening of the free substrate surface. It is argued that both phenomena, to be pursued experimentally, should be much more general than the specific NaCl surface case. Most metals in particular possessing one or more close packed non-melting surface, such as Pb, Al or Au(111), should behave quite similarly., Comment: 9 pages, 6 figures; High Temperature Capillarity 2007 (HTC-07) proceeding

We describe the effective Lorentz forces on the ions of a generic insulating system in an magnetic field, in the context of Born-Oppenheimer ab-initio molecular dynamics. The force on each ion includes an important contribution of electronic origin, which depends explicitly on the velocity of all other ions. It is formulated in terms of a Berry curvature, in a form directly suitable for future first principles classical dynamics simulations based {\it e.g.,} on density functional methods. As a preliminary analytical demonstration we present the dynamics of an H$_2$ molecule in a field of intermediate strength, approximately describing the electrons through Slater's variational wavefunction., Comment: 5 pages, 2 figures; to appear in Phys. Rev. B

The physics of sliding nanofriction at high temperature near the substrate melting point $\Tmelt$ is so far unexplored. We conducted simulations of hard tips sliding on a prototype non-melting surface, NaCl(100), revealing in this regime two distinct and opposite phenomena for plowing and for grazing friction. We found a frictional drop close to $\Tmelt$ for deep plowing and wear, but on the contrary a frictional rise for grazing, wearless sliding. For both phenomena we obtain a fresh microscopic understanding, relating the former to ``skating'' through a local liquid cloud, the latter to linear response properties of the free substrate surface. It is argued that both phenomena should be pursued experimentally, and much more general than the specific NaCl surface case. Most metals in particular possessing one or more close packed non-melting surface, such as Pb, Al or Au(111), that should behave quite similarly., Comment: 5 pages, 6 figures

The ground state of Sn/Si(111) and Sn/Ge(111) surface $\alpha$-phases is reexamined theoretically, based on $ab-initio$ calculations where correlations are approximately included through the orbital dependence of the Coulomb interaction (in the local density + Hubbard U approximation). The effect of correlations is to destabilize the vertical buckling in Sn/Ge(111) and to make the surface magnetic, with a metal-insulator transition for both systems. This signals the onset of a stable narrow gap Mott-Hubbard insulating state, in agreement with very recent experiments. Antiferromagnetic exchange is proposed to be responsible for the observed $\Gamma$-point photoemission intensity, as well asfor the partial metallization observed above above 60 K in Sn/Si(111). Extrinsic metallization of Sn/Si(111) by, $e.g.$ alkali doping, could lead to a novel 2D triangular superconducting state of this and similar surfaces., Comment: 4 pages, 4 figures

We performed first principles density functional calculations and simulations of magic-size neutral NaCl nanocubes, and computed the the extraction of a Na neutral corner atom after donating an electron. The atomic structure of the resulting Na corner vacancy is presented., Comment: 7 pages, 2 figures, 1 table; to appear on Surf. Sci

Alkali halide (100) surfaces are anomalously poorly wetted by their own melt at the triple point. We carried out simulations for NaCl(100) within a simple (BMHFT) model potential. Calculations of the solid-vapor, solid-liquid and liquid-vapor free energies showed that solid NaCl(100) is a nonmelting surface, and that the incomplete wetting can be traced to the conspiracy of three factors: surface anharmonicities stabilizing the solid surface; a large density jump causing bad liquid-solid adhesion; incipient NaCl molecular correlations destabilizing the liquid surface, reducing in particular its entropy much below that of solid NaCl(100). Presently, we are making use of the nonmelting properties of this surface to conduct case study simulations of hard tips sliding on a hot stable crystal surface. Preliminary results reveal novel phenomena whose applicability is likely of greater generality., Comment: 7 pages, 6 figure; minor typo corrections; to appear in: Surf. Sci

Electronic and transport properties of a short Ni nanowire suspended between two semi-infinite ferromagnetic Ni leads are explored in the framework of density-functional theory. The spin-dependent ballistic conductance of the nanowire is calculated using a scattering-based approach and the Landauer-Buttiker formula. The total calculated conductance in units of $G_0 = 2e^2/h$ is around 1.6, in fairly good agreement with the broad peak observed around 1.5 for the last conductance step in break junctions. Separating contributions from different spins, we find nearly 0.5 $G_0$ from the majority spin $s$-like channel, whereas the remaining minority spin conductance of 1.1 $G_0$ contains significant contributions from several $d$ states, but much less than 0.5 $G_0$ from $s$ states. The influence of the structural relaxation on the magnetic properties and the ballistic conductance of the nanowire is also studied., Comment: 9 pages, 7 figures, to appear in PRB

This paper presents a broad theoretical and simulation study of the high temperature behavior of crystalline alkali halide surfaces typified by NaCl(100), of the liquid NaCl surface near freezing, and of the very unusual partial wetting of the solid surface by the melt. Simulations are conducted using two-body rigid ion BMHFT potentials, with full treatment of long-range Coulomb forces. After a preliminary check of the description of bulk NaCl provided by these potentials, which seems generally good even at the melting point, we carry out a new investigation of solid and liquid surfaces. Solid NaCl(100) is found in this model to be very anharmonic and yet exceptionally stable when hot. It is predicted by a thermodynamic integration calculation of the surface free energy that NaCl(100) should be a well ordered, non-melting surface, metastable even well above the melting point. By contrast, the simulated liquid NaCl surface is found to exhibit large thermal fluctuations and no layering order. In spite of that, it is shown to possess a relatively large surface free energy. The latter is traced to a surface entropy deficit, reflecting some kind of surface short range order. Finally, the solid-liquid interface free energy is derived through Young's equation from direct simulation of partial wetting of NaCl(100) by a liquid droplet. It is concluded that three elements, namely the exceptional anharmonic stability of the solid (100) surface, the molecular short range order at the liquid surface, and the costly solid liquid interface, all conspire to cause the anomalously poor wetting of the (100) surface by its own melt in the BMHFT model of NaCl -- and most likely also in real alkali halide surfaces., Comment: modified version of JCP 123, 164701 15 pages, 25 figures

The high temperature surface properties of alkali halide crystals are very unusual. Through molecular dynamics simulations based on Tosi-Fumi potentials, we predict that crystalline NaCl(100) should remain stable without any precursor signals of melting up to and even above the bulk melting point $T_m$. In a metastable state, it should even be possible to overheat NaCl (100) by at least 50 K. The reasons leading to this lack of surface self-wetting are investigated. We will briefly discuss the results of calculations of the solid-vapor and liquid-vapor interface free energies, showing that the former is unusually low and the latter unusually high, and explaining why. Due to that the mutual interaction among solid-liquid and liquid-vapor interfaces, otherwise unknown, must be strongly attractive at short distance, leading to the collapse of any liquid film attempting to nucleate at the solid surface. This scenario naturally explains the large incomplete wetting angle of a drop of melt on NaCl(100)., Comment: 9 pages, 3 figures; to appear in: "Highlights in the quantum theory of condensed matter", a volume in honour of M. P. Tosi on his 72th birthday, ed. Scuola Normale Superiore di Pisa

Alkali halide (100) crystal surfaces are anomalous, being very poorly wetted by their own melt at the triple point. We present extensive simulations for NaCl, followed by calculations of the solid-vapor, solid-liquid, and liquid-vapor free energies showing that solid NaCl(100) is a nonmelting surface, and that its full behavior can quantitatively be accounted for within a simple Born-Meyer-Huggins-Fumi-Tosi model potential. The incomplete wetting is traced to the conspiracy of three factors: surface anharmonicities stabilizing the solid surface; a large density jump causing bad liquid-solid adhesion; incipient NaCl molecular correlations destabilizing the liquid surface. The latter is pursued in detail, and it is shown that surface short-range charge order acts to raise the surface tension because incipient NaCl molecular formation anomalously reduces the surface entropy of liquid NaCl much below that of solid NaCl(100)., Comment: 4 pages, 3 figures

We studied theoretically the behavior of an injected electron-hole pair in crystalline polyethylene. Time-dependent adiabatic evolution by ab-initio molecular dynamics simulations show that the pair will become self-trapped in the perfect crystal, with a trapping energy of about 0.38 eV, with formation of a pair of trans-gauche conformational defects, three C$_2$H$_4$ units apart on the same chain. The electron is confined in the inter-chain pocket created by a local, 120$^\circ$ rotation of the chain between the two defects, while the hole resides on the chain and is much less bound. Despite the large energy stored in the trapped excitation, there does not appear to be a direct non-radiative channel for electron-hole recombination. This suggests that intrinsic self-trapping of electron-hole pairs inside the ideal quasi-crystalline fraction of PE might not be directly relevant for electrical damage in high-voltage cables., Comment: 9 pages, 3 figs; using iopart.cls

We present an extensive but concise review of our present understanding, largely based on theory and simulation work from our group, on the equilibrium behavior of solid surfaces and nanosystems close to the bulk melting point. In the first part we define phenomena, in particular surface melting and nonmelting, and review some related theoretical approaches, from heuristic theories to computer simulation. In the second part we describe the surface melting/nonmelting behavior of several different classes of solids, ranging from van der Waals crystals, to valence semiconductors, to ionic crystals and metals. In the third part, we address special cases such as strained solids, the defreezing of glass surfaces, and rotational surface melting. Next, we digress briefly to surface layering of a liquid metal, possibly leading to solid-like or hexatic two dimensional phases floating on the liquid. In the final part, the relationship of surface melting to the premelting of nanoclusters and nanowires is reviewed., Comment: 54 pages, 26 figures

Rubber friction on wet rough substrates at low velocities is typically 20-30% smaller than for the corresponding dry surfaces. We show that this cannot be due to hydrodynamics and propose a novel explanation based on a sealing effect exerted by rubber on substrate "pools" filled with water. Water effectively smoothens the substrate, reducing the major friction contribution due to induced viscoelastic deformations of the rubber by surface asperities. The theory is illustrated with applications related to tire-road friction., Comment: Format Revtex 4; 8 pages, 11 figures (no color); Published on Phys. Rev. B (http://link.aps.org/abstract/PRB/v71/e035428); previous work on the same topic: cond-mat/0412045

Surface roughness has a huge impact on many important phenomena. The most important property of rough surfaces is the surface roughness power spectrum C(q). We present surface roughness power spectra of many surfaces of practical importance, obtained from the surface height profile measured using optical methods and the Atomic Force Microscope. We show how the power spectrum determines the contact area between two solids. We also present applications to sealing, rubber friction and adhesion for rough surfaces, where the power spectrum enters as an important input., Comment: Topical review; 82 pages, 61 figures; Format: Latex (iopart). Some figures are in Postscript Level 3

Loss of braking power and rubber skidding on a wet road is still an open physics problem, since neither the hydrodynamical effects nor the loss of surface adhesion that are sometimes blamed really manage to explain the 20-30% observed loss of low speed tire-road friction. Here we advance a novel mechanism based on sealing of water-filled substrate pools by the rubber. The sealed-in water effectively smoothens the substrate, thus reducing the viscoelastic dissipation in bulk rubber induced by surface asperities, well established as a major friction contribution. Starting with the measured spectrum of asperities one can calculate the water-smoothened spectrum and from that the predicted friction reduction, which is of the right magnitude. The theory is directly supported by fresh tire-asphalt friction data., Comment: 5 pages, 4 figures. Published on Nature Materials (November 7th 2004)

The self-wetting properties of ionic crystal surfaces are studied, using NaCl(100) as a prototype case. The anomalously large contact angle measured long ago by Mutaftschiev is well reproduced by realistic molecular dynamics simulations. Based on these results, and on independent determinations of the liquid-vapor and the solid-vapor interface free energy, an estimate of the solid-liquid interface free energy is extracted. The solid-vapor surface free energy turns out to be small and similar to the liquid-vapor one, providing a direct thermodynamic explanation of the reduced wetting ability of the ionic melt., Comment: 8 pages, 3 figures. Proceedings of the ECOSS22 Conference (Praha, 2003) Published on Surface Science (Elsevier)

A recent exciting experiment by Ghosh et al. reported that the flow of an ion-containing liquid such as water through bundles of single-walled carbon nanotubes induces a voltage in the nanotubes that grows logarithmically with the flow velocity v0. We propose an explanation for this observation. Assuming that the liquid molecules nearest the nanotube form a 2D solid-like monolayer pinned through the adsorbed ions to the nanotubes, the monolayer sliding will occur by elastic loading followed by local yield (stick-slip). The drifting adsorbed ions produce a voltage in the nanotube through electronic friction against free electrons inside the nanotube. Thermally excited jumps over force-biased barriers, well-known in stick-slip, can explain the logarithmic voltage growth with flow velocity. We estimate the short circuit current and the internal resistance of the nanotube voltage generator., Comment: 8 pages, 3 figures; published on PRB (http://link.aps.org/abstract/PRB/v69/e235410) and on the Virtual Journal of Nanoscale Science and Technology (http://www.vjnano.org, July 14, 2002, Vol. 10, Iss. 2)

The scattering-based approach for calculating the ballistic conductance of open quantum systems is generalized to deal with magnetic transition metals as described by ultrasoft pseudo-potentials. As an application we present quantum-mechanical conductance calculations for monatomic Co and Ni nanowires with a magnetization reversal. We find that in both Co and Ni nanowires, at the Fermi energy, the conductance of $d$ electrons is blocked by a magnetization reversal, while the $s$ states (one per spin) are perfectly transmitted. $d$ electrons have a non-vanishing transmission in a small energy window below the Fermi level. Here, transmission is larger in Ni than in Co., Comment: 9 pages, 6 figures, to appear in PRB

Near a Mott transition, strong electron correlations may enhance Cooper pairing. This is demonstrated in the Dynamical Mean Field Theory solution of a twofold-orbital degenerate Hubbard model with inverted Hund's rules on-site exchange, which favors local spin-singlet configurations. Close to the Mott insulator, which is a local version of a valence bond insulator, a pseudogap non-Fermi-liquid metal, a superconductor, and a normal metal appear, in striking similarity with the physics of the cuprates. The strongly correlated superconducting state has a larger Drude weight than the corresponding normal state. The role of the impurity Kondo problem is underscored., Comment: 4 pages, 5 figures

The approach proposed by Choi and Ihm for calculating the ballistic conductance of open quantum systems is generalized to deal with magnetic transition metals. The method has been implemented with ultrasoft pseudopotentials and plane wave basis set in a DFT-LSDA ab-initio scheme. We present the quantum-mechanical conductance calculations for monatomic Ni nanowire with a single spin reversal. We find that a spin reversal blocks the conductance of $d$ electrons at the Fermi energy of the Ni nanowire. On the other hand, two $s$ electrons (one per each spin) are perfectly transmitted in the whole energy window giving $2G_0$ for the total conductance. The relevance of these results in connection with recent experimental data is discussed., Comment: 4 pages, 1 figure, to be published in Surface Science