Applications of a tight-binding total-energy method for transition and noble metals: Elastic constants, vacancies, and surfaces of monatomic metals

Transferable orthogonal tight-binding parameters for ZnS and CdS

Calculations of Slater-Koster (SK) parameters appearing in the tight-binding method using sp(3)d(5) basis sets for both the cationic and anionic species are presented for ZnS and CdS. We have adjusted these parameters to match the band structures obtained from the full potential linear augmented plane wave method. This operation has been carried out for a variety of structures namely zinc blende, wurtzite, rocksalt, CsCl and for a wide range of near-neighbor distances. The SK parameters have slightly different values for the same near-neighbor distance in different structures. Therefore, a least-squares fitting has been performed separately for each parameter as a function of only the near-neighbor distance to guarantee the transferability of these parameters to different structural environments. The fitted parameters are then used to calculate the electronic structure of small-sized clusters of ZnS and CdS in given geometries and the results are compared with ab initio results. A fairly good agreement found in the one-electron energy spectrum and total energy confirms transferability of the parameters to different length scales. A detailed account of the calculation procedure and calibration results is given in the present paper. These parameters can be used to study the electronic structure of large-sized clusters where first-principles methods are computationally demanding. It may be mentioned that the SK parameters do not satisfy the R(-(l + l' + 1)) Harrison scaling law for larger values of the near-neighbor distance R.

Structure and magnetism of bulk Fe and Cr: from plane waves to LCAO methods

Magnetic, structural and energetic properties of bulk Fe and Cr were studied using first-principles calculations within density functional theory (DFT). We aimed to identify the dependence of these properties on key approximations of DFT, namely the exchange-correlation functional, the pseudopotential and the basis set. We found a smaller effect of pseudopotentials (PPs) on Fe than on Cr. For instance, the local magnetism of Cr was shown to be particularly sensitive to the potentials representing the core electrons, i.e. projector augmented wave and Vanderbilt ultrasoft PPs predict similar results, whereas standard norm-conserving PPs tend to overestimate the local magnetic moments of Cr in bcc Cr and in dilute bcc FeCr alloys. This drawback is suggested to be closely correlated to the overestimation of Cr solution energy in the latter system. On the other hand, we point out that DFT methods with very reduced localized basis sets (LCAO: linear combination of atomic orbitals) give satisfactory results compared with more robust plane-wave approaches. A minimal-basis representation of '3d' electrons comes to be sufficient to describe non-trivial magnetic phases including spin spirals in both fcc Fe and bcc Cr, as well as the experimental magnetic ground state of bcc Cr showing a spin density wave (SDW) state. In addition, a magnetic 'spd' tight binding model within the Stoner formalism was proposed and validated for Fe and Cr. The respective Stoner parameters were obtained by fitting to DFT data. This efficient semiempirical approach was shown to be accurate enough for studying various collinear and non-collinear phases of bulk Fe and Cr. It also enabled a detailed investigation of different polarization states of SDW in bcc Cr, where the longitudinal state was suggested to be the ground state, consistent with existing experimental data.

Domain boundary formation in helical multishell gold nanowires

Helical multishell gold nanowires are studied theoretically for the formation mechanism of the helical domain boundary. Nanowires with a wire length of more than 10 nm are relaxed by quantum mechanical molecular dynamics simulation with a tight-binding form Hamiltonian. In the results, non-helical nanowires are transformed into helical ones with the formation of atom pair defects at the domain boundary, where the defective atom pair is moved from an inner shell. Analysis of local electronic structure shows a competitive feature of the energy gain of reconstruction on the wire surface and the energy loss of the defect formation. A simple energy scaling theory gives a general explanation of domain boundary formation.

Silicon carbide nanostructures: a tight binding approach

A tight-binding model Hamiltonian is newly parametrized for silicon carbide based on fits to a database of energy points calculated within the density functional theory approach of the electronic energy surfaces of nanoclusters and the total energy of bulk 3C and 2H polytypes at different densities. This TB model includes s and p angular momentum symmetries with nonorthogonal atomic basis functions. With the aid of the new TB model, minima of silicon carbide cagelike clusters, nanotubes, ring-shaped ribbons, and nanowires are predicted. Energetics, structure, growth sequences, and stability patterns are reported for the nanoclusters and nanotubes. The band structure of SiC nanotubes and nanowires indicates that the band gap of the nanotubes ranges from 0.57 to 2.38 eV depending on the chirality, demonstrating that these nanotubes are semiconductors or insulators. One type of nanowire is metallic, another type is semiconductor, and the rest are insulators.

High thermal conductivity of a hydrogenated amorphous silicon film

We measured the thermal conductivity kappa of an 80 microm thick hydrogenated amorphous silicon film prepared by hot-wire chemical-vapor deposition with the 3omega (80-300 K) and the time-domain thermo-reflectance (300 K) methods. The kappa is higher than any of the previous temperature dependent measurements and shows a strong phonon mean free path dependence. We also applied a Kubo based theory using a tight-binding method on three 1000 atom continuous random network models. The theory gives higher kappa for more ordered models, but not high enough to explain our results, even after extrapolating to lower frequencies with a Boltzmann approach. Our results show that this material is more ordered than any amorphous silicon previously studied.

Helical [110] gold nanowires make longer linear atomic chains

Quantum mechanical molecular dynamics shows that gold nanowires formed along the [110] direction reconstruct upon stress to form helical nanowires. The mechanism for this formation is discussed. These helical nanowires evolve on stretching to form linear atomic chains. Because helical nanowires do not form symmetrical tips, a requirement to stop the growth of atomic chains, these nanowires produce longer atomic chains than other nanowires. These results are obtained resorting to the use of tight-binding molecular dynamics and ab initio electronic structure calculations.

Are conductance plateaus independent events in atomic point contact measurements? A statistical approach

Conductance-elongation curves of gold atomic wires are measured using a scanning tunneling microscope break junction technique at room temperature. Landauer's conductance plateaus are individually identified and statistically analyzed. Both the probabilities to observe and the lengths of the two last plateaus (at conductance values close to 2e(2)/h and 4e(2)/h) are studied. All results converge to show that the occurrences of these two conductance plateaus on a conductance-elongation curve are statistically independent events.

Two-stage formation model and helicity of gold nanowires

A model for the formation of helical multishell gold nanowires is proposed and is confirmed with quantum mechanical molecular dynamics simulations. The model can explain the magic number of the helical gold nanowires in the multishell structure. The reconstruction from ideal nonhelical to realistic helical nanowires consists of two stages: dissociations of atoms on the outermost shell from atoms on the inner shell and slip deformations of atom rows generating (111)-like structure on the outermost shell. The elementary processes are governed by competition between energy loss and gain by s and d electrons together with the width of the d band. The possibility for the helical nanowires of platinum, silver, and copper is discussed.

Tight-binding molecular dynamics study of the role of defects on carbon nanotube moduli and failure

We performed tight-binding molecular dynamics on single-walled carbon nanotubes with and without a variety of defects to study their effect on the nanotube modulus and failure through bond rupture. For a pristine (5,5) nanotube, Young's modulus was calculated to be approximately 1.1 TPa, and brittle rupture occurred at a strain of 17% under quasistatic loading. The predicted modulus is consistent with values from experimentally derived thermal vibration and pull test measurements. The defects studied consist of moving or removing one or two carbon atoms, and correspond to a 1.4% defect density. The occurrence of a Stone-Wales defect does not significantly affect Young's modulus, but failure occurs at 15% strain. The occurrence of a pair of separated vacancy defects lowers Young's modulus by approximately 160 GPa and the critical or rupture strain to 13%. These defects apparently act independently, since one of these defects alone was independently determined to lower Young's modulus by approximately 90 GPa, also with a critical strain of 13%. When the pair of vacancy defects adjacent, however, Young's modulus is lowered by only approximately 100 GPa, but with a lower critical strain of 11%. In all cases, there is noticeable strain softening, for instance, leading to an approximately 250 GPa drop in the apparent secant modulus at 10% strain. When a chiral (10,5) nanotube with a vacancy defect was subjected to tensile strain, failure occurred through a continuous spiral-tearing mechanism that maintained a high level of stress (2.5 GPa) even as the nanotube unraveled. Since the statistical likelihood of defects occurring near each other increases with nanotube length, these studies may have important implications for interpreting the experimental distribution of moduli and critical strains.

Nonequilibrium Green's function method for thermal transport in junctions

We present a detailed treatment of the nonequilibrium Green's function method for thermal transport due to atomic vibrations in nanostructures. Some of the key equations, such as self-energy and conductance with nonlinear effect, are derived. A self-consistent mean-field theory is proposed. Computational procedures are discussed. The method is applied to a number of systems including one-dimensional chains, a benzene ring junction, and carbon nanotubes. Mean-field calculations of the Fermi-Pasta-Ulam model are compared with classical molecular dynamics simulations using a generalized Langevin heat bath. We find that nonlinearity suppresses thermal transport even at moderately high temperatures.

Gold nanowires and the effect of impurities

Metal nanowires and in particular gold nanowires have received a great deal of attention in the past few years. Experiments on gold nanowires have prompted theory and simulation to help answer questions posed by these studies. Here we present results of computer simulations for the formation, evolution and breaking of very thin Au nanowires. We also discuss the influence of contaminants, such as atoms and small molecules, and their effect on the structural and mechanical properties of these nanowires.

Insights into the fracture mechanisms and strength of amorphous and nanocomposite carbon

Tight-binding molecular dynamics simulations shed light into the fracture mechanisms and the ideal strength of tetrahedral amorphous carbon and of nanocomposite carbon containing diamond crystallites, two of the hardest materials. It is found that fracture in the nanocomposites, under tensile or shear load, occurs intergrain and so their ideal strength is similar to the pure amorphous phase. The onset of fracture takes place at weakly bonded sites in the amorphous matrix. On the other hand, the nanodiamond inclusions significantly enhance the elastic moduli, which approach those of diamond.

Ab initio based tight-binding molecular dynamics simulation of the sticking and scattering of O2∕Pt(111)

The sticking and scattering of O(2)Pt(111) has been studied by tight-binding molecular dynamics simulations based on an ab initio potential energy surface. We focus, in particular, on the sticking probability as a function of the angle of incidence and the energy and angular distributions in scattering. Our simulations provide an explanation for the seemingly paradox experimental findings that adsorption experiments suggest that the O(2)Pt(111) interaction potential should be strongly corrugated while scattering experiments indicate a rather small corrugation. The potential energy surface is indeed strongly corrugated which leads to a pronounced dependence of the sticking probability on the angle of incidence. The scattered O(2) molecules, however, experience a rather flat surface due to the fact that they are predominantly scattered at the repulsive tail of the potential.

Approximate density functionals applied to molecular quantum dots

Recently, molecular quantum dots (MQDs) have been investigated experimentally and found to exhibit the Kondo effect. The Kondo effect leads to an enhancement of the zero-voltage conductance. Here, we study a finite cluster model of a MQD by means of Kohn-Sham density functional theory. Furthermore, employing an implementation of Landauer's formula, we calculate the conductance of the dot. We find that the electronic structure and the molecular conductance depend strongly on the exchange-correlation functional employed. While the local spin density approximation and the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation qualitatively reproduce certain features of the Kondo effect, PBE hybrid does not. Based on the MQD, we discuss the limitations of using density functional theory to model molecular electronic devices.

Gold/palladium and silver/palladium colloids as novel metallic substrates for surface-enhanced Raman scattering

New surface-enhanced Raman scattering (SERS) substrates, composed of gold or silver colloidal nanoparticles doped with palladium, were prepared. These novel colloids are stable and maintain a satisfactory SERS efficiency, even after long aging. The interest in doping the coinage metal nanoparticles with palladium is due to the well-known catalytic activity of this metal. Transmission electron microscopy (TEM) and ultraviolet-visible absorption spectroscopy were used to characterize the shape and size of the metal particles. It was found that these bimetallic colloidal nanoparticles have a core-shell structure, with gold or silver coated with palladium clusters.

Using molecular similarity to construct accurate semiempirical electronic structure theories

Ab initio electronic structure methods give accurate results for small systems, but do not scale well to large systems. Chemical insight tells us that molecular functional groups will behave approximately the same way in all molecules, large or small. This molecular similarity is exploited in semiempirical methods, which couple simple electronic structure theories with parameters for the transferable characteristics of functional groups. We propose that high-level calculations on small molecules provide a rich source of parametrization data. In principle, we can select a functional group, generate a large amount of ab initio data on the group in various small-molecule environments, and "mine" this data to build a sophisticated model for the group's behavior in large environments. This work details such a model for electron correlation: a semiempirical, subsystem-based correlation functional that predicts a subsystem's two-electron density matrix as a functional of its one-electron density matrix. This model is demonstrated on two small systems: chains of linear, minimal-basis (H-H)(5), treated as a sum of four overlapping (H-H)(2) subsystems; and the aldehyde group of a set of HOC-R molecules. The results provide an initial demonstration of the feasibility of the approach.

Zero-voltage conductance of short gold nanowires

Using the Landauer formula, the conductance of short gold wires is studied. The required electronic structure calculations are performed with a self-consistent tight-binding method. We consider gold wires of single-atom diameter with a variable number (N=1, em leader,5) of atoms. Depending on N, we find considerable conductance variations with one conductance quantum being the upper limit. The results are confirmed by means of Friedel's sum rule. Tip-shaped clusters are used to provide the contact-wire interfaces and the relation between various tip structures and the conductance is discussed. Our predictions about the conductance variations agree qualitatively with new experimental results.

New mechanism for the alpha to omega martensitic transformation in pure titanium

We propose a new direct mechanism for the pressure driven alpha-->omega martensitic transformation in pure titanium. A systematic algorithm enumerates all possible pathways whose energy barriers are evaluated. A new, homogeneous pathway emerges with a barrier at least 4 times lower than other pathways. The pathway is shown to be favorable in any nucleation model.

Disorder and suppression of quantum confinement effects in Pd nanoparticles

Size-selectable ligand-passivated crystalline and amorphous Pd nanoparticles (<4 nm) are synthesized by a novel two-phase process and verified by high-resolution transmission electron microscopy. Scanning tunneling spectroscopy preformed at 5 K on these two types of nanoparticles exhibits clear Coulomb blockade and Coulomb staircases. Size dependent multipeak spectral features in the differential conductance curve are observed for the crystalline Pd particles but not for the amorphous particles. Theoretical analysis shows that these spectral features are related to the quantized electronic states in the crystalline Pd particle. The suppression of the quantum confinement effect in the amorphous particle arises from the reduction of the degeneracy of the eigenstates and the level broadening due to the reduced lifetime of the electronic states.

Slow fluctuations in enhanced Raman scattering and surface roughness relaxation

We propose an explanation for the recently measured slow fluctuations and "blinking" in the surface-enhanced Raman scattering spectrum of single molecules adsorbed on a silver colloidal particle. We suggest that these fluctuations may be related to the dynamic relaxation of the surface roughness on the nanometer scale and show that there are two classes of roughness with qualitatively different dynamics. The predictions agree with the measurements of surface roughness relaxation. Using a theoretical model for the kinetics of surface roughness relaxation in the presence of charges and optical electrical fields, we predict that the high-frequency electromagnetic field increases both the effective surface tension and the surface diffusion constant and thus accelerates the surface smoothing kinetics and time scale of the Raman fluctuations in a manner that is linear with the laser power intensity, while the addition of salt retards the surface relaxation kinetics and increases the time scale of the fluctuations. These predictions are in qualitative agreement with the Raman experiments.

Bcc and Fcc transition metals and alloys: a central role for the Jahn-Teller effect in explaining their ideal and distorted structures

Transition metal elements, alloys, and intermetallic compounds often adopt the body centered cubic (bcc) and face centered cubic (fcc) structures. By comparing quantitative density functional with qualitative tight-binding calculations, we analyze the electronic factors which make the bcc and fcc structures energetically favorable. To do so, we develop a tight-binding function, DeltaE(star), a function that measures the energetic effects of transferring electrons within wave vector stars. This function allows one to connect distortions in solids to the Jahn-Teller effect in molecules and to provide an orbital perspective on structure determining deformations in alloys. We illustrate its use by considering first a two-dimensional square net. We then turn to three-dimensional fcc and bcc structures, and distortions of these. Using DeltaE(star), we rationalize the differences in energy of these structures. We are able to deduce which orbitals are responsible for instabilities in seven to nine valence electron per atom (e(-)/a) bcc systems and five and six e(-)/a fcc structures. Finally we demonstrate that these results account for the bcc and fcc type structures found in both the elements and binary intermetallic compounds of group 4 through 9 transition metal atoms. The outline of a theory of metal structure deformations based on loss of point group operation rather than translational symmetry is presented.

Anisotropic lattice distortions in random alloys from first-principles theory

Within the framework of the exact muffin-tin orbitals (EMTO) theory we have developed a new method to calculate the total energy for random substitutional alloys. The problem of disorder is treated within the coherent potential approximation (CPA), and the total energy is obtained using the full charge density (FCD) technique. The FCD-EMTO-CPA method is suitable for determination of energy changes due to anisotropic lattice distortions in random alloys. In particular, we calculate the elastic constants of the Cu-rich face centered cubic Cu-Zn alloys ( alpha-brass) and optimize the c/a ratio for the hexagonal Zn-rich alloys for both the epsilon and eta phases.

Analytic environment-dependent tight-binding bond integrals: application to MoSi2

We present the first derivation of explicit analytic expressions for the environmental dependence of the sigma, pi, and delta bond integrals within the orthogonal two-center tight-binding approximation by using the recently developed bond-order potential theory to invert the nonorthogonality matrix. We illustrate the power of this new formalism by showing that it not only captures the transferability of the bond integrals between elemental bcc Mo and Si and binary C11(b) MoSi2 but also predicts the absence of any discontinuity between first and second nearest neighbors for the ddsigma bond integral even though large discontinuities exist for ppsigma, pppi, and ddpi.