Electromigration of vacancies in copper

The total current-induced force on atoms in a Cu wire containing a vacancy are calculated using the self consistent one-electron density matrix in the presence of an electric current, without separation into electron-wind and direct forces. By integrating the total current-induced force, the change in vacancy migration energy due to the current is calculated. We use the change in migration energy with current to infer an effective electromigration driving force F-e. Finally, we calculate the proportionality constant rho* between F-e and the current density in the wire.

Magnetic tight binding and the iron-chromium enthalpy anomaly

We describe a self-consistent magnetic tight-binding theory based in an expansion of the Hohenberg-Kohn density functional to second order, about a non-spin-polarized reference density. We show how a first order expansion about a density having a trial input magnetic moment leads to a fixed moment model. We employ a simple set of tight-binding parameters that accurately describes electronic structure and energetics, and show these to be transferable between first row transition metals and their alloys. We make a number of calculations of the electronic structure of dilute Cr impurities in Fe, which we compare with results using the local spin density approximation. The fixed moment model provides a powerful means for interpreting complex magnetic configurations in alloys; using this approach, we are able to advance a simple and readily understood explanation for the observed anomaly in the enthalpy of mixing.

Structural changes in quasi-one-dimensional many-electron systems: From linear to zigzag and beyond

Many-electron systems confined to a quasi-one-dimensional geometry by a cylindrical distribution of positive charge have been investigated by density functional computations in the unrestricted local spin density approximation. Our investigations have been focused on the low-density regime, in which electrons are localized. The results reveal a wide variety of different charge and spin configurations, including linear and zig-zag chains, single-and double-strand helices, and twisted chains of dimers. The spin-spin coupling turns from weakly antiferromagnetic at relatively high density, to weakly ferromagnetic at the lowest densities considered in our computations. The stability of linear chains of localized charge has been investigated by analyzing the radial dependence of the self-consistent potential and by computing the dispersion relation of low-energy harmonic excitations.

Unusual decoherence in qubit measurements with a Bose-Einstein condensate

We consider an electrostatic qubit located near a Bose-Einstein condensate (BEC) of noninteracting bosons in a double-well potential, which is used for qubit measurements. Tracing out the BEC variables we obtain a simple analytical expression for the qubit's density matrix. The qubit's evolution exhibits a slow (proportional to 1/root t) damping of the qubit's coherence term, which however turns to be a Gaussian one in the case of static qubit. This is in contrast to the exponential damping produced by most classical detectors. The decoherence is, in general, incomplete and strongly depends on the initial state of the qubit.

Entanglement entropy dynamics of Heisenberg chains

By means of the time dependent density matrix renormalization group algorithm we study the zero-temperature dynamics of the Von Neumann entropy of a block of spins in a Heisenberg chain after a sudden quench in the anisotropy parameter. In the absence of any disorder the block entropy increases linearly with time and then saturates. We analyse the velocity of propagation of the entanglement as a function of the initial and final anisotropies and compare our results, wherever possible, with those obtained by means of conformal field theory. In the disordered case we find a slower ( logarithmic) evolution which may signal the onset of entanglement localization.

The synthesis and electronic structure of a novel [NiS4Fe2(CO)(6)] radical cluster: Implications for the active site of the [NiFe] hydrogenases

A novel [Ni'S-4'Fe-2(CO)(6)] cluster (1: 'S-4'=(CH3C6H3S2)(2)(CH2)(3)) has been synthesised, structurally characterised and has been shown to undergo a chemically reversible reduction process at -1.31 V versus Fc(+)/Fc to generate the EPR-active monoanion 1(-). Multifrequency Q-, X- and S-band EPR spectra of Ni-61-enriched 1(-) show a well-resolved quartet hyperfine splitting in the low-field region due to the interaction with a single Ni-61 (I = 3/2) nucleus. Simulations of the EPR spectra require the introduction of a single angle of non-coincidence between g, and A(1), and g(3) and A(3) to reproduce all of the features in the S- and X-band spectra. This behaviour provides a rare example of the detection and measurement of non-coincidence effects from frozen-solution EPR spectra without the need for single-crystal measurements, and in which the S-band experiment is sensitive to the non-coincidence. An analysis of the EPR spectra of 1(-) reveals a 24% Ni contribution to the SOMO in 1(-), supporting a delocalisation of the spin-density across the NiFe2 cluster. This observation is supported by IR spectroscopic results which show that the CO stretching frequencies, v(CO), shift to lower frequency by about 70 cm(-1) when 1 is reduced to 1(-). Density functional calculations provide a framework for the interpretation of the spectroscopic properties of 1(-) and suggest that the SOMO is delocalised over the whole cluster, but with little S-centre participation. This electronic structure contrasts with that of the Ni-A, -B, -C and -L forms of [NiFe] hydrogenase in which there is considerable S participation in the SOMO.

Angular correlations in radiative cascades following resonant electron capture by highly charged ions

We investigate the angular correlations between the photons emitted in the dielectronic recombination (DR) of initially hydrogenlike heavy ions. The theoretical analysis is performed based on a density-matrix approach and Dirac's relativistic theory. Special emphasis has been placed upon the effects of the higher-order, nondipole terms in the expansion of the electron-photon interaction. To illustrate these effects, we present and discuss detailed calculations for K-LL DR of initially hydrogenlike xenon, gold, and uranium. These computations show that the angular correlations are significantly affected by interference between the leading electric-dipole (E1) and the magnetic-quadrupole (M2) transitions.

Analog of Rabi oscillations in resonant electron-ion systems

Quantum coherence between electron and ion dynamics, observed in organic semiconductors by means of ultrafast spectroscopy, is the object of recent theoretical and computational studies. To simulate this kind of quantum coherent dynamics, we have introduced in a previous article [L. Stella, M. Meister, A. J. Fisher, and A. P. Horsfield, J. Chem. Phys. 127, 214104 (2007)] an improved computational scheme based on Correlated Electron-Ion Dynamics (CEID). In this article, we provide a generalization of that scheme to model several ionic degrees of freedom and many-body electronic states. To illustrate the capability of this extended CEID, we study a model system which displays the electron-ion analog of the Rabi oscillations. Finally, we discuss convergence and scaling properties of the extended CEID along with its applicability to more realistic problems. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3589165]

Full characterization of the quantum linear-zigzag transition in atomic chains

A string of repulsively interacting particles exhibits a phase transition to a zigzag structure, by reducing the transverse trap potential or the interparticle distance. Based on the emergent symmetry Z2 it has been argued that this instability is a quantum phase transition, which can be mapped to an Ising model in transverse field. An extensive Density Matrix Renormalization Group analysis is performed, resulting in an high-precision evaluation of the critical exponents and of the central charge of the system, confirming that the quantum linear-zigzag transition belongs to the critical Ising model universality class. Quantum corrections to the classical phase diagram are computed, and the range of experimental parameters where quantum effects play a role is provided. These results show that structural instabilities of one-dimensional interacting atomic arrays can simulate quantum critical phenomena typical of ferromagnetic systems.

Universal tight binding model for chemical reactions in solution and at surfaces. I. Organic molecules

As is now well established, a first order expansion of the Hohenberg-Kohn total energy density functional about a trial input density, namely, the Harris-Foulkes functional, can be used to rationalize a non self consistent tight binding model. If the expansion is taken to second order then the energy and electron density matrix need to be calculated self consistently and from this functional one can derive a charge self consistent tight binding theory. In this paper we have used this to describe a polarizable ion tight binding model which has the benefit of treating charge transfer in point multipoles. This admits a ready description of ionic polarizability and crystal field splitting. It is necessary in constructing such a model to find a number of parameters that mimic their more exact counterparts in the density functional theory. We describe in detail how this is done using a combination of intuition, exact analytical fitting, and a genetic optimization algorithm. Having obtained model parameters we show that this constitutes a transferable scheme that can be applied rather universally to small and medium sized organic molecules. We have shown that the model gives a good account of static structural and dynamic vibrational properties of a library of molecules, and finally we demonstrate the model's capability by showing a real time simulation of an enolization reaction in aqueous solution. In two subsequent papers, we show that the model is a great deal more general in that it will describe solvents and solid substrates and that therefore we have created a self consistent quantum mechanical scheme that may be applied to simulations in heterogeneous catalysis.

Ultra-fast spin dynamics in the 2s2p^2 configuration of C^+ studied by time-dependent R-matrix theory

We employ time-dependent R-matrix theory to study ultra-fast dynamics in the doublet 2s2p(2) configuration of C+ for a total magnetic quantum number M = 1. In contrast to the dynamics observed for M = 0, ultra-fast dynamics for M = 1 is governed by spin dynamics in which the 2s electron acts as a flag rather than a spectator electron. Under the assumption that m(S) = 1/2, m(2s) = 1/2 allows spin dynamics involving the two 2p electrons, whereas m(2s) = -1/2 prevents spin dynamics of the two 2p electrons. For a pump-probe pulse scheme with (h) over bar omega(pump) = 10.9 eV and (h) over bar omega(probe) = 16.3 eV and both pulses six cycles long, little sign of spin dynamics is observed in the total ionization probability. Signs of spin dynamics can be observed, however, in the ejected-electron momentum distributions. We demonstrate that the ejected-electron momentum distributions can be used for unaligned targets to separate the contributions of initial M = 0 and M = 1 levels. This would, in principle, allow unaligned target ions to be used to obtain information on the different dynamics in the 2s2p(2) configuration for the M = 0 and M = 1 levels from a single experime

Spin and spatial dynamics in electron-impact scattering off S-wave He using R-matrix with Time-Dependence theory

R-matrix with time-dependence theory is applied to electron-impact ionisation processes for He in the S-wave model. Cross sections for electron-impact excitation, ionisation and ionisation with excitation for impact energies between 25 and 225 eV are in excellent agreement with benchmark cross sections. Ultra-fast dynamics induced by a scattering event is observed through time-dependent signatures associated with autoionisation from doubly excited states. Further insight into dynamics can be obtained through examination of the spin components of the time-dependent wavefunction.

Tools for analysing configuration interaction wavefunctions

The configuration interaction (CI) approach to quantum chemical calculations is a well-established means of calculating accurately the solution to the Schrodinger equation for many-electron systems. It represents the many-body electron wavefunction as a sum of spin-projected Slater determinants of orthogonal one-body spin-orbitals. The CI wavefunction becomes the exact solution of the Schrodinger equation as the length of the expansion becomes infinite, however, it is a difficult quantity to visualise and analyse for many-electron problems. We describe a method for efficiently calculating the spin-averaged one- and two-body reduced density matrices rho(psi)((r) over bar; (r) over bar' ) and Gamma(psi)((r) over bar (1), (r) over bar (2); (r) over bar'(1), (r) over bar'(2)) of an arbitrary CI wavefunction Psi. These low-dimensional functions are helpful tools for analysing many-body wavefunctions; we illustrate this for the case of the electron-electron cusp. From rho and Gamma one can calculate the matrix elements of any one- or two-body spin-free operator (O) over cap. For example, if (O) over cap is an applied electric field, this field can be included into the CI Hamiltonian and polarisation or gating effects may be studied for finite electron systems. (C) 2003 Elsevier B.V. All rights reserved.

Current-driven magnetic rearrangements in spin-polarized point contacts

A method for investigating the dynamics of atomic magnetic moments in current-carrying magnetic point contacts under bias is presented. This combines the nonequilibrium Green's function (NEGF) method for evaluating the current and the charge density with a description of the dynamics of the magnetization in terms of quasistatic thermally activated transitions between stationary configurations. This method is then implemented in a tight-binding (TB) model with parameters chosen to simulate the main features of the electronic structures of magnetic transition metals. We investigate the domain wall (DW) migration in magnetic monoatomic chains sandwiched between magnetic leads, and for realistic parameters find that collinear arrangement of the magnetic moments of the chain is always favorable. Several stationary magnetic configurations are identified, corresponding to a different number of Bloch walls in the chain and to a different current. The relative stability of these configurations depends on the geometrical details of the junction and on the bias; however, we predict transitions between different configurations with activation barriers of the order of a few tens of meV. Since different magnetic configurations are associated with different resistances, this suggests an intrinsic random telegraph noise at microwave frequencies in the I-V curves of magnetic atomic point contacts at room temperature. Finally, we investigate whether or not current-induced torques are conservative.

Eigenvalues of the time-delay matrix in overlapping resonances

We discuss the properties of the lifetime or the time-delay matrix Q(E) for multichannel scattering, which is related to the scattering matrix S(E) by Q = i?S(dS†/dE). For two overlapping resonances occurring at energies E with widths G(? = 1, 2), with an energy-independent background, only two eigenvalues of Q(E) are proved to be different from zero and to show typical avoided-crossing behaviour. These eigenvalues are expressible in terms of the four resonance parameters (E , G) and a parameter representing the strength of the interaction of the resonances. An example of the strong and weak interaction in an overlapping double resonance is presented for the positronium negative ion. When more than two resonances overlap (? = 1, ..., N), no simple representation of each eigenvalue has been found. However, the formula for the trace of the Q-matrix leads to the expression d(E) = -?arctan[(G/2)/(E - E)] + d(E) for the eigenphase sum d(E) and the background eigenphase sum d(E), in agreement with the known form of the state density. The formulae presented in this paper are useful in a parameter fitting of overlapping resonances. © 2006 IOP Publishing Ltd.