Search

We present a first principle study of the electron-phonon (e-p) interaction at the Be(0001) surface. The real and imaginary part of the e-p self energy are calculated for the surface state in the binding energy range from the $\bar{\Gamma}$ point to the Fermi level. Our calculation shows an overall good agreement with several photoemission data measured at high and low temperatures. Additionally, we show that the energy derivative of real part of the self-energy presents a strong temperature and energy variation close to $E_{F}$, making it difficult to measure its value just at $E_{F}$., Comment: Accepted in Phys. Rev. Lett., 5 figures

Recent experiments showed that the last, single channel conductance step in monatomic gold contacts exhibits significant fluctuations as a function of stretching. From simulations of a stretched gold nanowire linked to deformable tips, we determine the distribution of the bond lengths between atoms forming the nanocontact and analyze its influence on the electronic conductance within a simplified single channel approach. We show that the inhomogeneous distribution of bond lengths can explain the occurrence and the 5% magnitude of conductance fluctuations below the quantum conductance unit $g_{o}= 2e^2/h$.

We have investigated Pd nanowires theoretically, and found that, unlike either metallic or free atomic Pd, they exhibit Hund's rule magnetism. In long, monoatomic wires, we find a spin moment of 0.7 Bohr magnetons per atom, whereas for short, monoatomic wires between bulk leads, the predicted moment is about 0.3 Bohr magnetons per wire atom. In contrast, a coaxial (6,1) wire was found to be nonmagnetic. The origin of the wire magnetism is analyzed., Comment: 6 pages, including 4 figures

We generalize to magnetic transition metals the approach proposed by Choi and Ihm for calculating the complex band structure of periodic systems, a key ingredient for future calculations of conductivity of an open quantum system within the Landauer-Buttiker theory. The method is implemented with ultrasoft pseudopotentials and plane wave basis set in a DFT-LSDA ab-initio scheme. As a first example, we present the complex band structure of bulk fcc Ni (which constitutes the tips of a Ni nanocontact) and monatomic Ni wire (the junction between two tips). Based on our results, we anticipate some features of the spin-dependent conductance in a Ni nanocontact., Comment: 5 pages, 2 figures, to be published in "Surface Science"

The expulsion dynamics of the last liquid monolayer of molecules confined between two surfaces has been analyzed by solving the two-dimensional (2D) Navier-Stokes equation for a compressible liquid. We find that the squeeze-out is characterized by the parameter g0 ~ P0/(rho c^2), where P0 is the average perpendicular (squeezing) pressure, rho the liquid (3D) density and c the longitudinal sound velocity in the monolayer film. When g0 << 1 the result of the earlier incompressible treatment is recovered. Numerical results for the squeeze-out time, and for the time-dependence of the radius of the squeezed-out region, indicate that compressibility effects may be non-negligible both in time and in space. In space, they dominate at the edge of the squeeze-out region. In time, they are strongest right at the onset of the squeeze-out process, and just before its completion., Comment: revtex4, 6 pages, 4 figures. Published on PRB on December 31, 2002

We investigate the effects of strain on a crystal surface close to the bulk melting temperature T_m, where surface melting usually sets in. Strain lowers the bulk melting point, so that at a fixed temperature below but close to T_m the thickness of the quasi-liquid film is expected to grow with strain, irrespective of sign. In addition, a strain-induced solid surface free energy increase/decrease takes place, favoring/disfavoring surface melting depending on the sign of strain relative to surface stress. In the latter case one can produce a strain-induced prewetting transition, where for increasing temperature the liquid film suddenly jumps from zero to a finite thickness. This phenomenology is illustrated by a realistic molecular dynamics simulation of strained Al(110)., Comment: Acceped for publication on Surface Science

It has been shown in recent experiments that electronic transport through a gold monatomic nanowire is dissipative above a threshold voltage due to excitation of phonons via the electron-phonon interaction. We address that data by computing, via density functional theory, the zone boundary longitudinal phonon frequency of a perfect monatomic nanowire during its mechanical elongation. The theoretical frequency that we find for an ideally strained nanowire is not compatible with experiment if a uniformly distributed stretch is assumed. With the help of a semi-empirical Au-Au potential, we model the realistic nanowire stretching as exerted by two tips. In this model we see that strain tends to concentrate in the junctions, so that the mean strain of the nanowire is roughly one half of the ideal value. With this reduced strain, the calculated phonon softening is in much better agreement with experiment., Comment: 9 pages,3 figures, Surface Science, in press

Distortions of the $\sqrt3\times\sqrt3$ Sn/Ge(111) and Sn/Si(111) surfaces are shown to reflect a disproportionation of an integer pseudocharge, $Q$, related to the surface band occupancy. A novel understanding of the $(3\times3)$-1U (``1 up, 2 down'') and 2U (``2 up, 1 down'') distortions of Sn/Ge(111) is obtained by a theoretical study of the phase diagram under strain. Positive strain keeps the unstrained value Q=3 but removes distorsions. Negative strain attracts pseudocharge from the valence band causing first a $(3\times3)$-2U distortion (Q=4) on both Sn/Ge and Sn/Si, and eventually a $(\sqrt3\times\sqrt3)$-3U (``all up'') state with Q=6. The possibility of a fluctuating phase in unstrained Sn/Si(111) is discussed., Comment: Revtex, 5 pages, 3 figures

High temperature superconductivity in doped Mott insulators such as the cuprates contradicts the conventional wisdom that electron repulsion is detrimental to superconductivity. Because doped fullerene conductors are also strongly correlated, the recent discovery of high-critical-temperature, presumably s-wave, superconductivity in C$_{60}$ field effect devices is even more puzzling. We examine a dynamical-mean-field solution of a model for electron doped fullerenes which shows how strong correlations can indeed enhance superconductivity close to the Mott transition. We argue that the mechanism responsible for this enhancement could be common to a wider class of strongly correlated models, including those for cuprate superconductors., Comment: 5 pages, 3 figures

Based on an accurate first principles description of the energetics in H-bonded KDP, we conduct a first study of nuclear quantum effects and of the changes brought about by deuteration. Cluster tunneling involving also heavy ions is allowed, the main effect of deuteration being a depletion of the proton probability density at the O-H-O bridge center, which in turn weakens its proton-mediated covalent bonding. The ensuing lattice expansion couples selfconsistently with the proton off-centering, thus explaining both the giant isotope effect, and its close connection with geometrical effects., Comment: 4 two-column pages, 4 figures

The lowest conductance step for a Ni nanocontact is anomalously small in comparison with the large expected number of conducting channels. We present electronic structure calculations for an extremely idealized Ni nanobridge consisting of just a monatomic nanowire. Our calculations show that no less than eight single spin bands cross the Fermi level in a nonmagnetic Ni monatomic wire, dropping marginally to seven in the more stable, fully ferromagnetic state. However, when we build in the wire a magnetization reversal, or domain wall, by forcing the net magnetization to be zero, we suddenly find that d electrons selectively cease to propagate across the wall. s electron propagation remains, and can account for the small observed conductance steps., Comment: 9 pages, 4 figures, Surface Science, to appear

Strain can affect the morphology of a crystal surface, and cause modifications of its reconstruction even when weak, as in the case of mechanical bending. We carried out calculations of strain-dependent surface free energy and direct bending simulations demonstrating the change of incommensurate reconstruction in Au(111) under strain, in good agreement with recent data. Time-dependent strain should cause a sliding of the topmost layer over the second, suggesting an interesting case of nanofriction. Bending strain could also be used to fine tune the spacing of selectively absorbed nanoclusters., Comment: 4 pages, 5 figures. Published on PRB Rapid Comm. Major revision

An unusual quasi-two-dimensional heavy band mass van Hove singularity (vHs) lies very near the Fermi energy in MgCNi3, recently reported to superconduct at 8.5 K. This compound is strongly exchange enhanced and is unstable to ferromagnetism upon hole doping with 12% Mg --> Na or Li. The 1/4-depleted fcc (frustrated) Ni sublattice and lack of Fermi surface nesting argues against competing antiferromagnetic and charge density wave instabilities. We identify an essentially infinite mass along the M-Gamma line, leading to quasi-two-dimensionality of this vHs may promote unconventional p-wave pairing that could coexist with superconductivity., Comment: 4 two-column pages, 4 figures

We study the influence of surface roughness on the adhesion of elastic solids. Most real surfaces have roughness on many different length scales, and this fact is taken into account in our analysis. We consider in detail the case when the surface roughness can be described as a self affine fractal, and show that when the fractal dimension D_f >2.5, the adhesion force may vanish, or be at least strongly reduced. We consider the block-substrate pull-off force as a function of roughness, and find a partial detachment transition preceding a full detachment one. The theory is in good qualitative agreement with experimental data., Comment: 15 pages, 16 figures. Submitted to J. Chem. Phys

We argue that a normal Fermi liquid and a singlet, spin gapped Mott insulator cannot be continuously connected, and that some intermediate phase must intrude between them. By explicitly working out a case study where the singlet insulator is stabilized by orbital degeneracy and an inverted Hund's rule coupling, mimicking a Jahn-Teller effect, we find that the intermediate phase is a superconductor., Comment: 4 pages, 3 figure, Changes to the figures. Accepted for publication on Physical Review Letters

We propose a two-stage mechanism for the rotational surface disordering phase transition of a molecular crystal, as realized in C$_{60}$ fullerite. Our study, based on Monte Carlo simulations, uncovers the existence of a new intermediate regime, between a low temperature ordered $(2 \times 2)$ state, and a high temperature $(1 \times 1)$ disordered phase. In the intermediate regime there is partial disorder, strongest for a subset of particularly frustrated surface molecules. These concepts and calculations provide a coherent understanding of experimental observations, with possible extension to other molecular crystal surfaces., Comment: 4 pages, 2 figures

We present a calculation of the change of free energy of a solid surface upon bending of the solid. It is based on extracting the surface stress through a molecular dynamics simulation of a bent slab by using a generalized stress theorem formula, and subsequent integration of the stress with respect to strain as a function of bending curvature. The method is exemplified by obtaining and comparing free energy changes with curvature of various reconstructed Au(001) surfaces., Comment: 14 pages, 2 figures, accepted for publication in Surface Science (ECOSS-19)

Due to particle conservation, Canonical Molecular Dynamics (MD) simulations fail in the description of surface phase transitions involving coverage or lateral density changes. However, a step on the surface can act effectively as a source or a sink of atoms, in the simulation as well as in real life. A single surface step can be introduced by suitably modifying planar Periodic Boundary Conditions (PBC), to accommodate the generally inequivalent stacking of two adjacent layers. We discuss here how, through the introduction of two orthogonal surface steps, particle number conservation may no longer represent a fatal constraint for the study of these surface transitions. As an example, we apply the method for estimating temperature-induced lateral density increase of the reconstructed Au (001) surface; the resulting anisotropic cell change is consistent with experimental observations. Moreover, we implement this kind of scheme in conjunction with the variable curvature MD method, recently introduced by our group., Comment: 9 pages, 5 figures, accepted for publication in Surface Science (ECOSS-19)

We discuss the spectroscopic fingerprints that a surface Mott-Hubbard insulator should show at the intra-atomic level. The test case considered is that of the Si-terminated SiC(0001) sqrt{3}xsqrt{3} surface, which is known experimentally to be insulating. We argue that, due to the Mott-Hubbard phenomenon, spin unpaired electrons in the Si adatom dangling bonds are expected to give rise to a Si-2p core level spectrum with a characteristic three-peaked structure, as seen experimentally. This structure results from the joint effect of intra-atomic exchange, spatial anisotropy, and spin-orbit coupling. Auger intensities are also discussed., Comment: 4 pages, 2 figures, ECOSS-18 conference

We have carried out first principles electronic structure and total energy calculations for a series of ultrathin aluminum nanowires, based on structures obtained by relaxing the model wires of Gulseren et al. The number of conducting channels is followed as the wires radius is increased. The results suggest that pentagonal wires should be detectable, as the only ones who can yield a channel number between 8 and 10., Comment: 9 pages + 3 figures, to appear on Surface Science

Off-lattice Grand Canonical Monte Carlo simulations of the clean diamond (111) surface, based on the effective many-body Brenner potential, yield the $(2\times1)$ Pandey reconstruction in agreement with \emph{ab-initio} calculations and predict the existence of new meta-stable states, very near in energy, with all surface atoms in three-fold graphite-like bonding. We believe that the long-standing debate on the structural and electronic properties of this surface could be solved by considering this type of carbon-specific configurations., Comment: 4 pages + 4 figures, Phys. Rev. B Rapid Comm., in press (15Apr00). For many additional details (animations, xyz files) see electronic supplement to this paper at http://www.sci.kun.nl/tvs/carbon/meta.html

The large surface density changes associated with the (100) noble metals surface hex-reconstruction suggest the use of non-particle conserving simulation methods. We present an example of a surface Grand Canonical Monte Carlo applied to the transformation of a square non reconstructed surface to the hexagonally covered low temperature stable Au(100). On the other hand, classical Molecular Dynamics allows to investigate microscopic details of the reconstruction dynamics, and we show, as an example, retraction of a step and its interplay with the surface reconstruction/deconstruction mechanism., Comment: 9 pages, 5 figures, accepted for publication in Surf. Rev. and Letters (ICSOS-6)

Deposition/removal of metal atoms on the hex reconstructed (100) surface of Au, Pt and Ir should present intriguing aspects, since a new island implies hex -> square deconstruction of the substrate, and a new crater the square -> hex reconstruction of the uncovered layer. To obtain a microscopic understanding of how islands/craters form in these conditions, we have conducted simulations of island and crater growth on Au(100), whose atomistic behavior, including the hex reconstruction on top of the square substrate, is well described by mean s of classical many-body forces. By increasing/decreasing the Au coverage on Au(100), we find that island/craters will not grow unless they exceed a critical size of about 8-10 atoms. This value is close to that which explains the nonlinear coverage dependence observed in molecular adsorption on the closely related surface Pt (100). This threshold size is rationalized in terms of a transverse step correlation length, measuring the spatial extent where reconstruction of a given plane is disturbed by the nearby step., Comment: 11 pages, 5 figures, accepted for publication in Surface Science (ECOSS-18)

The reentrant layering transitions found in rare gas adsorption on solid substrates have conflictually been explained either in terms of preroughening (PR), or of top layer melting-solidification phenomena. We obtain adsorption isotherms of Lennard-Jones particles on an attractive substrate by off lattice Grand Canonical Monte Carlo (GCMC) simulation, and reproduce reentrant layering. Microscopic analysis, including layer-by-layer occupancies, surface diffusion and pair correlations, confirms the switch of the top surface layer from solid to quasi-liquid across the transition temperature. At the same time, layer occupancy is found at each jump to switch from close to full to close to half, indicating a disordered flat (DOF) surface and establishing preroughening as the underlying mechanism. Our results suggest that top layer melting is essential in triggering preroughening, which thus represents the threshold transition to surface melting in rare gas solids., Comment: 4 pages, 4 figures; accepted for publication in Phys. Rev. Letters

We show that two distinct $3 \times 3$ ground states, one nonmagnetic, metallic, and distorted, the other magnetic, semimetallic (or insulating) and undistorted, compete in $\alpha$-phase adsorbates on semiconductor (111) surfaces. In Sn/Ge(111), LSDA/GGA calculations indicate, in agreement with experiment, that the distorted metallic ground state prevails. The reason for stability of this state is analysed, and is traced to a sort of bond density wave, specifically a modulation of the antibonding state filling between the adatom and a Ge-Ge bond directly underneath.

A glass surface may still flow below the bulk glass transition temperature, where the underlying bulk is frozen. Assuming the existence at T=T* of a bulk thermodynamical glass transition, we show that the glass-vapor interface is generally wetted by a liquid layer of thickness ~ -ln(T*-T) when T--> T*. Contrary to standard surface melting of crystals however, the integrated value of the diffusivity across the interface remains finite for T-->T*. Difference in shape induced by bulk and by surface flow is discussed as a possible means of experimental detection of surface defreezing., Comment: five pages, three figures

We discuss the possibility that, besides roughening, a vicinal surface could display preroughening (PR), and consider the possible mechanisms for its promotion. Within the framework of a terrace-step-kink model, it turns out that a PR transition is possible, and could be induced by a short-range repulsion between parallel kinks along the same step or on adjacent steps, or even by some kind of extended range step-step repulsion. We discuss the possible relevance of this phenomenon to the anomalous roughening behaviour recently reported for Ag(115)., Comment: 9 pages, 3 postscript figures, submitted to Surface Science

We study the quantum phase transition between a band (``ionic'') insulator and a Mott-Hubbard insulator, realized at a critical value U=Uc in a bipartite Hubbard model with two inequivalent sites, whose on-site energies differ by an offset Delta. The study is carried out both in D=1 and D=2 (square and honeycomb lattices), using exact Lanczos diagonalization, finite-size scaling, and Berry's phase calculations of the polarization. The Born effective charge jump from positive infinity to negative infinity previously discovered in D=1 by Resta and Sorella is confirmed to be directly connected with the transition from the band insulator to the Mott insulating state, in agreement with recent work of Ortiz et al.. In addition, symmetry is analysed, and the transition is found to be associated with a reversal of inversion symmetry in the ground state, of magnetic origin. We also study the D=1 excitation spectrum by Lanczos diagonalization and finite-size scaling. Not only the spin gap closes at the transition, consistent with the magnetic nature of the Mott state, but also the charge gap closes, so that the intermediate state between the two insulators appears to be metallic. This finding, rationalized within unrestricted Hartree-Fock as due to a sign change of the effective on-site energy offset Delta for the minority spin electrons, underlines the profound difference between the two insulators. The band-to-Mott insulator transition is also studied and found in the same model in D=2. There too we find an associated, although weaker, polarization anomaly, with some differences between square and honeycomb lattices. The honeycomb lattice, which does not possess an inversion symmetry, is used to demonstrate the possibility of an inverted piezoelectric effect in this kind of ionic Mott insulator., Comment: 27 pages, RevTex + 8 postscript figures

We study the unconventional insulating state in A_4C_{60} with a variety of approaches, including density functional calculations and dynamical mean-field theory. While the former predicts a metallic state, in disagreement with experiment, the latter yields a (paramagnetic) Mott-Jahn-Teller insulator. In that state, conduction between molecules is blocked by on-site Coulomb repulsion, magnetism is suppressed by intra-molecular Jahn-Teller effect, and important excitations (such as optical and spin gap) should be essentially intra-molecular. Experimental gaps of 0.5 eV and 0.1 eV respectively compare well with molecular ion values, in agreement with this picture., Comment: 4 pages, 2 postscript figures

We present ab-initio electronic structure calculations for the Si-terminated SiC(0001)$\sqrt{3}\times\sqrt{3}$ surface. While local density approximation (LDA) calculations predict a metallic ground state with a half-filled narrow band, Coulomb effects, included by the spin-polarized LDA+U method, result in a magnetic (Mott-Hubbard) insulator with a gap of 1.5 eV, comparable with the experimental value of 2.0 eV. The calculated value of the inter-site exchange parameter, J=30K, leads to the prediction of a paramagnetic Mott state, except at very low temperatures. The observed Si 2p surface core level doublet can naturally be explained as an on-site exchange splitting., Comment: RevTex, 4 pages, 4 eps-figures

We consider an exchange model describing two isotropic spin-1/2 Heisenberg antiferromagnets coupled by a quartic term on the square lattice. The model is relevant for systems with orbital degeneracy and strong electron-vibron coupling in the large Hubbard repulsion limit, and is known to show a spin-Peierls-like dimerization in one dimension. In two dimensions we calculate energy gaps, susceptibilities, and correlation functions with a Green's Function Monte Carlo. We find a finite spin gap and no evidence of any kind of order. We conclude that the ground state is, most likely, a spin liquid of resonating valence bonds., Comment: 4 pages, 4 figures, Revtex

We have examined theoretically the spontaneous thinning process of tip-suspended nanowires, and subsequently studied the structure and stability of the monatomic gold wires recently observed by Transmission Electron Microscopy (TEM). The methods used include thermodynamics, classical many-body force simulations, Local Density (LDA) and Generalized Gradient (GGA) electronic structure calculations as well as ab-initio simulations including the two tips. The wire thinning is well explained in terms of a thermodynamic tip suction driving migration of surface atoms from the wire to the tips. For the same reason the monatomic wire becomes progressively stretched. Surprisingly, however, all calculations so far indicate that the stretched monatomic gold wire should be unstable against breaking, contrary to the apparent experimental stability. The possible reasons for the observed stability are discussed., Comment: 4 figures

We consider the exchange Hamiltonian H_ST = -J Sum_{} (2 S_r S_r' - 1/2) (2 T_r T_r' - 1/2) describing two isotropic spin-1/2 Heisenberg antiferromagnets coupled by a quartic term on equivalent bonds. The model is relevant for systems with orbital degeneracy and strong electron-vibron coupling in the large Hubbard repulsion limit. To investigate the ground state properties we use a Green's Function Monte Carlo, calculating energy gaps and correlation functions, the latter through the forward walking technique. In one dimension we find that the ground state is a ``crystal'' of valence bond dimers. In two dimensions, the spin gap appears to remain finite in the thermodynamic limit, and, consistently, the staggered magnetization -- signal of Neel long range order -- seems to vanish. From the analysis of dimer-dimer correlation functions, however, we find no sign of a valence bond crystal. A spin liquid appears as a plausible scenario compatible with our findings., Comment: 3 pages, 3 figures

Suppose you look at today's stock prices and bet on the value of the first digit. One could guess that a fair bet should correspond to the frequency of $1/9 = 11.11%$ for each digit from 1 to 9. This is by no means the case, and one can easily observe a strong prevalence of the small values over the large ones. The first three integers 1,2 and 3 alone have globally a frequency of 60% while the other six values 4, 5, 6, 7, 8 and 9 appear only in 40% of the cases. This situation is actually much more general than the stock market and it occurs in a variety of number catalogs related to natural phenomena. The first observation of this property traces back to S. Newcomb in 1881 but a more precise account was given by F. Benford in 1938. In this note we illustrate these observations with the enlightening specific example of the stock market. We also identify the general mechanism for the origin of this uneven distribution in the multiplicative nature of fluctuations in economics and in many natural phenomena. This provides a natural explanation for the ubiquitous presence of the Benford's law in many different phenomena with the common element that their fluctuations refer to a fraction of their values. This brings us close to the problem of the spontaneous origin of scale invariant properties in various phenomena which is a debated question at the frontier of different fields., Comment: 9 pages, 2 figures; improved version with added references

In a recent paper [Phys. Rev. Lett. 80 (1998) 1690], Krim and coworkers have observed that the friction force, acting on a thin physisorbed layer of N_2 sliding on a lead film, abruptly decreases by a factor of ~2 when the lead film is cooled below its superconductivity transition temperature. We discuss the possible mechanisms for the abruptness of the sliding friction drop, and also discuss the relevance of these results to the problem of electronic friction., Comment: 5 pages, no figures

We present ab-initio electronic structure results on the surface of sqr3 adsorbates. In particular, we address the issue of metal-insulator instabilities, charge-density-waves (CDWs) or spin-density-waves (SDWs), driven by partly filled surface states and their 2D Fermi surface, and/or by the onset of magnetic instabilities. The focus is both on the newly discovered commensurate CDW transitions in the Pb/Ge(111) and Sn/Ge(111) structures, and on the puzzling semiconducting behavior of the Pb/Ge(111), K/Si(111):B and SiC(0001) surfaces. In all cases, the main factor driving the instability appears to be an extremely narrow surface state band. We have carried out so far preliminary calculations for the Si/Si(111) surface, chosen as our model system, within the gradient corrected local density (LDA+GC) and local spin density (LSD+GC) approximations, with the aim of understanding the possible interplay between 2D Fermi surface and electron correlations in the surface + adsorbate system. Our spin- unrestricted results show that the sqr3 paramagnetic surface is unstable towards a commensurate SDW with periodicity 3x3 and magnetization 1/3., Comment: 9 pages, 4 Postscript figures, to be published in Surf. Sci

Motivated by the recent experimental evidence of commensurate surface CDW in Pb/Ge(111) and Sn/Ge(111) $\sqrt{3}$-adlayer structures, as well as by the insulating states found on K/Si(111):B and SiC(0001), we have investigated the role of electron-electron interactions, and also of electron-phonon coupling, on the narrow surface state band originating from the dangling bond orbitals of the adsorbate. We model the problem by an extended two-dimensional Hubbard model at half-filling on a triangular lattice. We include an on-site Hubbard repulsion U and a nearest-neighbor V, plus a long-ranged Coulomb tail. The electron-phonon interaction is treated in the deformation potential approximation. We have explored the phase diagram of the model including the possibility of commensurate 3x3 phases, using mainly the Hartree-Fock approximation. For U larger than the bandwidth we find magnetic insulators, possibly corresponding to the situation in SiC and in K/Si. For smaller U, the inter-site repulsion V can stabilize metallic CDW phases, reminiscent of the 3x3 structures of Sn/Ge, and possibly of Pb/Ge., Comment: 10 pages, 3 figures, presented at ECOSS-17 sept 1997

In the context of perovskite quantum paraelectrics, we study the effects of a quadrupolar interaction $J_q$, in addition to the standard dipolar one $J_d$. We concentrate here on the nonlinear dielectric response $\chi_{P}^{(3)}$, as the main response function sensitive to quadrupolar (in our case antiquadrupolar) interactions. We employ a 3D quantum four-state lattice model and mean-field theory. The results show that inclusion of quadrupolar coupling of moderate strength ($J_q \sim {{1}\over{4}} J_d$) is clearly accompanied by a double change of sign of $\chi_{P}^{(3)}$ from negative to positive, near the quantum temperature $T_Q$ where the quantum paraelectric behaviour sets in. We fit our $\chi_{P}^{(3)}$ to recent experimental data for SrTiO$_3$, where the sign change is identified close to $T_Q \sim 37 K$., Comment: 22 pages

We discuss the results of a study of restricted solid-on-solid models for fcc (110) surfaces. These models are simple modifications of the exactly solvable BCSOS model, and are able to describe a $(2\times 1)$ missing-row reconstructed surface as well as an unreconstructed surface. They are studied in two different ways. The first is by mapping the problem onto a quantum spin-1/2 one-dimensional hamiltonian of the Heisenberg type, with competing $S^z_iS^z_j$ couplings. The second is by standard Monte Carlo simulations. We find phase diagrams with the following features, which we believe to be quite generic: (i) two flat, ordered phases (unreconstructed and missing-row reconstructed); a rough, disordered phase; an intermediate disordered flat (DF) phase, characterized by monoatomic steps, whose physics is shown to be akin to that of a dimer spin state. (ii) a transition line from the $(2\times 1)$ reconstructed phase to the DF phase showing exponents which appear to be close, within our numerical accuracy, to the 2D-Ising universality class. (iii) a critical (preroughening) line with variable exponents, separating the unreconstructed phase from the DF phase. Possible signatures and order parameters of the DF phase are investigated., Comment: Revtex (22 pages) + 15 figures (uuencoded file)

We have carried out an ab-initio study of $\alpha$-Sn (111), with the aim of predicting and understanding its structure, reconstructions, and electronic states. We consider a variety of structural possibilities, and optimize them by moving atoms according to Hellmann-Feynman forces. Our results indicate that the unreconstructed surface is highly unstable, while a variety of reconstructions compete for the true ground state. Extrapolated trends from diamond to Si to Ge are well borne out, with a $2\times 1 \pi$-bonded chain reconstruction prevailing in the absence of adatoms, and a $c(4\times 2)$ or $(2\times 2)$ basic adatom-restatom unit reconstruction otherwise. Accompanying surface bucklings are in both cases larger than in Si and Ge, with consequently large ionic charge transfers predicted. Search for a $\beta$-Sn-like metallic state of the surface turned out to be inconclusive., Comment: 27 pages,12 Postscript figures,submitted to Phy.Rev.B

We study the Jahn-Teller effect of positive fullerene ions $^2$C$_{60}^{+}$ and $^1$C$_{60}^{2+}$. The aim is to discover if this case, in analogy with the negative ion, possesses a Berry phase or not, and what are the consequences on dynamical Jahn-Teller quantization. Working in the linear and spherical approximation, we find no Berry phase in $^1$C$_{60}^{2+}$, and presence/absence of Berry phase for coupling of one $L=2$ hole to an $L=4$/$L=2$ vibration. We study in particular the special equal-coupling case ($g_2=g_4$), which is reduced to the motion of a particle on a 5-dimensional sphere. In the icosahedral molecule, the final outcome assesses the presence/absence of a Berry phase of $\pi$ for the $h_u$ hole coupled to $G_g$/$H_h$ vibrations. Some qualitative consequences on ground-state symmetry, low-lying excitations, and electron emission from C$_{60}$ are spelled out., Comment: 31 pages (RevTeX), 3 Postscript figures (uuencoded)

I present an informal discussion of various cases where two-dimensional surface metal-insulator structural and charge-density-wave instabilities driven by partly filled surface states have been advocated. These include reconstructions of clean semiconductor surfaces and of W(100) and Mo(100), as well as anomalies on the hydrogen-covered surfaces H/W(110) and H/Mo(110), and possibly alkali-covered surfaces such as K/Cu(111). In addition I will also discuss the opposite type of phenomena, namely surface insulator-metal transitions, which can be argued to occur on $\alpha -$% Ga(001), high-temperature Ge(111), and probably Be(0001)., Comment: 33 pages, revtex, figure available via fax+39-40-3787-528

Recent work has raised considerable interest on the nature of thin metallic wires. We have investigated the melting behavior of thin cylindrical Pb wires with the axis along a (110) direction, using molecular dynamics and a well-tested many-body potential. We find that---in analogy with cluster melting---the melting temperature $T_m (R)$ of a wire with radius $R$ is lower than that of a bulk solid, $T_m^b$, by $T_m (R) = T_m^b -c/R$. Surface melting effects, with formation of a thin skin of highly diffusive atoms at the wire surface, is observed. The diffusivity is lower where the wire surface has a flat, local (111) orientation, and higher at (110) and (100) rounded areas. The possible relevance to recent results on non-rupturing thin necks between an STM tip and a warm surface is addressed., Comment: 10 pages, 4 postscript figures are appended, RevTeX, SISSA Ref. 131/94/CM/SS

We show that electron hopping in a lattice of molecules possessing a Berry phase naturally leads to pairing. Our building block is a simple molecular site model inspired by C$_{60}$, but realized in closer similarity with Na$_3$. In the resulting model electron hopping must be accompanied by orbital operators, whose function is to switch on and off the Berry phase as the electron number changes. The effective hamiltonians (electron-rotor and electron-pseudospin) obtained in this way are then shown to exhibit a strong pairing phenomenon, by means of 1D linear chain case studies. This emerges naturally from numerical studies of small $N$-site rings, as well as from a BCS-like mean-field theory formulation. The pairing may be explained as resulting from the exchange of singlet pairs of orbital excitations, and is intimately connected with the extra degeneracy implied by the Berry phase when the electron number is odd. The relevance of this model to fullerides, to other molecular superconductors, as well as to present and future experiments, is discussed., Comment: 30 pages, RevTex

Available simulation methods, suitable to describe solid-solid phase transitions occurring upon increasing of presssure and/or temperature, are based on empirical interatomic potentials: this restriction reduces the predictive power, and thus the general usefulness of numeric simulations in this very relevant field. We present a new simulation scheme which allows, for the first time, the simulation of these phenomena with the correct quantum-mechanical description of interatomic forces and internal stress, along with the correct statistical mechanics of ionic degrees of freedom. The method is obtained by efficiently combining the Car-Parrinello method for ab- initio molecular dynamics with the Parrinello Rahman method to account for a variable cell shape. Within this scheme phase trasformations may spontaneously take place during the simulation with variation of external pressure and/or temperature. The validity of the method is demonstrated by simulating the metal-insulator transition in Silicon (from diamond structure to simple hexagonal structure) under high pressure., Comment: 11 pages REVTEX, 4 figures available upon request, SISSA-CM-93/XXX

We have considered the processes which lead to elastic scattering between two far ultraviolet or X--ray photons while they propagate inside a solid, modeled as a simple electron gas. The new ingredient, with respect to the standard theory of photon--photon scattering in vacuum, is the presence of low--energy, nonrelativistic electron--hole excitations. Owing to the existence of two--photon vertices, the scattering processes in the metal are predominantly of second order, as opposed to fourth order for the vacuum case. The main processes in second order are dominated by exchange of virtual plasmons between the two photons. For two photons of similar energy $\hbar \Omega$, this gives rise to a cross section rising like $\Omega^2$ up to maximum of around $10^{-32}$~cm$^2$, and then decreasing like $\Omega^{-6}$. The maximal cross section is found for the photon wavevector $k \sim k_{F}$, the Fermi surface size, which typically means a photon energy $\hbar \Omega$ in the keV range. Possible experiments aimed at checking the existence of these rare but seemingly measurable elastic photon--photon scattering processes are discussed, using in particular intense synchrotron sources., Comment: 33 pages, TeX, Version 3.1, S.I.S.S.A. preprint 35/93/CM

We have begun a study of quantum ferroelectrics and paraelectrics. Simple 2D short-range lattice model hamiltonians are constructed, keeping in mind the phenomenology of real perovskite systems, like $SrTiO_{3}$ and $KTaO_{3}$. Pertinent quantum tunneling phenomena, and the presence of an ice-like constraint are demonstrated. The two simplest models, namely a plain quantum four-state clock model, and a constrained one, are then studied in some detail. We show the equivalence of the former, but not of the latter, to a quantum Ising model. For the latter, we describe a very good analytical wavefunction valid in the special case of zero coupling ($J = 0$). In order to study the full quantum statistical mechanics of both models, a Path Integral Monte Carlo calculation is set up, and implemented with a technique, which even in the constrained case permits a good convergence for increasing time slice number $m$. The method is applied first to the unconstrained model, which serves as a check, and successively to the constrained quantum four-state clock model. It is found that in both cases, a quantum phase transition still takes place at finite coupling J, between a ferroelectric and a quantum paraelectric state, even when the constraint hinders disordering of the ferroelectric state. This model paraelectric state has a finite excitation gap, and no broken symmetry. The possible role of additional ("oxygen hopping") kinetic terms in making closer contact with the known phenomenology of $SrTiO_{3}$ is discussed., Comment: 47 pages, TeX, Version 3.1, S.I.S.S.A. preprint 109/93/CM/MB

We propose, based on ab initio calculations, that the ground state (001) surface of $\alpha$-Ga should be covered by two layers of GaIII, a denser phase which is stable in the bulk only at high pressure and temperature. We discuss how this novel prediction can explain recent STM observations, including anomalous thermal stability of this surface., Comment: 10 pages, figures available by fax upon request, REVTEX, ref. SISSA 30/93/CM/SS/SC