Phase diagram of a Heisenberg spin-Peierls model with quantum phonons

Using a new version of the density-matrix renormalization group we determine the phase diagram of a model of an antiferromagnetic Heisenberg spin chain where the spins interact with quantum phonons. A quantum phase transition from a gapless spin-fluid state to a gapped dimerized phase occurs at a nonzero value of the spin-phonon coupling. The transition is in the same universality class as that of a frustrated spin chain, to which the model maps in the diabatic limit. We argue that realistic modeling of known spin-Peierls materials should include the effects of quantum phonons.

Solution of a two-leg spin ladder system

A model for a spin-1/2 ladder system with two legs is introduced. It is demonstrated that this model is solvable via the Bethe ansatz method for arbitrary values of the rung coupling J. This is achieved by a suitable mapping from the Hubbard model with appropriate twisted boundary conditions. We determine that a phase transition between gapped and gapless spin excitations occurs at the critical value J(c) = 1/2 of the rung coupling.

Integrable impurity spin ladder systems

Two different types of integrable impurities in a spin ladder system are proposed. The impurities are introduced in such a way that the integrability of the models is not violated. The models are solved exactly with the Bethe ansatz equations as well as the energy eigenvalues obtained. We show for both models that a phase transition between gapped and gapless spin excitations occurs at a critical value of the rung coupling J. In addition, the dependence of the impurities on this phase transition is determined explicitly. In one of the models the spin gap decreases by increasing the impurity strength A. Moreover, for a fixed A, a reduction in the spin gap by increasing the impurity concentration is also observed.

Large-N solutions of the Heisenberg and Hubbard-Heisenberg models on the anisotropic triangular lattice: application to Cs2CuCl4 and to the layered organic superconductors kappa-(BEDT-TTF)(2)X (BEDT-TTF bis(ethylene-dithio)tetrathiafulvalene); X anion)

We solve the Sp(N) Heisenberg and SU(N) Hubbard-Heisenberg models on the anisotropic triangular lattice in the large-N limit. These two models may describe respectively the magnetic and electronic properties of the family of layered organic materials K-(BEDT-TTF)(2)X, The Heisenberg model is also relevant to the frustrated antiferromagnet, Cs2CuCl4. We find rich phase diagrams for each model. The Sp(N) :antiferromagnet is shown to have five different phases as a function of the size of the spin and the degree of anisotropy of the triangular lattice. The effects of fluctuations at finite N are also discussed. For parameters relevant to Cs2CuCl4 the ground state either exhibits incommensurate spin order, or is in a quantum disordered phase with deconfined spin-1/2 excitations and topological order. The SU(N) Hubbard-Heisenberg model exhibits an insulating dimer phase, an insulating box phase, a semi-metallic staggered flux phase (SFP), and a metallic uniform phase. The uniform and SFP phases exhibit a pseudogap, A metal-insulator transition occurs at intermediate values of the interaction strength.

Spin squeezing as a measure of entanglement in a two-qubit system

We show that the two definitions of spin squeezing extensively used in the literature [M. Kitagawa and M. Ueda, Phys. Rev. A 47, 5138 (1993) and D.J. Wineland , Phys. Rev. A 50, 67 (1994)] give different predictions of entanglement in the two-atom Dicke system. We analyze differences between the definitions and show that the spin squeezing parameter of Kitagawa and Ueda is a better measure of entanglement than the commonly used spectroscopic spin squeezing parameter. We illustrate this relation by examining different examples of a driven two-atom Dicke system in which spin squeezing and entanglement arise dynamically. We give an explanation of the source of the difference using the negativity criterion for entanglement.

Spin-detection in a quantum electromechanical shuttle system

We study the electrical transport of a harmonically bound, single-molecule C-60 shuttle operating in the Coulomb blockade regime, i.e. single electron shuttling. In particular, we examine the dependance of the tunnel current on an ultra-small stationary force exerted on the shuttle. As an example, we consider the force exerted on an endohedral N@C-60 by the magnetic field gradient generated by a nearby nanomagnet. We derive a Hamiltonian for the full shuttle system which includes the metallic contacts, the spatially dependent tunnel couplings to the shuttle, the electronic and motional degrees of freedom of the shuttle itself and a coupling of the shuttle's motion to a phonon bath. We analyse the resulting quantum master equation and find that, due to the exponential dependence of the tunnel probability on the shuttle-contact separation, the current is highly sensitive to very small forces. In particular, we predict that the spin state of the endohedral electrons of N@C-60 in a large magnetic gradient field can be distinguished from the resulting current signals within a few tens of nanoseconds. This effect could prove useful for the detection of the endohedral spin-state of individual paramagnetic molecules such as N@C-60 and P@C-60, or the detection of very small static forces acting on a C-60 shuttle.

Entanglement between a qubit and the environment in the spin-boson model

The quantitative description of the quantum entanglement between a qubit and its environment is considered. Specifically, for the ground state of the spin-boson model, the entropy of entanglement of the spin is calculated as a function of α, the strength of the ohmic coupling to the environment, and ɛ, the level asymmetry. This is done by a numerical renormalization group treatment of the related anisotropic Kondo model. For ɛ=0, the entanglement increases monotonically with α, until it becomes maximal for α→1-. For fixed ɛ>0, the entanglement is a maximum as a function of α for a value, α=αM

Matrix elements of U(2n) generators in a multishell spin-orbit basis. I. The del-operator MEs in a two-shell composite Gelfand-Paldus basis

This is the first in a series of three articles which aimed to derive the matrix elements of the U(2n) generators in a multishell spin-orbit basis. This is a basis appropriate to many-electron systems which have a natural partitioning of the orbital space and where also spin-dependent terms are included in the Hamiltonian. The method is based on a new spin-dependent unitary group approach to the many-electron correlation problem due to Gould and Paldus [M. D. Gould and J. Paldus, J. Chem. Phys. 92, 7394, (1990)]. In this approach, the matrix elements of the U(2n) generators in the U(n) x U(2)-adapted electronic Gelfand basis are determined by the matrix elements of a single Ll(n) adjoint tensor operator called the del-operator, denoted by Delta(j)(i) (1 less than or equal to i, j less than or equal to n). Delta or del is a polynomial of degree two in the U(n) matrix E = [E-j(i)]. The approach of Gould and Paldus is based on the transformation properties of the U(2n) generators as an adjoint tensor operator of U(n) x U(2) and application of the Wigner-Eckart theorem. Hence, to generalize this approach, we need to obtain formulas for the complete set of adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. The nonzero shift coefficients are uniquely determined and may he evaluated by the methods of Gould et al. [see the above reference]. In this article, we define zero-shift adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis which are appropriate to the many-electron problem. By definition, these are proportional to the corresponding two-shell del-operator matrix elements, and it is shown that the Racah factorization lemma applies. Formulas for these coefficients are then obtained by application of the Racah factorization lemma. The zero-shift adjoint reduced Wigner coefficients required for this procedure are evaluated first. All these coefficients are needed later for the multishell case, which leads directly to the two-shell del-operator matrix elements. Finally, we discuss an application to charge and spin densities in a two-shell molecular system. (C) 1998 John Wiley & Sons.

Matrix elements of U(2n) generators in a multishell spin-orbit basis. II. The two-shell nonzero shift ACCs and U(2n) generator MEs

This is the second in a series of articles whose ultimate goal is the evaluation of the matrix elements (MEs) of the U(2n) generators in a multishell spin-orbit basis. This extends the existing unitary group approach to spin-dependent configuration interaction (CI) and many-body perturbation theory calculations on molecules to systems where there is a natural partitioning of the electronic orbital space. As a necessary preliminary to obtaining the U(2n) generator MEs in a multishell spin-orbit basis, we must obtain a complete set of adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. The zero-shift coefficients were obtained in the first article of the series. in this article, we evaluate the nonzero shift adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. We then demonstrate that the one-shell versions of these coefficients may be obtained by taking the Gelfand-Tsetlin limit of the two-shell formulas. These coefficients,together with the zero-shift types, then enable us to write down formulas for the U(2n) generator matrix elements in a two-shell spin-orbit basis. Ultimately, the results of the series may be used to determine the many-electron density matrices for a partitioned system. (C) 1998 John Wiley & Sons, Inc.

Effect of disorder on quantum phase transitions in anisotropic XY spin chains in a transverse field

We present some exact results for the effect of disorder on the critical properties of an anisotropic XY spin chain in a transverse held. The continuum limit of the corresponding fermion model is taken and in various cases results in a Dirac equation with a random mass. Exact analytic techniques can then be used to evaluate the density of states and the localization length. In the presence of disorder the ferromagnetic-paramagnetic or Ising transition of the model is in the same universality class as the random transverse field Ising model solved by Fisher using a real-space renormalization-group decimation technique (RSRGDT). If there is only randomness in the anisotropy of the magnetic exchange then the anisotropy transition (from a ferromagnet in the x direction to a ferromagnet in the y direction) is also in this universality class. However, if there is randomness in the isotropic part of the exchange or in the transverse held then in a nonzero transverse field the anisotropy transition is destroyed by the disorder. We show that in the Griffiths' phase near the Ising transition that the ground-state energy has an essential singularity. The results obtained for the dynamical critical exponent, typical correlation length, and for the temperature dependence of the specific heat near the Ising transition agree with the results of the RSRODT and numerical work. [S0163-1829(99)07125-8].

Single-spin measurement using single-electron transistors to probe two-electron systems

We present a method for measuring single spins embedded in a solid by probing two-electron systems with a single-electron transistor (SET). Restrictions imposed by the Pauli principle on allowed two-electron states mean that the spin state of such systems has a profound impact on the orbital states (positions) of the electrons, a parameter which SET's are extremely well suited to measure. We focus on a particular system capable of being fabricated with current technology: a Te double donor in Si adjacent to a Si/SiO2, interface and lying directly beneath the SET island electrode, and we outline a measurement strategy capable of resolving single-electron and nuclear spins in this system. We discuss the limitations of the measurement imposed by spin scattering arising from fluctuations emanating from the SET and from lattice phonons. We conclude that measurement of single spins, a necessary requirement for several proposed quantum computer architectures, is feasible in Si using this strategy.

Exactly solvable quantum spin ladders associated with the orthogonal and symplectic Lie algebras

We extend the results of spin ladder models associated with the Lie algebras su(2(n)) to the case of the orthogonal and symplectic algebras o(2(n)), sp(2(n)) where n is the number of legs for the system. Two classes of models are found whose symmetry, either orthogonal or symplectic, has an explicit n dependence. Integrability of these models is shown for an arbitrary coupling of XX-type rung interactions and applied magnetic field term.

Quantum spin ladder systems associated with su(2 vertical bar 2)

Two integrable quantum spin ladder systems will be introduced associated with the fundamental su(2 \2) solution of the Yang-Baxter equation. The first model is a generalized quantum Ising system with Ising rung interactions. In the second model the addition of extra interactions allows us to impose Heisenberg rung interactions without violating integrability. The existence of a Bethe ansatz solution for both models allows us to investigate the elementary excitations for antiferromagnetic rung couplings. We find that the first model does not show a gap whilst in the second case there is a gap for all positive values of the rung coupling.

Molecular orbital calculations of two-electron states for P-donor solid-state spin qubits

We theoretically study the Hilbert space structure of two neighboring P-donor electrons in silicon-based quantum computer architectures. To use electron spins as qubits, a crucial condition is the isolation of the electron spins from their environment, including the electronic orbital degrees of freedom. We provide detailed electronic structure calculations of both the single donor electron wave function and the two-electron pair wave function. We adopted a molecular orbital method for the two-electron problem, forming a basis with the calculated single donor electron orbitals. Our two-electron basis contains many singlet and triplet orbital excited states, in addition to the two simple ground state singlet and triplet orbitals usually used in the Heitler-London approximation to describe the two-electron donor pair wave function. We determined the excitation spectrum of the two-donor system, and study its dependence on strain, lattice position, and interdonor separation. This allows us to determine how isolated the ground state singlet and triplet orbitals are from the rest of the excited state Hilbert space. In addition to calculating the energy spectrum, we are also able to evaluate the exchange coupling between the two donor electrons, and the double occupancy probability that both electrons will reside on the same P donor. These two quantities are very important for logical operations in solid-state quantum computing devices, as a large exchange coupling achieves faster gating times, while the magnitude of the double occupancy probability can affect the error rate.

Integrable anisotropic spin-ladder model

We present an integrable spin-ladder model, which possesses a free parameter besides the rung coupling J. Wang's system based on the SU(4) symmetry can be obtained as a special case. The model is exactly solvable by means of the Bethe ansatz method. We determine the dependence on the anisotropy parameter of the phase transition between gapped and gapless spin excitations and present the phase diagram. Finally, we show that the model is a special case of a more general Hamiltonian with three free parameters.