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.

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.

Stokes-operator-squeezed continuous-variable polarization states

We investigate nonclassical Stokes-operator variances in continuous-wave polarization-squeezed laser light generated from one and two optical parametric amplifiers. A general expression of how Stokes-operator variances decompose into two-mode quadrature operator variances is given. Stokes parameter variance spectra for four different polarization-squeezed states have been measured and compared with a coherent state. Our measurement results are visualized by three-dimensional Stokes-operator noise volumes mapped on the quantum Poincare sphere. We quantitatively compare the channel capacity of the different continuous-variable polarization states for communication protocols. It is shown that squeezed polarization states provide 33% higher channel capacities than the optimum coherent beam protocol.

Experimental investigation of continuous-variable quantum teleportation

We report the experimental demonstration of quantum teleportation of the quadrature amplitudes of a light field. Our experiment was stably locked for long periods, and was analyzed in terms of fidelity F and with signal transfer T-q=T++T- and noise correlation V-q=Vinparallel to out+Vinparallel to out-. We observed an optimum fidelity of 0.64+/-0.02, T-q=1.06+/-0.02, and V-q=0.96+/-0.10. We discuss the significance of both T-q>1 and V-q

Quantum theory of the far-off-resonance continuous-wave Raman laser: Heisenberg-Langevin approach

We present the quantum theory of the far-off-resonance continuous-wave Raman laser using the Heisenberg-Langevin approach. We show that the simplified quantum Langevin equations for this system are mathematically identical to those of the nondegenerate optical parametric oscillator in the time domain with the following associations: pump pump, Stokes signal, and Raman coherence idler. We derive analytical results for both the steady-state behavior and the time-dependent noise spectra, using standard linearization procedures. In the semiclassical limit, these results match with previous purely semiclassical treatments, which yield excellent agreement with experimental observations. The analytical time-dependent results predict perfect photon statistics conversion from the pump to the Stokes and nonclassical behavior under certain operational conditions.

A study of the radiation chemistry of poly(chlorotrifluoroethylene) by ESR spectroscopy

The ESR spectra of poly(chlorotrifluoroethylene) were recorded following gamma-radiolysis under vacuum at room temperature and 77 K. The very broad spectrum at 77 K revealed little fine structure with which to identity the radicals formed upon irradiation, but subsequent photobleaching and annealing studies, together with radiolytic studies at higher temperatures, afforded scope for making radical assignments. Both main-chain radicals and a range of chain-end radicals have been identified. The G-values for radical formation were 1.55, 0.36 and 0.32 at 77 K, 273 K and room temperature, respectively. (C) 2003 Elsevier Science Ltd. All rights reserved.

Tailoring teleportation to the quantum alphabet

We introduce a refinement of the standard continuous variable teleportation measurement and displacement strategies. This refinement makes use of prior knowledge about the target state and the partial information carried by the classical channel when entanglement is nonmaximal. This gives an improvement in the output quality of the protocol. The strategies we introduce could be used in current continuous variable teleportation experiments.

Microcanonical temperature for a classical field: Application to Bose-Einstein condensation

We show that the projected Gross-Pitaevskii equation (PGPE) can be mapped exactly onto Hamilton's equations of motion for classical position and momentum variables. Making use of this mapping, we adapt techniques developed in statistical mechanics to calculate the temperature and chemical potential of a classical Bose field in the microcanonical ensemble. We apply the method to simulations of the PGPE, which can be used to represent the highly occupied modes of Bose condensed gases at finite temperature. The method is rigorous, valid beyond the realms of perturbation theory, and agrees with an earlier method of temperature measurement for the same system. Using this method we show that the critical temperature for condensation in a homogeneous Bose gas on a lattice with a uv cutoff increases with the interaction strength. We discuss how to determine the temperature shift for the Bose gas in the continuum limit using this type of calculation, and obtain a result in agreement with more sophisticated Monte Carlo simulations. We also consider the behavior of the specific heat.

Quantum dynamics as a physical resource

How useful is a quantum dynamical operation for quantum information processing? Motivated by this question, we investigate several strength measures quantifying the resources intrinsic to a quantum operation. We develop a general theory of such strength measures, based on axiomatic considerations independent of state-based resources. The power of this theory is demonstrated with applications to quantum communication complexity, quantum computational complexity, and entanglement generation by unitary operations.