Full S matrix calculation via a single real-symmetric Lanczos recursion: The Lanczos artificial boundary inhomogeneity method

We present an efficient and robust method for the calculation of all S matrix elements (elastic, inelastic, and reactive) over an arbitrary energy range from a single real-symmetric Lanczos recursion. Our new method transforms the fundamental equations associated with Light's artificial boundary inhomogeneity approach [J. Chem. Phys. 102, 3262 (1995)] from the primary representation (original grid or basis representation of the Hamiltonian or its function) into a single tridiagonal Lanczos representation, thereby affording an iterative version of the original algorithm with greatly superior scaling properties. The method has important advantages over existing iterative quantum dynamical scattering methods: (a) the numerically intensive matrix propagation proceeds with real symmetric algebra, which is inherently more stable than its complex symmetric counterpart; (b) no complex absorbing potential or real damping operator is required, saving much of the exterior grid space which is commonly needed to support these operators and also removing the associated parameter dependence. Test calculations are presented for the collinear H+H-2 reaction, revealing excellent performance characteristics. (C) 2004 American Institute of Physics.

Fungible dynamics: There are only two types of entangling multiple-qubit interactions

What interactions are sufficient to simulate arbitrary quantum dynamics in a composite quantum system? It has been shown that all two-body Hamiltonian evolutions can be simulated using any fixed two-body entangling n-qubit Hamiltonian and fast local unitaries. By entangling we mean that every qubit is coupled to every other qubit, if not directly, then indirectly via intermediate qubits. We extend this study to the case where interactions may involve more than two qubits at a time. We find necessary and sufficient conditions for an arbitrary n-qubit Hamiltonian to be dynamically universal, that is, able to simulate any other Hamiltonian acting on n qubits, possibly in an inefficient manner. We prove that an entangling Hamiltonian is dynamically universal if and only if it contains at least one coupling term involving an even number of interacting qubits. For odd entangling Hamiltonians, i.e., Hamiltonians with couplings that involve only an odd number of qubits, we prove that dynamic universality is possible on an encoded set of n-1 logical qubits. We further prove that an odd entangling Hamiltonian can simulate any other odd Hamiltonian and classify the algebras that such Hamiltonians generate. Thus, our results show that up to local unitary operations, there are only two fundamentally different types of entangling Hamiltonian on n qubits. We also demonstrate that, provided the number of qubits directly coupled by the Hamiltonian is bounded above by a constant, our techniques can be made efficient.

Energy as an entanglement witness for quantum many-body systems

We investigate quantum many-body systems where all low-energy states are entangled. As a tool for quantifying such systems, we introduce the concept of the entanglement gap, which is the difference in energy between the ground-state energy and the minimum energy that a separable (unentangled) state may attain. If the energy of the system lies within the entanglement gap, the state of the system is guaranteed to be entangled. We find Hamiltonians that have the largest possible entanglement gap; for a system consisting of two interacting spin-1/2 subsystems, the Heisenberg antiferromagnet is one such example. We also introduce a related concept, the entanglement-gap temperature: the temperature below which the thermal state is certainly entangled, as witnessed by its energy. We give an example of a bipartite Hamiltonian with an arbitrarily high entanglement-gap temperature for fixed total energy range. For bipartite spin lattices we prove a theorem demonstrating that the entanglement gap necessarily decreases as the coordination number is increased. We investigate frustrated lattices and quantum phase transitions as physical phenomena that affect the entanglement gap.

Fast simulation of a quantum phase transition in an ion-trap realizable unitary map

We show that a specific implementation of a unitary map on multiple qubits in an ion trap is physically equivalent to a Hamiltonian evolution that belongs to the same universality class as the transverse Ising Hamiltonian. We suggest experimental signatures, and present numerical simulations for the case of four qubits.

Quantum dynamics of a model for two Josephson-coupled Bose-Einstein condensates

In this work, we investigate the quantum dynamics of a model for two singlemode Bose-Einstein condensates which are coupled via Josephson tunnelling. Using direct numerical diagonalization of the Hamiltonian, we compute the time evolution of the expectation value for the relative particle number across a wide range of couplings. Our analysis shows that the system exhibits rich and complex behaviours varying between harmonic and non-harmonic oscillations, particularly around the threshold coupling between the delocalized and selftrapping phases. We show that these behaviours are dependent on both the initial state of the system and regime of the coupling. In addition, a study of the dynamics for the variance of the relative particle number expectation and the entanglement for different initial states is presented in detail.

On the chaotic rotation of a liquid-filled gyrostat via the Melnikov-Holmes-Marsden integral

The non-linear motions of a gyrostat with an axisymmetrical, fluid-filled cavity are investigated. The cavity is considered to be completely filled with an ideal incompressible liquid performing uniform rotational motion. Helmholtz theorem, Euler's angular momentum theorem and Poisson equations are used to develop the disturbed Hamiltonian equations of the motions of the liquid-filled gyrostat subjected to small perturbing moments. The equations are established in terms of a set of canonical variables comprised of Euler angles and the conjugate angular momenta in order to facilitate the application of the Melnikov-Holmes-Marsden (MHM) method to investigate homoclinic/heteroclinic transversal intersections. In such a way, a criterion for the onset of chaotic oscillations is formulated for liquid-filled gyrostats with ellipsoidal and torus-shaped cavities and the results are confirmed via numerical simulations. (c) 2006 Elsevier Ltd. All rights reserved.

Analysis of chaotic instabilities in a rotating body with internal energy dissipation

Melnikov's method is used to analytically predict the onset of chaotic instability in a rotating body with internal energy dissipation. The model has been found to exhibit chaotic instability when a harmonic disturbance torque is applied to the system for a range of forcing amplitude and frequency. Such a model may be considered to be representative of the dynamical behavior of a number of physical systems such as a spinning spacecraft. In spacecraft, disturbance torques may arise under malfunction of the control system, from an unbalanced rotor, from vibrations in appendages or from orbital variations. Chaotic instabilities arising from such disturbances could introduce uncertainties and irregularities into the motion of the multibody system and consequently could have disastrous effects on its intended operation. A comprehensive stability analysis is performed and regions of nonlinear behavior are identified. Subsequently, the closed form analytical solution for the unperturbed system is obtained in order to identify homoclinic orbits. Melnikov's method is then applied on the system once transformed into Hamiltonian form. The resulting analytical criterion for the onset of chaotic instability is obtained in terms of critical system parameters. The sufficient criterion is shown to be a useful predictor of the phenomenon via comparisons with numerical results. Finally, for the purposes of providing a complete, self-contained investigation of this fundamental system, the control of chaotic instability is demonstated using Lyapunov's method.

On the chaotic instability of a nonsliding liquid-filled top with a small spheroidal base via Melnikov-Holmes-Marsden integrals

Chaotic orientations of a top containing a fluid filled cavity are investigated analytically and numerically under small perturbations. The top spins and rolls in nonsliding contact with a rough horizontal plane and the fluid in the ellipsoidal shaped cavity is considered to be ideal and describable by finite degrees of freedom. A Hamiltonian structure is established to facilitate the application of Melnikov-Holmes-Marsden (MHM) integrals. In particular, chaotic motion of the liquid-filled top is identified to be arisen from the transversal intersections between the stable and unstable manifolds of an approximated, disturbed flow of the liquid-filled top via the MHM integrals. The developed analytical criteria are crosschecked with numerical simulations via the 4th Runge-Kutta algorithms with adaptive time steps.

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.

T-Q relation and exact solution for the XYZ chain with general non-diagonal boundary terms

We propose that the Baxter's Q-operator for the quantum XYZ spin chain with open boundary conditions is given by the j -> infinity limit of the corresponding transfer matrix with spin-j (i.e., (2j + I)-dimensional) auxiliary space. The associated T-Q relation is derived from the fusion hierarchy of the model. We use this relation to determine the Bethe Ansatz solution of the eigenvalues of the fundamental transfer matrix. The solution yields the complete spectrum of the Hamiltonian. (c) 2006 Elsevier B.V. All rights reserved.

A Bethe ansatz solvable model for superpositions of Cooper pairs and condensed molecular bosons

We introduce a general Hamiltonian describing coherent superpositions of Cooper pairs and condensed molecular bosons. For particular choices of the coupling parameters, the model is integrable. One integrable manifold, as well as the Bethe ansatz solution, was found by Dukelsky et al. [J. Dukelsky, G.G. Dussel, C. Esebbag, S. Pittel, Phys. Rev. Lett. 93 (2004) 050403]. Here we show that there is a second integrable manifold, established using the boundary quantum inverse scattering method. In this manner we obtain the exact solution by means of the algebraic Bethe ansatz. In the case where the Cooper pair energies are degenerate we examine the relationship between the spectrum of these integrable Hamiltonians and the quasi-exactly solvable spectrum of particular Schrodinger operators. For the solution we derive here the potential of the Schrodinger operator is given in terms of hyperbolic functions. For the solution derived by Dukelsky et al., loc. cit. the potential is sextic and the wavefunctions obey PT-symmetric boundary conditions. This latter case provides a novel example of an integrable Hermitian Hamiltonian acting on a Fock space whose states map into a Hilbert space of PE-symmetric wavefunctions defined on a contour in the complex plane. (c) 2006 Elsevier B.V. All rights reserved.

Teleportation via Multi-Qubit Channels

We investigate the problem of teleporting an unknown qubit state to a recipient via a channel of 2L qubits. In this procedure a protocol is employed whereby L Bell state measurements are made and information based on these measurements is sent via a classical channel to the recipient. Upon receiving this information the recipient determines a local gate which is used to recover the original state. We find that the 2(2L)-dimensional Hilbert space of states available for the channel admits a decomposition into four subspaces. Every state within a given subspace is a perfect channel, and each sequence of Bell measurements projects 2L qubits of the system into one of the four subspaces. As a result, only two bits of classical information need be sent to the recipient for them to determine the gate. We note some connections between these four subspaces and ground states of many-body Hamiltonian systems, and discuss the implications of these results towards understanding entanglement in multi-qubit systems.

XSophe, a computer simulation software suite for the analysis of electron paramagnetic resonance spectra

The XSophe computer simulation software suite consisting of a daemon, the XSophe interface and the computational program Sophe is a state of the art package for the simulation of electron paramagnetic resonance spectra. The Sophe program performs the computer simulation and includes a number of new technologies including; the SOPHE partition and interpolation schemes, a field segmentation algorithm, homotopy, parallelisation and spectral optimisation. The SOPHE partition and interpolation scheme along with a field segmentation algorithm greatly increases the speed of simulations for most systems. Multidimensional homotopy provides an efficient method for accurately tracing energy levels and hence tracing transitions in the presence of energy level anticrossings and looping transitions and allowing computer simulations in frequency space. Recent enhancements to Sophe include the generalised treatment of distributions of orientational parameters, termed the mosaic misorientation linewidth model and a faster more efficient algorithm for the calculation of resonant field positions and transition probabilities. For complex systems the parallelisation enables the simulation of these systems on a parallel computer and the optimisation algorithms in the suite provide the experimentalist with the possibility of finding the spin Hamiltonian parameters in a systematic manner rather than a trial-and-error process. The XSophe software suite has been used to simulate multifrequency EPR spectra (200 MHz to 6 00 GHz) from isolated spin systems (S > ~½) and coupled centres (Si, Sj _> I/2). Griffin, M.; Muys, A.; Noble, C.; Wang, D.; Eldershaw, C.; Gates, K.E.; Burrage, K.; Hanson, G.R."XSophe, a Computer Simulation Software Suite for the Analysis of Electron Paramagnetic Resonance Spectra", 1999, Mol. Phys. Rep., 26, 60-84.

Dimethylsulfide dehydrogenase from rhodovulum sulfidophilum: EPR spectroscopy and biochemical analysis reveal its place in the DMSO reductase family of molybdenum enzymes

Dimethylsulfide (DMS) dehydrogenase catalyses the oxidation of DMS to dimethylsulfoxide. The purified enzyme has three subunits of Mr = 94, 38 and 32 kDa and has an optical spectrum dominated by a b-type cytochrome. The metal ion and nucleotide analysis revealed 0.5 g-atom Mo, 9.8 g-atom Fe and 1.96 mol GMP per tool of enzyme. Taken together, these data indicate that DMS dehydrogenase contains a bis(MGD)Mo cofactor. A comparison of the Nterminal amino acid sequence of DMS dehydrogenase revealed that the Mo-containing ct-subunit was most closely related to the c~-subunits of nitrate reductase (NarG) and selenate reductase (SerA). Similarly, the [~-subunit of DMS dehydrogenase was most closely related to the [3-subunits of nitrate reductase (NarH) and selenate reductase (SerB). Variable temperature X-band EPR spectra (120-2K) of 'as isolated' DMS dehydrogenase showed resonances arising from multiple redox centres, Mo(V), [3Fe-4S] +, [4Fe-4S] ÷. A pH dependent EPR study of the Mo(V) centre in lH20 and 2H20 reveals the presence of three Mo(V) species in equilibrium, Mo(V)-OH2, Mo(V)-X and Mo(V)-OH. Between pH6 and 8.2 the dominant species is Mo(V)-OH2 and Mo(V)-X is a minor component. X is probably the anion, chloride. Comparison of the rhombicity and anisotropy parameters for the Mo(V) species in DMS dehydrogenase with other Mo(V) centres in metalloproteins showed that it was most similar to the low pH nitrite spectrum of E. coli nitrate reductase (NarGHI). The spin Hamiltonian parameters (2.0158, 1.8870, 1.8620) for the [4Fe-4S] + cluster suggests the presence of histidine (N) coordination to iron in this cluster. It is suggested that this unusual [Fe-S] cluster may be associated with a histidine-cysteine rich sequence at the N-terminus of the ct-subunit of DMS dehydrogenase.