Universal simulation of Hamiltonian dynamics for quantum systems with finite-dimensional state spaces

What interactions are sufficient to simulate arbitrary quantum dynamics in a composite quantum system? Dodd [Phys. Rev. A 65, 040301(R) (2002)] provided a partial solution to this problem in the form of an efficient algorithm to simulate any desired two-body Hamiltonian evolution using any fixed two-body entangling N-qubit Hamiltonian, and local unitaries. We extend this result to the case where the component systems are qudits, that is, have D dimensions. As a consequence we explain how universal quantum computation can be performed with any fixed two-body entangling N-qudit Hamiltonian, and local unitaries.

The crystallography and morphology of Mo2C in ferrite

The orientation relationships between hexagonal Mo2C precipitates (H) in ferrite (B) have been determined by electron diffraction to an accuracy of +/-2degrees. With one exception, the 19 results are consistent with the previously reported Pitsch and Schrader (P/S) orientation relationship. However, these more accurate determinations show clearly that there is a systematic deviation of up to 5.5degrees from the exact P/S relationship and that this deviation consists of a small rotation about the parallel close packed directions-[100](B)//[2 (1) over bar(1) over bar0](H). The long direction of the Mo2C needles has been determined unequivocally in terms of the orientation relationship to be [100](B)//[2 (1) over bar(1) over bar0](H). Moire fringes between precipitate and matrix have been used to improve the accuracy of the orientation relationship results and to determine the lattice parameters of the carbide precipitates investigated. The Moire fringe analysis has shown small systematic departures from the exact parallelism between [100](B) and [2 (1) over bar(1) over bar0](H) along the length of Mo2C needles and a lowering of the carbide lattice parameter that is consistent with the replacement of Mo by Fe in the carbide. The orientation relationship results, including the observed systematic deviation from the exact P/S relationship, are shown to be consistent with the edge-to-edge model. (C) 2002 Kluwer Academic Publishers.

Preparing encoded states in an oscillator

Recently a scheme has been proposed for constructing quantum error-correcting codes that embed a finite-dimensional code space in the infinite-dimensional Hilbert space of a system described by continuous quantum variables. One of the difficult steps in this scheme is the preparation of the encoded states. We show how these states can be generated by coupling a continuous quantum variable to a single qubit. An ion trap quantum computer provides a natural setting for a continuous system coupled to a qubit. We discuss how encoded states may be generated in an ion trap.

Evidence from mitochondrial DNA that head lice and body lice of humans (Phthiraptera : Pediculidae) are conspecific

The specific status of the head and body lice of humans has been debated for more than 200 yr. To clarify the specific status of head and body lice, we sequenced 524 base pairs (bp) of the cytochrome oxidase I (COI) gene of 28 head and 28 body lice from nine countries. Ten haplotypes that differed by 1-5 bp at II nucleotide positions were identified. A phylogeny of these sequences indicates that these head and body lice are not from reciprocally monophyletic lineages. Indeed, head and body lice share three of the 10 haplotypes we found. F-ST values and exact tests of haplotype frequencies showed significant differences between head and body lice. However, the same tests also showed significant differences among lice from different countries. Indeed, more of the variation in haplotype frequencies was explained by differences among lice from different countries than by differences between head and body lice. Our results indicate the following: (1) bead and body lice do not represent reciprocally monophyletic lineages and are conspecific; (2) gene flow among populations of lice from different countries is limited; and (3) frequencies of COI haplotypes can be used to study maternal gene flow among populations of head and body lice and thus transmission of lice among their human hosts.

QCD thermal phase transition in the presence of a small chemical potential

We propose a new method to investigate the thermal properties of QCD with a small quark chemical potential mu. Derivatives of quark and gluonic observables with respect to mu are computed at mu=0 for two flavors of p4 improved staggered fermions with ma=0.1,0.2 on a 16(3)x4 lattice, and used to calculate the leading order Taylor expansion in mu of the location of the pseudocritical point about mu=0. This expansion should be well behaved for the small values of mu(q)/T(c)similar to0.1 relevant for BNL RHIC phenomenology, and predicts a critical curve T-c(mu) in reasonable agreement with estimates obtained using exact reweighting. In addition, we contrast the case of isoscalar and isovector chemical potentials, quantify the effect of munot equal0 on the equation of state, and comment on the complex phase of the fermion determinant in QCD with munot equal0.

Confined-unconfined flow in a horizontal confined aquifer

Mixed confined and unconfined groundwater flow occurs in a bounded initially dry aquifer when the hydraulic head at the side boundary suddenly rises above the elevation of the aquifer's top boundary. The flow problem as modelled by the Boussinesq equation is non-trivial because of the involvement of two moving boundaries. The transformed equation (based on a similarity transformation) can, however, be dealt with more easily. Here, we present an approximate analytical solution for this flow problem. The approximate solution is compared with an 'exact' numerical solution and found to be a very accurate description for describing the mixed confined and unconfined flow in the confined aquifer. (C) 2002 Elsevier Science B.V. All rights reserved.

A design method for superconducting MRI magnets with ferromagnetic material

The emphasis of this work is on the optimal design of MRI magnets with both superconducting coils and ferromagnetic rings. The work is directed to the automated design of MRI magnet systems containing superconducting wire and both `cold' and `warm' iron. Details of the optimization procedure are given and the results show combined superconducting and iron material MRI magnets with excellent field characteristics. Strong, homogeneous central magnetic fields are produced with little stray or external field leakage. The field calculations are performed using a semi-analytical method for both current coil and iron material sources. Design examples for symmetric, open and asymmetric clinical MRI magnets containing both superconducting coils and ferromagnetic material are presented.

Rapid field calculations for the effect of ferromagnetic material in MRI magnet design

Magnetic resonance imaging (MRI) magnets have very stringent constraints on the homogeneity of the static magnetic field that they generate over desired imaging regions. The magnet system also preferably generates very little stray field external to its structure, so that ease of siting and safety are assured. This work concentrates on deriving, means of rapidly computing the effect of 'cold' and 'warm' ferromagnetic material in or around the superconducting magnet system, so as to facilitate the automated design of hybrid material MR magnets. A complete scheme for the direct calculation of the spherical harmonics of the magnetic field generated by a circular ring of ferromagnetic material is derived under the conditions of arbitrary external magnetizing fields. The magnetic field produced by the superconducting coils in the system is computed using previously developed methods. The final, hybrid algorithm is fast enough for use in large-scale optimization methods. The resultant fields from a practical example of a 4 T, clinical MRI magnet containing both superconducting coils and magnetic material are presented.

Diverse synaptic mechanisms generate direction selectivity in the rabbit retina

The synaptic conductance of the On-Off direction-selective ganglion cells was measured during visual stimulation to determine whether the direction selectivity is a property of the circuitry presynaptic to the ganglion cells or is generated by postsynaptic interaction of excitatory and inhibitory inputs. Three synaptic asymmetries were identified that contribute to the generation of direction-selective responses: (1) a presynaptic mechanism producing stronger excitation in the preferred direction, (2) a presynaptic mechanism producing stronger inhibition in the opposite direction, and (3) postsynaptic interaction of excitation with spatially offset inhibition. Although the on- and off-responses showed the same directional tuning, the off-response was generated by all three mechanisms, whereas the on- response was generated primarily by the two presynaptic mechanisms. The results indicate that, within a single neuron, different strategies are used within distinct dendritic arbors to accomplish the same neural computation.

The use of stochastic dynamic programming in optimal landscape reconstruction for metapopulations

A decision theory framework can be a powerful technique to derive optimal management decisions for endangered species. We built a spatially realistic stochastic metapopulation model for the Mount Lofty Ranges Southern Emu-wren (Stipiturus malachurus intermedius), a critically endangered Australian bird. Using diserete-time Markov,chains to describe the dynamics of a metapopulation and stochastic dynamic programming (SDP) to find optimal solutions, we evaluated the following different management decisions: enlarging existing patches, linking patches via corridors, and creating a new patch. This is the first application of SDP to optimal landscape reconstruction and one of the few times that landscape reconstruction dynamics have been integrated with population dynamics. SDP is a powerful tool that has advantages over standard Monte Carlo simulation methods because it can give the exact optimal strategy for every landscape configuration (combination of patch areas and presence of corridors) and pattern of metapopulation occupancy, as well as a trajectory of strategies. It is useful when a sequence of management actions can be performed over a given time horizon, as is the case for many endangered species recovery programs, where only fixed amounts of resources are available in each time step. However, it is generally limited by computational constraints to rather small networks of patches. The model shows that optimal metapopulation, management decisions depend greatly on the current state of the metapopulation,. and there is no strategy that is universally the best. The extinction probability over 30 yr for the optimal state-dependent management actions is 50-80% better than no management, whereas the best fixed state-independent sets of strategies are only 30% better than no management. This highlights the advantages of using a decision theory tool to investigate conservation strategies for metapopulations. It is clear from these results that the sequence of management actions is critical, and this can only be effectively derived from stochastic dynamic programming. The model illustrates the underlying difficulty in determining simple rules of thumb for the sequence of management actions for a metapopulation. This use of a decision theory framework extends the capacity of population viability analysis (PVA) to manage threatened species.

One-dimensional quantum walk with unitary noise

The effect of unitary noise on the discrete one-dimensional quantum walk is studied using computer simulations. For the noiseless quantum walk, starting at the origin (n=0) at time t=0, the position distribution P-t(n) at time t is very different from the Gaussian distribution obtained for the classical random walk. Furthermore, its standard deviation, sigma(t) scales as sigma(t)similar tot, unlike the classical random walk for which sigma(t)similar toroott. It is shown that when the quantum walk is exposed to unitary noise, it exhibits a crossover from quantum behavior for short times to classical-like behavior for long times. The crossover time is found to be Tsimilar toalpha(-2), where alpha is the standard deviation of the noise.

Discrete element modelling of power draw of tumbling mills

The power required to operate large mills is typically 5-10 MW. Hence, optimisation of power consumption will have a significant impact on overall economic performance and environmental impact. Power draw modelling results using the discrete element code PFC3D have been compared with results derived from the widely used empirical Model of Morrell. This is achieved by calculating the power draw for a range of operating conditions for constant mill size and fill factor using two modelling approaches. fThe discrete element modelling results show that, apart from density, selection of the appropriate material damping ratio is critical for the accuracy of modelling of the mill power draw. The relative insensitivity of the power draw to the material stiffness allows selection of moderate stiffness values, which result in acceptable computation time. The results obtained confirm that modelling of the power draw for a vertical slice of the mill, of thickness 20% of the mill length, is a reliable substitute for modelling the full mill. The power draw predictions from PFC3D show good agreement with those obtained using the empirical model. Due to its inherent flexibility, power draw modelling using PFC3D appears to be a viable and attractive alternative to empirical models where necessary code and computer power are available.

Modeling diffraction in free-space optical interconnects by the mode expansion method

Free-space optical interconnects (FSOIs), made up of dense arrays of vertical-cavity surface-emitting lasers, photodetectors and microlenses can be used for implementing high-speed and high-density communication links, and hence replace the inferior electrical interconnects. A major concern in the design of FSOIs is minimization of the optical channel cross talk arising from laser beam diffraction. In this article we introduce modifications to the mode expansion method of Tanaka et al. [IEEE Trans. Microwave Theory Tech. MTT-20, 749 (1972)] to make it an efficient tool for modelling and design of FSOIs in the presence of diffraction. We demonstrate that our modified mode expansion method has accuracy similar to the exact solution of the Huygens-Kirchhoff diffraction integral in cases of both weak and strong beam clipping, and that it is much more accurate than the existing approximations. The strength of the method is twofold: first, it is applicable in the region of pronounced diffraction (strong beam clipping) where all other approximations fail and, second, unlike the exact-solution method, it can be efficiently used for modelling diffraction on multiple apertures. These features make the mode expansion method useful for design and optimization of free-space architectures containing multiple optical elements inclusive of optical interconnects and optical clock distribution systems. (C) 2003 Optical Society of America.

Calculation of bound and resonance states of HO2 for nonzero total angular momentum

Bound and resonance states of HO2 have been calculated quantum mechanically by the Lanczos homogeneous filter diagonalization method [Zhang and Smith, Phys. Chem. Chem. Phys. 3, 2282 (2001); J. Chem. Phys. 115, 5751 (2001)] for nonzero total angular momentum J = 1,2,3. For lower bound states, agreement between the results in this paper and previous work is quite satisfactory; while for high lying bound states and resonances these are the first reported results. A helicity quantum number V assignment (within the helicity conserving approximation) is performed and the results indicate that for lower bound states it is possible to assign the V quantum numbers unambiguously, but for resonances it is impossible to assign the V helicity quantum numbers due to strong mixing. In fact, for the high-lying bound states, the mixing has already appeared. These results indicate that the helicity conserving approximation is not good for the resonance state calculations and exact quantum calculations are needed to accurately describe the reaction dynamics for HO2 system. Analysis of the resonance widths shows that most of the resonances are overlapping and the interferences between them lead to large fluctuations from one resonance to another. In accord with the conclusions from earlier J = 0 calculations, this indicates that the dissociation of HO2 is essentially irregular. (C) 2003 American Institute of Physics.

Learning the dynamics of embedded clauses

Recent work by Siegelmann has shown that the computational power of recurrent neural networks matches that of Turing Machines. One important implication is that complex language classes (infinite languages with embedded clauses) can be represented in neural networks. Proofs are based on a fractal encoding of states to simulate the memory and operations of stacks. In the present work, it is shown that similar stack-like dynamics can be learned in recurrent neural networks from simple sequence prediction tasks. Two main types of network solutions are found and described qualitatively as dynamical systems: damped oscillation and entangled spiraling around fixed points. The potential and limitations of each solution type are established in terms of generalization on two different context-free languages. Both solution types constitute novel stack implementations - generally in line with Siegelmann's theoretical work - which supply insights into how embedded structures of languages can be handled in analog hardware.