Close packed transitions in slit-shaped pores: Density functional theory study of methane adsorption capacity in carbon

An important feature of improving lattice gas models and classical isotherms is the incorporation of a pore size dependent capacity, which has hitherto been overlooked. In this paper, we develop a model for predicting the temperature dependent variation in capacity with pore size. The model is based on the analysis of a lattice gas model using a density functional theory approach at the close packed limit. Fluid-fluid and solid-fluid interactions are modeled by the Lennard-Jones 12-6 potential and Steele's 10-4-3, potential respectively. The capacity of methane in a slit-shaped carbon pore is calculated from the characteristic parameters of the unit cell, which are extracted by minimizing the grand potential of the unit cell. The capacities predicted by the proposed model are in good agreement with those obtained from grand canonical Monte Carlo simulation, for pores that can accommodate up to three adsorbed layers. Single particle and pair distributions exhibit characteristic features that correspond to the sequence of buckling and rhombic transitions that occur as the slit pore width is increased. The model provides a useful tool to model continuous variation in the microstructure of an adsorbed phase, namely buckling and rhombic transitions, with increasing pore width. (C) 2002 American Institute of Physics.

Simulation of the micro-physics of rocks using LSMearth

The particle-based Lattice Solid Model (LSM) was developed to provide a basis to study the physics of rocks and the nonlinear dynamics of earthquakes (MORA and PLACE, 1994; PLACE and MORA, 1999). A new modular and flexible LSM approach has been developed that allows different microphysics to be easily included in or removed from the model. The approach provides a virtual laboratory where numerical experiments can easily be set up and all measurable quantities visualised. The proposed approach provides a means to simulate complex phenomena such as fracturing or localisation processes, and enables the effect of different micro-physics on macroscopic behaviour to be studied. The initial 2-D model is extended to allow three-dimensional simulations to be performed and particles of different sizes to be specified. Numerical bi-axial compression experiments under different confining pressure are used to calibrate the model. By tuning the different microscopic parameters (such as coefficient of friction, microscopic strength and distribution of grain sizes), the macroscopic strength of the material and can be adjusted to be in agreement with laboratory experiments, and the orientation of fractures is consistent with the theoretical value predicted based on Mohr-Coulomb diagram. Simulations indicate that 3-D numerical models have different macroscopic properties than in 2-D and, hence, the model must be recalibrated for 3-D simulations. These numerical experiments illustrate that the new approach is capable of simulating typical rock fracture behaviour. The new model provides a basis to investigate nucleation, rupture and slip pulse propagation in complex fault zones without the previous model limitations of a regular low-level surface geometry and being restricted to two-dimensions.

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.

Electron exchange coupling for single-donor solid-state spin qubits

Intervalley interference between degenerate conduction band minima has been shown to lead to oscillations in the exchange energy between neighboring phosphorus donor electron states in silicon [B. Koiller, X. Hu, and S. Das Sarma, Phys. Rev. Lett. 88, 027903 (2002); Phys. Rev. B 66, 115201 (2002)]. These same effects lead to an extreme sensitivity of the exchange energy on the relative orientation of the donor atoms, an issue of crucial importance in the construction of silicon-based spin quantum computers. In this article we calculate the donor electron exchange coupling as a function of donor position incorporating the full Bloch structure of the Kohn-Luttinger electron wave functions. It is found that due to the rapidly oscillating nature of the terms they produce, the periodic part of the Bloch functions can be safely ignored in the Heitler-London integrals as was done by Koiller, Hu, and Das Sarma, significantly reducing the complexity of calculations. We address issues of fabrication and calculate the expected exchange coupling between neighboring donors that have been implanted into the silicon substrate using an 15 keV ion beam in the so-called top down fabrication scheme for a Kane solid-state quantum computer. In addition, we calculate the exchange coupling as a function of the voltage bias on control gates used to manipulate the electron wave functions and implement quantum logic operations in the Kane proposal, and find that these gate biases can be used to both increase and decrease the magnitude of the exchange coupling between neighboring donor electrons. The zero-bias results reconfirm those previously obtained by Koiller, Hu, and Das Sarma.

Oxide ionic conductivity and microstructures of Sm- or La-doped CeO2-based systems

The concept of crystallographic index termed the effective index is suggested and applied to the design of ceria (CeO2)-based electrolytes to maximize oxide ionic conductivity. The suggested index considers the fluorite structure, and combines the expected oxygen vacancy level with the ionic radius mismatch between host and dopant cations. Using this approach, oxide ionic conductivity of Sm- or La-doped CeO2-based system has been optimized and tested under operating conditions of a solid oxide fuel cell. In the observation of microstructure in atomic scale, both Sm-doped CeO2 and La-doped CeO2 electrolytes had large micro-domains over 10 nm in the lattice. On the other hand, Sm or La and alkaline earth co-doped CeO2-based electrolytes with high effective index had small micro-domains around 1-3 nm in the microstructure. The large micro-domain would prevent oxide ion from passing through the lattice. Therefore, it is concluded that the improvement of ionic conductivity is reflected in changes of microstructure in atomic scale. (C) 2002 Elsevier Science B.V. All rights reserved.

Numerical study of hydrogenic effective mass theory for an impurity P donor in Si in the presence of an electric field and interfaces

In this paper we examine the effects of varying several experimental parameters in the Kane quantum computer architecture: A-gate voltage, the qubit depth below the silicon oxide barrier, and the back gate depth to explore how these variables affect the electron density of the donor electron. In particular, we calculate the resonance frequency of the donor nuclei as a function of these parameters. To do this we calculated the donor electron wave function variationally using an effective-mass Hamiltonian approach, using a basis of deformed hydrogenic orbitals. This approach was then extended to include the electric-field Hamiltonian and the silicon host geometry. We found that the phosphorous donor electron wave function was very sensitive to all the experimental variables studied in our work, and thus to optimize the operation of these devices it is necessary to control all parameters varied in this paper.

Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

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.

Actively shielded multi-layer gradient coil designs with improved cooling properties

In standard cylindrical gradient coils consisting of a single layer of wires, a limiting factor in achieving very large magnetic field gradients is the rapid increase in coil resistance with efficiency. This is a particular problem in small-bore scanners, such as those used for MR microscopy. By adopting a multi-layer design in which the coil wires are allowed to spread out into multiple layers wound at increasing radii, a more favourable scaling of resistance with efficiency is achieved, thus allowing the design of more powerful gradient coils with acceptable resistance values. Previously this approach has been applied to the design of unshielded, longitudinal, and transverse gradient coils. Here, the multi-layer approach has been extended to allow the design of actively shielded multi-layer gradient coils, and also to produce coils exhibiting enhanced cooling characteristics. An iterative approach to modelling the steady-state temperature distribution within the coil has also been developed. Results indicate that a good level of screening can be achieved in multi-layer coils, that small versions of such coils can yield higher efficiencies at fixed resistance than conventional two-layer (primary and screen) coils, and that performance improves as the number of layers of increases. Simulations show that by optimising multi-layer coils for cooling it is possible to achieve significantly higher gradient strengths at a fixed maximum operating temperature. A four-layer coil of 8 mm inner diameter has been constructed and used to test the steady-state temperature model. (C) 2003 Elsevier Inc. All rights reserved.

Calculation of electric fields induced by body and head motion in high-field MRI

In modern magnetic resonance imaging (MRI), patients are exposed to strong, nonuniform static magnetic fields outside the central imaging region, in which the movement of the body may be able to induce electric currents in tissues which could be possibly harmful. This paper presents theoretical investigations into the spatial distribution of induced electric fields and currents in the patient when moving into the MRI scanner and also for head motion at various positions in the magnet. The numerical calculations are based on an efficient, quasi-static, finite-difference scheme and an anatomically realistic, full-body, male model. 3D field profiles from an actively shielded 4T magnet system are used and the body model projected through the field profile with a range of velocities. The simulation shows that it possible to induce electric fields/currents near the level of physiological significance under some circumstances and provides insight into the spatial characteristics of the induced fields. The results are extrapolated to very high field strengths and tabulated data shows the expected induced currents and fields with both movement velocity and field strength. (C) 2003 Elsevier Science (USA). All rights reserved.

Multilayer transverse gradient coil design

In small, cylindrical gradient coils consisting of a single layer of wires, the limiting factor in achieving large magnetic field gradients is the rapid increase in coil resistance with efficiency. This behavior results from the decrease in the maximum usable wire diameter as the number of turns is increased. By adopting a multilayer design in which the coil wires are allowed to spread out into multiple layers wound at increasing radii, a more favorable scaling of resistance with efficiency is achieved, thus allowing the design of more powerful gradient coils with acceptable resistance values. By extending the theory used to design standard cylindrical gradient coils, mathematical expressions have been developed that allow the design of multilayer coils. These expressions have previously been applied to the design of a four-layer z-gradient coil. As a further development, the equations have now been modified to allow the design of multilayer transverse gradient coils. The variation in coil performance with the number of layers employed has been investigated for coils of a size suitable for use in NMR microscopy, and the effect of constructing the coil using wires or cuts in a continuous conducting surface has also been assessed. We find that at fixed resistance a small wire-wound two-layer coil offers an increase in efficiency of a factor of about 1.5 compared with a single-layer coil. In addition, a two-layer coil of 10-mm inner diameter has been designed and built. This coil had an efficiency of 0.41 Tm-1 A(-1), a resistance of 0.96 +/- 0.01 Omega, and an inductance of 22.3 +/- 0.2 muH. The coil produces a gradient that deviates from linearity by less than 5% over a central cylindrical region of interest of height and length 6.2 mm. (C) 2003 Wiley Periodicals, Inc.

New structure formation on gamma-irradiation of poly(chlorotrifluoroethylene)

The radiation chemistry of PCTFE at different temperatures has been studied. The polymer was irradiated under vacuum to absorbed doses of up to 1500 kGy. Three irradiation temperatures were chosen. These included ambient temperature, a temperature well above the T, and a temperature above the crystalline melting temperature. These were 298, 423 and 493 K, respectively. The formation of new structures was identified by solid-state FTIR and F-19 NMR. No branching was observed below the melting temperature, but branches were observed above the melting temperature. G-values for chain-end formation were 1.5 and 2.4 at room temperature and 423 K, respectively and the G-value for the formation of double bonds was found to be < 0.1. For the irradiations at 493 K, the G-values for the formation of chain ends, double bonds and branching points were 3.6, 0.2 and 0.5, respectively. The presence of long-chain branches within the polymer structure could not be proven for radiolysis at 493 K, but scission predominates and network formation does not occur upon irradiation. DSC studies of the polymers irradiated at ambient temperature were consistent with chain scission leading to an increase in the percentage crystallinity, as observed for other fluoropolymers. (C) 2003 Elsevier Science Ltd. All rights reserved.

Rates of electronic energy transfer in conformationally flexible bichromophoric macrocyclic complexes: a combined experimental and molecular modeling study

Electronic energy transfer (EET) rate constants between a naphthalene donor and anthracene acceptor in [ZnL4a](ClO4)(2) and [ZnL4b](ClO4)(2) were determined by time-resolved fluorescence where L-4a and L-4b are the trans and cis isomers of 6-((anthracen-9-yl-methyl)amino)-6,13-dimethyl-13-((naphthalen-1-yl-methyl)amino)-1,4,8,11-tetraazacyclotetradecane, respectively. These isomers differ in the relative disposition of the appended chromophores with respect to the macrocyclic plane. The trans isomer has an energy transfer rate constant (k(EET)) of 8.7 x 10(8) s(-1), whereas that of the cis isomer is significantly faster (2.3 x 10(9) s(-1)). Molecular modeling was used to determine the likely distribution of conformations in CH3CN solution for these complexes in an attempt to identify any distance or orientation dependency that may account for the differing rate constants observed. The calculated conformational distributions together with analysis by H-1 NMR for the [ZnL4a](2+) trans complex in the common trans-III N-based isomer gave a calculated Forster rate constant close to that observed experimentally. For the [ZnL4b](2+) cis complex, the experimentally determined rate constant may be attributed to a combination of trans-Ill and trans-I N-based isomeric forms of the complex in solution.

Photon-added detection

The production of conditional quantum states and quantum operations based on the result of measurement is now seen as a key tool in quantum information and metrology. We propose a different type of photon number detector. It functions nondeterministically, but when successful, it has high fidelity. The detector, which makes use of an n-photon auxiliary Fock state and high efficiency homodyne detection, allows a tunable trade-off between fidelity and probability. By sacrificing probability of operation, an excellent approximation to a photon-number detector is achieved.