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.

Utilization of the three-dimensional volcano surface to understand the chemistry of multiphase systems in heterogeneous catalysis

CO hydrogenation is used as a model system to understand why multiphase catalysts are chemically important in heterogeneous catalysis. By including both adsorption and subsequent surface reactions, kinetic equations are derived with two fundamental properties, the chemisorption energies of C and O (Delta H-C and Delta H-O, respectively). By plotting the activity against Delta H-C and Delta H-O, a 3-D volcano surface is obtained. Because of the constraint between Delta H-C and Delta H-O on monophase systems, a maximum can be achieved. However, if multiphase systems are used, such a constraint can be released and the global maximum may be achieved.

Mixed-dimensional coupling in finite element models

In many finite element analysis models it would be desirable to combine reduced- or lower-dimensional element types with higher-dimensional element types in a single model. In order to achieve compatibility of displacements and stress equilibrium at the junction or interface between the differing element types, it is important in such cases to integrate into the analysis some scheme for coupling the element types. A novel and effective scheme for establishing compatibility and equilibrium at the dimensional interface is described and its merits and capabilities are demonstrated. Copyright (C) 2000 John Wiley & Sons, Ltd.

Decoherence-based exploration of d-dimensional one-way quantum computation: Information transfer and basic gates

We study the effects of amplitude and phase damping decoherence in d-dimensional one-way quantum computation. We focus our attention on low dimensions and elementary unidimensional cluster state resources. Our investigation shows how information transfer and entangling gate simulations are affected for d >= 2. To understand motivations for extending the one-way model to higher dimensions, we describe how basic qudit cluster states deteriorate under environmental noise of experimental interest. In order to protect quantum information from the environment, we consider encoding logical qubits into qudits and compare entangled pairs of linear qubit-cluster states to single qudit clusters of equal length and total dimension. A significant reduction in the performance of cluster state resources for d > 2 is found when Markovian-type decoherence models are present.

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.

Boson pairs in a one-dimensional split trap

We describe the properties of a pair of ultracold bosonic atoms in a one-dimensional harmonic trapping potential with a tunable zero-ranged barrier at the trap center. The full characterization of the ground state is done by calculating the reduced single-particle density, the momentum distribution, and the two-particle entanglement. We derive several analytical expressions in the limit of infinite repulsion (Tonks-Girardeau limit) and extend the treatment to finite interparticle interactions by numerical solution. As pair interactions in double wells form a fundamental building block for many-body systems in periodic potentials, our results have implications for a wide range of problems.

Comparing the implementation of two-dimensional numerical quadrature on GPU, FPGA and ClearSpeed systems to study electron scattering by atoms

The use of accelerators, with compute architectures different and distinct from the CPU, has become a new research frontier in high-performance computing over the past ?ve years. This paper is a case study on how the instruction-level parallelism offered by three accelerator technologies, FPGA, GPU and ClearSpeed, can be exploited in atomic physics. The algorithm studied is the evaluation of two electron integrals, using direct numerical quadrature, a task that arises in the study of intermediate energy electron scattering by hydrogen atoms. The results of our ‘productivity’ study show that while each accelerator is viable, there are considerable differences in the implementation strategies that must be followed on each.

Entangling two distant oscillators with a quantum reservoir

The generation of entanglement between two oscillators that interact via a common reservoir is theoretically studied. The reservoir is modeled by a one-dimensional harmonic crystal initially in thermal equilibrium. Starting from a separable state, the oscillators can become entangled after a transient time, that is of the order of the thermalization time scale. This behaviour is observed at finite temperature even when the oscillators are at a distance significantly larger than the crystal's interparticle spacing. The underlying physical mechanisms can be explained by the dynamical properties of the collective variables of the two oscillators which may decouple from or be squeezed by the reservoir. Our predictions can be tested with an ion chain in a linear Paul trap. Copyright (C) EPLA, 2011

Scaling of the entanglement spectrum near quantum phase transitions

The entanglement spectrum describing quantum correlations in many-body systems has been recently recognized as a key tool to characterize different quantum phases, including topological ones. Here we derive its analytically scaling properties in the vicinity of some integrable quantum phase transitions and extend our studies also to nonintegrable quantum phase transitions in one-dimensional spin models numerically. Our analysis shows that, in all studied cases, the scaling of the difference between the two largest nondegenerate Schmidt eigenvalues yields with good accuracy critical points and mass scaling exponents.