Complex numbers and symmetries in quantum mechanics, and a nonlinear superposition principle for Wigner functions

Complex numbers appear in the Hilbert space formulation of quantum mechanics, but not in the formulation in phase space. Quantum symmetries are described by complex, unitary or antiunitary operators defining ray representations in Hilbert space, whereas in phase space they are described by real, true representations. Equivalence of the formulations requires that the former representations can be obtained from the latter and vice versa. Examples are given. Equivalence of the two formulations also requires that complex superpositions of state vectors can be described in the phase space formulation, and it is shown that this leads to a nonlinear superposition principle for orthogonal, pure-state Wigner functions. It is concluded that the use of complex numbers in quantum mechanics can be regarded as a computational device to simplify calculations, as in all other applications of mathematics to physical phenomena.

First-principles quantum dynamics in interacting Bose gases II: stochastic gauges

First principles simulations of the quantum dynamics of interacting Bose gases using the stochastic gauge representation are analysed. In a companion paper, we showed how the positive-P representation can be applied to these problems using stochastic differential equations. That method, however, is limited by increased sampling error as time evolves. Here, we show how the sampling error can be greatly reduced and the simulation time significantly extended using stochastic gauges. In particular, local stochastic gauges (a subset) are investigated. Improvements are confirmed in numerical calculations of single-, double- and multi-mode systems in the weak-mode coupling regime. Convergence issues are investigated, including the recognition of two modes by which stochastic equations produced by phase-space methods in general can diverge: movable singularities and a noise-weight relationship. The example calculated here displays wave-like behaviour in spatial correlation functions propagating in a uniform 1D gas after a sudden change in the coupling constant. This could in principle be tested experimentally using Feshbach resonance methods.

First-principles quantum dynamics in interacting Bose gases: I. The positive P representation

The performance of the positive P phase-space representation for exact many- body quantum dynamics is investigated. Gases of interacting bosons are considered, where the full quantum equations to simulate are of a Gross-Pitaevskii form with added Gaussian noise. This method gives tractable simulations of many-body systems because the number of variables scales linearly with the spatial lattice size. An expression for the useful simulation time is obtained, and checked in numerical simulations. The dynamics of first-, second- and third-order spatial correlations are calculated for a uniform interacting 1D Bose gas subjected to a change in scattering length. Propagation of correlations is seen. A comparison is made with other recent methods. The positive P method is particularly well suited to open systems as no conservation laws are hard-wired into the calculation. It also differs from most other recent approaches in that there is no truncation of any kind.

Gaussian operator bases for correlated fermions

We formulate a general multi-mode Gaussian operator basis for fermions, to enable a positive phase-space representation of correlated Fermi states. The Gaussian basis extends existing bosonic phase-space methods to Fermi systems and thus allows first-principles dynamical or equilibrium calculations in quantum many-body Fermi systems. We prove the completeness of the basis and derive differential forms for products with one- and two-body operators. Because the basis satisfies fermionic superselection rules, the resulting phase space involves only c-numbers, without requiring anticommuting Grassmann variables. Furthermore, because of the overcompleteness of the basis, the phase-space distribution can always be chosen positive. This has important consequences for the sign problem in fermion physics.

Many-body quantum dynamics of polarization squeezing in optical fibers

We report new experiments that test quantum dynamical predictions of polarization squeezing for ultrashort photonic pulses in a birefringent fiber, including all relevant dissipative effects. This exponentially complex many-body problem is solved by means of a stochastic phase-space method. The squeezing is calculated and compared to experimental data, resulting in excellent quantitative agreement. From the simulations, we identify the physical limits to quantum noise reduction in optical fibers. The research represents a significant experimental test of first-principles time-domain quantum dynamics in a one-dimensional interacting Bose gas coupled to dissipative reservoirs.

Phase diagram for a class of spin-1/2 Heisenberg models interpolating between the square-lattice, the triangular-lattice, and the linear-chain limits

We study the spin-1/2 Heisenberg models on an anisotropic two-dimensional lattice which interpolates between the square lattice at one end, a set of decoupled spin chains on the other end, and the triangular-lattice Heisenberg model in between. By series expansions around two different dimer ground states and around various commensurate and incommensurate magnetically ordered states, we establish the phase diagram for this model of a frustrated antiferromagnet. We find a particularly rich phase diagram due to the interplay of magnetic frustration, quantum fluctuations, and varying dimensionality. There is a large region of the usual two-sublattice Neel phase, a three-sublattice phase for the triangular-lattice model, a region of incommensurate magnetic order around the triangular-lattice model, and regions in parameter space where there is no magnetic order. We find that the incommensurate ordering wave vector is in general altered from its classical value by quantum fluctuations. The regime of weakly coupled chains is particularly interesting and appears to be nearly critical. [S0163-1829(99)10421-1].

Effect of disorder on quantum phase transitions in anisotropic XY spin chains in a transverse field

We present some exact results for the effect of disorder on the critical properties of an anisotropic XY spin chain in a transverse held. The continuum limit of the corresponding fermion model is taken and in various cases results in a Dirac equation with a random mass. Exact analytic techniques can then be used to evaluate the density of states and the localization length. In the presence of disorder the ferromagnetic-paramagnetic or Ising transition of the model is in the same universality class as the random transverse field Ising model solved by Fisher using a real-space renormalization-group decimation technique (RSRGDT). If there is only randomness in the anisotropy of the magnetic exchange then the anisotropy transition (from a ferromagnet in the x direction to a ferromagnet in the y direction) is also in this universality class. However, if there is randomness in the isotropic part of the exchange or in the transverse held then in a nonzero transverse field the anisotropy transition is destroyed by the disorder. We show that in the Griffiths' phase near the Ising transition that the ground-state energy has an essential singularity. The results obtained for the dynamical critical exponent, typical correlation length, and for the temperature dependence of the specific heat near the Ising transition agree with the results of the RSRODT and numerical work. [S0163-1829(99)07125-8].

Further results on the relationship between mu-invariant measures and quasi-stationary distributions for absorbing continuous-time Markov chains

This note considers continuous-time Markov chains whose state space consists of an irreducible class, C, and an absorbing state which is accessible from C. The purpose is to provide results on mu-invariant and mu-subinvariant measures where absorption occurs with probability less than one. In particular, the well-known premise that the mu-invariant measure, m, for the transition rates be finite is replaced by the more natural premise that m be finite with respect to the absorption probabilities. The relationship between mu-invariant measures and quasi-stationary distributions is discussed. (C) 2000 Elsevier Science Ltd. All rights reserved.

X-ray absorption and phase contrast imaging to study the interplay between plant roots and soil structure

Plant performance is, at least partly, linked to the location of roots with respect to soil structure features and the micro-environment surrounding roots. Measurements of root distributions from intact samples, using optical microscopy and field tracings have been partially successful but are imprecise and labour-intensive. Theoretically, X-ray computed micro-tomography represents an ideal solution for non-invasive imaging of plant roots and soil structure. However, before it becomes fast enough and affordable or easily accessible, there is still a need for a diagnostic tool to investigate root/soil interplay. Here, a method for detection of undisturbed plant roots and their immediate physical environment is presented. X-ray absorption and phase contrast imaging are combined to produce projection images of soil sections from which root distributions and soil structure can be analyzed. The clarity of roots on the X-ray film is sufficient to allow manual tracing on an acetate sheet fixed over the film. In its current version, the method suffers limitations mainly related to (i) the degree of subjectivity associated with manual tracing and (ii) the difficulty of separating live and dead roots. The method represents a simple and relatively inexpensive way to detect and quantify roots from intact samples and has scope for further improvements. In this paper, the main steps of the method, sampling, image acquisition and image processing are documented. The potential use of the method in an agronomic perspective is illustrated using surface and sub-surface soil samples from a controlled wheat trial. Quantitative characterization of root attributes, e.g. radius, length density, branching intensity and the complex interplay between roots and soil structure, is presented and discussed.

Structure-hepatic disposition relationships for cationic drugs in isolated perfused rat livers: Transmembrane exchange and cytoplasmic binding process

This work studied the structure-hepatic disposition relationships for cationic drugs of varying lipophilicity using a single-pass, in situ rat liver preparation. The lipophilicity among the cationic drugs studied in this work is in the following order: diltiazem. propranolol. labetalol. prazosin. antipyrine. atenolol. Parameters characterizing the hepatic distribution and elimination kinetics of the drugs were estimated using the multiple indicator dilution method. The kinetic model used to describe drug transport (the two-phase stochastic model) integrated cytoplasmic binding kinetics and belongs to the class of barrier-limited and space-distributed liver models. Hepatic extraction ratio (E) (0.30-0.92) increased with lipophilicity. The intracellular binding rate constant (k(on)) and the equilibrium amount ratios characterizing the slowly and rapidly equilibrating binding sites (K-S and K-R) increase with the lipophilicity of drug (k(on) : 0.05-0.35 s(-1); K-S : 0.61-16.67; K-R : 0.36-0.95), whereas the intracellular unbinding rate constant (k(off)) decreases with the lipophilicity of drug (0.081-0.021 s(-1)). The partition ratio of influx (k(in)) and efflux rate constant (k(out)), k(in)/k(out), increases with increasing pK(a) value of the drug [from 1.72 for antipyrine (pK(a) = 1.45) to 9.76 for propranolol (pK(a) = 9.45)], the differences in k(in/kout) for the different drugs mainly arising from ion trapping in the mitochondria and lysosomes. The value of intrinsic elimination clearance (CLint), permeation clearance (CLpT), and permeability-surface area product (PS) all increase with the lipophilicity of drug [CLint (ml . min(-1) . g(-1) of liver): 10.08-67.41; CLpT (ml . min(-1) . g(-1) of liver): 10.80-5.35; PS (ml . min(-1) . g(-1) of liver): 14.59-90.54]. It is concluded that cationic drug kinetics in the liver can be modeled using models that integrate the presence of cytoplasmic binding, a hepatocyte barrier, and a vascular transit density function.

Complex natural resonances of conducting planar objects buried in a dielectric half-space

A scheme is presented to incorporate a mixed potential integral equation (MPIE) using Michalski's formulation C with the method of moments (MoM) for analyzing the scattering of a plane wave from conducting planar objects buried in a dielectric half-space. The robust complex image method with a two-level approximation is used for the calculation of the Green's functions for the half-space. To further speed up the computation, an interpolation technique for filling the matrix is employed. While the induced current distributions on the object's surface are obtained in the frequency domain, the corresponding time domain responses are calculated via the inverse fast Fourier transform (FFT), The complex natural resonances of targets are then extracted from the late time response using the generalized pencil-of-function (GPOF) method. We investigate the pole trajectories as we vary the distance between strips and the depth and orientation of single, buried strips, The variation from the pole position of a single strip in a homogeneous dielectric medium was only a few percent for most of these parameter variations.

The use of liquid phase adsorption isotherms for characterization of activated carbons

The characterization of three commercial activated carbons was carried out using the adsorption of various compounds in the aqueous phase. For this purpose the generalized adsorption isotherm was employed, and a modification of the Dubinin-Radushkevich pore filling model, incorporating repulsive contributions to the pore potential as well as bulk liquid phase nonideality, was used as the local isotherm. Eight different flavor compounds were used as adsorbates, and the isotherms were jointly fitted to yield a common pore size distribution for each carbon. The bulk liquid phase nonideality was incorporated through the UNIFAC activity coefficient model, and the repulsive contribution to the pore potential was incorporated through the Steele 10-4-3 potential model. The mean micropore network coordination number for each carbon was also determined from the fitted saturation capacity based on percolation theory. Good agreement between the model and the experimental data was observed. In addition, excellent agreement between the bimodal gamma pore size distribution and density functional theory-cum-regularization-based pore size distribution obtained by argon adsorption was also observed, supporting the validity of the model. The results show that liquid phase adsorption, using adsorptive molecules of different sizes, can be an effective means of characterizing the pore size distribution as well as connectivity. Alternately, if the carbon pore size distribution is independently known, the method can be used to measure critical molecular sizes. (C) 2001 Elsevier Science.

Characterization of activated carbons using liquid phase adsorption

A modification of the Dubinin-Radushkevich pore filling model by incorporation of the repulsive contribution to the pore potential, and of bulk non-ideality, is proposed in this paper for characterization of activated carbon using liquid phase adsorption. For this purpose experiments have been performed using ethyl propionate, ethyl butyrate, and ethyl isovalerate as adsorbates and the microporous-mesoporous activated carbons Filtrasorb 400, Norit ROW 0.8 and Norit ROX 0.8 as adsorbents. The repulsive contribution to the pore potential is incorporated through a Lennard-Jones intermolecular potential model, and the bulk-liquid phase non-ideality through the UNIFAC activity coefficient model. For the characterization of activated carbons, the generalized adsorption isotherm is utilized with a bimodal gamma function as the pore size distribution function. It is found that the model can represent the experimental data very well, and significantly better than when the classical energy-size relationship is used, or when bulk non-ideality is neglected. Excellent agreement between the bimodal gamma pore size distribution and DFT-cum-regularization based pore size distribution is also observed, supporting the validity of the proposed model. (C) 2001 Elsevier Science Ltd. All rights reserved.

Reversible luminescence thermochromism in dipotassiumsodium tris[dicyanoargentate(I)] and the role of structural phase transitions

As a function of temperature, the layered compound K2Na[Ag(CN)213 displays dramatic variations in luminescence thermochromism with major trend changes occurring around 80 K. In order to understand these interesting optical properties, high-resolution neutron diffraction investigations were performed on a polycrystalline sample of this material in the temperature range from 1.5 to 300 K, and previous synchrotron X-ray data of Larochelle et al. (Solid State Commun. 114, 155 (2000)) were reinterpreted. The corresponding significant structural changes were found to be continuous with an anomalous increase of the monoclinic c-lattice parameter with decreasing temperature, associated with slight reorientations of two inequivalent, approximately linear N-C-Ag-C-N units. In the whole temperature range, the crystal structure is monoclinic with the space group C2/m. Based on the structural results, the major luminescence thermochromism changes around 80 K are attributed to the dominance of a back energy transfer process from low- to high-energy excitons at high temperatures. (E) 2002 Elsevier Science (USA).