Abstract: Two CoSi2/Si1-xGex/Si(100) heterostructures, with different Ge content, made by molecular beam epitaxy are characterized by high resolution electron microscopy. In general, the interface between the CoSi2 thin film and the Si1-xGex layer is of a high structural quality and the strained Si1-xGex layer exhibits few defects. For both samples, different interface structures are present, although the dominant interfacial configuration is similar to the unreconstructed interface present at the CoSi2/Si(100) interface. Only occasionally (2x1) reconstructed interface regions are found which are just a few nanometers in length. Phenomena such as Ge segregation and the introduction of defects are also observed in the Si1-xGex layer. We attribute the minimal presence of the reconstructed interface to both the (2x8):Si1-xGex(100) surface reconstruction and the Ge segregation that takes place. (C) 1999 American Institute of Physics. [S0021-8979(99)02104-0].

Abstract: We utilize an exact variational numerical procedure to calculate the ground-state properties of a polaron in the presence of Rashba and linear Dresselhaus spin-orbit coupling. We find that when the linear Dresselhaus spin-orbit coupling approaches the Rashba spin-orbit coupling, the Van Hove singularity in the density of states will be shifted away from the bottom of the band and finally disappear when the two spin-orbit couplings are tuned to be equal. The effective mass will be suppressed; the trend will become more significant for low phonon frequency. The presence of two dominant spin-orbit couplings will make it possible to tune the effective mass with more varied observables.

Abstract: Magnetotransport in bilayer graphene in the presence of a weak and periodic potential is investigated in the presence of a perpendicular magnetic field B. The modulation broadens the Landau levels into bands and for weak magnetic fields leads to the well-known Weiss oscillations in their bandwidth and their transport coefficients at very low B and to the Shubnikov-de Haas oscillations at larger B. The amplitude of the Weiss oscillations is severely reduced if the periodic potentials applied to the two layers oscillate out of phase. We also contrast some results with those corresponding to single-layer graphene. Relative to them the flat-band condition and the oscillation amplitude differ substantially, due to the interlayer coupling, and agree only when this coupling is extremely weak. We further show that the Hall conductivity exhibits the well-known steps at half-integer and integer multiples of 4e(2)/h in single-layer and bilayer graphene, respectively, even for very weak magnetic fields. The results are pertinent to weak and periodic corrugations when the potential modulation dominates the strain-induced magnetic modulation.

Abstract: We performed a detailed investigation of the oxidation of the technologically relevant GaAs(001)-beta 2(2x4) surface via density functional calculations. The purpose is to gain insights on the atomistic mechanisms and local bondings that underlie the degradation of the surface properties once exposed to oxygen. The study comprises the adsorption of single O atoms, through the sampling of several adsorption sites, and the subsequent formation of the O adsorbate at increasing coverage by taking into account multiple-atom adsorption. Based on the evaluation of the energetics and the structural properties of the atomistic models generated, the results here reported delineate a consistent picture of the initial stage of the surface oxidation: (i) at low coverage, in the limit of single O insertions, oxygen is incorporated on the surface forming a twofold-bridging Ga-O-As bond; (ii) at increasing coverage, as multiple O atoms are involved, this is accompanied by the formation of a threefold-coordinated bond (with two Ga and one As atoms); (iii) the latter has important implications regarding the electronic properties of the adsorbate since this O bonding may result in the formation of As dangling bonds. Moreover, a clear trend of increased energy gain for the incorporation of neighboring O atoms compared to single O insertions indicates that the formation of oxide clusters is favored over a regime of uniform oxidation. Our findings provide a detailed description of the O bonding and stress the importance of modeling the adsorption of multiple O atoms for an accurate description of the surface oxidation.

Abstract: The Coulomb interaction and the correlation of a remote electron with a single layer of graphene is investigated in the presence of a magnetic field applied perpendicular to the graphene layer. The remote electron polarizes the electron gas in the graphene layer, which we describe in terms of excitations of virtual plasmons in graphene. The composite quasiparticle formed by electron plus polarization is called a plasmon polaron. The ground-state energy of this quasiparticle is calculated within perturbation theory for remote electrons in different environments.

Abstract: We present a comprehensive study of the crystal structure, magnetic structure, and microscopic magnetic model of (CuBr)LaNb2O7, the Br analog of the spin-gap quantum magnet (CuCl) LaNb2O7. Despite similar crystal structures and spin lattices, the magnetic behavior and even peculiarities of the atomic arrangement in the Cl and Br compounds are very different. The high- resolution x-ray and neutron data reveal a split position of Br atoms in (CuBr) LaNb2O7. This splitting originates from two possible configurations developed by [CuBr] zigzag ribbons. While the Br atoms are locally ordered in the ab plane, their arrangement along the c direction remains partially disordered. The predominant and energetically more favorable configuration features an additional doubling of the c lattice parameter that was not observed in (CuCl) LaNb2O7. (CuBr) LaNb2O7 undergoes long-range antiferromagnetic ordering at T-N = 32 K, which is nearly 70% of the leading exchange coupling J4 similar or equal to 48 K. The Br compound does not show any experimental signatures of low-dimensional magnetism because the underlying spin lattice is three-dimensional. The coupling along the c direction is comparable to the couplings in the ab plane, even though the shortest Cu-Cu distance along c (11.69 angstrom) is three times larger than nearest-neighbor distances in the ab plane (3.55 angstrom). The stripe antiferromagnetic long-range order featuring columns of parallel spins in the ab plane and antiparallel spins along c is verified experimentally and confirmed by the microscopic analysis.

Abstract: Iron carbides play a crucial role in steel manufacturing and processing and to a large extent determine the physical properties of steel products. The modified embedded atom method (MEAM) in combination with Lee's Fe-C potential is a good candidate for molecular dynamics simulations on larger Fe-C systems. Here, we investigate the stability and crystal structures of pure iron and binary iron carbides using MEAM and compare them with the experimental data and quantum-mechanical density functional theory calculations. The analysis shows that the Fe-C potential gives reasonable results for the relative stability of iron and iron carbides. The performance of MEAM for the prediction of the potential energy and the calculated lattice parameters at elevated temperature for pure iron phases and cementite are investigated as well. The conclusion is that Lee's MEAM Fe-C potential provides a promising basis for further molecular dynamics simulations of Fe-C alloys and steels at lower temperatures (up to 800 K).

Abstract: Using atomistic simulations we investigate the morphological properties of graphene deposited on top of a nanostructured substrate. Sinusoidally corrugated surfaces, steps, elongated trenches, one-dimensional and cubic barriers, spherical bubbles, Gaussian bumps, and Gaussian depressions are considered as support structures for graphene. The graphene-substrate interaction is governed by van der Waals forces and the profile of the graphene layer is determined by minimizing the energy using molecular dynamics simulations. Based on the obtained optimum configurations, we found that (i) for graphene placed over sinusoidally corrugated substrates with corrugation wavelengths longer than 2 nm, the graphene sheet follows the substrate pattern while for supported graphene it is always suspended across the peaks of the substrate, (ii) the conformation of graphene to the substrate topography is enhanced when increasing the energy parameter in the van der Waals model, (iii) the adhesion of graphene into the trenches depends on the width of the trench and on the graphene's orientation, i. e., in contrast to a small-width (3 nm) nanoribbon with armchair edges, the one with zigzag edges follows the substrate profile, (iv) atomic-scale graphene follows a Gaussian bump substrate but not the substrate with a Gaussian depression, and (v) the adhesion energy due to van der Waals interaction varies in the range [0.1-0.4] J/m(2).

Abstract: We present a tight-binding investigation of strained bilayer graphene within linear elasticity theory, focusing on the different environments experienced by the A and B carbon atoms of the different sublattices. We find that the inequivalence of the A and B atoms is enhanced by the application of perpendicular strain epsilon(zz), which provides a physical mechanism for opening a band gap, most effectively obtained when pulling the two graphene layers apart. In addition, perpendicular strain introduces electron-hole asymmetry and can result in linear electronic dispersion near the K point. Our findings suggest experimental means for strain-engineered band gaps in bilayer graphene.

Abstract: To date, several experimental groups reported measurements of the thickness dependence of T-c of atomically uniform single-crystalline Pb nanofilms. The reported amplitude of the T-c oscillations varies significantly from one experiment to another. Here we propose that the reason for this unresolved issue is an interplay of the quantum-size variations in the single-electron density of states with thickness-dependent oscillations in the phonon-mediated electron-electron coupling. Such oscillations in the coupling depend on the substrate material, the quality of the interface, the protection cover, and other details of the fabrication process, changing from one experiment to another. This explains why the available data do not exhibit one-voice consistency about the amplitude of the T-c oscillations. Our analyses are based on a numerical solution of the Bogoliubov-de Gennes equations for a superconducting slab.

Abstract: The lattice dynamics of low-frequency rigid-plane modes in metallic (graphene multilayers, GML) and in insulating (hexagonal boron-nitride multilayers, BNML) layered crystals is investigated. The frequencies of shearing and compression (stretching) modes depend on the layer number N and are presented in the form of fan diagrams. The results for GML and BNML are very similar. In both cases, only the interactions (van der Waals and Coulomb) between nearest-neighbor planes are effective, while the interactions between more distant planes are screened. A comparison with recent Raman scattering results on low-frequency shear modes in GML [Tan et al., Nat. Mater., in press, doi: 10.1038/nmat3245, (2012)] is made. Relations with the low-lying rigid-plane phonon dispersions in the bulk materials are established. Master curves, which connect the fan diagram frequencies for any given N, are derived. Static and dynamic thermal correlation functions for rigid-layer shear and compression modes are calculated. The results might be of use for the interpretation of friction force experiments on multilayer crystals.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

Abstract: We consider two possible mechanisms for single-vortex fluctuative entry/exit through the surface barrier in ultrathin superconducting disk-shaped nanoislands made of Pb and consisting of just a few monoatomic layers, which can be fabricated using modern techniques. We estimate tunneling probabilities and establish criteria for the crossover between these two mechanisms depending on magnetic field and system sizes. For the case of vortex entry, quantum tunneling dominates on the major part of the temperature/flux phase diagram. For the case of vortex exit, thermal activation turns out to be more probable. This nontrivial result is due to the subtle balance between the barrier height and width, which determine rates of the thermal activation and quantum tunneling, respectively.