title: Six-Dimensional Quantum Dynamics of Adsorption and Desorption of H_2 at Pd(100): Steering and Steric Effects

abstract: We report the first six-dimensional quantum dynamical calculations of dissociative adsorption and associative desorption. Using a potential energy surface obtained by density functional theory calculations, we show that the initial decrease of the sticking probability with increasing kinetic energy in the system H_2/Pd(100), which is usually attributed to the existence of a molecular adsorption state, is due to dynamical steering. In addition, we examine the influence of rotational motion and orientation of the hydrogen molecule on adsorption and desorption.


title: The Ab-Initio Simulation of the Liquid Ga-Se System

abstract: Ab-initio dynamical simulation is used to study the liquid Ga-Se system at the three concentrations Ga$_2$Se, GaSe and Ga$_2$Se$_3$ at the temperature 1300~K. The simulations are based on the density functional pseudopotential technique, with the system maintained on the Born-Oppenheimer surface by conjugate gradients minimization. We present results for the partial structure factors and radial distribution functions, which reveal how the liquid structure depends on the composition. Our calculations of the electrical conductivity $\sigma$ using the Kubo-Greenwood approximation show that $\sigma$ depends very strongly on the composition. We show how this variation of $\sigma$ is related to the calculated electronic density of states. Comparisons with recent experimental determinations of the structure and conductivity are also presented.


title: Localization and band-gap pinning in semiconductor superlattices with layer-thickness fluctuations

abstract: We consider (AlAs)_n/(GaAs)_n superlattices with random thickness fluctuations Delta-n around the nominal period n. Using three-dimensional pseudopotential plane-wave band theory, we show that (i) any amount Delta-n/n of thickness fluctuations leads to band-edge wavefunction localization, (ii) for small Delta-n/n the SL band gap is pinned at the gap level produced by a single layer with ``wrong'' thickness n + Delta-n, (iii) the bound states due to monolayer thickness fluctuations lead to significant band-gap reductions, (iv) <001> AlAs/GaAs SL's with monolayer thickness fluctuations have a direct band gap, while the ideal <001> SL's are indirect for n<4.


title: Water Chemisorption and Reconstruction of the MgO Surface

abstract: The observed reactivity of MgO with water is in apparent conflict with theoretical calculations which show that molecular dissociation does not occur on a perfect (001) surface. We have performed ab-initio total energy calculations which show that a chemisorption reaction involving a reconstruction to form a (111) hydroxyl surface is strongly preferred with Delta E = -90.2kJ/mol. We conclude that protonation stabilizes the otherwise unstable (111) surface and that this, not the bare (001), is the most stable surface of MgO under ambient conditions.


title: Interlayer Coupling in Magnetic/Pd Multilayers

abstract: The Anderson model of local-state conduction electron mixing is applied to the p roblem of interlayer magnetic coupling in metallic multilayered structures with palladium (Pd) spacer layers. An oscillation period of $5$ spacer monolayers and the tendency towards ferromag netic bias of the interlayer magnetic coupling that we obtain are consistent with the experimental data.


title: Spin-dependent resonant tunneling through semimetallic ErAs quantum wells

abstract: Resonant tunneling through semimetallic ErAs quantum wells embedded in GaAs structures with AlAs barriers was recently found to exhibit an intriguing behavior in magnetic fields which is explained in terms of tunneling selection rules and the spin-polarized band structure including spin-orbit coupling.


title: Band structure analysis of the conduction-band mass anisotropy in 6H and 4H SiC

abstract: The band structures of 6H and 4H SiC calculated by means of the FP-LMTO method are used to determine the effective mass tensors for their conduction-band minima. The results are shown to be consistent with recent optically detected cyclotron resonance measurements and predict an unusual band filling dependence for 6H-SiC.


title: Comparison of the Electronic Structures and Energetics of Ferroelectric LiNbO3 and LiTaO3

abstract: This paper explains the origin of the ferroelectric instability in LiNbO3 and LiTaO3 and compares the electronic structures and energetics of the two materials.


title: Density-Functional Theory of the nonlinear optical susceptibility: application to cubic semiconductors

abstract: We present a general scheme for the computation of the time dependent (TD) quadratic susceptibility ($\chi^{(2)}$) of an extended insulator obtained by applying the `$2n+1$' theorem to the action functional as defined in TD density functional theory. The resulting expression for $\chi^{(2)}$ includes self-consistent local-field effects, and is a simple function of the linear response of the system. We compute the static $\chi^{(2)}$ of nine III-V and five II-VI semiconductors using the local density approximation(LDA) obtaining good agreement with experiment. For GaP we also evaluate the TD $\chi^{(2)}$ for second harmonic generation using TD-LDA.


title: The Effect of Strain on Thermodynamics of the Weakly First-Order Phase Transition

abstract: Elastic matrix distortion around a growing inclusion of a new phase is analyzed and the associated contribution to the Gibbs free energy is considered. The constant-composition transformation from the parent to product phase is considered within the frame of Landau theory of phase transitions. The volume misfit between the inclusion and matrix is assumed to originate from the transformation volume change coupled with the phenomenological order parameter. The minimization of free energy with respect to the volume change and order parameter gives the dependence of Gibbs energy on the volume fraction of the product phase. The transformation proceeds in a finite temperature region with the equilibrium volume fraction dependent on temperature rather than at a fixed temperature as it would be expected for the first-order transition. The activation processes are shown to be irrelevant and the transformation kinetics is found to be fluctuationless.


title: Towards a Linear-Scaling DFT Technique: The Density Matrix Approach

abstract: A recently proposed linear-scaling scheme for density-functional pseudopotential calculations is described in detail. The method is based on a formulation of density functional theory in which the ground state energy is determined by minimization with respect to the density matrix, subject to the condition that the eigenvalues of the latter lie in the range [0,1]. Linear-scaling behavior is achieved by requiring that the density matrix should vanish when the separation of its arguments exceeds a chosen cutoff. The limitation on the eigenvalue range is imposed by the method of Li, Nunes and Vanderbilt. The scheme is implemented by calculating all terms in the energy on a uniform real-space grid, and minimization is performed using the conjugate-gradient method. Tests on a 512-atom Si system show that the total energy converges rapidly as the range of the density matrix is increased. A discussion of the relation between the present method and other linear-scaling methods is given, and some problems that still require solution are indicated.


title: The Adsorption of H2O on TiO2 and SnO2(110) Studied by First-Principles Calculations

abstract: First-principles calculations based on density functional theory and the pseudopotential method have been used to investigate the energetics of H$_2$O adsorption on the (110) surface of TiO$_2$ and SnO$_2$. Full relaxation of all atomic positions is performed on slab systems with periodic boundary conditions, and the cases of full and half coverage are studied. Both molecular and dissociative (H$_2$O $\rightarrow$ OH$^-$ + H$^+$) adsorption are treated, and allowance is made for relaxation of the adsorbed species to unsymmetrical configurations. It is found that for both TiO$_2$ and SnO$_2$ an unsymmetrical dissociated configuration is the most stable. The symmetrical molecularly adsorbed configuration is unstable with respect to lowering of symmetry, and is separated from the fully dissociated configuration by at most a very small energy barrier. The calculated dissociative adsorption energies for TiO$_2$ and SnO$_2$ are in reasonable agreement with the results of thermal desorption experiments. Calculated total and local electronic densities of states for dissociatively and molecularly adsorbed configurations are presented and their relation with experimental UPS spectra is discussed.


title: The Influence of Gradient Corrections on Bulk and Surface Properties of TiO2 and SnO2

abstract: First-principles calculations based on density functional theory and the pseudo\-potential method have been used to investigate the influence of gradient corrections to the standard LDA technique on the equilibrium structure and energetics of rutile TiO$_2$ and SnO$_2$ perfect crystals and their (110) surfaces. We find that gradient corrections increase the calculated lattice parameters by roughly 3~\%, as has been found for other types of material. Gradient corrections give only very minor changes to the equilibrium surface structure, but reduce the surface energies by about 30~\%.


title: Theoretical Study of One-dimensional Chains of Metal Atoms in Nanotubes

abstract: Using first-principles total-energy pseudopotential calculations, we have studied the properties of chains of potassium and aluminum in nanotubes. For BN tubes, there is little interaction between the metal chains and the tubes, and the conductivity of these tubes is through carriers located at the inner part of the tube. In contrast, for small radius carbon nanotubes, there are two types of interactions: charge-transfer (dominant for alkali atoms) leading to strong ionic cohesion, and hybridization (for multivalent metal atoms) resulting in a smaller cohesion. For Al-atomic chains in carbon tubes, we show that both effects contribute. New electronic properties related to these confined atomic chains of metal are analyzed.


title: Electronic Structure and Valence Band Spectra of Bi4Ti3O12

abstract: The x-ray photoelectron valence band spectrum and x-ray emission valence-band spectra (Ti K _beta_5, Ti L_alpha, O K_alpha) of Bi4Ti3O12 are presented (analyzed in the common energy scale) and interpreted on the basis of a band-structure calculation for an idealized I4/mmm structure of this material.


title: Ab initio study of the volume dependence of dynamical and thermodynamical properties of silicon

abstract: Motivated by the negative thermal expansion observed for silicon between 20 K and 120 K, we present first an ab initio study of the volume dependence of interatomic force constants, phonon frequencies of TA(X) and TA(L) modes, and of the associated mode Gruneisen parameters. The influence of successive nearest neighbors shells is analysed. Analytical formulas, taking into account interactions up to second nearest neighbors, are developped for phonon frequencies of TA(X) and TA(L) modes and the corresponding mode Gruneisen parameters. We also analyze the volume and pressure dependence of various thermodynamic properties (specific heat, bulk modulus, thermal expansion), and point out the effect of the negative mode Gruneisen parameters of the acoustic branches on these properties. Finally, we present the evolution of the mean square atomic displacement and of the atomic temperature factor with the temperature for different volumes, for which the anomalous effects are even greater.


title: Stability of Liquid Crystalline Macroemulsions

abstract: The thermodynamic stability of emulsions of liquid crystal in water (glycerol) matrices is demonstrated for a wide range of materials and concentrations. Coalescence is prevented by an energy barrier for a topological ring defect formation in a neck between the two merging droplets. There is a characteristic size of emulsion droplets, typically tens of microns or more, controlled by the balance of elastic and anchoring energies of the liquid crystal. On removal of liquid crystallinity (by raising the temperature above $T_{ni}$ in thermotropic nematic materials, for example) the energy barriers for coalescence disappear and emulsion droplets can merge quickly, controlled only by the traditional kinetic effects. Practical applications of this effect, as well as some wider theoretical implications are discussed in the end.


title: The de Haas-van Alphen effect in canonical and grand canonical multiband Fermi liquid

abstract: A qualitatively different character of dHvA oscillations has been found in a multiband (quasi)two dimensional Fermi liquid with a fixed fermion density $n_{e}$ (canonical ensemble) compared with an open system where the chemical potential $\mu$ is kept fixed (grand canonical ensemble). A new fundamental period $P_{f}$ appears when $n_{e}$ is fixed, a damping of the Landau levels is relatively small and a background density of states is negligible. $P_{f}$ is determined by the total density rather than by the partial densities of carriers in different bands: $P_{f}=1/(2n_{e}\phi)$ for spin-split Landau levels and $P_{f}=1/(n_{e}\phi)$ in the case of spin degenerate levels where $\phi$ is the flux quantum.


title: Effect of quantum fluctuations on structural phase transitions in SrTiO_3 and BaTiO_3

abstract: Using path-integral Monte Carol simulations and an ab initio effective Hamiltonian, we study the effects of quantum fluctuations on structural phase transitions in the cubic perovskite compounds SrTiO3 and BaTiO3. We find quantum fluctuations affect ferroelectric (FE) transitions more strongly than antiferrodistortive (AFD) ones, even though the effective mass of a single FE local mode is larger. For SrTiO3 we find that the quantum fluctuations suppress the FE transition completely, and reduce the AFD transition temperature from 130K to 110K. For BaTiO3, quantum fluctuations do not affect the order of the transition, but do reduce the transition temperature by 35-50 K. The implications of the calculations are discussed.


title: First-principles investigation of 180-degree domain walls in BaTiO_3

abstract: We present a first-principles study of 180-degree ferroelectric domain walls in tetragonal barium titanate. The theory is based on an effective Hamiltonian that has previously been determined from first-principles ultrasoft-pseudopotential calculations. Statistical properties are investigated using Monte Carlo simulations. We compute the domain-wall energy, free energy, and thickness, analyze the behavior of the ferroelectric order parameter in the interior of the domain wall, and study its spatial fluctuations. An abrupt reversal of the polarization is found, unlike the gradual rotation typical of the ferromagnetic case.


title: Composition Dependence of the Structure and Electronic Properties of Liquid Ga-Se Alloys Studied by Ab Initio Molecular Dynamics Simulation

abstract: Ab initio molecular dynamics simulation is used to study the structure and electronic properties of the liquid Ga-Se system at the three compositions Ga$_2$Se, GaSe and Ga$_2$Se$_3$, and of the GaSe and Ga$_2$Se$_3$ crystals. The calculated equilibrium structure of GaSe crystal agrees well with available experimental data. The neutron-weighted liquid structure factors calculated from the simulations are in reasonable agreement with recent neutron diffraction measurements. Simulation results for the partial radial distribution functions show that the liquid structure is closely related to that of the crystals. A close similarity between solid and liquid is also found for the electronic density of states and charge density. The calculated electronic conductivity decreases strongly with increasing Se content, in accord with experimental measurements.


title: Pressure Dependence of Born Effective Charges, Dielectric Constant and Lattice Dynamics in SiC

abstract: The pressure dependence of the Born effective charge, dielectric constant and zone-center LO and TO phonons have been determined for $3C$-SiC by a linear response method based on the linearized augmented plane wave calculations within the local density approximation. The Born effective charges are found to increase nearly linearly with decreasing volume down to the smallest volume studied, $V/V_0=0.78$, corresponding to a pressure of about 0.8 Mbar. This seems to be in contradiction with the conclusion of the turnover behavior recently reported by Liu and Vohra [Phys.\ Rev.\ Lett.\ {\bf 72}, 4105 (1994)] for $6H$-SiC. Reanalyzing their procedure to extract the pressure dependence of the Born effective charges, we suggest that the turnover behavior they obtained is due to approximations in the assumed pressure dependence of the dielectric constant $\varepsilon_\infty$, the use of a singular set of experimental data for the equation of state, and the uncertainty in measured phonon frequencies, especially at high pressure.


title: Molecular Dynamics Simulation of Electron Trapping in the Sapphire Lattice

abstract: Energy storage and release in dielectric materials can be described on the basis of the charge trapping mechanism. Most phenomenological aspects have been recently rationalized in terms of the space charge model~\cite{blaise,blaise1}. Dynamical aspects are studied here by performing Molecular Dynamics simulations. We show that an excess electron introduced into the sapphire lattice (\alumina) can be trapped only at a limited number of sites. The energy gained by allowing the electron to localize in these sites is of the order of 4-5 eV, in good agreement with the results of the space charge model. Displacements of the neighboring ions due to the implanted charge are shown to be localized in a small region of about 5~\AA. Detrapping is observed at 250 $K$. The ionic displacements turn out to play an important role in modifying the potential landscape by lowering, in a dynamical way, the barriers that cause localization at low temperature.


title: Enhanced electron-phonon coupling at the Mo and W(110) surfaces induced by adsorbed hydrogen

abstract: The possible occurrence of either a charge-density-wave or a Kohn anomaly is governed by the presence of Fermi-surface nesting and the subtle interaction of electrons and phonons. Recent experimental and theoretical investigations suggest such an effect for the hydrogen covered Mo and W(110) surfaces. Using density-functional theory we examine the electronic structure and the electron-phonon coupling of these systems. Besides good agreement with the experimental phonon frequencies our study provides a characterization and quantitative analysis of an unusual scenario determining the electronic, vibrational, and structural properties of these surfaces.


title: Ab initio pseudopotential study of Fe, Co, and Ni employing the spin-polarized LAPW approach

abstract: The ground-state properties of Fe, Co, and Ni are studied with the linear-augmented-plane-wave (LAPW) method and norm-conserving pseudopotentials. The calculated lattice constant, bulk modulus, and magnetic moment with both the local-spin-density approximation (LSDA) and the generalized gradient approximation (GGA) are in good agreement with those of all-electron calculations, respectively. The GGA results show a substantial improvement over the LSDA results, i.e., better agreement with experiment. The accurate treatment of the nonlinear core-valence exchange and correlation interaction is found to be essential for the determination of the magnetic properties of 3d transition metals. The present study demonstrates the successful application of the LAPW pseudopotential approach to the calculation of ground-state properties of magnetic 3d transition metals.


title: Theoretical study of O adlayers on Ru(0001)

abstract: Recent experiments performed at high pressures indicate that ruthenium can support unusually high concentrations of oxygen at the surface. To investigate the structure and stability of high coverage oxygen structures, we performed density functional theory calculations, within the generalized gradient approximation, for O adlayers on Ru(0001) from low coverage up to a full monolayer. We achieve quantitative agreement with previous low energy electron diffraction intensity analyses for the (2x2) and (2x1) phases and predict that an O adlayer with a (1x1) periodicity and coverage of 1 monolayer can form on Ru(0001), where the O adatoms occupy hcp-hollow sites.


title: Unusually Large Thermal Expansion of Ag(111)

abstract: We investigate the thermal behavior of the (111) surface of silver, using phonon frequencies obtained from ab initio total energy calculations, and anharmonic effects treated within a quasiharmonic approximation. Our results reproduce the experimental observation of a large and anomalous increase in the surface thermal expansion at high temperatures [P. Statiris et al., Phys. Rev. Lett. 72, 3574 (1994)]. Surprisingly, we find that this increase can be attributed to a rapid softening of the in-plane phonon frequencies, rather than due to the anharmonicity of the out-of-plane surface phonon modes. This provides evidence for a new mechanism for the enhancement of surface anharmonicity. A comparison with Al(111) shows that the two surfaces behave quite differently, with no evidence for such anomalous behavior on Al(111).


title: Force calculation and atomic-structure optimization for the full-potential linearized augmented plane-wave code WIEN

abstract: Following the approach of Yu, Singh, and Krakauer [Phys. Rev. B 43 (1991) 6411] we extended the linearized augmented plane wave code WIEN of Blaha, Schwarz, and coworkers by the evaluation of forces. In this paper we describe the approach, demonstrate the high accuracy of the force calculation, and use them for an efficient geometry optimization of poly-atomic systems.


title: Parallelisation of algorithms for ab initio computation of material properties

abstract: Parallel algorithms for ab initio calculations of vibrations modes of solids are presented and implemented under PVM. Load balancing and communication problems are dealt with in order to increase parallelism efficiency. For accurate time measurements, synchronisation of processes must be achieved. The results obtained by working with 1,2,4 and 8 processors of a Convex MetaSeries computer are presented, showing a very good efficiency. Further parallelisation of the codes (involving a three-dimensional Fast Fourier Transform) with compiler directives on a Convex Exemplar, is discussed.


title: Potential Energy Surface for H_2 Dissociation over Pd(100)

abstract: The potential energy surface (PES) of dissociative adsorption of H_2 on Pd(100) is investigated using density functional theory and the full-potential linear augmented plane wave (FP-LAPW) method. Several dissociation pathways are identified which have a vanishing energy barrier. A pronounced dependence of the potential energy on ``cartwheel'' rotations of the molecular axis is found. The calculated PES shows no indication of the presence of a precursor state in front of the surface. Both results indicate that steering effects determine the observed decrease of the sticking coefficient at low energies of the H_2 molecules. We show that the topology of the PES is related to the dependence of the covalent H(s)-Pd(d) interactions on the orientation of the H_2 molecule.


title: A model of linear chain submonolayer structures. Application to Li/W(112) and Li/Mo(112)

abstract: We propose a lattice gas model to account for linear chain structures adsorbed on (112) faces of W and Mo. The model includes a dipole-dipole interaction as well as a long-ranged indirect interaction. We have explicitly demonstrated that the periodic ground states depend on a competition between dipole-dipole and indirect interaction. The effect of temperature is studied within the molecular-field approximation. The numerical results show that for dipole-dipole interaction only, all long periodic linear chain phases are suppressed to low temperatures. However, when the long-range indirect interaction becomes important, the long-periodic linear chain phases start to fill up the phase diagram and develop a high thermal stability. Model parameters are chosen to reconstruct a sequence of long-periodic phases as observed experimentally for Li/Mo(112) and Li/W(112).


title: Adlayer core-level shifts of admetal monolayers on transition metal substrates and their relation to the surface chemical reactivity

abstract: Using density-functional-theory we study the electronic and structural properties of a monolayer of Cu on the fcc (100) and (111) surfaces of the late 4d transition metals, as well as a monolayer of Pd on Mo bcc(110). We calculate the ground states of these systems, as well as the difference of the ionization energies of an adlayer core electron and a core electron of the clean surface of the adlayer metal. The theoretical results are compared to available experimental data and discussed in a simple physical picture; it is shown why and how adlayer core-level binding energy shifts can be used to deduce information on the adlayer's chemical reactivity.


title: Scattering of rare-gas atoms at a metal surface: evidence of anticorrugation of the helium-atom potential-energy surface and the surface electron density

abstract: Recent measurements of the scattering of He and Ne atoms at Rh(110) suggest that these two rare-gas atoms measure a qualitatively different surface corrugation: While Ne atom scattering seemingly reflects the electron-density undulation of the substrate surface, the scattering potential of He atoms appears to be anticorrugated. An understanding of this perplexing result is lacking. In this paper we present density functional theory calculations of the interaction potentials of He and Ne with Rh(110). We find that, and explain why, the nature of the interaction of the two probe particles is qualitatively different, which implies that the topographies of their scattering potentials are indeed anticorrugated.


title: The Effective Particle-Hole Interaction and the Optical Response of Simple Metal Clusters

abstract: Following Sham and Rice [L. J. Sham, T. M. Rice, Phys. Rev. 144 (1966) 708] the correlated motion of particle-hole pairs is studied, starting from the general two-particle Greens function. In this way we derive a matrix equation for eigenvalues and wave functions, respectively, of the general type of collective excitation of a N-particle system. The interplay between excitons and plasmons is fully described by this new set of equations. As a by-product we obtain - at least a-posteriori - a justification for the use of the TDLDA for simple-metal clusters.


title: Role of Self-Interaction Effects in the Geometry Optimization of Small Metal Clusters

abstract: By combining the Self-Interaction Correction (SIC) with pseudopotential perturbation theory, the role of self-interaction errors inherent to the Local Density Approximation (LDA) to Density Functional Theory is estimated in the determination of ground state and low energy isomeric structures of small metallic clusters. Its application to neutral sodium clusters with 8 and 20 atoms shows that the SIC provides sizeable effects in Na_8, leading to a different ordering of the low lying isomeric states compared with ab-initio LDA predictions, whereas for Na_20, the SIC effects are less pronounced, such that a quantitative agreement is achieved between the present method and ab-initio LDA calculations.


title: Formation Energy, Stress, and Relaxations of Low-Index Rhodium Surfaces

abstract: Ab initio local-density-functional-theory calculations of formation energies, surface stress, and multilayer relaxations are reported for the (111), (100), and (110) surfaces of Rh. The study is performed using ultrasoft pseudopotentials and plane waves in a parallel implementation. Calculated values are in good agreement with previous studies where they exist. In particular, we confirm recent results on surface stress anisotropy of transition metals.


title: Wetting of the Au(110) substrate: Homoepitaxial islands and layers

abstract: The morphology of epitaxial structures formed in solidification of liquid Au films on the Au(110) surface was studied by molecular dynamics simulation based on a many--body interatomic potential. The (1x2) reconstructed and smooth phase at the temperature T=500 K, as well as a deconstructed and rough phase of the Au(110) substrate at T=900 K, were investigated. The three--dimensional islands with the (111) oriented facets were formed in solidification of thinner liquid films. At the same time, the substrate below these islands underwent a weak faceting process. Conversely, the solidification of thicker liquid films resulted in the flat solid films. The two surface phases of the Au(110) substrate induced different structure of these epitaxial solid films. The films were studied via construction of the Voronoi polyhedra, density profiles, and surface stress calculations.


title: First-Principles Calculation of the Structure of Mercury

abstract: Mercury has perhaps the strangest behavior of any of the metals. Although the other metals in column IIB have an $hcp$ ground state, mercury's ground state is the body centered tetragonal $\beta$Hg phase. The most common phase of mercury is the rhombohedral $\alpha$Hg phase, which is stable from 79K to the melting point and meta-stable below 79K. Another rhombohedral phase, $\gamma$Hg, is believed to exist at low temperatures. First-principles calculations are used to study the energetics of the various phases of mercury. Even when partial spin-orbit effects are included, the calculations indicate that the hexagonal close packed structure is the ground state. It is suggested that a better treatment of the spin-orbit interaction might alter this result.


title: Separable Dual Space Gaussian Pseudo-potentials

abstract: We present pseudo-potential coefficients for the first two rows of the periodic table. The pseudo potential is of a novel analytic form, that gives optimal efficiency in numerical calculations using plane waves as basis set. At most 7 coefficients are necessary to specify its analytic form. It is separable and has optimal decay properties in both real and Fourier space. Because of this property, the application of the nonlocal part of the pseudo-potential to a wave-function can be done in an efficient way on a grid in real space. Real space integration is much faster for large systems than ordinary multiplication in Fourier space since it shows only quadratic scaling with respect to the size of the system. We systematically verify the high accuracy of these pseudo-potentials by extensive atomic and molecular test calculations.


title: Linear-Response Calculation of Electron-Phonon Coupling Parameters

abstract: An ab initio method for calculating electron-phonon coupling parameters is presented. The method is an extension of the plane-wave-based linear-response method for the calculation of lattice dynamics. Results for the mass enhancement parameter and the electron-phonon spectral function for Al, Pb and Li are presented. Comparisons are made to available experimental data.


title: Coulomb interaction and ferroelectric instability of BaTiO3

abstract: Using first-principles calculations, the phonon frequencies at the $\Gamma$ point and the dielectric tensor are determined and analysed for the cubic and rhombohedral phases of BaTiO$_{3}$. The dipole-dipole interaction is then separated \`a la Cochran from the remaining short-range forces, in order to investigate their respective influence on lattice dynamics. This analysis highlights the delicate balance of forces leading to an unstable phonon in the cubic phase and demonstrates the extreme sensitivity of this close compensation to minute effective charge changes. Within our decomposition, the stabilization of the unstable mode in the rhombohedral phase or under isotropic pressure has a different origin.


title: Semiconductor effective charges from tight-binding theory

abstract: We calculate the transverse effective charges of zincblende compound semiconductors using Harrison's tight-binding model to describe the electronic structure. Our results, which are essentially exact within the model, are found to be in much better agreement with experiment than previous perturbation-theory estimates. Efforts to improve the results by using more sophisticated variants of the tight-binding model were actually less successful. The results underline the importance of including quantities that are sensitive to the electronic wavefunctions, such as the effective charges, in the fitting of tight-binding models.


title: Development of a Kohn-Sham like potential in the Self-Consistent Atomic Deformation Model

abstract: This is a brief description of how to derive the local ``atomic'' potentials from the Self-Consistent Atomic Deformation (SCAD) model density function. Particular attention is paid to the spherically averaged case.


title: The Structure, Dynamics and Electronic Structure of Liquid Ag-Se Alloys Investigated by Ab Initio Simulation

abstract: Ab initio molecular-dynamics simulations have been used to investigate the structure, dynamics and electronic properties of the liquid alloy Ag(1-x)Se(x) at 1350 K and at the three compositions x=0.33, 0.42 and 0.65. The calculations are based on density-functional theory in the local density approximation and on the pseudopotential plane-wave method. The reliability of the simulations is confirmed by detailed comparisons with very recent neutron diffraction results for the partial structure factors and radial distribution functions (RDF) of the stoichiometric liquid Ag2Se. The simulations show a dramatic change of the Se-Se RDF with increasing Se content. This change is due to the formation of Se clusters bound by covalent bonds, the Se-Se bond length being almost the same as in pure c-Se and l-Se. The clusters are predominantly chain-like, but for higher x a large fraction of 3-fold coordinated Se atoms is also found. It is shown that the equilibrium fractions of Se present as isolated atoms and in clusters can be understood on a simple charge-balance model based on an ionic interpretation. The Ag and Se diffusion coefficients both increase with Se content, in spite of the Se clustering. An analysis of the Se-Se bond dynamics reveals surprisingly short bond lifetimes of less than 1 ps. The changes in the density of states with composition arise directly from the formation of Se-Se covalent bonds. Results for the electronic conductivity obtained using the Kubo-Greenwood approximation are in adequate agreement with experiment for l-Ag2Se, but not for the high Se contents. Possible reasons for this are discussed.


title: Third order dielectric susceptibility in a model quantum paraelectric

abstract: 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$.


title: Applications of a tight-binding total energy method for transition and noble metals: Elastic Constants, Vacancies, and Surfaces of Monatomic Metals

abstract: A recent tight-binding scheme provides a method for extending the results of first principles calculations to regimes involving $10^2 - 10^3$ atoms in a unit cell. The method uses an analytic set of two-center, non-orthogonal tight-binding parameters, on-site terms which change with the local environment, and no pair potential. The free parameters in this method are chosen to simultaneously fit band structures and total energies from a set of first-principles calculations for monatomic fcc and bcc crystals. To check the accuracy of this method we evaluate structural energy differences, elastic constants, vacancy formation energies, and surface energies, comparing to first-principles calculations and experiment. In most cases there is good agreement between this theory and experiment. We present a detailed account of the method, a complete set of tight-binding parameters, and results for twenty-nine of the alkaline earth, transition and noble metals.


title: Statistical Properties of the Dense Hydrogen Plasma: an ab initio Molecular Dynamics Investigation

abstract: The metallic regime of the hydrogen plasma is studied by ab initio Molecular Dynamics simulations, for classical protons and fully degenerate electrons, in the strong coupling regime of the protons. The breakdown of linear screening observed for decreasing density gives rise to a surprisingly rich low-temperature phase diagram, showing in particular a dramatic drop of the melting temperature of the proton crystal. Extensive dynamical simulations reveal the remarkable persistence of a weakly damped high-frequency ion acoustic mode (plasmon-like), even under conditions of strong electron screening. This collective mode should disappear in the molecular phase, thus providing a probe for the metal-insulator transition in a region of parameters where experiment is difficult to achieve. Finite-size effects arising in the simulation of liquid metals are discussed. The status of the dense Hydrogen matter is extensively reviewed in the introduction.


title: Density functional theory study of Na at Al(111) and O at Ru(0001)

abstract: The success of density functional theory for the description of the adsorption of atoms on surfaces is well established, and based on recent calculations using gradient corrections, it has been shown that it also describes well the dissociative adsorption of molecules at surfaces - admittedly however, the data base for reactions at surfaces is still somewhat small. In the present paper the power of density functional theory calculations is demonstrated by investigations for two different adsorption systems, namely, one with a strongly electropositive adsorbate [Na on Al(111)] and one with a strongly electronegative adsorbate [O on Ru(0001)]. In each case, new hitherto not expected adsorbate phases have been predicted by the theory: For Na on Al(111) the stability of a "four-layer" surface alloy was identified while for O on Ru(0001) it was predicted that the formation of a (1 x 1)-O adlayer should be possible which implies that the apparent saturation coverage of 1/2 is due to kinetic hindering.


title: Dynamical Quantum Processes of Molecular Beams at Surfaces: Dissociative Adsorption of Hydrogen on Metal Surfaces

abstract: Due to the improvement of computer power and the development of efficient algorithms it is now possible to combine high-dimensional quantum dynamical calculations of the dissociative adsorption of molecular beams with reliable ab-initio potential energy surfaces (PES). In this brief review two recent examples of such studies of the systems H_2/Cu(111), where adsorption is hindered by a noticeable energy barrier, and H_2/Pd(100), where activated as well as non-activated paths to adsorption exist, are presented. The effect of lateral surface corrugations on the sticking probability in the tunneling and the classical regime and the role of additional parallel momentum are discussed in the context of the H_2/Cu(111) results. For the system H_2/Pd(100) it is shown that the initial decrease of the sticking probability with increasing kinetic energy, which is usually attributed to a precursor mechanism, can be explained by dynamical steering. In addition, the influence of rotation on the adsorption and desorption dynamics is examined.


title: Six-Dimensional Quantum Dynamics of Adsorption and Desorption of H_2 at Pd(100): No Need for a Molecular Precursor Adsorption State

abstract: We report six-dimensional quantum dynamical calculations of dissociative adsorption and associative desorption of the system H_2/Pd(100) using an ab-initio potential energy surface. We focus on rotational effects in the steering mechanism, which is responsible for the initial decrease of the sticking probability with kinetic energy. In addition, steric effects are briefly discussed.


title: Vibrational Spectra of Defects in Silicon: An Orbital Radii Approach

abstract: A phenomenological approach to the stretching mode vibrational frequencies of defects in semiconductors is proposed. A novel quantum scale is defined in terms of the first principles pseudopotential based orbital radius and the principal quantum number of the element concerned. A universal linear relationship between the Sanderson electronegativity and this quantum scale is established. Next, we show that the stretching mode vibrational frequencies of hydrogen and chlorine in the silicon network scale linearly with this quantum scale. Predictions and identifications of defect environments around the Si-H and Si-Cl are possible. The assignments of vibrational modes in porous silicon are critically examined. We discuss our proposed scale in the context of Mendeleveyan scales in general, and suggest justifications for it. We believe that our approach can be gainfully extended to the vibrational spectra of other semiconductors.

categories: mtrl-th atom-ph chem-ph cond-mat.mtrl-sci


title: Deformation of Small Compressed Droplets

abstract: We investigate the elastic properties of small droplets under compression. The compression of a bubble by two parallel plates is solved exactly and it is shown that a lowest-order expansion of the solution reduces to a form similar to that obtained by Morse and Witten. Other systems are studied numerically and results for configurations involving between 2 and 20 compressing planes are presented. It is found that the response to compression depends on the number of planes. The shear modulus is also calculated for common lattices and the stability crossover between f.c.c.\ and b.c.c.\ is discussed.


title: A Model for the Elasticity of Compressed Emulsions

abstract: We present a new model to describe the unusual elastic properties of compressed emulsions. The response of a single droplet under compression is investigated numerically for different Wigner-Seitz cells. The response is softer than harmonic, and depends on the coordination number of the droplet. Using these results, we propose a new effective inter-droplet potential which is used to determine the elastic response of a monodisperse collection of disordered droplets as a function of volume fraction. Our results are in excellent agreement with recent experiments. This suggests that anharmonicity, together with disorder, are responsible for the quasi-linear increase of $G$ and $\Pi$ observed at $\varphi_c$.


title: First-principles study of stability and vibrational properties of tetragonal PbTiO_3

abstract: A first-principles study of the vibrational modes of PbTiO_3 in the ferroelectric tetragonal phase has been performed at all the main symmetry points of the Brillouin zone (BZ). The calculations use the local-density approximation and ultrasoft pseudopotentials with a plane-wave basis, and reproduce well the available experimental information on the modes at the Gamma point, including the LO-TO splittings. The work was motivated in part by a previously reported transition to an orthorhombic phase at low temperatures [(J. Kobayashi, Y. Uesu, and Y. Sakemi, Phys. Rev. B {\bf 28}, 3866 (1983)]. We show that a linear coupling of orthorhombic strain to one of the modes at Gamma plays a role in the discussion of the possibility of this phase transition. However, no mechanical instabilities (soft modes) are found, either at Gamma or at any of the other high-symmetry points of the BZ.


title: Mechanism of Poisoning the Catalytic Activity of Pd(100) by a Sulfur Adlayer

abstract: The modification of the potential-energy surface (PES) of H_2 dissociation over Pd(100) as induced by the presence of a (2x2) S adlayer is investigated by density-functional theory and the linear augmented plane wave method. It is shown that the poisoning effect of S originates from the formation of energy barriers hindering the dissociation of H_2. The barriers are in the entrance channel of the PES and their magnitude strongly depends on the lateral distance of the H_2 molecule from the S adatoms.


title: On metallic clusters squeezed in atomic cages

abstract: The stability of metallic clusters of sodium (Na) in the octahedral cages of Na-doped fullerites Na6C60 and Na11C60 is discussed within a Thomas-Fermi model. It is shown that the tetrahedral Na4-cluster in Na6C60 has an electric charge of cca. +2.7 (in electron charge units), while the body-centered cubic Na9-cluster in Na11C60 is almost electrically neutral.


title: Modeling of Dislocation Structures in Materials

abstract: A phenomenological model of the evolution of an ensemble of interacting dislocations in an isotropic elastic medium is formulated. The line-defect microstructure is described in terms of a spatially coarse-grained order parameter, the dislocation density tensor. The tensor field satisfies a conservation law that derives from the conservation of Burgers vector. Dislocation motion is entirely dissipative and is assumed to be locally driven by the minimization of plastic free energy. We first outline the method and resulting equations of motion to linear order in the dislocation density tensor, obtain various stationary solutions, and give their geometric interpretation. The coupling of the dislocation density to an externally imposed stress field is also addressed, as well as the impact of the field on the stationary solutions.


title: On The Low-Frequency Vibrational Modes of C$_{60}$

abstract: The vibrational spectrum of C$_{60}$ is compared to the spectrum of a classical isotropic elastic spherical shell. We show correlations between the low frequency modes of C$_{60}$ and those of the spherical shell. We find the spherical model gives the approximate frequency ordering for the low frequency modes. We estimate a Poisson ratio of $\sigma\approx 0.30$ and a transverse speed of sound of $v_s\approx 1800$ m/s for the equivalent elastic shell. We also find that $\omega({\rm M_1})/\omega({\rm M_0})=\sqrt{3\over 2}$ for the shell modes ${\rm M_0}$ and ${\rm M_1}$, independent of elastic constants. We find that this ratio compares favorably with an experimental value of 1.17.


title: Inelastic Collapse of Three Particles

abstract: A system of three particles undergoing inelastic collisions in arbitrary spatial dimensions is studied with the aim of establishing the domain of ``inelastic collapse''---an infinite number of collisions which take place in a finite time. Analytic and simulation results show that for a sufficiently small restitution coefficient, $0\leq r<7-4\sqrt{3}\approx 0.072$, collapse can occur. In one dimension, such a collapse is stable against small perturbations within this entire range. In higher dimensions, the collapse can be stable against small variations of initial conditions, within a smaller $r$ range, $0\leq r<9-4\sqrt{5}\approx 0.056$.


title: Structure, barriers and relaxation mechanisms of kinks in the 90-degree partial dislocation in silicon

abstract: Kink defects in the 90-degree partial dislocation in silicon are studied using a linear-scaling density-matrix technique. The asymmetric core reconstruction plays a crucial role, generating at least four distinct kink species as well as soliton defects. The energies and migration barriers of these entities are calculated and compared with experiment. As a result of certain low-energy kinks, a peculiar alternation of the core reconstruction is predicted. We find the solitons to be remarkably mobile even at very low temperature, and propose that they mediate the kink relaxation dynamics.


title: Modeling of Covalent Bonding in Solids by Inversion of Cohesive Energy Curves

abstract: We provide a systematic test of empirical theories of covalent bonding in solids using an exact procedure to invert ab initio cohesive energy curves. By considering multiple structures of the same material, it is possible for the first time to test competing angular functions, expose inconsistencies in the basic assumption of a cluster expansion, and extract general features of covalent bonding. We test our methods on silicon, and provide the direct evidence that the Tersoff-type bond order formalism correctly describes coordination dependence. For bond-bending forces, we obtain skewed angular functions that favor small angles, unlike existing models. As a proof-of-principle demonstration, we derive a Si interatomic potential which exhibits comparable accuracy to existing models.


title: Generalized Stacking Fault Energy Surfaces and Dislocation Properties of Silicon: A First-Principles Theoretical Study

abstract: The generalized stacking fault (GSF) energy surfaces have received considerable attention due to their close relation to the mechanical properties of solids. We present a detailed study of the GSF energy surfaces of silicon within the framework of density functional theory. We have calculated the GSF energy surfaces for the shuffle and glide set of the (111) plane, and that of the (100) plane of silicon, paying particular attention to the effects of the relaxation of atomic coordinates. Based on the calculated GSF energy surfaces and the Peierls-Nabarro model, we obtain estimates for the dislocation profiles, core energies, Peierls energies, and the corresponding stresses for various planar dislocations of silicon.


title: Slip energy barriers in aluminum and implications for ductile versus brittle behavior

abstract: We conisder the brittle versus ductile behavior of aluminum in the framework of the Peierls-model analysis of dislocation emission from a crack tip. To this end, we perform first-principles quantum mechanical calculations for the unstable stacking energy $\gamma_{us}$ of aluminum along the Shockley partial slip route. Our calculations are based on density functional theory and the local density approximation and include full atomic and volume relaxation. We find that in aluminum $\gamma_{us} = 0.224$ J/m$^2$. Within the Peierls-model analysis, this value would predict a brittle solid which poses an interesting problem since aluminum is typically considered ductile. The resolution may be given by one of three possibilites: (a) Aluminum is indeed brittle at zero temperature, and becomes ductile at a finite temperature due to motion of pre-existing dislocations which relax the stress concentration at the crack tip. (b) Dislocation emission at the crack tip is itself a thermally activated process. (c) Aluminum is actually ductile at all temperatures and the theoretical model employed needs to be significantly improved in order to resolve the apparent contradiction.


title: Magnetic susceptibility of insulators from first principles

abstract: We present an {\it ab initio} approach for the computation of the magnetic susceptibility $\chi$ of insulators. The approach is applied to compute $\chi$ in diamond and in solid neon using density functional theory in the local density approximation, obtaining good agreement with experimental data. In solid neon, we predict an observable dependence of $\chi$ upon pressure.


title: The Melting Curve and Premelting of MgO

abstract: The melting curve for MgO was obtained using molecular dynamics and a non-empirical, many-body potential. We also studied premelting effects by computing the dynamical structure factor in the crystal on approach to melting. The melting curve simulations were performed with periodic boundary conditions with cells up to 512 atoms using the ab-initio Variational Induced Breathing (VIB) model. The melting curve was obtained by computing $% \Delta H_m$ and $\Delta V_m$ and integrating the Clapeyron equation. Our $% \Delta H_m$ is in agreement with previous estimates and we obtain a reasonable $\Delta V_m$, but our melting slope dT/dP (114 K/GPa) is three times greater than that of Zerr and Boehler [1994] (35 K/GPa), suggesting a problem with the experimental melting curve, or an indication of exotic, non-ionic behavior of MgO liquid. We computed $S(q,\omega)$ from simulations of 1000 atom clusters using the Potential Induced Breathing (PIB) model. A low frequency peak in the dynamical structure factor $% S(q,\omega)$ arises below the melting point which appears to be related to the onset of bulk many-atom diffusive exchanges. These exchanges may help destabilize the crystalline state and be related to intrinsic crystalline instability suggested in earlier simulations.


title: Anomalies in He atom scattering spectra of the H-covered Mo(110) and W(110) surfaces

abstract: Helium atom scattering (HAS) studies of the H-covered Mo(110) and W(110) surfaces reveal a twofold anomaly in the respective dispersion curves. In order to explain this unusual behavior we performed density-functional theory calculations of the atomic and electronic structure, the vibrational properties, and the electronic susceptibility of those surfaces. Our work provides evidence for hydrogen adsorption induced Fermi-{}surface nesting. The respective nesting vectors are in excellent agreement with the HAS data and recent angle resolved photoemission experiments of the H-covered alloy system Mo_0.95Re_0.05(110). Also, we investigated the electron-phonon coupling and discovered that the Rayleigh phonon frequency is lowered for those critical wave vectors compared to the clean surfaces. Moreover, the smaller indentation in the HAS spectra can be clearly identified as a Kohn anomaly. Based on our results for the susceptibility and the recently improved understanding of the He scattering mechanism we argue that the larger anomalous dip is due to electron-hole excitations by the He scattering.


title: Semiempirical Hartree-Fock calculations for KNbO3

abstract: In applying the semiempirical intermediate neglect of differential overlap (INDO) method based on the Hartree-Fock formalism to a cubic perovskite-based ferroelectric material KNbO3, it was demonstrated that the accuracy of the method is sufficient for adequately describing the small energy differences related to the ferroelectric instability. The choice of INDO parameters has been done for a system containing Nb. Based on the parametrization proposed, the electronic structure, equilibrium ground state structure of the orthorhombic and rhombohedral phases, and Gamma-TO phonon frequencies in cubic and rhombohedral phases of KNbO3 were calculated and found to be in good agreement with the experimental data and with the first-principles calculations available.


title: Multigrid Methods in Electronic Structure Calculations

abstract: We describe a set of techniques for performing large scale ab initio calculations using multigrid accelerations and a real-space grid as a basis. The multigrid methods provide effective convergence acceleration and preconditioning on all length scales, thereby permitting efficient calculations for ill-conditioned systems with long length scales or high energy cut-offs. We discuss specific implementations of multigrid and real-space algorithms for electronic structure calculations, including an efficient multigrid-accelerated solver for Kohn-Sham equations, compact yet accurate discretization schemes for the Kohn-Sham and Poisson equations, optimized pseudo\-potentials for real-space calculations, efficacious computation of ionic forces, and a complex-wavefunction implementation for arbitrary sampling of the Brillioun zone. A particular strength of a real-space multigrid approach is its ready adaptability to massively parallel computer architectures, and we present an implementation for the Cray-T3D with essentially linear scaling of the execution time with the number of processors. The method has been applied to a variety of periodic and non-periodic systems, including disordered Si, a N impurity in diamond, AlN in the wurtzite structure, and bulk Al. The high accuracy of the atomic forces allows for large step molecular dynamics; e.g., in a 1 ps simulation of Si at 1100 K with an ionic step of 80 a.u., the total energy was conserved within 27 microeV per atom.


title: Phonon softening and superconductivity in tellurium under pressure

abstract: The phonon dispersion and the electron-phonon interaction for the $\beta$-Po and the bcc high pressure phases of tellurium are computed with density-functional perturbation theory. Our calculations reproduce and explain the experimentally observed pressure dependence of the superconducting critical temperature (T$_{\rm c}$) and confirm the connection between the jump in T$_{\rm c}$ and the structural phase transition. The phonon contribution to the free energy is shown to be responsible for the difference in the structural transition pressure observed in low and room temperature experiments.


title: Effects of crack tip geometry on dislocation emission and cleavage: A possible path to enhanced ductility

abstract: We present a systematic study of the effect of crack blunting on subsequent crack propagation and dislocation emission. We show that the stress intensity factor required to propagate the crack is increased as the crack is blunted by up to thirteen atomic layers, but only by a relatively modest amount for a crack with a sharp 60$^\circ$ corner. The effect of the blunting is far less than would be expected from a smoothly blunted crack; the sharp corners preserve the stress concentration, reducing the effect of the blunting. However, for some material parameters blunting changes the preferred deformation mode from brittle cleavage to dislocation emission. In such materials, the absorption of preexisting dislocations by the crack tip can cause the crack tip to be locally arrested, causing a significant increase in the microscopic toughness of the crack tip. Continuum plasticity models have shown that even a moderate increase in the microscopic toughness can lead to an increase in the macroscopic fracture toughness of the material by several orders of magnitude. We thus propose an atomic-scale mechanism at the crack tip, that ultimately may lead to a high fracture toughness in some materials where a sharp crack would seem to be able to propagate in a brittle manner. Results for blunt cracks loaded in mode II are also presented.


title: Structure and Stability of Si(114)-(2x1)

abstract: We describe a recently discovered stable planar surface of silicon, Si(114). This high-index surface, oriented 19.5 degrees away from (001) toward (111), undergoes a 2x1 reconstruction. We propose a complete model for the reconstructed surface based on scanning tunneling microscopy images and first-principles total-energy calculations. The structure and stability of Si(114)-(2x1) arises from a balance between surface dangling bond reduction and surface stress relief, and provides a key to understanding the morphology of a family of surfaces oriented between (001) and (114).


title: Point-charge electrostatics in disordered alloys

abstract: A simple analytic model of point-ion electrostatics has been previously proposed in which the magnitude of the net charge q_i on each atom in an ordered or random alloy depends linearly on the number N_i^(1) of unlike neighbors in its first coordination shell. Point charges extracted from recent large supercell (256-432 atom) local density approximation (LDA) calculations of Cu-Zn random alloys now enable an assessment of the physical validity and accuracy of the simple model. We find that this model accurately describes (i) the trends in q_i vs. N_i^(1), particularly for fcc alloys, (ii) the magnitudes of total electrostatic energies in random alloys, (iii) the relationships between constant-occupation-averaged charges <q_i> and Coulomb shifts <V_i> (i.e., the average over all sites occupied by either $A$ or $B$ atoms) in the random alloy, and (iv) the linear relation between the site charge q_i and the constant- charge-averaged Coulomb shift (i.e., the average over all sites with the same charge) for fcc alloys. However, for bcc alloys the fluctuations predicted by the model in the q_i vs. V_i relation exceed those found in the LDA supercell calculations. We find that (a) the fluctuations present in the model have a vanishing contribution to the electrostatic energy. (b) Generalizing the model to include a dependence of the charge on the atoms in the first three (two) shells in bcc (fcc) - rather than the first shell only - removes the fluctuations, in complete agreement with the LDA data. We also demonstrate an efficient way to extract charge transfer parameters of the generalized model from LDA calculations on small unit cells.


title: Inhibited Al diffusion and growth roughening on Ga-coated Al (100)

abstract: Ab initio calculations indicate that the ground state for Ga adsorption on Al (100) is on-surface with local unit coverage. On Ga-coated Al (100), the bridge diffusion barrier for Al is large, but the Al$\rightarrow$Ga {\it exchange barrier is zero}: the ensuing incorporation of randomly deposited Al's into the Ga overlayer realizes a percolation network, efficiently recoated by Ga atoms. Based on calculated energetics, we predict rough surface growth at all temperatures; modeling the growth by a random deposition model with partial relaxation, we find a power-law divergent roughness $w\sim t^{\,0.07\pm0.02}$.


title: Gross and Scheffler Reply

abstract: Among other things (e.g. steering and steric effects in dissociative adsorption) we had predicted that the initial sticking probability of H_2 molecules impinging at clean Pd(100) exhibits oscillations, reflecting the quantum nature of the scattering process. In the preceding comment Rettner and Auerbach (RA) analyze experimental results and conclude that these oscillations are not detectable and thus either not existing or at least very small. In this reply we argue that the experimental study of RA is not conclusive to rule out the existence of quantum oscillations in the scattering of H_2 and note several problems and incongruities in their study.


title: Potential, core-level and d band shifts at transition metal surfaces

abstract: We have extended the validity of the correlation between the surface 3d-core-level shift (SCLS) and the surface d band shift (SDBS) to the entire 4d transition metal series and to the neighboring elements Sr and Ag via accurate first-principles calculations. We find that the correlation is quasilinear and robust with respect to the differencies both between initial and final-state calculations of the SCLS's and two distinct measures of the SDBS's. We show that despite the complex spatial dependence of the surface potential shift (SPS) and the location of the 3d and 4d orbitals in different regions of space, the correlation exists because the sampling of the SPS by the 3d and 4d orbitals remains similar. We show further that the sign change of the SCLS's across the transition series does indeed arise from the d band-narrowing mechanism previously proposed. However, while in the heavier transition metals the predicted increase of d electrons in the surface layer relative to the bulk arises primarily from transfers from s and p states to d states within the surface layer, in the lighter transition metals the predicted decrease of surface d electrons arises primarily from flow out into the vacuum.


title: The GaAs Equilibrium Crystal Shape from First-Principles

abstract: Surface energies for different GaAs surface orientations have been calculated as a function of the chemical potential. We use an energy density formalism within the first-principles pseudopotential density-functional approach. The equilibrium crystal shape (ECS) has been derived from the surface energies for the (110), (100), (111), and (-1-1-1) orientations. Under As-rich conditions all four considered surface orientations exist in thermodynamic equilibrium, in agreement with experimental observations. Moreover, our calculations allow us to decide on previous contradictory theoretical values for the surface energies of the (111) and (-1-1-1) facets.


title: Influence of molecular vibrations on dissociative adsorption

abstract: The influence of molecular vibrations on dissociative adsorption is studied by six-dimensional quantum dynamical calculations. For the system H_2 at Pd(100), which possesses non-activated pathways, it is shown that large vibrational effects exist and that they are not due to a strongly curved reaction path and a late dissociation-hindering minimum barrier, as was previously assumed. Instead, they are caused by the lowering of the H-H vibrational frequency during the dissociation and the multi-dimensionality of the potential energy surface. Still there are quantitative discrepancies between theory and experiment identified.


title: Anisotropy of Growth of the Close-Packed Surfaces of Silver

abstract: The growth morphology of clean silver exhibits a profound anisotropy: The growing surface of Ag(111) is typically very rough while that of Ag(100) is smooth and flat. This serious and important difference is unexpected, not understood, and hitherto not observed for any other metal. Using density functional theory calculations of self-diffusion on flat and stepped Ag(100) we find, for example, that at flat regions a hopping mechanism is favored, while across step edges diffusion proceeds by an exchange process. The calculated microscopic parameters explain the experimentally reported growth properties.


title: Reconstructions and faceting of H-covered Al(111)

abstract: First principles calculations show that H overlayers induce a honeycomb vacancy reconstruction on Al(111). Above one ML H coverage, Al(111) facets into surfaces vicinal to (111), with {100} faceted single and double steps and vacancies on the {111} terraces. These H induced reconstructions are stable because H binds more strongly to {100} than to {111} facets. Faceting of Al(111) explains step bunching in H-mediated epitaxy of Al(111). H-covered Al(100) is stable against faceting.


title: Variational quantum Monte Carlo ground state of GaAs

abstract: Variational quantum Monte Carlo calculations are reported for the bulk GaAs semiconductor in order to present values for the ground-state energy, the lattice constant, the bulk modulus, and some derived properties. The statistical accuracy is significantly higher than the remaining differences to the experimental values, especially for the total-energy upper bound. The agreement with experiment is satisfactory. The results are also compared with those of density functional calculations. The accuracy of our results is comparable to the best of those calculations.


title: Scattering of hydrogen molecules from a reactive surface: Strong off-specular and rotationally inelastic diffraction

abstract: Six-dimensional quantum dynamical calculations of the scattering of H_2 from a Pd(100) surface using a potential energy surface derived from density-functional theory calculations are presented. Due to the corrugation and anisotropy of the PES strong off-specular and rotationally inelastic diffraction is found. The dependence of the diffraction intensitities on the incident kinetic energy is closely examined. In particular we focus on the quantum oscillations for normal and off-normal incidence.


title: Ab initio Molecular Dynamics Study of D_2 Desorption from Si(100)

abstract: Ab initio molecular dynamics calculations of deuterium desorbing from Si(100) have been performed in order to monitor the energy redistribution among the hydrogen and silicon degrees of freedom during the desorption process. The calculations show that part of the potential energy at the transition state to desorption is transferred to the silicon lattice. The deuterium molecules leave the surface vibrationally hot and rotationally cold, in agreement with experiments; the mean kinetic energy, however, is larger than found in experiments.


title: The equilibrium shape of InAs quantum dots grown on a GaAs(001) substrate

abstract: The equilibrium shape of strained InAs quantum dots grown epitaxially on a GaAs(001) substrate is derived as a function of volume. InAs surface energies are calculated within density-functional theory, and a continuum approach is applied for the elastic relaxation energies.


title: Local Isoelectronic Reactivity of Solid Surfaces

abstract: The quantity w^N(r) = ( 1/ k^2 T_el)[partial n(r, T_el) / partial T_el]_(v(r),N) is introduced as a convenient measure of the local isoelectronic reactivity of surfaces. It characterizes the local polarizability of the surface and it can be calculated easily. The quantity w^N(r) supplements the charge transfer reactivity measured e.g. by the local softness to which it is closely related. We demonstrate the applicability and virtues of the function w^N(r) for the example of hydrogen dissociation and adsorption on Pd(100).


title: Structure and stability of a high-coverage (1x1) oxygen phase on Ru(0001)

abstract: The formation of chemisorbed O-phases on Ru(0001) by exposure to O_2 at low pressures is apparently limited to coverages Theta <= 0.5. Using low-energy electron diffraction and density functional theory we show that this restriction is caused by kinetic hindering and that a dense O overlayer (Theta = 1) can be formed with a (1x1) periodicity. The structural and energetic properties of this new adsorbate phase are analyzed and discussed in view of attempts to bridge the so-called "pressure gap" in heterogeneous catalysis. It is argued that the identified system actuates the unusually high rate of oxidizing reactions at Ru surfaces under high oxygen pressure conditions.


title: Microscopic Estimates for Electromigration Velocities of Intragranular Voids in Aluminum Lines

abstract: We explore the effect of faceting on possible mechanisms for mass transport around electromigration voids in aluminum interconnects. Motivated by linear response estimates which suggest that particle flux would be much higher along steps than across terraces on a clean aluminum surface, we study step nucleation in the presence of a small driving force along a surface. We find that step nucleation, even on a nearly defect-free void surface, would be slow if the step energy is equal to that calculated for a clean aluminum surface. In the presence of a uniform electromigration force, the creation of new steps between existing ones should not occur unless the free energy cost of a step is much less than thermal energies. We conclude that voids cannot move intragranularly at $\mu$m/hr rates without help from other factors such as local heating and impurities.


title: Ab initio statistical mechanics of the ferroelectric phase transition in PbTiO3

abstract: An effective Hamiltonian for the ferroelectric transition in PbTiO3 is constructed from first-principles density-functional-theory total-energy and linear-response calculations through the use of a localized, symmetrized basis set of ``lattice Wannier functions.'' Explicit parametrization of the polar lattice Wannier functions is used for subspace projection, addressing the issues of LO-TO splitting and coupling to the complementary subspace. In contrast with ferroelectric BaTiO3 and KNbO3, we find significant involvement of the Pb atom in the lattice instability. Monte Carlo simulations for this Hamiltonian show a first-order cubic-tetragonal transition at 660 K. Resulting temperature dependence of spontaneous polarization, c/a ratio and unit-cell volume near the transition are in good agreement with experiment. Comparison of Monte Carlo results with mean field theory analysis shows that both strain and fluctuations are necessary to produce the first order character of this transition.


title: Ab-initio study of the anomalies in the He atom scattering spectra of H/Mo(110) and H/W(110)

abstract: Helium atom scattering (HAS) studies of the H-covered Mo(110) and W(110) surfaces reveal a twofold anomaly in the respective dispersion curves. In order to explain this unusual behavior we performed density functional theory calculations of the atomic and electronic structure, the vibrational properties, and the spectrum of electron-hole excitations of those surfaces. Our work provides evidence for hydrogen adsorption induced Fermi surface nesting. The respective nesting vectors are in excellent agreement with the HAS data and recent angle resolved photoemission experiments of the H-covered alloy system Mo_0.95Re_0.05(110). Also, we investigated the electron-phonon coupling and discovered that the Rayleigh phonon frequency is lowered for those critical wave vectors. Moreover, the smaller indentation in the HAS spectra can be clearly identified as a Kohn anomaly. Based on our results for the susceptibility and the recently improved understanding of the He scattering mechanism we argue that the larger anomalous dip is due to a direct interaction of the He atoms with electron-hole excitations at the Fermi level.


title: Polymers in Curved Boxes

abstract: We apply results derived in other contexts for the spectrum of the Laplace operator in curved geometries to the study of an ideal polymer chain confined to a spherical annulus in arbitrary space dimension D and conclude that the free energy compared to its value for an uncurved box of the same thickness and volume, is lower when $D < 3$, stays the same when $D = 3$, and is higher when \mbox{$D > 3$}. Thus confining an ideal polymer chain to a cylindrical shell, lowers the effective bending elasticity of the walls, and might induce spontaneous symmetry breaking, i.e. bending. (Actually, the above mentioned results show that {\em {any}} shell in $D = 3$ induces this effect, except for a spherical shell). We compute the contribution of this effect to the bending rigidities in the Helfrich free energy expression.


title: Anomalous effective charges and far IR optical absorption of Al$_2$Ru from first principles

abstract: For the orthorhombic intermetallic semiconductor Al$_2$Ru, the bandstructure, valence charge density, zone center optical phonon frequencies, and Born effective charge and electronic dielectric tensors are calculated using variational density functional perturbation theory with {\em ab initio} pseudopotentials and a plane wave basis set. Good agreement is obtained with recent measurements on polycrystalline samples which showed anomalously strong far IR absorption by optical phonons, while analysis of the valence charge density shows that the static ionic charges of Al and Ru are negligible. Hybridization is proposed as the single origin both of the semiconducting gap and the anomalous Born effective charges. Analogous behavior is expected in related compounds such as NiSnZr, PbTe, skutterudites, and Al-transition-metal quasicrystals.


title: Reciprocity Theorems for Ab Initio Force Calculations

abstract: We present a method for calculating ab initio interatomic forces which scales quadratically with the size of the system and provides a physically transparent representation of the force in terms of the spatial variation of the electronic charge density. The method is based on a reciprocity theorem for evaluating an effective potential acting on a charged ion in the core of each atom. We illustrate the method with calculations for diatomic molecules.


title: Band Structure and Transport Properties of CrO_2

abstract: Local Spin Density Approximation (LSDA) is used to calculate the energy bands of both the ferromagnetic and paramagnetic phases of metallic CrO_2. The Fermi level lies in a peak in the paramagnetic density of states, and the ferromagnetic phase is more stable. As first predicted by Schwarz, the magnetic moment is 2 \mu_B per Cr atom, with the Fermi level for minority spins lying in an insulating gap between oxygen p and chromium d states ("half-metallic" behavior.) The A_1g Raman frequency is predicted to be 587 cm^{-1}. Drude plasma frequencies are of order 2eV, as seen experimentally by Chase. The measured resistivity is used to find the electron mean-free path l, which is only a few angstroms at 600K, but nevertheless, resistivity continues to rise as temperature increases. This puts CrO_2 into the category of "bad metals" in common with the high T_c superconductors, the high T metallic phase of VO_2, and the ferromagnet SrRuO_3. In common with both SrRuO_3 and Sr_2RuO_4, the measured specific heat \gamma is higher than band theory by a renormalization factor close to 4.


title: Reconstruction and thermal stability of the cubic SiC(001) surfaces

abstract: The (001) surfaces of cubic SiC were investigated with ab-initio molecular dynamics simulations. We show that C-terminated surfaces can have different c(2x2) and p(2x1) reconstructions, depending on preparation conditions and thermal treatment, and we suggest experimental probes to identify the various reconstructed geometries. Furthermore we show that Si-terminated surfaces exhibit a p(2x1) reconstruction at T=0, whereas above room temperature they oscillate between a dimer row and an ideal geometry below 500 K, and sample several patterns including a c(4x2) above 500 K.


title: Ab initio simulations of liquid systems: Concentration dependence of the electric conductivity of NaSn alloys

abstract: Liquid NaSn alloys in five different compositions (20, 40, 50, 57 and 80% sodium) are studied using density functional calculations combined with molecular dynamics(Car-Parrinello method). The frequency-dependent electric conductivities for the systems are calculated by means of the Kubo-Greenwood formula. The extrapolated DC conductivities are in good agreement with the experimental data and reproduce the strong variation with the concentration. The maximum of conductivity is obtained, in agreement with experiment, near the equimolar composition. The strong variation of conductivity, ranging from almost semiconducting up to metallic behaviour, can be understood by an analysis of the densities-of-states.


title: Correlation effects in MgO and CaO: Cohesive energies and lattice constants

abstract: A recently proposed computational scheme based on local increments has been applied to the calculation of correlation contributions to the cohesive energy of the CaO crystal. Using ab-initio quantum chemical methods for evaluating individual increments, we obtain 80% of the difference between the experimental and Hartree-Fock cohesive energies. Lattice constants corrected for correlation effects deviate by less than 1% from experimental values, in the case of MgO and CaO.


title: Surface Effects in Ferroelectrics: Periodic Slab Computations for BaTiO_3

abstract: Total energies, electronic structure, surface energies, polarization, potentials and charge densities were studied for slabs of BaTiO_3 using the Linearized Augmented Plane Wave (LAPW) method. The depolarization field inhibits ferroelectricity in the slabs, and the macroscopic field set up across a ferroelectric slab is sufficient to cause electronic states to span the gap and give a metallic band structure, but the band shifts are not rigid and O p states tend to pile up at the Fermi level. There are electronic surface states, especially evident on TiO_2 surfaces. The dangling bonds bond back to the surface Ti's and make the surface stable and reactive. The BaO surfaces are more ionic than the bulk.


title: Basis Functions for Linear-Scaling First-Principles Calculations

abstract: In the framework of a recently reported linear-scaling method for density-functional-pseudopotential calculations, we investigate the use of localized basis functions for such work. We propose a basis set in which each local orbital is represented in terms of an array of `blip functions'' on the points of a grid. We analyze the relation between blip-function basis sets and the plane-wave basis used in standard pseudopotential methods, derive criteria for the approximate equivalence of the two, and describe practical tests of these criteria. Techniques are presented for using blip-function basis sets in linear-scaling calculations, and numerical tests of these techniques are reported for Si crystal using both local and non-local pseudopotentials. We find rapid convergence of the total energy to the values given by standard plane-wave calculations as the radius of the linear-scaling localized orbitals is increased.


title: Substrate-adsorbate coupling in CO-adsorbed copper

abstract: The vibrational properties of carbon monoxide adsorbed to the copper (100) surface are explored within density functional theory. Atoms of the substrate and adsorbate are treated on an equal footing in order to examine the effect of substrate--adsorbate coupling. This coupling is found to have a significant effect on the vibrational modes, particularly the in-plane frustrated translation, which mixes strongly with substrate phonons and broadens into a resonance. The predicted lifetime due to this harmonic decay mechanism is in excellent quantitative agreement with experiment.


title: Ab initio molecular dynamics study of the desorption of D_2 from Si(100)

abstract: Ab initio molecular dynamics calculations of deuterium desorbing from Si(100) have been performed in order to monitor the energy redistribution among the various D$_2$ and silicon degrees of freedom during the desorption process. The calculations show that a considerable part of the potential energy at the transition state to desorption is transferred to the silicon lattice. The deuterium molecules leave the surface vibrationally hot and rotationally cold, in agreement with thermal desorption experiments; the mean kinetic energy, however, is larger than found in a laser-induced desorption experiment. We discuss possible reasons for this discrepancy.


title: Role of defects and impurities in doping of GaN

abstract: We have calculated formation energies and position of the defect levels for all native defects and for a variety of donor and acceptor impurities employing first-principles total-energy calculations. An analysis of the numerical results gives direct insight into defect concentrations and impurity solubility with respect to growth parameters (temperature, chemical potentials) and into the mechanisms limiting the doping levels in GaN. We show how compensation and passivation by native defects or impurities, solubility issues, and incorporation of dopants on other sites influence the acceptor doping levels.


title: Reconstructions of Ir (110) and (100): an ab initio study

abstract: Prediction criteria for surface reconstructions are discussed, with reference to ab initio calculations of the (110)-$1\times 2$ missing-row and (100)-$5\times 1$ quasi-hexagonal reconstructions of Ir and Rh.


title: Strained tetragonal states and Bain paths in metals

abstract: Paths of tetragonal states between two phases of a material, such as bcc and fcc, are called Bain paths. Two simple Bain paths can be defined in terms of special imposed stresses, one of which applies directly to strained epitaxial films. Each path goes far into the range of nonlinear elasticity and reaches a range of structural parameters in which the structure is inherently unstable. In this paper we identify and analyze the general properties of these paths by density functional theory. Special examples include vanadium, cobalt and copper, and the epitaxial path is used to identify an epitaxial film as related uniquely to a bulk phase.


title: Steering and isotope effects in the dissociative adsorption of H_2/Pd(100)

abstract: The interaction of hydrogen with many transition metal surfaces is characterized by a coexistence of activated with non-activated paths to adsorption with a broad distribution of barrier heights. By performing six-dimensional quantum dynamical and classical molecular dynamics calculations using the same potential energy surface derived from ``ab initio'' calculations for the system H_2/Pd(100) we show that these features of the potential energy surface lead to strong steering effects in the dissociative adsorption dynamics. The adsorption dynamics shows only a small isotope effect which is purely due to the quantum nature of hydrogen.


title: Diffusivity of Ga and Al adatoms on GaAs(001)

abstract: The diffusivity of Ga and Al adatoms on the (2x4) reconstructed GaAs(001) surface are evaluated using detailed ab initio total energy calculations of the potential energy surface together with transition state theory. A strong diffusion anisotropy is found, with the direction of fastest diffusion being parallel to the surface As-dimer orientation. In contrast to previous calculations we identify a short--bridge position between the two As atoms of a surface dimer as the adsorption site for Al and Ga adatoms.


title: A Model for the Thermal Expansion of Ag(111) and other Metal Surfaces

abstract: We develop a model to study the thermal expansion of surfaces, wherein phonon frequencies are obtained from ab initio total energy calculations. Anharmonic effects are treated exactly in the direction normal to the surface, and within a quasiharmonic approximation in the plane of the surface. We apply this model to the Ag(111) and Al(111) surfaces, and find that our calculations reproduce the experimental observation of a large and anomalous increase in the surface thermal expansion of Ag(111) at high temperatures [P. Statiris, H.C. Lu and T. Gustafsson, Phys. Rev. Lett. 72, 3574 (1994)]. Surprisingly, we find that this increase can be attributed to a rapid softening of the in-plane phonon frequencies, rather than due to the anharmonicity of the out-of-plane surface phonon modes. This provides evidence for a new mechanism for the enhancement of surface anharmonicity. A comparison with Al(111) shows that the two surfaces behave quite differently, with no evidence for such anomalous behavior on Al(111).


title: Theory of Alkali Induced Reconstruction of the Cu(100) Surface

abstract: LEED experiments show that Li adsorbed at Cu(100) surfaces at room temperature induces a (2x1) missing row substrate reconstruction while adsorption at lower temperatures, T=180 K, results in an unreconstructed Cu(100)+c(2x2)--Li overlayer structure. Substrate reconstruction has not been observed for Na nor for K adsorption. In order to study the specific reconstruction behavior of the Li adsorbate ab initio DFT calculations have been performed on Cu(100)+Ad, Ad = Li, Na, K systems at coverages Theta_Ad=0.25-0.5 with and without reconstruction. The calculations show that the (2x1) MR reconstructed surface lies energetically above the ideal (1x1) surface by 0.2 eV per unit cell. However, alkali binding is stronger in the MR geometry as compared to that of the ideal surface where the increase in bond strength becomes smaller in going from Li to Na to K. As a result, the MR reconstructed and the overlayer adsorbate systems are energetically very close for Cu(100)+Li while for Na and K the overlayer geometry is always favored.


title: Fermionic representation for basic hypergeometric functions related to Schur polynomials

abstract: We present the fermionic representation for the q-deformed hypergeometric functions related to Schur polynomials considered by S.Milne \cite{Milne}. For $q=1$ these functions are also known as hypergeometric functions of matrix argument which are related to zonal spherical polynomials for $GL(N,C)/U(N)$ symmetric space. We show that these multivariable hypergeometric functions are tau-functions of the KP hierarchy. At the same time they are the ratios of Toda lattice tau-functions considered by Takasaki in \cite{Tinit}, \cite{T} evaluated at certain values of higher Toda lattice times. The variables of the hypergeometric functions are related to the higher times of those hierarchies via Miwa change of variables. The discrete Toda lattice variable shifts parameters of hypergeometric functions. Hypergeometric functions of type ${}_pF_s$ can be also viewed as group 2-cocycle for the $\Psi$DO on the circle of the order $p-s \leq 1$ (the group times are higher times of TL hierarchy and the arguments of hypergeometric function). We get the determinant representation and the integral representation of special type of KP tau-functions, these results generalize some of Milne's results in \cite{Milne}. We write down a system of linear differential and difference equations for these tau-functions (string equations). We present also fermionic representation for special type of Gelfand-Graev hypergeometric functions.


title: Local estimates for entropy densities in coupled map lattices

abstract: We present a method to derive an upper bound for the entropy density of coupled map lattices with local interactions from local observations. To do this, we use an embedding technique being a combination of time delay and spatial embedding. This embedding allows us to identify the local character of the equations of motion. Based on this method we present an approximate estimate of the entropy density by the correlation integral.


title: On Soliton-type Solutions of Equations Associated with N-component Systems

abstract: The algebraic geometric approach to $N$-component systems of nonlinear integrable PDE's is used to obtain and analyze explicit solutions of the coupled KdV and Dym equations. Detailed analysis of soliton fission, kink to anti-kink transitions and multi-peaked soliton solutions is carried out. Transformations are used to connect these solutions to several other equations that model physical phenomena in fluid dynamics and nonlinear optics.


title: On Peakon Solutions of the Shallow Water Equation

abstract: A new parameterization of the Jacobi inversion problem is used along with the dynamics of the peaks to describe finite time interaction of peakon weak solutions of the Shallow Water equation.


title: One-Two Dimensional Nonlinear Pulse Interaction

abstract: The peculiar intergrability of the Davey-Stewartson equation allows us to find analytically solutions describing the simultaneous formation and interaction of one-dimensional and two-dimensional localized coherent structures. The predicted phenomenology allows us to address the issue of interaction of solitons of different dimensionality that may serve as a starting point for the understanding of hybrido-dimensional collisions recently observed in nonlinear optical media.


title: Q-deformed solitons and quantum solitons of the Maxwell-Bloch lattice

abstract: We report for the first time exact solutions of a completely integrable nonlinear lattice system for which the dynamical variables satisfy a q-deformed Lie algebra - the Lie-Poisson algebra su_q(2). The system considered is a q-deformed lattice for which in continuum limit the equations of motion become the envelope Maxwell-Bloch (or SIT) equations describing the resonant interaction of light with a nonlinear dielectric. Thus the N-soliton solutions we here report are the natural q-deformations, necessary for a lattice, of the well-known multi-soliton and breather solutions of self-induced transparency (SIT). The method we use to find these solutions is a generalization of the Darboux-Backlund dressing method. The extension of these results to quantum solitons is sketched.


title: Integrable Structure of Interface Dynamics

abstract: We establish the equivalence of a 2D contour dynamics to the dispersionless limit of the integrable Toda hierarchy constrained by a string equation. Remarkably, the same hierarchy underlies 2D quantum gravity.

categories: nlin.SI cond-mat.stat-mech hep-th nlin.PS


title: Relaxation oscillations and negative strain rate sensitivity in the Portevin - Le Chatelier effect

abstract: A characteristic feature of the Portevin - Le Chatelier effect or the jerky flow is the stick-slip nature of stress-strain curves which is believed to result from the negative strain rate dependence of the flow stress. The latter is assumed to result from the competition of a few relevant time scales controlling the dynamics of jerky flow. We address the issue of time scales and its connection to the negative strain rate sensitivity of the flow stress within the framework of a model for the jerky flow which is known to reproduce several experimentally observed features including the negative strain rate sensitivity of the flow stress. We attempt to understand the above issues by analyzing the geometry of the slow manifold underlying the relaxational oscillations in the model. We show that the nature of the relaxational oscillations is a result of the atypical bent geometry of the slow manifold. The analysis of the slow manifold structure helps us to understand the time scales operating in different regions of the slow manifold. Using this information we are able to establish connection with the strain rate sensitivity of the flow stress. The analysis also helps us to provide a proper dynamical interpretation for the negative branch of the strain rate sensitivity.

categories: nlin.CD cond-mat.mtrl-sci math-ph math.MP


title: Some Techniques for the Measurement of Complexity in Tierra

abstract: Recently, Adami and coworkers have been able to measure the information content of digital organisms living in their {\em Avida} artificial life system. They show that over time, the organisms behave like Maxwell's demon, accreting information (or complexity) as they evolve. In {\em Avida} the organisms don't interact with each other, merely reproduce at a particular rate (their fitness), and attempt to evaluate an externally given arithmetic function in order win bonus fitness points. Measuring the information content of a digital organism is essentially a process of counting the number of genotypes that give rise to the same phenotype. Whilst Avidan organisms have a particularly simple phenotype, Tierran organisms interact with each other, giving rise to an ecology of phenotypes. In this paper, I discuss techniques for comparing pairs of Tierran organisms to determine if they are phenotypically equivalent. I then discuss a method for computing an estimate of the number of phenotypically equivalent genotypes that is more accurate than the ``hot site'' estimate used by Adami's group. Finally, I report on an experimental analysis of a Tierra run.


title: A new method to introduce additional separated variables for high-order binary constrained flows

abstract: Degrees of freedom for high-order binary constrained flows of soliton equations admitting $2\times 2$ Lax matrices are $2N+k_0$. It is known that $N+k_0$ pairs of canonical separated variables for their separation of variables can be introduced directly via their Lax matrices. In present paper we propose a new method to introduce the additional $N$ pairs of canonical separated variables and $N$ additional separated equations. The Jacobi inversion problems for high-order binary constrained flows and for soliton equations are also established. This new method can be applied to all high-order binary constrained flows admitting $2\times 2$ Lax matrices.


title: Lyapunov Instability for a hard-disk fluid in equilibrium and nonequilibrium thermostated by deterministic scattering

abstract: We compute the full Lyapunov spectra for a hard-disk fluid under temperature gradient and shear. The system is thermalized by deterministic and time-reversible scattering at the boundary. This thermostating mechanism allows for energy fluctuations around a mean value which is reflected by only two vanishing Lyapunov exponents in equilibrium and nonequilibrium. The Lyapunov exponents are calculated with a recently developed formalism for systems with elastic hard collisions. In a nonequilibrium steady state the average phase space volume is contracted onto a fractal attractor leading to a negative sum of Lyapunov exponents. Since the system is driven inhomogeneously we do not expect the conjugate pairing rule to hold which is confirmed numerically.


title: (2+0)-Dimensional Integrable Equations and Exact Solutions

abstract: We propose a nonlinear $\sigma$-model in a curved space as a general integrable elliptic model. We construct its exact solutions and obtain energy estimates near the critical point. We consider the Pohlmeyer transformation in Euclidean space and investigate the gauge equivalence conditions for a broad class of elliptic equations. We develop the inverse scattering transform method for the $\sinh$-Gordon equation and evaluate its exact and asymptotic solutions.


title: Frobenius-Perron Resonances for Maps with a Mixed Phase Space

abstract: Resonances of the time evolution (Frobenius-Perron) operator P for phase space densities have recently been shown to play a key role for the interrelations of classical, semiclassical and quantum dynamics. Efficient methods to determine resonances are thus in demand, in particular for Hamiltonian systems displaying a mix of chaotic and regular behavior. We present a powerful method based on truncating P to a finite matrix which not only allows to identify resonances but also the associated phase space structures. It is demonstrated to work well for a prototypical dynamical system.


title: Asymmetric Squares as Standing Waves in Rayleigh-Benard Convection

abstract: Possibility of asymmetric square convection is investigated numerically using a few mode Lorenz-like model for thermal convection in Boussinesq fluids confined between two stress free and conducting flat boundaries. For relatively large value of Rayleigh number, the stationary rolls become unstable and asymmetric squares appear as standing waves at the onset of secondary instability. Asymmetric squares, two dimensional rolls and again asymmetric squares with their corners shifted by half a wavelength form a stable limit cycle.


title: The statistical properties of the city transport in Cuernavaca (Mexico) and Random matrix ensembles

abstract: We analyze statistical properties of the city bus transport in Cuernavaca (Mexico) and show that the bus arrivals display probability distributions conforming those given by the Unitary Ensemble of random matrices.