Identifying vibrations that destabilize crystals and characterize the glassy state

Greaves, George; Meneau, F.; Majerus, O.; Jones, D.G., (2005) 'Identifying vibrations that destabilize crystals and characterize the glassy state', Science 308(5726) pp.1299-1302 RAE2008

Quantum Monte Carlo in the Interaction Representation --- Application to a Spin-Peierls Model

A quantum Monte Carlo algorithm is constructed starting from the standard perturbation expansion in the interaction representation. The resulting configuration space is strongly related to that of the Stochastic Series Expansion (SSE) method, which is based on a direct power series expansion of exp(-beta*H). Sampling procedures previously developed for the SSE method can therefore be used also in the interaction representation formulation. The new method is first tested on the S=1/2 Heisenberg chain. Then, as an application to a model of great current interest, a Heisenberg chain including phonon degrees of freedom is studied. Einstein phonons are coupled to the spins via a linear modulation of the nearest-neighbor exchange. The simulation algorithm is implemented in the phonon occupation number basis, without Hilbert space truncations, and is exact. Results are presented for the magnetic properties of the system in a wide temperature regime, including the T-->0 limit where the chain undergoes a spin-Peierls transition. Some aspects of the phonon dynamics are also discussed. The results suggest that the effects of dynamic phonons in spin-Peierls compounds such as GeCuO3 and NaV2O5 must be included in order to obtain a correct quantitative description of their magnetic properties, both above and below the dimerization temperature.

Fabrication and characterization of doped and porous silicon nanowires as anodes for lithium ion batteries

By using Si(100) with different dopant type (n++-type (As) or p-type (B)), it is shown how metal-assisted chemically (MAC) etched silicon nanowires (Si NWs) can form with rough outer surfaces around a solid NW core for p-type NWs, and a unique, defined mesoporous structure for highly doped n-type NWs. High resolution electron microscopy techniques were used to define the characteristic roughening and mesoporous structure within the NWs and how such structures can form due to a judicious choice of carrier concentration and dopant type. Control of roughness and internal mesoporosity is demonstrated during the formation of Si NWs from highly doped n-type Si(100) during electroless etching through a systematic investigation of etching parameters (etching time, AgNO3 concentration, %HF and temperature). Raman scattering measurements of the transverse optical phonon confirm quantum size effects and phonon scattering in mesoporous wires associated with the etching condition, including quantum confinement effects for the nanocrystallites of Si comprising the internal structure of the mesoporous NWs. Laser power heating of NWs confirms phonon confinement and scattering from internal mesoporosity causing reduced thermal conductivity. The Li+ insertion and extraction characteristics at n-type and p-type Si(100) electrodes with different carrier density and doping type are investigated by cyclic voltammetry and constant current measurements. The insertion and extraction potentials are demonstrated to vary with cycling and the occurrence of an activation effect is shown in n-type electrodes where the charge capacity and voltammetric currents are found to be much higher than p-type electrodes. X-ray photo-electron spectroscopy (XPS) and Raman scattering demonstrate that highly doped n-type Si(100) retains Li as a silicide and converts to an amorphous phase as a two-step phase conversion process. The findings show the succinct dependence of Li insertion and extraction processes for uniformly doped Si(100) single crystals and how the doping type and its effect on the semiconductor-solution interface dominate Li insertion and extraction, composition, crystallinity changes and charge capacity. The effect of dopant, doping density and porosity of MAC etched Si NWs are investigated. The CV response is shown to change in area (current density) with increasing NW length and in profile shape with a changing porosity of the Si NWs. The CV response also changes with scan rate indicative of a transition from intercalation or alloying reactions, to pseudocapactive charge storage at higher scan rates and for p-type NWs. SEM and TEM show a change in structure of the NWs after Li insertion and extraction due to expansion and contraction of the Si NWs. Galvanostatic measurements show the cycling behavior and the Coulombic efficiency of the Si NWs in comparison to their bulk counterparts.

Modeling non-harmonic behavior of materials from experimental inelastic neutron scattering and thermal expansion measurements

Based on thermodynamic principles, we derive expressions quantifying the non-harmonic vibrational behavior of materials, which are rigorous yet easily evaluated from experimentally available data for the thermal expansion coefficient and the phonon density of states. These experimentally- derived quantities are valuable to benchmark first-principles theoretical predictions of harmonic and non-harmonic thermal behaviors using perturbation theory, ab initio molecular-dynamics, or Monte-Carlo simulations. We illustrate this analysis by computing the harmonic, dilational, and anharmonic contributions to the entropy, internal energy, and free energy of elemental aluminum and the ordered compound FeSi over a wide range of temperature. Results agree well with previous data in the literature and provide an efficient approach to estimate anharmonic effects in materials.

Optical Characterization of Indium Gallium Nitride for Application in High-Efficiency Solar Photovoltaic Cells

The semiconductor alloy indium gallium nitride (InxGa1-xN) offers substantial potential in the development of high-efficiency multi-junction photovoltaic devices due to its wide range of direct band gaps, strong absorption and other optoelectronic properties. This work uses a variety of characterization techniques to examine the properties of InxGa1-xN thin films deposited in a range of compositions by a novel plasma-enhanced evaporation deposition system. Due to the high vapour pressure and low dissociation temperature of indium, the indium incorporation and, ultimately, control of the InxGa1-xN composition was found to be influenced to a greater degree by deposition temperature than variations in the In:Ga source rates in the investigated region of deposition condition space. Under specific deposition conditions, crystalline films were grown in an advantageous nano-columnar microstructure with deposition temperature influencing column size and density. The InxGa1-xN films were determined to have very strong absorption coefficients with band gaps indirectly related to indium content. However, the films also suffer from compositional inhomogeneity and In-related defect complexes with strong phonon coupling that dominates the emission mechanism. This, in addition to the presence of metal impurities, harms the alloy’s electronic properties as no significant photoresponse was observed. This research has demonstrated the material properties that make the InxGa1-xN alloy attractive for multi-junction solar cells and the benefits/drawbacks of the plasma-enhanced evaporation deposition system. Future work is needed to overcome significant challenges relating to crystalline quality, compositional homogeneity and the optoelectronic properties of In-rich InxGa1-xN films in order to develop high-performance photovoltaic devices.

Measurement and interpretation of the mid-infrared properties of single crystal and polycrystalline gold

The contribution of electron-phonon scattering and grain boundary scattering to the mid-IR (lambda = 3.392 mum) properties of An has been assessed by examining both bulk, single crystal samples-Au(1 1 1) and Au(1 1 0)-and thin film, polycrystalline An samples at 300 K and 100 K by means of surface plasmon polariton excitation. The investigation constitutes a stringent test for the in-vacuo Otto-configuration prism coupler used to perform the measurements, illustrating its strengths and limitations. Analysis of the optical response is guided by a physically based interpretation of the Drude model. Relative to the reference case of single crystal Au at 100 K (epsilon = - 568 + i17.5), raising the temperature to 300 K causes increased electron-phonon scattering that accounts for a reduction of similar to40 nm in the electron mean free path. Comparison of a polycrystalline sample to the reference case determines a mean free path due to grain boundary scattering of similar to 17 nm, corresponding to about half the mean grain size as determined from atomic force microscopy and indicating a high reflectance coefficient for the An grain boundaries. An analysis combining consideration of grain boundary scattering and the inclusion of a small percentage of voids in the polycrystalline film by means of an effective medium model indicates a value for the grain boundary reflection coefficient in the range 0.55-0.71. (C) 2005 Elsevier B.V. All rights reserved.

Nonlinear transport and localization in single-walled carbon nanotubes

We have measured the electrical transport properties of mats of single-walled carbon nanotubes (SWNT) as a function of applied electric and magnetic fields. We find that at low temperatures the resistance as a function of temperature R(T) follows the Mott variable range hopping (VRH) formula for hopping in three dimensions. Measurement of the electric field dependence of the resistance R(E) allows for the determination of the Bohr radius of a localized state a = 700nm. The magnetoresistance (MR) of SWNT mat samples is large and negative at all temperatures and fields studied. The low field negative MR is proportional to H2, in agreement with variable range hopping in two or three dimensions. 3D VRH indicates good intertube contacts, implying that the localization is due to the disorder experienced by the individual tubes. The 3D localization radius gives a measure of the ID localization length on the individual tubes, which we estimate to be >700 nm. Implications for the electron-phonon mean free path are discussed.

Correlated electron-ion dynamics: the excitation of atomic motion by energetic electrons

Correlated electron-ion dynamics (CEID) is an extension of molecular dynamics that allows us to introduce in a correct manner the exchange of energy between electrons and ions. The formalism is based on a systematic approximation: small amplitude moment expansion. This formalism is extended here to include the explicit quantum spread of the ions and a generalization of the Hartree-Fock approximation for incoherent sums of Slater determinants. We demonstrate that the resultant dynamical equations reproduce analytically the selection rules for inelastic electron-phonon scattering from perturbation theory, which control the mutually driven excitations of the two interacting subsystems. We then use CEID to make direct numerical simulations of inelastic current-voltage spectroscopy in atomic wires, and to exhibit the crossover from ionic cooling to heating as a function of the relative degree of excitation of the electronic and ionic subsystems.

Inelastic current-voltage spectroscopy of atomic wires

A tight-binding model is developed to describe the electron-phonon coupling in atomic wires under an applied voltage and to model, their inelastic current-voltage spectroscopy. Particular longitudinal phonons are found to have greatly enhanced coupling to the electronic states of the system. This leads to a large drop in differential conductance at threshold energies associated with these phonons. It is found that with increasing tension these energies decrease, while the size of the conductance drops increases, in agreement with experiment.

Relative energetics and structural properties of zirconia using a self-consistent tight-binding model

We describe an empirical, self-consistent, orthogonal tight-binding model for zirconia, which allows for the polarizability of the anions at dipole and quadrupole levels and for crystal field splitting of the cation d orbitals, This is achieved by mixing the orbitals of different symmetry on a site with coupling coefficients driven by the Coulomb potentials up to octapole level. The additional forces on atoms due to the self-consistency and polarizabilities are exactly obtained by straightforward electrostatics, by analogy with the Hellmann-Feynman theorem as applied in first-principles calculations. The model correctly orders the zero temperature energies of all zirconia polymorphs. The Zr-O matrix elements of the Hamiltonian, which measure covalency, make a greater contribution than the polarizability to the energy differences between phases. Results for elastic constants of the cubic and tetragonal phases and phonon frequencies of the cubic phase are also presented and compared with some experimental data and first-principles calculations. We suggest that the model will be useful for studying finite temperature effects by means of molecular dynamics.

Inelastic quantum transport in nanostructures: The self-consistent Born approximation and correlated electron-ion dynamics

A dynamical method for inelastic transport simulations in nanostructures is compared to a steady-state method based on nonequilibrium Green's functions. A simplified form of the dynamical method produces, in the steady state in the weak-coupling limit, effective self-energies analogous to those in the Born approximation due to electron-phonon coupling. The two methods are then compared numerically on a resonant system consisting of a linear trimer weakly embedded between metal electrodes. This system exhibits an enhanced heating at high biases and long phonon equilibration times. Despite the differences in their formulation, the static and dynamical methods capture local current-induced heating and inelastic corrections to the current with good agreement over a wide range of conditions, except in the limit of very high vibrational excitations where differences begin to emerge.

Nonconservative current-induced forces: A physical interpretation

We give a physical interpretation of the recently demonstrated non-conservative nature of interatomic forces in current-carrying nanostructures. We start from the analytical expression for the curl of these forces, and evaluate it for a point defect in a current-carrying system. We obtain a general definition of the capacity of electrical current flow to exert a non-conservative force, and thus do net work around closed paths, by a formal non-invasive test procedure. Second, we show that the gain in atomic kinetic energy in time, generated by non-conservative current-induced forces, is equivalent to the uncompensated stimulated emission of directional phonons. This connection with electron-phonon interactions quantifies explicitly the intuitive notion that non-conservative forces work by angular momentum transfer.

Temperature-dependent gap edge in strong-coupling superconductors determined using the Eliashberg-Nambu formalism

Using the theory of Eliashberg and Nambu for strong-coupling superconductors, we have calculated the gap function for a model superconductor and a selection of real superconductors includong the elements Al, Sn, Tl, Nb, In, Pb and Hg and one alloy, Bi2Tl. We have determined thetemperature-dependent gap edge in each and found that in materials with weak electron-phonon ($\lambda 1.20$), not only is the gap edge double valued but it also departs significantly from the BCS form and develops a shoulderlike structure which may, in some cases, denote a gap edge exceeding the $T = 0$ value. These computational results support the insights obtained by Leavens in an analytic consideration of the general problem. Both the shoulder and double value arise from a common origin seated in the form of the gap function in strong coupled materials at finite temperatures. From the calculated gap function, we can determine the densities of states in the materials and the form of the tunneling current-voltage characteristics for junctions with these materials as electroddes. By way of illustration, results are shown for the contrasting cases of Sn ($\lambda=0.74$) and Hg ($\lambad=1.63$). The reported results are distinct in several ways from BCS predictions and provide an incentive determinative experimental studies with techniques such as tunneling and far infrared absorption.

Extracting a more realistic pseudopotential for aluminum, lead, niobium and tantalum from superconductor electron tunnelling spectroscopy data

Electron tunnelling spectroscopy, developed to extract superconductive metals the electron-phonon spectral density, $\alpha^2F(\nu)$, is found to be a powerful tool also for extracting a more realistic pseudopotential from such metals. The pseudopotential so extracted has a range of surprising but physically reasonable properties and regenerates $\alpha^2F(\nu)$ accurately. Free from most of its long-standing uncertainties, thie pseudopotential may be useful in a number of active fields.

Optical Coalescence and Applications in Cavity Quantum Electrodynamics

We investigate a hitherto largely unexplored regime of cavity quantum electrodynamics in which a highly-reflective element positioned between the end-mirrors of a typical Fabry--P\'erot resonator strongly modifies the cavity response function, such that two longitudinal modes with different spatial parity are brought close to frequency degeneracy. We examine applications of this generic `optical coalescence' phenomenon for the generation of enhanced photon--phonon nonlinearities in optomechanics and atom--photon nonlinearities in cavity quantum electrodynamics with strongly-coupled emitters.