Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

Quantum technology: the second quantum revolution

We are currently in the midst of a second quantum revolution. The first quantum revolution gave us new rules that govern physical reality. The second quantum revolution will take these rules and use them to develop new technologies. In this review we discuss the principles upon which quantum technology is based and the tools required to develop it. We discuss a number of examples of research programs that could deliver quantum technologies in coming decades including: quantum information technology, quantum electromechanical systems, coherent quantum electronics, quantum optics and coherent matter technology.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

Theoretical and numerical analyses of convective instability in porous media with temperature-dependent viscosity

Exact analytical solutions of the critical Rayleigh numbers have been obtained for a hydrothermal system consisting of a horizontal porous layer with temperature-dependent viscosity. The boundary conditions considered are constant temperature and zero vertical Darcy velocity at both the top and bottom of the layer. Not only can the derived analytical solutions be readily used to examine the effect of the temperature-dependent viscosity on the temperature-gradient driven convective flow, but also they can be used to validate the numerical methods such as the finite-element method and finite-difference method for dealing with the same kind of problem. The related analytical and numerical results demonstrated that the temperature-dependent viscosity destabilizes the temperature-gradient driven convective flow and therefore, may affect the ore body formation and mineralization in the upper crust of the Earth. Copyright (C) 2003 John Wiley Sons, Ltd.

Measuring the decoherence rate in a semiconductor charge qubit

We describe a method by which the decoherence time of a solid-state qubit may be measured. The qubit is coded in the orbital degree of freedom of a single electron bound to a pair of donor impurities in a semiconductor host. The qubit is manipulated by adiabatically varying an external electric field. We show that by measuring the total probability of a successful qubit rotation as a function of the control field parameters, the decoherence rate may be determined. We estimate various system parameters, including the decoherence rates due to electromagnetic fluctuations and acoustic phonons. We find that, for reasonable physical parameters, the experiment is possible with existing technology. In particular, the use of adiabatic control fields implies that the experiment can be performed with control electronics with a time resolution of tens of nanoseconds.

Quantum process tomography of a controlled-NOT gate

We demonstrate complete characterization of a two-qubit entangling process-a linear optics controlled-NOT gate operating with coincident detection-by quantum process tomography. We use a maximum-likelihood estimation to convert the experimental data into a physical process matrix. The process matrix allows an accurate prediction of the operation of the gate for arbitrary input states and a calculation of gate performance measures such as the average gate fidelity, average purity, and entangling capability of our gate, which are 0.90, 0.83, and 0.73, respectively.

Utilizing encoding in scalable linear optics quantum computing

We present a scheme which offers a significant reduction in the resources required to implement linear optics quantum computing. The scheme is a variation of the proposal of Knill, Laflamme and Milburn, and makes use of an incremental approach to the error encoding to boost probability of success.

Signature of mott-insulator transition with ultracold fermions in a one-dimensional optical lattice

Using the exact Bethe ansatz solution of the Hubbard model and Luttinger liquid theory, we investigate the density profiles and collective modes of one-dimensional ultracold fermions confined in an optical lattice with a harmonic trapping potential. We determine a generic phase diagram in terms of a characteristic filling factor and a dimensionless coupling constant. The collective oscillations of the atomic mass density, a technique that is commonly used in experiments, provide a signature of the quantum phase transition from the metallic phase to the Mott-insulator phase. A detailed experimental implementation is proposed.

Finite-temperature correlations and density profiles of an inhomogeneous interacting one-dimensional Bose gas

We calculate the density profiles and density correlation functions of the one-dimensional Bose gas in a harmonic trap, using the exact finite-temperature solutions for the uniform case, and applying a local density approximation. The results are valid for a trapping potential that is slowly varying relative to a correlation length. They allow a direct experimental test of the transition from the weak-coupling Gross-Pitaevskii regime to the strong-coupling, fermionic Tonks-Girardeau regime. We also calculate the average two-particle correlation which characterizes the bulk properties of the sample, and find that it can be well approximated by the value of the local pair correlation in the trap center.

Double diffusion-driven convective instability of three-dimensional fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the double diffusion-driven convective instability of three-dimensional fluid-saturated geological fault zones when they are heated uniformly from below. The fault zone is assumed to be more permeable than its surrounding rocks. In particular, we have derived exact analytical solutions to the total critical Rayleigh numbers of the double diffusion-driven convective flow. Using the corresponding total critical Rayleigh numbers, the double diffusion-driven convective instability of a fluid-saturated three-dimensional geological fault zone system has been investigated. The related theoretical analysis demonstrates that: (1) The relative higher concentration of the chemical species at the top of the three-dimensional geological fault zone system can destabilize the convective flow of the system, while the relative lower concentration of the chemical species at the top of the three-dimensional geological fault zone system can stabilize the convective flow of the system. (2) The double diffusion-driven convective flow modes of the three-dimensional geological fault zone system are very close each other and therefore, the system may have the similar chance to pick up different double diffusion-driven convective flow modes, especially in the case of the fault thickness to height ratio approaching 0. (3) The significant influence of the chemical species diffusion on the convective instability of the three-dimensional geological fault zone system implies that the seawater intrusion into the surface of the Earth is a potential mechanism to trigger the convective flow in the shallow three-dimensional geological fault zone system.

Mechanical loading modulates glutamate receptor subunit expression in bone

The cellular mechanisms coupling mechanical loading with bone remodeling remain unclear. In the CNS, the excitatory amino acid glutamate (Glu) serves as a potent neurotransmitter exerting its effects via various membrane Glu receptors (GluR). Nerves containing Glu exist close to bone cells expressing functional GluRs. Demonstration of a mechanically sensitive glutamate/aspartate transporter protein and the ability of glutamate to stimulate bone resorption in vitro suggest a role for glutamate linking mechanical load and bone remodeling. We used immunohistochemical techniques to identify the expression of N-methyl-D-aspartate acid (NMDA) and non-NMDA (AMPA or kainate) ionotropic GluR subunits on bone cells in vivo. In bone sections from young adult rats, osteoclasts expressed numerous GluR subunits including AMPA (GluR2/3 and GluR4), kainic acid (GluR567) and NMDA (NMDAR2A, NMDAR2B and NMDAR2C) receptor subtypes. Bone lining cells demonstrated immunoexpression for NMDAR2A, NMDAR2B, NMDAR2C, GluR567, GluR23, GuR2 and GluR4 subunits. Immunoexpression was not evident on osteocytes, chondrocytes or vascular channels. To investigate the effects of mechanical loading on GluR expression, we used a Materials Testing System (MTS) to apply 10 N sinusoidal axial compressive loads percutaneously to the right limbs (radius/ulna, tibia/fibula) of rats. Each limb underwent 300-load cycles/day (cycle rate, 1 Hz) for 4 consecutive days. Contralateral, non-loaded limbs served as controls. Mechanically loaded limbs revealed a load-induced loss of immunoexpression for GluR2/3, GluR4, GluR567 and NMDAR2A on osteoclasts and NMDAR2A, NMDAR2B, GluR2/3 and GluR4 on bone lining cells. Both neonatal rabbit and rat osteoclasts were cultured on bone slices to investigate the effect of the NMDA receptor antagonist, MK801, and the AMPA/kainic acid receptor antagonist, NBQX, on osteoclast resorptive activity in vitro. The inhibition of resorptive function seen suggested that both NMDAR and kainic acid receptor function are required for normal osteoclast function. While the exact role of ionotropic GluRs in skeletal tissue remains unclear, the modulation of GluR subunit expression by mechanical loading lends further support for participation of Glu as a mechanical loading effector. These ionotropic receptors appear to be functionally relevant to normal osteoclast resorptive activity. (C) 2005 Elsevier Inc. All rights reserved.

Periodic electric field enhanced transport through membranes

We examine the mean flux across a homogeneous membrane of a charged tracer subject to an alternating, symmetric voltage waveform. The analysis is based on the Nernst-Planck flux equation, with electric field subject to time dependence only. For low frequency electric fields the quasi steady-state flux can be approximated using the Goldman model, which has exact analytical solutions for tracer concentration and flux. No such closed form solutions can be found for arbitrary frequencies, however we find approximations for high frequency. An approximation formula for the average flux at all frequencies is also obtained from the two limiting approximations. Numerical integration of the governing equation is accomplished by use of the numerical method of lines and is performed for four different voltage waveforms. For the different voltage profiles, comparisons are made with the approximate analytical solutions which demonstrates their applicability. (c) 2005 Elsevier B.V. All rights reserved.

The Qu-Prolog unification algorithm: Formalisation and correctness

Qu-Prolog is an extension of Prolog which performs meta-level computations over object languages, such as predicate calculi and lambda-calculi, which have object-level variables, and quantifier or binding symbols creating local scopes for those variables. As in Prolog, the instantiable (meta-level) variables of Qu-Prolog range over object-level terms, and in addition other Qu-Prolog syntax denotes the various components of the object-level syntax, including object-level variables. Further, the meta-level operation of substitution into object-level terms is directly represented by appropriate Qu-Prolog syntax. Again as in Prolog, the driving mechanism in Qu-Prolog computation is a form of unification, but this is substantially more complex than for Prolog because of Qu-Prolog's greater generality, and especially because substitution operations are evaluated during unification. In this paper, the Qu-Prolog unification algorithm is specified, formalised and proved correct. Further, the analysis of the algorithm is carried out in a frame-work which straightforwardly allows the 'completeness' of the algorithm to be proved: though fully explicit answers to unification problems are not always provided, no information is lost in the unification process.

Paramagnetism of caesium titanium alum

Caesium titanium alum, CsTi(SO4)(2) . 12H(2)O, is a beta alum and exhibits a large trigonal field and a dynamic Jahn-Teller effect. Exact calculations of the linear (2)T(2)xe Jahn-Teller coupling show that in the strict S-6 Site symmetry the ground multiplet consists of a Kramers doublet 2 Gamma(6) with magnetic splitting factors g(parallel to)=1.1 and g perpendicular to=0, a Gamma(4) Gamma(5) doublet at similar to 60 cm(-1) with g(parallel to)=2.51 and g(perpendicular to)=0.06 and another Gamma(4) Gamma(5) doublet at similar to 270 cm(-1) with g(parallel to)=1.67 and g(perpendicular to)=1.83. The controversial g values observed below 4.2 K, g(parallel to)=1.25 and g(perpendicular to)=1.14, are shown to arise from low symmetry distortions. These distortions couple the vibronic levels and induce into the ground state the off-diagonal axial Zeeman interaction that exists between the first excited and the ground vibronic levels. (C) 1997 American Institute of Physics.

Quantum Lie algebras, their existence, uniqueness and q-antisymmetry

Quantum Lie algebras are generalizations of Lie algebras which have the quantum parameter h built into their structure. They have been defined concretely as certain submodules L-h(g) of the quantized enveloping algebras U-h(g). On them the quantum Lie product is given by the quantum adjoint action. Here we define for any finite-dimensional simple complex Lie algebra g an abstract quantum Lie algebra g(h) independent of any concrete realization. Its h-dependent structure constants are given in terms of inverse quantum Clebsch-Gordan coefficients. We then show that all concrete quantum Lie algebras L-h(g) are isomorphic to an abstract quantum Lie algebra g(h). In this way we prove two important properties of quantum Lie algebras: 1) all quantum Lie algebras L-h(g) associated to the same g are isomorphic, 2) the quantum Lie product of any Ch(B) is q-antisymmetric. We also describe a construction of L-h(g) which establishes their existence.